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.

SIGNALS, the Star formation, Ionized Gas, and Nebular Abundances Legacy Survey, is a large observing program designed to investigate massive star formation and HII regions in a sample of local extended galaxies. The program will use the imaging Fourier transform spectrograph SITELLE at the Canada-France-Hawaii Telescope. Over 355 hours (54.7 nights) have been allocated beginning in fall 2018 for eight consecutive semesters. Once completed, SIGNALS will provide a statistically reliable laboratory to investigate massive star formation, including over 50 000 resolved HII regions : the largest, most complete, and homogeneous database of spectroscopically and spatially resolved extragalactic HII regions ever assembled. For each field observed, three datacubes covering the spectral bands of the filters SN1 (363 -386 nm), SN2 (482 - 513 nm), and SN3 (647 - 685 nm) are gathered. The spectral resolution selected for each spectral band is 1000, 1000, and 5000, respectively. As defined, the project sample will facilitate the study of small-scale nebular physics and many other phenomena linked to star formation at a mean spatial resolution of 20 pc. This survey also has considerable legacy value for additional topics including planetary nebulae, diffuse ionized gas, andsupernova remnants. The purpose of this paper is to present a general outlook of the survey, notably the observing strategy, galaxy sample, and science requirements.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Recent observations of column densities in molecular clouds find lognormal distributions with power-law high-density tails. These results are often interpreted as indications that supersonic turbulence dominates the dynamics of the observed clouds. We calculate and present the column-density distributions of three clouds, modeled with very different techniques, none of which is dominated by supersonic turbulence. The first star-forming cloud is simulated using smoothed particle hydrodynamics (SPH); in this case gravity, opposed only by thermal-pressure forces, drives the evolution. The second cloud is magnetically subcritical with subsonic turbulence, simulated using nonideal MHD; in this case the evolution is due to gravitationally-driven ambipolar diffusion. The third cloud is isothermal, self-gravitating, and has a smooth density distribution analytically approximated with a uniform inner region and an r^-2 profile at larger radii. We show that in all three cases the column-density distributions are lognormal. Power-law tails develop only at late times (or, in the case of the smooth analytic profile, for strongly centrally concentrated configurations), when gravity dominates all opposing forces. It therefore follows that lognormal column-density distributions are generic features of diverse model clouds, and should not be interpreted as being a consequence of supersonic turbulence.

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.

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.

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.

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.

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.

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.

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<z<8, we show that the effect of lensing can become significant. Hence, if a population of very-high-z supernovae was ever discovered, it would be crucial to fully understand the effect of lensing, before these SNe could be used to constrain cosmological models. We show that the distance moduli m-M for an open CDM universe and a LCDM universe are comparable at z>2. Therefore if supernovae up to these redshifts were ever discovered, it is still the ones in the range 0.3<z<1 that would distinguish these two models.

The difference in density profiles of the contributions from different density peaks to dark matter halos results in certain expectations about the Milky Way's stellar halo. We cut our simulated halo stars into two populations: those forming before/during the last major merger, and those accreted after the last major merger. The former population are more centrally located at z=0, while stars forming in low mass late forming proto-galaxies are spread through the halo. A difference in observed binding energy distinguishes these two populations. We look at possible chemical abundance signatures of the two populations. We also show that galaxies forming in isolated low sigma peaks will form from primordial material. Thus, even though the oldest stars are centrally concentrated as they originated in the early collapsing, densest regions, primordial stars would be found distributed throughout the halo. Thus, the lack of observed metal free stars can be taken as directly constraining the Population III IMF, and the lowest metallicity observed stars can be interpreted as holding clues to the chemical yields of Pop III stars.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Increasingly large populations of disk galaxies are now being observed at increasingly high redshifts, providing new constraints on our knowledge of how such galaxies evolve. Are these observations consistent with a cosmology in which structures form hierarchically? To probe this question, we employ SPH/N-body galaxy scale simulations of late-type galaxies. We examine the evolution of these simulated disk galaxies from redshift 1 to 0, looking at the mass-size and luminosity-size relations, and the thickness parameter, defined as the ratio of scale-height to scale-length. The structural parameters of our simulated disks settle down quickly, and after redshift z=1 the galaxies evolve to become only slightly flatter. Our present day simulated galaxies are larger, more massive, less bright, and redder than at z=1. The inside-out nature of the growth of our simulated galaxies reduces, and perhaps eliminates, expectations of evolution in the size-mass relation.

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.

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.

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.

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.

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<z<8, we show that the effect of lensing can become significant. Hence, if a population of very-high-z supernovae was ever discovered, it would be crucial to fully understand the effect of lensing, before these SNe could be used to constrain cosmological models.

