Showing 151–198 of 198 results for author: Padilla, N

Future virialized structures: An analysis of superstructures in SDSS-DR7

Authors: Heliana E. Luparello, Marcelo Lares, Diego G. Lambas, Nelson D. Padilla

Abstract: We construct catalogues of present superstructures that, according to a LCDM scenario, will evolve into isolated, virialized structures in the future. We use a smoothed luminosity density map derived from galaxies in SDSS-DR7 data and separate high luminosity density peaks. The luminosity density map is obtained from a volume-limited sample of galaxies in the spectroscopic galaxy catalogue, within… ▽ More We construct catalogues of present superstructures that, according to a LCDM scenario, will evolve into isolated, virialized structures in the future. We use a smoothed luminosity density map derived from galaxies in SDSS-DR7 data and separate high luminosity density peaks. The luminosity density map is obtained from a volume-limited sample of galaxies in the spectroscopic galaxy catalogue, within the SDSS-DR7 footprint area and in the redshift range 0.04 < z < 0.12. Other two samples are constructed for calibration and testing purposes, up to z = 0.10 and z = 0.15. The luminosity of each galaxy is spread using an Epanechnikov kernel of 8Mpc/h radius, and the map is constructed on a 1 Mpc/h cubic cells grid. Future virialized structures (FVS) are identified as regions with overdensity above a given threshold, calibrated using a LCDM numerical simulation, and the criteria presented by Dünner et al. (2006). We assume a constant mass-to-luminosity ratio and impose the further condition of a minimum luminosity of 10^{12}Lsol. According to our calibrations with a numerical simulation, these criteria lead to a negligible contamination by less overdense (non FVS) superstructures.We present a catalogue of superstructures in the SDSS-DR7 area within redshift 0.04 < z < 0.12 and test the reliability of our method by studying different subsamples as well as a mock catalogue.We compute the luminosity and volume distributions of the superstructures finding that about 10% of the luminosity (mass) will end up in future virialized structures. The fraction of groups and X-ray clusters in these superstructures is higher for groups/clusters of higher mass, suggesting that future cluster mergers will involve the most massive systems. We also analyse known structures in the present Universe and compare with our catalogue of FVS. △ Less

Submitted 15 April, 2011; v1 submitted 10 January, 2011; originally announced January 2011.

Comments: 14 pages, 11 figures, modified to match accepted version in MNRAS. PDF with high resolution colour figures is available at http://www.oac.uncor.edu/apache2-default/adminweb/html/WEB/preprints/2011.01/FVS-DR7.pdf

Host galaxy-active galactic nucleus alignments in the SDSS-DR7

Authors: Claudia del P. Lagos, Nelson D. Padilla, Michael A. Strauss, Sofia A. Cora, Lei Hao

Abstract: We determine the intrinsic shapes and orientations of 27,450 type I and II active galactic nucleus (AGN) galaxies in the spectroscopic sample of the SDSS DR7, by studying the distribution of projected axis ratios of AGN hosts. Our aim is to study possible alignments between the AGN and host galaxy systems (e.g. the accretion disc and the galaxy angular momentum) and the effect of dust obscuration… ▽ More We determine the intrinsic shapes and orientations of 27,450 type I and II active galactic nucleus (AGN) galaxies in the spectroscopic sample of the SDSS DR7, by studying the distribution of projected axis ratios of AGN hosts. Our aim is to study possible alignments between the AGN and host galaxy systems (e.g. the accretion disc and the galaxy angular momentum) and the effect of dust obscuration geometry on the AGN type. We define control samples of non-AGN galaxies that mimic the morphology, colour, luminosity and concentration distributions of the AGN population, taking into account the effects of dust extinction and reddening. Assuming that AGN galaxies have the same underlying three-dimensional shape distribution as their corresponding control samples, we find that the spiral and elliptical type I AGN populations are strongly skewed toward face-on galaxies, while ellipticals and spirals type II AGN are skewed toward edge-on orientations. These findings rule out random orientations for AGN hosts at high confidence for type I spirals (delta chi^2~230) and type II ellipticals (delta chi^2~15), while the signal for type I ellipticals and type II spirals is weaker (delta chi^2~3 and delta chi^2~6, respectively). We obtain a much stronger tendency for the type II spirals to be edge-on when just high [OIII] equivalent width (EW) AGN are considered, suggesting that >20% of low [OIII] EW edge-on type II AGN may be missing from the optical sample. Galactic dust absorption of the broad-line region alone cannot explain the observed inclination angle and projected axis ratio distributions of type I and II Seyfert types, implying that obscuration by a small-scale circumnuclear torus is necessary. These results favour a scenario in which the angular momentum of the material which feeds the black hole retains a memory of its original gas source at least to some small, non-negligible degree. △ Less

Submitted 11 February, 2011; v1 submitted 19 November, 2010; originally announced November 2010.

Comments: 17 pages, 11 figures. Accepted in MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.414:2148-2162,2011

Large scale anisotropies on halo infall

Authors: Laura Ceccarelli, Dante J. Paz, Nelson Padilla, Diego G. Lambas

Abstract: We perform a statistical analysis of the peculiar velocity field around dark matter haloes in numerical simulations. We examine different properties of the infall of material onto haloes and its relation to central halo shapes and the shape of the large scale surrounding regions (LSSR). We find that the amplitude of the infall velocity field along the halo shape minor axis is larger than that alon… ▽ More We perform a statistical analysis of the peculiar velocity field around dark matter haloes in numerical simulations. We examine different properties of the infall of material onto haloes and its relation to central halo shapes and the shape of the large scale surrounding regions (LSSR). We find that the amplitude of the infall velocity field along the halo shape minor axis is larger than that along the major axis. This is consistent for general triaxial haloes, and for both prolate and oblate systems. We also report a strong anisotropy of the velocity field along the principal axes of the LSSR. The infall velocity field around dark matter haloes reaches a maximum value along the direction of the minor axis of the LSSR, whereas along the direction of its major axis, it exhibits the smallest velocities. We also analyse the dependence of the matter velocity field on the local environment. The amplitude of the infall velocity at high local density regions is larger than at low local density regions. The velocity field tends to be more laminar along the direction towards the minor axis of the LSSR, where the mean ratio between flow velocity and velocity dispersion is of order unity and nearly constant up to scales of 15 Mpc/h. We also detect anisotropies in the outflowing component of the velocity field, showing a maximum amplitude along the surrounding LSSR major axis. △ Less

Submitted 10 November, 2010; originally announced November 2010.

Comments: 12 pages, 9 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society

The nature of assembly bias - I. Clues from a LCDM cosmology

Authors: Ivan Lacerna, Nelson Padilla

Abstract: We present a new proxy for the overdensity peak height for which the large-scale clustering of haloes of a given mass does not vary significantly with the assembly history. The peak height, usually taken to be well represented by the virial mass, can instead be approximated by the mass inside spheres of different radii, which in some cases can be larger than the virial radius and therefore include… ▽ More We present a new proxy for the overdensity peak height for which the large-scale clustering of haloes of a given mass does not vary significantly with the assembly history. The peak height, usually taken to be well represented by the virial mass, can instead be approximated by the mass inside spheres of different radii, which in some cases can be larger than the virial radius and therefore include mass outside the individual host halo. The sphere radii are defined as r = $a$ delta_t + $b$ log_10(M_vir/M_nl), where delta_t is the age relative to the typical age of galaxies hosted by haloes with virial mass M_vir, M_nl is the non-linear mass, and $a$=0.2 and $b$=-0.02 are the free parameters adjusted to trace the assembly bias effect. Note that $r$ depends on both halo mass and age. In this new approach, some of the objects which were initially considered low-mass peaks belong to regions with higher overdensities. At large scales, i.e. in the two-halo regime, this model properly recovers the simple prescription where the bias responds to the height of the mass peak alone, in contrast to the usual definition (virial mass) that shows a strong dependence on additional halo properties such as formation time. The dependence on the age in the one-halo term is also remarkably reduced. The population of galaxies whose "peak height" changes with this new definition consists mainly of old stellar populations and are preferentially hosted by low-mass haloes located near more massive objects. The latter is in agreement with recent results which indicate that old, low-mass haloes would suffer truncation of mass accretion by nearby larger haloes or simply due to the high density of their surroundings, thus showing an assembly bias effect. The change in mass is small enough that the Sheth et al. (2001) mass function is still a good fit to the resulting distribution of new masses. △ Less

Submitted 5 November, 2010; originally announced November 2010.

Comments: 13 pages, 10 figures, submitted to MNRAS, comments welcome

Low Surface Brightness Galaxies in the SDSS: the link between environment, star-forming properties and AGN

Authors: Gaspar Galaz, Rodrigo Herrera-Camus, Diego Garcia-Lambas, Nelson Padilla

Abstract: Using the Sloan Digital Sky Survey (SDSS) data release 4 (DR 4), we investigate the spatial distribution of low and high surface brightness galaxies (LSBGs and HSBGs, respectively). In particular, we focus our attention on the influence of interactions between galaxies on the star formation strength in the redshift range $0.01 < z < 0.1$. With cylinder counts and projected distance to the first an… ▽ More Using the Sloan Digital Sky Survey (SDSS) data release 4 (DR 4), we investigate the spatial distribution of low and high surface brightness galaxies (LSBGs and HSBGs, respectively). In particular, we focus our attention on the influence of interactions between galaxies on the star formation strength in the redshift range $0.01 < z < 0.1$. With cylinder counts and projected distance to the first and fifth-nearest neighbor as environment tracers, we find that LSBGs tend to have a lack of companions compared to HSBGs at small scales ($<2$ Mpc). Regarding the interactions, we have evidence that the fraction of LSBGs with strong star formation activity increases when the distance between pairs of galaxies ($r_{p}$) is smaller than about four times the Petrosian radius ($r_{90}$) of one of the components. Our results suggest that, rather than being a condition for their formation, the isolation of LSBGs is more connected with their survival and evolution. The effect of the interaction on the star formation strength, measured by the average value of the birthrate parameter $b$, seems to be stronger for HSBGs than for LSBGs. The analysis of our population of LSBGs and HSBGs hosting an AGN show that, regardless of the mass range, the fraction of LSBGs having an AGN is lower than the corresponding fraction of HSBGs with an AGN. Also, we observe that the fraction of HSBGs and LSBGs having an AGN increases with the bulge luminosity. These results, and those concerning the star-forming properties of LSBGs as a function of the environment, fit with the scenario proposed by some authors where, below a given threshold of surface mass density, low surface brightness disks are unable to propagate instabilities, preventing the formation and evolution of massive black holes in the centers of LSBGs. △ Less

Submitted 12 December, 2010; v1 submitted 22 July, 2010; originally announced July 2010.

Comments: 33 pages, 13 Figures, 2 Tables. Accepted for publication in The Astrophysical Journal (January 2011 Issue)

Clustering and descendants of MUSYC galaxies at z<1.5

Authors: Nelson Padilla, Daniel Christlein, Eric Gawiser, Roberto González, Lucia Guaita, Leopoldo Infante

Abstract: We measure the evolution of galaxy clustering out to a redshift of z~1.5 using data from two MUSYC fields, the Extended Hubble Deep Field South (EHDF-S) and the Extended Chandra Deep Field South (ECDF-S). We use photometric redshift information to calculate the projected-angular correlation function, omega(sigma), from which we infer the projected correlation function Xi(sigma). We demonstrate tha… ▽ More We measure the evolution of galaxy clustering out to a redshift of z~1.5 using data from two MUSYC fields, the Extended Hubble Deep Field South (EHDF-S) and the Extended Chandra Deep Field South (ECDF-S). We use photometric redshift information to calculate the projected-angular correlation function, omega(sigma), from which we infer the projected correlation function Xi(sigma). We demonstrate that this technique delivers accurate measurements of clustering even when large redshift measurement errors affect the data. To this aim we use two mock MUSYC fields extracted from a LambdaCDM simulation populated with GALFORM semi-analytic galaxies which allow us to assess the degree of accuracy of our estimates of Xi(sigma) and to identify and correct for systematic effects in our measurements. We study the evolution of clustering for volume limited subsamples of galaxies selected using their photometric redshifts and rest-frame r-band absolute magnitudes. We find that the real-space correlation length r_0 of bright galaxies, M_r<-21 (rest-frame) can be accurately recovered out to z~1.5, particularly for ECDF-S given its near-infrared photometric coverage. There is mild evidence for a luminosity dependent clustering in both fields at the low redshift samples (up to <z>=0.57), where the correlation length is higher for brighter galaxies by up to 1Mpc/h between median rest-frame r-band absolute magnitudes of -18 to -21.5. As a result of the photometric redshift measurement, each galaxy is assigned a best-fit template; we restrict to E and E+20%Sbc types to construct subsamples of early type galaxies (ETGs). Our ETG samples show a strong increase in r_0 as the redshift increases, making it unlikely (95% level) that ETGs at median redshift z_med=1.15 are the direct progenitors of ETGs at z_med=0.37 with equivalent passively evolved luminosities. (ABRIDGED) △ Less

