Midplane Pressure Research Articles

Metallicity Dependence of Pressure-regulated Feedback-modulated Star Formation in the TIGRESS-NCR Simulation Suite

We present a new suite of numerical simulations of the star-forming interstellar medium (ISM) in galactic disks using the TIGRESS-NCR framework. Distinctive aspects of our simulation suite are (1) sophisticated and comprehensive numerical treatments of essential physical processes including magnetohydrodynamics, self-gravity, and galactic differential rotation, as well as photochemistry, cooling, and heating coupled with direct ray-tracing UV radiation transfer and resolved supernova feedback and (2) wide parameter coverage including the variation in metallicity over Z′≡Z/Z⊙∼0.1-3 , gas surface density Σgas ∼ 5–150 M ⊙ pc−2, and stellar surface density Σstar ∼ 1–50 M ⊙ pc−2. The range of emergent star formation rate surface density is ΣSFR ∼ 10−4–0.5 M ⊙ kpc−2 yr−1, and ISM total midplane pressure is P tot/k B = 103–106 cm−3 K, with P tot equal to the ISM weight W . For given Σgas and Σstar, we find ΣSFR∝Z′0.3 . We provide an interpretation based on the pressure-regulated feedback-modulated (PRFM) star formation theory. The total midplane pressure consists of thermal, turbulent, and magnetic stresses. We characterize feedback modulation in terms of the yield ϒ, defined as the ratio of each stress to ΣSFR. The thermal feedback yield varies sensitively with both weight and metallicity as ϒth∝W−0.46Z′−0.53 , while the combined turbulent and magnetic feedback yield shows weaker dependence ϒturb+mag∝W−0.22Z′−0.18 . The reduction in ΣSFR at low metallicity is due mainly to enhanced thermal feedback yield, resulting from reduced attenuation of UV radiation. With the metallicity-dependent calibrations we provide, PRFM theory can be used for a new subgrid star formation prescription in cosmological simulations where the ISM is unresolved.

On the distribution of the cold neutral medium in galaxy discs

ABSTRACT The cold neutral medium (CNM) is an important part of the galactic gas cycle and a precondition for the formation of molecular and star-forming gas, yet its distribution is still not fully understood. In this work, we present extremely high resolution simulations of spiral galaxies with time-dependent chemistry such that we can track the formation of the CNM, its distribution within the galaxy, and its correlation with star formation. We find no strong radial dependence between the CNM fraction and total neutral atomic hydrogen (H i) due to the decreasing interstellar radiation field counterbalancing the decreasing gas column density at larger galactic radii. However, the CNM fraction does increase in spiral arms where the CNM distribution is clumpy, rather than continuous, overlapping more closely with H2. The CNM does not extend out radially as far as H i, and the vertical scale height is smaller in the outer galaxy compared to H i with no flaring. The CNM column density scales with total mid-plane pressure and disappears from the gas phase below values of PT/kB = 1000 K cm−3. We find that the star formation rate density follows a similar scaling law with CNM column density to the total gas Kennicutt–Schmidt law. In the outer galaxy, we produce realistic vertical velocity dispersions in the H i purely from galactic dynamics, but our models do not predict CNM at the extremely large radii observed in H i absorption studies of the Milky Way. We suggest that extended spiral arms might produce isolated clumps of CNM at these radii.

Saturation of the magnetorotational instability and the origin of magnetically elevated accretion discs

ABSTRACT We propose that the strength of angular momentum transport in accretion discs threaded by net vertical magnetic field is determined by a self-regulation mechanism: the magnetorotational instability (MRI) grows until its own turbulent resistivity damps the fastest growing mode on the scale of the disc thickness. Given weak assumptions as to the structure of MRI-derived turbulence, supported by prior simulation evidence, the proposed mechanism reproduces the known scaling of the viscous α-parameter, $\alpha \propto \beta _z^{-1/2}$. Here, $\beta _z = 8\pi p_{\rm g}/B_{z0}^2$ is the initial plasma β-parameter on the disc mid-plane, Bz0 is the net field, and pg is the mid-plane gas pressure. We generalize the argument to discs with strong suprathermal toroidal magnetic fields, where the MRI growth rate is modified from the weak-field limit. Additional sources of turbulence are required if such discs are to become magnetically elevated, with the increased scale heights near the mid-plane that are seen in simulations. We speculate that tearing modes, associated with current sheets broadened by the effective resistivity, are a possible source of enhanced turbulence in elevated discs.

