Adaptive Mesh Refinement Code RAMSES Research Articles

Discs are born eccentric

Context. Recent observations have begun probing the early phases of disc formation, but little data yet exists on disc structure and morphology of Class 0 objects. Using simulations, we are able to lay out predictions of disc morphologies expected in future surveys of young discs. Based on detailed simulations of ab initio star formation by core collapse, we predict that early discs must be eccentric. Aims. In this Letter, we study the morphology and, in particular, the eccentricity of discs formed in non-ideal magnetohydrodynamic (MHD) collapse simulations. We attempt to show that discs formed by cloud collapse are likely to be eccentric. Methods. We ran non-ideal MHD collapse simulations in the adaptive mesh refinement code RAMSES with radiative transfer. We used state-of-the-art analysis methods to measure the disc eccentricity. Results. We find that despite no asymmetry in the initial conditions, the discs formed are eccentric, with eccentricities on the order of 0.1. Conclusions. These results may have important implications for protoplanetary disc dynamics and planet formation. The presence of eccentricity in young discs that is not seen at later stages of disc evolution is in tension with current viscous eccentricity damping models. This implies that there may be an as-yet undiscovered circularisation mechanism in circumstellar discs.

The contribution of binary star formation via core fragmentation on protostellar multiplicity

Context. Observations of young multiple star systems find a bimodal distribution in companion frequency and separation. The origin of these peaks has often been attributed to binary formation via core and disc fragmentation. However, theory and simulations suggest that young stellar systems that form via core fragmentation undergo significant orbital evolution. Aims. We investigate the influence of the environment on the formation and orbital evolution of multiple star systems, and how core fragmentation contributes to the formation of close (20 − 100 AU) binaries. We use multiple simulations of star formation in giant molecular clouds and compare them to the multiplicity statistics of the Perseus star-forming region. Methods. Simulations were run with the adaptive mesh refinement code RAMSES with sufficient resolution to resolve core fragmentation beyond 400 AU and dynamical evolution down to 16.6 AU, but without the possibility of resolving disc fragmentation. The evolution of the resulting stellar systems was followed over millions of years. Results. We find that star formation in lower gas density environments is more clustered; however, despite this, the fractions of systems that form via dynamical capture and core fragmentation are broadly consistent at ∼40% and ∼60%, respectively. In all gas density environments, we find that the typical scale at which systems form via core fragmentation is 103 − 3.5 AU. After formation, we find that systems that form via core fragmentation have slightly lower inspiral rates (∼10−1.68 AU yr−1 measured over the first 10 000 yr) compared to dynamical capture (∼10−1.32 AU yr−1). We then compared the simulation with the conditions most similar to the Perseus star-forming region to determine whether the observed bimodal distribution can be replicated. We find that it can be replicated, but it is sensitive to the evolutionary state of the simulation. Conclusions. Our results indicate that a significant number of low-mass close binaries with separations from 20 − 100 AU can be produced via core fragmentation or dynamical capture due to efficient inspiral, without the need for a further contribution from disc fragmentation.

Gasflows in Barred Galaxies with Big Orbital Loops—A Comparative Study of Two Hydrocodes

We study the flow of gas in a barred-galaxy model, in which a considerable part of the underlying stable periodic orbits have loops where, close to the ends of the bar, several orbital families coexist and chaos dominates. Such conditions are typically encountered in a zone between the 4:1 resonance and corotation. The purpose of our study is to understand the gaseous flow in the aforementioned environment and trace the morphology of the shocks that form. We use two conceptually different hydrodynamic schemes for our calculations, namely, the mesh-free Lagrangian SPH method and the adaptive mesh refinement code RAMSES. This allows us to compare responses by means of the two algorithms. We find that the big loops of the orbits, mainly belonging to the x1 stable periodic orbits, do not help the shock loci to approach corotation. They deviate away from the regions occupied by the loops, bypass them and form extensions at an angle with the straight-line shocks. Roughly at the distance from the center at which we start to observe the big loops, we find characteristic “tails” of dense gas streaming towards the straight-line shocks. The two codes give complementary information for understanding the hydrodynamics of the models.

ACACIA: a new method to produce on-the-fly merger trees in the ramses code

