Search

Spectral data cubes of the interacting pair of galaxies NGC 2535 and NGC 2536 (the Arp 82 system) targeting bright emission lines in the visible band, obtained with the imaging Fourier transform spectrometer (iFTS) SITELLE attached to the Canada-France-Hawaii Telescope (CFHT), are presented. Analysis of H$\upalpha$ velocity maps reveals a bar in $\rm NGC\,2536$. In $\rm NGC\,2535$, we find strong non-circular motions outside the ocular ring, in the elliptical arc and tidal tails of $\rm NGC\,2535$ and a misalignment between the kinematic and photometric position angles. We detect 155 HII region complexes in the interacting pair of galaxies and determine oxygen abundances for 66 of them using different calibrators. We find, regardless of the indicator used, that the oxygen abundance distribution in $\rm NGC\,2536$ is shallow whereas, in $\rm NGC\,2535$, it is best fitted by two slopes, the break occurring beyond the ocular ring. The inner slope is comparable to the one observed in isolated normal star-forming galaxies but the outer slope is shallow. We present a numerical simulation of the interaction that reproduces the observed tidal features, kinematics, and metallicity distribution, to investigate the effect of the interaction on the galaxies. The model indicates that the galaxies have undergone a close encounter, strongly prograde for the primary, and are half way in their course to a second close encounter., Comment: 25 pages, 18 figures, accepted for publication in mnras

We present a series of chemodynamical simulations of Magellanic-like systems consisting of two interacting, equal-mass dwarf galaxies orbiting a massive host galaxy, including feedback and star formation, tides, and ram pressure. We study the star formation and chemical enrichment history of the dwarfs, and the production of a Magellanic Stream analogue. The dwarfs interact with each other through tidal forces, distorting their morphologies and triggering star formation. A stream is naturally produced as outflows, induced by feedback and interactions, are stretched by tidal forces. Counter to some recent simulations, we find that the leading arm persists even in the presence of ram pressure from the host galaxy. Interactions between the dwarfs and the host galaxies produce multiple kinematic components in the stream, as observed. A combination of ongoing star-formation and entrained low-metallicity gas causes the stream to have a complex chemical structure, with an average metallicity that is generally lower than that of the dwarfs., Comment: Accepted to ApJ

The formation and evolution of local brightest cluster galaxies (BCGs) is investigated by determining the stellar populations and dynamics from the galaxy core, though the outskirts and into the intracluster light (ICL). Integral spectroscopy of 23 BCGs observed out to 4 r_e is collected and high signal-to-noise regions are identified. Stellar population synthesis codes are used to determine the age, metallicity, velocity, and velocity dispersion of stars within each region. The intracluster light (ICL) spectra are best modeled with populations that are younger and less metal-rich than those of the BCG cores. The average BCG core age of the sample is 13.3$\pm$ 2.8 Gyr and the average metallicity is [Fe/H] = 0.30 $\pm$ 0.09, whereas for the ICL the average age is 9.2$\pm$3.5 Gyr and the average metallicity is [Fe/H] = 0.18$\pm$0.16. The velocity dispersion profile is seen to be rising or flat in most of the sample (17/23), and those with rising values reach the value of the host cluster's velocity dispersion in several cases. The most extended BCGs are closest to the peak of the cluster's X-ray luminosity. The results are consistent with the idea that the BCG cores and inner regions formed quickly and long ago, with the outer regions and ICL forming more recently, and continuing to assemble through minor merging. Any recent star formation in the BCGs is a minor component, and is associated with the cluster cool core status., Comment: 22 pages, 21 figures, MNRAS, accepted

By implementing a dynamic wind-tunnel model in a smoothed-particle chemodynamic/hydrodynamic simulation suite, we have investigated the effects of ram pressure and tidal forces on dwarf galaxies similar to the Magellanic Clouds, within host galaxies with gas and dark matter halos that are varied, to compare the relative effects of tides and ram pressure. We concentrate on how the distributions of metals are affected by interactions. We find that while ram pressure and tidal forces have some effect on dwarf galaxy outflows, these effects do not produce large differences in the metal distributions of the dwarf disks other than truncation in the outer regions in some cases, and that confinement from the host galaxy gas halo appears to be more significant than ram pressure stripping. We find that stochastic variations in the star formation rate can explain the remaining variations in disk metal properties. This raises questions on the cause of low metallicities in dwarf galaxies., Comment: Submitted to ApJ, under 2nd review (very minor revisions)

To investigate the role of bars in the chemical evolution of isolated disc galaxies, we performed a series of 39 gas dynamical simulations of isolated barred and unbarred galaxies with various masses, initial gas fractions, and AGN feedback models. The presence of a bar drives a substantial amount of gas toward the central region of the galaxy. In the most massive galaxies, this results in a violent starburst, followed by a drop in star formation resulting from gas exhaustion. The time delay between Type Ia and Type II supernovae explosions means that barred galaxies experience a rapid increase in [O/H] in the central region, and a much more gradual increase in [Fe/H]. In unbarred galaxies, star formation proceeds at a slow and steady rate, and oxygen and iron are produced at steady rates which are similar except for a time offset. Comparing the abundance ratios in barred and unbarred galaxies with the same central stellar mass $M_*$, we find in barred galaxies an enhancement of 0.07 dex in [O/H], 0.05 dex in [Fe/H], and 0.05 dex in [O/Fe]. The [O/H] enhancement is in excellent agreement with observations from the SDSS. The initial gas fraction has very little effect on the abundance ratios in barred and unbarred galaxies, unless the galaxies experience a starburst. We considered AGN-host galaxies located near the bottom of the AGN regime, $M_*>3\times10^{10}M_\odot$, where AGN feedback dominates over supernovae feedback. We found that the impact of AGN feedback on the central abundances is marginal., Comment: 16 pages, 11 figures, accepted for publication in MNRAS

We present a numerical study of the impact of AGN accretion and feedback on the star formation history of barred disc galaxies. Our goal is to determine whether the effect of feedback is positive (enhanced star formation) or negative (quenched star formation), and to what extent. We performed a series of 12 hydrodynamical simulations of disc galaxies, 10 barred and 2 unbarred, with various initial gas fractions and AGN feedback prescriptions. In barred galaxies, gas is driven toward the centre of the galaxy and causes a starburst, followed by a slow decay, while in unbarred galaxies the SFR increases slowly and steadily. AGN feedback suppresses star formation near the central black hole. Gas is pushed away from the black hole, and collides head-on with inflowing gas, forming a dense ring at a finite radius where star formation is enhanced. We conclude that both negative and positive feedback are present, and these effects mostly cancel out. There is no net quenching or enhancement in star formation, but rather a displacement of the star formation sites to larger radii. In unbarred galaxies, where the density of the central gas is lower, quenching of star formation near the black hole is more efficient, and enhancement of star formation at larger radii is less efficient. As a result, negative feedback dominates. Lowering the gas fraction reduces the star formation rate at all radii, whether or not there is a bar or an AGN., Comment: 18 pages, 17 figures. Accepted for publication in MNRAS