Submitted 7 July, 2010; v1 submitted 29 June, 2010; originally announced June 2010.

Comments: 16 pages, 12 figures, 2 tables, accepted for publication in MNRAS

Journal ref: 2010, MNRAS, 409, 184

VLT LBG Redshift Survey II: Interactions between galaxies and the IGM at z ~3

Authors: N. H. M. Crighton, R. Bielby, T. Shanks, L. Infante, C. G. Bornancini, N. Bouche, D. G. Lambas, J. D. Lowenthal, D. Minniti, S. L. Morris, N. Padilla, C. Peroux, P. Petitjean, T. Theuns, P. Tummuangpak, P. M. Weilbacher, L. Wisotzki, G. Worseck

Abstract: We have measured redshifts for 243 z ~3 quasars in nine VLT VIMOS LBG redshift survey areas, each of which is centred on a known bright quasar. Using spectra of these quasars, we measure the cross-correlation between neutral hydrogen gas causing the Lya forest and 1020 Lyman-break galaxies at z ~3. We find an increase in neutral hydrogen absorption within 5 h^-1 Mpc of a galaxy in agreement with t… ▽ More We have measured redshifts for 243 z ~3 quasars in nine VLT VIMOS LBG redshift survey areas, each of which is centred on a known bright quasar. Using spectra of these quasars, we measure the cross-correlation between neutral hydrogen gas causing the Lya forest and 1020 Lyman-break galaxies at z ~3. We find an increase in neutral hydrogen absorption within 5 h^-1 Mpc of a galaxy in agreement with the results of Adelberger et al. (2003, 2005). The Lya-LBG cross-correlation can be described by a power-law on scales larger than 3 h^-1 Mpc. When galaxy velocity dispersions are taken into account our results at smaller scales (<2 h^-1 Mpc) are also in good agreement with the results of Adelberger et al. (2005). There is little immediate indication of a region with a transmission spike above the mean IGM value which might indicate the presence of star-formation feedback. To measure the galaxy velocity dispersions, which include both intrinsic LBG velocity dispersion and redshift errors, we have used the LBG-LBG redshift space distortion measurements of Bielby et al. (2010). We find that the redshift-space transmission spike implied in the results of Adelberger et al. (2003) is too narrow to be physical in the presence of the likely LBG velocity dispersion and is likely to be a statistical fluke. Nevertheless, neither our nor previous data can rule out the presence of a narrow, real-space transmission spike, given the evidence of the increased Lya absorption surrounding LBGs which can mask the spike's presence when convolved with a realistic LBG velocity dispersion. Finally, we identify 176 CIV systems in the quasar spectra and find an LBG-CIV correlation strength on scales of 10 h^-1 Mpc consistent with the relation measured at ~Mpc scales. △ Less

Submitted 22 June, 2010; originally announced June 2010.

Comments: 33 pages, 16 figures, accepted by MNRAS

The LABOCA survey of the Extended Chandra Deep Field South: A photometric redshift survey of submillimetre galaxies

Authors: J. L. Wardlow, Ian Smail, K. E. K. Coppin, D. M. Alexander, W. N. Brandt, A. L. R. Danielson, B. Luo, A. M. Swinbank, F. Walter, A. Weiss, Y. Q. Xue, S. Zibetti, F. Bertoldi, A. D. Biggs, S. C. Chapman, H. Dannerbauer, J. S. Dunlop, E. Gawiser, R. J. Ivison, K. K. Knudsen, A. Kovacs, C. G. Lacey, K. M. Menten, N. Padilla, H. -W. Rix , et al. (1 additional authors not shown)

Abstract: [abridged] We derive photometric redshifts from 17-band optical to mid-IR photometry of 74 robust counterparts to 68 of the 126 submillimetre galaxies (SMGs) selected at 870um by LABOCA observations in the ECDFS. The median photometric redshift of identified SMGs is z=2.2\pm0.1, the interquartile range is z=1.8-2.7 and we identify 10 (~15%) high-redshift (z>3) SMGs. We derive a simple redshift est… ▽ More [abridged] We derive photometric redshifts from 17-band optical to mid-IR photometry of 74 robust counterparts to 68 of the 126 submillimetre galaxies (SMGs) selected at 870um by LABOCA observations in the ECDFS. The median photometric redshift of identified SMGs is z=2.2\pm0.1, the interquartile range is z=1.8-2.7 and we identify 10 (~15%) high-redshift (z>3) SMGs. We derive a simple redshift estimator for SMGs based on the 3.6 and 8um fluxes, which is accurate to Delta_z~0.4 for SMGs at z<4. A statistical analysis of sources around unidentified SMGs identifies a population of likely counterparts with a redshift distribution peaking at z=2.5\pm0.3, which likely comprises ~60% of the unidentified SMGs. This confirms that the bulk of the undetected SMGs are coeval with those detected in the radio/mid-IR. We conclude that ~15% of all the SMGs are below the flux limits of our survey and lie at z>3 and hence ~30% of all SMGs have z>3. We estimate that the full S_870um>4mJy SMG population has a median redshift of 2.5\pm0.6. In contrast to previous suggestions we find no significant correlation between S_870um and redshift. The median stellar mass of the SMGs derived from SED fitting is (9.2\pm0.9)x10^10Msun and the interquartile range is (4.7-14)x10^10Msun, although we caution that uncertainty in the star-formation histories results in a factor of ~5 uncertainty in these stellar masses. The median characteristic dust temperature of SMGs is 35.9\pm1.4K and the interquartile range is 28.5-43.3K. The infrared luminosity function shows that SMGs at z=2-3 typically have higher far-IR luminosities and luminosity density than those at z=1-2. This is mirrored in the evolution of the star-formation rate density (SFRD) for SMGs which peaks at z~2. The maximum contribution of bright SMGs to the global SFRD (~5% for SMGs with S_870um>4mJy; ~50% for SMGs with S_870um>1mJy) also occurs at z~2. △ Less

Submitted 10 June, 2010; originally announced June 2010.

Comments: Submitted for publication in MNRAS: 31 pages, 15 figures, 4 tables

Journal ref: 2011, MNRAS, 415, 1497

The VLT LBG Redshift Survey I: Clustering and Dynamics of ~1000 Galaxies at z~3

Authors: Rich Bielby, T. Shanks, P. M. Weilbacher, L. Infante, N. H. M. Crighton, Carlos Bornancini, N. Bouché, P. Héraudeau, D. G. Lambas, J. Lowenthal, D. Minniti, N. Padilla, P. Petitjean, T. Theuns

Abstract: We present the initial imaging and spectroscopic data acquired as part of the VLT VIMOS Lyman-break galaxy Survey. UBR (or UBVI) imaging covers five 36'x36' fields centred on bright z>3 QSOs, allowing ~21,000 2<z<3.5 galaxy candidates to be selected using the Lyman-break technique. We performed spectroscopic follow-up using VIMOS, measuring redshifts for 1020 z>2 LBGs and 10 z>2 QSOs from a total… ▽ More We present the initial imaging and spectroscopic data acquired as part of the VLT VIMOS Lyman-break galaxy Survey. UBR (or UBVI) imaging covers five 36'x36' fields centred on bright z>3 QSOs, allowing ~21,000 2<z<3.5 galaxy candidates to be selected using the Lyman-break technique. We performed spectroscopic follow-up using VIMOS, measuring redshifts for 1020 z>2 LBGs and 10 z>2 QSOs from a total of 19 VIMOS pointings. From the galaxy spectra, we observe a 625+-510 km/s velocity offset between the ISM absorption and Ly-alpha emission line redshifts. Using the photometric and spectroscopic catalogues, we have analysed the galaxy clustering at z~3. In the photometric case, the angular correlation function, w(theta), is well fit by a double power-law with clustering scale-length, r_0 = 3.19+0.32-0.54 Mpc/h for r < 1 Mpc/h and r_0 = 4.59+0.31-0.33 Mpc/h at larger scales. Using the redshift sample we estimate the semi-projected correlation function, w_p(sigma) and find r_0 = 3.67+0.23-0.24 Mpc/h for the VLT sample and r_0 = 3.98+0.14-0.15 Mpc/h for a combined VLT+Keck sample. From the z-space correlation functions and assuming the above xi(r) models, we find that the combined VLT and Keck surveys require a galaxy velocity dispersion, <w_z^2>^1/2 ~ 700 km/s, higher than the ~400 km/s found by previous authors. We also measure a value for the gravitational growth rate parameter of beta(z=3) = 0.48+-0.17, implying a low value for the bias of b = 2.06+1.1-0.5. This value is consistent with the galaxy clustering amplitude which gives b = 2.22+-0.16, assuming the standard cosmology, implying that the evolution of the gravitational growth rate is also consistent with Einstein gravity. We have compared our LBG clustering amplitudes with lower redshift measurements and find that the clustering strength is not inconsistent with that of low-redshift L* spirals for simple 'long-lived' galaxy models. △ Less

Submitted 4 April, 2011; v1 submitted 17 May, 2010; originally announced May 2010.

Comments: 27 pages, 25 figures, accepted in MNRAS

Unveiling the Nature of Submillimeter Galaxy SXDF850.6

Authors: B. Hatsukade, D. Iono, T. Yoshikawa, M. Akiyama, J. S. Dunlop, R. J. Ivison, A. B. Peck, S. Ikarashi, A. Biggs, H. Ezawa, H. Hanami, P. Ho, D. H. Hughes, R. Kawabe, K. Kohno, S. Matsushita, K. Nakanishi, N. Padilla, G. Petitpas, Y. Tamura, J. Wagg, D. J. Wilner, G. W. Wilson, T. Yamada, M. S. Yun

Abstract: We present an 880 micron Submillimeter Array (SMA) detection of the submillimeter galaxy SXDF850.6. SXDF850.6 is a bright source (S(850 micron) = 8 mJy) detected in the SCUBA Half Degree Extragalactic Survey (SHADES), and has multiple possible radio counterparts in its deep radio image obtained at the VLA. Our new SMA detection finds that the submm emission coincides with the brightest radio emi… ▽ More We present an 880 micron Submillimeter Array (SMA) detection of the submillimeter galaxy SXDF850.6. SXDF850.6 is a bright source (S(850 micron) = 8 mJy) detected in the SCUBA Half Degree Extragalactic Survey (SHADES), and has multiple possible radio counterparts in its deep radio image obtained at the VLA. Our new SMA detection finds that the submm emission coincides with the brightest radio emission that is found ~8" north of the coordinates determined from SCUBA. Despite the lack of detectable counterparts in deep UV/optical images, we find a source at the SMA position in near-infrared and longer wavelength images. We perform SED model fits to UV-optical-IR photometry (u, B, V, R, i', z', J, H, K, 3.6 micron, 4.5 micron, 5.8 micron, and 8.0 micron) and to submm-radio photometry (850 micron, 880 micron, 1100 micron, and 21 cm) independently, and we find both are well described by starburst templates at a redshift of z ~= 2.2 (+/- 0.3). The best-fit parameters from the UV-optical-IR SED fit are a redshift of z = 1.87 (+0.15/-0.07), a stellar mass of M_star = 2.5 +2.2/-0.3 x 10^11 M_sun, an extinction of A_V = 3.0 (+0.3/-1.0) mag, and an age of 720 (+1880/-210) Myr. The submm-radio SED fit provides a consistent redshift of z ~ 1.8-2.5, an IR luminosity of L_IR = (7-26) x 10^12 L_sun, and a star formation rate of 1300-4500 M_sun/yr. These results suggest that SXDF850.6 is a mature system already having a massive amount of old stellar population constructed before its submm bright phase and is experiencing a dusty starburst, possibly induced by major mergers. △ Less

