Showing 1–50 of 78 results for author: Sommer-Larsen, J

Chasing Lyman alpha-emitting galaxies at z = 8.8

Authors: Peter Laursen, Jesper Sommer-Larsen, Bo Milvang-Jensen, Johan P. U. Fynbo, Alexei O. Razoumov

Abstract: With a total integration time of 168 hours and a narrowband (NB) filter tuned to Lyman alpha at z = 8.8, the UltraVISTA survey has set out to find some of the most distant galaxies, on the verge of the Epoch of Reionization. Previous calculations of the expected number of detected Lya-emitting galaxies (LAEs) at this redshift did not explicitly take into account the radiative transfer (RT) of Lya.… ▽ More With a total integration time of 168 hours and a narrowband (NB) filter tuned to Lyman alpha at z = 8.8, the UltraVISTA survey has set out to find some of the most distant galaxies, on the verge of the Epoch of Reionization. Previous calculations of the expected number of detected Lya-emitting galaxies (LAEs) at this redshift did not explicitly take into account the radiative transfer (RT) of Lya. In this work we combine a theoretical model for the halo mass function with numerical results from high-res cosmological hydrosimulations with LyC+Lya RT, assessing the visibility of LAEs residing in these halos. Uncertainties such as cosmic variance and the anisotropic escape of Lya are taken into account, and it is predicted that once the survey has finished, the probabilities of detecting none, one, or more than one are ~90%, ~10%, and ~1%; a significantly smaller success rate compared to earlier predictions, due to the combined effect of a highly neutral IGM scattering Lya to such large distances from the galaxy that they fall outside the observational aperture, and to the actual depth of the survey being less than predicted. Because the IGM affects NB and broadband (BB) magnitudes differently, we argue for a relaxed color selection criterion of NB - BB ~ +0.85. But since the flux is continuum-dominated, even if a galaxy is detectable in the NB its probability of being selected as a NB excess object is <~35%. Various properties of galaxies at this redshift are predicted, e.g. UV and Lya LFs, M*-Mh relation, spectral shape, optimal aperture, and the anisotropic escape of Lya through both a dusty ISM and a partly neutral IGM. Finally, we describe and publish a fast numerical code for adding numbers with asymmetric uncertainties ("x_{-sigma_1}^{+sigma_2}") proving to be significantly better than the standard, but wrong, way of adding upper and lower uncertainties in quadrature separately. △ Less

Submitted 19 June, 2018; originally announced June 2018.

Comments: Submitted to A&A, comments are welcome

Journal ref: A&A 627, A84 (2019)

The high-redshift evolution of the Red Sequence scatter from joint simulations and HAWK-I Cluster Survey

Authors: A. D. Romeo, P. Cerulo, Kang Xi, E. Contini, J. Sommer-Larsen, I. Gavignaud

Abstract: We study the evolution of the Red Sequence (RS) scatter in galaxy clusters and groups simultaneously using predictions from our simulations (cosmological hydrodynamic + semi-analytical) as well as observational data from the HAWK-I Cluster Survey (HCS), a sample of galaxy clusters at redshifts $0.8 < z < 1.5$. We analyze the intrinsic scatter of the RS to investigate whether the stellar age can be… ▽ More We study the evolution of the Red Sequence (RS) scatter in galaxy clusters and groups simultaneously using predictions from our simulations (cosmological hydrodynamic + semi-analytical) as well as observational data from the HAWK-I Cluster Survey (HCS), a sample of galaxy clusters at redshifts $0.8 < z < 1.5$. We analyze the intrinsic scatter of the RS to investigate whether the stellar age can be its main effective driver, at the same time assessing the role of metallicity variations in shaping the RS building at around epoch $z\sim$1 and beyond. To this purpose we rely on various methods to derive the average age and age spread from the RS colour scatter of the HCS sample, with the aid of population synthesis models. The RS scatter predicted by the models at $z< 0.7$ is found to not depend on the star formation history adopted, whilst at $z>the correlation found between age and rest-frame colour scatters is quite robust, although all age scatter estimations ultimately depend on the definition of RS as well as on the completeness limits adopted. We find that the age spread of RS galaxies predicted by both hydrodynamical simulations and SAM increases with cosmic epoch, while the ratio between the age spread and the average age remains approximately constant. Both trends are in agreement with observational results from both the HCS and other literature samples. △ Less

Submitted 14 November, 2016; originally announced November 2016.

Comments: submitted to ApJ

The galaxy counterpart of the high-metallicity and 16 kpc impact parameter DLA towards Q0918+1636 - a challenge to galaxy formation models?

Authors: J. Sommer-Larsen, J. P. U. Fynbo

Abstract: The quasar Q0918+1636 (z=3.07) has an intervening high-metallicity Damped Lyman-alpha Absorber (DLA) along the line of sight, at a redshift of z=2.58. The DLA is located at a large impact parameter of 16.2 kpc, and has an almost solar metallicity. It is shown, that a novel type of cosmological galaxy formation models, invoking a new SNII feedback prescription, the Haardt & Madau (2012) UVB field a… ▽ More The quasar Q0918+1636 (z=3.07) has an intervening high-metallicity Damped Lyman-alpha Absorber (DLA) along the line of sight, at a redshift of z=2.58. The DLA is located at a large impact parameter of 16.2 kpc, and has an almost solar metallicity. It is shown, that a novel type of cosmological galaxy formation models, invoking a new SNII feedback prescription, the Haardt & Madau (2012) UVB field and explicit treatment of UVB self-shielding, can reproduce the observed characteristics of the DLA. UV radiation from young stellar populations in the galaxy, in particular in the photon energy range 10.36-13.61 eV (relating to Sulfur II abundance), are also considered in the analysis. It is found that a) for L~L* galaxies (at z=2.58), about 10% of the sight-lines through the galaxies at impact parameter 16.2 kpc will display a Sulfur II column density N(SII)$>$ 10$^{15.82}$ cm$^{-2}$ (the observed value for the DLA), and b) considering only cases where a near-solar metallicity will be detected at 16.2 kpc impact parameter, the probability distribution of galaxy SFR peaks near the value observed for the DLA galaxy counterpart of ~27 Msun/yr. It is argued, that the bulk of the alpha-elements, like Sulfur, traced by the high metal column density, b=16.2 kpc absorption lines, were produced by evolving young stars in the inner galaxy, and later transported outward by galactic winds. △ Less

Submitted 25 October, 2016; originally announced October 2016.

Comments: 22 pages, 24 figures, MNRAS in press

Journal ref: MNRAS, 464, 2441 (2017)

SImulator of GAlaxy Millimetre/submillimetre Emission (SIGAME): CO emission from massive z=2 main-sequence galaxies

Authors: Karen P. Olsen, Thomas R. Greve, Christian Brinch, Jesper Sommer-Larsen, Jesper Rasmussen, Sune Toft, Andrew Zirm

Abstract: We present SIGAME (SImulator of GAlaxy Millimetre/submillimetre Emission), a new numerical code designed to simulate the 12CO rotational line emission spectrum of galaxies. Using sub-grid physics recipes to post-process the outputs of smoothed particle hydrodynamics (SPH) simulations, a molecular gas phase is condensed out of the hot and partly ionized SPH gas. The gas is subjected to far-UV radia… ▽ More We present SIGAME (SImulator of GAlaxy Millimetre/submillimetre Emission), a new numerical code designed to simulate the 12CO rotational line emission spectrum of galaxies. Using sub-grid physics recipes to post-process the outputs of smoothed particle hydrodynamics (SPH) simulations, a molecular gas phase is condensed out of the hot and partly ionized SPH gas. The gas is subjected to far-UV radiation fields and cosmic ray ionization rates which are set to scale with the local star formation rate volume density. Level populations and radiative transport of the CO lines are solved with the 3-D radiative transfer code LIME. We have applied SIGAME to cosmological SPH simulations of three disc galaxies at z=2 with stellar masses in the range ~(0.5-2)x10^11 Msun and star formation rates ~40-140 Msun/yr. Global CO luminosities and line ratios are in agreement with observations of disc galaxies at z~2 up to and including J=3-2 but falling short of the few existing J=5-4 observations. The central 5 kpc regions of our galaxies have CO 3-2/1-0 and 7-6/1-0 brightness temperature ratios of ~0.55-0.65 and ~0.02-0.08, respectively, while further out in the disc the ratios drop to more quiescent values of ~0.5 and <0.01. Global CO-to-H2 conversion (alpha_CO) factors are ~=1.5 Msun*pc^2/(K km s/1), i.e. ~2-3 times below typically adopted values for disc galaxies, and alpha_CO increases with radius, in agreement with observations of nearby galaxies. Adopting a top-heavy Giant Molecular Cloud (GMC) mass spectrum does not significantly change the results. Steepening the GMC density profile leads to higher global line ratios for J_up>=3 and CO-to-H2 conversion factors [~=3.6 Msun*pc^2/(K km/s)]. △ Less

Submitted 26 February, 2016; v1 submitted 30 June, 2015; originally announced July 2015.

Comments: 28 pages, 20 figures. Accepted for Publication in MNRAS. Substantial revisions from the previous version, including tests with model galaxies similar to the Milky Way. Improved figures and added tables

Journal ref: MNRAS 457, 3306-3333 (2016)

A study on the multicolour evolution of Red Sequence galaxy populations: insights from hydrodynamical simulations and semi-analytical models

Authors: A. D. Romeo, Xi Kang, E. Contini, J. Sommer-Larsen, R. Fassbender, N. R. Napolitano, V. Antonuccio-Delogu, I. Gavignaud

Abstract: By means of our own cosmological-hydrodynamical simulation and semi-analytical model we studied galaxy population properties in clusters and groups, spanning over 10 different bands from UV to NIR, and their evolution since redshift z=2. We compare our results in terms of galaxy red/blue fractions and luminous-to-faint ratio (LFR) on the Red Sequence (RS) with recent observational data reaching be… ▽ More By means of our own cosmological-hydrodynamical simulation and semi-analytical model we studied galaxy population properties in clusters and groups, spanning over 10 different bands from UV to NIR, and their evolution since redshift z=2. We compare our results in terms of galaxy red/blue fractions and luminous-to-faint ratio (LFR) on the Red Sequence (RS) with recent observational data reaching beyond z=1.5. Different selection criteria were tested in order to retrieve galaxies belonging to the RS: either by their quiescence degree measured from their specific SFR ("Dead Sequence"), or by their position in a colour-colour plane which is also a function of sSFR. In both cases, the colour cut and the limiting magnitude threshold were let evolving with redshift, in order to follow the natural shift of the characteristic luminosity in the LF. We find that the Butcher-Oemler effect is wavelength-dependent, with the fraction of blue galaxies increasing steeper in optical colours than in NIR. Besides, only when applying a lower limit in terms of fixed absolute magnitude, a steep BO effect can be reproduced, while the blue fraction results less evolving when selecting samples by stellar mass or an evolving magnitude limit. We then find that also the RS-LFR behaviour, highly debated in the literature, is strongly dependent on the galaxy selection function: in particular its very mild evolution recovered when measured in terms of stellar mass, is in agreement with values reported for some of the highest redshift confirmed (proto)clusters. As to differences through environments, we find that normal groups and (to a lesser extent) cluster outskirts present the highest values of both star forming fraction and LFR at low z, while fossil groups and cluster cores the lowest: this separation among groups begins after z~0.5, while earlier all group star forming properties are undistinguishable. △ Less

Submitted 30 June, 2015; v1 submitted 22 May, 2015; originally announced May 2015.

Comments: revised version, A&A accepted (11 pages, 6 figures)

Journal ref: A&A 581, A50 (2015)

Galactic winds - How to launch galactic outflows in typical Lyman-break galaxies

Authors: Wolfgang von Glasow, Martin Krause, Jesper Sommer-Larsen, Andreas Burkert