We performed a series of 29 gasdynamical simulations of disc galaxies, barred and unbarred, with various stellar masses, to study the impact of the bar on star formation history. Unbarred galaxies evolve very smoothly, with a star formation rate (SFR) that varies by at most a factor of three over a period of 2 Gyr. The evolution of barred galaxies is much more irregular, especially at high stellar masses. In these galaxies, the bar drives a substantial amount of gas toward the centre, resulting in a high SFR, and producing a starburst in the most massive galaxies. Most of the gas is converted into stars, and gas exhaustion leads to a rapid drop of star formation after the starburst. In massive barred galaxies (stellar mass M* > 2x10^10 Msun) the large amount of gas funnelled toward the centre is completely consumed by the starburst, while in lower-mass barred galaxies it is only partially consumed. Gas concentration is thus higher in lower-mass barred galaxies than it is in higher-mass ones. Even though unbarred galaxies funnelled less gas toward their centre, the lower SFR allows this gas to accumulate. At late times, the star formation efficiency is higher in barred galaxies than unbarred ones, enabling these galaxies to maintain a higher SFR with a smaller gas supply. Several properties, such as the global SFR, central SFR, or central gas concentration, vary monotonically with time for unbarred galaxies, but not for barred galaxies. Therefore one must be careful when comparing barred and unbarred galaxies that share one observational property, since these galaxies might be at very different stages of their respective evolution., Comment: 15 pages, 14 figures

The mass-metallicity relation shows that the galaxies with the lowest mass have the lowest metallicities. As most dwarf galaxies are in group environments, interaction effects such as tides could contribute to this trend. We perform a series of smoothed particle hydrodynamics (SPH) simulations of dwarf galaxies in external tidal fields to examine the effects of tides on their metallicities and metallicity gradients. In our simulated galaxies, gravitational instabilities drive gas inwards and produce centralized star formation and a significant metallicity gradient. Strong tides can contribute to these instabilities, but their primary effect is to strip the outer low-metallicity gas, producing a truncated gas disk with a large metallicity. This suggests that the role of tides on the mass-metallicity relation is to move dwarf galaxies to higher metallicities., Comment: Accepted to ApJ

We perform a series of smoothed particle hydrodynamics simulations of isolated dwarf galaxies to compare different metal mixing models. In particular, we examine the role of diffusion in the production of enriched outflows, and in determining the metallicity distributions of gas and stars. We investigate different diffusion strengths, by changing the pre-factor of the diffusion coefficient, by varying how the diffusion coefficient is calculated from the local velocity distribution, and by varying whether the speed of sound is included as a velocity term. Stronger diffusion produces a tighter [O/Fe]-[Fe/H] distribution in the gas, and cuts off the gas metallicity distribution function at lower metallicities. Diffusion suppresses the formation of low-metallicity stars, even with weak diffusion, and also strips metals from enriched outflows. This produces a remarkably tight correlation between "metal mass-loading" (mean metal outflow rate divided by mean metal production rate) and the strength of diffusion, even when the diffusion coefficient is calculated in different ways. The effectiveness of outflows at removing metals from dwarf galaxies and the metal distribution of the gas is thus dependent on the strength of diffusion. By contrast, we show that the metallicities of stars are not strongly dependent on the strength of diffusion, provided that some diffusion is present., Comment: Submitted to ApJ

Using a large N-body cosmological simulation combined with a subgrid treatment of galaxy formation, we study the formation and evolution of the galaxy and cluster population in a comoving volume (100 Mpc)^3 in a LCDM universe. At z = 0, our computational volume contains 1788 clusters with mass M_cl > 1.1x10^12 Msun, including 18 massive clusters with M_cl > 10^14 Msun. It also contains 1 088 797 galaxies with mass M_gal > 2x10^9 Msun and luminosity L > 9.5x10^5 Lsun. For each cluster, we identified the brightest cluster galaxy (BCG). We then computed the fraction f_BNC of clusters in which the BCG is not the closest galaxy to the center of the cluster in projection, and the ratio Dv/s, where Dv is the difference in radial velocity between the BCG and the whole cluster, and s is the radial velocity dispersion of the cluster. f_BNC increases from 0.05 for low-mass clusters (M_cl ~ 10^12 Msun) to 0.5 for high-mass ones (M_cl > 10^14 Msun), with no dependence on cluster redshift. The values of Dv/s vary from 0 to 1.8. These results are consistent with previous observational studies, and indicate that the central galaxy paradigm, which states that the BCG should be at rest at the center of the cluster, is usually valid, but exceptions are too common to be ignored. Analysis of the merger trees for the 18 most massive clusters in the simulation reveals that 16 of these clusters have experienced major mergers in the past. These mergers leave each cluster in a non-equilibrium state, but eventually the cluster settles into an equilibrium configuration, unless it is disturbed by another major merger. We found evidence that these mergers are responsible for the off-center positions and peculiar velocities of some BCGs. Our results thus support the merging-group scenario, in which some clusters form by the merger of smaller groups in which the galaxies have already formed., Comment: 32 pages, 11 figures (one in color). Accepted for publication in The Astrophysical Journal

We have modified our particle-by-particle adaptation of the made-to-measure (M2M) method, with the aim of modelling the Galactic disc from upcoming Galactic stellar survey data. In our new particle-by-particle M2M algorithm, PRIMAL, the observables of the target system are compared with those of the model galaxy at the position of the target stars, i.e. particles. The mass of the model particles are adjusted to reproduce the observables of the target system, and the gravitational potential is automatically adjusted by the changing mass of the particles. This paper builds upon our previous work, introducing likelihood-based velocity constraints in PRIMAL. In this paper we apply PRIMAL to barred disc galaxies created by a N-body simulation in a known dark matter potential, with no error in the observables. This paper demonstrates that PRIMAL can recover the radial profiles of the surface density, velocity dispersion in the radial and perpendicular directions, and the rotational velocity of the target discs, along with the apparent bar structure and pattern speed of the bar, especially when the reference frame is adjusted so that the bar angle of the target galaxy is aligned to that of the model galaxy at every timestep., Comment: 13 pages, 10 figures, Accepted 2013 April 16. Received 2013 April 8; in original form 2013 February 4