Introducing TIGRESS-NCR. I. Coregulation of the Multiphase Interstellar Medium and Star Formation Rates

Massive, young stars are the main source of energy that maintains multiphase structure and turbulence in the interstellar medium (ISM), and without this “feedback” the star formation rate (SFR) would be much higher than is observed. Rapid energy loss in the ISM and efficient energy recovery by stellar feedback lead to coregulation of SFRs and the ISM state. Realistic approaches to this problem should solve for the dynamical evolution of the ISM, including star formation and the input of feedback energy self-consistently and accurately. Here, we present the TIGRESS-NCR numerical framework, in which UV radiation, supernovae, cooling and heating processes, and gravitational collapse are modeled explicitly. We use an adaptive ray-tracing method for UV radiation transfer from star clusters represented by sink particles, accounting for attenuation by dust and gas. We solve photon-driven chemical equations to determine the abundances of hydrogen (time dependent) and carbon/oxygen-bearing species (steady state), which then set cooling and heating rates self-consistently. Applying these methods, we present high-resolution magnetohydrodynamics simulations of differentially rotating local galactic disks representing typical conditions of nearby star-forming galaxies. We analyze ISM properties and phase distributions and show good agreement with existing multiwavelength galactic observations. We measure midplane pressure components (turbulent, thermal, and magnetic) and the weight, demonstrating that vertical dynamical equilibrium holds. We quantify the ratios of pressure components to the SFR surface density, which we call the feedback yields. The TIGRESS-NCR framework will allow for a wide range of parameter exploration, including in low-metallicity systems.

Ultra-diffuse Galaxies as Extreme Star-forming Environments. II. Star Formation and Pressure Balance in H i-rich UDGs

In addition to occupying the extreme, diffuse tail of the dwarf galaxy population, ultra-diffuse galaxies (UDGs) are themselves a key laboratory in which to study star formation in extreme low-density environments. In the second paper of this series, we compare the spatially resolved star formation activity of 22 H i-selected UDGs and 21 “normal” dwarf galaxies within 120 Mpc to predictions within the pressure-regulated, feedback-modulated (PRFM) theory of star formation. To do so, we employ a joint spectral energy distribution fitting method that allows us to estimate star formation rate and stellar mass surface density from UV-optical imaging. We find that the PRFM framework extends successfully to the UDG regime—although the UDGs in our sample show unusually low star formation rate surface densities given their H i content, this low star formation efficiency can be naturally explained by the diffuse structure of the UDGs. In fact, when cast in the PRFM framework, the relationship between midplane pressure and star formation in the UDG sample is in good agreement not only with the “normal” dwarf reference sample, but also with measurements from more massive galaxies. Our results suggest that despite their low star formation efficiencies, the H i-rich UDGs need not be forming stars in an exotic manner. We also find that the UDGs are likely H2 poor compared even to the overall dwarf population.

Pressure-regulated, Feedback-modulated Star Formation in Disk Galaxies

The star formation rate (SFR) in galactic disks depends on both the quantity of the available interstellar medium (ISM) gas and its physical state. Conversely, the ISM’s physical state depends on the SFR, because the “feedback” energy and momentum injected by recently formed massive stars is crucial to offsetting losses from turbulent dissipation and radiative cooling. The ISM’s physical state also responds to the gravitational field that confines it, with increased weight driving higher pressure. In a quasi-steady state, it is expected that the mean total pressure of different thermal phases will match each other, that the component pressures and total pressure will satisfy thermal and dynamical equilibrium requirements, and that the SFR will adjust as needed to provide the requisite stellar radiation and supernova feedback. The pressure-regulated, feedback-modulated (PRFM) theory of the star-forming ISM formalizes these ideas, leading to a prediction that the SFR per unit area, ΣSFR, will scale nearly linearly with ISM weight . In terms of the large-scale gas surface density Σgas, stellar plus dark matter density ρ sd, and effective ISM velocity dispersion σ eff, an observable weight estimator is , and this is predicted to match the total midplane pressure P tot. Using a suite of multiphase magnetohydrodynamic simulations run with the TIGRESS computational framework, we test the principles of the PRFM model and calibrate the total feedback yield ϒtot = P tot/ΣSFR ∼ 1000 km s−1, as well as its components. We compare the results from TIGRESS to theory, previous numerical simulations, and observations, finding excellent agreement.