Submitted 20 January, 2010; originally announced January 2010.

Comments: 7 pages, 5 figures, Accepted for publication in Astrophysical Journal

VISTA Variables in the Via Lactea (VVV): The public ESO near-IR variability survey of the Milky Way

Authors: D. Minniti, P. W. Lucas, J. P. Emerson, R. K. Saito, M. Hempel, P. Pietrukowicz, A. V. Ahumada, M. V. Alonso, J. Alonso-García, J. I. Arias, R. M. Bandyopadhyay, R. H. Barbá, B. Barbuy, L. R. Bedin, E. Bica, J. Borissova, L. Bronfman, G. Carraro, M. Catelan, J. J. Clariá, N. Cross, R. de Grijs, I. Dékány, J. E. Drew, C. Fariña , et al. (41 additional authors not shown)

Abstract: We describe the public ESO near-IR variability survey (VVV) scanning the Milky Way bulge and an adjacent section of the mid-plane where star formation activity is high. The survey will take 1929 hours of observations with the 4-metre VISTA telescope during five years (2010-2014), covering ~10^9 point sources across an area of 520 deg^2, including 33 known globular clusters and ~350 open clusters… ▽ More We describe the public ESO near-IR variability survey (VVV) scanning the Milky Way bulge and an adjacent section of the mid-plane where star formation activity is high. The survey will take 1929 hours of observations with the 4-metre VISTA telescope during five years (2010-2014), covering ~10^9 point sources across an area of 520 deg^2, including 33 known globular clusters and ~350 open clusters. The final product will be a deep near-IR atlas in five passbands (0.9-2.5 microns) and a catalogue of more than 10^6 variable point sources. Unlike single-epoch surveys that, in most cases, only produce 2-D maps, the VVV variable star survey will enable the construction of a 3-D map of the surveyed region using well-understood distance indicators such as RR Lyrae stars, and Cepheids. It will yield important information on the ages of the populations. The observations will be combined with data from MACHO, OGLE, EROS, VST, Spitzer, HST, Chandra, INTEGRAL, WISE, Fermi LAT, XMM-Newton, GAIA and ALMA for a complete understanding of the variable sources in the inner Milky Way. This public survey will provide data available to the whole community and therefore will enable further studies of the history of the Milky Way, its globular cluster evolution, and the population census of the Galactic Bulge and center, as well as the investigations of the star forming regions in the disk. The combined variable star catalogues will have important implications for theoretical investigations of pulsation properties of stars. △ Less

Submitted 11 December, 2009; v1 submitted 5 December, 2009; originally announced December 2009.

Comments: accepted for publication in New Astronomy

Journal ref: New Astron.15:433-443,2010

LSST Science Book, Version 2.0

Authors: LSST Science Collaboration, Paul A. Abell, Julius Allison, Scott F. Anderson, John R. Andrew, J. Roger P. Angel, Lee Armus, David Arnett, S. J. Asztalos, Tim S. Axelrod, Stephen Bailey, D. R. Ballantyne, Justin R. Bankert, Wayne A. Barkhouse, Jeffrey D. Barr, L. Felipe Barrientos, Aaron J. Barth, James G. Bartlett, Andrew C. Becker, Jacek Becla, Timothy C. Beers, Joseph P. Bernstein, Rahul Biswas, Michael R. Blanton, Joshua S. Bloom , et al. (223 additional authors not shown)

Abstract: A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south… ▽ More A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy. △ Less

Submitted 1 December, 2009; originally announced December 2009.

Comments: 596 pages. Also available at full resolution at http://www.lsst.org/lsst/scibook

Automated detection of filaments in the large scale structure of the universe

Authors: Roberto E. Gonzalez, Nelson E. Padilla

Abstract: We present a new method to identify large scale filaments and apply it to a cosmological simulation. Using positions of haloes above a given mass as node tracers, we look for filaments between them using the positions and masses of all the remaining dark-matter haloes. In order to detect a filament, the first step consists in the construction of a backbone linking two nodes, which is given by a… ▽ More We present a new method to identify large scale filaments and apply it to a cosmological simulation. Using positions of haloes above a given mass as node tracers, we look for filaments between them using the positions and masses of all the remaining dark-matter haloes. In order to detect a filament, the first step consists in the construction of a backbone linking two nodes, which is given by a skeleton-like path connecting the highest local dark matter (DM) density traced by non-node haloes. The filament quality is defined by a density and gap parameters characterising its skeleton, and filament members are selected by their binding energy in the plane perpendicular to the filament. This membership condition is associated to characteristic orbital times; however if one assumes a fixed orbital timescale for all the filaments, the resulting filament properties show only marginal changes, indicating that the use of dynamical information is not critical for the method. We test the method in the simulation using massive haloes($M>10^{14}$h$^{-1}M_{\odot}$) as filament nodes. The main properties of the resulting high-quality filaments (which corresponds to $\simeq33%$ of the detected filaments) are, i) their lengths cover a wide range of values of up to $150 $h$^{-1}$Mpc, but are mostly concentrated below 50h$^{-1}$Mpc; ii) their distribution of thickness peaks at $d=3.0$h$^{-1}$Mpc and increases slightly with the filament length; iii) their nodes are connected on average to $1.87\pm0.18$ filaments for $\simeq 10^{14.1}M_{\odot}$ nodes; this number increases with the node mass to $\simeq 2.49\pm0.28$ filaments for $\simeq 10^{14.9}M_{\odot}$ nodes. △ Less

Submitted 13 May, 2010; v1 submitted 1 December, 2009; originally announced December 2009.

Comments: 17 pages, 13 figures, MNRAS Accepted

Journal ref: Mon.Not.Roy.Astron.Soc.407:1449-1463,2010

Local and global environmental effects on galaxies and active galactic nuclei

Authors: Nelson Padilla, Diego Garcia Lambas, Roberto Gonzalez

Abstract: We study the properties of SDSS galaxies with and without AGN detection as a function of the local and global environment measured via the local density, the mass of the galaxy host group (parameterised by the group luminosity) and distance to massive clusters. Our results can be divided in two main subjects, the environments of galaxies and their relation to the assembly of their host haloes, and… ▽ More We study the properties of SDSS galaxies with and without AGN detection as a function of the local and global environment measured via the local density, the mass of the galaxy host group (parameterised by the group luminosity) and distance to massive clusters. Our results can be divided in two main subjects, the environments of galaxies and their relation to the assembly of their host haloes, and the environments of AGN. (i) For the full SDSS sample, we find indications that the local galaxy density is the most efficient parameter to separate galaxy populations, but we also find that galaxies at fixed local density show some remaining variation of their properties as a function of the distance to the nearest cluster of galaxies (in a range of 0 to 10 cluster virial radii). These differences seem to become less significant if the galaxy samples are additionally constrained to be hosted by groups of similar total luminosity. (ii) In AGN host galaxies, the morphology-density relation is much less noticeable when compared to the behaviour of the full SDSS sample. In order to interpret this result we analyse control samples constructed using galaxies with no detected AGN activity with matching distributions of redshifts, stellar masses, r-band luminosities, g-r colours, concentrations, local densities, host group luminosities, and fractions of central and satellite galaxies. The control samples also show a similar small dependence on the local density indicating an influence from the AGN selection, but their colours are slightly bluer compared to the AGN hosts regardless of local density. Furthermore, even when the local density is held fixed at intermediate or high values, and the distance to the closest cluster of galaxies is allowed to vary, AGN control galaxies away from clusters tend to be bluer than the AGN hosts. (ABRIDGED) △ Less

Submitted 17 November, 2010; v1 submitted 27 November, 2009; originally announced November 2009.

Comments: 18 pages, 7 figures, MNRAS, in press

Journal ref: Mon.Not.Roy.Astron.Soc.409:936,2010

Lyman-Alpha-Emitting Galaxies at z = 2.1 in ECDF-S: Building Blocks of Typical Present-day Galaxies?

Authors: L. Guaita, E. Gawiser, N. Padilla, H. Francke, N. A. Bond, C. Gronwall, R. Ciardullo, J. J. Feldmeier, S. Sinawa, G. A. Blanc, S. Virani

Abstract: We discovered a sample of 250 Ly-Alpha emitting (LAE) galaxies at z=2.1 in an ultra-deep 3727 A narrow-band MUSYC image of the Extended Chandra Deep Field-South. LAEs were selected to have rest-frame equivalent widths (EW) > 20 A and emission line fluxes > 2.0 x 10^(-17)erg /cm^2/s, after carefully subtracting the continuum contributions from narrow band photometry. The median flux of our sample… ▽ More We discovered a sample of 250 Ly-Alpha emitting (LAE) galaxies at z=2.1 in an ultra-deep 3727 A narrow-band MUSYC image of the Extended Chandra Deep Field-South. LAEs were selected to have rest-frame equivalent widths (EW) > 20 A and emission line fluxes > 2.0 x 10^(-17)erg /cm^2/s, after carefully subtracting the continuum contributions from narrow band photometry. The median flux of our sample is 4.2 x 10^(-17)erg/cm^2/s, corresponding to a median Lya luminosity = 1.3 x 10^(42) erg/s at z=2.1. At this flux our sample is > 90% complete. Approximately 4% of the original NB-selected candidates were detected in X-rays by Chandra, and 7% were detected in the rest-frame far-UV by GALEX. At luminosity>1.3 x 10^42 erg/s, the equivalent width distribution is unbiased and is represented by an exponential with scale-length of 83+/-10 A. Above this same luminosity threshold, we find a number density of 1.5+/-0.5 x 10^-3 Mpc^-3. Neither the number density of LAEs nor the scale-length of their EW distribution show significant evolution from z=3 to z=2. We used the rest frame UV luminosity to estimate a median star formation rate of 4 M_(sun) /yr. The median rest frame UV slope, parametrized by B-R, is that typical of dust-free, 0.5-1 Gyr old or moderately dusty, 300-500 Myr old populations. Approximately 40% of the sample occupies the z~2 star-forming galaxy locus in the UVR two color diagram. Clustering analysis reveals that LAEs at z=2.1 have r_0=4.8+/-0.9 Mpc and a bias factor b=1.8+/-0.3. This implies that z=2.1 LAEs reside in dark matter halos with median masses Log(M/M_(sun))=11.5^(+0.4)_(-0.5), which are among of the lowest-mass halos yet probed at this redshift. We used the Sheth-Tormen conditional mass function to study the descendants of these LAEs and found that their typical present-day descendants are local galaxies with L* properties, like the Milky Way. △ Less

Submitted 9 March, 2010; v1 submitted 12 October, 2009; originally announced October 2009.

Comments: 35 pages, 9 figures, ApJ, in press

Journal ref: Astrophys.J.714:255-269,2010

Introducing the Photometric Maximum Likelihood Method: Galaxy Luminosity Functions at z<1.2 in MUSYC-ECDFS

Authors: Daniel Christlein, Eric Gawiser, Danilo Marchesini, Nelson Padilla

Abstract: We present a new maximum likelihood method for the calculation of galaxy luminosity functions from multi-band photometric surveys without spectroscopic data. The method evaluates the likelihood of a trial luminosity function by directly comparing the predicted distribution of fluxes in a multi-dimensional photometric space to the observations, and thus does not require the intermediate step of c… ▽ More We present a new maximum likelihood method for the calculation of galaxy luminosity functions from multi-band photometric surveys without spectroscopic data. The method evaluates the likelihood of a trial luminosity function by directly comparing the predicted distribution of fluxes in a multi-dimensional photometric space to the observations, and thus does not require the intermediate step of calculating photometric redshifts. We apply this algorithm to ~27,000 galaxies with m_R<=25 in the MUSYC-ECDFS field, with a focus on recovering the luminosity function of field galaxies at z<1.2. Our deepest LFs reach M_r=-14 and show that the field galaxy LF deviates from a Schechter function, exhibiting a steep upturn at intermediate magnitudes that is due to galaxies of late spectral types. △ Less

Submitted 10 September, 2009; originally announced September 2009.

Comments: 24 pages, 11 figures, accepted for publication by MNRAS

Global environmental effects versus galaxy interactions

Authors: Josefa Perez, Patricia Tissera, Nelson Padilla, M. Sol Alonso, Diego G. Lambas, .