Despite having entered the era of "precision cosmology," the formation of galaxies within the favored CDM cosmological paradigm remains problematic. By relating our N-body/smooth particle hydro-dynamical simulation to an extensive range of Galactic and extragalactic observations, we shed light on the formation and evolution of the Milky Way and other late type galaxies. In light of recent observations of the stellar populations of extragalactic thick disks, we examine the proposal that the thick disk formed in a high redshift period characterized by gas rich merging. We show such a scenario to be consistent with color observations. We then follow the evolution of structural parameters of the subsequently formed disk galaxy, from redshift ~1 to the present. Consistent with observation, little evolution in the ratio of scale-height to scale-length is found in our simulated galaxy, despite its somewhat chaotic origins. The simulated galaxy in this report forms part of a larger study of a suite of galaxies, with which these issues are to be studied in detail.

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.

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.

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)

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.

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.

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.

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.

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)

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.

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.

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.

Using an analytical model, we compute the distribution of image separations resulting from gravitational lensing of distant sources, for 7 COBE-normalized CDM models with various combinations of Omega_0 and lambda_0. Our model assumes that multiple imaging results from strong lensing by individual galaxies. We model galaxies as nonsingular isothermal spheres, and take into account the finite angular size of the sources. Our model neglects the contribution of the background matter distribution, and assumes that lensing is entirely caused by galaxies. To test the validity of this assumption, we performed a series of ray-tracing experiments to study the effect of the background matter on the distribution of image separations. The analytical model predicts that the distributions of image separations are virtually indistinguishable for flat, cosmological constant models with different values of Omega_0. For models with no cosmological constant, the distributions of image separations do depend upon Omega_0, but this dependence is weak. We conclude that while the number of multiple-imaged sources can put strong constraints on the cosmological parameters, the distribution of image separations does not constrain the cosmological models in any significant way, and mostly provides constraints on the structure of the galaxies responsible for lensing.

We summarize the ideas that led to the Adaptive Smoothed Particle Hydrodynamics (ASPH) algorithm, with anisotropic smoothing and shock-tracking. We then identify a serious new problem for SPH simulations with shocks and radiative cooling --- false cooling --- and discuss a possible solution based on the shock-tracking ability of ASPH.

We consider an explosion at the center of a halo which forms at the intersection of filaments inside a cosmological pancake, a convenient test-bed model for galaxy formation. ASPH/P3M simulations reveal that such explosions are anisotropic. The energy and metals are channeled into the low density regions, away from the pancake. The pancake remains essentially undisturbed, even if the explosion is strong enough to blow away all the gas located inside the halo and reheat the IGM surrounding the pancake. Infall quickly replenishes this ejected gas and gradually restores the gas fraction as the halo continues to grow. Estimates of the collapse epoch and SN energy-release for galaxies of different mass in the CDM model can relate these results to scale-dependent questions of blow-out and blow-away and their implication for early IGM heating and metal enrichment and the creation of gas-poor dwarf galaxies.

Adaptive SPH and N-body simulations were carried out to study the evolution of the equilibrium structure of dark matter halos that result from the gravitational instability and fragmentation of cosmological pancakes. Such halos resemble those formed by hierarchical clustering from realistic initial conditions in a CDM universe and, therefore, serve as a test-bed model for studying halo dynamics. The dark matter density profile is close to the universal halo profile identified previously from N-body simulations of structure formation in CDM, with a total mass and concentration parameter which grow linearly with scale factor a. When gas is included, this concentration parameter is slightly larger than the pure N-body result. We also find that the dark matter velocity distribution is less isotropic and more radial than found by N-body simulations of CDM.

