Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “explainability”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 937 records · Page 52

The Utility of Macroecological Rules for Microbial Biogeography

Macroecological rules have been developed for plants and animals that describe large-scale distributional patterns and attempt to explain the underlying physiological and ecological processes behind them. Similarly, microorganisms exhibit patterns in relative abundance, distribution, diversity, and traits across space and time, yet it remains unclear the extent to which microorganisms follow macroecological rules initially developed for macroorganisms. Additionally, the usefulness of these rules as a null hypothesis when surveying microorganisms has yet to be fully evaluated. With rapid advancements in sequencing technology, we have seen a recent increase in microbial studies that utilize macroecological frameworks. Here, we review and synthesize these macroecological microbial studies with two main objectives: (1) to determine to what extent macroecological rules explain the distribution of host-associated and free-living microorganisms, and (2) to understand which environmental factors and stochastic processes may explain these patterns among microbial clades (archaea, bacteria, fungi, and protists) and habitats (host-associated and free living; terrestrial and aquatic). Overall, 78% of microbial macroecology studies focused on free living, aquatic organisms. In addition, most studies examined macroecological rules at the community level with only 35% of studies surveying organismal patterns across space. At the community level microorganisms often tracked patterns of macroorganisms for island biogeography (74% confirm) but rarely followed Latitudinal Diversity Gradients (LDGs) of macroorganisms (only 32% confirm). However, when microorganisms and macroorganisms shared the same macroecological patterns, underlying environmental drivers (e.g., temperature) were the same. Because we found a lack of studies for many microbial groups and habitats, we conclude our review by outlining several outstanding questions and creating recommendations for future studies in microbial ecology.

59 BASIC BIOLOGICAL SCIENCES↗

Genome-Wide Association Study for Identifying Loci that Affect Fillet Yield, Carcass, and Body Weight Traits in Rainbow Trout (Oncorhynchus mykiss)

Fillet yield (FY, %) is an economically-important trait in rainbow trout aquaculture that affects production efficiency. Despite that, FY has received little attention in breeding programs because it is difficult to measure on a large number of fish and cannot be directly measured on breeding candidates. The recent development of a high-density SNP array for rainbow trout has provided the needed tool for studying the underlying genetic architecture of this trait. A genome-wide association study (GWAS) was conducted for FY, body weight at 10 (BW10) and 13 (BW13) months post-hatching, head-off carcass weight (CAR), and fillet weight (FW) in a pedigreed rainbow trout population selectively bred for improved growth performance. The GWAS analysis was performed using the weighted single-step GBLUP method (wssGWAS). Phenotypic records of 1447 fish (1.5 kg at harvest) from 299 full-sib families in three successive generations, of which 875 fish from 196 full-sib families were genotyped, were used in the GWAS analysis. A total of 38,107 polymorphic SNPs were analyzed in a univariate model with hatch year and harvest group as fixed effects, harvest weight as a continuous covariate, and animal and common environment as random effects. A new linkage map was developed to create windows of 20 adjacent SNPs for use in the GWAS. The two windows with largest effect for FY and FW were located on chromosome Omy9 and explained only 1.0–1.5% of genetic variance, thus suggesting a polygenic architecture affected by multiple loci with small effects in this population. One window on Omy5 explained 1.4 and 1.0% of the genetic variance for BW10 and BW13, respectively. Three windows located on Omy27, Omy17, and Omy9 (same window detected for FY) explained 1.7, 1.7, and 1.0%, respectively, of genetic variance for CAR. Among the detected 100 SNPs, 55% were located directly in genes (intron and exons). Nucleotide sequences of intragenic SNPs were blasted to the Mus musculus genome to create a putative gene network. The network suggests that differences in the ability to maintain a proliferative and renewable population of myogenic precursor cells may affect variation in growth and fillet yield in rainbow trout.

54 ENVIRONMENTAL SCIENCES↗

Genomic Diversity Evaluation of Populus trichocarpa Germplasm for Rare Variant Genetic Association Studies