Abstract: We explore properties of close galaxy pairs and merging systems selected from the SDSS-DR4 in different environments with the aim to assess the relative importance of the role of interactions over global environmental processes. For this purpose, we perform a comparative study of galaxies with and without close companions as a function of local density and host-halo mass, carefully removing sour… ▽ More We explore properties of close galaxy pairs and merging systems selected from the SDSS-DR4 in different environments with the aim to assess the relative importance of the role of interactions over global environmental processes. For this purpose, we perform a comparative study of galaxies with and without close companions as a function of local density and host-halo mass, carefully removing sources of possible biases. We find that at low and high local density environments, colours and morphologies of close galaxy pairs are very similar to those of isolated galaxies. At intermediate densities, we detect significant differences, indicating that close pairs could have experienced a more rapid transition onto the red sequence than isolated galaxies. The presence of a correlation between colours and morphologies indicates that the physical mechanism responsible for the colour transformation also operates changing galaxy morphologies. Regardless of dark matter halo mass, we show that the percentage of red galaxies in close pairs and in the control sample are comparable at low and high local density environments. However, at intermediate local densities, the gap in the red fraction between close pairs and the control galaxies increases from ~10% in low mass haloes up to ~50% in the most massive ones. Our findings suggest that in intermediate density environments galaxies are efficiently pre-processed by close encounters and mergers before entering higher local density regions. (Abridge) △ Less

Submitted 21 April, 2009; v1 submitted 18 April, 2009; originally announced April 2009.

Comments: submitted to MNRAS, 12 pages, 9 figures (For people who have tried to download this paper, I've already changed the wrong version)

Galaxy Clusters in the Line of Sight to Background Quasars: II. Environmental effects on the sizes of baryonic halo sizes

Authors: N. Padilla, I. Lacerna, S. Lopez, L. F. Barrientos, P. Lira, H. Andrews, N. Tejos

Abstract: Based on recent results on the frequency of MgII absorption line systems in the "QSO behind RCS clusters" survey (QbC), we analyse the effects of the cluster environment on the sizes of baryonic haloes around galaxies. We use two independent models, i) an empirical halo occupation model which fits current measurements of the clustering and luminosity function of galaxies at low and high redshift… ▽ More Based on recent results on the frequency of MgII absorption line systems in the "QSO behind RCS clusters" survey (QbC), we analyse the effects of the cluster environment on the sizes of baryonic haloes around galaxies. We use two independent models, i) an empirical halo occupation model which fits current measurements of the clustering and luminosity function of galaxies at low and high redshifts, and ii) the GALFORM semi-analytic model of galaxy formation, which follows the evolution of the galaxy population from first principles, adjusted to match the statistics of low and high redshift galaxies. In both models we constrain the MgII halo sizes of field and cluster galaxies using observational results on the observed MgII statistics. Our results for the field are in good agreement with previous works, indicating a typical \mgii\ halo size of $r_MgII ~ 50h_71^-1kpc in the semi-analytic model, and slightly lower in the halo occupation number approach. For the cluster environment, we find that both models require a median MgII halo size of r_MgII< 10h_71^-1kpc in order to reproduce the observed statistics on absorption line systems in clusters of galaxies. Based on the Chen & Tinker (2008) result that stronger systems occur closer to the MgII halo centre, we find that strong absorption systems in clusters of galaxies occur at roughly a fixed fraction of the cold-warm halo size out to 1h_71^-1Mpc from the cluster centres. In contrast, weaker absorption systems appear to occur at progressively shorter relative fractions of this halo as the distance to the cluster centre decreases. △ Less

Submitted 11 February, 2009; originally announced February 2009.

Comments: 12 pages, 8 figures, accepted for publication in MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.395:1135,2009

Cosmic queuing: galaxy satellites, building blocks and the hierarchical clustering paradigm

Authors: Claudia del P. Lagos, Nelson D. Padilla, Sofia A. Cora

Abstract: We study the properties of building blocks (BBs, i.e. accreted satellites) and surviving satellites of present-day galaxies using the SAG semi-analytic model of galaxy formation in the context of a concordance Lambda Cold Dark Matter (LCDM) cosmology. We consider large numbers of DM halo merger trees spanning a wide range of masses (~1x10^10 - 2.14x10^15 Msun). We find higher metallicities for B… ▽ More We study the properties of building blocks (BBs, i.e. accreted satellites) and surviving satellites of present-day galaxies using the SAG semi-analytic model of galaxy formation in the context of a concordance Lambda Cold Dark Matter (LCDM) cosmology. We consider large numbers of DM halo merger trees spanning a wide range of masses (~1x10^10 - 2.14x10^15 Msun). We find higher metallicities for BBs with respect to surviving satellites, an effect produced by the same processes behind the build-up of the mass-metallicity relation. We prove that these metallicity differences arise from the higher peak height in the density fluctuation field occupied by BBs and central galaxies which have collapsed into a single object earlier than surviving satellites. BBs start to form stars earlier, during the peak of the merger activity in LCDM, and build-up half of their final stellar mass (measured at the moment of disruption) up to four times faster than surviving satellites. Surviving satellites keep increasing their stellar masses rather quiescently down to z~1. The difference between the metallicities of satellites, BBs and central galaxies depends on the host DM halo mass, in a way that can be used as a further test for the concordance cosmology. △ Less

Submitted 23 April, 2009; v1 submitted 13 January, 2009; originally announced January 2009.

Comments: 5 pages, 4 figures. Accepted for publication in MNRAS Letters

Black hole spin and radio loudness in a LCDM universe

Authors: Claudia del P. Lagos, Nelson D. Padilla, Sofia A. Cora

Abstract: We use a combination of a cosmological N-body simulation of the concordance Lambda cold dark matter (LCDM) paradigm and a semi-analytic model of galaxy formation to investigate the spin development of central supermassive black holes (BHs) and its relation to the BH host galaxy properties. In order to compute BH spins, we use the alpha-model of Shakura & Sunyaev and consider the King et al. warp… ▽ More We use a combination of a cosmological N-body simulation of the concordance Lambda cold dark matter (LCDM) paradigm and a semi-analytic model of galaxy formation to investigate the spin development of central supermassive black holes (BHs) and its relation to the BH host galaxy properties. In order to compute BH spins, we use the alpha-model of Shakura & Sunyaev and consider the King et al. warped disc alignment criterion. The orientation of the accretion disc is inferred from the angular momentum of the source of accreted material, which bears a close relationship to the large-scale structure in the simulation. We find that the final BH spin depends almost exclusively on the accretion history and only weakly on the warped disc alignment. The main mechanisms of BH spin-up are found to be gas cooling processes and disc instabilities, a result that is only partially compatible with Monte-Carlo models where the main spin-up mechanisms are major mergers and disc instabilities; the latter results are reproduced when implementing randomly oriented accretion discs in our model. Regarding the BH population, we find that more massive BHs, which are hosted by massive ellipticals, have higher spin values than less-massive BHs, hosted by spiral galaxies. We analyse whether gas accretion rates and BH spins can be used as tracers of the radio loudness of active galactic nuclei (AGN). We find that the current observational indications of an increasing trend of radio-loud AGN fractions with stellar and BH mass can be easily obtained when placing lower limits on the BH spin, with a minimum influence from limits on the accretion rates; a model with random accretion disc orientations is unable to reproduce this trend. (ABRIDGED) △ Less

Submitted 5 January, 2009; originally announced January 2009.

Comments: 13 pages, 8 figures. Accepted for publication in MNRAS

Galaxy star formation in different environments

Authors: R. E. Gonzalez, N. D. Padilla

Abstract: We use a semi-analytic model of galaxy formation to study signatures of large-scale modulations in the star formation (SF) activity in galaxies. In order to do this we carefully define local and global estimators of the density around galaxies. The former are computed using a voronoi tessellation technique and the latter are parameterised by the normalised distance to haloes and voids, in terms… ▽ More We use a semi-analytic model of galaxy formation to study signatures of large-scale modulations in the star formation (SF) activity in galaxies. In order to do this we carefully define local and global estimators of the density around galaxies. The former are computed using a voronoi tessellation technique and the latter are parameterised by the normalised distance to haloes and voids, in terms of the virial and void radii, respectively. As a function of local density, galaxies show a strong modulation in their SF, a result that is in agreement with those from several authors. When taking subsamples of equal local density at different large-scale environments, we find relevant global effects whereby the fraction of red galaxies diminishes for galaxies in equal local density environments farther away from clusters and closer to voids. In general, the semianalytic simulation is in good agreement with the available observational results, and offers the possibility to disentangle many of the processes responsible for the variation of galaxy properties with the environment; we find that the changes found in samples of galaxies with equal local environment but different distances to haloes or voids come from the variations in the underlying mass function of dark-matter haloes. There is an additional possible effect coming from the host dark-matter halo ages, indicating that halo assembly also plays a small but significant role ($1.14σ$) in shaping the properties of galaxies, and in particular, hints at a possible spatial correlation in halo/stellar mass ages. △ Less

Submitted 30 April, 2009; v1 submitted 8 November, 2008; originally announced November 2008.

Comments: 14 pages, 11 figures, MNRAS Accepted

On the accuracy of the high redshift cluster luminosity function

Authors: Roberto P. Muñoz, Nelson D. Padilla, Luis F. Barrientos

Abstract: We study the reliability of the statistical background subtraction method for computing the Ks-band luminosity function of cluster galaxies at z~1 using mock Red-sequence Cluster Survey cluster catalogues constructed from GALFORM semi-analytic galaxies. The underlying cluster luminosity function in the mocks are compatible with recent estimates at z~1 by several authors. We simulate different sa… ▽ More We study the reliability of the statistical background subtraction method for computing the Ks-band luminosity function of cluster galaxies at z~1 using mock Red-sequence Cluster Survey cluster catalogues constructed from GALFORM semi-analytic galaxies. The underlying cluster luminosity function in the mocks are compatible with recent estimates at z~1 by several authors. We simulate different samples where the number of clusters with Ks-band photometry goes from 5 to a maximum of 50, in order to find the most suitable observational sample to carry out this study; the current observational status in the nIR wavelength range has been reached using 5 real clusters at z~1. We compute the composite luminosity function for several samples of galaxy clusters with masses 1.5x10^14 M_sun assuming a flux limited, complete sample of galaxies down to Ks=21.0 magnitudes. We find that the Schechter fit parameters Ks* and alpha for a sample of galaxies with no redshift information are rather poorly constrained if both parameters are allowed to vary freely; if alpha is fixed at a fiducial value, then Ks* shows significantly improved stochastic uncertainties but can be influenced by systematic deviations. We find a significantly improved accuracy in the luminosity function parameters when adding photometric redshift information for bright cluster galaxies. The dwarf-to-giant ratios inferred from the luminosity functions of red-sequence galaxies in the mock catalogue agree very well with the underlying values. Finally, we find that in order to use estimates of Ks* to study the formation redshift of cluster galaxies at z=1, the sample would need to contain 520 z~1 clusters, for an accuracy of 2 Gyr at the 68 per cent confidence level. △ Less

Submitted 8 October, 2008; originally announced October 2008.

Comments: 12 pages, 8 figures. Accepted for publication in MNRAS

Evolution of environment dependent galaxy properties in the Sloan Digital Sky Survey

Authors: Ana Laura O'Mill, Nelson Padilla, Diego Garcia Lambas

Abstract: We use photometric redshifts to analyse the effect of local environment on galaxy colours at redshifts z < 0.63 in the SDSS data release 6. We construct mock SDSS-DR6 catalogues using semi-analytic galaxies to study possible systematic effects on the characterisation of environment and colour statistics due to the uncertainty in the determination of redshifts. We use the projected galaxy density… ▽ More We use photometric redshifts to analyse the effect of local environment on galaxy colours at redshifts z < 0.63 in the SDSS data release 6. We construct mock SDSS-DR6 catalogues using semi-analytic galaxies to study possible systematic effects on the characterisation of environment and colour statistics due to the uncertainty in the determination of redshifts. We use the projected galaxy density derived from the distance to the nearest neighbours with a suitable radial velocity threshold to take into account the uncertainties in the photometric redshift estimates. Our findings indicate that the use of photometric redshifts greatly improves estimates of projected local galaxy density when galaxy spectra are not available. We find a tight relationship between spectroscopic and photometric derived densities, both in the SDSS-DR6 data (up to z=0.3) and mock catalogues (up to z=0.63). At z=0, faint galaxies show a clear increase of the red galaxy fraction as the local density increases. Bright galaxies, on the other hand, show a constant red galaxy fraction. We are able to track the evolution of this fraction to z=0.55 for galaxies brighter than M_r=-21.5 and find that the fraction of blue galaxies with respect to the total population progressively becomes higher as the redshift increases, at a rate of 15% Gyr. Also, at any given redshift, bright galaxies show a larger red population, indicating that the star-formation activity shifts towards smaller objects as the redshift decreases. △ Less

