Search for: All records

Turbulent radiative mixing layers (TRMLs) form at the interface of cold, dense gas and hot, diffuse gas in motion with each other. TRMLs are ubiquitous in and around galaxies on a variety of scales, including galactic winds and the circumgalactic medium. They host the intermediate-temperature gases that are efficient in radiative cooling, thus playing a crucial role in controlling the cold gas supply, phase structure, and spectral features of galaxies. In this work, we develop an intuitive analytic 1.5-dimensional model for TRMLs that includes a simple parameterization of the effective turbulent conductivity and viscosity and a piecewise power-law cooling curve. Our analytic model reproduces the mass flux, total cooling, and phase structure of 3D simulations of TRMLs at a fraction of the computational cost. It also reveals essential insights into the physics of TRMLs, particularly the importance of the viscous dissipation of relative kinetic energy in balancing radiative cooling as the shear Mach number approaches unity. This dissipation takes place both in the intermediate-temperature phase, which reduces the enthalpy flux from the hot phase, and in the cold phase, which enhances radiative cooling. Additionally, our model provides a fast and easy way of computing the column density and surface brightness of TRMLs, which can be directly linked to observations.

 

The driving of turbulence in galaxies is deeply connected with the physics of feedback, star formation, outflows, accretion, and radial transport in disks. The velocity dispersion of gas in galaxies therefore offers a promising observational window into these processes. However, the relative importance of each of these mechanisms remains controversial. In this work we revisit the possibility that turbulence on galactic scales is driven by the direct impact of accreting gaseous material on the disk. We measure this effect in a disk-like star-forming galaxy in IllustrisTNG, using the high-resolution cosmological magnetohydrodynamical simulation TNG50. We employ Lagrangian tracer particles with a high time cadence of only a few million years to identify accretion and other events. The energies of particles are measured by stacking the events in bins of time around the event. The average effect of each event is measured by fitting explicit models for the kinetic and turbulent energies as a function of time. These measurements are corroborated by cross-correlating the turbulent energy with other time series and searching for signals of causality, i.e., asymmetries across zero time lag. We find that accretion contributes to the large-scale turbulent kinetic energy even if it does not dominate in this ∼5 × 10⁹M_⊙stellar mass galaxy. Extrapolating this finding to a range of galaxy masses, we find that there are regimes where energy from direct accretion may dominate the turbulent energy budget, particularly in disk outskirts, galaxies less massive than the Milky Way, and at redshift ∼2.

 

Galactic outflows driven by supernovae (SNe) are thought to be a powerful regulator of a galaxy’s star-forming efficiency. Mass, energy, and metal outflows (η_M,η_E, andη_Z, here normalized by the star formation rate, the SNe energy, and metal production rates, respectively) shape galaxy properties by both ejecting gas and metals out of the galaxy and by heating the circumgalactic medium (CGM), preventing future accretion. Traditionally, models have assumed that galaxies self-regulate by ejecting a large fraction of the gas, which enters the interstellar medium (ISM), although whether such high mass loadings agree with observations is still unclear. To better understand how the relative importance of ejective (i.e., high mass loading) versus preventative (i.e., high energy loading) feedback affects the present-day properties of galaxies, we develop a simple gas-regulator model of galaxy evolution, where the stellar mass, ISM, and CGM are modeled as distinct reservoirs which exchange mass, metals, and energy at different rates within a growing halo. Focusing on the halo mass range from 10¹⁰to 10¹²M_⊙, we demonstrate that, with reasonable parameter choices, we can reproduce the stellar-to-halo mass relation and the ISM-to-stellar mass relation with low-mass-loaded (η_M∼ 0.1–10) but high-energy-loaded (η_E∼ 0.1–1) winds, with self-regulation occurring primarily through heating and cooling of the CGM. We show that the model predictions are robust against changes to the mass loading of outflows but are quite sensitive to our choice of the energy loading, preferringη_E∼ 1 for the lowest-mass halos and ∼0.1 for Milky Way–like halos.

 