Predicting HCN, HCO+, multi-transition CO, and dust emission of star-forming galaxies

The interstellar medium is a turbulent, multiphase, and multi-scale medium that follows scaling relations that link the surface density, volume density, and velocity dispersion with the cloud size. Galactic clouds range from below 1 pc to about 100 pc in size. Extragalactic clouds appear to follow the same range, although they are only now becoming observable in atomic and molecular lines. Analytical models of galactic gaseous disks need to take the multi-scale and multiphase nature of the interstellar medium into account. They can be described as clumpy star-forming accretion disks in vertical hydrostatic equilibrium, with the mid-plane pressure balancing the gravity of the gaseous and stellar disk. Interstellar medium turbulence is taken into account by applying Galactic scaling relations to the cold atomic and molecular gas phases. Turbulence is maintained through energy injection by supernovae. With the determination of the gas mass fraction at a given spatial scale, the equilibrium gas temperature between turbulent heating and line cooling, the molecular abundances, and the molecular line emission can be calculated. The resulting model radial profiles of infrared, HI, CO, HCN, and HCO+ emission are compared to THINGS, HERACLES, EMPIRE, SINGS, and GALEX observations of 17 local spiral galaxies. The model free parameters were constrained for each galactic radius independently. The Toomre parameter, which measures the stability against star formation (cloud collapse), exceeds unity in the inner disk of a significant number of galaxies. In two galaxies it also exceeds unity in the outer disk. Therefore, in spiral galaxies Qtot = 1 is not ubiquitous. The model gas velocity dispersion is consistent with the observed HI velocity dispersion where available. Within our model, HCN and HCO+ is already detectable in relatively low-density gas (∼1000 cm−3). We derive CO and HCN conversion factors and molecular gas depletion times. Both conversion factors are consistent with values found in the literature. Whereas in the massive galaxies the viscous timescale greatly exceeds the star-formation timescale, the viscous timescale is smaller than the star-formation timescale within R ∼ 2 Rd, the disk scale length, in the low-mass galaxies. We suggest that massive spiral galaxies undergo starvation in the absence of gas accretion from the halo, whereas in low-mass galaxies the fuel for star formation reaches R ∼ 2 Rd from outside via a thick gas disk component with a high radial infall velocity observable in the HI line.

Effects of Varying Mass Inflows on Star Formation in Nuclear Rings of Barred Galaxies

Abstract Observations indicate that the star formation rate (SFR) of nuclear rings varies considerably with time and is sometimes asymmetric rather than being uniform across a ring. To understand what controls temporal and spatial distributions of ring star formation, we run semiglobal, hydrodynamic simulations of nuclear rings subject to time-varying and/or asymmetric mass inflow rates. These controlled variations in the inflow lead to variations in the star formation, while the ring orbital period (18 Myr) and radius (600 pc) remain approximately constant. We find that both the mass inflow rate and supernova feedback affect the ring SFR. An oscillating inflow rate with period Δτ in and amplitude 20 causes large-amplitude (a factor of ≳5), quasi-periodic variations of the SFR when Δτ in ≳ 50 Myr. We find that the time-varying interstellar medium (ISM) weight and midplane pressure track each other closely, establishing an instantaneous vertical equilibrium. The measured time-varying depletion time is consistent with the prediction from self-regulation theory provided the time delay between star formation and supernova feedback is taken into account. The supernova feedback is responsible only for small-amplitude (a factor of ∼2) fluctuations of the SFR with a timescale ≲40 Myr. Asymmetry in the inflow rate does not necessarily lead to asymmetric star formation in nuclear rings. Only when the inflow rate from one dust lane is suddenly increased by a large factor do the rings undergo a transient period of lopsided star formation.