ABSTRACT The implementation of ACACIA, a new algorithm to generate dark matter halo merger trees with the Adaptive Mesh Refinement code RAMSES, is presented. The algorithm is fully parallel and based on the Message Passing Interface. As opposed to most available merger tree tools, it works on the fly during the course of the N-body simulation. It can track dark matter substructures individually using the index of the most bound particle in the clump. Once a halo (or a sub-halo) merges into another one, the algorithm still tracks it through the last identified most bound particle in the clump, allowing to check at later snapshots whether the merging event was definitive, or whether it was only temporary, with the clump only traversing another one. The same technique can be used to track orphan galaxies that are not assigned to a parent clump anymore because the clump dissolved due to numerical overmerging. We study in detail the impact of various parameters on the resulting halo catalogues and corresponding merger histories. We then compare the performance of our method using standard validation diagnostics, demonstrating that we reach a quality similar to the best available and commonly used merger tree tools. As a proof of concept, we use our merger tree algorithm together with a parametrized stellar-mass-to-halo-mass relation and generate a mock galaxy catalogue that shows good agreement with observational data.

Protoplanetary Disk Birth in Massive Star-forming Clumps: The Essential Role of the Magnetic Field

Abstract Protoplanetary disks form through angular momentum conservation in collapsing dense cores. In this work, we perform the first simulations with a maximal resolution down to the astronomical unit (au) of protoplanetary disk formation, through the collapse of 1000 M ⊙ clumps, treating self-consistently both non-ideal magnetohydrodynamics with ambipolar diffusion as well as radiative transfer in the flux-limited diffusion approximation including stellar feedback. Using the adaptive mesh-refinement code RAMSES, we investigate the influence of the magnetic field on the disks properties with three models. We show that, without magnetic fields, a population dominated by large disks is formed that is not consistent with Class 0 disk properties as estimated from observations. The inclusion of magnetic field leads, through magnetic braking, to a very different evolution. When it is included, small <50 au disks represent about half the population. In addition, about 70% of the stars have no disk in this case, which suggests that our resolution is still insufficient to preserve the smaller disks. With ambipolar diffusion, the proportion of small disks is also prominent and we report a flat mass distribution around 0.01–0.1M ⊙ and a typical disk-to-star mass ratios of ∼10−2–10−1. This work shows that the magnetic field and its evolution plays a prominent role in setting the initial properties of disk populations.

Cosmic ray feedback from supernovae in dwarf galaxies

The regulation of the baryonic content in dwarf galaxies is a long-standing problem. Supernovae (SNe) are supposed to play a key role in forming large-scale galactic winds by removing important amounts of gas from galaxies. SNe are efficient accelerators of non-thermal particles, so-called cosmic rays (CRs), which can substantially modify the dynamics of the gas and conditions to form large-scale galactic winds. We investigate how CR injection by SNe impacts the star formation and the formation of large-scale winds in dwarf galaxies, and whether it can produce galaxy star-formation rates (SFR) and wind properties closer to observations. We ran CR magneto-hydrodynamical simulations of dwarf galaxies at high resolution (9 pc) with the adaptive mesh refinement codeRAMSES. Those disc galaxies are embedded in isolated halos of mass of 1010and 1011 M⊙, and CRs are injected by SNe. We included CR isotropic and anisotropic diffusion with various diffusion coefficients, CR radiative losses, and CR streaming. The injection of CR energy into the interstellar medium smooths out the highest gas densities, which reduces the SFR by a factor of 2–3. Mass outflow rates are significantly greater with CR diffusion, by 2 orders of magnitudes for the higher diffusion coefficients. Without diffusion and streaming, CRs are inefficient at generating winds. CR streaming alone allows for the formation of winds but which are too weak to match observations. The formation of galactic winds strongly depends on the diffusion coefficient: for low coefficients, CR energy stays confined in high density regions where CR energy losses are highest, and higher coefficients, which allow for a more efficient leaking of CRs out of dense gas, produce stronger winds. CR diffusion leads to colder and denser winds than without CRs, and brings outflow rates and mass loading factors much closer to observations.

The impact of AGN feedback on the 1D power spectra from the Ly α forest using the Horizon-AGN suite of simulations