Using a multiple-lens plane algorithm, we study light propagation in inhomogeneous universes, for 43 different COBE-normalized Cold Dark Matter models, with various values of the density parameter Omega_0, cosmological constant lambda_0, Hubble constant H_0, and rms density fluctuation sigma_8. This is the largest cosmological parameter survey ever done in this field. We performed a total of 3,798 experiments, each experiment consisting of propagating a square beam of angular size 21.9'' x 21.9'' composed of 116,281 light rays from the observer up to redshift z=3. These experiments provide statistics of the magnification, shear, and multiple imaging of distant sources. The results of these experiments can be compared with observations, and eventually help constraining the possible values of the cosmological parameters. Additionally, they provide insight into the gravitational lensing process and its complex relationship with the various cosmological parameters.

We introduce a set of four new publicly available N-body simulations, the most recent additions to the Texas P^3M Database. Our models probe the less studied parameter space region of moderate volume (100/h Mpc box) combined with fine mass resolution (~10^12 M_Sol, roughly comparable to a L_* galaxy), making these simulations especially suitable for study of major large-scale structure features such as voids, and for comparison with the largest three-dimensional redshift surveys currently available. Our cosmological models (LCDM, TOCDM, OCDM, TCDM) are all COBE-normalized, and when possible (LCDM and TOCDM) also cluster-normalized, based on the X-ray cluster M--T relation. The COBE- and cluster-normalized LCDM model reiterates the attractiveness of this currently favored model which does not require the introduction of tilt in order to fit the constraints imposed by observations of other cosmological parameters.

As an idealized model of the effects of energy release by supernovae during galaxy formation, we consider an explosion at the center of a halo which forms at the intersection of filaments in the plane of a cosmological pancake by gravitational instability during pancake collapse. Such halos resemble the virialized objects found in N-body simulations in a CDM universe and, therefore, serve as a convenient, scale-free test-bed model for galaxy formation. ASPH/P3M simulations reveal that such explosions are anisotropic. The energy and metals are channeled into the low density regions, away from the pancake plane. The pancake remains essentially undisturbed, even if the explosion is strong enough to blow away all the gas located inside the halo at the onset of the explosion and reheat the IGM surrounding the pancake. Infall quickly replenishes this ejected gas and gradually restores the gas fraction as the halo mass continues to grow. Estimates of the collapse epoch and SN energy-release for galaxies of different mass in the CDM model can relate these results to scale-dependent questions of blow-out and blow-away and their implication for early IGM heating and metal enrichment and the creation of dark-matter-dominated dwarf galaxies.

Adaptive SPH and N-body simulations were carried out to study the effect of gasdynamics on the structure of dark matter halos that result from the gravitational instability and fragmentation of cosmological pancakes. Such halos resemble those formed in a hierarchically clustering CDM universe and serve as a test-bed model for studying halo dynamics. With no gas, the density profile is close to the universal profile identified previously from N-body simulations of structure formation in CDM. When gas is included, the gas in the halo is approximately isothermal, and both the dark matter and the gas have singular central density profiles which are steeper than that of the dark matter with no gas. This worsens the disagreement between observations of constant density cores in cosmological halos and the singular ones found in simulations. We also find that the dark matter velocity distribution is less isotropic than found by N-body simulations of CDM, because of the strongly filamentary substructure.

When density fluctuations collapse gravitationally out of the expanding cosmological background universe to form galaxies, the secondary energy release which results can affect their subsequent evolution profoundly. We focus here on the effects of one form of such energy release - explosions, such as might result from the supernovae which end the lives of the first generation of massive stars to form inside protogalaxies. We are particularly interested in the consequences of the nonspherical geometry and continuous infall which are characteristic of galaxy formation from realistic initial and boundary conditions. As an idealized model which serves to illustrate and quantify the importance of these effects, we study the effect of explosions on the quasi-spherical objects which form at the intersections of filaments in the plane of a cosmological pancake, as a result of gravitational instability and fragmentation of the pancake. We study the formation and evolution of these "galaxies," subject to the explosive injection of energy at their centers, by numerical gas dynamical simulation in 3D utilizing our new, anisotropic version of Smoothed Particle Hydrodynamics, Adaptive SPH ("ASPH"), with a P3M gravity solver.