In a novel approach employing implicit likelihood inference (ILI), also known as likelihood-free inference, we calibrate the parameters of cosmological hydrodynamic simulations against observations, which has previously been unfeasible due to the high computational cost of these simulations. For computational efficiency, we train neural networks as emulators on ∼1000 cosmological simulations from the CAMELS project to estimate simulated observables, taking as input the cosmological and astrophysical parameters, and use these emulators as surrogates for the cosmological simulations. Using the cosmic star formation rate density (SFRD) and, separately, the stellar mass functions (SMFs) at different redshifts, we perform ILI on selected cosmological and astrophysical parameters (Ω_m,σ₈, stellar wind feedback, and kinetic black hole feedback) and obtain full six-dimensional posterior distributions. In the performance test, the ILI from the emulated SFRD (SMFs) can recover the target observables with a relative error of 0.17% (0.4%). We find that degeneracies exist between the parameters inferred from the emulated SFRD, confirmed with new full cosmological simulations. We also find that the SMFs can break the degeneracy in the SFRD, which indicates that the SMFs provide complementary constraints for the parameters. Further, we find that a parameter combination inferred from an observationally inferred SFRD reproduces the target observed SFRD very well, whereas, in the case of the SMFs, the inferred and observed SMFs show significant discrepancies that indicate potential limitations of the current galaxy formation modeling and calibration framework, and/or systematic differences and inconsistencies between observations of the SMFs.

 

The dividing line between galaxies that are quenched by reionization (“relics”) and galaxies that survive reionization (i.e., continue forming stars) is commonly discussed in terms of a halo mass threshold. We probe this threshold in a physically more complete and accurate way than has been possible to date, using five extremely high resolution (M_target= 4M_⊙) cosmological zoom-in simulations of dwarf galaxies within the halo mass range (1–4) × 10⁹M_⊙. The employed LYRA simulation model features resolved interstellar medium physics and individual, resolved supernova explosions. Interestingly, two out of five of the simulated dwarf galaxies lie close to the threshold mass but are neither full reionization relics nor full reionization survivors. These galaxies initially quench at the time of reionization but merely remain quiescent for ∼500 Myr. Atz∼ 5 they recommence star formation in a synchronous way and remain star-forming until the present day. The parallel timing indicates consistent sound-crossing and cooling times between the halos. While the star formation histories we find are diverse, we show that they are directly related to the ability of a given halo to retain and cool gas. Whereas the latter is most strongly dependent on the mass (or virial temperature) of the host halo at the time of reionization, it also depends on its growth history, the UV background (and its decrease at late times), and the amount of metals retained within the halo.

 

Cluster spiral galaxies suffer catastrophic losses of the cool, neutral gas component of their interstellar medium due to ram pressure stripping, contributing to the observed quenching of star formation in the disc compared to galaxies in lower density environments. However, the short-term effects of ram pressure on the star formation rate and active galactic nucleus (AGN) activity of galaxies undergoing stripping remain unclear. Numerical studies have recently demonstrated cosmic rays can dramatically influence galaxy evolution for isolated galaxies, yet their influence on ram pressure stripping remains poorly constrained. We perform the first cosmic ray magnetohydrodynamic simulations of an L* galaxy undergoing ram pressure stripping, including radiative cooling, self-gravity of the gas, star formation, and stellar feedback. We find the microscopic transport of cosmic rays plays a key role in modulating the star formation enhancement experienced by spirals at the outskirts of clusters compared to isolated spirals. Moreover, we find that galaxies undergoing ram pressure stripping exhibit enhanced gas accretion on to their centres, which may explain the prevalence of AGNs in these objects. In agreement with observations, we find cosmic rays significantly boost the global radio emission of cluster spirals. Although the gas removal rate is relatively insensitive to cosmic ray physics, we find that cosmic rays significantly modify the phase distribution of the remaining gas disc. These results suggest observations of galaxies undergoing ram pressure stripping may place novel constraints on cosmic ray calorimetry and transport.

 

Galaxy mergers are known to host abundant young massive cluster (YMC) populations, whose formation mechanism is still not well-understood. Here, we present a high-resolution galaxy merger simulation with explicit star formation and stellar feedback prescriptions to investigate how mergers affect the properties of the interstellar medium and YMCs. Compared with a controlled simulation of an isolated galaxy, the mass fraction of dense and high-pressure gas is much higher in mergers. Consequently, the mass function of both molecular clouds and YMCs becomes shallower and extends to higher masses. Moreover, cluster formation efficiency is significantly enhanced and correlates positively with the star formation rate surface density and gas pressure. We track the orbits of YMCs and investigate the time evolution of tidal fields during the course of the merger. At an early stage of the merger, the tidal field strength correlates positively with YMC mass, λtid ∝ M0.71, which systematically affects the shape of the mass function and age distribution of the YMCs. At later times, most YMCs closely follow the orbits of their host galaxies, gradually sinking into the centre of the merger remnant due to dynamical friction, and are quickly dissolved via efficient tidal disruption. Interestingly, YMCs formed during the first passage, mostly in tidal tails and bridges, are distributed over a wide range of galactocentric radii, greatly increasing their survivability because of the much weaker tidal field in the outskirts of the merger system. These YMCs are promising candidates for globular clusters that survive to the present day.

 

