Simulations Of Isolated Galaxies Research Articles

Radial properties of dust in galaxies: Comparison between observations and isolated galaxy simulations

We study the importance of several processes that influence the evolution of dust and its grain size distribution on spatially resolved scales in nearby galaxies. Here, we compiled several multi-wavelength observations for the nearby galaxies NGC\,628 (M74), NGC\,5457 (M101), NGC\,598 (M33), and NGC 300. We applied spatially resolved spectral energy distribution (SED) fitting to the latest iteration of infrared data to get constraints on the galaxy dust masses and the small-to-large grain abundance ratio (SLR). We separated each galaxy into radial rings and obtained the radial profiles of the properties mentioned above. For comparison, we took the radial profiles of the stellar mass and gas mass surface density for NGC\,628 combined with its metallicity gradient in the literature to calibrate a single-galaxy simulation using the GADGET4-OSAKA code. The simulations include a parametrization to separate the dense and diffuse phases of the ISM where different dust-evolution mechanisms are in action. We find that our simulation can reproduce the radial profile of dust mass surface density but overestimates the SLR in NGC\,628. Changing the dust-accretion timescale has little impact on the dust mass or SLR, as most of the available metals are accreted onto dust grains at early times ($<3$\,Gyr), except in the outer regions of the galaxy where the metallicity is below $2 $. This suggests we can only constrain the accretion timescale of galaxies at extremely low metallicities where accretion still competes with other mechanisms controlling the dust budget. The overestimation of the SLR likely results from (i) overly efficient shattering processes in the diffuse interstellar medium (ISM), which were calibrated to reproduce Milky Way-type galaxies and/or (ii) our use of a diffuse and dense gas density subgrid model that does not entirely capture the intricacies of the small-scale structure present in NGC\,628. We conclude that future modeling efforts will need to focus on improving the subgrid recipes to mimic the multi-phase gas distribution in galaxies before the efficiency of dust evolution processes can be calibrated for galaxies other than the Milky Way.

Observational signatures of the dust size evolution in isolated galaxy simulations

We aim to provide observational signatures of the dust size evolution in the interstellar medium (ISM). In particular, we explore indicators of the polycyclic aromatic hydrocarbon (PAH) mass fraction ($q_ PAH $), defined as the mass fraction of PAHs relative to the total dust grains. In addition, we validate our dust evolution model by comparing the observational signatures from our simulations to observations. We used the hydrodynamic simulation code, GADGET4-OSAKA to model the dust properties of Milky Way-like and NGC 628-like galaxies representing star-forming galaxies. This code incorporates the evolution of grain size distributions driven by dust production and interstellar processing. Furthermore, we performed post-processing dust radiative transfer calculations with SKIRT based on the hydrodynamic simulations to predict the observational properties of the simulations. We find that the intensity ratio between 8 um and 24 um ($I_ nu um )/I_ nu um )$) is correlated with $q_ PAH $ and can be used as an indicator of the PAH mass fraction. However, this ratio is influenced by the local radiation field. We suggest the 8 um -to-total infrared intensity ratio ($ I_ nu um )/ I_ TIR $) as another indicator of the PAH mass fraction, since it is tightly correlated with the PAH mass fraction. Furthermore, we explored the spatially resolved evolutionary properties of the PAH mass fraction in the simulated Milky Way-like galaxy using $ I_ nu um )/ I_ TIR $. We find that the spatially resolved PAH mass fraction increases with metallicity at $Z due to the interplay between accretion and shattering, whereas it decreases at $Z because of coagulation. Also, coagulation decreases the PAH mass fraction in regions with a high hydrogen surface density. Finally, we compared the above indicators in the NGC 628-like simulation with those observed in NGC 628 by Herschel Spitzer and JWST . Consequently, we find that our simulation underestimates the PAH mass fraction throughout the entire galaxy by a factor of $ 8$ on average. This could be due to the efficient loss of PAHs by coagulation in our model, suggesting that our treatment of PAHs in dense regions needs to be improved.

Entropy-conserving Scheme for Modeling Nonthermal Energies in Fluid Dynamics Simulations