We calculate the maximum fraction of matter which is able to condense out of the expanding background universe, for any universe which will expand forever. We use a simple spherical model for the growth of density fluctuations in the universe. This model includes open, matter-dominated universe and universes in which there is an uniform background component (e.g. the cosmological constant or "quintessence"). In these background universes, Omega eventually drops significantly below unity. When this happens, gravitational instability is suppressed and, with it, so is the growth of the collapsed fraction. We identify a limitation of the Press-Schechter (PS) approximation. In this approximation, the mass function determined from the predicted collapse of overdense regions is multiplied by a correction factor of 2, assuming that each bound fluctuation will accrete an equal share of mass from nearby underdense regions. Our model determines the actual value of the correction factor. We show that, while the factor of 2 adopted by the PS approximation is correct for an Einstein-de Sitter universe, it is not correct when the freeze-out of fluctuation growth inherent in the more general class of background universes described above occurs. When freeze-out occurs, the correction factor reduces to unity and the PS approximation must overestimate the collapsed fraction. We apply our model to open CDM and flat Lambda-CDM models. For H0=70 and Omega0=0.3, these models yield asymptotic collapsed fractions of 0.0361 and 0.0562, respectively, on the galaxy cluster mass-scale, only 55% of the values determined by the PS approximation. These results have implications for the use of the latter approximation to compare the observed space density of X-ray clusters today with that predicted by cosmological models.

We have simulated the formation and evolution of large-scale structure in the universe, for 68 different COBE-normalized cosmological models. For each cosmological model, we have performed between 1 and 3 simulations, for a total of 160 simulations. This constitutes the largest database of cosmological simulations ever assembled, and the largest cosmological parameter space ever covered by such simulations. We are making this database available to the astronomical community. We provide instructions for accessing the database and for converting the data from computational units to physical units. The database includes Tilted Cold Dark Matter (TCDM) models, Tilted Open Cold Dark Matter (TOCDM) models, and Tilted Lambda Cold Dark Matter (TLCDM) models. (For several simulations, the primordial exponent n of the power spectrum is near unity, hence these simulations can be considered as "untilted.") The simulations cover a 4-dimensional cosmological parameter phase space, the parameters being the present density parameter Omega_0, cosmological constant lambda_0, and Hubble constant H_0, and the rms density fluctuation sigma_8 at scale 8/h Mpc. All simulations were performed using a P3M algorithm with 64^3 particles on a 128^3 mesh, in a cubic volume of comoving size 128 Mpc. Each simulation starts at a redshift of 24, and is carried up to the present. More simulations will be added to the database in the future.

We compute statistical properties of weak gravitational lensing by large-scale structure in three Cold Dark Matter models. We use a P$^3$M $N$-body code to simulate the formation and evolution of large-scale structure in the universe. We perform $1.1\times10^7$ ray-tracing experiments for each model using the multiple lens-plane algorithm. From the results of these experiments, we calculate the probability distribution functions (PDF) of the convergences, shears, and magnifications, and their root-mean-square (rms) values. We find that the rms values of the convergence and shear agree with the predictions of a nonlinear analytical model. We also find that the PDFs of the magnifications $\mu$ have a peak at values slightly smaller than $\mu=1$, and are strongly skewed toward large magnifications. In particular, for the high-density model, a power-law tail appears in the magnification distribution at large magnifications for sources at redshifts $z_s>2$. The rms values of the magnifications essentially agree with the nonlinear analytical predictions for sources at low redshift, but exceed these predictions for high redshift sources, once the power-law tail appears. We study the effect of magnification bias on the luminosity functions of high-redshift quasars, using the calculated PDFs of the magnifications. We show that the magnification bias is moderate in the absence of the power-law tail in the magnification distribution, but depends strongly on the value of the density parameter. In presence of the power-law tail, the bias becomes considerable, especially at the bright end of the luminosity functions.

Using a multi-plane lensing method that we have developed, we follow the evolution of light beams as they propagate through inhomogeneous universes. We use a P3M code to simulate the formation and evolution of large-scale structure. The resolution of the simulations is increased to sub-Megaparsec scales by using a Monte Carlo method to locate galaxies inside the computational volume according to the underlying particle distribution. We consider cold dark matter models normalized to COBE, and perform a large parameter survey by varying the cosmological parameters Omega_0, lambda_0, H_0, and n (the tilt of the primordial power spectrum).