Abstract: We perform hydrodynamical simulations of a young galactic disc embedded in a hot gaseous halo using parameters typical for Lyman break galaxies (LBGs). We take into account the (static) gravitational potentials due to a dark matter halo, a stellar bulge and a disc of stars and gas. Star formation is treated by a local Kennicutt-Schmidt law. We simplify the structure of the interstellar medium by r… ▽ More We perform hydrodynamical simulations of a young galactic disc embedded in a hot gaseous halo using parameters typical for Lyman break galaxies (LBGs). We take into account the (static) gravitational potentials due to a dark matter halo, a stellar bulge and a disc of stars and gas. Star formation is treated by a local Kennicutt-Schmidt law. We simplify the structure of the interstellar medium by restricting the computational domain to a 25th of the full azimuthal angle, effectively assuming large-scale axisymmetry and neglecting any effects of spiral structure, and focus on the large-scale ISM drivers, the superbubbles. Supernovae are triggered randomly and have preset event sizes of several tens to hundreds. We further investigate different halo gas pressures and energy injection methods. Many of our simulated galaxies, but not all, develop bipolar outflows. We characterise the strength of the outflow by mass and energy outflow rates, and investigate the effect of changes to the details of the model. We find that supernovae are more effective if comprised into larger superbubbles. The weight and the pressure of the halo gas is able to quench galactic outflows. The wind emerges from a series of superbubbles in regions where a critical star formation density is exceeded. The superbubbles expand into the gaseous halo at slightly supersonic speed, producing radiative shock waves with similar characteristics as the absorptions systems observed around LBGs. △ Less

Submitted 13 June, 2013; originally announced June 2013.

Comments: 20 pages, 21 figures, accepted by MNRAS

Evolution of the Mass-Metallicity relations in passive and star-forming galaxies from SPH-cosmological simulations

Authors: A. D. Romeo Velonà, J. Sommer-Larsen, N. R. Napolitano, V. Antonuccio-Delogu, S. Cielo, I. Gavignaud

Abstract: We present results from SPH-cosmological simulations, including self-consistent modelling of SN feedback and chemical evolution, of galaxies belonging to two clusters and twelve groups. We reproduce the mass-metallicity (ZM) relation of galaxies classified in two samples according to their star-forming activity, as parametrized by their sSFR, across a redshift range up to z=2. Its slope shows ir… ▽ More We present results from SPH-cosmological simulations, including self-consistent modelling of SN feedback and chemical evolution, of galaxies belonging to two clusters and twelve groups. We reproduce the mass-metallicity (ZM) relation of galaxies classified in two samples according to their star-forming activity, as parametrized by their sSFR, across a redshift range up to z=2. Its slope shows irrelevant evolution in the passive sample, being steeper in groups than in clusters. However, the sub-sample of high-mass passive galaxies only is characterized by a steep increase of the slope with redshift, from which it can be inferred that the bulk of the slope evolution of the ZM relation is driven by the more massive passive objects. (...ABRIDGED...) The ZM relation for the star-forming sample reveals an increasing scatter with redshift, indicating that it is still being built at early epochs. The star-forming galaxies make up a tight sequence in the SFR-M_* plane at high redshift, whose scatter increases with time alongside with the consolidation of the passive sequence. We also confirm the anti-correlation between sSFR and stellar mass, pointing at a key role of the former in determining the galaxy downsizing, as the most significant means of diagnostics of the star formation efficiency. Likewise, an anti-correlation between sSFR and metallicity can be established for the star-forming galaxies, while on the contrary more active galaxies in terms of simple SFR are also metal-richer. We discuss these results in terms of the mechanisms driving the evolution within the high- and low-mass regimes at different epochs: mergers, feedback-driven outflows and the intrinsic variation of the star formation efficiency. △ Less

Submitted 4 September, 2014; v1 submitted 7 May, 2013; originally announced May 2013.

Comments: Emended list of authors

Journal ref: ApJ 770, 155 (2014)

On the Use of Ly-alpha Emitters as Probes of Reionization

Authors: Hannes Jensen, Peter Laursen, Garrelt Mellema, Ilian T. Iliev, Jesper Sommer-Larsen, Paul R. Shapiro

Abstract: We use numerical simulations to study the effects of the patchiness of a partly reionized intergalactic medium (IGM) on the observability of Ly-alpha emitters (LAEs) at high redshifts (z ~ 6). We present a new model that divides the Ly-alpha radiative transfer into a (circum-)galactic and an extragalactic (IGM) part, and investigate how the choice of intrinsic line model affects the IGM transmissi… ▽ More We use numerical simulations to study the effects of the patchiness of a partly reionized intergalactic medium (IGM) on the observability of Ly-alpha emitters (LAEs) at high redshifts (z ~ 6). We present a new model that divides the Ly-alpha radiative transfer into a (circum-)galactic and an extragalactic (IGM) part, and investigate how the choice of intrinsic line model affects the IGM transmission results. We use our model to study the impact of neutral hydrogen on statistical observables such as the Ly-alpha restframe equivalent width (REW) distribution, the LAE luminosity function and the two-point correlation function. We find that if the observed changes in LAE luminosity functions and equivalent width distributions between z ~ 6 and z ~ 7 are to be explained by an increased IGM neutral fraction alone, we require an extremely late and rapid reionization scenario, where the Universe was ~ 40 % ionized at z = 7, ~ 50 % ionized at z = 6.5 and ~ 100 % ionized at z = 6. This is in conflict with other observations, suggesting that intrinsic LAE evolution at z > 6 cannot be completely neglected. We show how the two-point correlation function can provide more robust constraints once future observations obtain larger LAE samples, and provide predictions for the sample sizes needed to tell different reionization scenarios apart. △ Less

Submitted 18 October, 2012; v1 submitted 18 June, 2012; originally announced June 2012.

Comments: Accepted for publication in MNRAS

Stellar population gradients from cosmological simulations: dependence on mass and environment in local galaxies

Authors: C. Tortora, A. D. Romeo, N. R. Napolitano, V. Antonuccio-Delogu, A. Meza, J. Sommer-Larsen, M. Capaccioli

Abstract: The age and metallicity gradients for a sample of group and cluster galaxies from N-body+hydrodynamical simulation are analyzed in terms of galaxy stellar mass. Dwarf galaxies show null age gradient with a tail of high and positive values for systems in groups and cluster outskirts. Massive systems have generally zero age gradients which turn to positive for the most massive ones. Metallicity grad… ▽ More The age and metallicity gradients for a sample of group and cluster galaxies from N-body+hydrodynamical simulation are analyzed in terms of galaxy stellar mass. Dwarf galaxies show null age gradient with a tail of high and positive values for systems in groups and cluster outskirts. Massive systems have generally zero age gradients which turn to positive for the most massive ones. Metallicity gradients are distributed around zero in dwarf galaxies and become more negative with mass; massive galaxies have steeper negative metallicity gradients, but the trend flatten with mass. In particular, fossil groups are characterized by a tighter distribution of both age and metallicity gradients. We find a good agreement with both local observations and independent simulations. The results are also discussed in terms of the central age and metallicity, as well as the total colour, specific star formation and velocity dispersion. △ Less

Submitted 13 September, 2010; originally announced September 2010.

Comments: 9 pages, 5 figures, accepted for publication on MNRAS

Intergalactic Transmission and its Impact on the Lyα Line

Authors: Peter Laursen, Jesper Sommer-Larsen, Alexei O. Razoumov

Abstract: We study the intergalactic transmission of radiation in the vicinity of the Lyα wavelength. Simulating sightlines through the intergalactic medium (IGM) in detailed cosmological hydrosimulations, the impact of the IGM on the shape of the line profile from Lyα emitting galaxies at redshifts 2.5 to 6.5 is investigated. In particular we show that taking into account the correlation of the density and… ▽ More We study the intergalactic transmission of radiation in the vicinity of the Lyα wavelength. Simulating sightlines through the intergalactic medium (IGM) in detailed cosmological hydrosimulations, the impact of the IGM on the shape of the line profile from Lyα emitting galaxies at redshifts 2.5 to 6.5 is investigated. In particular we show that taking into account the correlation of the density and velocity fields of the IGM with the galaxies, the blue part of the spectrum may be appreciably reduced, even at relatively low redshifts. This may in some cases provide an alternative to the often-invoked outflow scenario, although it is concluded that this model is still a plausible explanation of the many asymmetric Lyα profiles observed. Applying the calculated wavelength dependent transmission to simulated spectra from Lyα emitting galaxies, we derive the fraction of photons that are lost in the IGM, in addition to what is absorbed internally in the galaxies due to dust. Moreover, by comparing the calculated transmission of radiation blueward of the Lyα line, the total optical depth to Thomson scattering of cosmic microwave background, with corresponding observations, we are able to constrain the epoch when the Universe was reionized to z <~ 8.5. △ Less

Submitted 6 December, 2010; v1 submitted 7 September, 2010; originally announced September 2010.

Comments: Substantially extended, ~30 references added, 1.5 page extra (article style) in particular on the impact of the IGM at z~5.8 and z~6.5, 2 extra figures, unnecessary fluff cut out, accepted for publication in ApJ

Hot gas haloes around disc galaxies: O VII column densities from galaxy formation simulations

Authors: Evangelia Ntormousi, Jesper Sommer-Larsen

Abstract: Numerical models of disc galaxy formation predict the existence of extended, hot ( T~10^6 K) gas haloes around present day spirals. The X-ray luminosity of these haloes is predicted to increase strongly with galaxy mass. However, searches for their X-ray emission have not been successful so far. We calculate the all sky O VII column density distributions for the haloes of three Milky Way like disc… ▽ More Numerical models of disc galaxy formation predict the existence of extended, hot ( T~10^6 K) gas haloes around present day spirals. The X-ray luminosity of these haloes is predicted to increase strongly with galaxy mass. However, searches for their X-ray emission have not been successful so far. We calculate the all sky O VII column density distributions for the haloes of three Milky Way like disc galaxies, resulting from cosmological high-resolution, N-body/gasdynamical simulations. We perform calculations both including the disc gas and without it, so the disc contribution to the column density is quantified. It is found that the column densities estimated for Milky Way-like galaxies are just below the observational upper limit, making a test of the hot halo paradigm likely within observational reach. △ Less

Submitted 23 July, 2010; originally announced July 2010.

Comments: accepted for publication in MNRAS

Cosmological Simulations of Massive Compact High-z Galaxies

Authors: Jesper Sommer-Larsen, Sune Toft

Abstract: In order to investigate the structure and dynamics of the recently discovered massive (M_* > 10^11 M_sun) compact z~2 galaxies, cosmological hydrodynamical/N-body simulations of a proto-cluster region have been undertaken. At z=2, the highest resolution simulation contains ~5800 resolved galaxies, of which 509, 27 and 5 have M_* > 10^10 M_sun, > 10^11 M_sun and > 4x10^11 M_sun, respectively. Eff… ▽ More In order to investigate the structure and dynamics of the recently discovered massive (M_* > 10^11 M_sun) compact z~2 galaxies, cosmological hydrodynamical/N-body simulations of a proto-cluster region have been undertaken. At z=2, the highest resolution simulation contains ~5800 resolved galaxies, of which 509, 27 and 5 have M_* > 10^10 M_sun, > 10^11 M_sun and > 4x10^11 M_sun, respectively. Effective radii and characteristic stellar densities have been determined for all galaxies. At z=2, for the definitely well resolved mass range of M_* > 10^11 Msun, the mass-size relation is consistent with observational findings for the most compact z~2 galaxies. The very high velocity dispersion recently measured for a compact z~2 galaxy (~510 km/s; van Dokkum et al 2009) can be matched at about the 1-sigma level, although a somewhat larger mass than the estimated M_* ~ 2 x 10^11 M_sun is indicated. For the above mass range, the galaxies have an average axial ratio <b/a> = 0.64 +/- 0.02 with a dispersion of 0.1, an average rotation to 1D velocity dispersion ratio <v/sigma> = 0.46 +/- 0.06 with a dispersion of 0.3, and a maximum value of v/sigma ~ 1.1. Rotation and velocity anisotropy both contribute in flattening the compact galaxies. Some of the observed compact galaxies appear flatter than any of the simulated galaxies. Finally, it is found that the massive compact galaxies are strongly baryon dominated in their inner parts, with typical dark matter mass fractions of order only 20% inside of r=2R_eff. △ Less

Submitted 4 September, 2009; originally announced September 2009.

Comments: 10 pages, 8 figures, submitted to ApJ

Journal ref: Astrophys.J.721:1755-1764,2010

High-Redshift Star-Forming Galaxies: Angular Momentum and Baryon Fraction, Turbulent Pressure Effects and the Origin of Turbulence

