Search

We investigate the impact of AGN feedback, on the entropy and characteristic temperature measurements of galaxy groups detected in the SRG/eROSITA's first All-Sky Survey (eRASS1) to shed light on the characteristics of the feedback mechanisms. We analyze deeper eROSITA observations of 1178 galaxy groups detected in eRASS1. We divide the sample into 271 subsamples and extract average thermodynamic properties, including electron density, temperature, and entropy at three characteristic radii along with the integrated temperature by jointly analyzing X-ray images and spectra following a Bayesian approach. We present the tightest constraints on the impact of AGN feedback through our average entropy and characteristic temperature measurements of the largest group sample used in X-ray studies, incorporating major systematics in our analysis. We find that entropy shows an increasing trend with temperature in the form of a power-law-like relation at the higher intra-group medium temperatures, while for the low mass groups, a slight flattening is observed on the average entropy. Overall, the observed entropy measurements agree well with the earlier measurements in the literature. The comparisons with the state-of-the-art cosmological hydrodynamic simulations (MillenniumTNG, Magneticum, OWL simulations) after the applications of the selection function calibrated for our galaxy groups reveal that observed entropy profiles in the cores are below the predictions of simulations. At the mid-region, the entropy measurements agree well with the Magneticum simulations, whereas the predictions of MillenniumTNG and OWL simulations fall below observations. At the outskirts, the overall agreement between the observations and simulations improves, with Magneticum simulations reproducing the observations the best. Our measurements will pave the way for more realistic AGN feedback implementations in simulations., Comment: 25 pages, 16 figures, 3 tables, Submitted to A&A

New observational facilities are probing astrophysical transients such as stellar explosions and gravitational wave (GW) sources at ever increasing redshifts, while also revealing new features in source property distributions. To interpret these observations, we need to compare them to predictions from stellar population models. Such models require the metallicity-dependent cosmic star formation history ($\mathcal{S}(Z,z)$) as an input. Large uncertainties remain in the shape and evolution of this function. In this work, we propose a simple analytical function for $\mathcal{S}(Z,z)$. Variations of this function can be easily interpreted, because the parameters link to its shape in an intuitive way. We fit our analytical function to the star-forming gas of the cosmological TNG100 simulation and find that it is able to capture the main behaviour well. As an example application, we investigate the effect of systematic variations in the $\mathcal{S}(Z,z)$ parameters on the predicted mass distribution of locally merging binary black holes (BBH). Our main findings are: I) the locations of features are remarkably robust against variations in the metallicity-dependent cosmic star formation history, and II) the low mass end is least affected by these variations. This is promising as it increases our chances to constrain the physics that governs the formation of these objects., Comment: Accepted for publication in ApJ, 6 Figures, code available at https://github.com/LiekeVanSon/SFRD_fit, made using showyourwork

The visibility of high-redshift Lyman-alpha emitting galaxies (LAEs) provides important constraints on galaxy formation processes and the Epoch of Reionization (EoR). However, predicting realistic and representative statistics for comparison with observations represents a significant challenge in the context of large-volume cosmological simulations. The THESAN project offers a unique framework for addressing such limitations by combining state-of-the-art galaxy formation (IllustrisTNG) and dust models with the Arepo-RT radiation-magneto-hydrodynamics solver. In this initial study we present Lyman-alpha centric analysis for the flagship simulation that resolves atomic cooling haloes throughout a (95.5 cMpc)^3 region of the Universe. To avoid numerical artifacts we devise a novel method for accurate frequency-dependent line radiative transfer in the presence of continuous Hubble flow, transferable to broader astrophysical applications as well. Our scalable approach highlights the utility of LAEs and red damping-wing transmission as probes of reionization, which reveal nontrivial trends across different galaxies, sightlines, and frequency bands that can be modelled in the framework of covering fractions. In fact, after accounting for environmental factors influencing large-scale ionized bubble formation such as redshift and UV magnitude, the variation across galaxies and sightlines mainly depends on random processes including peculiar velocities and self-shielded systems that strongly impact unfortunate rays more than others. Throughout the EoR local and cosmological optical depths are often greater than or less than unity such that the exp(-tau) behavior leads to anisotropic and bimodal transmissivity. Future surveys will benefit by targeting both rare bright objects and Goldilocks zone LAEs to infer the presence of these (un)predictable (dis)advantages., Comment: 23 pages, 23 figures, MNRAS, in press. Please visit www.thesan-project.com for more details

We introduce the THESAN project, a suite of large volume (L = 95.5 cMpc) radiation-magnetohydrodynamic simulations that simultaneously model the large-scale statistical properties of the IGM during reionization and the resolved characteristics of the galaxies responsible for it. The flagship simulation has dark matter (DM) and baryonic mass resolutions of $3.1 \times 10^6$ $M_\odot$ and $5.8 \times 10^5 $ $M_\odot$, respectively. The gravitational forces are softened on scales of 2.2 ckpc with the smallest cell sizes reaching 10 pc at z=5.5, enabling predictions down to the atomic cooling limit. The simulations use an efficient radiation hydrodynamics solver (AREPO-RT) that precisely captures the interaction between ionizing photons and gas, coupled to well-tested galaxy formation (IllustrisTNG) and dust models to accurately predict the properties of galaxies. Through a complementary set of medium resolution simulations we investigate the changes to reionization introduced by different assumptions for ionizing escape fractions, varying DM models, and numerical convergence. The simulations produce realistic reionization histories that match the observed evolution of the global neutral hydrogen fraction and electron scattering optical depth to reionization. They also match a wealth of high-redshift observationally inferred data, including the stellar-to-halo-mass relation, stellar mass function, star formation rate density, and the mass-metallicity relation, despite the galaxy formation model being mainly calibrated at z=0. We demonstrate that different reionization models give rise to varied bubble size distributions that imprint unique signatures on the 21cm emission, especially on the slope of the power spectrum at large spatial scales, enabling current and upcoming 21cm experiments to accurately characterise the sources that dominate the ionizing photon budget. [abridged], Comment: MNRAS, in press. Please visit www.thesan-project.com for more details