Genome-wide association studies are powerful tools to elucidate the genome-to-phenome relationship. In order to explain most of the observed heritability of a phenotypic trait, a sufficient number of individuals and a large set of genetic variants must be examined. The development of high-throughput technologies and cost-efficient resequencing of complete genomes have enabled the genome-wide identification of genetic variation at large scale. As such, almost all existing genetic variation becomes available, and it is now possible to identify rare genetic variants in a population sample. Rare genetic variants that were usually filtered out in most genetic association studies are the most numerous genetic variations across genomes and hold great potential to explain a significant part of the missing heritability observed in association studies. Rare genetic variants must be identified with high confidence, as they can easily be confounded with sequencing errors. In this study, we used a pre-filtered data set of 1,014 pure Populus trichocarpa entire genomes to identify rare and common small genetic variants across individual genomes. We compared variant calls between Platypus and HaplotypeCaller pipelines, and we further applied strict quality filters for improved genetic variant identification. Finally, we only retained genetic variants that were identified by both variant callers increasing calling confidence. Based on these shared variants and after stringent quality filtering, we found high genomic diversity in P. trichocarpa germplasm, with 7.4 million small genetic variants. Importantly, 377k non-synonymous variants (5% of the total) were uncovered. We highlight the importance of genomic diversity and the potential of rare defective genetic variants in explaining a significant portion of P. trichocarpa's phenotypic variability in association genetics. The ultimate goal is to associate both rare and common alleles with poplar's wood quality traits to support selective breeding for an improved bioenergy feedstock.

Genetics & Heredity↗

Modeling HIV-1 Within-Host Dynamics After Passive Infusion of the Broadly Neutralizing Antibody VRC01

VRC01 is a broadly neutralizing antibody that targets the CD4 binding site of HIV-1 gp120. Passive administration of VRC01 in humans has assessed the safety and the effect on plasma viremia of this monoclonal antibody (mAb) in a phase 1 clinical trial. After VRC01 infusion, the plasma viral load in most of the participants was reduced but had particular dynamics not observed during antiretroviral therapy. In this paper, we introduce different mathematical models to explain the observed dynamics and fit them to the plasma viral load data. Based on the fitting results we argue that a model containing reversible Ab binding to virions and clearance of virus-VRC01 complexes by a two-step process that includes (1) saturable capture followed by (2) internalization/degradation by phagocytes, best explains the data. This model predicts that VRC01 may enhance the clearance of Ab-virus complexes, explaining the initial viral decay observed immediately after antibody infusion in some participants. Because Ab-virus complexes are assumed to be unable to infect cells, i.e., contain neutralized virus, the model predicts a longer-term viral decay consistent with that observed in the VRC01 treated participants. By assuming a homogeneous viral population sensitive to VRC01, the model provides good fits to all of the participant data. However, the fits are improved by assuming that there were two populations of virus, one more susceptible to antibody-mediated neutralization than the other.

60 APPLIED LIFE SCIENCES↗

GRAIN SIZE CONSTRAINTS ON HL TAU WITH POLARIZATION SIGNATURE

The millimeter-wave polarization of the protoplanetary disk around HL Tau has been interpreted as the emission from elongated dust grains aligned with the magnetic field in the disk. However, the self-scattering of thermal dust emission may also explain the observed millimeter-wave polarization. In this paper, we report a modeling of the millimeter-wave polarization of the HL Tau disk with the self-polarization. Dust grains are assumed to be spherical and to have a power-law size distribution. We change the maximum grain size with a fixed dust composition in a fixed disk model to find the grain size to reproduce the observed signature. We find that the direction of the polarization vectors and the polarization degree can be explained with the self-scattering. Moreover, the polarization degree can be explained only if the maximum grain size is ∼150 μm. The obtained grain size from the polarization is different from that which has been previously expected from the spectral index of the dust opacity coefficient (a millimeter or larger) if the emission is optically thin. We discuss that porous dust aggregates may solve the inconsistency of the maximum grain size between the two constraints.

79 ASTRONOMY AND ASTROPHYSICS↗

THE INFLUENCE OF NONUNIFORM CLOUD COVER ON TRANSIT TRANSMISSION SPECTRA

We model the impact of nonuniform cloud cover on transit transmission spectra. Patchy clouds exist in nearly every solar system atmosphere, brown dwarfs, and transiting exoplanets. Our major findings suggest that fractional cloud coverage can exactly mimic high mean molecular weight atmospheres and vice versa over certain wavelength regions, in particular, over the Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) bandpass (1.1–1.7 μm). We also find that patchy cloud coverage exhibits a signature that is different from uniform global clouds. Furthermore, we explain analytically why the “patchy cloud-high mean molecular weight” degeneracy exists. We also explore the degeneracy of nonuniform cloud coverage in atmospheric retrievals on both synthetic and real planets. We find from retrievals on a synthetic solar composition hot Jupiter with patchy clouds and a cloud-free high mean molecular weight warm Neptune that both cloud-free high mean molecular weight atmospheres and partially cloudy atmospheres can explain the data equally well. Another key finding is that the HST WFC3 transit transmission spectra of two well-observed objects, the hot Jupiter HD 189733b and the warm Neptune HAT-P-11b, can be explained well by solar composition atmospheres with patchy clouds without the need to invoke high mean molecular weight or global clouds. The degeneracy between high molecular weight and solar composition partially cloudy atmospheres can be broken by observing the molecular Rayleigh scattering differences between the two. Furthermore, the signature of partially cloudy limbs also appears as a ∼100 ppm residual in the ingress and egress of the transit light curves, provided that the transit timing is known to seconds.