Authors: A. Burkert, R. Genzel, N. Bouche, G. Cresci, S. Khochfar, J. Sommer-Larsen, A. Sternberg, T. Naab, N. Foerster-Schreiber, L. Tacconi, K. Shapiro, E. Hicks, D. Lutz, R. Davies, P. Buschkamp, S. Genel

Abstract: The structure of a sample of high-redshift (z=2), rotating galaxies with high star formation rates and turbulent gas velocities of sigma=40-80 km/s is investigated. Fitting the observed disk rotational velocities and radii with a Mo, Mao, White (1998) (MMW) model requires unusually large disk spin parameters lambda_d>0.1 and disk-to-dark halo mass fraction m_d=0.2, close to the cosmic baryon fract… ▽ More The structure of a sample of high-redshift (z=2), rotating galaxies with high star formation rates and turbulent gas velocities of sigma=40-80 km/s is investigated. Fitting the observed disk rotational velocities and radii with a Mo, Mao, White (1998) (MMW) model requires unusually large disk spin parameters lambda_d>0.1 and disk-to-dark halo mass fraction m_d=0.2, close to the cosmic baryon fraction. The galaxies segregate into dispersion-dominated systems with 1<vmax/sigma<3, maximum rotational velocities vmax<200 km/s and disk half-light radii rd=1-3 kpc and rotation-dominated systems with vmax>200 km/s, vmax/sigma>3 and rd=4-8 kpc. For the dispersion-dominated sample, radial pressure gradients partly compensate the gravitational force, reducing the rotational velocities. Including this pressure effect in the MMW model, dispersion-dominated galaxies can be fitted well with spin parameters lf lambda_d=0.03-0.05 for high disk mass fractions of m_d=0.2 and with lambda_d=0.01-0.03 for m_d=0.05. These values are in good agreement with cosmological expectations. For the rotation-dominated sample however pressure effects are small and better agreement with theoretically expected disk spin parameters can only be achieved if the dark halo mass contribution in the visible disk regime (2-3*rd) is smaller than predicted by the MMW model. We argue that these galaxies can still be embedded in standard cold dark matter halos if the halos did not contract adiabatically in response to disk formation. It is shown that the observed high turbulent gas motions of the galaxies are consistent with a Toomre instability parameter Q=1 which is equal to the critical value, expected for gravitational disk instability to be the major driver of turbulence. The dominant energy source of turbulence is then the potential energy of the gas in the disk. △ Less

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

Comments: 23 pages, 4 figures, ApJ, in press

Lyman alpha Radiative Transfer with Dust: Escape Fractions from Simulated High-Redshift Galaxies

Authors: Peter Laursen, Jesper Sommer-Larsen, Anja C. Andersen

Abstract: The Lyman alpha emission line is an essential diagnostic tool for probing galaxy formation and evolution. Not only is it commonly the strongest observable line from high-redshift galaxies but from its shape detailed information about its host galaxy can be revealed. However, due to the scattering nature of Lya photons increasing their path length in a non-trivial way, if dust is present in the g… ▽ More The Lyman alpha emission line is an essential diagnostic tool for probing galaxy formation and evolution. Not only is it commonly the strongest observable line from high-redshift galaxies but from its shape detailed information about its host galaxy can be revealed. However, due to the scattering nature of Lya photons increasing their path length in a non-trivial way, if dust is present in the galaxy the line may be severely suppressed and its shape altered. In order to interpret observations correctly, it is thus of crucial significance to know how much of the emitted light actually escapes the galaxy. In the present work, using a combination of high-resolution cosmological hydro-simulations and an adaptively refinable Monte Carlo Lya radiative transfer code including an advanced model of dust, the escape fractions f_esc of Lya radiation from high-redshift (z = 3.6) galaxies are calculated. In addition to the average escape fraction, the variation of f_esc in different directions and from different parts of the galaxies is investigated, as well as the effect on the emergent spectrum. Escape fractions from a sample of simulated galaxies of representative physical properties are found to decrease for increasing galaxy virial mass M_vir, from f_esc approaching unity for M_vir ~ 10^9 M_sun to f_esc less than 10% for M_vir ~ 10^12 M_sun. In spite of the dust being almost grey, it is found that the emergent spectrum is affected non-uniformly, with the escape fraction of photons close to the line center being much higher than of those in the wings, thus effectively narrowing the Lya line. △ Less

Submitted 16 July, 2009; originally announced July 2009.

Comments: Submitted to ApJ; comments are very welcome

Moderate Steepening of Galaxy Cluster Dark Matter Profiles by Baryonic Pinching

Authors: Jesper Sommer-Larsen, Marceau Limousin

Abstract: To assess the effect of baryonic ``pinching'' of galaxy cluster dark matter (DM) haloes, cosmological (LCDM) TreeSPH simulations of the formation and evolution of two galaxy clusters have been performed, with and without baryons included. The simulations with baryons invoke star formation, chemical evolution with non-instantaneous recycling, metallicity dependent radiative cooling, strong star… ▽ More To assess the effect of baryonic ``pinching'' of galaxy cluster dark matter (DM) haloes, cosmological (LCDM) TreeSPH simulations of the formation and evolution of two galaxy clusters have been performed, with and without baryons included. The simulations with baryons invoke star formation, chemical evolution with non-instantaneous recycling, metallicity dependent radiative cooling, strong star-burst, driven galactic super-winds and the effects of a meta-galactic UV field, including simplified radiative transfer. The two clusters have T_X~3 and 6 keV, respectively, and, at z~0, both host a prominent, central cD galaxy. Comparing the simulations without and with baryons, it is found for the latter that the inner DM density profiles, r<50-100 kpc, steepen considerably: Delta(alpha)~0.5-0.6, where -alpha is the logarithmic DM density gradient. This is mainly due to the central stellar cDs becoming very massive, as a consequence of the onset of late time cooling flows and related star formation. Once these spurious cooling flows have been corrected for, and the cluster gravitational potentials dynamically adjusted, much smaller pinching effects are found: Delta(alpha)~0.1. Including the effects of baryonic pinching, central slopes of alpha~1.0 and 1.2 are found for the DM in the two clusters, interestingly close to recent observational findings. For the simulations with baryons, the inner density profile of DM+ICM gas combined is found to be only very marginally steeper than that of the DM, Delta(alpha)<0.05. However, the total matter inner density profiles are found to be Delta(alpha)~0.5 steeper than the inner profiles in the dark matter only simulations. △ Less

Submitted 2 June, 2009; originally announced June 2009.

Comments: 11 Pages, 8 Figures, submitted to MNRAS, printing in colour recommended

Ionizing radiation from z=4-10 galaxies

Authors: Alexei O. Razoumov, Jesper Sommer-Larsen

Abstract: We compute the escape of ionizing radiation from galaxies in the redshift interval z=4-10, i.e., during and after the epoch of reionization, using a high-resolution set of galaxies, formed in fully cosmological simulations. The simulations invoke early, energetic feedback, and the galaxies evolve into a realistic population at z=0. Our galaxies cover nearly four orders of magnitude in masses (10… ▽ More We compute the escape of ionizing radiation from galaxies in the redshift interval z=4-10, i.e., during and after the epoch of reionization, using a high-resolution set of galaxies, formed in fully cosmological simulations. The simulations invoke early, energetic feedback, and the galaxies evolve into a realistic population at z=0. Our galaxies cover nearly four orders of magnitude in masses (10^{7.8}-10^{11.5}\msun) and more than five orders in star formation rates (10^{-3.5}-10^{1.7}\msun\yr^{-1}), and we include an approximate treatment of dust absorption. We show that the source-averaged Lyman-limit escape fraction at z=10.4 is close to 80% declining monotonically with time as more massive objects build up at lower redshifts. Although the amount of dust absorption is uncertain to 1-1.5 dex, it is tightly correlated with metallicity; we find that dust is unlikely to significantly impact the observed UV output. These results support reionization by stellar radiation from low-luminosity dwarf galaxies and are also compatible with Lyman continuum observations and theoretical predictions at z\sim3-4. △ Less

Submitted 21 January, 2010; v1 submitted 12 March, 2009; originally announced March 2009.

Comments: ApJ accepted, 9 figures, printing in colour recommended, comments welcome; replaced version: 2 new figures, updated content

Journal ref: Astrophys.J.710:1239-1246,2010

Hot gas halos around disk galaxies: Confronting cosmological simulations with observations

Authors: Jesper Rasmussen, Jesper Sommer-Larsen, Kristian Pedersen, Sune Toft, Andrew Benson, Richard G. Bower, Lisbeth F. Grove

Abstract: Models of disk galaxy formation commonly predict the existence of an extended reservoir of accreted hot gas surrounding massive spirals at low redshift. As a test of these models, we use X-ray and H-alpha data of the two massive, quiescent edge-on spirals NGC 5746 and NGC 5170 to investigate the amount and origin of any hot gas in their halos. Contrary to our earlier claim, the Chandra analysis… ▽ More Models of disk galaxy formation commonly predict the existence of an extended reservoir of accreted hot gas surrounding massive spirals at low redshift. As a test of these models, we use X-ray and H-alpha data of the two massive, quiescent edge-on spirals NGC 5746 and NGC 5170 to investigate the amount and origin of any hot gas in their halos. Contrary to our earlier claim, the Chandra analysis of NGC 5746, employing more recent calibration data, does not reveal any significant evidence for diffuse X-ray emission outside the optical disk, with a 3-sigma upper limit to the halo X-ray luminosity of 4e39 erg/s. An identical study of the less massive NGC 5170 also fails to detect any extraplanar X-ray emission. By extracting hot halo properties of disk galaxies formed in cosmological hydrodynamical simulations, we compare these results to expectations for cosmological accretion of hot gas by spirals. For Milky Way-sized galaxies, these high-resolution simulations predict hot halo X-ray luminosities which are lower by a factor of ~2 compared to our earlier results reported by Toft et al. (2002). We find the new simulation predictions to be consistent with our observational constraints for both NGC 5746 and NGC 5170, while also confirming that the hot gas detected so far around more actively star-forming spirals is in general probably associated with stellar activity in the disk. Observational results on quiescent disk galaxies at the high-mass end are nevertheless providing powerful constraints on theoretical predictions, and hence on the assumed input physics in numerical studies of disk galaxy formation and evolution. △ Less

Submitted 3 March, 2009; originally announced March 2009.

Comments: 16 pages, 8 figures, ApJ accepted. Supersedes arXiv:astro-ph/0610893

Journal ref: Astrophys.J.697:79-93,2009

The SINS survey: modeling the dynamics of z~2 galaxies and the high-z Tully-Fisher relation

Authors: G. Cresci, E. K. S. Hicks, R. Genzel, N. M. Foerster Schreiber, R. Davies, N. Bouche', P. Buschkamp, S. Genel, K. Shapiro, L. Tacconi, J. Sommer-Larsen, A. Burkert, F. Eisenhauer, O. Gerhard, D. Lutz, T. Naab, A. Sternberg, A. Cimatti, E. Daddi, D. K. Erb, J. Kurk, S. L. Lilly, A. Renzini, A. Shapley, C. C. Steidel , et al. (1 additional authors not shown)

Abstract: We present the modeling of SINFONI integral field dynamics of 18 star forming galaxies at z ~ 2 from Halpha line emission. The galaxies are selected from the larger sample of the SINS survey, based on the prominence of ordered rotational motions with respect to more complex merger induced dynamics. The quality of the data allows us to carefully select systems with kinematics dominated by rotatio… ▽ More We present the modeling of SINFONI integral field dynamics of 18 star forming galaxies at z ~ 2 from Halpha line emission. The galaxies are selected from the larger sample of the SINS survey, based on the prominence of ordered rotational motions with respect to more complex merger induced dynamics. The quality of the data allows us to carefully select systems with kinematics dominated by rotation, and to model the gas dynamics across the whole galaxy using suitable exponential disk models. We obtain a good correlation between the dynamical mass and the stellar mass, finding that large gas fractions Mgas~M*) are required to explain the difference between the two quantities. We use the derived stellar mass and maximum rotational velocity Vmax from the modeling to construct for the first time the stellar mass Tully-Fisher relation at z ~ 2.2. The relation obtained shows a slope similar to what is observed at lower redshift, but we detect an evolution of the zero point. We find that at z ~ 2.2 there is an offset in log(M*) for a given rotational velocity of 0.41+-0.11 with respect to the local Universe. This result is consistent with the predictions of the latest N-body/hydrodynamical simulations of disk formation and evolution, which invoke gas accretion onto the forming disk in filaments and cooling flows. This scenario is in agreement with other dynamical evidence from SINS, where gas accretion from the halo is required to reproduce the observed properties of a large fraction of the z ~ 2 galaxies. △ Less