Abstract We present a suite of high-resolution simulations of an isolated dwarf galaxy using four different hydrodynamical codes: Gizmo , Arepo , Gadget , and Ramses . All codes adopt the same physical model, which includes radiative cooling, photoelectric heating, star formation, and supernova (SN) feedback. Individual SN explosions are directly resolved without resorting to subgrid models, eliminating one of the major uncertainties in cosmological simulations. We find reasonable agreement on the time-averaged star formation rates as well as the joint density–temperature distributions between all codes. However, the Lagrangian codes show significantly burstier star formation, larger SN-driven bubbles, and stronger galactic outflows compared to the Eulerian code. This is caused by the behavior in the dense, collapsing gas clouds when the Jeans length becomes unresolved: Gas in Lagrangian codes collapses to much higher densities than that in Eulerian codes, as the latter is stabilized by the minimal cell size. Therefore, more of the gas cloud is converted to stars and SNe are much more clustered in the Lagrangian models, amplifying their dynamical impact. The differences between Lagrangian and Eulerian codes can be reduced by adopting a higher star formation efficiency in Eulerian codes, which significantly enhances SN clustering in the latter. Adopting a zero SN delay time reduces burstiness in all codes, resulting in vanishing outflows as SN clustering is suppressed. 

The MeerKAT Fornax Survey maps the distribution and kinematics of atomic neutral hydrogen gas (H  I ) in the nearby Fornax galaxy cluster using the MeerKAT telescope. The 12 deg 2 survey footprint covers the central region of the cluster out to ∼ R vir and stretches south-west out to ∼2 R vir to include the NGC 1316 galaxy group. The H  I column density sensitivity (3 σ over 25 km s −1 ) ranges from 5 × 10 19 cm −2 at a resolution of ∼10″ (∼1 kpc at the 20 Mpc distance of Fornax) down to ∼10 18 cm −2 at ∼1′ (∼6 kpc), and slightly below this level at the lowest resolution of ∼100″ (∼10 kpc). The H  I mass sensitivity (3 σ over 50 km s −1 ) is 6 × 10 5 M ⊙ . The H  I velocity resolution is 1.4 km s −1 . In this paper, we describe the survey design and H  I data processing, and we present a sample of six galaxies with long, one-sided, starless H  I tails (only one of which was previously known) radially oriented within the cluster and with measurable internal velocity gradients. We argue that the joint properties of the H  I tails represent the first unambiguous evidence of ram pressure shaping the distribution of H  I in the Fornax cluster. The disturbed optical morphology of all host galaxies supports the idea that the tails consist of H  I that was initially pulled out of the galaxies’ stellar body by tidal forces. Ram pressure was then able to further displace the weakly bound H  I and give the tails their current direction, length, and velocity gradient. 

Abstract Physical and chemical properties of the interstellar medium (ISM) at subgalactic (∼kiloparsec) scales play an indispensable role in controlling the ability of gas to form stars. In this paper, we use the TNG50 cosmological simulation to explore the physical parameter space of eight resolved ISM properties in star-forming regions to constrain the areas of this hyperspace where most star-forming environments exist. We deconstruct our simulated galaxies spanning a wide range of mass ( M ⋆ = 10 7 –10 11 M ⊙ ) and redshift (0 ≤ z ≤ 3) into kiloparsec-sized regions and statistically analyze the gas/stellar surface densities, gas metallicity, vertical stellar velocity dispersion, epicyclic frequency, and dark-matter volumetric density representative of each region in the context of their star formation activity and environment (radial galactocentric location). By examining the star formation rate (SFR) weighted distributions of these properties, we show that stars primarily form in two distinct environmental regimes, which are brought about by an underlying bicomponent radial SFR profile in galaxies. We examine how the relative prominence of these regimes depends on galaxy mass and cosmic time. We also compare our findings with those from integral field spectroscopy observations and find similarities as well as departures. Further, using dimensionality reduction, we characterize the aforementioned hyperspace to reveal a high degree of multicollinearity in relationships among ISM properties that drive the distribution of star formation at kiloparsec scales. Based on this, we show that a reduced 3D representation underpinned by a multivariate radius relationship is sufficient to capture most of the variance in the original 8D space.