79 ASTRONOMY AND ASTROPHYSICS↗

THE B/C AND SUB-IRON/IRON COSMIC RAY RATIOS—FURTHER EVIDENCE IN FAVOR OF THE SPIRAL-ARM DIFFUSION MODEL

The boron to carbon (B/C) and sub-Fe/Fe ratios provide an important clue on cosmic ray (CR) propagation within the Galaxy. These ratios estimate the grammage that the CRs traverse as they propagate from their sources to Earth. Attempts to explain these ratios within the standard CR propagation models require ad hoc modifications and even with those these models necessitate inconsistent grammages to explain both ratios. As an alternative, physically motivated model, we have proposed that CRs originate preferably within the galactic spiral arms. CR propagation from dynamic spiral arms has important imprints on various secondary to primary ratios, such as the B/C ratio and the positron fraction. We use our spiral-arm diffusion model with the spallation network extended up to nickel to calculate the sub-Fe/Fe ratio. We show that without any additional parameters the spiral-arm model consistently explains both ratios with the same grammage, providing further evidence in favor of this model.

79 ASTRONOMY AND ASTROPHYSICS↗

TEMPORAL VARIABILITY FROM THE TWO-COMPONENT ADVECTIVE FLOW SOLUTION AND ITS OBSERVATIONAL EVIDENCE

In the propagating oscillatory shock model, the oscillation of the post-shock region, i.e., the Compton cloud, causes the observed low-frequency quasi-periodic oscillations (QPOs). The evolution of QPO frequency is explained by the systematic variation of the Compton cloud size, i.e., the steady radial movement of the shock front, which is triggered by the cooling of the post-shock region. Thus, analysis of the energy-dependent temporal properties in different variability timescales can diagnose the dynamics and geometry of accretion flows around black holes. We study these properties for the high-inclination black hole source XTE J1550-564 during its 1998 outburst and the low-inclination black hole source GX 339-4 during its 2006–07 outburst using RXTE /PCA data, and we find that they can satisfactorily explain the time lags associated with the QPOs from these systems. We find a smooth decrease of the time lag as a function of time in the rising phase of both sources. In the declining phase, the time lag increases with time. We find a systematic evolution of QPO frequency and hard lags in these outbursts. In XTE J1550-564, the lag changes from hard to soft (i.e., from a positive to a negative value) at a crossing frequency (ν {sub c}) of ∼3.4 Hz. We present possible mechanisms to explain the lag behavior of high and low-inclination sources within the framework of a single two-component advective flow model.

79 ASTRONOMY AND ASTROPHYSICS↗

LITHIUM DEPLETION IS A STRONG TEST OF CORE-ENVELOPE RECOUPLING

Rotational mixing is a prime candidate for explaining the gradual depletion of lithium from the photospheres of cool stars during the main sequence. However, previous mixing calculations have relied primarily on treatments of angular momentum transport in stellar interiors incompatible with solar and stellar data in the sense that they overestimate the internal differential rotation. Instead, recent studies suggest that stars are strongly differentially rotating at young ages but approach a solid body rotation during their lifetimes. We modify our rotating stellar evolution code to include an additional source of angular momentum transport, a necessary ingredient for explaining the open cluster rotation pattern, and examine the consequences for mixing. We confirm that core-envelope recoupling with a ∼20 Myr timescale is required to explain the evolution of the mean rotation pattern along the main sequence, and demonstrate that it also provides a more accurate description of the Li depletion pattern seen in open clusters. Recoupling produces a characteristic pattern of efficient mixing at early ages and little mixing at late ages, thus predicting a flattening of Li depletion at a few Gyr, in agreement with the observed late-time evolution. Using Li abundances we argue that the timescale for core-envelope recoupling during the main sequence decreases sharply with increasing mass. We discuss the implications of this finding for stellar physics, including the viability of gravity waves and magnetic fields as agents of angular momentum transport. We also raise the possibility of intrinsic differences in initial conditions in star clusters using M67 as an example.