ABSTRACT The Lyman-α forest is a powerful probe for cosmology, but it is also strongly impacted by galaxy evolution and baryonic processes such as active galactic nucleus (AGN) feedback, which can redistribute mass and energy on large scales. We constrain the signatures of AGN feedback on the 1D power spectrum of the Lyman-α forest using a series of eight hydro-cosmological simulations performed with the adaptive mesh refinement code ramses. This series starts from the Horizon-AGN simulation and varies the subgrid parameters for AGN feeding, feedback, and stochasticity. These simulations cover the whole plausible range of feedback and feeding parameters according to the resulting galaxy properties. AGNs globally suppress the Lyman-α power at all scales. On large scales, the energy injection and ionization dominate over the supply of gas mass from AGN-driven galactic winds, thus suppressing power. On small scales, faster cooling of denser gas mitigates the suppression. This effect increases with decreasing redshift. We provide lower and upper limits of this signature at nine redshifts between z = 4.25 and 2.0, making it possible to account for it at post-processing stage in future work given that running simulations without AGN feedback can save considerable amounts of computing resources. Ignoring AGN feedback in cosmological inference analyses leads to strong biases with 2 per cent shift on σ8 and 1 per cent shift on ns, which represents twice the standards deviation of the current constraints on ns.

What determines the formation and characteristics of protoplanetary discs?

Context. Planets form in protoplanetary discs. Their masses, distribution, and orbits sensitively depend on the structure of the protoplanetary discs. However, what sets the initial structure of the discs in terms of mass, radius and accretion rate is still unknown. Aims. It is therefore of great importance to understand exactly how protoplanetary discs form and what determines their physical properties. We aim to quantify the role of the initial dense core magnetisation, rotation, turbulence, and misalignment between rotation and magnetic field axis as well as the role of the accretion scheme onto the central object. Methods. We performed non-ideal magnetohydrodynamics numerical simulations using the adaptive mesh refinement code Ramses of a collapsing, one solar mass molecular core to study the disc formation and early, up to 100 kyr, evolution. We paid particular attention to the impact of numerical resolution and accretion scheme. Results. We found that the mass of the central object is almost independent of the numerical parameters such as the resolution and the accretion scheme onto the sink particle. The disc mass and to a lower extent its size, however heavily depend on the accretion scheme, which we found is itself resolution dependent. This implies that the accretion onto the star and through the disc are largely decoupled. For a relatively large domain of initial conditions (except at low magnetisation), we found that the properties of the disc do not change too significantly. In particular both the level of initial rotation and turbulence do not influence the disc properties provide the core is sufficiently magnetised. After a short relaxation phase, the disc settles in a stationary state. It then slowly grows in size but not in mass. The disc itself is weakly magnetised but its immediate surrounding on the contrary is highly magnetised. Conclusions. Our results show that the disc properties directly depend on the inner boundary condition, i.e. the accretion scheme onto the central object. This suggests that the disc mass is eventually controlled by a small-scale accretion process, possibly the star-disc interaction. Because of ambipolar diffusion and its significant resistivity, the disc diversity remains limited and except for low magnetisation, their properties are weakly sensitive to initial conditions such as rotation and turbulence.

Global simulations of galactic discs: violent feedback from clustered supernovae during bursts of star formation

ABSTRACT A suite of idealized, global, gravitationally unstable, star-forming galactic disc simulations with 2 pc spatial resolution, performed with the adaptive mesh refinement code ramses, is used in this paper to predict the emergent effects of supernova feedback. The simulations include a simplified prescription for the formation of single stellar populations of mass $\sim 100 \, {\rm M}_{\odot }$, radiative cooling, photoelectric heating, an external gravitational potential for a dark matter halo and an old stellar disc, self-gravity, and a novel implementation of supernova feedback. The results of these simulations show that gravitationally unstable discs can generate violent supersonic winds with mass-loading factors η ≳ 10, followed by a galactic fountain phase. These violent winds are generated by highly clustered supernovae exploding in dense environments created by gravitational instability, and they are not produced in simulation without self-gravity. The violent winds significantly perturb the vertical structure of the disc, which is later re-established during the galactic fountain phase. Gas resettles into a quasi-steady, highly turbulent disc with volume-weighted velocity dispersion $\sigma \gt 50 \, {\rm km\, s}^{-1}$. The new configuration drives weaker galactic winds with a mass-loading factor η ≤ 0.1. The whole cycle takes place in ≤10 dynamical times. Such high time variability needs to be taken into account when interpreting observations of galactic winds from starburst and post-starburst galaxies.