Stellar halos in early-type galaxies (ETGs) are shaped by their accretion and merger histories. We use a sample of 1114 ETGs in the TNG100 simulation with stellar masses $10^{10.3}\leq M_{*}/M_\odot\leq 10^{12}$, selected at z=0 within the range of g-r colour and lambda-ellipticity diagram populated by observed ETGs. We study how the rotational support and intrinsic shapes of the stellar halos depend on the fraction of accreted stars, overall and separately by major, minor, and mini mergers. Accretion histories in TNG100 ETGs as well as the radial distributions of ex-situ stars $f_{ex}(R)$ strongly correlate with stellar mass. Low-mass ETGs have characteristic peaked rotation profiles and near-oblate shapes with rounder halos that are completely driven by the in-situ stars. At high $f_{ex}$ major mergers decrease the in-situ peak in rotation velocity, flatten the $V_{*}/\sigma_{*}(R)$ profiles, and increase the triaxiality of the stellar halos. Kinematic transition radii do not trace the transition between in-situ and ex-situ dominated regions, but for systems with $M_{*}>10^{10.6}M_\odot$ the local rotational support decreases with the local ex-situ fraction $f_{ex}(R)$ and triaxiality increases with $f_{ex}$. These correlations are followed by fast and slow rotators alike with a continuous and overlapping sequence of properties. Merger events dynamically couple stars and dark matter: in high mass ETGs and at large radii where $f_{ex}\gtrsim0.5$, both components tend to have similar intrinsic shapes and rotational support, and nearly aligned principal axes and spin directions. Based on these results we suggest that extended photometry and kinematics of massive ETGs ($M_{*}>10^{10.6}M_\odot$) can be used to estimate the local fraction of ex-situ stars and to approximate the intrinsic shapes and rotational support of the co-spatial dark matter component. [abridged], Comment: 22 pages, 17 figures, accepted for publication on A&A update version including a new section on high-redshift compact progenitor galaxies and a new Appendix with with tables reporting relevant physical relations derived throughout this paper

We characterize the photometric and kinematic properties of simulated early-type galaxy (ETG) stellar halos, and compare them to observations. We select a sample of ~1200 ETGs in the TNG100 and TNG50 simulations, spanning a stellar mass range of $10^{10.3}-10^{12}M_{\odot}$ and within the range of (g-r) colour and lambda-ellipticity diagram populated by observed ETGs. We determine photometric parameters, intrinsic shapes, and kinematic observables in their extended stellar halos. We study the variation in kinematics from center to halo and connect it to a change in the intrinsic shape of the galaxies. We find that the simulated galaxy sample reproduces the diversity of kinematic properties observed in ETG halos. Simulated fast rotators (FRs) divide almost evenly in one third having flat lambda profiles and high halo rotational support, a third with gently decreasing profiles, and another third with low halo rotation. Slow rotators (SRs) tend to have increased rotation in the outskirts, with half of them exceeding lambda=0.2. For $M_{*}>10^{11.5}M_{\odot}$ halo rotation is unimportant. A similar variety of properties is found for the stellar halo intrinsic shapes. Rotational support and shape are deeply related: the kinematic transition to lower rotational support is accompanied by a change towards rounder intrinsic shape. Triaxiality in the halos of FRs increases outwards and with stellar mass. Simulated SRs have relatively constant triaxiality profiles. Simulated stellar halos show a large variety of structural properties, with quantitative but no clear qualitative differences between FRs and SRs. At the same stellar mass, stellar halo properties show a gradual transition and significant overlap between the two families, despite the clear bimodality in the central regions. This is in agreement with observations of extended photometry and kinematics. [abridged], Comment: accepted for publication in A&A, 25 pages, 22 figures

Over the past century, major advances in astronomy and astrophysics have been largely driven by improvements in instrumentation and data collection. With the amassing of high quality data from new telescopes, and especially with the advent of deep and large astronomical surveys, it is becoming clear that future advances will also rely heavily on how those data are analyzed and interpreted. New methodologies derived from advances in statistics, computer science, and machine learning are beginning to be employed in sophisticated investigations that are not only bringing forth new discoveries, but are placing them on a solid footing. Progress in wide-field sky surveys, interferometric imaging, precision cosmology, exoplanet detection and characterization, and many subfields of stellar, Galactic and extragalactic astronomy, has resulted in complex data analysis challenges that must be solved to perform scientific inference. Research in astrostatistics and astroinformatics will be necessary to develop the state-of-the-art methodology needed in astronomy. Overcoming these challenges requires dedicated, interdisciplinary research. We recommend: (1) increasing funding for interdisciplinary projects in astrostatistics and astroinformatics; (2) dedicating space and time at conferences for interdisciplinary research and promotion; (3) developing sustainable funding for long-term astrostatisics appointments; and (4) funding infrastructure development for data archives and archive support, state-of-the-art algorithms, and efficient computing., Comment: Submitted to the Astro2020 Decadal Survey call for science white papers

We write in response to the call from the 2020 Decadal Survey to submit white papers illustrating the most pressing scientific questions in astrophysics for the coming decade. We propose exploration as the central question for the Decadal Committee's discussions.The history of astronomy shows that paradigm changing discoveries are not driven by well formulated scientific questions, based on the knowledge of the time. They were instead the result of the increase in discovery space fostered by new telescopes and instruments. An additional tool for increasing the discovery space is provided by the analysis and mining of the increasingly larger amount of archival data available to astronomers. Revolutionary observing facilities, and the state of the art astronomy archives needed to support these facilities, will open up the universe to new discovery. Here we focus on exploration for compact objects and multi messenger science. This white paper includes science examples of the power of the discovery approach, encompassing all the areas of astrophysics covered by the 2020 Decadal Survey.

We present a machine-learning approach for estimating galaxy cluster masses from Chandra mock images. We utilize a Convolutional Neural Network (CNN), a deep machine learning tool commonly used in image recognition tasks. The CNN is trained and tested on our sample of 7,896 Chandra X-ray mock observations, which are based on 329 massive clusters from the IllustrisTNG simulation. Our CNN learns from a low resolution spatial distribution of photon counts and does not use spectral information. Despite our simplifying assumption to neglect spectral information, the resulting mass values estimated by the CNN exhibit small bias in comparison to the true masses of the simulated clusters (-0.02 dex) and reproduce the cluster masses with low intrinsic scatter, 8% in our best fold and 12% averaging over all. In contrast, a more standard core-excised luminosity method achieves 15-18% scatter. We interpret the results with an approach inspired by Google DeepDream and find that the CNN ignores the central regions of clusters, which are known to have high scatter with mass., Comment: 10 pages, 6 figures, accepted for publication in The Astrophysical Journal