79 ASTRONOMY AND ASTROPHYSICS↗

Resolving the discrepancy of galaxy merger fraction measurements at z ∼ 0-3

We measure the merger fraction of massive galaxies using the UltraVISTA/COSMOS catalog, complemented with the deeper, higher resolution 3DHST+CANDELS catalog, presenting the largest mass-complete photometric merger sample up to z∼3. We find that the variation in the mass ratio probe can explain the discrepant redshift evolution of the merger fraction in the literature: selecting mergers using the H {sub 160}-band flux ratio leads to an increasing merger fraction with redshift, while selecting mergers using the stellar mass ratio reveals a merger fraction with little redshift dependence at z=1--3. Defining major and minor mergers as having stellar mass ratios of 1:1–4:1 and 4:1–10:1, respectively, the results imply ∼1 major merger and ∼0.7 minor merger on average for a massive (log (M{sub ⋆}/M{sub ⊙})⩾10.8) galaxy during z=0.1--2.5. There may be an additional ∼0.5(0.4) major (minor) merger if we use the H-band flux ratio selection. The observed amount of major merging alone is sufficient to explain the observed number density evolution for the very massive (log (M{sub ⋆}/M{sub ⊙})⩾11.1) galaxies. The observed number of major and minor mergers can increase the size of a massive quiescent galaxy by a factor of two at most. This amount of merging is enough to bring the compact quiescent galaxies formed at z>2 to lie at 1.5σ below the mean of the stellar mass–size relation as measured in some works (e.g., Newman et al.), but additional mechanisms are needed to fully explain the evolution, and to be consistent with works suggesting stronger evolution.

79 ASTRONOMY AND ASTROPHYSICS↗

Dynamical model for spindown of solar-type stars

After their formation, stars slow down their rotation rates by the removal of angular momentum from their surfaces, e.g., via stellar winds. Explaining how this rotation of solar-type stars evolves in time is currently an interesting but difficult problem in astrophysics. Despite the complexity of the processes involved, a traditional model, where the removal of angular momentum by magnetic fields is prescribed, has provided a useful framework to understand observational relations between stellar rotation, age, and magnetic field strength. Here, for the first time, a spindown model is proposed where loss of angular momentum by magnetic fields evolves dynamically, instead of being prescibed kinematically. To this end, we evolve the stellar rotation and magnetic field simultaneously over stellar evolution time by extending our previous work on a dynamo model which incorporates nonlinear feedback mechanisms on rotation and magnetic fields. We show that our extended model reproduces key observations and is capable of explaining the presence of the two branches of (fast and slow rotating) stars which have different relations between rotation rate Ω versus time (age), magnetic field strength |B| versus rotation rate, and frequency of magnetic field ω{sub cyc} versus rotation rate. For fast rotating stars we find that: (i) there is an exponential spindown Ω∝e{sup −1.35t}, with t measured in Gyr; (ii) magnetic activity saturates for higher rotation rate; (iii) ω{sub cyc}∝Ω{sup 0.83}. For slow rotating stars we find: (i) a power-law spindown Ω∝t{sup −0.52}; (ii) that magnetic activity scales roughly linearly with rotation rate; (iii) ω{sub cyc}∝Ω{sup 1.16}. The results obtained from our investigations are in good agreement with observations. The Vaughan–Preston gap is consistently explained in our model by the shortest spindown timescale in this transition from fast to slow rotators. Our results highlight the importance of self-regulation of magnetic fields and rotation by direct and indirect interactions involving nonlinear feedback in stellar evolution.

79 ASTRONOMY AND ASTROPHYSICS↗

TOI-1634 b: An Ultra-short-period Keystone Planet Sitting inside the M-dwarf Radius Valley