Star Formation in Nuclear Rings with the TIGRESS Framework

Abstract Nuclear rings are sites of intense star formation at the centers of barred galaxies. To understand what determines the structure and star formation rate (SFR; ) of nuclear rings, we run semi-global, hydrodynamic simulations of nuclear rings subject to constant mass inflow rates . We adopt the TIGRESS framework of Kim & Ostriker to handle radiative heating and cooling, star formation, and related supernova (SN) feedback. We find that the SN feedback is never strong enough to destroy the ring or quench star formation everywhere in the ring. Under the constant , the ring star formation is very steady and persistent, with the SFR exhibiting only mild temporal fluctuations. The ring SFR is tightly correlated with the inflow rate as , for a range of . Within the ring, vertical dynamical equilibrium is maintained, with the midplane pressure (powered by SN feedback) balancing the weight of the overlying gas. The SFR surface density is correlated nearly linearly with the midplane pressure, as predicted by the pressure-regulated, feedback-modulated star formation theory. Based on our results, we argue that the ring SFR is causally controlled by , while the ring gas mass adapts to the SFR to maintain the vertical dynamical equilibrium under the gravitational field arising from both gas and stars.

Cosmic-Ray Diffusion Suppression in Star-forming Regions Inhibits Clump Formation in Gas-rich Galaxies

Abstract Observations of the γ-ray emission around star clusters, isolated supernova remnants, and pulsar wind nebulae indicate that the cosmic-ray (CR) diffusion coefficient near acceleration sites can be suppressed by a large factor compared to the Galaxy average. We explore the effects of such local suppression of CR diffusion on galaxy evolution using simulations of isolated disk galaxies with regular and high gas fractions. Our results show that while CR propagation with constant diffusivity can make gaseous disks more stable by increasing the midplane pressure, large-scale CR pressure gradients cannot prevent local fragmentation when the disk is unstable. In contrast, when CR diffusivity is suppressed in star-forming regions, the accumulation of CRs in these regions results in strong local pressure gradients that prevent the formation of massive gaseous clumps. As a result, the distribution of dense gas and star formation changes qualitatively: a globally unstable gaseous disk does not violently fragment into massive star-forming clumps but maintains a regular grand-design spiral structure. This effect regulates star formation and disk structure and is qualitatively different from and complementary to the global role of CRs in vertical hydrostatic support of the gaseous disk and in driving galactic winds.

The EDGE-CALIFA survey: self-regulation of star formation at kpc scales

ABSTRACT The processes that regulate star formation are essential to understand how galaxies evolve. We present the relation between star formation rate density, ΣSFR , and hydrostatic mid-plane pressure,Ph, for 4260 star-forming regions of kpc size located in 96 galaxies included in the EDGE-CALIFA survey covering a wide range of stellar masses and morphologies. We find that these two parameters are tightly correlated, showing a smaller scatter in comparison to other star-forming relations. A power law, with a slightly sublinear index, is a good representation of this relation. Its residuals show a significant anticorrelation with both stellar age and metallicity whereas the total stellar mass may also play a secondary role in shaping the ΣSFR–Ph relation. For actively star-forming regions, we find that the effective feedback momentum per unit stellar mass (p*/m*), measured from thePh/ΣSFR ratio increases withPh. The median value of this ratio for all the sampled regions is larger than the expected momentum just from supernovae explosions. Morphology of the galaxies, including bars, does not seem to have a significant impact in the ΣSFR–Ph relation. Our analysis indicates that local ΣSFR self-regulation comes mainly from momentum injection to the interstellar medium from supernovae explosions. However, other mechanisms in disc galaxies may also play a significant role in shaping the ΣSFR at kpc scales. Our results also suggest thatPh is the main parameter that modulates star formation at kpc scales, rather than individual components of the baryonic mass.

First Results from SMAUG: Characterization of Multiphase Galactic Outflows from a Suite of Local Star-forming Galactic Disk Simulations