Long dynamical timescales in the outskirts of galaxies preserve the information content of their accretion histories, for example in the form of stellar population gradients. We present a detailed analysis of the stellar halo properties of a statistically representative sample of early-type galaxies from the Illustris simulation and show that stellar population gradients at large radii can indeed be used to infer basic properties of galactic accretion histories. We measure metallicity, age, and surface-brightness profiles in quiescent Illustris galaxies ranging from $\mathrm{M}_\star = 10^{10} - 2\times 10^{12}\;\mathrm{M}_\odot$ and show that they are in reasonable agreement with observations. At fixed mass, galaxies that accreted little of their stellar halo material tend to have steeper metallicity and surface-brightness profiles between 2 - 4 effective radii (R$_e$) than those with larger accreted fractions. Profiles of metallicity and surface-brightness in the stellar halo typically flatten from z = 1 to the present. This suggests that the accretion of stars into the stellar halo tends to flatten metallicity and surface-brightness profiles, a picture which is supported by the tight correlation between the two gradients in the stellar halo. We find no statistical evidence of additional information content related to accretion histories in stellar halo metallicity profiles beyond what is contained in surface-brightness profiles. Age gradients in the stellar halo do not appear to be sensitive to galactic accretion histories, and none of the stellar population gradients studied are strongly correlated with the mean merger mass-ratio. Future observations that reach large radii outside galaxies will have the best potential to constrain galactic accretion histories., Comment: 18 pages, 12 figures. Resubmitted to the ApJ after referee comments

We use a set of high-resolution N-body simulations of the Galactic disk to study its interactions with the population of satellites predicted cosmologically. One simulation illustrates that multiple passages of massive satellites with different velocities through the disk generate a wobble, having the appearance of rings in face-on projections of the stellar disk. They also produce flares in the disk outer parts and gradually heat the disk through bending waves. A different numerical experiment shows that an individual satellite as massive as the Sagittarius dwarf galaxy passing through the disk will drive coupled horizontal and vertical oscillations of stars in underdense regions, with small significant associated heating. This experiment shows that vertical excursions of stars in these low-density regions can exceed 1 kpc in the Solar neighborhood, resembling the coherent vertical oscillations recently detected locally. They can also induce non-zero vertical streaming motions as large as 10-20 km s$^{-1}$, consistent with recent observations in the Galactic disk. This phenomenon appears as a local ring, with modest associated disk heating., Comment: Accepted for publication on Astrophysical Journal, minor revision, 1 table added

The influence of both major and minor mergers is expected to significantly affect gradients of stellar ages and metallicities in the outskirts of galaxies. Measurements of observed gradients are beginning to reach large radii in galaxies, but a theoretical framework for connecting the findings to a picture of galactic build-up is still in its infancy. We analyze stellar populations of a statistically representative sample of quiescent galaxies over a wide mass range from the Illustris simulation. We measure metallicity and age profiles in the stellar halos of quiescent Illustris galaxies ranging in stellar mass from $10^{10}$ to $10^{12} M_\odot$, accounting for observational projection and luminosity-weighting effects. We find wide variance in stellar population gradients between galaxies of similar mass, with typical gradients agreeing with observed galaxies. We show that, at fixed mass, the fraction of stars born in-situ within galaxies is correlated with the metallicity gradient in the halo, confirming that stellar halos contain unique information about the build-up and merger histories of galaxies., Comment: 6 pages, 2 figures, proceedings of IAU Symposium 317: "The General Assembly of Galaxy Halos: Structure, Origin and Evolution"

We confirm our recent prediction of the "pitchfork" foreground signature in power spectra of high-redshift 21 cm measurements where the interferometer is sensitive to large-scale structure on all baselines. This is due to the inherent response of a wide-field instrument and is characterized by enhanced power from foreground emission in Fourier modes adjacent to those considered to be the most sensitive to the cosmological H I signal. In our recent paper, many signatures from the simulation that predicted this feature were validated against Murchison Widefield Array (MWA) data, but this key pitchfork signature was close to the noise level. In this paper, we improve the data sensitivity through the coherent averaging of 12 independent snapshots with identical instrument settings and provide the first confirmation of the prediction with a signal-to-noise ratio > 10. This wide-field effect can be mitigated by careful antenna designs that suppress sensitivity near the horizon. Simple models for antenna apertures that have been proposed for future instruments such as the Hydrogen Epoch of Reionization Array and the Square Kilometre Array indicate they should suppress foreground leakage from the pitchfork by ~40 dB relative to the MWA and significantly increase the likelihood of cosmological signal detection in these critical Fourier modes in the three-dimensional power spectrum., Comment: Accepted for publication in ApJ Letters, 6 pages, 2 figures

The separation of the faint cosmological background signal from bright astrophysical foregrounds remains one of the most daunting challenges of mapping the high-redshift intergalactic medium with the redshifted 21 cm line of neutral hydrogen. Advances in mapping and modeling of diffuse and point source foregrounds have improved subtraction accuracy, but no subtraction scheme is perfect. Precisely quantifying the errors and error correlations due to missubtracted foregrounds allows for both the rigorous analysis of the 21 cm power spectrum and for the maximal isolation of the "EoR window" from foreground contamination. We present a method to infer the covariance of foreground residuals from the data itself in contrast to previous attempts at a priori modeling. We demonstrate our method by setting limits on the power spectrum using a 3 h integration from the 128-tile Murchison Widefield Array. Observing between 167 and 198 MHz, we find at 95% confidence a best limit of Delta^2(k) < 3.7 x 10^4 mK^2 at comoving scale k = 0.18 hMpc^-1 and at z = 6.8, consistent with existing limits., Comment: 18 pages, 9 figures. Accepted to Phys. Rev. D