We have performed a series of chemodynamical simulations of barred disc galaxies. Our goal is to determine the physical processes responsible for the increase in the central gas-phase metallicity and of the central SFR observed in SDSS. We define a 2kpc diameter central aperture to approximate the integrated spectroscopic fibre measurements from the SDSS. The chemical evolution observed within this central region depends critically upon the relative size of the bar and the aperture, which evolves strongly with time. At t~0.5Gyr, a strong bar forms, whose length is considerably longer than the 2kpc aperture. The stars and gas lose angular momentum and follow elongated orbits that cause an intense mixing of the gas between the central region and its surroundings. During the next 1.5Gyr, the orbits of the gas contract significantly until the entire gas bar is contained in the 2kpc aperture, resulting in a net flux of gas into the central region. During this period, the metallicity in the central region increases steadily, and this enrichment is dominated by metal-rich gas that is flowing into the central region. The main result of this work is therefore that the observed enrichment in the centres of barred galaxies is not dominated by in-situ enrichment by stars formed in the centre. Rather, star formation occurs along the full length of the bar, much of which occurs initially outside the 2kpc aperture. About 50% of the metals that end up in the central region originate from this extended bar-long star formation, but flow into the central region due to loss of angular momentum. We conclude that there is no direct connection between central SFR and central metallicity. The central metallicity does not originate exclusively from central stars. Instead, the global SFR (especially along the bar) and the large-scale flow of enriched gas play a major role., Comment: Accepted for publication in MNRAS

Dust and gas energetics are incorporated into a cluster-scale simulation of star formation in order to study the effect of heating and cooling on the star formation process. We build on our previous work by calculating separately the dust and gas temperatures. The dust temperature is set by radiative equilibrium between heating by embedded stars and radiation from dust. The gas temperature is determined using an energy-rate balance algorithm which includes molecular cooling, dust-gas collisional energy transfer, and cosmic-ray ionization. The fragmentation proceeds roughly similarly to simulations in which the gas temperature is set to the dust temperature, but there are differences. The structure of regions around sink particles have properties similar to those of Class 0 objects, but the infall speeds and mass accretion rates were, on average, higher than those seen for regions forming only low-mass stars. The gas and dust temperature have complex distributions not well modeled by approximations that ignore the detailed thermal physics. There is no simple relationship between density and kinetic temperature. In particular, high density regions have a large range of temperatures, determined by their location relative to heating sources. The total luminosity underestimates the star formation rate at these early stages, before ionizing sources are included, by an order of magnitude. As predicted in our previous work, a larger number of intermediate mass objects form when improved thermal physics is included, but the resulting IMF still has too few low mass stars. However, if we consider recent evidence on core-to-star efficiencies, the match to the IMF is improved., Comment: 10 figures, 3 colored. Accepted for publication in ApJ

We combine a N-body simulation algorithm with a subgrid treatment of galaxy formation, mergers, and tidal destruction, and an observed conditional luminosity function Phi(L|M), to study the origin and evolution of galactic and extragalactic light inside a cosmological volume of size (100 Mpc)^3, in a concordance LCDM model. This algorithm simulates the growth of large-scale structures and the formation of clusters, the evolution of the galaxy population in clusters, the destruction of galaxies by mergers and tides, and the evolution of the intracluster light. We find that destruction of galaxies by mergers dominates over destruction by tides by about an order of magnitude at all redshifts. However, tidal destruction is sufficient to produce intracluster light fractions f_ICL that are sufficiently high to match observations. The bulk of the intracluster light (60%) is provided by intermediate galaxies of total masses 10^11 Msun-10^12 Msun and stellar masses 6x10^8 Msun-3x10^10 Msun that were tidally destroyed by even more massive galaxies. The contribution of low-mass galaxies to the intracluster light is small and the contribution of dwarf galaxies is negligible, even though, by numbers, most galaxies that are tidally destroyed are dwarfs. Tracking clusters back in time, we find that their values of f_ICL tend to increase over time, but can experience sudden changes that are sometimes non-monotonic. These changes occur during major mergers involving clusters of comparable masses but very different intracluster luminosities. Most of the tidal destruction events take place in the central regions of clusters. As a result, the intracluster light is more centrally concentrated than the galactic light. Our results support tidal destruction of intermediate-mass galaxies as a plausible scenario for the origin of the intracluster light., Comment: 13 figures, one in color. Accepted for publication in ApJ

We present a model to self-consistently describe the joint evolution of starburst galaxies and the galactic wind resulting from this evolution. We combine the population synthesis code Starburst99 with a semi-analytical model of galactic outflows and a model for the distribution and abundances of chemical elements inside the outflows. Starting with a galaxy mass, formation redshift, and adopting a particular form for the star formation rate, we describe the evolution of the stellar populations in the galaxy, the evolution of the metallicity and chemical composition of the interstellar medium (ISM), the propagation of the galactic wind, and the metal-enrichment of the intergalactic medium (IGM). In this paper, we study the properties of the model, by varying the mass of the galaxy, the star formation rate, and the efficiency of star formation. Our main results are the following: (1) For a given star formation efficiency f*, a more extended period of active star formation tends to produce a galactic wind that reaches a larger extent. If f* is sufficiently large, the energy deposited by the stars completely expels the ISM. Eventually, the ISM is being replenished by mass loss from supernovae and stellar winds. (2) For galaxies with masses above 10^11 Msun, the material ejected in the IGM always falls back onto the galaxy. Hence lower-mass galaxies are the ones responsible for enriching the IGM. (3) Stellar winds play a minor role in the dynamical evolution of the galactic wind, because their energy input is small compared to supernovae. However, they contribute significantly to the chemical composition of the galactic wind. We conclude that the history of the ISM enrichment plays a determinant role in the chemical composition and extent of the galactic wind, and therefore its ability to enrich the IGM., Comment: 17 pages, 14 color figures. Accepted for publication in MNRAS

Galactic outflows play a major role in the evolution of galaxies and the intergalactic medium (IGM). The energy deposited into the interstellar medium by supernovae and active galactic nuclei can accelerate the gas past the escape velocity, and eject it into the IGM. This will affect the subsequent evolution of the galaxy, by reducing or eliminating star formation, and quenching the accretion of matter onto the central AGN. Galactic outflows is the main process by which energy and processed interstellar matter is transported into the IGM. This affects the subsequent formation of other galaxies. The energy carried by outflows can strip protogalactic halos of their gas, preventing galaxies from forming. Conversely, the metals carried by outflows can modify the composition and cooling rates of the gas in protogalactic halos, favoring the formation of galaxies. In this paper, I review the various techniques used to simulate galactic outflows and their impact on galaxy and IGM evolution., Comment: 8 pages, 4 figures. Invited talk presented at the 11th Asian-Pacific Regional IAU Meeting, Chiang Mai, Thailand, July 28, 2011

We performed numerical simulations of mergers between gas-rich disc galaxies, which result in the formation of late-type galaxies. Stars formed during the merger end up in a thick disc that is partially supported by velocity dispersion and has high [alpha/Fe] ratios at all metallicities. Stars formed later end up in a thin, rotationally supported disc which has lower [alpha/Fe] ratios. While the structural and kinematical properties of the merger remnants depend strongly upon the orbital parameters of the mergers, we find a clear chemical signature of gas-rich mergers., Comment: 5 pages, 5 figures, IAU symposium No. 277 : Tracing the Ancestry of Galaxies