Abstract Large-scale outflows in star-forming galaxies are observed to be ubiquitous and are a key aspect of theoretical modeling of galactic evolution, the focus of the Simulating Multiscale Astrophysics to Understand Galaxies (SMAUG) project. Gas blown out from galactic disks, similar to gas within galaxies, consists of multiple phases with large contrasts of density, temperature, and other properties. To study multiphase outflows as emergent phenomena, we run a suite of rougly parsec-resolution local galactic disk simulations using the TIGRESS framework. Explicit modeling of the interstellar medium (ISM), including star formation and self-consistent radiative heating plus supernova feedback, regulates ISM properties and drives the outflow. We investigate the scaling of outflow mass, momentum, energy, and metal loading factors with galactic disk properties, including star formation rate (SFR) surface density (ΣSFR ∼ 10−4 − 1 M ⊙ kpc−2 yr−1), gas surface density ( ), and total midplane pressure (or weight; ). The main components of outflowing gas are mass-delivering cool gas (T ∼ 104 K) and energy/metal-delivering hot gas (T ≳ 106 K). Cool mass outflow rates measured at outflow launch points (one or two scale heights ) are 1–100 times the SFR (decreasing with ΣSFR), although in massive galaxies most mass falls back owing to insufficient outflow velocity. The hot galactic outflow carries mass comparable to 10% of the SFR, together with 10%–20% of the energy and 30%–60% of the metal mass injected by SN feedback. Importantly, our analysis demonstrates that in any physically motivated cosmological wind model it is crucial to include at least two distinct thermal wind components.

Pressure balance in the multiphase ISM of cosmologically simulated disc galaxies

ABSTRACT Pressure balance plays a central role in models of the interstellar medium (ISM), but whether and how pressure balance is realized in a realistic multiphase ISM is not yet well understood. We address this question by using a set of FIRE-2 cosmological zoom-in simulations of Milky Way-mass disc galaxies, in which a multiphase ISM is self-consistently shaped by gravity, cooling, and stellar feedback. We analyse how gravity determines the vertical pressure profile as well as how the total ISM pressure is partitioned between different phases and components (thermal, dispersion/turbulence, and bulk flows). We show that, on average and consistent with previous more idealized simulations, the total ISM pressure balances the weight of the overlying gas. Deviations from vertical pressure balance increase with increasing galactocentric radius and with decreasing averaging scale. The different phases are in rough total pressure equilibrium with one another, but with large deviations from thermal pressure equilibrium owing to kinetic support in the cold and warm phases, which dominate the total pressure near the mid-plane. Bulk flows (e.g. inflows and fountains) are important at a few disc scale heights, while thermal pressure from hot gas dominates at larger heights. Overall, the total mid-plane pressure is well-predicted by the weight of the disc gas and we show that it also scales linearly with the star formation rate surface density (ΣSFR). These results support the notion that the Kennicutt–Schmidt relation arises because ΣSFR and the gas surface density (Σg) are connected via the ISM mid-plane pressure.

Simulations of Unsteady Three-Dimensional Hypersonic Double-Wedge Flow Experiments

Hypersonic flow over a three-dimensional 30–55 deg double wedge is investigated for three different conditions of varying enthalpy and Reynolds number with nitrogen as the test gas. Streamwise and spanwise wall heat flux are compared to the experimental data that were obtained at the University of Illinois at Urbana–Champaign and the California Institute of Technology. In previous studies, there were concerns that portions of the flowfield were rarefied; however, after an in-depth analysis, we find the flowfield to be firmly in the continuum regime. Comparisons of heat flux at an instant of time in the simulations show good agreement with the experiments. However, averaged heat flux over a long period of time relative to the experimental test time are underpredicted and develop low-frequency periodic motion. To better understand the unsteadiness and three-dimensionality of the flowfield, we analyze the mass flux through a spanwise plane, and we perform sparsity-promoting dynamic mode decomposition for the wall heat flux and midplane pressure. Finally, we analyze the effect of vibrational nonequilibrium on the flowfield. It is found that the flowfield dynamics vary with increasing enthalpy and Reynolds number. Furthermore, the flowfield is found to be three dimensional, unsteady, and asymmetric for all cases.