Recent studies have shown that galaxies accrete most of their baryons via the cold mode, from streams with temperatures T~10^4-10^5 K. At these temperatures, the streams should radiate primarily in the Lya line and have therefore been proposed as a model to power the extended, high-redshift objects known as Lya blobs and other high-redshift Lya sources. We introduce a new Lya radiative transfer code, aRT, and apply it to cosmological hydrodynamical simulations. We address physical and numerical issues that are critical to making accurate predictions for the cooling luminosity, but that have been mostly neglected or treated simplistically so far. We highlight the importance of self-shielding and of properly treating sub-resolution models in simulations. Most existing simulations do not self-consistently incorporate these effects, which can lead to order-of-magnitude errors in the predicted cooling luminosity. Using a combination of post-processing ionizing radiative transfer and re-simulation techniques, we develop an approximation to the consistent evolution of the self-shielded gas. We quantify the dependence of the Lya cooling luminosity on halo mass at z=3 for the simplified problem of pure gas accretion. While cooling in massive halos (without additional energy input from star formation and AGN) is in principle sufficient to produce L_alpha~10^43-10^44 erg s^-1 blobs, this requires including energy released in gas of density sufficient to form stars, but which is kept 100% gaseous in our optimistic estimates. Excluding emission from such dense gas yields lower luminosities by up to one to two orders of magnitude at high masses, making it difficult to explain the observed Lya blobs with pure cooling. Resonant scattering produces diffuse Lya halos, even for centrally concentrated emission, and broad double peaked line profiles. [Abridged], Comment: 29 pages, 13 figures. ApJ, in press

We present a merger-driven evolutionary model for the production of luminous (LIRGs) and ultraluminous infrared galaxies (ULIRGs) with warm IR colours. Our results show that simulations of gas--rich major mergers including star formation, black hole growth, and feedback can produce warm (U)LIRGs. We also find that while the warm evolutionary phase is associated with increased AGN activity, star formation alone may be sufficient to produce warm IR colours. However, the transition can be suppressed entirely - even when there is a significant AGN contribution - when we assume a single-phase ISM, which maximizes the attenuation. Finally, our evolutionary models are consistent with the 25-to-60 micron flux density ratio vs. L_HX/L_IR relation for local LIRGs and ULIRGs, and predict the observed scatter in IR color at fixed L_HX/L_IR. Therefore, our models suggest a cautionary note in the interpretation of warm IR colours: while associated with periods of active black hole growth, they are probably produced by a complex mix of star formation and AGN activity intermediate between the cold star formation dominated phase and the birth of a bright, unobscured quasar., Comment: 5 pages, accepted for publication in MNRAS Letters

The ionizing background determines the ionization balance and the thermodynamics of the cosmic gas. It is therefore a fundamental ingredient to theoretical and empirical studies of both the IGM and galaxy formation. We present here a new calculation of its spectrum that satisfies the empirical constraints we recently obtained by combining state-of-the-art luminosity functions and intergalactic opacity measurements. In our preferred model, star-forming galaxies and quasars each contribute substantially to the HI ionizing field at z<3, with galaxies rapidly overtaking quasars at higher redshifts as quasars become rarer. In addition to our fiducial model, we explore the physical dependences of the calculated background and clarify how recombination emission contributes to the ionization rates. We find that recombinations do not simply boost the ionization rates by the number of reemitted ionizing photons as many of these rapidly redshift below the ionization edges and have a distribution of energies. A simple analytic model that captures the main effects seen in our numerical radiative transfer calculations is given. Finally, we discuss the effects of HeII reionization by quasars on both the spectrum of the ionizing background and on the thermal history of the IGM. In regions that have yet to be reionized, the spectrum is expected to be almost completely suppressed immediately above 54.4 eV while a background of higher-energy (>~0.5 keV) photons permeates the entire universe owing to the frequency-dependence of the photoionization cross section. We provide an analytic model of the heat input during HeII reionization and its effects on the temperature-density relation., Comment: 36 pages, including 16 figures, accepted by ApJ. Data in electronic form can be obtained at http://www.cfa.harvard.edu/~cgiguere/uvbkg.html

We investigate the implications of the intergalactic opacity for the evolution of the cosmic UV luminosity density and its sources. Our main constraint is our measurement of the Lya forest opacity at redshifts 2~3. Combined with direct measurements of the galaxy UV luminosity function, this requires only a small fraction f_esc~0.5% of galactic hydrogen ionizing photons to escape their source for galaxies to solely account for the entire ionizing background. Under the assumption that the galactic UV emissivity traces the star formation rate, current state-of-the-art observational estimates of the star formation rate density appear to underestimate the total photoionization rate at z~4 by a factor ~4, are in tension with recent determinations of the UV luminosity function, and fail to reionize the Universe by z~6 if extrapolated to arbitrarily high redshift. A theoretical star formation history peaking earlier fits the Lya forest photoionization rate well, reionizes the Universe in time, and is in better agreement with the rate of z~4 gamma-ray bursts observed by Swift. Quasars suffice to doubly ionize helium by z~3 and likely contribute a non-negligible and perhaps dominant fraction of the hydrogen ionizing background at their z~2 peak. [Abridged], Comment: 26 pages, including 14 figures, submitted to ApJ and revised following the referee's report

We measure the width of the MgII $\lambda2799$ line in quasar spectra from the SDSS, 2QZ and 2SLAQ surveys and, by invoking an unnormalised virial mass estimator, relate the scatter in line width to the scatter in mass in the underlying black hole population. We find conclusive evidence for a trend such that there is less scatter in line width, and hence black hole mass, in more luminous objects. However, the most luminous objects in our sample show such a low degree of scatter in line width that, when combined with measures for the intrinsic scatter in the radius-luminosity relation for the broad-line region in active galaxies, an inconsistency arises in the virial technique for estimating black hole masses. This analysis implies that, at least for the most luminous quasars, either there is little-to-no intrinsic scatter in the radius-luminosity relation or the MgII broad emission line region is not totally dominated by virial velocities. Finally we exploit the measured scatter in line widths to constrain models for the velocity field of the broad-line region. We show that the lack of scatter in broad line-widths for luminous quasars is inconsistent with a pure planar/disk-like geometry for the broad-line region... (abridged), Comment: 19 pages 16 figures, accepted MNRAS. High res version available at: http://www.physics.usyd.edu.au/~sfine/qso_bl_scatter.ps

We investigate the implications of our measurement of the Lyman-alpha forest opacity at redshifts 2~3. Our measurement argues against a star formation rate density declining beyond z~3, in contrast with existing state-of-the-art determinations of the cosmic star formation history from direct galaxy counts. Stellar emission from galaxies therefore likely reionized the Universe., Comment: 5 pages, including 1 figure, published by ApJ