Submitted 26 June, 2008; originally announced June 2008.

Comments: 9 pages, 6 figures. Accepted by MNRAS

Effects of AGN feedback on LCDM galaxies

Authors: Claudia del P. Lagos, Sofia A. Cora, Nelson D. Padilla, .

Abstract: We study the effects of Active Galactic Nuclei (AGN) feedback on the formation and evolution of galaxies in a semi-analytic model of galaxy formation. This model is an improved version of the one described by Cora (2006), which now considers the growth of black holes (BHs) as driven by (i) gas accretion during merger-driven starbursts and mergers with other BHs, (ii) accretion during starbursts… ▽ More We study the effects of Active Galactic Nuclei (AGN) feedback on the formation and evolution of galaxies in a semi-analytic model of galaxy formation. This model is an improved version of the one described by Cora (2006), which now considers the growth of black holes (BHs) as driven by (i) gas accretion during merger-driven starbursts and mergers with other BHs, (ii) accretion during starbursts triggered by disc instabilities, and (iii) accretion of gas cooled from quasi-hydrostatic hot gas haloes. It is assumed that feedback from AGN operates in the later case. The model has been calibrated in order to reproduce observational correlations between BH mass and mass, velocity dispersion, and absolute magnitudes of the galaxy bulge. AGN feedback has a strong impact on reducing or even suppressing gas cooling, an effect that becomes important at lower redshifts. This phenomenon helps to reproduce the observed galaxy luminosity function (LF) in the optical and near IR bands at z=0, and the cosmic star formation rate and stellar mass functions over a wide redshift range (0<z<5). It also allows to have a population of massive galaxies already in place at z>1, which are mostly early-type and have older and redder stellar populations than lower mass galaxies, reproducing the observed bimodality in the galaxy colour distribution, and the morphological fractions. The evolution of the optical QSO LF is also reproduced, provided that the presence of a significant fraction of obscured QSOs is assumed. We explore the effects of AGN feedback during starbursts and new recent prescriptions for dynamical friction time-scales. (ABRIDGED) △ Less

Submitted 13 May, 2008; originally announced May 2008.

Comments: 18 pages, 15 figures. Accepted for publication in MNRAS. High resolution figures available at: http://www.astro.puc.cl/~clagos/AGNOnLCDMGals.pdf

Voids in the 2dFGRS and LCDM simulations: spatial and dynamical properties

Authors: Laura Ceccarelli, Nelson Padilla, Carlos Valotto, Diego Garcia Lambas

Abstract: We perform a statistical study on the distribution and dynamics of voids in the 2dFGRS. Our statistics are tested and calibrated using mock 2dFGRS catalogues. We analyse the dynamics around voids in mock and real 2dFGRS surveys. The void-galaxy cross-correlation redshift-space distortions show evidence of the continuing growth of voids, confirming another prediction of the hierarchical clusterin… ▽ More We perform a statistical study on the distribution and dynamics of voids in the 2dFGRS. Our statistics are tested and calibrated using mock 2dFGRS catalogues. We analyse the dynamics around voids in mock and real 2dFGRS surveys. The void-galaxy cross-correlation redshift-space distortions show evidence of the continuing growth of voids, confirming another prediction of the hierarchical clustering scenario. A non-linear outflow model can be used to provide quantitative estimates of the outflow velocities around 2dFGRS voids. These results are consistent with maximum outflows of 110km/s, 210km/s and 270km/s for voids of <r_void>=7.5, 12.5 and 17.5h^-1Mpc, assuming a galaxy bias of b=1. As an application for future surveys, our study of the mock catalogues shows that direct measurements of the expansion of voids can be obtained using peculiar velocity data. We find that it would also be possible to detect differences in the velocity dispersion of galaxies in the directions parallel and perpendicular to the void walls. △ Less

Submitted 6 May, 2008; originally announced May 2008.

Comments: 13 pages, 13 figures, published in MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.373:1440,2006

Large-scale modulation of star formation in void walls

Authors: Laura Ceccarelli, Nelson Padilla, Diego Garcia Lambas

Abstract: We perform a statistical study of the characteristics of galaxies in voids and void walls in the SDSS and 2dFGRS catalogues. We investigate dependencies of the distribution of galaxy spectral types and colours as a function of the relative position to the void centres for different luminosity and local density ranges. We find a trend towards bluer, star forming galaxies in void walls beyond the… ▽ More We perform a statistical study of the characteristics of galaxies in voids and void walls in the SDSS and 2dFGRS catalogues. We investigate dependencies of the distribution of galaxy spectral types and colours as a function of the relative position to the void centres for different luminosity and local density ranges. We find a trend towards bluer, star forming galaxies in void walls beyond the local density dependence. These results indicate that luminosity and local density do not entirely determine the distribution of galaxy properties such as colours and spectral types, and point towards a large scale modulation of star formation. We argue that this effect is due to the lower accretion and merger history of galaxies arriving at void walls from the emptier inner void regions. △ Less

Submitted 6 May, 2008; originally announced May 2008.

Comments: 5 pages, 6 figures, submitted to MRNAS Letters

The influence of halo assembly on galaxies and galaxy groups

Authors: Tatiana Zapata, Josefa Perez, Nelson Padilla, Patricia Tissera

Abstract: In this paper, we study the variations of group galaxy properties according to the assembly history in SDSS-DR6 selected groups. Using mock SDSS group catalogues, we find two suitable indicators of group formation time: i) the isolation of the group, defined as the distance to the nearest neighbor in terms of its virial radius, and ii) the concentration, measured as the group inner density calcu… ▽ More In this paper, we study the variations of group galaxy properties according to the assembly history in SDSS-DR6 selected groups. Using mock SDSS group catalogues, we find two suitable indicators of group formation time: i) the isolation of the group, defined as the distance to the nearest neighbor in terms of its virial radius, and ii) the concentration, measured as the group inner density calculated using the fifth nearest bright galaxy to the group centre. Groups within narrow ranges of mass in the mock catalogue show increasing group age with isolation and concentration. However, in the observational data the stellar age, as indicated by the spectral type, only shows a correlation with concentration. We study groups of similar mass and different assembly history, finding important differences in their galaxy population. Particularly, in high mass SDSS groups, the number of members, mass-to-light ratios, red galaxy fractions and the magnitude difference between the brightest and second brightest group galaxies, show different trends as a function of isolation and concentration, even when it is expected that the latter two quantities correlate with group age. Conversely, low mass SDSS groups appear to be less sensitive to their assembly history. The correlations detected in the SDSS are not consistent with the trends measured in the mock catalogues. However, discrepancies can be explained in terms of the disagreement found in the age-isolation trends, suggesting that the model might be overestimating the effects of environment.We discuss how the modeling of the cold gas in satellite galaxies could be responsible for this problem. These results can be used to improve our understanding of the evolution of galaxies in high-density environments. △ Less

Submitted 12 January, 2009; v1 submitted 28 April, 2008; originally announced April 2008.

Comments: 11 pages, 7 figures, accepted for publication in MNRAS. Send comments to taz@astro.iag.usp.br

Angular momentum-Large-scale structure alignments in LCDM models and the SDSS

Authors: Dante Paz, Federico Stasyszyn, Nelson Padilla

Abstract: We study the alignments between the angular momentum of individual objects and the large-scale structure in cosmological numerical simulations and real data from the Sloan Digital Sky Survey, Data Release 6. To this end we measure anisotropies in the two point cross-correlation function around simulated halos and observed galaxies, studying separately the 1- and 2-halo regimes. The alignment of… ▽ More We study the alignments between the angular momentum of individual objects and the large-scale structure in cosmological numerical simulations and real data from the Sloan Digital Sky Survey, Data Release 6. To this end we measure anisotropies in the two point cross-correlation function around simulated halos and observed galaxies, studying separately the 1- and 2-halo regimes. The alignment of the angular momentum of dark-matter haloes in LCDM simulations is found to be dependent on scale and halo mass. At large distances (2-halo regime), the spins of high mass haloes are preferentially oriented in the direction perpendicular to the distribution of matter; lower mass systems show a weaker trend that may even reverse to show an angular momentum in the plane of the matter distribution. In the 1-halo term regime, the angular momentum is aligned in the direction perpendicular to the matter distribution; the effect is stronger than for the 1-halo term and increases for higher mass systems. On the observational side, we focus our study on galaxies in the Sloan Digital Sky Survey, Data Release 6 (SDSS-DR6) with elongated apparent shapes, and study alignments with respect to the major semi-axis. We find an excess of structure in the direction of the major semi-axis for all samples; the red sample shows the highest alignment (2.7+-0.08%) and indicates that the angular momentum of flattened spheroidals tends to be perpendicular to the large-scale structure. (Abridged) △ Less

Submitted 27 June, 2008; v1 submitted 28 April, 2008; originally announced April 2008.

Comments: 10 pages, 6 figures, accepted for publication in MNRAS; the definitive version is available at www.blackwell-synergy.com

Galaxy Clusters in the Line of Sight to Background Quasars: I. Survey Design and Incidence of MgII Absorbers at Cluster Redshifts

Authors: S. Lopez, L. F. Barrientos, P. Lira, N. Padilla, D. G. Gilbank, M. D. Gladders, J. Maza, N. Tejos, M. Vidal, H. K. C. Yee

Abstract: We describe the first optical survey of absorption systems associated with galaxy clusters at z= 0.3-0.9. We have cross-correlated SDSS DR3 quasars with high-redshift cluster/group candidates from the Red-Sequence Cluster Survey. We have found 442 quasar-cluster pairs for which the MgII doublet might be detected at a transverse (physical) distance d<2 Mpc from the cluster centers. To investigate… ▽ More We describe the first optical survey of absorption systems associated with galaxy clusters at z= 0.3-0.9. We have cross-correlated SDSS DR3 quasars with high-redshift cluster/group candidates from the Red-Sequence Cluster Survey. We have found 442 quasar-cluster pairs for which the MgII doublet might be detected at a transverse (physical) distance d<2 Mpc from the cluster centers. To investigate the incidence (dN/dz) and equivalent-width distribution n(W) of MgII systems at cluster redshifts, two statistical samples were drawn out of these pairs: one made of high-resolution spectroscopic quasar observations (46 pairs), and one made of quasars used in MgII searches found in the literature (375 pairs). The results are: (1) the population of strong MgII systems (W_0>2.0 Ang.) near cluster redshifts shows a significant (>3 sigma) overabundance (up to a factor of 15) when compared with the 'field' population; (2) the overabundance is more evident at smaller distances (d<1 Mpc) than larger distances (d<2 Mpc) from the cluster center; and, (3) the population of weak MgII systems (W_0<0.3 Ang.) near cluster redshifts conform to the field statistics. Unlike in the field, this dichotomy makes n(W) in clusters appear flat and well fitted by a power-law in the entire W-range. A sub-sample of the most massive clusters yields a stronger and still significant signal. Since either the absorber number density or filling-factor/cross-section affects the absorber statistics, an interesting possibility is that we have detected the signature of truncated halos due to environmental effects. Thus, we argue that the excess of strong systems is due to a population of absorbers in an overdense galaxy environment, and the lack of weak systems to a different population, that got destroyed in the cluster environment. (Abridged) △ Less

Submitted 5 March, 2008; originally announced March 2008.

Comments: Accepted for publication in the Astrophysical Journal

The shapes of galaxies in the Sloan Digital Sky Survey

Authors: Nelson D. Padilla, Michael A. Strauss