We investigate the large-scale influence of outflows from AGNs in enriching the IGM with metals in a cosmological context. We combine cosmological simulations of large scale structure formation with a detailed model of metal enrichment, in which outflows expand anisotropically along the direction of least resistance, distributing metals into the IGM. The metals carried by the outflows are generated by two separate stellar populations: stars located near the central AGN, and stars located in the greater galaxy. Using this algorithm, we performed a series of 5 simulations of the propagation of AGN-driven outflows in a cosmological volume of size (128/h Mpc)^3 in a Lambda-CDM universe, and analyze the resulting metal enrichment of the IGM. We found that the metallicity induced in the IGM is greatly dominated by AGNs having bolometric luminosity L > 10^9 L_sun, sources with 10^8 < L / L_sun < 10^9 having a negligible contribution. Our simulations produced an average IGM metallicity of [O/H] = -5 at z = 5.5, which then rises gradually, and remains relatively flat at a value [O/H] = -2.8 between z = 2 and z = 0. The ejection of metals from AGN host galaxies by AGN-driven outflows is found to enrich the IGM to > 10 - 20% of the observed values, the number dependent on redshift. The enriched IGM volume fractions are small at z > 3, then rise rapidly to the following values at z = 0: 6 - 10% of the volume enriched to [O/H] > -2.5, 14 - 24% volume to [O/H] > -3, and 34 - 45% volume to [O/H] > -4. At z > 2, there is a gradient of the induced enrichment, the metallicity decreasing with increasing IGM density, enriching the underdense IGM to higher metallicities, a trend more prominent with increasing anisotropy of the outflows. This can explain observations of metal-enriched low-density IGM at z = 3 - 4., Comment: 18 pages, 12 figures, submitted

We combine an analytic model for anisotropic outflows and galaxy formation with numerical simulations of large-scale structure and halo formation to study the impact of galactic outflows on the evolution of the IGM. We have simulated the evolution of a comoving volume (15 Mpc)^3 in the LCDM universe. We follow the formation of 20000-60000 galaxies and simulate the galactic outflows produced by these galaxies, for five outflow opening angles, alpha=60, 90, 120, 150, and 180 degrees (isotropic outflows). Anisotropic outflows follow the path of least resistance and thus travel preferentially into low-density regions, away from cosmological structures where galaxies form. These anisotropic outflows are less likely to overlap with one another, or to hit pre-galactic collapsing halos and strip them of their gas, preventing a galaxy from forming. Going from 180 deg to 60 deg, the number of galaxies that actually form doubles, producing twice as many outflows, and these outflows overlap to a lesser extent. As a result, the metal volume filling factor of the IGM goes from 8% for isotropic outflows up to 28% for anisotropic ones. High density regions are more efficiently enriched than low density ones (~80% compared to ~20% by volume), even though most enriched regions are low densities. Increasing the anisotropy of outflows increases the extent of enrichment at all densities, low and high. This is in part because anisotropic outflows are more numerous. When this effect is factored-out, we find that the probability a galaxy will enrich systems at densities up to 10 rho_mean is higher for increasingly anisotropic outflows. This is an effect of the dynamical evolution of the IGM. Anisotropic outflows expand preferentially into underdense gas, but that gas can later accrete onto overdense structures., Comment: Submitted to ApJ. Figures 4, 5 and 6 are not included because of their large sizes. They can be downloaded from: http://www.astro.phy.ulaval.ca/staff/hugo/web_out/figures456.tar.gz

We investigate the effect of heating by luminosity sources in a simulation of clustered star formation. Our heating method involves a simplified continuum radiative transfer method that calculates the dust temperature. The gas temperature is set by the dust temperature. We present the results of four simulations, two simulations assume an isothermal equation of state and the two other simulations include dust heating. We investigate two mass regimes, i.e., 84 Msun and 671 Msun, using these two different energetics algorithms. The mass functions for the isothermal simulations and simulations which include dust heating are drastically different. In the isothermal simulation, we do not form any objects with masses above 1 Msun. However, the simulation with dust heating, while missing some of the low-mass objects, forms high-mass objects (~20 Msun) which have a distribution similar to the Salpeter IMF. The envelope density profiles around the stars formed in our simulation match observed values around isolated, low-mass star-forming cores. We find the accretion rates to be highly variable and, on average, increasing with final stellar mass. By including radiative feedback from stars in a cluster-scale simulation, we have determined that it is a very important effect which drastically affects the mass function and yields important insights into the formation of massive stars., Comment: 19 pages, 28 figures. See http://www.astro.phy.ulaval.ca/staff/hugo/dust/ms_dust.big.pdf for high resolution version of document

We used an N-body smoothed particle hydrodynamics algorithm, with a detailed treatment of star formation, supernovae feedback, and chemical enrichment, to perform eight simulations of mergers between gas-rich disc galaxies. We vary the mass ratio of the progenitors, their rotation axes, and their orbital parameters and analyze the kinematic, structural, and chemical properties of the remnants. Six of these simulations result in the formation of a merger remnant with a disc morphology as a result of the large gas-fraction of the remnants. We show that stars formed during the merger (a sudden starburst occur in our simulation and last for 0.2-0.3 Gyr) and those formed after the merger have different kinematical and chemical properties. The first ones are located in thick disc or the halo. They are partially supported by velocity dispersion and have high [alpha/Fe] ratios even at metallicities as high as [Fe/H]=-0.5. The former ones -- the young component -- are located in a thin disc rotationally supported and have lower [alpha/Fe] ratios. The difference in the rotational support of both components results in the rotation of the thick disc lagging that of the thin disc by as much as a factor of two, as recently observed.We find that, while the kinematic and structural properties of the merger remnant depends strongly upon the orbital parameters of the mergers, there is a remarkable uniformity in the chemical properties of the mergers. This suggests that general conclusions about the chemical signature of gas-rich mergers can be drawn., Comment: 17 pages, 20 figures, to be published in MNRAS