The universe is permeated by a network of filaments, sheets, and knots collectively forming a "cosmic web.'' The discovery of the cosmic web, especially through its signature of absorption of light from distant sources by neutral hydrogen in the intergalactic medium, exemplifies the interplay between theory and experiment that drives science, and is one of the great examples in which numerical simulations have played a key and decisive role. We recount the milestones in our understanding of cosmic structure, summarize its impact on astronomy, cosmology, and physics, and look ahead by outlining the challenges faced as we prepare to probe the cosmic web at new wavelengths., Comment: 10 pages, 2 figures. Appeared as a solicited Perspective article in the January 4, 2008 special issue of Science on the cosmic web

We directly measure the evolution of the intergalactic Lyman-alpha effective optical depth, tau_eff, over the redshift range 2 is <1% at z=2, 4% at z=3, and 12% at z=4. We provide estimates of the level of absorption arising from metals in the Ly-alpha forest based on both direct and statistical metal removal results in the literature, finding that this contribution is ~6-9% at z=3 and decreases monotonically with redshift. The high precision of our measurement indicates significant departures from the best-fit power-law redshift evolution, particularly near z=3.2., Comment: 4 pages, including 2 figures, to appear in the proceedings of "A Century of Cosmology", August 27-31, 2007, San Servolo, Venice, Italy; for full details, see arXiv:0709.2382

We directly measure the evolution of the intergalactic Lya effective optical depth, tau_eff, over the redshift range 2~3 based on directly fitting the quasar continua in the Lya forest, which have generally neglected this effect, are therefore likely biased low. We provide estimates of the level of absorption arising from metals in the Lya forest based on both direct and statistical metal removal results in the literature, finding that this contribution is ~6-9% at z=3 and decreases monotonically with redshift. The high precision of our measurement, attaining 3% in redshift bins of width Delta z=0.2 around z=3, indicates significant departures from the best-fit power-law redshift evolution (tau_eff=0.0018(1+z)^3.92, when metals are left in), particularly near z=3.2. The observed downward departure is statistically consistent with a similar feature detected in a precision statistical measurement using Sloan Digital Sky Survey spectra by Bernardi and coworkers using an independent approach., Comment: 27 pages, including 18 figures, published in ApJ

We show that the Lyman-alpha optical depth statistics in the proximity regions of quasar spectra depend on the mass of the dark matter halos hosting the quasars. This is owing to both the overdensity around the quasars and the associated infall of gas toward them. For a fiducial quasar host halo mass of (3.0+/-1.6) h^-1 x 10^12 Msun, as inferred by Croom et al. from clustering in the 2dF QSO Redshift Survey, we show that estimates of the ionizing background (Gamma^bkg) from proximity effect measurements could be biased high by a factor of ~2.5 at z=3 owing to neglecting these effects alone. The clustering of galaxies and other active galactic nuclei around the proximity effect quasars enhances the local background, but is not expected to skew measurements by more than a few percent. Assuming the measurements of Gamma^bkg based on the mean flux decrement in the Ly-alpha forest to be free of bias, we demonstrate how the proximity effect analysis can be inverted to measure the mass of the dark matter halos hosting quasars. In ideal conditions, such a measurement could be made with a precision comparable to the best clustering constraints to date from a modest sample of only about 100 spectra. We discuss observational difficulties, including continuum flux estimation, quasar systematic redshift determination, and quasar variability, which make accurate proximity effect measurements challenging in practice. These are also likely to contribute to the discrepancies between existing proximity effect and flux decrement measurements of Gamma^bkg., Comment: 25 pages, including 14 figures, accepted by ApJ

(ABRIDGED) We examine the fundamental scaling relations of elliptical galaxies formed through mergers. Using hundreds of simulations to judge the impact of progenitor galaxy properties on merger remnants, we find that gas dissipation provides an important contribution to tilt in the Fundamental Plane relation. Dissipationless mergers of disks produce remnants that occupy the virial plane. As the gas content of disk galaxies is increased, the tilt of the Fundamental Plane relation increases and the slope of the Re-M_* relation steepens. For gas fractions fgas > 30%, the simulated Fundamental Plane scalings approach those observed in the K-band. In our simulations, feedback from supermassive black hole growth has only a minor influence on the stellar-mass scaling relations of spheroidal galaxies, but may play a role in maintaining the observed Fundamental Plane tilt at optical wavelengths by suppressing residual star formation in merger remnants. We estimate that \approx 40-100% of the Fundamental Plane tilt induced by structural properties owes to trends in the central total-to-stellar mass ratio M_total/M_* produced by dissipation. Lower mass systems obtain greater phase- space densities than higher mass systems, producing a galaxy mass-dependent central M_total/M_* and a corresponding tilt in the Fundamental Plane., Comment: Version accepted by ApJ, 20 pages, 18 figures, resolution reduced for size

We study the physical and photometric properties of galaxies at z=4 in cosmological hydrodynamic simulations of a lambda-CDM universe. We focus on galaxies satisfying the GOODS "B-dropout" criteria. Our goals are: (1) to study the nature of high-redshift galaxies; (2) to test the simulations against published measurements of high-redshift galaxies; (3) to find relations between photometric measurements by HST/ACS (0.4 -- 1 micron) and Spitzer/IRAC (3.6 -- 8 micron) and the intrinsic physical properties of GOODS "B-dropouts" such as stellar mass, stellar age, dust reddening, and star-formation rate; and (4) to assess how representative the GOODS survey is at this epoch. Our simulations predict that high-redshift galaxies show strong correlations in star formation rate versus stellar mass, and weaker correlations versus environment and age, such that GOODS galaxies are predicted to be the most massive, most rapidly star-forming galaxies at that epoch, living preferentially in dense regions. The simulated rest-frame UV luminosity function (LF) and integrated luminosity density are in broad agreement with observations at z~4. The predicted rest-frame optical (observed 3.6 micron) LF is similar to the rest-frame UV LF, shifted roughly one magnitude brighter. We predict that GOODS detects less than 50% of the total stellar mass density formed in galaxies more massive than 10^8.7 M_sun by z=4, mainly because of brightness limits in the HST/ACS bands. The most rapidly star forming galaxies in our simulations have rates exceeding 1000 M_sun yr^-1, similar to observed sub-mm galaxies. The star formation rates of these galaxies show at most a mild excess (2--3x) over the rates that would be expected for their stellar mass. Whether these bright galaxies would be observable as LBGs depends on the uncertain effects of dust reddening., Comment: 27 pages, 15 figures, submitted to ApJ