Using a multi-plane lensing method that we have developed, we follow the evolution of light beams as they propagate through inhomogeneous universes. We use a P3M code to simulate the formation and evolution of large-scale structure. The resolution of the simulations is increased to sub-Megaparsec scales by using a Monte Carlo method to locate galaxies inside the computational volume according to the underlying particle distribution. The galaxies are approximated by isothermal spheres, with each morphological type having its own distribution of masses and core radii. The morphological types are chosen in order to reproduce the observed morphology-density relation. This algorithm has an effective resolution of 9 orders of magnitudes in length, from the size of superclusters down to the core radii of the smallest galaxies. We consider cold dark matter models normalized to COBE, and perform a large parameter survey by varying the cosmological parameters Omega_0, lambda_0, H_0, and n (the tilt of the primordial power spectrum). The values of n are chosen by imposing particular values or sigma_8, the rms mass fluctuation at a scale of 8/h Mpc. We use the power spectrum given by Bunn & White. This is the largest parameter survey ever done is this field.

This paper presents an alternative formulation of the ASPH algorithm for evolving anisotropic smoothing kernels, in which the geometric approach of Shapiro et al. (1996; Paper I) is replaced by an approach involving a local transformation of coordinates to those in which the underlying anisotropic volume changes appear to be isotropic. The ASPH method is presented in 2D and 3D, including a number of details not previously included in Paper I, some of which represent either advances or different choices with respect to Paper I. Among the advances included here are an asynchronous time-integration scheme with different time steps for different particles and the generalization of the ASPH method to 3D. The shock-tracking algorithm described in Paper I for locally adapting the artificial viscosity to restrict viscous heating just to particles encountering shocks, is not included here. Instead, we adopt a different interpolation kernel for use with the artificial viscosity, which has the effect of spatially localizing effects of the artificial viscosity. This version of the ASPH method in 2D and 3D is then applied to a series of 1D, 2D, and 3D test problems, and the results are compared to those of standard SPH applied to the same problems. These include the problem of cosmological pancake collapse, the Riemann shock tube, cylindrical and spherical Sedov blast waves, the collision of two strong shocks, and problems involving shearing disks intended to test the angular momentum conservation properties of the method. These results further support the idea that ASPH has significantly better resolving power than standard SPH for a wide range of problems, including that of cosmological structure formation. (Abridged)

We present a set of cosmological variables, called "supercomoving variables," which are particularly useful for describing the gas dynamics of cosmic structure formation. For ideal gas with gamma=5/3, the supercomoving position, velocity, density, temperature, and pressure are constant in time in a uniform, isotropic, adiabatically expanding universe. Expressed in terms of these supercomoving variables, the cosmological fluid conservation equations and the Poisson equation closely resemble their noncosmological counterparts. This makes it possible to generalize noncosmological results and techniques to cosmological problems, for a wide range of cosmological models. These variables were initially introduced by Shandarin for matter-dominated models only. We generalize supercomoving variables to models with a uniform component corresponding to a nonzero cosmological constant, domain walls, cosmic strings, a nonclumping form of nonrelativistic matter (e.g. massive nettrinos), or radiation. Each model is characterized by the value of the density parameter Omega0 of the nonrelativistic matter component in which density fluctuation is possible, and the density parameter OmegaX of the additional, nonclumping component. For each type of nonclumping background, we identify FAMILIES within which different values of Omega0 and OmegaX lead to fluid equations and solutions in supercomoving variables which are independent of Omega0 and OmegaX. We also include the effects of heating, radiative cooling, thermal conduction, viscosity, and magnetic fields. As an illustration, we describe 3 familiar cosmological problems in supercomoving variables: the growth of linear density fluctuations, the nonlinear collapse of a 1D plane-wave density fluctuation leading to pancake formation, and the Zel'dovich approximation.

We simulate the growth of large-scale structure in the universe using a N-body code. By combining these simulations with a Monte-Carlo method, we generate galaxy distributions at present that reproduces the observed morphology-density relation, with most ellipticals concentrated in the densest regions. By "tying" each galaxy to its nearest particle, we trace the trajectory of that galaxy back in time. This enables us to reconstruct the distribution of galaxies at high redshift, and the trajectory of each galaxy from its formation epoch to the present. Our goal is to determine whether the morphological type of galaxies is primarily determined by the initial conditions in which these galaxies form, or by evolutionary processes occurring later. We compare the environment in which galaxies are at the epoch of galaxy formation (z=3) and at the present. Making the null hypothesis that morphological types do not evolve, we compare the galaxies that form in low density environments but end up later in high density environments to the ones that form also in low density environment but remain in low density environment. The first group contains a larger proportion of E and S0 galaxies than the second group. We assume that the galaxy formation process cannot distinguish a low density environment that will always remain low density from one that will eventually become high density. Hence, these results force us to discard the null hypothesis of no morphological evolution. Our study suggests that 75% of the E and S0 galaxies observed at present formed as such, while the remaining 25% formed as spirals, and underwent morphological evolution. We conclude that most galaxies did not undergo morphological evolution, but a non-negligible fraction did.