We investigate the cosmological-scale influence of outflows driven by AGNs on metal enrichment of the intergalactic medium. AGNs are located in dense cosmological structures which tend to be anisotropic. We designed a semi-analytical model for anisotropic AGN outflows which expand away along the direction of least resistance. This model was implemented into a cosmological numerical simulation algorithm for simulating the growth of large-scale structure in the universe. Using this modified algorithm, we perform a series of 9 simulations inside cosmological volumes of size $(128 h^{-1}{\rm Mpc})^3$, in a concordance $\Lambda$CDM universe, varying the opening angle of the outflows, the lifetimes of the AGNs, their kinetic fractions, and their level of clustering. For each simulation, we compute the volume fraction of the IGM enriched in metals by the outflows. The resulting enriched volume fractions are relatively small at $z \gtrsim 2.5$, and then grow rapidly afterward up to $z = 0$. We find that AGN outflows enrich from 65% to 100% of the entire universe at the present epoch, for different values of the model parameters. The enriched volume fraction depends weakly on the opening angle of the outflows. However, increasingly anisotropic outflows preferentially enrich underdense regions, a trend found more prominent at higher redshifts and decreasing at lower redshifts. The enriched volume fraction increases with increasing kinetic fraction and decreasing AGN lifetime and level of clustering., Comment: 19 pages, 16 figures, submitted. The version uploaded here does not contain Figs 5, 6 & 7, because of their large sizes. Those can be found along with the full paper at: http://www.astro.phy.ulaval.ca/staff/paramita/AllPages/Talks-Posters/Papers_Thesis/ms_AGNoutflow.pdf

We simulate anisotropic outflows of AGN, and investigate the large-scale impact of the cosmological population of AGN outflows over the Hubble time by performing N-body LambdaCDM simulations. Using the observed quasar luminosity function to get the redshift and luminosity distribution, and analytical models for the outflow expansion, AGNs are allowed to evolve in a cosmological volume. By the present epoch, 13 - 25% of the total volume is found to be pervaded by AGN outflows, with 10^{-9} G magnetic field., Comment: 4 pages, 2 figures. Contributed talk at conference "The Central Kiloparsec: Active Galactic Nuclei and Their Hosts", Ierapetra, Crete, 4-6 June 2008. To appear in Volume 79 of the Memorie della Societa Astronomica Italiana

We use a series of ray-tracing experiments to determine the magnification distribution of high-redshift sources by gravitational lensing. We determine empirically the relation between magnification and redshift, for various cosmological models. We then use this relation to estimate the effect of lensing on the determination of the cosmological parameters from observations of high-z supernovae. We found that, for supernovae at redshifts z<1.8 the effect of lensing is negligible compared to the intrinsic uncertainty in the measurements. Using mock data in the range 1.82. Therefore if supernovae up to these redshifts were ever discovered, it is still the ones in the range 0.3

The main goal of this paper is to compare the relative importance of destruction by tides, vs. destruction by mergers, in order to assess if tidal destruction of dwarf galaxies in clusters is a viable scenario for explaining the origin of intracluster stars. We have designed a simple algorithm for simulating the evolution of isolated clusters. The distribution of galaxies in the cluster is evolved using a direct gravitational N-body algorithm combined with a subgrid treatment of physical processes such as mergers, tidal disruption, and galaxy harassment. Using this algorithm, we have performed a total of 227 simulations. Our main results are (1) destruction of dwarf galaxies by mergers dominates over destruction by tides, and (2) the destruction of dwarf galaxies by tides is sufficient to explain the observed intracluster light in clusters., Comment: 21 pages, submitted to Journal of Astrophysics and Astronomy

We show that radiative feedback due to reionization has a pronounced effect on the extent of mechanical feedback due to galactic outflows. The photoionization of the Intergalactic Medium (IGM) suppresses low-mass galaxy formation by photoheating the gas and limiting atomic line cooling. The number of low-mass galaxies is central for the enrichment of the IGM as these objects have the capacity to enrich a significant fraction (by volume) of the Universe. We use a modified version of our galactic outflow model, combined with a simple criterion for suppression, to investigate the potential impact upon the IGM. We find that this suppression strongly reduces the enrichment of the IGM and is sensitive to the reionization history. We also investigate the contribution of halos of different masses with varying degrees of suppression., Comment: 4 pages, 3 figures, class file included, accepted by ApJ Letters, minor changes and expanded Figure 3

We use cosmological, chemo-dynamical, smoothed particle hydrodynamical simulations of Milky-Way-analogue galaxies to find the expected present-day distributions of both metal-free stars that formed from primordial gas and the oldest star populations. We find that metal-free stars continue to form until z~4 in halos that are chemically isolated and located far away from the biggest progenitor of the final system. As a result, if the Population III initial mass function allows stars with low enough mass to survive until z=0 (< 0.8 Msol), they would be distributed throughout the Galactic halo. On the other hand, the oldest stars form in halos that collapsed close to the highest density peak of the final system, and at z=0 they are located preferentially in the central region of the Galaxy, i.e., in the bulge. According to our models, these trends are not sensitive to the merger histories of the disk galaxies or the implementation of supernova feedback. Furthermore, these full hydrodynamics results are consistent with our N-body results in Paper I, and lend further weight to the conclusion that surveys of low-metallicity stars in the Galactic halo can be used to directly constrain the properties of primordial stars. In particular, they suggest that the current lack of detections of metal-free stars implies that their lifetimes were shorter than a Hubble time, placing constraints on the metal-free initial mass function., Comment: Accepted by ApJ. Emulate ApJ style

We employ N-body, smoothed particle hydrodynamical simulations, including detailed treatment of chemical enrichment, to follow a gas-rich merger which results in a galaxy with disk morphology. We trace the kinematic, structural and chemical properties of stars formed before, during, and after the merger. We show that such a merger produces two exponential disk components, with the older, hotter component having a scale-length 20% larger than the later-forming, cold disk. Rapid star formation during the merger quickly enriches the protogalactic gas reservoir, resulting in high metallicities of the forming stars. These stars form from gas largely polluted by Type II supernovae, which form rapidly in the merger-induced starburst. After the merger, a thin disk forms from gas which has had time to be polluted by Type Ia supernovae. Abundance trends are plotted, and we examine the proposal that increased star formation during gas-rich mergers may explain the high alpha-to-iron abundance ratios which exist in the relatively high-metallicity thick disk component of the Milky Way., Comment: 10 pages, 1 color figure. Accepted for publication in ApJ

We have designed an analytical model for the evolution of anisotropic galactic outflows. These outflows follow the path of least resistance, and thus travel preferentially into low-density regions, away from cosmological structures where galaxies form. We show that anisotropic outflows can significantly enrich low-density systems with metals., Comment: Proceedings of Chemodynamics 2006, Lyon, 2 pages, 1 figure, style file included

We have developed an analytical model for the evolution of anisotropic galactic outflows. These outflows follow the path of least resistance, and thus travel preferentially into low-density regions, away from cosmological structures where galaxies form. We show that anisotropic outflows can significantly enrich low-density systems, while reducing the enrichment of overdense regions., Comment: Proceedings of Cosmic Frontiers 2006, Durham, 2 pages, 1 figure, style file included