Abstract: We determine the underlying shapes of spiral and elliptical galaxies in the Sloan Digital Sky Survey Data Release 6 from the observed distribution of projected galaxy shapes, taking into account the effects of dust extinction and reddening. We assume that the underlying shapes of spirals and ellipticals are well approximated by triaxial ellipsoids. The elliptical galaxy data are consistent with… ▽ More We determine the underlying shapes of spiral and elliptical galaxies in the Sloan Digital Sky Survey Data Release 6 from the observed distribution of projected galaxy shapes, taking into account the effects of dust extinction and reddening. We assume that the underlying shapes of spirals and ellipticals are well approximated by triaxial ellipsoids. The elliptical galaxy data are consistent with oblate spheroids, with a correlation between luminosity and ellipticity: the mean values of minor to middle axis ratios are 0.41+-0.03 for Mr ~ -18 ellipticals, and 0.76+-0.04 for Mr ~-22.5 ellipticals. Ellipticals show almost no dependence of axial ratio on galaxy colour, implying a negligible dust optical depth. There is a strong variation of spiral galaxy shapes with colour indicating the presence of dust. The intrinsic shapes of spiral galaxies in the SDSS-DR6 are consistent with flat disks with a mean and dispersion of thickness to diameter ratio of (21+-2)%, and a face-on ellipticity, e, of ln(e)=-2.33+-0.79. Not including the effects of dust in the model leads to disks that are systematically rounder by up to 60%. More luminous spiral galaxies tend to have thicker and rounder disks than lower-luminosity spirals. Both elliptical and spiral galaxies tend to be rounder for larger galaxies. The marginalised value of the edge-on r-band dust extinction E_0 in spiral galaxies is E_0 ~ 0.45 magnitudes for galaxies of median colours, increasing to E_0=1 magnitudes for g-r>0.9 and E_0=1.9 for the luminous and most compact galaxies, with half-light radii <2kpc/h. △ Less

Submitted 19 May, 2008; v1 submitted 6 February, 2008; originally announced February 2008.

Comments: 18 pages, 14 figures, accepted for publication in MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.388:1321-1334,2008

Ly Alpha-Emitting Galaxies at z=3.1: L* Progenitors Experiencing Rapid Star Formation

Authors: Eric Gawiser, Harold Francke, Kamson Lai, Kevin Schawinski, Caryl Gronwall, Robin Ciardullo, Ryan Quadri, Alvaro Orsi, L. Felipe Barrientos, Guillermo A. Blanc, Giovanni Fazio, John J. Feldmeier, Jia-Sheng Huang, Leopoldo Infante, Paulina Lira, Nelson Padilla, Edward N. Taylor, Ezequiel Treister, C. Megan Urry, Pieter G. van Dokkum, Shanil N. Virani

Abstract: We studied the clustering properties and multiwavelength spectral energy distributions of a complete sample of 162 Ly Alpha-Emitting (LAE) galaxies at z=3.1 discovered in deep narrow-band MUSYC imaging of the Extended Chandra Deep Field South. LAEs were selected to have observed frame equivalent widths >80A and emission line fluxes >1.5E-17 erg/cm^2/s. Only 1% of our LAE sample appears to host A… ▽ More We studied the clustering properties and multiwavelength spectral energy distributions of a complete sample of 162 Ly Alpha-Emitting (LAE) galaxies at z=3.1 discovered in deep narrow-band MUSYC imaging of the Extended Chandra Deep Field South. LAEs were selected to have observed frame equivalent widths >80A and emission line fluxes >1.5E-17 erg/cm^2/s. Only 1% of our LAE sample appears to host AGN. The LAEs exhibit a moderate spatial correlation length of r_0=3.6+0.8-1.0 Mpc, corresponding to a bias factor b=1.7+0.3-0.4, which implies median dark matter halo masses of log10(M_med) = 10.9+0.5-0.9 M_sun. Comparing the number density of LAEs, (1.5+-0.3)E-3/Mpc^3, with the number density of these halos finds a mean halo occupation ~1-10%. The evolution of galaxy bias with redshift implies that most z=3.1 LAEs evolve into present-day galaxies with L<2.5L*, whereas other z>3 galaxy populations typically evolve into more massive galaxies. Halo merger trees show that z=0 descendants occupy halos with a wide range of masses, with a median descendant mass close to that of L*. Only 30% of LAEs have sufficient stellar mass (>~3E9 M_sun) to yield detections in deep Spitzer-IRAC imaging. A two-population SED fit to the stacked UBVRIzJK+[3.6,4.5,5.6,8.0]micron fluxes of the IRAC-undetected objects finds that the typical LAE has low stellar mass (1.0+0.6-0.4 E9 M_sun), moderate star formation rate (2+-1 M_sun/yr), a young component age of 20+30-10 Myr, and little dust (A_V<0.2). The best fit model has 20% of the mass in the young stellar component, but models without evolved stars are also allowed. △ Less

Submitted 14 October, 2007; originally announced October 2007.

Comments: ApJ, in press, 7 pages including 4 color figures

Journal ref: Astrophysical Journal 671, 278 (2007)

Spatial clustering of USS sources and galaxies

Authors: Carlos G. Bornancini, Nelson Padilla, Diego G. Lambas, Carlos De Breuck

Abstract: We present measurements of the clustering properties of galaxies in the field of redshift range $0.5 \lesssim z \lesssim 1.5$ Ultra Steep Spectrum (USS) radio sources selected from Sydney University Molonglo Sky Survey and NRAO VLA Sky Survey. Galaxies in these USS fields were identified in deep near-IR observations, complete down to $K_s=20$, using {\tt IRIS2} instrument at the AAT telescope. W… ▽ More We present measurements of the clustering properties of galaxies in the field of redshift range $0.5 \lesssim z \lesssim 1.5$ Ultra Steep Spectrum (USS) radio sources selected from Sydney University Molonglo Sky Survey and NRAO VLA Sky Survey. Galaxies in these USS fields were identified in deep near-IR observations, complete down to $K_s=20$, using {\tt IRIS2} instrument at the AAT telescope. We used the redshift distribution of $K_{s} < 20$ galaxies taken from Cimatti et al. (2002) to constrain the correlation length $r_0$. We find a strong correlation signal of galaxies with $K_{s} < 20$ around our USS sample. A comoving correlation length $r_{0}=14.0\pm2.8$ $h^{-1}$ Mpc and $γ=1.98\pm0.15$ are derived in a flat cosmological model Universe. We compare our findings with those obtained in a cosmological N--body simulation populated with GALFORM semi-analytic galaxies. We find that clusters of galaxies with masses in the range $M=10^{13.4-14.2}$ $h^{-1}$ M$_{\sun}$ have a cluster--galaxy cross--correlation amplitude comparable to those found between USS hosts and galaxies. These results suggest that distant radio galaxies are excellent tracers of galaxy overdensities and pinpoint the progenitors of present day rich clusters of galaxies. △ Less

Submitted 1 February, 2006; originally announced February 2006.

Comments: Accepted for publication in MNRAS, 5 pages including 3 PostScript figures and 1 table

Journal ref: Mon.Not.Roy.Astron.Soc.368:619-622,2006

Cluster infall in the concordance LCDM model

Authors: Maximiliano C. Pivato, Nelson D. Padilla, Diego G. Lambas

Abstract: We perform statistical analyses of the infall of dark-matter onto clusters in numerical simulations within the concordance LCDM model. By studying the infall profile around clusters of different mass, we find a linear relation between the maximum infall velocity and mass which reach 900km/s for the most massive groups. The maximum infall velocity and the group mass follow a suitable power law fi… ▽ More We perform statistical analyses of the infall of dark-matter onto clusters in numerical simulations within the concordance LCDM model. By studying the infall profile around clusters of different mass, we find a linear relation between the maximum infall velocity and mass which reach 900km/s for the most massive groups. The maximum infall velocity and the group mass follow a suitable power law fit of the form, V_{inf}^{max} = (M/m_0)^{gamma}. By comparing the measured infall velocity to the linear infall model with an exponential cutoff introduced by Croft et al., we find that the best agreement is obtained for a critical overdensity delta_c = 45. We study the dependence of the direction of infall with respect to the cluster centres, and find that in the case of massive groups, the maximum alignment occurs at scales r ~ 6Mpc/h. We obtain a logarithmic power-law relation between the average infall angle and the group mass. We also study the dependence of the results on the local dark-matter density, finding a remarkable difference in the dynamical behaviour of low- and high-density particles. △ Less

Submitted 29 August, 2006; v1 submitted 6 December, 2005; originally announced December 2005.

Comments: 8 pages, 10 figures. Accepted to publication in MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.373:1409-1415,2006

Galaxy Peculiar Velocities and Infall onto Groups

Authors: M. L. Ceccarelli, C. Valotto, D. G. Lambas, N. Padilla, R. Giovanelli, M. Haynes

Abstract: We perform statistical analyses to study the infall of galaxies onto groups and clusters in the nearby Universe. The study is based on the UZC and SSRS2 group catalogs and peculiar velocity samples. We find a clear signature of infall of galaxies onto groups over a wide range of scales 5 h^{-1} Mpc<r<30 h^{-1} Mpc, with an infall amplitude on the order of a few hundred kilometers per second. We… ▽ More We perform statistical analyses to study the infall of galaxies onto groups and clusters in the nearby Universe. The study is based on the UZC and SSRS2 group catalogs and peculiar velocity samples. We find a clear signature of infall of galaxies onto groups over a wide range of scales 5 h^{-1} Mpc<r<30 h^{-1} Mpc, with an infall amplitude on the order of a few hundred kilometers per second. We obtain a significant increase in the infall amplitude with group virial mass (M_{V}) and luminosity of group member galaxies (L_{g}). Groups with M_{V}<10^{13} M_{\odot} show infall velocities V_{infall} \simeq 150 km s^{-1} whereas for M_{V}>10^{13} M_{\odot} a larger infall is observed, V_{infall} \simeq 200 km s^{-1}. Similarly, we find that galaxies surrounding groups with L_{g}<10^{15} L_{\odot} have V_{infall} \simeq 100 km s^{-1}, whereas for L_{g}>10^{15} L_{\odot} groups, the amplitude of the galaxy infall can be as large as V_{infall} \simeq 250 km s^{-1}. The observational results are compared with the results obtained from mock group and galaxy samples constructed from numerical simulations, which include galaxy formation through semianalytical models. We obtain a general agreement between the results from the mock catalogs and the observations. The infall of galaxies onto groups is suitably reproduced in the simulations and, as in the observations, larger virial mass and luminosity groups exhibit the largest galaxy infall amplitudes. We derive estimates of the integrated mass overdensities associated with groups by applying linear theory to the infall velocities after correcting for the effects of distance uncertainties obtained using the mock catalogs. The resulting overdensities are consistent with a power law with δ\sim 1 at r \sim 10 h^{-1}Mpc. △ Less

Submitted 16 September, 2005; originally announced September 2005.

Comments: 25 pages, 10 figures

Journal ref: Astrophys.J.622:853-861,2005

Shapes of clusters and groups of galaxies: Comparison of model predictions with observations

Authors: D. J. Paz, D. G. Lambas, N. Padilla, M. Merchan

Abstract: We study the properties of the 3-dimensional and projected shapes of haloes using high resolution numerical simulations and observational data where the latter comes from the 2PIGG (Eke et al. 2004) and SDSS-DR3GC group catalogues (Merchan & Zandivarez 2005). We investigate the dependence of halo shape on characteristics such as mass and number of members. In the 3-dimensional case, we find a si… ▽ More We study the properties of the 3-dimensional and projected shapes of haloes using high resolution numerical simulations and observational data where the latter comes from the 2PIGG (Eke et al. 2004) and SDSS-DR3GC group catalogues (Merchan & Zandivarez 2005). We investigate the dependence of halo shape on characteristics such as mass and number of members. In the 3-dimensional case, we find a significant correlation between the mass and halo shape; massive systems are more prolate than small haloes. We detect a source of strong systematics in estimates of the triaxiality of a halo, which is found to be a strong function of the number of members; LCDM haloes usually characterised by triaxial shapes, slightly bent toward prolate forms, appear more oblate when taking only a small subset of the halo particles. The ellipticities of observed 2PIGG and SDSS-DR3GC groups are found to be strongly dependent on the number of group members, so that poor groups appear more elongated than rich ones. However, this is again an artifact caused by poor statistics and not an intrinsic property of the galaxy groups, nor an effect from observational biases. We interpret these results with the aid of a GALFORM mock 2PIGG catalogue. When comparing the group ellipticities in mock and real catalogues, we find an excellent agreement between the trends of shapes with number of group members. When carefully taking into account the effects of low number statistics, we find that more massive groups are consistent with more elongated shapes. Finally, our studies find no significant correlations between the shape of observed 2PIGG or SDSS-DR3GC groups with the properties of galaxy members such as colour or spectral type index. △ Less

