7,152 research outputs found
A scheme for radiation pressure and photon diffusion with the M1 closure in RAMSES-RT
We describe and test an updated version of radiation-hydrodynamics (RHD) in
the RAMSES code, that includes three new features: i) radiation pressure on
gas, ii) accurate treatment of radiation diffusion in an unresolved optically
thick medium, and iii) relativistic corrections that account for Doppler
effects and work done by the radiation to first order in v/c. We validate the
implementation in a series of tests, which include a morphological assessment
of the M1 closure for the Eddington tensor in an astronomically relevant
setting, dust absorption in a optically semi-thick medium, direct pressure on
gas from ionising radiation, convergence of our radiation diffusion scheme
towards resolved optical depths, correct diffusion of a radiation flash and a
constant luminosity radiation, and finally, an experiment from Davis et al. of
the competition between gravity and radiation pressure in a dusty atmosphere,
and the formation of radiative Rayleigh-Taylor instabilities. With the new
features, RAMSES-RT can be used for state-of-the-art simulations of radiation
feedback from first principles, on galactic and cosmological scales, including
not only direct radiation pressure from ionising photons, but also indirect
pressure via dust from multi-scattered IR photons reprocessed from
higher-energy radiation, both in the optically thin and thick limits.Comment: 25 pages, 13 figures, accepted for publication in MNRAS. Revised to
match published versio
Globular Cluster Formation in the Virgo Cluster
Metal poor globular clusters (MPGCs) are a unique probe of the early
universe, in particular the reionization era. Systems of globular clusters in
galaxy clusters are particularly interesting as it is in the progenitors of
galaxy clusters that the earliest reionizing sources first formed. Although the
exact physical origin of globular clusters is still debated, it is generally
admitted that globular clusters form in early, rare dark matter peaks (Moore et
al. 2006; Boley et al. 2009). We provide a fully numerical analysis of the
Virgo cluster globular cluster system by identifying the present day globular
cluster system with exactly such early, rare dark matter peaks. A popular
hypothesis is that that the observed truncation of blue metal poor globular
cluster formation is due to reionization (Spitler et al. 2012; Boley et al.
2009; Brodie & Strader 2006); adopting this view, constraining the formation
epoch of MPGCs provides a complementary constraint on the epoch of
reionization. By analyzing both the line of sight velocity dispersion and the
surface density distribution of the present day distribution we are able to
constrain the redshift and mass of the dark matter peaks. We find and quantify
a dependence on the chosen line of sight of these quantities, whose strength
varies with redshift, and coupled with star formation efficiency arguments find
a best fitting formation mass and redshift of and . We predict intracluster MPGCs in
the Virgo cluster. Our results confirm the techniques pioneered by Moore et al.
(2006) when applied to the the Virgo cluster and extend and refine the analytic
results of Spitler et al. (2012) numerically.Comment: 13 Pages, 13 Figures, submitted to MNRA
Infall near clusters of galaxies: comparing gas and dark matter velocity profiles
We consider the dynamics in and near galaxy clusters. Gas, dark matter and
galaxies are presently falling into the clusters between approximately 1 and 5
virial radii. At very large distances, beyond 10 virial radii, all matter is
following the Hubble flow, and inside the virial radius the matter particles
have on average zero radial velocity. The cosmological parameters are imprinted
on the infall profile of the gas, however, no method exists, which allows a
measurement of it. We consider the results of two cosmological simulations
(using the numerical codes RAMSES and Gadget) and find that the gas and dark
matter radial velocities are very similar. We derive the relevant dynamical
equations, in particular the generalized hydrostatic equilibrium equation,
including both the expansion of the Universe and the cosmological background.
This generalized gas equation is the main new contribution of this paper. We
combine these generalized equations with the results of the numerical
simulations to estimate the contribution to the measured cluster masses from
the radial velocity: inside the virial radius it is negligible, and inside two
virial radii the effect is below 40%, in agreement the earlier analyses for DM.
We point out how the infall velocity in principle may be observable, by
measuring the gas properties to distance of about two virial radii, however,
this is practically not possible today.Comment: 7 pages, 3 figures, to appear in MNRA
Modelling CO emission from hydrodynamic simulations of nearby spirals, starbursting mergers, and high-redshift galaxies
We model the intensity of emission lines from the CO molecule, based on
hydrodynamic simulations of spirals, mergers, and high-redshift galaxies with
very high resolutions (3pc and 10^3 Msun) and detailed models for the
phase-space structure of the interstellar gas including shock heating, stellar
feedback processes and galactic winds. The simulations are analyzed with a
Large Velocity Gradient (LVG) model to compute the local emission in various
molecular lines in each resolution element, radiation transfer and opacity
effects, and the intensity emerging from galaxies, to generate synthetic
spectra for various transitions of the CO molecule. This model reproduces the
known properties of CO spectra and CO-to-H2 conversion factors in nearby
spirals and starbursting major mergers. The high excitation of CO lines in
mergers is dominated by an excess of high-density gas, and the high turbulent
velocities and compression that create this dense gas excess result in broad
linewidths and low CO intensity-to-H2 mass ratios. When applied to
high-redshift gas-rich disks galaxies, the same model predicts that their
CO-to-H2 conversion factor is almost as high as in nearby spirals, and much
higher than in starbursting mergers. High-redshift disk galaxies contain giant
star-forming clumps that host a high-excitation component associated to gas
warmed by the spatially-concentrated stellar feedback sources, although CO(1-0)
to CO(3-2) emission is overall dominated by low-excitation gas around the
densest clumps. These results overall highlight a strong dependence of CO
excitation and the CO-to-H2 conversion factor on galaxy type, even at similar
star formation rates or densities. The underlying processes are driven by the
interstellar medium structure and turbulence and its response to stellar
feedback, which depend on global galaxy structure and in turn impact the CO
emission properties.Comment: A&A in pres
A scheme for radiation pressure and photon diffusion with the M1 closure in ramses-rt
We describe and test an updated version of radiation-hydrodynamics in the ramses code, that includes three new features: (i) radiation pressure on gas, (ii) accurate treatment of radiation diffusion in an unresolved optically thick medium, and (iii) relativistic corrections that account for Doppler effects and work done by the radiation to first order in v/c. We validate the implementation in a series of tests, which include a morphological assessment of the M1 closure for the Eddington tensor in an astronomically relevant setting, dust absorption in an optically semithick medium, direct pressure on gas from ionizing radiation, convergence of our radiation diffusion scheme towards resolved optical depths, correct diffusion of a radiation flash and a constant luminosity radiation, and finally, an experiment from Davis etal. of the competition between gravity and radiation pressure in a dusty atmosphere, and the formation of radiative Rayleigh-Taylor instabilities. With the new features, ramses-rt can be used for state-of-the-art simulations of radiation feedback from first principles, on galactic and cosmological scales, including not only direct radiation pressure from ionizing photons, but also indirect pressure via dust from multiscattered IR photons reprocessed from higher-energy radiation, both in the optically thin and thick limit
Collapse, outflows and fragmentation of massive, turbulent and magnetized prestellar barotropic cores
Stars and more particularly massive stars, have a drastic impact on galaxy
evolution. Yet the conditions in which they form and collapse are still not
fully understood. In particular, the influence of the magnetic field on the
collapse of massive clumps is relatively unexplored, it is thus of great
relevance in the context of the formation of massive stars to investigate its
impact. We perform high resolution, MHD simulations of the collapse of hundred
solar masses, turbulent and magnetized clouds, using the adaptive mesh
refinement code RAMSES. We compute various quantities such as mass
distribution, magnetic field and angular momentum within the collapsing core
and study the episodic outflows and the fragmentation that occurs during the
collapse. The magnetic field has a drastic impact on the cloud evolution. We
find that magnetic braking is able to substantially reduce the angular momentum
in the inner part of the collapsing cloud. Fast and episodic outflows are being
launched with typical velocities of the order of 3-5 km s although the
highest velocities can be as high as 30-40 km s. The fragmentation in
several objects, is reduced in substantially magnetized clouds with respect to
hydrodynamical ones by a factor of the order of 1.5-2. We conclude that
magnetic fields have a significant impact on the evolution of massive clumps.
In combination with radiation, magnetic fields largely determine the outcome of
massive core collapse. We stress that numerical convergence of MHD collapse is
a challenging issue. In particular, numerical diffusion appears to be important
at high density therefore possibly leading to an over-estimation of the number
of fragments.Comment: accepted for publication in A&
- …