We have developed an analytical model to describe the evolution of anisotropic galactic outflows. With it, we investigate the impact of varying opening angle on galaxy formation and the evolution of the IGM. We have implemented this model in a Monte Carlo algorithm to simulate galaxy formation and outflows in a cosmological context. Using this algorithm, we have simulated the evolution of a comoving volume of size [12h^(-1)Mpc]^3 in the LCDM universe. Starting from a Gaussian density field at redshift z=24, we follow the formation of ~20,000 galaxies, and simulate the galactic outflows produced by these galaxies. When these outflows collide with density peaks, ram pressure stripping of the gas inside the peak may result. This occurs in around half the cases and prevents the formation of galaxies. Anisotropic outflows follow the path of least resistance, and thus travel preferentially into low-density regions, away from cosmological structures (filaments and pancakes) where galaxies form. As a result, the number of collisions is reduced, leading to the formation of a larger number of galaxies. Anisotropic outflows can significantly enrich low-density systems with metals. Conversely, the cross-pollution in metals of objects located in a common cosmological structure, like a filament, is significantly reduced. Highly anisotropic outflows can travel across cosmological voids and deposit metals in other, unrelated cosmological structures., Comment: 32 pages, 9 figures (2 color). Revised version accepted in ApJ

The H atoms inside minihalos (i.e. halos with virial temperatures T_vir < 10^4 K, in the mass range roughly from 10^4 M_sun to 10^8 M_sun) during the cosmic dark ages in a LambdaCDM universe produce a redshifted background of collisionally-pumped 21-cm line radiation which can be seen in emission relative to the cosmic microwave background (CMB). Previously, we used semi-analytical calculations of the 21-cm signal from individual halos of different mass and redshift and the evolving mass function of minihalos to predict the mean brightness temperature of this 21-cm background and its angular fluctuations. Here we use high-resolution cosmological N-body and hydrodynamic simulations of structure formation at high redshift (z>8) to compute the mean brightness temperature of this background from both minihalos and the intergalactic medium (IGM) prior to the onset of Ly-alpha radiative pumping. We find that the 21-cm signal from gas in collapsed, virialized minihalos dominates over that from the diffuse shocked gas in the IGM., Comment: 10 pages, 8 figures; Accepted for publication in ApJ, with minor revision in text

In order to explain disk galaxy formation within the hierarchical structure formation, it seems that gas rich mergers must play an important role. We review here our previous studies which have shown the importance of mergers at high redshift being gas rich, in the formation of both the stellar halo and thick disk components of disk galaxies. Regulation of star formation in the building blocks of our galaxy is required to form a low mass low metallicity stellar halo. This regulation results in high redshift, gas rich mergers during which the thick disk forms. In these proceedings, we categorise stars from our simulated disk galaxy into thin and thick disk components by using the Toomre diagram., Comment: To appear in the proceedings of "Island Universes: Structure and Evolution of Disk Galaxies" (Terschelling, the Netherlands)

N-body simulations find a universal structure for the halos which result from the nonlinear growth of Gaussian-random-noise density fluctuations in the CDM universe. This talk summarized our attempts to derive and explain this universal structure by analytical approximation and simplified models. As an example, we show here that a 1D spherical infall model involving a fluid approximation derived from the Boltzmann equation can explain not only the halo density profile but its phase-space density profile, as well., Comment: To appear in the proceedings of 21st IAP colloquium, "Mass Profiles and Shapes of Comological Structures"

The H atoms inside minihalos (i.e. halos with virial temperatures T_vir < 10^4 K, in the mass range roughly from 10^4 M_sun to 10^8 M_sun) during the cosmic dark ages in a LCDM universe produce a redshifted background of collisionally-pumped 21-cm line radiation which can be seen in emission relative to the cosmic microwave background (CMB). Previously, we used semi-analytical calculations of the 21-cm signal from individual halos of different mass and redshift and the evolving mass function of minihalos to predict the mean brightness temperature of this 21-cm background and its angular fluctuations. Here we use high-resolution cosmological N-body and hydrodynamic simulations of structure formation at high redshift (z > 8) to compute the mean brightness temperature of this background from both minihalos and the intergalactic medium (IGM) prior to the onset of Ly-alpha radiative pumping. We find that the 21-cm signal from gas in collapsed, virialized minihalos dominates over that from the diffuse shocked gas in the IGM., Comment: 8 pages, 5 figures. To appear in proceedings of UC Irvine May 2005 workshop on "First Light & Reionization", eds. E. Barton & A. Cooray, New Astronomy Reviews, in press

Observations suggest that the structural parameters of disk galaxies have not changed greatly since redshift 1. We examine whether these observations are consistent with a cosmology in which structures form hierarchically. We use SPH/N-body galaxy-scale simulations to simulate the formation and evolution of Milky-Way-like disk galaxies by fragmentation, followed by hierarchical merging. The simulated galaxies have a thick disk, that forms in a period of chaotic merging at high redshift, during which a large amount of alpha-elements are produced, and a thin disk, that forms later and has a higher metallicity. Our simulated disks settle down quickly and do not evolve much since redshift z~1, mostly because no major mergers take place between z=1 and z=0. During this period, the disk radius increases (inside-out growth) while its thickness remains constant. These results are consistent with observations of disk galaxies at low and high redshift., Comment: Will appear in Proceedings of the 9th Asian Pacific Regional IAU Meeting

This document describes analytical and numerical techniques for the generation of Gaussian density fields, which represent cosmological density perturbations. The mathematical techniques involved in the generation of density harmonics in k-space, the filtering of the density fields, and the normalization of the power spectrum to the measured temperature fluctuations of the Cosmic Microwave Background, are presented in details. These techniques are well-known amongst experts, but the current literature lacks a formal description. I hope that this technical report will prove useful to new researchers moving into this field, sparing them the task of reinventing the wheel., Comment: Technical report, 24 pages, 7 figures. Revised version

We present a series of simulations of the fragmentation of a molecular cloud, leading to the formation of a cluster of protostellar cores. The purpose of these simulations is to address a specific numerical problem called artificial fragmentation, that plagues SPH simulations of cloud fragmentation. We argue that this is a serious problem that needs to be addressed, and that the only reasonable and practical way to address it is to use a relatively new technique called particle splitting. Our largest simulation has an effective resolution of 256^3 particles (much higher than most previous SPH simulations of cloud fragmentation) and results in the formation of a dense cluster containing ~3000 protostellar cores. It is the first simulation of this kind to properly resolve the Jeans mass throughout the entire system, at all times, thus preventing artificial fragmentation., Comment: 47 pages, 15 figures (2 grayscale, one color), ApJ Suppl, in press

We use a series of ray-tracing experiments to determine the magnification distribution of high-redshift sources by gravitational lensing. We determine empirically the relation between magnification and redshift, for various cosmological models. We then use this relation to estimate the effect of lensing on the determination of the cosmological parameters from observations of high-z supernovae. We found that, for supernovae at redshifts z<1.8, the effect of lensing is negligible compared to the intrinsic uncertainty in the measurements. Using mock data in the range 1.8