Studies of close-in planets orbiting M dwarfs have suggested that the M-dwarf radius valley may be well explained by distinct formation timescales between enveloped terrestrials and rocky planets that form at late times in a gas-depleted environment. This scenario is at odds with the picture that close-in rocky planets form with a primordial gaseous envelope that is subsequently stripped away by some thermally driven mass-loss process. These two physical scenarios make unique predictions of the rocky/enveloped transition’s dependence on orbital separation such that studying the compositions of planets within the M-dwarf radius valley may be able to establish the dominant physics. Here, we present the discovery of one such keystone planet: the ultra-short-period planet TOI-1634 b (P = 0.989 days, F=121F{sub ⊕}, r{sub p}=1.790{sub -0.081}{sup +0.080} R {sub ⊕}) orbiting a nearby M2 dwarf (K {sub s} = 8.7, R {sub s} = 0.450 R {sub ⊙}, M {sub s} = 0.502 M {sub ⊙}) and whose size and orbital period sit within the M-dwarf radius valley. We confirm the TESS-discovered planet candidate using extensive ground-based follow-up campaigns, including a set of 32 precise radial velocity measurements from HARPS-N. We measure a planetary mass of 4.91{sub -0.70}{sup +0.68} M {sub ⊕}, which makes TOI-1634 b inconsistent with an Earth-like composition at 5.9σ and thus requires either an extended gaseous envelope, a large volatile-rich layer, or a rocky composition that is not dominated by iron and silicates to explain its mass and radius. The discovery that the bulk composition of TOI-1634 b is inconsistent with that of Earth supports the gas-depleted formation mechanism to explain the emergence of the radius valley around M dwarfs with M{sub s}≲0.5 M {sub ⊙}.

36 MATERIALS SCIENCE↗

H I-to-H{sub 2} transition layers in the star-forming region W43

The process of atomic-to-molecular (H i-to-H{sub 2}) gas conversion is fundamental for molecular-cloud formation and star formation. 21 cm observations of the star-forming region W43 revealed extremely high H i column densities, of 120–180 M{sub ⊙} pc{sup −2}, a factor of 10–20 larger than predicted by H i-to-H{sub 2} transition theories. We analyze the observed H i with a theoretical model of the H i-to-H{sub 2} transition, and show that the discrepancy between theory and observation cannot be explained by the intense radiation in W43, nor be explained by variations of the assumed volume density or H{sub 2} formation rate coefficient. We show that the large observed H i columns are naturally explained by several (9–22) H i-to-H{sub 2} transition layers, superimposed along the sightlines of W43. We discuss other possible interpretations such as a non-steady-state scenario and inefficient dust absorption. The case of W43 suggests that H i thresholds reported in extragalactic observations are probably not associated with a single H i-to-H{sub 2} transition, but are rather a result of several transition layers (clouds) along the sightlines, beam-diluted with diffuse intercloud gas.

79 ASTRONOMY AND ASTROPHYSICS↗

Gamma–Gamma Absorption in the γ-ray Binary System PSR B1259-63/LS 2883

The observed TeV light curve from the γ -ray binary PSR B1259-63/LS 2883 shows a decrease in the flux at periastron that has not been fully explained by emission mechanisms alone. This observed decrease can, however, be explained by γγ absorption due to the stellar and disk photons. We calculate the γγ absorption in PSR B1259-63/LS 2883 taking into account photons from both the circumstellar disk and star, assuming that the γ -rays originate at the position of the pulsar. The γγ absorption due to the circumstellar disk photons produces a ≈14% decrease in the flux, and there is a total decrease of ≈52% (>1 TeV) within a few days before periastron, accompanied by a hardening of the γ -ray photon index. While the γγ absorption alone is not sufficient to explain the full complexity of the H.E.S.S. γ -ray light curve, it results in a significant decrease in the predicted flux, which is coincident with the observed decrease. In addition, we have calculated an upper limit on the γγ absorption, assuming that the emission is produced at the apex of the bow shock. Future observations with CTA during the 2021 periastron passage may be able to confine the location of the emission based on the degree of γγ absorption, as well as measure the hardening of the spectrum around periastron.

79 ASTRONOMY AND ASTROPHYSICS↗

General Model for Light Curves of Chromospherically Active Binary Stars

The starspots on the surface of many chromospherically active binary stars concentrate on long-lived active longitudes separated by 180°. Shifts in activity between these two longitudes, the “flip-flop” events, have been observed in single stars like FK Comae and binary stars like σ Geminorum. Recently, interferometry has revealed that ellipticity may at least partly explain the flip-flop events in σ Geminorum. This idea was supported by the double-peaked shape of the long-term mean light curve of this star. Here we show that the long-term mean light curves of 14 chromospherically active binaries follow a general model that explains the connection between orbital motion, changes in starspot distribution, ellipticity, and flip-flop events. Surface differential rotation is probably weak in these stars, because the interference of two constant period waves may explain the observed light curve changes. These two constant periods are the active longitude period (P{sub act}) and the orbital period (P{sub orb}). We also show how to apply the same model to single stars, where only the value of P {sub act} is known. Finally, we present a tentative interference hypothesis about the origin of magnetic fields in all spectral types of stars.