Submitted 26 February, 2009; originally announced February 2009.

Comments: 10 pages, 6 figures, submitted to ApJ

Journal ref: Astrophys.J.697:115-132,2009

Lyman alpha Radiative Transfer in Cosmological Simulations using Adaptive Mesh Refinement

Authors: Peter Laursen, Alexei O. Razoumov, Jesper Sommer-Larsen, ;

Abstract: A numerical code for solving various Lyman alpha (Lya) radiative transfer (RT) problems is presented. The code is suitable for an arbitrary, three-dimensional distribution of Lya emissivity, gas temperature, density, and velocity field. Capable of handling Lya RT in an adaptively refined grid-based structure, it enables detailed investigation of the effects of clumpiness of the interstellar (or… ▽ More A numerical code for solving various Lyman alpha (Lya) radiative transfer (RT) problems is presented. The code is suitable for an arbitrary, three-dimensional distribution of Lya emissivity, gas temperature, density, and velocity field. Capable of handling Lya RT in an adaptively refined grid-based structure, it enables detailed investigation of the effects of clumpiness of the interstellar (or intergalactic) medium. The code is tested against various geometrically and physically idealized configurations for which analytical solutions exist, and subsequently applied to three "Lyman-break galaxies", extracted from high-resolution cosmological simulations at redshift z = 3.6. Proper treatment of the Lya scattering reveals a diversity of surface brightness (SB) and line profiles. Specifically, for a given galaxy the maximum observed SB can vary by an order of magnitude, and the total flux by a factor of 3 - 6, depending on the viewing angle. This may provide an explanation for differences in observed properties of high-redshift galaxies, and in particular a possible physical link between Lyman-break galaxies and regular Lya emitters. △ Less

Submitted 16 July, 2009; v1 submitted 21 May, 2008; originally announced May 2008.

Comments: Minor corrections to match version published in ApJ, plus one typo concerning erroneously swapped "<" and ">" immediately after Eq. 17

Journal ref: Astrophys.J.696:853-869,2009; Erratum-ibid.702:824,2009

The evolution of the galaxy Red Sequence in simulated clusters and groups

Authors: Alessio D. Romeo, N. R. Napolitano, G. Covone, J. Sommer-Larsen, V. Antonuccio-Delogu, M. Capaccioli

Abstract: N-body + hydrodynamical simulations of the formation and evolution of galaxy groups and clusters in a LambdaCDM cosmology are used in order to follow the building-up of the colour-magnitude relation in two clusters and in 12 groups. We have found that galaxies, starting from the more massive, move to the Red Sequence (RS) as they get aged over times and eventually set upon a ``dead sequence'' (D… ▽ More N-body + hydrodynamical simulations of the formation and evolution of galaxy groups and clusters in a LambdaCDM cosmology are used in order to follow the building-up of the colour-magnitude relation in two clusters and in 12 groups. We have found that galaxies, starting from the more massive, move to the Red Sequence (RS) as they get aged over times and eventually set upon a ``dead sequence'' (DS) once they have stopped their bulk star formation activity. Fainter galaxies keep having significant star formation out to very recent epochs and lie broader around the RS. Environment plays a role as galaxies in groups and cluster outskirts hold star formation activity longer than the central cluster regions. However galaxies experiencing infall from the outskirts to the central parts keep star formation on until they settle on to the DS of the core galaxies. Merging contributes to mass assembly until z~1, after which major events only involve the brightest cluster galaxies. The emerging scenario is that the evolution of the colour-magnitude properties of galaxies within the hierarchical framework is mainly driven by star formation activity during dark matter halos assembly. Galaxies progressively quenching their star formation settle to a very sharp ``red and dead'' sequence, which turns out to be universal, its slope and scatter being almost independent of the redshift (since at least z~1.5) and environment. Differently from the DS, the operatively defined RS evolves more evidently with z, the epoch when it changes its slope being closely corresponding to that at which the passive galaxies population takes over the star forming one: this goes from z~1 in clusters down to 0.4 in normal groups. △ Less

Submitted 14 July, 2008; v1 submitted 9 April, 2008; originally announced April 2008.

Comments: MNRAS in press

Journal ref: 2008 MNRAS 389, 13

Strong Lensing in Abell 1703: Constraints on the Slope of the Inner Dark Matter Distribution

Authors: M. Limousin, J. Richard, J. -P. Kneib, H. Brink, R. Pello, E. Jullo, H. Tu, J. Sommer-Larsen, E. Egami, M. J. Michalowski, R. Cabanac, D. P. Stark

Abstract: In this article, we apply strong lensing techniques in Abell 1703, a massive X-ray luminous galaxy cluster at z=0.28. Our analysis is based on imaging data both from space and ground in 8 bands, complemented with a spectroscopic survey. Abell 1703 looks rather circular from the general shape of its multiply imaged systems and present a dominant giant elliptical cD galaxy in its centre. This clus… ▽ More In this article, we apply strong lensing techniques in Abell 1703, a massive X-ray luminous galaxy cluster at z=0.28. Our analysis is based on imaging data both from space and ground in 8 bands, complemented with a spectroscopic survey. Abell 1703 looks rather circular from the general shape of its multiply imaged systems and present a dominant giant elliptical cD galaxy in its centre. This cluster exhibits a remarkable bright 'central ring' formed by 4 bright images at z_{spec}=0.888 located very close to the cD galaxy, providing observational constraints that are potentially very interesting to probe the central mass distribution. The stellar contribution from the cD galaxy (~1.25 10^{12} M_{sun} within 7") is accounted for in our parametric mass modelling, and the underlying smooth dark matter component distribution is described using a generalized NFW profile parametrized with a central logarithmic slope α. We find that within the range where observational constraints are present (from ~5" to ~50"), the slope of the dark matter distribution in Abell1703 is equal to 1.09^{+0.05}_{-0.11} (3σconfidence level). The concentration parameter is equal to c_{200} ~ 3.5, and the scale radius is constrained to be larger than the region where observational constraints are available. Within this radius, the 2D mass is equal to M(50")=2.4 10^{14} M_{\sun}. We cannot draw any conclusions on cosmological models at this point since we lack results from realistic numerical simulations containing baryons to make a proper comparison. We advocate the need for a sample of observed and simulated unimodal relaxed galaxy clusters in order to make reliable comparisons, and potentially provide a test of cosmological models. △ Less

Submitted 17 July, 2008; v1 submitted 28 February, 2008; originally announced February 2008.

Comments: Accepted for publication, conclusion unchanged but strengthened by a serie of tests suggested by a constructive referee report. Higher resolution available at http://www.dark-cosmology.dk/~marceau/1703.php.html

Reconciling the Metallicity Distributions of Gamma-ray Burst, Damped Lyman-alpha, and Lyman-break Galaxies at z=3

Authors: Johan P. U. Fynbo, J. Xavier Prochaska, Jesper Sommer-Larsen, Miroslava Dessauges-Zavadsky, Palle Moller

Abstract: We test the hypothesis that the host galaxies of long-duration gamma-ray bursts (GRBs) as well as quasar-selected damped Lyman-alpha (DLA) systems are drawn from the population of UV-selected star-forming, high-z galaxies (generally referred to as Lyman-break galaxies). Specifically, we compare the metallicity distributions of the GRB and DLA populations to simple models where these galaxies are… ▽ More We test the hypothesis that the host galaxies of long-duration gamma-ray bursts (GRBs) as well as quasar-selected damped Lyman-alpha (DLA) systems are drawn from the population of UV-selected star-forming, high-z galaxies (generally referred to as Lyman-break galaxies). Specifically, we compare the metallicity distributions of the GRB and DLA populations to simple models where these galaxies are drawn randomly from the distribution of star-forming galaxies according to their star-formation rate and HI cross-section respectively. We find that it is possible to match both observational distributions assuming very simple and constrained relations between luminosity, metallicity and HI sizes. The simple model can be tested by observing the luminosity distribution of GRB host galaxies and by measuring the luminosity and impact parameters of DLA selected galaxies as a function of metallicity. Our results support the expectation that GRB and DLA samples, in contrast to magnitude limited surveys, provide an almost complete census of z=3 star-forming galaxies that are not heavily dust-obscured. △ Less

Submitted 13 April, 2008; v1 submitted 22 January, 2008; originally announced January 2008.

Comments: 31 pages, 6 figures. Accepted for publication in ApJ

Journal ref: Astrophys. J. 683 (2008) 321-328

Sunyaev-Zeldovich profiles for clusters and groups of galaxies

Authors: Troels Haugboelle, Jesper Sommer-Larsen, Kristian Pedersen

Abstract: The Sunyaev-Zeldovich (SZ) effect gives a measure of the thermal energy and electron pressure in groups and clusters of galaxies. In the near future SZ surveys will map hundreds of systems, shedding light on the pressure distribution in the systems. The thermal energy is related to the total mass of a system of galaxies, but it is only a projection that is observed through the SZ effect. A model… ▽ More The Sunyaev-Zeldovich (SZ) effect gives a measure of the thermal energy and electron pressure in groups and clusters of galaxies. In the near future SZ surveys will map hundreds of systems, shedding light on the pressure distribution in the systems. The thermal energy is related to the total mass of a system of galaxies, but it is only a projection that is observed through the SZ effect. A model for the 3D distribution of pressure is needed to link the SZ signal to the total mass of the system. In this work we construct an empirical model for the 2D and 3D SZ profile, and compare it to a set of realistic high resolution SPH simulations of galaxy clusters and groups, and to a stacked SZ profile for massive clusters derived from WMAP data. Furthermore, we combine observed temperature profiles with dark matter potentials to yield an additional constraint, under the assumption of hydrostatic equilibrium. We find a very tight correlation between the characteristic scale in the model, the integrated SZ signal, and the total mass in the systems with a scatter of only 4%. The model only contains two free parameters, making it readily applicable even to low resolution SZ observations of galaxy clusters. A fitting routine for the model that can be applied to observed or simulated data can be found at http://www.phys.au.dk/~haugboel/software.shtml △ Less

Submitted 17 December, 2007; originally announced December 2007.

Comments: 7 pages, 9 figures

Can gravitational infall energy lead to the observed velocity dispersion in DLAs?

Authors: Alexei O. Razoumov, Michael L. Norman, Jason X. Prochaska, Jesper Sommer-Larsen, Arthur M. Wolfe, Yi-Jung Yang

Abstract: The median observed velocity width v_90 of low-ionization species in damped Ly-alpha systems is close to 90 km/s, with approximately 10% of all systems showing v_90 > 210 km/s at z=3. We show that a relative shortage of such high-velocity neutral gas absorbers in state-of-the-art galaxy formation models is a fundamental problem, present both in grid-based and particle-based numerical simulations… ▽ More The median observed velocity width v_90 of low-ionization species in damped Ly-alpha systems is close to 90 km/s, with approximately 10% of all systems showing v_90 > 210 km/s at z=3. We show that a relative shortage of such high-velocity neutral gas absorbers in state-of-the-art galaxy formation models is a fundamental problem, present both in grid-based and particle-based numerical simulations. Using a series of numerical simulations of varying resolution and box size to cover a wide range of halo masses, we demonstrate that energy from gravitational infall alone is insufficient to produce the velocity dispersion observed in damped Ly-alpha systems, nor does this dispersion arise from an implementation of star formation and feedback in our highest resolution (~ 45 pc) models, if we do not put any galactic winds into our models by hand. We argue that these numerical experiments highlight the need to separate dynamics of different components of the multiphase interstellar medium at z=3. △ Less

Submitted 5 April, 2008; v1 submitted 22 October, 2007; originally announced October 2007.

Comments: 12 Pages, 9 Figures, accepted to ApJ, printing in colour recommended

Probing the truncation of galaxy dark matter halos in high density environments from hydrodynamical N-body simulations

Authors: Marceau Limousin, Jesper Sommer-Larsen, Priyamvada Natarajan, Bo Milvang-Jensen

Abstract: We analyze high resolution, N-body hydrodynamical simulations of fiducial galaxy clusters to probe tidal stripping of the dark matter subhalos. These simulations include a prescription for star formation allowing us to track the fate of the stellar component as well. We investigate the effect of tidal stripping on cluster galaxies hosted in these dark matter subhalos as a function of cluster-cen… ▽ More We analyze high resolution, N-body hydrodynamical simulations of fiducial galaxy clusters to probe tidal stripping of the dark matter subhalos. These simulations include a prescription for star formation allowing us to track the fate of the stellar component as well. We investigate the effect of tidal stripping on cluster galaxies hosted in these dark matter subhalos as a function of cluster-centric radius. To quantify the extent of the dark matter halos of cluster galaxies, we introduce the half mass radius r_half as a diagnostic, and study its evolution with projected cluster-centric distance R as a function of redshift. We find a well defined trend for (r_half,R): the closer the galaxies are to the center of the cluster, the smaller the half mass radius. Interestingly, this trend is inferred in all redshift frames examined in this work ranging from z=0 to z=0.7. At z=0, galaxy halos in the central regions of clusters are found to be highly truncated, with the most compact half mass radius of 10 kpc. We also find that r_half depends on luminosity and we present scaling relations of r_half with galaxy luminosity. The corresponding total mass of the cluster galaxies is also found to increase with projected cluster-centric distance and luminosity, but with more scatter than the (r_half,R) trend. Comparing the distribution of stellar mass to total mass for cluster galaxies, we find that the dark matter component is preferentially stripped, whereas the stellar component remains protected by the halo and is much less affected by tidal forces. We compare these results with galaxy-galaxy lensing probes of r_half and find qualitative agreement. (Abridged) △ Less

Submitted 3 March, 2009; v1 submitted 21 June, 2007; originally announced June 2007.

Comments: Accepted for publication in ApJ

Journal ref: Astrophys.J.696:1771-1779,2009

Ongoing Galactic Accretion: Simulations and Observations of Condensed Gas in Hot Halos

Authors: J. E. G. Peek, M. E. Putman, Jesper Sommer-Larsen

Abstract: Ongoing accretion onto galactic disks has been recently theorized to progress via the unstable cooling of the baryonic halo into condensed clouds. These clouds have been identified as analogous to the High-Velocity Clouds (HVCs) observed in HI in our Galaxy. Here we compare the distribution of HVCs observed around our own Galaxy and extra-planar gas around the Andromeda galaxy to these possible… ▽ More Ongoing accretion onto galactic disks has been recently theorized to progress via the unstable cooling of the baryonic halo into condensed clouds. These clouds have been identified as analogous to the High-Velocity Clouds (HVCs) observed in HI in our Galaxy. Here we compare the distribution of HVCs observed around our own Galaxy and extra-planar gas around the Andromeda galaxy to these possible HVC analogs in a simulation of galaxy formation that naturally generates these condensed clouds. We find a very good correspondence between these observations and the simulation, in terms of number, angular size, velocity distribution, overall flux and flux distribution of the clouds. We show that condensed cloud accretion only accounts for ~ 0.2 M_solar / year of the current overall Galactic accretion in the simulations. We also find that the simulated halo clouds accelerate and become more massive as they fall toward the disk. The parameter space of the simulated clouds is consistent with all of the observed HVC complexes that have distance constraints, except the Magellanic Stream which is known to have a different origin. We also find that nearly half of these simulated halo clouds would be indistinguishable from lower-velocity gas and that this effect is strongest further from the disk of the galaxy, thus indicating a possible missing population of HVCs. These results indicate that the majority of HVCs are consistent with being infalling, condensed clouds that are a remnant of Galaxy formation. △ Less

Submitted 11 September, 2007; v1 submitted 2 May, 2007; originally announced May 2007.

Comments: 10 pages, 6 figures, ApJ Accepted. Some changes to technique

Coincident, 100 kpc-scale damped Lyman alpha absorption towards a binary QSO: how large are galaxies at z ~ 3?

Authors: Sara L. Ellison, Joseph F. Hennawi, Crystal L. Martin, Jesper Sommer-Larsen

Abstract: We report coincident damped Lyman alpha (DLA) and sub-DLA absorption at z = 2.66 and z = 2.94 towards the z ~ 3 13.8 arcsecond separation binary quasar SDSS 1116+4118 AB. At the redshifts of the absorbers, this angular separation corresponds to a proper transverse separation of ~ 110 kpc. A third absorber, a sub-DLA at z = 2.47, is detected towards SDSS 1116+4118 B, but no corresponding high col… ▽ More We report coincident damped Lyman alpha (DLA) and sub-DLA absorption at z = 2.66 and z = 2.94 towards the z ~ 3 13.8 arcsecond separation binary quasar SDSS 1116+4118 AB. At the redshifts of the absorbers, this angular separation corresponds to a proper transverse separation of ~ 110 kpc. A third absorber, a sub-DLA at z = 2.47, is detected towards SDSS 1116+4118 B, but no corresponding high column density absorber is present towards SDSS 1116+4118 A. We use high resolution galaxy simulations and a clustering analysis to interpret the coincident absorption and its implications for galaxy structure at z ~ 3. We conclude that the common absorption in the two lines of sight is unlikely to arise from a single galaxy, or a galaxy plus satellite system, and is more feasibly explained by a group of two or more galaxies with separations ~ 100 kpc. The impact of these findings on single line of sight observations is also discussed; we show that abundances of DLAs may be affected by up to a few tenths of a dex by line of sight DLA blending. From a Keck ESI spectrum of the two quasars, we measure metal column densities for all five absorbers and determine abundances for the three absorbers with log N(HI) > 20. For the two highest N(HI) absorbers, we determine high levels of metal enrichment, corresponding to 1/3 and 1/5 solar. These metallicities are amongst the highest measured for DLAs at any redshift and are consistent with values measured in Lyman break galaxies at 2 < z < 3. For the DLA at z = 2.94 we also infer an approximately solar ratio of alpha-to-Fe peak elements from [S/Zn] = +0.05, and measure an upper limit for the molecular fraction in this particular line of sight of log f(H_2)< -5.5. △ Less

Submitted 13 April, 2007; originally announced April 2007.

Comments: Accepted for publication in MNRAS. Full resolution simulation images available in pdf copy of the manuscript at http://www.astro.uvic.ca/~sara/1116.pdf

Journal ref: Mon.Not.Roy.Astron.Soc.378:801-818,2007

Where are the Cosmic Metals at z~3?

Authors: Jesper Sommer-Larsen, Johan P. U. Fynbo

Abstract: The temperature distribution of the cosmic gas-phase oxygen at z~3 is determined by combining high resolution cosmological simulations with the observed, (rest-frame) V-band galaxy LF. The simulations invoke three different stellar IMFs, a Kroupa (K98), a Salpeter (S) and an Arimoto-Yoshii (AY), spanning a range of a factor of five in chemical yield and specific SNII energy feedback. Oxygen is… ▽ More The temperature distribution of the cosmic gas-phase oxygen at z~3 is determined by combining high resolution cosmological simulations with the observed, (rest-frame) V-band galaxy LF. The simulations invoke three different stellar IMFs, a Kroupa (K98), a Salpeter (S) and an Arimoto-Yoshii (AY), spanning a range of a factor of five in chemical yield and specific SNII energy feedback. Oxygen is found in all T phases (cold: logT<~4, warm: logT~4.5, hot: logT>~5), in particular for the (top-heavy) AY IMF. For the K98 and S IMFs, the most important phase is the cold one. Moreover, the cold phase alone contains 1-3 times the mass of oxygen in galactic stars for the three IMFs. The implications of this in relation to DLA studies are discussed. In relation to ``missing metals'' it is found that a significant fraction of the oxygen is located in a warm/hot phase that may be very difficult to detect. Moreover, it is found that less than about 20-25% of the cosmic oxygen is associated with galaxies brighter than M_V~-22, the faintest LBGs probed by current metallicity determinations. From the LBG based, lambda~1500 AA luminosity density history at z>3, the mean oxygen density at z=3 is inferred. Comparing this to the models, it is found that the (solar neighbourhood type) K98 IMF is strongly excluded (too small yield), the S is marginally excluded, and the AY matches the constraint well. The optimal IMF has a yield between the S and AY. A similar conclusion is reached on the basis of the C IV abundance of the moderate density IGM. △ Less

Submitted 23 October, 2007; v1 submitted 28 March, 2007; originally announced March 2007.

Comments: 22 Pages, 16 Figures, MNRAS, Accepted version, Printing in colour recommended

Journal ref: Mon.Not.Roy.Astron.Soc.385:3-22,2008

Modeling Lyman continuum emission from young galaxies

Authors: Alexei O. Razoumov, Jesper Sommer-Larsen

Abstract: Based on cosmological simulations, we model Lyman continuum emission from a sample of 11 high-redshift star forming galaxies spanning a mass range of a factor 20. Each of the 11 galaxies has been simulated both with a Salpeter and a Kroupa initial mass function (IMF). We find that the Lyman continuum (LyC) luminosity of an average star forming galaxy in our sample declines from z=3.6 to 2.4 due… ▽ More Based on cosmological simulations, we model Lyman continuum emission from a sample of 11 high-redshift star forming galaxies spanning a mass range of a factor 20. Each of the 11 galaxies has been simulated both with a Salpeter and a Kroupa initial mass function (IMF). We find that the Lyman continuum (LyC) luminosity of an average star forming galaxy in our sample declines from z=3.6 to 2.4 due to the steady gas infall and higher gas clumping at lower redshifts, increasingly hampering the escape of ionizing radiation. The galaxy-to-galaxy variation of apparent LyC emission at a fixed redshift is caused in approximately equal parts by the intrinsic variations in the LyC emission and by orientation effects. The combined scatter of an order of magnitude can explain the variance in the far-UV spectra of high-redshift galaxies detected by Shapley et al. (2006). Our results imply that the cosmic galactic ionizing UV luminosity would be monotonically decreasing from z=3.6 to 2.4, curiously anti-correlated with the star formation rate in the smaller galaxies, which on average rises during this redshift interval. △ Less

Submitted 21 June, 2007; v1 submitted 12 February, 2007; originally announced February 2007.

Comments: 8 pages, 12 figures, ApJ, in press

Evolution of the Red Sequence in simulated galaxy groups

Authors: A. D. Romeo, N. R. Napolitano, G. Covone, J. Sommer-Larsen

Abstract: N--body + hydrodynamical simulations of formation and evolution of galaxy groups in a LambdaCDM cosmology have been performed. The properties of the galaxy populations in 12 groups are here discussed, with focus on the colour-magnitude relation in both normal and fossil groups. N--body + hydrodynamical simulations of formation and evolution of galaxy groups in a LambdaCDM cosmology have been performed. The properties of the galaxy populations in 12 groups are here discussed, with focus on the colour-magnitude relation in both normal and fossil groups. △ Less

Submitted 15 December, 2006; originally announced December 2006.

Comments: Proceedings of the Conference "DEEP06 - At the edge of the Universe", Sintra, 9-13 October 2006; 5 pages, 3 colour figures

Investigating hot gas in the halos of two massive spirals: Observations and cosmological simulations

Authors: Jesper Rasmussen, Jesper Sommer-Larsen, Kristian Pedersen, Sune Toft, Andrew Benson, Richard G. Bower, Lisbeth F. Olsen

Abstract: Models of disk galaxy formation commonly predict the existence of an extended reservoir of hot gas surrounding massive spirals at low redshift. As a test of these models, we have obtained X-ray and optical data of the two massive edge-on spirals NGC 5746 and NGC 5170, in order to investigate the amount and origin of hot gas in their disks and halos. Chandra observations of NGC 5746 reveal eviden… ▽ More Models of disk galaxy formation commonly predict the existence of an extended reservoir of hot gas surrounding massive spirals at low redshift. As a test of these models, we have obtained X-ray and optical data of the two massive edge-on spirals NGC 5746 and NGC 5170, in order to investigate the amount and origin of hot gas in their disks and halos. Chandra observations of NGC 5746 reveal evidence for diffuse X-ray emission with a total luminosity of ~7 x 10^39 erg/s surrounding this galaxy out to at least ~20 kpc from the disk, whereas an identical study of the less massive NGC 5170 fails to detect any extraplanar X-ray emission. Unlike the case for other disk galaxies with detected X-ray halos, the halo emission around NGC 5746 is not accompanied by extraplanar H-alpha or radio emission, and there is no evidence for significant nuclear or starburst activity in the disk. In contrast to these other cases, the emission around NGC 5746 therefore appears to arise from the cooling of externally accreted material rather than from disk outflows. To verify this idea, we present results of cosmological simulations of galaxy formation and evolution, showing our observations to be in good agreement with expectations for cosmological accretion, while also confirming that the X-ray halos of other spirals do not fit well into an accretion scenario. We find that the estimated cooling rate of hot halo gas around NGC 5746 would provide sufficient material for star formation in the disk to proceed at its present rate. This lends support to the idea that a supply of hot ambient gas is potentially available as fuel for star formation in massive, nearby spirals, and suggests that accretion of hot gas could be important for maintaining the stellar disks of such galaxies. △ Less

Submitted 30 October, 2006; originally announced October 2006.

Comments: 20 pages, 8 figures, submitted to ApJ

Lyman alpha Resonant Scattering in Young Galaxies - Predictions from Cosmological Simulations

Authors: P. Laursen, J. Sommer-Larsen

Abstract: We present results obtained with a 3D, Ly alpha radiative transfer code, applied to a fully cosmological galaxy formation simulation. The developed Monte Carlo code is capable of treating an arbitrary distribution of source Ly alpha emission, neutral hydrogen density, temperature, and peculiar velocity of the interstellar medium. We investigate the influence of resonant scattering on the appeara… ▽ More We present results obtained with a 3D, Ly alpha radiative transfer code, applied to a fully cosmological galaxy formation simulation. The developed Monte Carlo code is capable of treating an arbitrary distribution of source Ly alpha emission, neutral hydrogen density, temperature, and peculiar velocity of the interstellar medium. We investigate the influence of resonant scattering on the appearance and properties of young galaxies by applying the code to a simulated "Lyman Break Galaxy" at redshift z = 3.6, and of star formation rate 22 M_sun/yr and total Ly alpha luminosity 2.0 X 10^43 erg/s. It is found that resonant scattering of Ly alpha radiation can explain that young galaxies frequently are observed to be more extended on the sky in Ly alpha than in the optical. Moreover, it is shown that, for the system investigated, due to the anisotropic escape of the photons, the appearent maximum surface brightness can differ by a factor of ~15, and the total derived luminosity by a factor of ~4, depending on the orientation of the system relative to the observer. △ Less

Submitted 18 July, 2009; v1 submitted 25 October, 2006; originally announced October 2006.

Comments: Letter updated to match version published in ApJ

Journal ref: Astrophys.J.Lett.657:L69-L72,2007

Broadband UBVRI Photometry of Horizontal-Branch and Metal-Poor Candidates from the HK and Hamburg/ESO Surveys. I

Authors: Timothy C. Beers, Chris Flynn, Silvia Rossi, Jesper Sommer-Larsen, Ronald Wilhelm, Brian Marsteller, YoungSun Lee, Nathan De Lee, Julie Krugler, Constantine P. Deliyannis, Franz-Josef Zickgraf, Johan Holmberg, Anna Onehag, Anders Eriksson, Donald M. Terndrup, Samir Salim, Johannes Andersen, Birgitta Nordstrom, Norbert Christlieb, Anna Frebel

Abstract: We report broadband UBV and/or BVRI CCD photometry for a total of 1857 stars in the thick-disk and halo populations of the Galaxy. The majority of our targets were selected as candidate field horizontal-branch or other A-type stars (FHB/A, N = 576), or candidate low-metallicity stars (N = 1221), from the HK and Hamburg/ESO objective-prism surveys. Similar data for a small number of additional st… ▽ More We report broadband UBV and/or BVRI CCD photometry for a total of 1857 stars in the thick-disk and halo populations of the Galaxy. The majority of our targets were selected as candidate field horizontal-branch or other A-type stars (FHB/A, N = 576), or candidate low-metallicity stars (N = 1221), from the HK and Hamburg/ESO objective-prism surveys. Similar data for a small number of additional stars from other samples are also reported. These data are being used for several purposes. In the case of the FHB/A candidates they are used to accurately separate the lower-gravity FHB stars from various higher-gravity A-type stars, a subsample that includes the so-called Blue Metal Poor stars, halo and thick-disk blue stragglers, main-sequence A-type dwarfs, and Am and Ap stars. These data are also being used to derive photometric distance estimates to high-velocity hydrogen clouds in the Galaxy and for improved measurements of the mass of the Galaxy. Photometric data for the metal-poor candidates are being used to refine estimates of stellar metallicity for objects with available medium-resolution spectroscopy, to obtain distance estimates for kinematic analyses, and to establish initial estimates of effective temperature for analysis of high-resolution spectroscopy of the stars for which this information now exists. △ Less

Submitted 1 October, 2006; originally announced October 2006.

Comments: 22 pages, including 3 figures, 5 tables, and two ascii files of full data, accepted for publication in the Astrophysical Journal (Supplements)

Escape of ionizing radiation from star forming regions in young galaxies

Authors: Alexei O. Razoumov, Jesper Sommer-Larsen

Abstract: Using results from high-resolution galaxy formation simulations in a standard Lambda-CDM cosmology and a fully conservative multi-resolution radiative transfer code around point sources, we compute the energy-dependent escape fraction of ionizing photons from a large number of star forming regions in two galaxies at five different redshifts from z=3.8 to 2.39. All escape fractions show a monoton… ▽ More Using results from high-resolution galaxy formation simulations in a standard Lambda-CDM cosmology and a fully conservative multi-resolution radiative transfer code around point sources, we compute the energy-dependent escape fraction of ionizing photons from a large number of star forming regions in two galaxies at five different redshifts from z=3.8 to 2.39. All escape fractions show a monotonic decline with time, from (at the Lyman-limit) ~6-10% at z=3.6 to ~1-2% at z=2.39, due to higher gas clumping at lower redshifts. It appears that increased feedback can lead to higher f_esc at z>3.4 via evacuation of gas from the vicinity of star forming regions and to lower f_esc at z<2.39 through accumulation of swept-up shells in denser environments. Our results agree well with the observational findings of \citet{inoue..06} on redshift evolution of f_esc in the redshift interval z=2-3.6. △ Less

Submitted 19 September, 2006; originally announced September 2006.

Comments: four pages, four figures, submitted to ApJL

Journal ref: Astrophys.J.651:L89-L92,2006

CO line emission from Lyman Break Galaxies. Cosmological simulations and predictions for ALMA

Authors: T. R. Greve, J. Sommer-Larsen

Abstract: [Abridged] The detection of the rotational lines of CO in proto-galaxies in the early Universe provides one of the most promising ways of probing the fundamental physical properties of a galaxy, such as its size, dynamical mass, gas density, and temperature. While such observations are currently limited to the most luminous galaxies, the advent of ALMA will change the situation dramatically, res… ▽ More [Abridged] The detection of the rotational lines of CO in proto-galaxies in the early Universe provides one of the most promising ways of probing the fundamental physical properties of a galaxy, such as its size, dynamical mass, gas density, and temperature. While such observations are currently limited to the most luminous galaxies, the advent of ALMA will change the situation dramatically, resulting in the detection of numerous normal galaxies at high redshifts. Maps and spectra of rotational CO line emission were calculated from a cosmological N-body/hydrodynamical TreeSPH simulation of a z ~ 3 "Lyman break galaxy'' of UV luminosity about one order of magnitude below L*. To simulate a typical observation of our system with ALMA, we imposed characteristic noise, angular, and spectral resolution constraints. The CO line properties predicted by our simulation are in good agreement with the two Lyman break systems detected in CO to date. We find that while supernovae explosions from the ongoing star formation carve out large cavities in the molecular ISM, they do not generate large enough gas outflows to make a substantial imprint on the CO line profile. This implies that for most proto-galaxies - except possibly the most extreme cases - stellar feedback effects do not affect CO as a dynamical mass tracer. Detecting CO in sub-L* galaxies at z ~3 will push ALMA to the limits of its cababilities, and whether a source is detected or not may depend critically on its inclination angle. Both these effects (sensitivity and inclination) will severely impair the ability of ALMA to infer the gas kinematics and dynamical masses using line observations. △ Less

Submitted 27 March, 2008; v1 submitted 30 August, 2006; originally announced August 2006.

Comments: 4 pages (letter format), 1 figure. 2008, A&A, 480, 335

Baryonic Pinching of Galactic Dark Matter Halos

Authors: Michael Gustafsson, Malcolm Fairbairn, Jesper Sommer-Larsen

Abstract: High resolution cosmological N-body simulations of four galaxy-scale dark matter halos are compared to corresponding N-body/hydrodynamical simulations containing dark matter, stars and gas. The simulations without baryons share features with others described in the literature in that the dark matter density slope continuously decreases towards the center, with a density ρ~r^{-1.3+/-0.2}, at abou… ▽ More High resolution cosmological N-body simulations of four galaxy-scale dark matter halos are compared to corresponding N-body/hydrodynamical simulations containing dark matter, stars and gas. The simulations without baryons share features with others described in the literature in that the dark matter density slope continuously decreases towards the center, with a density ρ~r^{-1.3+/-0.2}, at about 1% of the virial radius for our Milky Way sized galaxies. The central cusps in the simulations which also contain baryons steepen significantly, to ρ~r^{-1.9+/-0.2}, with an indication of the inner logarithmic slope converging. Models of adiabatic contraction of dark matter halos due to the central build-up of stellar/gaseous galaxies are examined. The simplest and most commonly used model, by Blumenthal et al., is shown to overestimate the central dark matter density considerably. A modified model proposed by Gnedin et al. is shown to be a considerable improvement, but not perfect. Moreover it is found that the contraction parameters not only depend on the orbital structure of the dark-matter-only halos but also on the stellar feedback prescription which is most relevant for the baryonic distribution. Implications for dark matter annihilation at the galactic center are discussed and it is found that although our simulations show a considerable reduced halo contraction as compared to the Blumenthal et al. model, the fluxes from dark matter annihilation is still expected to be enhanced by at least a factor of a hundred as compared to dark-matter-only halos. Finally, it is shown that while dark-matter-only halos are typically prolate, the dark matter halos containing baryons are mildly oblate with minor-to-major axis ratios of c/a=0.73+/-0.11, with their flattening aligned with the central baryonic disks. △ Less

Submitted 3 January, 2007; v1 submitted 29 August, 2006; originally announced August 2006.

Comments: 16 pages, 14 figures and 8 tables. Minor text improvements, typo corrections and an updated reference list to match the published version

Journal ref: Phys.Rev. D74 (2006) 123522

The Tully-Fisher relation and its evolution with redshift in cosmological simulations of disc galaxy formation

Authors: Laura Portinari, Jesper Sommer-Larsen

Abstract: We present predictions on the evolution of the Tully-Fisher (TF) relation with redshift, based on cosmological N-body/hydrodynamical simulations of disc galaxy formation and evolution. The simulations invoke star formation and stellar feedback, chemical evolution with non-instantaneous recycling, metallicity dependent radiative cooling and effects of a meta-galactic UV field, including simplifie… ▽ More We present predictions on the evolution of the Tully-Fisher (TF) relation with redshift, based on cosmological N-body/hydrodynamical simulations of disc galaxy formation and evolution. The simulations invoke star formation and stellar feedback, chemical evolution with non-instantaneous recycling, metallicity dependent radiative cooling and effects of a meta-galactic UV field, including simplified radiative transfer. At z=0, the simulated and empirical TF relations are offset by about 0.4 magnitudes (1 sigma) in the B and I bands. The origin of these offsets is somewhat unclear, but it may not necessarily be a problem of the simulations only. As to evolution, we find a brightening of the TF relation between z=0 and z=1 of about 0.85 mag in rest-frame B band, with a non-evolving slope. The brightening we predict is intermediate between the (still quite discrepant) observational estimates. This evolution is primarily a luminosity effect, while the stellar mass TF relation shows negligible evolution. The individual galaxies do gain stellar mass between z=1 and z=0, by a 50-100%; but they also correspondingly increase their characteristic circular speed. As a consequence, individually they mainly evolve ALONG the stellar mass TF relation, while the relation as such does not show any significant evolution. △ Less

Submitted 29 November, 2006; v1 submitted 21 June, 2006; originally announced June 2006.

Comments: 14 pages, 6 figures; minor revisions and added references; version accepted for publication on MNRAS. Fig. 1,2,3,4 available also in colour version from the source files

Journal ref: Mon.Not.Roy.Astron.Soc.375:913-924,2007

The density profiles of hot galactic halo gas

Authors: Steen H. Hansen, Jesper Sommer-Larsen

Abstract: Extended gas haloes around galaxies are a ubiquitous prediction of galaxy formation scenarios. However, the density profiles of this hot halo gas is virtually unknown, although various profiles have been suggested on theoretical grounds. In order to quantitatively address the gas profile, we compare galaxies from direct cosmological simulations with analytical solutions of the underlying gas equ… ▽ More Extended gas haloes around galaxies are a ubiquitous prediction of galaxy formation scenarios. However, the density profiles of this hot halo gas is virtually unknown, although various profiles have been suggested on theoretical grounds. In order to quantitatively address the gas profile, we compare galaxies from direct cosmological simulations with analytical solutions of the underlying gas equations. We find remarkable agreement between simulations and theoretical predictions. We present an expression for this gas profile with a non-trivial dependence on the total mass profile. This expression is useful when setting up equilibrium galaxy models for numerical experiments. △ Less

Submitted 13 June, 2006; originally announced June 2006.

Comments: 5 pages, 2 figures

Journal ref: Astrophys.J.653:L17-L20,2006

Where are the ``Missing'' Galactic Baryons?

Authors: Jesper Sommer-Larsen

Abstract: Based on 19 high-resolution N-body/gas-dynamical galaxy formation simulations in the LCDM cosmology it is shown, that for a galaxy like the Milky Way, in addition to the baryonic mass of the galaxy itself, about 70% extra baryonic mass should reside around the galaxy (inside of the virial radius), chiefly in the form of hot gas. Averaging over the entire field galaxy population, this ``external'… ▽ More Based on 19 high-resolution N-body/gas-dynamical galaxy formation simulations in the LCDM cosmology it is shown, that for a galaxy like the Milky Way, in addition to the baryonic mass of the galaxy itself, about 70% extra baryonic mass should reside around the galaxy (inside of the virial radius), chiefly in the form of hot gas. Averaging over the entire field galaxy population, this ``external'' component amounts to 64-85% of the baryonic mass of the population itself. These results are supported by the recent detection of very extended, soft X-ray emission from the halo of the quiescent, massive disk galaxy NGC 5746. Some of the hot gas may, by thermal instability, have condensed into mainly pressure supported, warm clouds, similar to the Galactic High Velocity Clouds (HVCs). Based on an ultra-high resolution cosmological test simulation of a Milky Way like galaxy (with a gas particle mass and gravity softening length of only 7600 h^-1 Msun and 83 h^-1 pc,respectively), it is argued, that the hot gas phase dominates over the warm gas phase, in the halo. Finally, an origin of HVCs as ``leftovers'' from filamentary, ``cold'' accretion events, mainly occurring early in the history of galaxies, is proposed. △ Less

Submitted 27 February, 2006; originally announced February 2006.

Comments: 4 Pages, 4 Figures, Submitted to ApJL, Printing in colour recommended

Journal ref: Astrophys.J. 644 (2006) L1-L4

Soft X-rays from disk galaxy halos, Ly-alpha emission from forming galaxies, and the z~1 Tully-Fisher relation

Authors: Jesper Sommer-Larsen

Abstract: Extended, soft X-ray emission from the halo of a very large disk galaxy has been detected out to r~20-40 kpc. The luminosity and surface brightness distribution is in very good agreement with predictions by recent, cosmological galaxy formation models. Predictions of Ly-alpha emission, associated with ``cold'' accretion of filamentary gas onto galaxies, are presented. In particular these predict… ▽ More Extended, soft X-ray emission from the halo of a very large disk galaxy has been detected out to r~20-40 kpc. The luminosity and surface brightness distribution is in very good agreement with predictions by recent, cosmological galaxy formation models. Predictions of Ly-alpha emission, associated with ``cold'' accretion of filamentary gas onto galaxies, are presented. In particular these predictions are compared to the properties of a recently detected Ly-alpha ``blob'', which shows no obvious continuum counter-part in any of 12 broad-bands, from X-rays (Chandra) to infrared (Spitzer). Finally, the predicted evolution of the Tully-Fisher relation, going from z=0 to 1, is discussed in relation to recent observations. △ Less

Submitted 19 December, 2005; originally announced December 2005.

Comments: 6 pages, To appear in the proceedings of the "Island Universes: Structure and Evolution of Disk Galaxies" conference held in Terschelling, Netherlands, July 2005, ed. R. de Jong (Springer: Dordrecht)

A Lyman-alpha blob in the GOODS South field: evidence for cold accretion onto a dark matter halo

Authors: Kim Nilsson, Johan P. U. Fynbo, Palle Moller, Jesper Sommer-Larsen, Cedric Ledoux

Abstract: We report on the discovery of a z = 3.16 Lyman-alpha emitting blob in the GOODS South field. The blob has a total Ly-alpha luminosity of ~ 10^(43) erg s^(-1) and a diameter larger than 60 kpc. The available multi-wavelength data in the GOODS field consists of 13 bands from X-rays (Chandra) to infrared (Spitzer). Unlike other discovered Ly-alpha blobs, this blob shows no obvious continuum counter… ▽ More We report on the discovery of a z = 3.16 Lyman-alpha emitting blob in the GOODS South field. The blob has a total Ly-alpha luminosity of ~ 10^(43) erg s^(-1) and a diameter larger than 60 kpc. The available multi-wavelength data in the GOODS field consists of 13 bands from X-rays (Chandra) to infrared (Spitzer). Unlike other discovered Ly-alpha blobs, this blob shows no obvious continuum counter-part in any of the broad-bands. In particular, no optical counter-parts are found in the deep HST/ACS imaging available. For previously published blobs, AGN (Active Galactic Nuclei) or 'superwind' models have been found to provide the best match with the data. We here argue that the most probable origin of the extended Ly-alpha emission from the blob in the GOODS South field is cold accretion onto a dark matter halo. △ Less

Submitted 24 April, 2006; v1 submitted 15 December, 2005; originally announced December 2005.

Comments: 4 pages, 2 tables, 2 figures, Accepted to A&A Letters, minor changes to text

Journal ref: Astron.Astrophys. 452 (2006) L23-L26

Discovery of a very extended X-ray halo around a quiescent spiral galaxy - the "missing link" of galaxy formation

Authors: Kristian Pedersen, Jesper Rasmussen, Jesper Sommer-Larsen, Sune Toft, Andrew J. Benson, Richard G. Bower

Abstract: Hot gaseous haloes surrounding galaxies and extending well beyond the distribution of stars are a ubiquitous prediction of galaxy formation scenarios. The haloes are believed to consist of gravitationally trapped gas with a temperature of millions of Kelvin. The existence of such hot haloes around massive elliptical galaxies has been established through their X-ray emission. While gas out-flowin… ▽ More Hot gaseous haloes surrounding galaxies and extending well beyond the distribution of stars are a ubiquitous prediction of galaxy formation scenarios. The haloes are believed to consist of gravitationally trapped gas with a temperature of millions of Kelvin. The existence of such hot haloes around massive elliptical galaxies has been established through their X-ray emission. While gas out-flowing from starburst spiral galaxies has been detected, searches for hot haloes around normal, quiescent spiral galaxies have so far failed, casting doubts on the fundamental physics in galaxy formation models. Here we present the first detection of a hot, large-scale gaseous halo surrounding a normal, quiescent spiral galaxy, NGC 5746, alleviating a long-standing problem for galaxy formation models. In contrast to starburst galaxies, where the X-ray halo can be powered by the supernova energy, there is no such power source in NGC 5746. The only compelling explanation is that we are here witnessing a galaxy forming from gradually in-flowing hot and dilute halo gas. △ Less

Submitted 13 January, 2006; v1 submitted 23 November, 2005; originally announced November 2005.

Comments: New Astronomy, in press

Journal ref: New Astron. 11 (2006) 465-470

Properties of Intra-group Stars and Galaxies in Galaxy Groups: "Normal" versus "Fossil" Groups

Authors: Jesper Sommer-Larsen

Abstract: Cosmological LCDM simulations of 12 M_vir~10^14 Msun galaxy groups have been performed, invoking star formation, chemical evolution with non-instantaneous recycling, metallicity dependent radiative cooling, strong star-burst driven galactic super-winds and effects of a meta-galactic UV field. At z=0, intra-group light (IGL) fractions are found to be 12-45%. Low values refer to groups with only a… ▽ More Cosmological LCDM simulations of 12 M_vir~10^14 Msun galaxy groups have been performed, invoking star formation, chemical evolution with non-instantaneous recycling, metallicity dependent radiative cooling, strong star-burst driven galactic super-winds and effects of a meta-galactic UV field. At z=0, intra-group light (IGL) fractions are found to be 12-45%. Low values refer to groups with only a small difference between the R-band magnitudes of the first and second ranked group galaxy, large are typical of "fossil" groups (FGs). The IG stars in the 4 FGs are 0.3-0.5 Gyr older than in the 8 nonFGs. For the IGL, B-R=~1.4, in good agreement with observations. For FGs/nonFGs the iron abundance of the IG stars is slightly sub-solar in the central parts (r~100 kpc) decreasing to about 40% solar at about 0.5 r_vir The IG stars are alpha-element enhanced with [O/Fe] increasing with r, and an overall [O/Fe]~0.45, indicating predominant SNII enrichment. The velocity distributions of the IG stars and group galaxies are, at r>~30 kpc, significantly more radially anisotropic for FGs than for nonFGs. So a characteristic of FG formation, apart from formation time (D'Onghia et al.), may be the "initial" velocity distribution of the group galaxies. For FGs one can dynamically infer the (dark matter dominated) mass distribution of the groups all the way to r_vir, from the kinematics of the IG stars or group galaxies. For the nonFGs this method overestimates the group mass at r>~200 kpc, by up to a factor of two at r_vir. This is interpreted as FGs being, in general, more relaxed than nonFGs. Finally, FGs of the above M_vir should host ~500 planetary nebulae at projected distances between 100 and 1000 kpc from the first ranked galaxy. All results appear consistent with the FG formation scenario of D'Onghia et al. △ Less

Submitted 27 March, 2006; v1 submitted 28 September, 2005; originally announced September 2005.

Comments: 12 pages, 15 figures, Accepted for MNRAS, Printing in colour recommended

Journal ref: Mon.Not.Roy.Astron.Soc.369:958-968,2006

Simulating galaxy clusters -- I. Thermal and chemical properties of the intra-cluster medium

Authors: A. D. Romeo, J. Sommer-Larsen, L. Portinari, V. Antonuccio-Delogu

Abstract: We have performed a series of N-body/hydrodynamical (TreeSPH) simulations of clusters and groups of galaxies, selected from cosmological N-body simulations within a $Λ$CDM framework: these objects have been re-simulated at higher resolution to $z$=0, in order to follow also the dynamical, thermal and chemical input on to the ICM from stellar populations within galaxies. The simulations include m… ▽ More We have performed a series of N-body/hydrodynamical (TreeSPH) simulations of clusters and groups of galaxies, selected from cosmological N-body simulations within a $Λ$CDM framework: these objects have been re-simulated at higher resolution to $z$=0, in order to follow also the dynamical, thermal and chemical input on to the ICM from stellar populations within galaxies. The simulations include metal dependent radiative cooling, star formation according to different IMFs, energy feedback as strong starburst-driven galactic super-winds, chemical evolution with non-instantaneous recycling of gas and heavy elements, effects of a meta-galactic UV field and thermal conduction in the ICM. In this Paper I of a series of three, we derive results, mainly at $z=0$, on the temperature and entropy profiles of the ICM, its X-ray luminosity, the cluster cold components (cold fraction as well as mass--to--light ratio) and the metal distribution between ICM and stars. In general, models with efficient super-winds, along with a top-heavy stellar IMF, are able to reproduce fairly well the observed $L_X-T$ relation, the entropy profiles and the cold fraction. Observed radial ICM temperature profiles can be matched, except for the gradual decline in temperature inside of $r\sim$~0.1$R_{\rm{vir}}$. Metal enrichment of the ICM gives rise to somewhat steep inner iron gradients; yet, the global level of enrichment compares well to observational estimates after correcting for the stars formed at late times at the base of the cooling flows; also the metal partition between stars and ICM gets into good agreement with observations. △ Less

Submitted 8 September, 2006; v1 submitted 16 September, 2005; originally announced September 2005.

Comments: 23 pages, 20 colour figures; final version accepted

Journal ref: Mon.Not.Roy.Astron.Soc.371:548-568,2006; Erratum-ibid.373:1648,2006

The Formation of Fossil Galaxy Groups in the hierarchical Universe

Authors: E. D'Onghia, J. Sommer-Larsen, A. D. Romeo, A. Burkert, K. Pedersen, L. Portinari, J. Rasmussen

Abstract: We use a set of twelve high-resolution N-body/hydrodynamical simulations in the $Λ$CDM cosmology to investigate the origin and formation rate of fossil groups (FGs), which are X-ray bright galaxy groups dominated by a large elliptical galaxy, with the second brightest galaxy being at least two magnitudes fainter. The simulations invoke star formation, chemical evolution with non-instantaneous re… ▽ More We use a set of twelve high-resolution N-body/hydrodynamical simulations in the $Λ$CDM cosmology to investigate the origin and formation rate of fossil groups (FGs), which are X-ray bright galaxy groups dominated by a large elliptical galaxy, with the second brightest galaxy being at least two magnitudes fainter. The simulations invoke star formation, chemical evolution with non-instantaneous recycling, metal dependent radiative cooling, strong star burst driven galactic super winds, effects of a meta-galactic UV field and full stellar population synthesis. We find an interesting correlation between the magnitude gap between the first and second brightest galaxy and the formation time of the group. It is found that FGs have assembled half of their final dark matter mass already at $z\ga1$, and subsequently typically grow by minor merging only, wheras non-FGs on average form later. The early assembly of FGs leaves sufficient time for galaxies of $L \sim L_*$ to merge into the central one by dynamical friction, resulting in the large magnitude gap at $z=0$. A fraction of 33$\pm$16% of the groups simulated are found to be fossil, whereas the observational estimate is $\sim$10-20%. The FGs are found to be X-ray over-luminous relative to non-FGs of the same optical luminosity, in qualitative agreement with observations. Finally, from a dynamical friction analysis is found that only because infall of $L \sim L_*$ galaxies happens along filaments with small impact parameters do FGs exist at all. △ Less

Submitted 31 August, 2005; v1 submitted 26 May, 2005; originally announced May 2005.

Comments: 4 pages, 3 figures, one figure removed. Accepted for publication in ApJ Letter

Journal ref: Astrophys.J. 630 (2005) L109-L112

Simulating galaxy clusters : the ICM and the galaxy populations

Authors: Alessio D. Romeo, J. Sommer-Larsen, L. Portinari

Abstract: Cosmological LambdaCDM TreeSPH simulations of the formation and evolution of galaxy groups and clusters have been performed. The simulations include: star formation, chemical evolution with non-instantaneous recycling, metal dependent radiative cooling, strong star burst and (optionally) AGN driven galactic super winds, effects of a meta-galactic UV field and thermal conduction. We report result… ▽ More Cosmological LambdaCDM TreeSPH simulations of the formation and evolution of galaxy groups and clusters have been performed. The simulations include: star formation, chemical evolution with non-instantaneous recycling, metal dependent radiative cooling, strong star burst and (optionally) AGN driven galactic super winds, effects of a meta-galactic UV field and thermal conduction. We report results on the temperature and entropy profiles of the ICM, the X-ray luminosity, cold fraction, M/L and IMLR ratios in gas and stars, metal abundances and gradients. Besides, the properties of the galaxy populations in the two richest clusters are discussed: global star formation rates of the cluster galaxies, the total K-band luminosity, the galaxy luminosity functions at z=0 and their redshift evolution, the colour-magnitude relation (``red sequence'') as resulting from metallicity effects, and the role of the IMF in reproducing colours and abundances of the stellar populations. Finally the contribute (20-40 %) to cluster light from the intra-cluster stars and the cD galaxy has been investigated: surface brightness profiles, mean colours and chemical abundances, kinematics (velocity dispersion and distributions) of the IC stellar populations. △ Less

Submitted 15 April, 2005; originally announced April 2005.

Comments: Conference proceedings, 8 pags., 7 figs

Journal ref: Proc. of Conf. "Baryons in Dark Matter Halos", Novigrad, 5-9 Oct. 2004, PoS SISSA

Simulating galaxy Clusters -II: global star formation histories and galaxy populations

Authors: Alessio D. Romeo, Laura Portinari, Jesper Sommer-Larsen

Abstract: Cosmological (LambdaCDM) TreeSPH simulations of the formation and evolution of galaxy groups and clusters have been performed. The simulations invoke star formation, chemical evolution with non-instantaneous recycling, metal dependent radiative cooling, strong star burst and (optionally) AGN driven galactic super winds, effects of a meta-galactic UV field and thermal conduction. The properties o… ▽ More Cosmological (LambdaCDM) TreeSPH simulations of the formation and evolution of galaxy groups and clusters have been performed. The simulations invoke star formation, chemical evolution with non-instantaneous recycling, metal dependent radiative cooling, strong star burst and (optionally) AGN driven galactic super winds, effects of a meta-galactic UV field and thermal conduction. The properties of the galaxy populations in two clusters, one Virgo-like (T~3 keV) and one (sub) Coma-like (T~6 keV), are discussed. The global star formation rates of the cluster galaxies are found to decrease very significantly with time from redshift z=2 to 0, in agreement with observations. The total K-band luminosity of the cluster galaxies correlates tightly with total cluster mass, and for models without additional AGN feedback, the zero point of the relation matches the observed one fairly well. The match to observed galaxy luminosity functions is reasonable, except for a deficiency of bright galaxies (M_B < -20), which becomes increasingly significant with super-wind strength. Results of a high resolution test indicate that this deficiency is not due to ``over--merging''. The redshift evolution of the luminosity functions from z=1 to 0 is mainly driven by luminosity evolution, but also by merging of bright galaxies with the cD. The colour--magnitude relation of the cluster galaxies matches the observed "red sequence" very well and, on average, galaxy metallicity increases with luminosity. As the brighter galaxies are essentially coeval, the colour--magnitude relation results from metallicity rather than age effects, as observed. △ Less

Submitted 1 June, 2005; v1 submitted 22 April, 2004; originally announced April 2004.

Comments: 15 pages, 10 figures. Final version accepted by MNRAS, presenting new simulations and major changes. Printing in colour recommended

Journal ref: Mon.Not.Roy.Astron.Soc. 361 (2005) 983-996

Simulating Galaxies and Galaxy Clusters

Authors: Jesper Sommer-Larsen

Abstract: The X-ray luminosities of the hot halo gas around simulated, Milky Way like disk galaxies have been determined, as a function of redshift. The X-ray luminosity increases significantly with redshift, in some cases as much as a factor 30 going from z=0 to 2. Consequently, the optimal detection redshift can be >~0.5. Results of fully cosmological simulations of galaxy groups and clusters, incorpo… ▽ More The X-ray luminosities of the hot halo gas around simulated, Milky Way like disk galaxies have been determined, as a function of redshift. The X-ray luminosity increases significantly with redshift, in some cases as much as a factor 30 going from z=0 to 2. Consequently, the optimal detection redshift can be >~0.5. Results of fully cosmological simulations of galaxy groups and clusters, incorporating star formation, chemical evolution with non-instantaneous recycling, metal-dependent radiative cooling, galactic super-winds and thermal conduction are presented. X-ray luminosities at z=0 are somewhat high and central entropies somewhat low compared to observations. This is likely a combined effect of chemical evolution and metal-dependent, radiative cooling. Central ICM abundance profiles are somewhat steep, and the observed level of ICM enrichment can only be reproduced with IMFs more top-heavy than the Salpeter IMF. In agreement with observations it is found that the iron in the ICM is in place already at z~1. The [Si/Fe] of the ICM decreases with time, and increases slightly with radius at a given time. The cluster galaxies match the observed "red sequence" very well, and the metallicity of cluster galaxies increases with galaxy mass, as observed. △ Less

Submitted 1 April, 2004; originally announced April 2004.

Comments: 6 pages, 2 B/W + 5 colour figures. Proceedings of the Vulcano conference on "Modelling the intergalactic and intracluster media", V. Antonuccio et al. (eds.), to be published online on MemSAIt - Supplementi. Printing in colour recommended

Report number: Nordita preprint #25 (2004)

Simulating Galaxy Clusters - III: Properties of the Intracluster Stars

Authors: J. Sommer-Larsen, A. D. Romeo, L. Portinari

Abstract: Cosmological (LCDM) TreeSPH simulations of galaxy clusters have been performed, invoking star-formation, chemical evolution with non-instantaneous recycling, metal dependent radiative cooling, strong star-burst and/or AGN driven super-winds, effects of a meta-galactic UV field and thermal conduction. Results for two clusters, one Virgo-like (T~3 keV) and one (sub) Coma-like (T~6 keV), are presen… ▽ More Cosmological (LCDM) TreeSPH simulations of galaxy clusters have been performed, invoking star-formation, chemical evolution with non-instantaneous recycling, metal dependent radiative cooling, strong star-burst and/or AGN driven super-winds, effects of a meta-galactic UV field and thermal conduction. Results for two clusters, one Virgo-like (T~3 keV) and one (sub) Coma-like (T~6 keV), are presented. At z=0, the stellar contents of both clusters consist of a central dominant (cD) galaxy surrounded by galaxies and intracluster (IC) stars. The IC stars are found to contribute 20-40% of the total cluster B-band light and to form at a mean redshift z_f~3, on average about 0.5 Gyr prior to the stars in galaxies. UBVRIJHK surface brightness profiles of the IC star populations are presented; the profile of the larger cluster matches the observed V-band profile of the cD in Abell 1413 (T~8 keV). The IC stellar populations have B-R=1.4-1.5, comparable to larger E galaxies. The mean [Fe/H]~0.0 in the central part of the cluster (r~100 kpc) decreasing to about -0.3 at r_vir. The IC stars are alpha-element enhanced, with a weak trend of [O/Fe] increasing with r, and an overall [O/Fe]~0.4 dex. The IC stars are kinematically colder than the galaxies: The velocity dispersions of the IC stars are at r~100-500 kpc only about half of those of the galaxies increasing to about 70% at r=1-2 Mpc. The projected velocity dispersion in the Virgo-like cluster at R>~50 kpc is 300-600 km/s depending on orientation and R. Rotation is dynamically insignificant for the IC stars. The velocity distributions of IC stars and galaxies are in one cluster highly radially anisotropic, in the other close to isotropic. △ Less

Submitted 4 March, 2006; v1 submitted 11 March, 2004; originally announced March 2004.

Comments: 12 pages, 14 figures, Final version accepted by MNRAS, Printing in colour recommended

Journal ref: Mon.Not.Roy.Astron.Soc. 357 (2005) 478-488

Mass-to-Light ratio, Initial Mass Function and chemical evolution in disc galaxies

Authors: L. Portinari, J. Sommer-Larsen, R. Tantalo

Abstract: Cosmological simulations of disc galaxy formation, when compared to the observed Tully-Fisher relation, suggest a low Mass-to-Light (M/L) ratio for the stellar component in spirals. We show that a number of "bottom-light" Initial Mass Functions (IMFs) suggested independently in literature, do imply M/L ratios as low as required, at least for late type spirals (Sbc/Sc). However the typical M/L ra… ▽ More Cosmological simulations of disc galaxy formation, when compared to the observed Tully-Fisher relation, suggest a low Mass-to-Light (M/L) ratio for the stellar component in spirals. We show that a number of "bottom-light" Initial Mass Functions (IMFs) suggested independently in literature, do imply M/L ratios as low as required, at least for late type spirals (Sbc/Sc). However the typical M/L ratio, and correspondingly the zero-point of the Tully-Fisher relation, is expected to vary considerably with Hubble type. Bottom-light IMFs tend to have a metal production in excess of what is tipically estimated for spiral galaxies. Suitable tuning of the IMF slope and mass limits, post-supernova fallback of metals onto black holes or metal outflows must then be invoked, to reproduce the observed chemical properties of disc galaxies. △ Less

Submitted 18 February, 2004; originally announced February 2004.

Comments: 4 pages, 6 colour figures. To appear on PASA (refereed proceedings of Galactic ChemoDynamics V, B.K. Gibson and D. Kawata eds.)