Shock-accelerated cosmic rays and streaming instability in the adaptive mesh refinement code Ramses

Cosmic rays (CRs) are thought to play a dynamically important role in several key aspects of galaxy evolution, including the structure of the interstellar medium, the formation of galactic winds, and the non-thermal pressure support of halos. We introduce a numerical model solving for the CR streaming instability and acceleration of CRs at shocks with a fluid approach in the adaptive mesh refinement code RAMSES. CR streaming is solved with a diffusion approach and its anisotropic nature is naturally captured. We introduce a shock finder for the RAMSES code that automatically detects shock discontinuities in the flow. Shocks are the loci for CR injection, and their efficiency of CR acceleration is made dependent on the upstream magnetic obliquity according to the diffuse shock acceleration mechanism. We show that the shock finder accurately captures shock locations and estimates the shock Mach number for several problems. The obliquity-dependent injection of CRs in the Sedov solution leads to situations where the supernova bubble exhibits large polar caps (homogeneous background magnetic field), or a patchy structure of the CR distribution (inhomogeneous background magnetic field). Finally, we combine both accelerated CRs with streaming in a simple turbulent interstellar medium box, and show that the presence of CRs significantly modifies the structure of the gas.

Dense gas formation and destruction in a simulated Perseus-like galaxy cluster with spin-driven black hole feedback

Context.Extended filamentary Hαemission nebulae are a striking feature of nearby galaxy clusters but the formation mechanism of the filaments, and the processes which shape their morphology remain unclear.Aims.We conduct an investigation into the formation, evolution and destruction of dense gas in the centre of a simulated, Perseus-like, cluster under the influence of a spin-driven jet. The jet is powered by the supermassive black hole (SMBH) located in the cluster’s brightest cluster galaxy. We particularly study the role played by condensation of dense gas from the diffuse intracluster medium, and the impact of direct uplifting of existing dense gas by the jets, in determining the spatial distribution and kinematics of the dense gas.Methods.We present a hydrodynamical simulation of an idealised Perseus-like cluster using the adaptive mesh refinement codeRAMSES. Our simulation includes a SMBH that self-consistently tracks its spin evolution via its local accretion, and in turn drives a large-scale jet whose direction is based on the black hole’s spin evolution. The simulation also includes a live dark matter (DM) halo, a SMBH free to move in the DM potential, star formation and stellar feedback.Results.We show that the formation and destruction of dense gas is closely linked to the SMBH’s feedback cycle, and that its morphology is highly variable throughout the simulation. While extended filamentary structures readily condense from the hot intra-cluster medium, they are easily shattered into an overly clumpy distribution of gas during their interaction with the jet driven outflows. Condensation occurs predominantly onto infalling gas located 5−15 kpc from the centre during quiescent phases of the central AGN, when the local ratio of the cooling time to free fall time falls below 20, i.e. whentcool/tff < 20.Conclusions.We find evidence for both condensation and uplifting of dense gas, but caution that purely hydrodynamical simulations struggle to effectively regulate the cluster cooling cycle and produce overly clumpy distributions of dense gas morphologies, compared to observation.

Emission from the circumgalactic medium: from cosmological zoom-in simulations to multiwavelength observables