79 ASTRONOMY AND ASTROPHYSICS↗

Age Spreads and the Temperature Dependence of Age Estimates in Upper Sco

Past estimates for the age of the Upper Sco Association are typically 11–13 Myr for intermediate-mass stars and 4–5 Myr for low-mass stars. In this study, we simulate populations of young stars to investigate whether this apparent dependence of estimated age on spectral type may be explained by the star formation history of the association. Solar and intermediate mass stars begin their pre-main sequence evolution on the Hayashi track, with fully convective interiors and cool photospheres. Intermediate-mass stars quickly heat up and transition onto the radiative Henyey track. As a consequence, for clusters in which star formation occurs on a timescale similar to that of the transition from a convective to a radiative interior, discrepancies in ages will arise when ages are calculated as a function of temperature instead of mass. Simple simulations of a cluster with constant star formation over several Myr may explain about half of the difference in inferred ages versus photospheric temperature; speculative constructions that consist of a constant star formation followed by a large supernova-driven burst could fully explain the differences, including those between F and G stars where evolutionary tracks may be more accurate. The age spreads of low-mass stars predicted from these prescriptions for star formation are consistent with the observed luminosity spread of Upper Sco. The conclusion that a lengthy star formation history will yield a temperature dependence in ages is expected from the basic physics of pre-main sequence evolution, and is qualitatively robust to the large uncertainties in pre-main sequence evolutionary models.

79 ASTRONOMY AND ASTROPHYSICS↗

Structure and Evolution of Internally Heated Hot Jupiters

Hot Jupiters receive strong stellar irradiation, producing equilibrium temperatures of 1000--2500 K. Incoming irradiation directly heats just their thin outer layer, down to pressures of ∼0.1 bars. In standard irradiated evolution models of hot Jupiters, predicted transit radii are too small. Previous studies have shown that deeper heating—at a small fraction of the heating rate from irradiation—can explain observed radii. Here we present a suite of evolution models for HD 209458b, where we systematically vary both the depth and intensity of internal heating, without specifying the uncertain heating mechanism(s). Our models start with a hot, high-entropy planet whose radius decreases as the convective interior cools. The applied heating suppresses this cooling. We find that very shallow heating—at pressures of 1--10 bars—does not significantly suppress cooling, unless the total heating rate is ≳10% of the incident stellar power. Deeper heating, at 100 bars, requires heating at only 1% of the stellar irradiation to explain the observed transit radius of 1.4R{sub Jup} after 5 Gyr of cooling. In general, more intense and deeper heating results in larger hot-Jupiter radii. Surprisingly, we find that heat deposited at 10{sup 4} bars—which is exterior to ≈99% of the planet’s mass—suppresses planetary cooling as effectively as heating at the center. In summary, we find that relatively shallow heating is required to explain the radii of most hot Jupiters, provided that this heat is applied early and persists throughout their evolution.

79 ASTRONOMY AND ASTROPHYSICS↗

The Physical Origin of Long Gas Depletion Times in Galaxies

We present a model that explains why galaxies form stars on a timescale significantly longer than the timescales of processes governing the evolution of interstellar gas. We show that gas evolves from a non-star-forming to a star-forming state on a relatively short timescale, and thus the rate of this evolution does not limit the star formation rate (SFR). Instead, the SFR is limited because only a small fraction of star-forming gas is converted into stars before star-forming regions are dispersed by feedback and dynamical processes. Thus, gas cycles into and out of a star-forming state multiple times, which results in a long timescale on which galaxies convert gas into stars. Our model does not rely on the assumption of equilibrium and can be used to interpret trends of depletion times with the properties of observed galaxies and the parameters of star formation and feedback recipes in simulations. In particular, the model explains how feedback self-regulates the SFR in simulations and makes it insensitive to the local star formation efficiency. We illustrate our model using the results of an isolated L {sub *}-sized galaxy simulation that reproduces the observed Kennicutt–Schmidt relation for both molecular and atomic gas. Interestingly, the relation for molecular gas is almost linear on kiloparsec scales, although a nonlinear relation is adopted in simulation cells. We discuss how a linear relation emerges from non-self-similar scaling of the gas density PDF with the average gas surface density.

79 ASTRONOMY AND ASTROPHYSICS↗