We compare the performance of energy-based and entropy-conserving schemes for modeling nonthermal energy components, such as unresolved turbulence and cosmic rays, using idealized fluid dynamics tests and isolated galaxy simulations. While both methods are aimed to model advection and adiabatic compression or expansion of different energy components, the energy-based scheme numerically solves the nonconservative equation for the energy density evolution, while the entropy-conserving scheme uses a conservative equation for modified entropy. Using the standard shock tube and Zel’dovich pancake tests, we show that the energy-based scheme results in a spurious generation of nonthermal energy on shocks, while the entropy-conserving method evolves the energy adiabatically to machine precision. We also show that, in simulations of an isolated L ⋆ galaxy, switching between the schemes results in ≈20%–30% changes of the total star formation rate and a significant difference in morphology, particularly near the galaxy center. We also outline and test a simple method that can be used in conjunction with the entropy-conserving scheme to model the injection of nonthermal energies on shocks. Finally, we discuss how the entropy-conserving scheme can be used to capture the kinetic energy dissipated by numerical viscosity into the subgrid turbulent energy implicitly, without explicit source terms that require calibration and can be rather uncertain. Our results indicate that the entropy-conserving scheme is the preferred choice for modeling nonthermal energy components, a conclusion that is equally relevant for Eulerian and moving-mesh fluid dynamics codes.

RAMSES-RTZ: non-equilibrium metal chemistry and cooling coupled to on-the-fly radiation hydrodynamics

ABSTRACT Emission and absorption lines from elements heavier than helium (metals) represent one of our strongest probes of galaxy formation physics across nearly all redshifts accessible to observations. The vast majority of simulations that model these metal lines often assume either collisional or photoionization equilibrium, or a combination of the two. For the few simulations that have relaxed these assumptions, a redshift-dependent meta-galactic UV background or fixed spectrum is often used in the non-equilibrium photoionization calculation, which is unlikely to be accurate in the interstellar medium where the gas can self-shield, as well as in the high-redshift circumgalactic medium, where locally emitted radiation may dominate over the UV background. In this work, we relax this final assumption by coupling the ionization states of individual metals to the radiation hydrodynamics solver present in ramses-rt. Our chemical network follows radiative recombination, dielectronic recombination, collisional ionization, photoionization, and charge transfer, and we use the ionization states to compute non-equilibrium optically thin metal-line cooling. The fiducial model solves for the ionization states of C, N, O, Mg, Si, S, Fe, and Ne in addition to H, He, and H2, but can be easily extended for other ions. We provide interfaces to two different ODE solvers that are competitive in both speed and accuracy. The code has been benchmarked across a variety of gas conditions to reproduce results from cloudy when equilibrium is reached. We show an example isolated galaxy simulation with on-the-fly radiative transfer that demonstrates the utility of our code for translating between simulations and observations without the use of idealized photoionization models.

The global star formation law by supernova feedback

ABSTRACT We address a simple model where the Kennicutt–Schmidt (KS) relation between the macroscopic densities of star formation rate (SFR, ρsfr) and gas (n) in galactic discs emerges from self-regulation of the SFR via supernova feedback. It arises from the physics of supernova bubbles, insensitive to the microscopic SFR recipe and not explicitly dependent on gravity. The key is that the filling factor of SFR-suppressed supernova bubbles self-regulates to a constant, f ∼ 0.5. Expressing the bubble fading radius and time in terms of n, the filling factor is $f\propto S\, n^{-s}$ with s ≃ 1.5, where S is the supernova rate density. A constant f thus refers to ρsfr ∝ n1.5, with a density-independent SFR efficiency per free-fall time ∼0.01. The self-regulation to f ∼ 0.5 and the convergence to a KS relation independent of the local SFR recipe are demonstrated in cosmological and isolated-galaxy simulations using different codes and recipes. In parallel, the spherical analysis of bubble evolution is generalized to clustered supernovae, analytically and via simulations, yielding s ≃ 1.5 ± 0.5. An analysis of photoionized bubbles about pre-supernova stars yields a range of KS slopes but the KS relation is dominated by the supernova bubbles. Superbubble blowouts may lead to an alternative self-regulation by outflows and recycling. While the model is oversimplified, its simplicity and validity in the simulations may argue that it captures the origin of the KS relation.

The echo of the bar buckling: Phase-space spirals in Gaia Data Release 2