Fallback disks are expected to form around new-born neutron stars following a supernova explosion. In almost all cases, the disk will pass through a propeller stage. If the neutron star is spinning rapidly (initial period $\sim 10$ ms) and has an ordinary magnetic moment ($\sim 10^{30}$ G cm$^3$), the rotational power transferred to the disk by the magnetic field of the neutron star will exceed the Eddington limit by many orders of magnitude, and the disk will be rapidly disrupted. Fallback disks can thus survive only around slow-born neutron stars and around black holes, assuming the latter do not torque their surrounding disks as strongly as do neutron stars. This might explain the apparent rarity of fallback disks around young compact objects., Comment: Submitted to Astrophysical Journal Letters

We use the LUQAS sample of 27 high resolution, high signal-to-noise QSO absorption spectra (Kim et al. 2004) and the results from Croft et al. (2002) together with a suite of high-resolution hydro-dynamical simulations run with the GADGET-II code, to infer the linear dark matter power spectrum on scales of \~ 0.3-30 Mpc at z=2.125 and at z=2.72., Comment: 4 pages, 2 figures; Talk given by M. Viel at the XXXIXth Rencontres de Moriond on "Exploring the Universe", La Thuile, Italy, March 28 - April 4, 2004

We describe hydrodynamical simulations of galaxy formation in a Lambda cold dark matter (CDM) cosmology performed using a subresolution model for star formation and feedback in a multiphase interstellar medium (ISM). In particular, we demonstrate the formation of a well-resolved disk galaxy. The surface brightness profile of the galaxy is exponential, with a B-band central surface brightness of 21.0 mag arcsec^-2 and a scale-length of R_d = 2.0 h^-1 kpc. We find no evidence for a significant bulge component. The simulated galaxy falls within the I-band Tully-Fisher relation, with an absolute magnitude of I = -21.2 and a peak stellar rotation velocity of V_rot=121.3 km s^-1. While the total specific angular momentum of the stars in the galaxy agrees with observations, the angular momentum in the inner regions appears to be low by a factor of ~2. The star formation rate of the galaxy peaks at ~7 M_sun yr^-1 between redshifts z=2-4, with the mean stellar age decreasing from \~10 Gyrs in the outer regions of the disk to ~7.5 Gyrs in the center, indicating that the disk did not simply form inside-out. The stars exhibit a metallicity gradient from 0.7 Z_sun at the edge of the disk to 1.3 Z_sun in the center. Using a suite of idealized galaxy formation simulations with different models for the ISM, we show that the effective pressure support provided by star formation and feedback in our multiphase model is instrumental in allowing the formation of large, stable disk galaxies. If ISM gas is instead modeled with an isothermal equation of state, or if star formation is suppressed entirely, growing gaseous disks quickly violate the Toomre stability criterion and undergo catastrophic fragmentation., Comment: 14 pages, 12 figures, LaTex (emulateapj.cls), submitted to ApJ, high resolution images available at http://www-cfa.harvard.edu/~brobertson/papers/galaxy/

We study the soft X-ray emission (0.1 - 1 keV) from the Warm-Hot Intergalactic Medium (WHIM) in a hydrodynamic simulation of a Cold Dark Matter universe. Our main goal is to investigate how such emission can be explored with a combination of imaging and spectroscopy, and to motivate future X-ray missions. We first present high resolution images of the X-ray emission in several energy bands, in which emission from different ion species dominates. We pick three different areas to study the high resolution spectra of X-rays from the warm-hot IGM: (a) a galaxy group; (b) a filament and (c) an underluminous region. By taking into account the background X-ray emission from AGNs and foreground emission from the Galaxy, we compute composite X-ray spectra of the selected regions. We briefly investigate angular clustering of the soft-X-ray emission, finding a strong signal. Most interestingly, the combination of high spectral resolution and angular information allows us to map the emission from the WHIM in 3 dimensions. We cross-correlate the positions of galaxies in the simulation with this redshift map of emission and detect the presence of six different ion species (Ne IX, Fe XVII, O VII, O VIII, N VII, C VI) in the large-scale structure traced by the galaxies. Finally we show how such emission can be detected and studied with future X-ray satellites, with particular attention to a proposed mission, the Missing Baryon Explorer, or MBE. We present simulated observations of the WHIM gas with MBE., Comment: 17 pages, 17 figures. Revised version, accepted by ApJ

(abridged) An outburst of more than 80 individual bursts, similar to those seen from Soft Gamma Repeaters (SGRs), was detected from the Anomalous X-ray Pulsar (AXP) 1E 2259+586 in 2002 June. Coincident with this burst activity were gross changes in the pulsed flux, persistent flux, energy spectrum, pulse profile and spin down of the underlying X-ray source. We present RXTE and XMM-Newton observations of 1E 2259+586 that show the evolution of the aforementioned source parameters during and following this episode. Specifically, we observe an X-ray flux increase by more than an order of magnitude having two distinct components. The first component is linked to the burst activity and decays within ~2 days during which the energy spectrum is considerably harder than during the quiescent state of the source. The second component decays over the year following the glitch according to a power law in time with an exponent -0.22 +/- 0.01. The pulsed fraction decreased initially to ~15% RMS, and the pulse profile changed significantly during the outburst. A glitch was observed in 1E 2259+586 that preceded the observed burst activity. A fraction of the glitch (~19%) recovered, although the recovery was not purely exponential. An exponential rise of ~20% of the frequency jump with a time scale of ~14 days results in a significantly better fit to the data, however, contamination from a systematic drift in the phase of the pulse profile cannot be excluded. The long-term post-glitch spin-down rate decreased in magnitude relative to the pre-glitch value. A comparison with SGR outburst properties, a physical interpretation of the results, and implications on the number of magnetar candidates in our Galaxy are discussed., Comment: Accepted for publication in ApJ. 51 pages and 13 figures. Figure 1 is modified significantly. Several other minor changes made to the text