ABSTRACT We simulate the flux emitted from galaxy haloes in order to quantify the brightness of the circumgalactic medium (CGM). We use dedicated zoom-in cosmological simulations with the hydrodynamical adaptive mesh refinement code ramses, which are evolved down to z = 0 and reach a maximum spatial resolution of 380 h−1 pc and a gas mass resolution up to $1.8\times 10^{5} \, h^{-1}\, \rm {M}_{\odot }$ in the densest regions. We compute the expected emission from the gas in the CGM using cloudy emissivity models for different lines (e.g. Lyα, C iv, O vi, C vi, O viii) considering UV background fluorescence, gravitational cooling and continuum emission. In the case of Lyα, we additionally consider the scattering of continuum photons. We compare our predictions to current observations and find them to be in good agreement at any redshift after adjusting the Lyα escape fraction. We combine our mock observations with instrument models for Faint Intergalactic Redshifted Emission Balloon-2 (FIREBall-2; UV balloon spectrograph) and HARMONI (visible and NIR IFU on the ELT) to predict CGM observations with either instrument and optimize target selections and observing strategies. Our results show that Lyα emission from the CGM at a redshift of 0.7 will be observable with FIREBall-2 for bright galaxies (NUV∼18 mag), while metal lines like O vi and C iv will remain challenging to detect. HARMONI is found to be well suited to study the CGM at different redshifts with various tracers.

Accurate tracer particles of baryon dynamics in the adaptive mesh refinement code Ramses

We present a new implementation of the tracer particles algorithm based on a Monte Carlo approach for the Eulerian adaptive mesh refinement code RAMSES. The purpose of tracer particles is to keep track of where fluid elements originate in Eulerian mesh codes, so as to follow their Lagrangian trajectories and re-processing history. We provide a comparison to the more commonly used velocity-based tracer particles, and show that the Monte Carlo approach reproduces the gas distribution much more accurately. We present a detailed statistical analysis of the properties of the distribution of tracer particles in the gas and report that it follows a Poisson law. We extend these Monte Carlo gas tracer particles to tracer particles for the stars and black holes, so that they can exchange mass back and forth between themselves. With such a scheme, we can follow the full cycle of baryons, that is, from gas-forming stars to the release of mass back to the surrounding gas multiple times, or accretion of gas onto black holes. The overall impact on computation time is ∼3% per tracer per initial cell. As a proof of concept, we study an astrophysical science case – the dual accretion modes of galaxies at high redshifts –, which highlights how the scheme yields information hitherto unavailable. These tracer particles will allow us to study complex astrophysical systems where both efficiency of shock-capturing Godunov schemes and a Lagrangian follow-up of the fluid are required simultaneously.

Protostellar birth with ambipolar and ohmic diffusion

Context. The transport of angular momentum is fundamental during the formation of low-mass stars; too little removal and rotation ensures stellar densities are never reached, too much and the absence of rotation means no protoplanetary disks can form. Magnetic diffusion is seen as a pathway to resolving this long-standing problem. Aims. We aim to investigate the impact of including resistive magnetohydrodynamics (MHD) in simulations of the gravitational collapse of a 1 M⊙ gas sphere, from molecular cloud densities to the formation of the protostellar seed; the second Larson core. Methods. We used the adaptive mesh refinement code RAMSES to perform two 3D simulations of collapsing magnetised gas spheres, including self-gravity, radiative transfer in the form of flux-limited diffusion, and a non-ideal gas equation of state to describe H2 dissociation which leads to the second collapse. The first run was carried out under the ideal MHD approximation, while ambipolar and ohmic diffusion was incorporated in the second calculation using resistivities computed from an equilibrium chemical network. Results. In the ideal MHD simulation, the magnetic field dominates the energy budget everywhere inside and around the first hydrostatic core, fueling interchange instabilities and driving a low-velocity outflow above and below the equatorial plane of the system. High magnetic braking removes essentially all angular momentum from the second core. On the other hand, ambipolar and ohmic diffusion create a barrier which prevents amplification of the magnetic field beyond 0.1 G in the first Larson core which is now fully thermally supported. A significant amount of rotation is preserved and a small Keplerian-like disk forms around the second core. The ambipolar and ohmic diffusions are effective at radii below 10 AU, indicating that a least ~1 AU is necessary to investigate the angular momentum transfer and the formation of rotationally supported disks. Finally, when studying the radiative efficiency of the first and second core accretion shocks, we found that it can vary by several orders of magnitude over the 3D surface of the cores. Conclusions. This proves that magnetic diffusion is a prerequisite to star formation. Not only does it enable the formation of protoplanetary disks in which planets will eventually form, it also plays a determinant role in the formation of the protostar itself.

Episodic accretion: the interplay of infall and disc instabilities