Submitted 2 September, 2005; originally announced September 2005.

Comments: 9 pages, 10 figures, submitted to MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.366:1503-1510,2006

Spatial and Dynamical Properties of Voids in a LCDM Universe

Authors: N. D. Padilla, L. Ceccarelli, D. G. Lambas

Abstract: We study statistical properties of voids in the distribution of mass, dark-matter haloes and galaxies (B_J<-16) in a LCDM numerical simulation populated with galaxies using a semi-analytic galaxy formation model(GALFORM, Cole et al. 2000). We find that the properties of voids selected from GALFORM galaxies are compatible with those of voids identified from a population of haloes with mass M>10^1… ▽ More We study statistical properties of voids in the distribution of mass, dark-matter haloes and galaxies (B_J<-16) in a LCDM numerical simulation populated with galaxies using a semi-analytic galaxy formation model(GALFORM, Cole et al. 2000). We find that the properties of voids selected from GALFORM galaxies are compatible with those of voids identified from a population of haloes with mass M>10^11.5 M_sun/h, similar to the median halo mass, M_med=10^11.3 M_sun/h. We also find that the number density of galaxy- and halo-defined voids can be up to two orders of magnitude higher than mass-defined voids for large void radii. As expected, there are outflow velocities which show their maximum at larger void-centric distances for larger voids (well described by a linear relation). The void-centric distance where this maximum occurs, follows a suitable power law fit of the form, log(d_vmax)=(r_void/A)^B. At sufficiently large distances, we find mild infall motions onto the sub-dense regions. We find that a similar analysis in redshift space would make both outflows and infalls to appear with a lower amplitude. We also find that the velocity dispersion of galaxies and haloes is larger in the direction parallel to the void walls by ~10-20%. We analyse the void-tracer cross-correlation functions in real- and redshift-space as a function of separation, and also as a function of separations parallel and perpendicular to the line of sight. The distortion pattern observed in xi(sigma,pi) is that of an elongation along the line of sight which extends out to large separations. Positive xi contours evidence finger-of-god motions at the void walls. Elongations along the line of sight are roughly comparable between galaxy-, halo- and mass-defined voids. (Abridged) △ Less

Submitted 12 August, 2005; originally announced August 2005.

Comments: 16 pages, 11 figures, accepted for Publication in MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.363:977-990,2005

The galaxy luminosity function around groups

Authors: R. E. Gonzalez, N. D. Padilla, G. Galaz, L. Infante