Chandra images present evidence for a non-uniform spatial distribution of discrete X-ray sources in the elliptical galaxy NGC 4261. This non-uniform distribution is inconsistent with the optical morphology of NGC 4261 at greater than a 99.9% confidence level. Similar evidence is seen in one more elliptical galaxy (NGC 4697; 98% confidence level) out of five cases we investigated. NGC 4261 and NGC 4697 have old stellar populations (9-15 Gyrs) and fine structure parameters of 1 and 0 respectively, suggesting no recent merging activity. On the basis of simulations of galaxy interactions, we propose that the X-ray sources responsible for the non-uniform distribution are associated with young stellar populations, related to the rejuvenating fall-back of material in tidal tails onto a relaxed merger remnant, or shock induced star-formation along the tidal tails., Comment: Accepted for publication in Ap.J. Letters

We present a determination of the bispectrum of the flux in the Lyman-alpha forest of QSO absorption spectra obtained from LUQAS which consists of spectra observed with the high resolution Echelle spectrograph UVES. Typical errors on the observed bispectrum as obtained from a jack-knife estimator are ~ 50%. For wavenumbers in the range 0.03 s/km 2., Comment: 6 pages, 3 figures, 1 table. A table can be found at http://www.ast.cam.ac.uk/~rtnigm/luqas.htm To be published in MNRAS Letters

We present a pulsar-atmosphere (PA) model for modulated thermal X-ray emission from cooling magnetized neutron stars. The model synthesizes the spectral properties of detailed stellar atmosphere calculations with the non-uniform surface properties anticipated for isolated, aging radio pulsars, and general relativistic effects on photon trajectories. We analyze the archival Chandra observations of the middle-aged radio pulsar PSR 0656+14 with the PA model and find it is an excellent representation of the phase-averaged X-ray spectrum for a broad range of polar effective temperature T_p and column density N_H ~ 10^{20} cm^{-2}. The spectral fits favor a sub-solar neutron star mass M < 1.0 M_{\sun}, a large radius R=15-16 km for a source distance d ~ 190 pc, while retaining consistency with theoretical neutron star equations-of-state. The modulated spectrum constrains the angular displacement of the pulsar spin axis to 30 +/- 2 \arcdeg with respect to the line of sight., Comment: 21 pages, 7 figures

We describe techniques for comparing spectra extracted from cosmological simulations and observational data, using the same methodology to link Lyman-alpha properties derived from the simulations with properties derived from observational data. The eventual goal is to measure the coherence or clustering properties of Lyman-alpha absorbers using observations of quasar pairs and groups. We quantify the systematic underestimate in opacity that is inherent in the continuum fitting process of observed spectra over a range of resolution and signal-to-noise ratio. We present an automated process for detecting and selecting absorption features over the range of resolution and signal-to-noise of typical observational data on the Lyman-alpha "forest". Using these techniques, we detect coherence over transverse scales out to 500 h^{-1}_{50} kpc in spectra extracted from a cosmological simulation at z = 2., Comment: 52 pages, includes 14 figures, to appear in ApJ v566 Feb 2002

This paper discusses a semi-numerical method of investigating the enrichment of the intergalactic medium by galactic winds. We find that most galaxies at z >~ 3 should be driving winds, and that (if these winds are similar to those at low-z) these winds should escape to large distances. Our calculations -- which permit exploration of a large region of model parameter space -- indicate that the wind velocity, the mass of the wind-driving galaxies, the fraction of ambient material entrained, and the available time (between wind launch and the observed redshift) all affect wind propagation significantly; other physical effects can be important but are sub-dominant. We find that under reasonable assumptions, the enrichment by 3 <~ z <~ 6 galaxies could account for the quantity of metals seen in the Ly-alpha forest, though it is presently unclear whether this enrichment is compatible with the intergalactic medium's detailed metal distribution or relative quiescence., Comment: 7 pages, 4 ps figures; to appear in "Chemical Enrichment of Intracluster and Intergalactic Medium", Vulcano, Italy, 14-18 May 2001, ASP Conference Series

WFPC2 images and STIS spectroscopic observations are presented of the double nucleus in the merger system NGC 6240. We find that: (a) the kinematics of the ionized gas is similar to that of the molecular gas, despite a different morphology; (b) the gaseous and stellar kinematics are quite different, suggesting an early merger stage; (c) neither the gaseous nor the stellar kinematics show an obvious sign of the supermassive black hole believed to be responsible for the X-ray emission of NGC 6240; and (d) the steep off-nuclear velocity gradient is not due to a 10E11 solar mass black hole, in contrast to earlier suggestions., Comment: To appear in: "The central kpc of starbursts and AGN: the La Palma connection", Eds. J.H. Knapen, J.E. Beckman, I. Shlosman and T.J. Mahoney, ASP conf. series, in press (2001)