Using zoom-simulations carried out with the adaptive mesh-refinement code RAMSES with a dynamic range of up to $2^{27} \approx 1.34 \times 10^8$ we investigate the accretion profiles around six stars embedded in different environments inside a (40 pc)$^3$ giant molecular cloud, the role of mass infall and disc instabilities on the accretion profile, and thus on the luminosity of the forming protostar. Our results show that the environment in which the protostar is embedded determines the overall accretion profile of the protostar. Infall on to the circumstellar disc may trigger gravitational disc instabilities in the disc at distances of around ~10 to ~50 au leading to rapid transport of angular momentum and strong accretion bursts. These bursts typically last for about ~10 to a ~100 yr, consistent with typical orbital times at the location of the instability, and enhance the luminosity of the protostar. Calculations with the stellar evolution code mesa show that the accretion bursts induce significant changes in the protostellar proper- ties, such as the stellar temperature and radius. We apply the obtained protostellar properties to produce synthetic observables with RADMC3D and predict that accretion bursts lead to ob- servable enhancements around 20 to 200 $\mu$m in the spectral energy distribution of Class 0 type young stellar objects.

The fate of the gaseous discs of galaxies that fall into clusters

Galaxy clusters are known to induce gas loss in infalling galaxies due to the ram pressure exerted by the intracluster medium over their gas content. In this paper, we investigate this process through a set of simulations of Milky Way like galaxies falling inside idealised clusters of 10$^{14}$ M$_\odot$ and 10$^{15}$ M$_\odot$, containing a cool-core or not, using the adaptive mesh refinement code RAMSES. We use these simulations to constrain how much of the initial mass contained in the gaseous disk of the galaxy will be converted into stars and how much of it will be lost, after a single crossing of the entire cluster. We find that, if the galaxy reaches the central region of a cool-core cluster, it is expected to lose all its gas, independently of its entry conditions and of the cluster's mass. On the other hand, it is expected to never lose all its gas after crossing a cluster without a cool-core just once. Before reaching the centre of the cluster, the SFR of the galaxy is always enhanced, by a factor of 1.5 to 3. If the galaxy crosses the cluster without being completely stripped, its final amount of gas is on average two times smaller after crossing the 10$^{15}$ M$_\odot$ cluster, relative to the 10$^{14}$ M$_\odot$ cluster. This is reflected in the final SFR of the galaxy, which is also two times smaller in the former, ranging from 0.5 $-$ 1 M$_\odot$ yr$^{-1}$, compared to 1 $-$ 2 M$_\odot$ yr$^{-1}$ for the latter.

Blossoms from black hole seeds: properties and early growth regulated by supernova feedback

Massive black holes (BHs) inhabit local galaxies, including the Milky Way and some dwarf galaxies. BH formation, occurring at early cosmic times, must account for the properties of BHs in today's galaxies, notably why some galaxies host a BH, and others do not. We investigate the formation, distribution and growth of BH `seeds' by using the adaptive mesh refinement code Ramses. We develop an implementation of BH formation in dense, low-metallicity environments, as advocated by models invoking the collapse of the first generation of stars, or of dense nuclear star clusters. The seed masses are computed one-by-one on-the-fly, based on the star formation rate and the stellar initial mass function. This self-consistent method to seed BHs allows us to study the distribution of BHs in a cosmological context and their evolution over cosmic time. We find that all high-mass galaxies tend to a host a BH, whereas low-mass counterparts have a lower probability of hosting a BH. After the end of the epoch of BH formation, this probability is modulated by the growth of the galaxy. The simulated BHs connect to low-redshift observational samples, and span a similar range in accretion properties as Lyman-Break Analogs. The growth of BHs in low-mass galaxies is stunted by strong supernova feedback. The properties of BHs in dwarf galaxies thus remain a testbed for BH formation. Simulations with strong supernova feedback, which is able to quench BH accretion in shallow potential wells, produce galaxies and BHs in better agreement with observational constraints.

FOLLOWING THE COSMIC EVOLUTION OF PRISTINE GAS. I. IMPLICATIONS FOR MILKY WAY HALO STARS