The role of galactic dynamics in shaping the physical properties of giant molecular clouds in Milky Way-like galaxies

ABSTRACT We examine the role of the large-scale galactic-dynamical environment in setting the properties of giant molecular clouds in Milky Way-like galaxies. We perform three high-resolution simulations of Milky Way-like discs with the moving-mesh hydrodynamics code arepo, yielding a statistical sample of ${\sim}80\, 000$ giant molecular clouds and ${\sim}55\, 000$ H i clouds. We account for the self-gravity of the gas, momentum, and thermal energy injection from supernovae and H ii regions, mass injection from stellar winds, and the non-equilibrium chemistry of hydrogen, carbon, and oxygen. By varying the external gravitational potential, we probe galactic-dynamical environments spanning an order of magnitude in the orbital angular velocity, gravitational stability, mid-plane pressure, and the gradient of the galactic rotation curve. The simulated molecular clouds are highly overdense (∼100×) and overpressured (∼25×) relative to the ambient interstellar medium. Their gravoturbulent and star-forming properties are decoupled from the dynamics of the galactic mid-plane, so that the kpc-scale star formation rate surface density is related only to the number of molecular clouds per unit area of the galactic mid-plane. Despite this, the clouds display clear, statistically significant correlations of their rotational properties with the rates of galactic shearing and gravitational free-fall. We find that galactic rotation and gravitational instability can influence their elongation, angular momenta, and tangential velocity dispersions. The lower pressures and densities of the H i clouds allow for a greater range of significant dynamical correlations, mirroring the rotational properties of the molecular clouds, while also displaying a coupling of their gravitational and turbulent properties to the galactic-dynamical environment.

Local Simulations of Spiral Galaxies with the TIGRESS Framework. I. Star Formation and Arm Spurs/Feathers

Abstract Spiral arms greatly affect gas flows and star formation in disk galaxies. We use local three-dimensional simulations of the vertically stratified, self-gravitating, differentially rotating interstellar medium (ISM) subject to a stellar spiral potential to study the effects of spiral arms on star formation and formation of arm spurs/feathers. We adopt the TIGRESS framework of Kim & Ostriker to handle radiative heating and cooling, star formation, and ensuing supernova (SN) feedback. We find that more than 90% of star formation takes place in spiral arms, but the global star formation rate (SFR) in models with spiral arms is enhanced by less than a factor of 2 compared to the no-arm counterpart. This results from a quasi-linear relationship between the SFR surface density and the gas surface density Σ, and supports the picture that spiral arms do not trigger star formation but rather concentrate star-forming regions. Correlated SN feedback produces gaseous spurs/feathers downstream from arms in both magnetized and unmagnetized models. These spurs/feathers are short lived and have magnetic fields parallel to their length, in contrast to the longer-lived features with perpendicular magnetic fields induced by gravitational instability. SN feedback drives the turbulent component of magnetic fields, with the total magnetic field strength sublinearly proportional to Σ. The total midplane pressure varies by a factor of ∼10 between arm and interarm regions but agrees locally with the total vertical ISM weight, while is locally consistent with the prediction of pressure-regulated, feedback-modulated theory.

Testing Feedback-regulated Star Formation in Gas-rich, Turbulent Disk Galaxies

Abstract In this paper we compare the molecular gas depletion times and midplane hydrostatic pressure in turbulent, star-forming disk galaxies to internal properties of these galaxies. For this analysis we use 17 galaxies from the DYNAMO sample of nearby (z ∼ 0.1) turbulent disks. We find a strong correlation, such that galaxies with lower molecular gas depletion time (t dep) have higher gas velocity dispersion (σ). Within the scatter of our data, our observations are consistent with the prediction that made in theories of feedback-regulated star formation. We also show a strong, single power-law correlation between midplane pressure (P) and star formation rate surface density (ΣSFR), which extends for 6 orders of magnitude in pressure. Disk galaxies with lower pressure are found to be roughly in agreement with theoretical predictions. However, in galaxies with high pressure we find P/ΣSFR values that are significantly larger than theoretical predictions. Our observations could be explained with any of the following: (1) the correlation of ΣSFR−P is significantly sublinear; (2) the momentum injected from star formation feedback (p */m *) is not a single, universal value; or (3) alternate sources of pressure support are important in gas-rich disk galaxies. Finally, using published survey results, we find that our results are consistent with the cosmic evolution of t dep(z) and σ(z). Our interpretation of these results is that the cosmic evolution of t dep may be regulated not just by the supply of gas but also by the internal regulation of star formation via feedback.