We present the analysis of a suite of simulations of a Virgo mass galaxy cluster. Undertaken within the framework of standard cold dark matter cosmology, these simulations were performed at differing resolutions and with increasingly complex physical processes, with the goal of identifying the effects of each on the evolution of the cluster. We focus on the cluster at the present epoch and examine properties including the radial distributions of density, temperature, entropy and velocity. We also map `observable' projected properties such as the surface mass density, X-ray surface brightness and SZ signature. We identify significant differences between the simulations, which highlights the need for caution when comparing numerical simulations to observations of galaxy clusters. While resolution affects the inner density profile in dark matter simulations, the addition of a gaseous component, especially one that cools and forms stars, affects the entire cluster. We conclude that both resolution and included physical processes play an important role in simulating the formation and evolution of galaxy clusters. Therefore, physical inferences drawn from simulations that do not include a gaseous component that can cool and form stars present a poor representation of reality. (Abridged), Comment: Accepted for publication in the Astrophysical Journal. Several changes from previous version, including new material

We propose a new method of analysis for the \lya forest, namely to measure the 1-point and 2-point joint probability distribution of the transmitted flux. The results for a sample of seven observed quasars and from two simulations of structure formation are shown and compared. Statistically significant differences in the 2-point function between the results of the numerical simulations and the observations are easily found. The analysis we suggest is very simple to apply to observed data sets, and we discuss its superiority over the traditional Voigt-profile fitting algorithms for accurate comparison to the predictions of theoretical models., Comment: 7 pages, 3 postscript figures, uses conf_iap.sty. To appear in Proceedings of 13th IAP Colloquium: Structure and Evolution of the IGM from QSO Absorption Line Systems, eds. P. Petitjean, S. Charlot

We study the stability of stellar dynamical equilibrium models for M32. Kinematic observations show that M32 has a central black hole of 3x10^6 solar masses, and a phase-space distribution function that is close to the `two-integral' form f=f(E,L_z). M32 is also rapidly rotating; 85-90% of the stars have the same sense of rotation around the symmetry axis. Previous work has shown that flattened, rapidly rotating two-integral models can be bar-unstable. We have performed N-body simulations to test whether this is the case for M32. Particle realizations with N=512,000 were studied for two representative inclinations, i=90 (edge-on) and i=55, corresponding to intrinsic axial ratios of q=0.73 and q=0.55, respectively. The time evolution of the models was calculated with a `self-consistent field' code on a Cray T3D parallel supercomputer. We find both models to be dynamically stable. This implies that they provide a physically meaningful description of M32, and that the inclination of M32 (and hence its intrinsic flattening) cannot be strongly constrained through stability arguments. Previous work on the stability of f(E,L_z) models has shown that the bar-mode is the only possibly unstable mode for systems rounder than q=0.3, and that the likelihood for this mode to be unstable increases with flattening and rotation rate. The f(E,L_z) models studied for M32 are stable, and M32 has a higher rotation rate than nearly all other elliptical galaxies. This suggests that f(E,L_z) models constructed to fit data for real elliptical galaxies will generally be stable for q>0.55, and possibly for flatter systems as well., Comment: 23 pages, LaTeX, with 8 PostScript figures; submitted to the Astrophysical Journal; mpeg animation of the stellar motions in M32 available at http://oposite.stsci.edu/pubinfo/Anim.html

We have measured the distribution function of the flux decrement caused by Lyman alpha forest absorption in a new sample of high resolution QSO spectra. The observations are compared to the results from two simulations of the Lya forest: an Eulerian Lambda-CDM model, and an SPH standard CDM model. Good agreement between the shapes of simulated and observed distributions is achieved by globally scaling the optical depth to match the mean flux decrements. This procedure amounts to a measurement of the parameter Omega_b^2 h^3 / Gamma (where Omega_b is the baryonic matter density and Gamma is the HI ionization rate). Estimating a lower limit Gamma > 7 10^{-13} s^{-1} from the abundance of known QSOs, we derive a lower limit to the baryon density, Omega_b h^2>0.021(0.017) for the Lambda-CDM (SCDM) model. In both cases the large values are inconsistent with some recent D/H determinations (Rugers & Hogan 1996a,b), favoring a low deuterium abundance as reported by Tytler, Fan & Burles (1996). Adopting a fixed Omega_b, we can determine the evolution of the ionizing radiation field. Our models predict the intensity to be approximately constant with redshift, consistent with the assumption that the ionizing background is produced by known quasars for z < 3. However, additional sources of ionizing photons are required at higher redshift., Comment: 35 pages latex (uses aaspp4 and psfig.sty), 10 postscript figures; submitted to ApJ. The complete paper can also be retrieved at http://astro.caltech.edu/~mr/preprints/fluxdec.ps.gz

Perturbation theory makes it possible to calculate the probability distribution function (PDF) of the large scale density field in the small variance limit. For top hat smoothing and scale-free Gaussian initial fluctuations, the result depends only on the linear variance, sigma_linear, and its logarithmic derivative with respect to the filtering scale -(n_linear+3)=dlog sigma_linear^2/dlog L (Bernardeau 1994). In this paper, we measure the PDF and its low-order moments in scale-free simulations evolved well into the nonlinear regime and compare the results with the above predictions, assuming that the spectral index and the variance are adjustable parameters, n_eff and sigma_eff=sigma, where sigma is the true, nonlinear variance. With these additional degrees of freedom, results from perturbation theory provide a good fit of the PDFs, even in the highly nonlinear regime. The value of n_eff is of course equal to n_linear when sigma << 1, and it decreases with increasing sigma. A nearly flat plateau is reached when sigma >> 1. In this regime, the difference between n_eff and n_linear increases when n_linear decreases. For initial power-spectra with n_linear=-2,-1,0,+1, we find n_eff ~ -9,-3,-1,-0.5 when sigma^2 ~ 100., Comment: 13 pages, 6 figures, Latex (MN format), submitted to MNRAS

We measure the scaling properties of the probability distribution of the smoothed density field in $N$-body simulations of expanding universes with scale-free initial power-spectra, with particular attention to the predictions of the stable clustering hypothesis. We concentrate our analysis on the ratios $S_Q(\ell)\equiv {\bar \xi}_Q/{\bar \xi}_2^{Q-1}$, $Q \leq 5$, where ${\bar \xi}_Q$ is the averaged $Q$-body correlation function over a cell of radius $\ell$. The behavior of the higher order correlations is studied through that of the void probability distribution function. As functions of ${\bar \xi}_2$, the quantities $S_Q$, $3 \leq Q \leq 5$, exhibit two plateaus separated by a smooth transition around ${\bar \xi}_2 \sim 1$. In the weakly nonlinear regime, ${\bar \xi}_2 \la 1$, the results are in reasonable agreement with the predictions of perturbation theory. In the nonlinear regime, ${\bar \xi}_2 > 1$, the function $S_Q({\bar \xi}_2)$ is larger than in the weakly nonlinear regime, and increasingly so with $-n$. It is well-fitted by the expression $S_Q= ({\bar \xi}_2/100)^{0.045(Q-2)}\ {\widetilde S}_Q$ for all $n$. This weak dependence on scale proves {\em a small, but significant departure from the stable clustering predictions} at least for $n=0$ and $n=+1$. The analysis of $P_0$ confirms that the expected scale-invariance of the functions $S_Q$ is not exactly attained in the part of the nonlinear regime we probe, except possibly for $n=-2$ and marginally for $n=-1$. In these two cases, our measurements are not accurate enough to be discriminant., Comment: 31 pages, postscript file, figure 1 missing. Postscript file including figure 1 available at ftp://ftp-astro-theory.fnal.gov:/pub/Publications/Pub-95-256-A