The recently emerging conviction that thick disks are prevalent in disk galaxies, and their seemingly ubiquitous old ages, means that the formation of the thick disk, perhaps more than any other component, holds the key to unravelling the evolution of the Milky Way, and indeed all disk galaxies. In Paper I, we proposed that the thick disk was formed in an epoch of gas rich mergers, at high redshift. This hypothesis was based on comparing N-body/SPH simulations to a variety of Galactic and extragalactic observations, including stellar kinematics, ages and chemical properties.Here examine our thick disk formation scenario in light of the most recent observations of extragalactic thick disks. In agreement, our simulted thick disks are old and relatively metal rich, with V-I colors that do not vary significantly with distance from the plane. Further, we show that our proposal results in an enhancement of alpha-elements in thick disk stars as compared with thin disk stars, consistent with observations of the relevant populations of the Milky Way. We also find that our scenario naturally leads to the formation of an old metal weak stellar halo population with high alpha-element abundances., Comment: submitted to ApJ

We present the summary of the theoretical aspects and algorithms used in an undergraduate (JRG) summer project based on numerical N-body simulations of collisionless systems. First, we review the importance of numerical N-body simulations in astrophysics. We introduce the different codes used and their performances. We then introduce four famous density profiles : Hernquist, NFW, truncated isothermal sphere, and Plummer. The history of these profiles and their dynamical properties are discussed in the third section. In the fourth section, we present the Barnes & Hut tree code, its features and performances. We then explain how to build and incorporate a multiple time stepping scheme in a tree code. The fifth section is dedicated to the different physical measurements we used to characterize the dynamics of the N-body system. Finally, we describe the future work that will be done with these codes, mainly the study of the adiabatic growth of a black hole at the center of a spherical distribution of stars., Comment: 18 pages, 6 figures (3 in color). Technical report

Dark-matter halos are the scaffolding around which galaxies and clusters are built. They form when the gravitational instability of primordial density fluctuations causes regions which are denser than average to slow their cosmic expansion, recollapse, and virialize. Understanding the equilibrium structure of these halos is thus a prerequisite for understanding galaxy and cluster formation. Numerical N-body simulations of structure formation from Gaussian-random-noise initial conditions in the CDM universe find a universal internal structure for halos. Objects as different in size and mass as dwarf spheroidal galaxies and galaxy clusters are predicted to have halos with the same basic structure when properly rescaled, independent of halo mass, of the shape of the power spectrum of primordial density fluctuations, and of the cosmological background parameters. This remarkable universality is a fundamental prediction of the CDM model, but our knowledge is limited to the ``empirical'' N-body simulation results, with little analytical understanding. We summarize here our attempts to fill this gap, in an effort to derive and give physical insight to the numerical results and extend them beyond the range of numerical simulation: (1) Simulated halos which form from highly simplified initial conditions involving gravitational instability in a cosmological pancake show that many of the universal properties of CDM halos are generic to cosmological gravitational collapse and do not require Gaussian-random-noise density fluctuations or hierarchical clustering. (2) A fluid approximation derived from the Boltzmann equation yields an analytical theory of halo dynamics which can explain many of the N-body results ... (abridged), Comment: invited article to appear in Progress in Dark Matter Research (Nova Science Publishers, Inc), 76 pages, 50 figures

We study the gravitational lensing of high-redshift sources in a LCDM universe. We have performed a series of ray-tracing experiments, and selected a subsample of cases of strong lensing (multiple images, arcs, and Einstein rings). For each case, we identify a massive galaxy that is primarily responsible for lensing, and studied how the various density inhomogeneities along the line of sight (other galaxies, background matter) affect the properties of the image. The matter located near the lensing galaxy, and physically associated with it, has a small effect. The background matter increases the magnification by a few percents at most, while nearby galaxies can increase it by up to about 10 percent. The effect on the image separation is even smaller. The only significant effect results from the random alignment of physically unassociated galaxies, which can increase the magnification by factors of several, create additional images, and turn arcs into rings. We conclude that the effect of environment on strong lensing in negligible in general, and might be important only in rare cases. We show that our conclusion does not depend on the radial density profile of the galaxies responsible for lensing., Comment: 23 pages, 7 figures (one in color). Accepted for publication in The Astrophysical Journal. Minor typos corrected

The gravitational lensing properties of cosmological halos depend upon the mass distribution within each halo. The description of halos as nonsingular, truncated isothermal spheres, a particular solution of the isothermal Lane-Emden equation (suitably modified for Lambda nonzero), has proved to be a useful approximation for the halos which form from realistic initial conditions in a CDM universe. The nonsingular TIS model reproduces many of the quantitative features of the N-body results for CDM halos, except in the very center, where CDM N-body halos show density profiles which vary as rho ~ r^(-alpha), alpha>1, instead of a small flat core. Possible discrepancies between these cuspy halo predictions of the CDM N-body simulations and observations of the inner mass profiles of dwarf and LSB disk galaxies based upon their rotation curves and of clusters based upon strong lensing measurements have led to a search for other diagnostics. A description of the lensing by TIS halos would be useful in this regard, as a self-consistent model for CDM halos in a proper cosmological context, nonsingular but otherwise consistent with the CDM N-body results. We derive here the basic lensing properties of individual TIS halos. For comparison, we also consider three singular profiles: the Navarro-Frenk-White density profile, the singular isothermal sphere, and the Schwarzschild lens. For all profiles, we compute the critical curves and caustics, the image separation, the magnification and brightness ratio, the shear, the time delay, and the average shear inside the tidal radius. This provides tools for studying the statistical properties of lensing by TIS and other lenses in the context of a theory of cosmological structure formation., Comment: 23 pages, 20 figures, submitted to MNRAS. Identical to version 1, except for typo in equation (17) corrected

We present the first gas dynamical simulations of the photoevaporation of cosmological minihalos overtaken by the ionization fronts which swept through the IGM during reionization in a LCDM universe, including the effects of radiative transfer. We demonstrate the phenomenon of I-front trapping inside minihalos, in which the weak, R-type fronts which traveled supersonically across the IGM decelerated when they encountered the dense, neutral gas inside minihalos, becoming D-type I-fronts, preceded by shock waves. For a minihalo with virial temperature T_vir < 10,000 K, the I-front gradually burned its way through the minihalo which trapped it, removing all of its baryonic gas by causing a supersonic, evaporative wind to blow backwards into the IGM, away from the exposed layers of minihalo gas just behind the advancing I-front. Such hitherto neglected feedback effects were widespread during reionization. N-body simulations and analytical estimates of halo formation suggest that sub-kpc minihalos such as these, with T_vir < 10,000 K, were so common as to dominate the absorption of ionizing photons. This means that previous estimates of the number of ionizing photons per H atom required to complete reionization which neglected this effect may be too low. Regardless of their effect on the progress of reionization, however, the minihalos were so abundant that random lines of sight thru the high-z universe should encounter many of them, which suggests that it may be possible to observe the processes described here in the absorption spectra of distant sources., Comment: 4 pages, 3 figures, to appear in, The Emergence of Cosmic Structure, the 13th Annual October Astrophysics Conference in Maryland, S. Holt & C. Reynolds, eds (AIP)