Breathing oscillations in a global simulation of a thin accretion disc

We study the oscillations of an axisymmetric, viscous, radiative, general relativistic hydrodynamical simulation of a geometrically thin disk around a non-rotating, $6.62\,M_\odot$ black hole. The numerical setup is initialized with a Novikov-Thorne, gas-pressure-dominated accretion disk, with an initial mass accretion rate of $\dot{m} = 0.01\,L_\mathrm{Edd}/c^2$ (where $L_\mathrm{Edd}$ is the Eddington luminosity and $c$ is the speed of light). Viscosity is treated with the $\alpha$-prescription. The simulation was evolved for about $1000$ Keplerian orbital periods at three Schwarzschild radii (ISCO radius). Power density spectra of the radial and vertical fluid velocity components, the total (gas $+$ radiation) midplane pressure, and the vertical component of radiative flux from the photosphere, all reveal strong power at the local breathing oscillation frequency. The first, second and third harmonics of the breathing oscillation are also clearly seen in the data. We quantify the properties of these oscillations by extracting eigenfunctions of the radial and vertical velocity components and total pressure. This confirms that these oscillations are associated with breathing motion.

Modeling the Atomic-to-molecular Transition in Cosmological Simulations of Galaxy Formation

Abstract Large-scale cosmological simulations of galaxy formation currently do not resolve the densities at which molecular hydrogen forms, implying that the atomic-to-molecular transition must be modeled either on the fly or in postprocessing. We present an improved postprocessing framework to estimate the abundance of atomic and molecular hydrogen and apply it to the IllustrisTNG simulations. We compare five different models for the atomic-to-molecular transition, including empirical, simulation-based, and theoretical prescriptions. Most of these models rely on the surface density of neutral hydrogen and the ultraviolet (UV) flux in the Lyman–Werner band as input parameters. Computing these quantities on the kiloparsec scale resolved by the simulations emerges as the main challenge. We show that the commonly used Jeans length approximation to the column density of a system can be biased and exhibits large cell-to-cell scatter. Instead, we propose to compute all surface quantities in face-on projections and perform the modeling in two dimensions. In general, the two methods agree on average, but their predictions diverge for individual galaxies and for models based on the observed midplane pressure of galaxies. We model the UV radiation from young stars by assuming a constant escape fraction and optically thin propagation throughout the galaxy. With these improvements, we find that the five models for the atomic-to-molecular transition roughly agree on average but that the details of the modeling matter for individual galaxies and the spatial distribution of molecular hydrogen. We emphasize that the estimated molecular fractions are approximate due to the significant systematic uncertainties.

Investigations of the static and dynamic characteristics of water-lubricated hydrodynamic journal bearing considering turbulent, thermohydrodynamic and misaligned effects

Abstract A misaligned thermohydrodynamic (THD) model considering turbulent-flow is presented to predict the performance of water-lubricated hydrodynamic journal bearing (WHJB). The load capacity, frictional power loss, flow rate, mid-plane pressure distribution, translational, tilting and coupled translational-tilting dynamic coefficients versus the rotary speed, eccentricity ratio and tilting angle are analyzed; quantitative comparisons between models considering and without considering turbulent, THD and misaligned effects are given. Results show that the turbulence improves the load capacity while the THD effect decreases it; the THD effect reduces the frictional power loss significantly while the turbulent effect is weak; at large rotary speed and eccentricity ratio, the THD and turbulent effects on the dynamic coefficients should be fully considered at both of aligned and misaligned conditions.