We present a high-resolution numerical study of the phase-space diversity in an isolated Milky Way-type galaxy. Using a single N-body simulation (N ≈ 0.14 × 109) we explore the formation, evolution, and spatial variation of the phase-space spirals similar to those recently discovered by Antoja et al. in the Milky Way disk with Gaia Data Release 2 (DR2). For the first time in the literature we use a self-consistent N-body simulation of an isolated Milky Way-type galaxy to show that the phase-space spirals develop naturally from vertical oscillations driven by the buckling of the stellar bar. Thus, we claim that the physical mechanism standing behind the observed incomplete phase-space mixing process can be internal and not necessarily due to the perturbation induced by a massive satellite. In our model, the bending oscillations propagate outward and produce axisymmetric variations of the mean vertical coordinate and vertical velocity component of about 100 − 200 pc and 1 − 2 km s−1, respectively. As a consequence, the phase-space wrapping results in the formation of patterns with various morphologies across the disk, depending on the bar orientation, distance to the galactic center, and time elapsed since the bar buckling. Once bending waves appear, they are supported for a long time via disk self-gravity. Such vertical oscillations trigger the formation of various time-dependent phase-space spirals in the entire disk. The underlying physical mechanism implies the link between in-plane and vertical motion that leads directly to phase-space structures whose amplitude and shape are in remarkable agreement with those of the phase-space spirals observed in the Milky Way disk. In our isolated galaxy simulation, phase-space spirals are still distinguishable at the solar neighborhood 3 Gyr after the buckling phase. The long-lived character of the phase-space spirals generated by the bar buckling instability cast doubts on the timing argument used so far to get back to the time of the onset of the perturbation: phase-space spirals may have been caused by perturbations originated several gigayearrs ago, and not as recent as suggested so far.

Osaka feedback model: isolated disc galaxy simulations