In the standard Cold Dark Matter (CDM) theory of structure formation, virialized minihalos (with T_vir = 10,000 K or less) form in abundance at high redshift (z > 6), during the cosmic "dark ages." The hydrogen in these minihalos, the first nonlinear baryonic structures to form in the universe, is mostly neutral and sufficiently hot and dense to emit strongly at the 21-cm line. We calculate the emission from individual minihalos and the radiation background contributed by their combined effect. Minihalos create a "21-cm forest" of emission lines. We predict that the angular fluctuations in this 21-cm background should be detectable with the planned LOFAR and SKA radio arrays, thus providing a direct probe of structure formation during the "dark ages." Such a detection will serve to confirm the basic CDM paradigm while constraining the background cosmology parameters, the shape of the power-spectrum of primordial density fluctuations, the onset and duration of the reionization epoch, and the conditions which led to the first stars and quasars. We present results here for the currently-favored, flat LCDM model, for different tilts of the primordial power spectrum. These minihalos will also cause a "21-cm forest" of absorption lines, as well, in the spectrum of radio continuum sources at high redshift, if the latter came into existence before the end of reionization., Comment: 4 pages, 2 figures, to appear in The Emergence of Cosmic Structure, the 13th Annual October Astrophysics Conference in Maryland, S. Holt & C. Reynolds, eds (AIP)

Hydrogen atoms inside virialized minihalos generate a radiation background from redshifted 21-cm line emission whose angular fluctuations reflect clustering during before and during reionization. We have shown elsewhere that this emission may be detectable with the planned Low Frequency Array (LOFAR) and Square Kilometer Array (SKA) in a flat LambdaCDM. This is a direct probe of structure during the ``Dark Ages'' at redshifts z>6 and down to smaller scales than have previously been constrained. In our original calculation, we used a standard approximation known as the ``linear bias'' [e.g. Mo & White (1996)]. Here we improve upon that treatment by considering the effect of nonlinear clustering. To accomplish this, we develop a new analytical method for calculating the nonlinear Eulerian bias of halos, which should be useful for other applications as well. Predictions of this method are compared with the results of LambdaCDM N-body simulations, showing significantly better agreement than the standard linear bias approximation. When applied to the 21-cm background from minihalos, our formalism predicts fluctuations that differ from our original predictions by up to 30% at low frequencies and small scales. However, within the range of frequencies and angular scales at which the signal could be observable by LOFAR and SKA as currently planned, the differences are small and our original predictions prove robust. Our results indicate that while a smaller frequency bandwidth of observation leads to a higher signal that is more sensitive to nonlinear effects, this effect is counteracted by the lowered sensitivity of the radio arrays. We calculate the best frequency bandwidth for these observations to be Delta nu~2 MHz.(abbreviated), Comment: 10 pages, 5 figures, submitted to MNRAS, changed to accomodate the comments of the referee, Figure 1 is with degraded resolution due to archive file size restrictions, full-resolution figure available from Authors upon request

Adaptive SPH and N-body simulations were carried out to study the collapse and evolution of dark matter halos that result from the gravitational instability and fragmentation of cosmological pancakes. Such halos resemble those formed by hierarchical clustering in a CDM universe and serve as a convenient test-bed model for studying halo dynamics. Our halos are in approximate virial equilibrium and roughly isothermal, as in CDM simulations. Their density profiles agree quite well with the fit to N-body results for CDM halos by Navarro, Frenk, & White (NFW). This test-bed model enables us to study the evolution of individual halos. The masses of our halos evolve in three stages: an initial collapse, continuous infall, and a final stage in which infall tapers off as a result of finite mass supply. In the continuous infall stage, halo mass grows at the rate expected for self-similar spherical infall, with M(a) proportional to the scale factor a. After the end of initial collapse at a=a_0, the concentration parameter grows linearly with a, c(a)~4a/a_0. The virial ratio 2T/|W| just after virialization is about 1.35, as predicted by the truncated isothermal sphere model and consistent with the value expected for a virialized halo in which mass infall contributes an effective surface pressure. Thereafter, the virial ratio evolves towards the value expected for an isolated halo, 2T/|W|~1. This mass accretion history and evolution of concentration parameter are very similar to those reported recently in N-body simulations of CDM. We therefore conclude that the fundamental properties of halo formation and evolution are generic to the formation of cosmological halos by gravitational instability and are not limited to hierarchical collapse scenarios or even to Gaussian-random-noise initial conditions., Comment: 2 pages, 4 figures, to appear in proceedings of "Galaxy Evolution: Theory and Observations" eds. V. Avila-Reese, C. Firmani, C. Frenk, & C. Allen, RevMexAA SC (2002)

The gravitational lensing properties of cosmological halos depend upon the mass distribution within each halo. The description of halos as nonsingular, truncated isothermal spheres, a particular solution of the isothermal Lane- Emden equation (suitably modified for Lambda non equal 0), has proved to be a useful approximation for the halos which form from realistic initial conditions in a CDM universe. We derive here the basic lensing properties of such halos, including the image separation, magnification, shear, and time-delay. We also provide analytical expressions for the critical curves and caustics. We show how the scale-free results we derive yield scale-dependent lensing properties which depend upon the cosmological background universe and the mass and collapse redshift of the lensing halos, according to the truncated isothermal sphere (TIS) model of CDM halos derived elsewhere. We briefly describe the application of these results to the currently-favored Lambda-CDM universe., Comment: 12 pages, 13 figures, submitted to MNRAS

In the standard Cold Dark Matter (CDM) theory of structure formation, virialized minihalos (with T_vir}<10^4 K) form in abundance at high redshift (z>6), during the cosmic ``dark ages.'' The hydrogen in these minihalos, the first nonlinear baryonic structures to form in the universe, is mostly neutral and sufficiently hot and dense to emit strongly at the 21-cm line. We calculate the emission from individual minihalos and the radiation background contributed by their combined effect. Minihalos create a ``21-cm forest'' of emission lines. We predict that the angular fluctuations in this 21-cm background should be detectable with the planned LOFAR and SKA radio arrays, thus providing a direct probe of structure formation during the ``dark ages.'' Such a detection will serve to confirm the basic CDM paradigm while constraining the shape of the power-spectrum of primordial density fluctuations down to much smaller scales than have previously been constrained, the onset and duration of the reionization epoch, and the conditions which led to the first stars and quasars. We present results here for the currently-favored, flat $\Lambda$CDM model, for different tilts of the primordial power spectrum., Comment: Revised version. 10 pages, 4 figures, re-submitted to ApJ Letters