Abstract: We present a study on the variations of the luminosity function of galaxies around clusters in a numerical simulation with semi-analytic galaxies, attempting to detect these variations in the 2dF Galaxy Redshift Survey. We subdivide the simulation box in equal-density regions around clusters, which we assume can be achieved by selecting objects at a given normalised distance ($r/r_{rms}$, where… ▽ More We present a study on the variations of the luminosity function of galaxies around clusters in a numerical simulation with semi-analytic galaxies, attempting to detect these variations in the 2dF Galaxy Redshift Survey. We subdivide the simulation box in equal-density regions around clusters, which we assume can be achieved by selecting objects at a given normalised distance ($r/r_{rms}$, where $r_{rms}$ is an estimate of the halo radius) from the group centre. The semi-analytic model predicts important variations in the luminosity function out to $r/r_{rms}\simeq5$. In brief, variations in the mass function of haloes around clusters (large dark-matter haloes with $M>10^{12}$h$^{-1}$M$_{\sun}$) lead to cluster central regions that present a high abundance of bright galaxies (high $M^*$ values) as well as low luminosity galaxies (high $α$); at $r/r_{rms}\simeq 3$ there is a lack of bright galaxies, which shows the depletion of galaxies in the regions surrounding clusters (minimum in $M^*$ and $α$), and a tendency to constant luminosity function parameters at larger cluster-centric distances. We take into account the observational biases present in the real data by reproducing the peculiar velocity effect on the redshifts of galaxies in the simulation box, and also by producing mock catalogues. We find that excluding from the analysis galaxies wich in projection are close to the centres of the groups provides results that are qualitatively consistent with the full simulation box results. △ Less

Submitted 9 August, 2005; originally announced August 2005.

Comments: 11 pages, 10 figures, MNRAS Accepted

Journal ref: Mon.Not.Roy.Astron.Soc.363:1008-1018,2005

Cosmological parameters from CMB measurements and the final 2dFGRS power spectrum

Authors: Ariel G. Sanchez, C. M. Baugh, W. J. Percival, J. A. Peacock, N. D. Padilla, S. Cole, C. S. Frenk, P. Norberg

Abstract: We derive constraints on cosmological parameters using the power spectrum of galaxy clustering measured from the final two-degree field galaxy redshift survey (2dFGRS) and a compilation of measurements of the temperature power spectrum and temperature-polarization cross-correlation of the cosmic microwave background radiation. We analyse a range of parameter sets and priors, allowing for massive… ▽ More We derive constraints on cosmological parameters using the power spectrum of galaxy clustering measured from the final two-degree field galaxy redshift survey (2dFGRS) and a compilation of measurements of the temperature power spectrum and temperature-polarization cross-correlation of the cosmic microwave background radiation. We analyse a range of parameter sets and priors, allowing for massive neutrinos, curvature, tensors and general dark energy models. In all cases, the combination of datasets tightens the constraints, with the most dramatic improvements found for the density of dark matter and the energy-density of dark energy. If we assume a flat universe, we find a matter density parameter of $Ω_{\rm m}=0.237 \pm 0.020$, a baryon density parameter of $Ω_{\rm b} = 0.041 \pm 0.002$, a Hubble constant of $H_{0}=74\pm2 {\rm kms}^{-1}{\rm Mpc}^{-1}$, a linear theory matter fluctuation amplitude of $σ_{8}=0.77\pm0.05$ and a scalar spectral index of $n_{\rm s}=0.954 \pm 0.023$ (all errors show the 68% interval). Our estimate of $n_{\rm s}$ is only marginally consistent with the scale invariant value $n_{\rm s}=1$; this spectrum is formally excluded at the 95% confidence level. However, the detection of a tilt in the spectrum is sensitive to the choice of parameter space. If we allow the equation of state of the dark energy to float, we find $w_{\rm DE}= -0.85_{-0.17}^{+0.18}$, consistent with a cosmological constant. We also place new limits on the mass fraction of massive neutrinos: $f_ν < 0.105$ at the 95% level, corresponding to $\sum m_ν < 1.2$ eV. △ Less

Submitted 24 March, 2006; v1 submitted 26 July, 2005; originally announced July 2005.

Comments: MNRAS, in press. Minor revision after referee's report. 22 pages, 18 plots. Colour figures for talks (including additional plots which do not appear in the paper) can be downloaded from http://star-www.dur.ac.uk/~cmb/Params/

Journal ref: Mon.Not.Roy.Astron.Soc.366:189-207,2006

The 2dF Galaxy Redshift Survey: The clustering of galaxy groups

Authors: Nelson D. Padilla, C. M. Baugh, V. R. Eke, P. Norberg, S. Cole, C. S. Frenk, D. J. Croton, I. K. Baldry, J. Bland-Hawthorn, T. Bridges, R. Cannon, M. Colless, C. Collins, W. Couch, G. Dalton, R. De Propris, S. P. Driver, G. Efstathiou, R. S. Ellis, K. Glazebrook, C. Jackson, O. Lahav, I. Lewis, S. Lumsden, S. Maddox , et al. (5 additional authors not shown)

Abstract: We measure the clustering of galaxy groups in the 2dFGRS Percolation-Inferred Galaxy Group (2PIGG) catalogue. The 2PIGG sample has 29,000 groups with at least two members. The clustering amplitude of the full 2PIGG catalogue is weaker than that of 2dFGRS galaxies, in agreement with theoretical predictions. We have subdivided the 2PIGG catalogue into samples that span a factor of 25 in median tot… ▽ More We measure the clustering of galaxy groups in the 2dFGRS Percolation-Inferred Galaxy Group (2PIGG) catalogue. The 2PIGG sample has 29,000 groups with at least two members. The clustering amplitude of the full 2PIGG catalogue is weaker than that of 2dFGRS galaxies, in agreement with theoretical predictions. We have subdivided the 2PIGG catalogue into samples that span a factor of 25 in median total luminosity. Our correlation function measurements span an unprecedented range of clustering strengths, connecting the regimes probed by groups fainter than L* galaxies and rich clusters. There is a steady increase in clustering strength with group luminosity; the most luminous groups are ten times more strongly clustered than the full 2PIGG catalogue. We demonstrate that the 2PIGG results are in very good agreement with the clustering of groups expected in the LCDM model. △ Less

Submitted 12 July, 2004; v1 submitted 25 February, 2004; originally announced February 2004.

Comments: 16 pages. Submitted to MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.352:211,2004

Galaxy groups in the 2dF galaxy redshift survey: Large Scale Structure with Groups

Authors: Ariel Zandivarez, Manuel E. Merchan, Nelson D. Padilla

Abstract: We use the 2dF Galaxy Group Catalogue constructed by Merchán & Zandivarez to study the large scale structure of the Universe traced by galaxy groups. We concentrate on the computation of the power spectrum and the two point correlation function. The resulting group power spectrum shows a similar shape to the galaxy power spectrum obtained from the 2dF Galaxy Redshift Survey by Percival et al., b… ▽ More We use the 2dF Galaxy Group Catalogue constructed by Merchán & Zandivarez to study the large scale structure of the Universe traced by galaxy groups. We concentrate on the computation of the power spectrum and the two point correlation function. The resulting group power spectrum shows a similar shape to the galaxy power spectrum obtained from the 2dF Galaxy Redshift Survey by Percival et al., but with a higher amplitude quantified by a relative bias in redshift space of $b_s(k) \sim 1.5$ . The group two point correlation function for the total sample is well described by a power law with correlation length $s_0=8.9 \pm 0.3 \mpc$ and slope $γ=-1.6 \pm 0.1$ on scales $s < 20 \mpc$. In order to study the dependence of the clustering properties on group mass we split the catalogue in four subsamples defined by different ranges of group virial masses. These computations allow a fair estimate of the relation described by the correlation length $s_0$ and the mean intergroup separation $d_c$ for galaxy systems of low mass. We also extend our study to the redshift space distortions of galaxy groups, where we find that the anisotropies in the clustering pattern of the 2dF group catalogue are consistent with gravitational instability, with a flattening of the redshift-space correlation function contours in the direction of the line of sight. △ Less

Submitted 19 March, 2003; originally announced March 2003.

Comments: 11 pages, 9 figures, resubmitted to MNRAS after revision

Journal ref: Mon.Not.Roy.Astron.Soc. 344 (2003) 247

Anisotropies in the redshift-space correlations of galaxy groups and clusters II: Analysis of observational cluster samples

Authors: N. D. Padilla, D. G. Lambas

Abstract: We study the redshift space correlation function of galaxy clusters for observational samples constructed in different surveys. We explore correlation amplitudes, pairwise velocity distributions and bias factors. Systematics in cluster identification procedures are the main source of biased estimates of the correlation amplitude and inferred velocity dispersions. By comparing our results with th… ▽ More We study the redshift space correlation function of galaxy clusters for observational samples constructed in different surveys. We explore correlation amplitudes, pairwise velocity distributions and bias factors. Systematics in cluster identification procedures are the main source of biased estimates of the correlation amplitude and inferred velocity dispersions. By comparing our results with those of numerical simulations we explore the strong influence on the clustering distortion pattern in redshift space from effects associated with the cluster identification procedure from two dimensional surveys as is the case of Abell and APM clusters. The identification of clusters in X-ray surveys improves this situation, although there are still systematic effects which are probably due to identification of optical sources in the determination of cluster redshifts. These systematics are particularly strong for the most luminous X-ray selected clusters in the Extended Bright Cluster Survey, which exhibits very large anisotropies, comparable to those present in the Abell catalogue. Our results demonstrate that forthcoming large redshift surveys will be extremely important for the construction of new samples of groups and clusters as well as improving the determination of optical and X-ray selected cluster distances, essential for reliable analyses of the large scale structure. △ Less

Submitted 16 February, 2003; originally announced February 2003.

Comments: 12 pages, 10 figures. Revised version re-submitted to MNRAS. Info and figures can be found at http://star-www.dur.ac.uk/~nelsonp/anisotropies

Anisotropies in the redshift-space correlations of galaxy groups and clusters I: Simulated catalogues

Authors: N. D. Padilla, D. G. Lambas

Abstract: We analyse the correlation function of mock galaxy clusters in redshift space. We constructed several mock catalogues designed to mimic the selection biases inherent in a variety of observational surveys. We explore different effects that contribute to the distortion of the clustering pattern; the pairwise velocity distribution of galaxy systems, coherent bulk motions, redshift errors and system… ▽ More We analyse the correlation function of mock galaxy clusters in redshift space. We constructed several mock catalogues designed to mimic the selection biases inherent in a variety of observational surveys. We explore different effects that contribute to the distortion of the clustering pattern; the pairwise velocity distribution of galaxy systems, coherent bulk motions, redshift errors and systematics in cluster identification. Our tests show that the redshift-space clustering pattern of galaxy systems is highly influenced by effects associated with the identification procedure from two dimensional surveys. These systems show a spuriously large correlation amplitude, an effect that is present and even stronger in a subsample whose angular positions coincide with 3-dimensional identified clusters. The use of a small number of redshift measurements in the determination of the cluster distance has a smaller impact on the measured correlation function. We show that an angular incompleteness in redshift surveys, such as that present in the 2dFGRS 100k public release, has no significant effect in the results. We suggest that the nature of projection effects arise mainly from structures along the line of sight in the filamentary large-scale clustering pattern. Thus spectroscopic surveys are the only means of providing unbiased cluster samples. △ Less

Submitted 16 February, 2003; originally announced February 2003.

Comments: 14 pages, 13 figures. Revised version re-submitted to MNRAS. Info and high resolution figures can be found at http://star-www.dur.ac.uk/~nelsonp/anisotropies

Journal ref: Mon.Not.Roy.Astron.Soc. 342 (2003) 519

The power spectrum of galaxy clustering in the APM survey

Authors: N. D. Padilla, C. M. Baugh

Abstract: We measure the power spectrum of galaxy clustering in real space from the APM Galaxy Survey. We present an improved technique for the numerical inversion of Limber's equation that relates the angular clustering of galaxies to an integral over the power spectrum in three dimensions. Our approach is underpinned by a large ensemble of mock galaxy catalogues constructed from the Hubble Volume N-body… ▽ More We measure the power spectrum of galaxy clustering in real space from the APM Galaxy Survey. We present an improved technique for the numerical inversion of Limber's equation that relates the angular clustering of galaxies to an integral over the power spectrum in three dimensions. Our approach is underpinned by a large ensemble of mock galaxy catalogues constructed from the Hubble Volume N-body simulations. The mock catalogues are used to test for systematic effects in the inversion algorithm and to estimate the errors on our measurement. We find that we can recover the power spectrum to an accuracy of better than 15% over three decades in wavenumber. A key advantage of the use of mock catalogues to infer errors is that we can apply our technique on scales for which the density fluctuations are not Gaussian, thus probing the regime that offers the best constraints on models of galaxy formation. On large scales, our measurement of the power spectrum is consistent with the shape of the mass power spectrum in the popular ``concordance'' cold dark matter model. The galaxy power spectrum on small scales is strongly affected by nonlinear evolution of density fluctuations, and, to a lesser degree, by galaxy bias. The rms variance in the galaxy distribution, when smoothed in spheres of radius 8 h^(-1)Mpc, is sigma^g_8=0.96^(+0.17)_(-0.20) and the shape of the power spectrum on large scales is described by a simple fitting formula with parameter Gamma= 0.19^(+0.13)_(-0.04) (these errors are the 1-sigma ranges for a two parameter fit). We use our measurement of the power spectrum to estimate the galaxy two point correlation function; the results are well described by a power law with correlation length r_0 = 5.9 +- 0.7 h^(-1)Mpc and slope gamma = 1.61 +- 0.06 for pair separations in the range 0.1 < r/[h^(-1)Mpc] < 20. △ Less

Submitted 28 July, 2003; v1 submitted 6 January, 2003; originally announced January 2003.

Comments: 20 pages, 19 postscript figures, final version accepted for publication in MNRAS. Info and figures at http://star-www.dur.ac.uk/~nelsonp/apm/APM.html

Journal ref: Mon.Not.Roy.Astron.Soc. 343 (2003) 796

Determination of the Linear Mass Power Spectrum From the Mass Function of Galaxy Clusters

Authors: A. G. Sanchez, N. D. Padilla, D. G. Lambas

Abstract: We develop a new method to determine the linear mass power spectrum using the mass function of galaxy clusters. We obtain the rms mass fluctuation sigma (M) using the expression for the mass function in the Press & Schechter (1974), Sheth, Mo & Tormen (2001) and Jenkins et al. (2001) formalisms. We apply different techniques to recover the adimensional power spectrum Delta ^{2}(k) from sigma (M)… ▽ More We develop a new method to determine the linear mass power spectrum using the mass function of galaxy clusters. We obtain the rms mass fluctuation sigma (M) using the expression for the mass function in the Press & Schechter (1974), Sheth, Mo & Tormen (2001) and Jenkins et al. (2001) formalisms. We apply different techniques to recover the adimensional power spectrum Delta ^{2}(k) from sigma (M) namely the k_{eff} approximation, the singular value decomposition and the linear regularization method. The application of these techniques to the tCDM and LCDM GIF simulations shows a high efficiency in recovering the theoretical power spectrum over a wide range of scales. We compare our results with those derived from the power spectrum of the spatial distribution of the same sample of clusters in the simulations obtained by application of the classical Feldman, Kaiser & Peacock (1994), FKP, method. We find that the mass function based method presented here can provide a very accurate estimate of the linear power spectrum, particularly for low values of k. This estimate is comparable, or even better behaved, than the FKP solution. The principal advantage of our method is that it allows the determination of the linear mass power spectrum using the joint information of objects of a wide range of masses without dealing with specific assumptions on the bias relative to the underlying mass distribution. △ Less

Submitted 23 July, 2002; originally announced July 2002.

Comments: 12 pages, 13 poscript figures, accepted for publication in MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc. 337 (2002) 161

Cluster correlations in redshift space

Authors: Nelson D. Padilla, Carlton M. Baugh

Abstract: We test an analytic model for the two-point correlations of galaxy clusters in redshift space using the Hubble Volume N-body simulations. The correlation function of clusters shows no enhancement along the line of sight, due to the lack of any virialised structures in the cluster distribution. However, the distortion of the clustering pattern due to coherent bulk motions is clearly visible. The… ▽ More We test an analytic model for the two-point correlations of galaxy clusters in redshift space using the Hubble Volume N-body simulations. The correlation function of clusters shows no enhancement along the line of sight, due to the lack of any virialised structures in the cluster distribution. However, the distortion of the clustering pattern due to coherent bulk motions is clearly visible. The distribution of cluster peculiar motions is well described by a Gaussian, except in the extreme high velocity tails. The simulations produce a small but significant number of clusters with large peculiar motions. The form of the redshift space power spectrum is strongly influenced by errors in measured cluster redshifts in extant surveys. When these errors are taken into account, the model reproduces the power spectrum recovered from the simulation to an accuracy of 15% or better over a decade in wavenumber. We compare our analytic predictions with the power spectrum measured from t! he! APM cluster redshift survey. The cluster power spectrum constrains the amplitude of density fluctuations, as measured by the linear rms variance in spheres of radius 8 h^(-1) Mpc, denoted by sigma_8. When combined with the constraints on sigma_8 and the density parameter Omega derived from the local abundance of clusters, we find a best fitting cold dark matter model with sigma_8~1.25 and Omega~0.2, for a power spectrum shape that matches that measured for galaxies. However, for the best fitting value of Omega and given the value of Hubble's constant from recent measurements, the assumed shape of the power spectrum is incompatible with the most readily motivated predictions from the cold dark matter paradigm. △ Less

Submitted 17 December, 2001; v1 submitted 18 April, 2001; originally announced April 2001.

Comments: 2002, MNRAS, 329, 431-444

Journal ref: Mon.Not.Roy.Astron.Soc. 329 (2002) 431

Redshift-Space Distortions of Group and Galaxy Correlations

Authors: Nelson D. Padilla, Manuel E. Merchán, Carlos A. Valotto, Diego G. Lambas, Marcio A. G. Maia

Abstract: The two-point correlation function is computed for galaxies and groups of galaxies selected using 3-dimensional information from the Updated Zwicky Galaxy Catalog - (UZC). The redshift space distortion of the correlation function $ξ(σ,π)$ in the directions parallel and perpendicular to the line of sight, induced by pairwise group peculiar velocities is evaluated. Two methods are used to characte… ▽ More The two-point correlation function is computed for galaxies and groups of galaxies selected using 3-dimensional information from the Updated Zwicky Galaxy Catalog - (UZC). The redshift space distortion of the correlation function $ξ(σ,π)$ in the directions parallel and perpendicular to the line of sight, induced by pairwise group peculiar velocities is evaluated. Two methods are used to characterize the pairwise velocity field of groups and galaxies. The best fitting parameters are $γ=-1.89 \pm 0.17$ and $r_0=9.7 \pm 4.5$ $h^{-1}$ Mpc for groups, and $γ=-2.00 \pm 0.03$ and $r_0=5.29 \pm 0.21$ $h^{-1}$ Mpc for galaxies. We find $β_{gx}=0.51 \pm 0.15$ for galaxies, while for groups we obtain a noisy estimate $β< 1.5$. △ Less

Submitted 21 February, 2001; originally announced February 2001.

Peculiar Velocities and the Mean Density Parameter

Authors: Nelson Padilla, Diego G. Lambas

Abstract: We study the peculiar velocity field inferred from the Mark III spirals using a new method of analysis. We estimate optimal values of Tully-Fisher scatter and zero-point offset, and we derive the 3-dimensional rms peculiar velocity ($σ_v$) of the galaxies in the samples analysed. We check our statistical analysis using mock catalogs derived from numerical simulations of CDM models considering me… ▽ More We study the peculiar velocity field inferred from the Mark III spirals using a new method of analysis. We estimate optimal values of Tully-Fisher scatter and zero-point offset, and we derive the 3-dimensional rms peculiar velocity ($σ_v$) of the galaxies in the samples analysed. We check our statistical analysis using mock catalogs derived from numerical simulations of CDM models considering measurement uncertainties and sampling variations. Our best determination for the observations is $σ_v= (660\pm50) km/s$. We use the linear theory relation between $σ_v$, the density parameter $Ω$, and the galaxy correlation function $ξ(r)$ to infer the quantity $β=Ω^{0.6}/b = 0.60^{+0.13}_{-0.11}$ where $b$ is the linear bias parameter of optical galaxies and the uncertainties correspond to bootstrap resampling and an estimated cosmic variance added in quadrature. Our findings are consistent with the results of cluster abundances and redshift space distortion of the two-point correlation function. These statistical measurements suggest a low value of the density parameter $Ω\sim 0.4$ if optical galaxies are not strongly biased tracers of mass. △ Less

Submitted 19 October, 1999; originally announced October 1999.

Comments: Accepted for publication in MNRAS. 8 pages latex (mn.sty), including 7 figures

Effects of Distance Uncertainties on Determinations of the Galaxy Peculiar Velocity Field

Authors: Nelson Padilla, Manuel Merchán, Diego G. Lambas

Abstract: Our statistical studies comprise the mean relative velocity and velocity correlation for pairs of galaxies as a function of separation, and a bulk flow analysis determined on spheres of $10 h^{-1}$ $Mpc$ radius. In order to test the statistical analysis we use COBE normalized CDM numerical simulations with different density parameters and cosmological constant and take into account the Tully-Fis… ▽ More Our statistical studies comprise the mean relative velocity and velocity correlation for pairs of galaxies as a function of separation, and a bulk flow analysis determined on spheres of $10 h^{-1}$ $Mpc$ radius. In order to test the statistical analysis we use COBE normalized CDM numerical simulations with different density parameters and cosmological constant and take into account the Tully-Fisher (TF) scatter and a possible TF zero-point offset, as well as variations in the results of the simulations arising from different observer positions. We compare the results of the mock catalogs with samples of spiral galaxies taken from the Mark III catalog. △ Less

Submitted 29 April, 1998; originally announced April 1998.

Comments: 18 pages, 6 figures, to appear in ApJ