ABSTRACT We make use of a new subgrid model of turbulent mixing to accurately follow the cosmological evolution of the first stars, the mixing of their supernova (SN) ejecta, and the impact on the chemical composition of the Galactic Halo. Using the cosmological adaptive mesh refinement code ramses, we implement a model for the pollution of pristine gas as described in Pan et al. Tracking the metallicity of Pop III stars with metallicities below a critical value allows us to account for the fraction of stars formed even in regions in which the gas’s average metallicity is well above We demonstrate that such partially mixed regions account for 0.5 to 0.7 of all Pop III stars formed up to z = 5. Additionally, we track the creation and transport of “primordial metals” (PM) generated by Pop III SNe. These neutron-capture deficient metals are taken up by second-generation stars and likely lead to unique abundance signatures characteristic of carbon-enhanced, metal-poor (CEMP-no) stars. As an illustrative example, we associate primordial metals with abundance ratios used by Keller et al. to explain the source of metals in the star SMSS J031300.36-670839.3, finding good agreement with the observed [Fe/H], [C/H], [O/H], and [Mg/Ca] ratios in CEMP-no Milky Way halo stars. Similar future simulations will aid in further constraining the properties of Pop III stars using CEMP observations, as well as improve predictions of the spatial distribution of Pop III stars, as will be explored by the next generation of ground- and space-based telescopes.

Comparison of black hole growth in galaxy mergers with gasoline and ramses

Supermassive black hole dynamics during galaxy mergers is crucial in determining the rate of black hole mergers and cosmic black hole growth. As simulations achieve higher resolution, it becomes important to assess whether the black hole dynamics is influenced by the treatment of the interstellar medium in different simulation codes. We here compare simulations of black hole growth in galaxy mergers with two codes: the Smoothed Particle Hydrodynamics code Gasoline, and the Adaptive Mesh Refinement code Ramses. We seek to identify predictions of these models that are robust despite differences in hydrodynamic methods and implementations of sub-grid physics. We find that the general behavior is consistent between codes. Black hole accretion is minimal while the galaxies are well-separated (and even as they "fly-by" within 10 kpc at first pericenter). At late stages, when the galaxies pass within a few kpc, tidal torques drive nuclear gas inflow that triggers bursts of black hole accretion accompanied by star formation. We also note quantitative discrepancies that are model-dependent: our Ramses simulations show less star formation and black hole growth, and a smoother gas distribution with larger clumps and filaments, than our Gasoline simulations. We attribute these differences primarily to the sub-grid models for black hole fueling and feedback and gas thermodynamics. The main conclusion is that differences exist quantitatively between codes, and this should be kept in mind when making comparisons with observations. However, reassuringly, both codes capture the same dynamical behaviors in terms of triggering of black hole accretion, star formation, and black hole dynamics.

Ambipolar diffusion in low-mass star formation

In this paper, we provide a more accurate description of the evolution of the magnetic flux redistribution during prestellar core collapse by including resistive terms in the magnetohydrodynamics (MHD) equations. We focus more particularly on the impact of ambipolar diffusion. We use the adaptive mesh refinement code RAMSES to carry out such calculations. The resistivities required to calculate the ambipolar diffusion terms were computed using a reduced chemical network of charged, neutral and grain species. The inclusion of ambipolar diffusion leads to the formation of a magnetic diffusion barrier in the vicinity of the core, preventing accumulation of magnetic flux in and around the core and amplification of the field above 0.1G. The mass and radius of the first Larson core remain similar between ideal and non-ideal MHD models. This diffusion plateau has crucial consequences on magnetic braking processes, allowing the formation of disk structures. Magnetically supported outflows launched in ideal MHD models are weakened when using non-ideal MHD. Contrary to ideal MHD misalignment between the initial rotation axis and the magnetic field direction does not significantly affect the results for a given mu, showing that the physical dissipation truly dominate over numerical diffusion. We demonstrate severe limits of the ideal MHD formalism, which yield unphysical behaviours in the long-term evolution of the system. This includes counter rotation inside the outflow, interchange instabilities, and flux redistribution triggered by numerical diffusion, none observed in non-ideal MHD. Disks with Keplerian velocity profiles form in all our non-ideal MHD simulations, with final mass and size which depend on the initial magnetisation. This ranges from a few 0.01 solar masses and 20-30 au for the most magnetised case (mu=2) to 0.2 solar masses and 40-80 au for a lower magnetisation (mu=5).