We describe a numerical algorithm which simulates the propagation of light in inhomogeneous universes, using the multiple lens-plane method. The deformation and deflection of light beams as they interact with each lens plane are computed using the filled-beam approximation. We use a N-body code to simulate large-scale structure formation in the universe. We extend the length resolution of the simulations to sub-Megaparsec scales by using a Monte-Carlo method for locating galaxies according to the underlying distribution of background matter. The observed galaxy 2-point correlation function and morphology-density relation are reproduced. This algorithm constitutes a major improvement over previous methods. We test this algorithm for three different CDM models: an Einstein-de Sitter model, an open model, and a flat model with a nonzero cosmological constant. In each case, we compute the shear and magnification caused by the presence of inhomogeneities. Our results are the following: (1) The magnification is totally dominated by the convergence; the shear contributes less than one part in 10^4. (2) Most of the cumulative shear and magnification is contributed by matter located at intermediate redshifts z=1-2. (3) The redshift where the largest contribution to shear and magnification occurs is model-dependent. (4) The bulk of the shear and magnification does not originate from direct hits with galaxies, but rather from the tidal influence of nearby and more distant galaxies, and background matter. (5) The average contributions of background matter and nearby galaxies to the shear is comparable for models with small Omega0. For the Einstein-de Sitter model, the contribution of the background matter exceeds the one of nearby galaxies by nearly one order of magnitude.

In theories in which the cosmological constant Lambda takes a variety of values in different ``subuniverses,'' the probability distribution of its observed values is conditioned by the requirement that there must be someone to measure it. This probability is proportional to the fraction of matter which is destined to condense out of the background into mass concentrations large enough to form observers. We calculate this ``collapsed fraction'' by a simple, pressure-free, spherically symmetric, nonlinear model for the growth of density fluctuations in a flat universe with arbitrary value of the cosmological constant, applied in a statistical way to the observed spectrum of density fluctuations at recombination. From this, the probability distribution for the vacuum energy density rho_V=Lambda/8pi G for Gaussian random density fluctuations is derived analytically. It is shown that the results depend on only one quantity, sigma^3 RHO, where sigma^2 and RHO are the variance and mean value of the fluctuating matter density field at recombination, respectively. To calculate sigma, we adopt the flat CDM model with nonzero cosmological constant and fix the amplitude and shape of the primordial power spectrum in accordance with data on cosmic microwave background anisotropy from the COBE satellite DMR experiment. A comparison of the results of this calculation of the likely values of rho_V with present observational bounds on the cosmological constant indicates that the small, positive value of rho_V (up to 3 times greater than the present cosmic mass density) suggested recently by several lines of evidence is a reasonably likely value to observe, even if all values of rho_V are equally likely a priori.

The gravitational instability of cosmological pancakes composed of collisionless dark matter in an Einstein-de Sitter universe is investigated numerically to demonstrate that pancakes are unstable with respect to fragmentation and the formation of filaments. A ``pancake'' is defined here as the nonlinear outcome of the growth of a 1D, sinusoidal, plane-wave, adiabatic density perturbation. We have used high resolution, 2D, N-body simulations by the Particle-Mesh (PM) method to study the response of pancakes to perturbation by either symmetric (density) or antisymmetric (bending or rippling) modes, with corresponding wavevectors k_s and k_a transverse to the wavevector k_p of the unperturbed pancake plane-wave. We consider dark matter which is initially ``cold'' (i.e. with no random thermal velocity in the initial conditions). We also investigate the effect of a finite, random, isotropic, initial velocity dispersion (i.e. initial thermal velocity) on the fate of pancake collapse and instability. Pancakes are shown to be gravitationally unstable with respect to all perturbations of wavelength l<l_p (where l_p= 2pi/k_p). These results are in contradiction with the expectations of an approximate, thin-sheet energy argument.