We study various implementations of supernova feedback model and present the results of our `Osaka feedback model' using isolated galaxy simulations performed by the smoothed particle hydrodynamics (SPH) code {\small GADGET-3}. Our model is a modified version of Stinson et al.'s work, and we newly add the momentum kick for SN feedback rather than only thermal feedback. We incorporate the physical properties of SN remnants from the results of Chevalier and McKee \& Ostriker, such as the effective radius of SN bubble and the remnant life-time, in the form of Sedov-Taylor (ST)-like solutions with the effect of radiative cooling. Our model utilizes the local, physical parameters such as density and temperature of the ISM rather than galactic or halo properties to determine the galactic wind velocity or mass-loading factor. The Osaka model succeeds in self-regulating star formation, and naturally produces galactic outflow with variable velocities depending on the local environment and available SN energy as a function of time.An important addition to our previous work by Aoyama et al. is the implementation of the {\small CELib} chemistry library which allows us to deal with the time-dependent input of energy and metal yields for type Ia \& II supernovae (SNe) and asymptotic giant branch (AGB) stars. As initial tests of our model, we apply it to isolated galaxy simulations, and examine various galactic properties and compare with observational data including metal abundances.

The impact of galactic disc environment on star-forming clouds

We explore the effect of different galactic disc environments on the properties of star-forming clouds through variations in the background potential in a set of isolated galaxy simulations. Rising, falling and flat rotation curves expected in halo dominated, disc dominated and Milky Way-like galaxies were considered, with and without an additional two-arm spiral potential. The evolution of each disc displayed notable variations that are attributed to different regimes of stability, determined by shear and gravitational collapse. The properties of a typical cloud were largely unaffected by the changes in rotation curve, but the production of small and large cloud associations was strongly dependent on this environment. This suggests that while differing rotation curves can influence where clouds are initially formed, the average bulk properties are effectively independent of the global environment. The addition of a spiral perturbation made the greatest difference to cloud properties, successfully sweeping the gas into larger, seemingly unbound, extended structures and creating large arm-interarm contrasts.

Black hole feeding and feedback: the physics inside the ‘sub-grid’

Black holes (BHs) are believed to be a key ingredient of galaxy formation. However, the galaxy-BH interplay is challenging to study due to the large dynamical range and complex physics involved. As a consequence, hydrodynamical cosmological simulations normally adopt sub-grid models to track the unresolved physical processes, in particular BH accretion; usually the spatial scale where the BH dominates the hydrodynamical processes (the Bondi radius) is unresolved, and an approximate Bondi-Hoyle accretion rate is used to estimate the growth of the BH. By comparing hydrodynamical simulations at different resolutions (300, 30, 3 pc) using a Bondi-Hoyle approximation to sub-parsec runs with non-parameterized accretion, our aim is to probe how well an approximated Bondi accretion is able to capture the BH accretion physics and the subsequent feedback on the galaxy. We analyse an isolated galaxy simulation that includes cooling, star formation, Type Ia and Type II supernovae, BH accretion and AGN feedback (radiation pressure, Compton heating/cooling) where mass, momentum, and energy are deposited in the interstellar medium through conical winds. We find that on average the approximated Bondi formalism can lead to both over- and under-estimations of the BH growth, depending on resolution and on how the variables entering into the Bondi-Hoyle formalism are calculated.

Black hole growth and AGN feedback under clumpy accretion

High-resolution simulations of supermassive black holes in isolated galaxies have suggested the importance of short (~10 Myr) episodes of rapid accretion caused by interactions between the black hole and massive dense clouds within the host. Accretion of such clouds could potentially provide the dominant source for black hole growth in high-z galaxies, but it remains unresolved in cosmological simulations. Using a stochastic subgrid model calibrated by high-resolution isolated galaxy simulations, we investigate the impact that variability in black hole accretion rates has on black hole growth and the evolution of the host galaxy. We find this clumpy accretion to more efficiently fuel high-redshift black hole growth. This increased mass allows for more rapid accretion even in the absence of high-density clumps, compounding the effect and resulting in substantially faster overall black hole growth. This increased growth allows the black hole to efficiently evacuate gas from the central region of the galaxy, driving strong winds up to ~2500 km/s, producing outflows ~10x stronger than the smooth accretion case, suppressing the inflow of gas onto the host galaxy, and suppressing the star formation within the galaxy by as much as a factor of two. This suggests that the proper incorporation of variability is a key factor in the co-evolution between black holes and their hosts.

The anatomy of a star-forming galaxy: pressure-driven regulation of star formation in simulated galaxies

We explore the regulation of star formation in star-forming galaxies through a suite of high-resolution isolated galaxy simulations. We use the SPH code GASOLINE, including photoelectric heating and metal cooling, which produces a multi-phase interstellar medium. We show that representative star formation and feedback sub-grid models naturally lead to a weak, sub-linear dependence between the amount of star formation and changes to star formation parameters. We incorporate these sub-grid models into an equilibrium pressure-driven regulation framework. We show that the sub-linear scaling arises as a consequence of the non-linear relationship between scale height and the effective pressure generated by stellar feedback. Thus, simulated star-formation regulation is sensitive to how well vertical structure in the ISM is resolved. Full galaxy disks experience density waves which drive locally time-dependent star formation. We develop a simple time-dependent, pressure-driven model that reproduces the response extremely well.

BLACK HOLE GROWTH AND ACTIVE GALACTIC NUCLEI OBSCURATION BY INSTABILITY-DRIVEN INFLOWS IN HIGH-REDSHIFT DISK GALAXIES FED BY COLD STREAMS

Disk galaxies at high redshift have been predicted to maintain high gas surface densities due to continuous feeding by intense cold streams leading to violent gravitational instability, transient features and giant clumps. Gravitational torques between the perturbations drive angular momentum out and mass in, and the inflow provides the energy for keeping strong turbulence. We use analytic estimates of the inflow for a self-regulated unstable disk at a Toomre stability parameter Q~1, and isolated galaxy simulations capable of resolving the nuclear inflow down to the central parsec. We predict an average inflow rate ~10 Msun/yr through the disk of a 10^11 Msun galaxy, with conditions representative of z~2 stream-fed disks. The inflow rate scales with disk mass and (1+z)^{3/2}. It includes clump migration and inflow of the smoother component, valid even if clumps disrupt. This inflow grows the bulge, while only a fraction ~ 10^-3 of it needs to accrete onto a central black hole (BH), in order to obey the observed BH-bulge relation. A galaxy of 10^11 Msun at z~2 is expected to host a BH of ~10^8 Msun, accreting on average with moderate sub-Eddington luminosity L_X ~ 10^42-43 erg/s, accompanied by brighter episodes when dense clumps coalesce. We note that in rare massive galaxies at z~6, the same process may feed 10^9 Msun BH at the Eddington rate. High central gas column densities can severely obscure AGN in high-redshift disks, possibly hindering their detection in deep X-ray surveys.