SEARCH · Engineering Papers

Results for “R”

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.

SourceResearch connectionFrom yearThrough year

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

At least 901 records · Page 50

All sources NASA NTRS DOE OSTI

DOE OSTI2022

From lattice QCD to in-medium heavy-quark interactions via deep learning

Bottomonium states are key probes for experimental studies of the quark-gluon plasma (QGP) created in high-energy nuclear collisions. Theoretical models of bottomonium productions in high-energy nuclear collisions rely on the in-medium interactions between the bottom and antibottom quarks, which can be characterized by real (V R (T, r)) and imaginary (V I (T, r)) potentials, as functions of temperature and spatial separation. Recently, the masses and thermal widths of up to 3S and 2P bottomonium states in QGP were calculated using lattice quantum chromodynamics (LQCD). Starting from these LQCD results and through a novel application of deep neural network (DNN), here, we obtain model-independent results for V R (T, r) and V I (T, r). The temperature dependence of V R (T, r) was found to be very mild between T ≈ 0 - 330 MeV. Meanwhile, VI(T, r) shows rapid increase with T and r, which is much larger than the perturbation theory based expectations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Dynamical electron correlation and the chemical bond. II. Recoupled pair bonds in the $a$ 4 Σ - states of CH and CF

We extended our studies of the effect of dynamical electron correlation on the covalent bonds in the AH and AF series (A = B–F) to the recoupled pair bonds in the excited $a$ 4 Σ - states of the CH and CF molecules. Dynamical correlation is energetically less important in the $a$ 4 Σ - states than in the corresponding X 2 Π states for both molecules, which is reflected in smaller changes in bond energies ($D$ e ). Changes in the equilibrium bond distance ($R$ e ) and vibrational frequency ($ω$ e ), on the other hand, are influenced by the changes in the slope and curvature of the dynamical electron correlation energy as a function of the internuclear distance $R$, $E$ DEC ($R$). In the CH($a$ 4 Σ - ) state, these changes are much smaller than in the CH(X 2 Π) state, but in the CF($a$ 4 Σ - ) state, they are larger, reflecting a significant difference in the shapes of $E$ DEC ($R$) curves. At large $R$, the shape of $E$ DEC ($R$) curves for covalent and recoupled pair bonds is similar although different in magnitude. For the CH($a$ 4 Σ - ) state, $E$ DEC ($R$) has a minimum at $R$ = $R$ e + 0.72 Å as the orbitals associated with the formation of the recoupled pair bond switch places. $E$ DEC (R) for the CF($a$ 4 Σ - ) state decreases continuously throughout the bound region of the potential energy curve because the dynamical correlation energy associated with the electrons in the lone pair orbitals is increasing. In conclusion, these results support our earlier conclusion that we still have much to learn about the nature of dynamical electron correlation in molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Self-interacting dark scalar spikes around black holes via relativistic Bondi accretion

For this work, we consider the spike mass density profile in a dark halo by self-consistently solving the relativistic Bondi accretion of dark matter onto a non-spining black hole of mass M. We assume that the dominant component of the dark matter in the halo is a Standard model gauge-singlet scalar. Its mass m ≃ 10 -5 eV and quartic self-coupling λ ≲ 10 19 are constrained to be compatible with the properties of galactic dark halos. In the hydrodynamic limit, we find that the accretion rate is bounded from below, M min = 96πG 2 M 2 m 4 /λ$\hbar$ 3 . Therefore, for M = 10 6 M ⊙ we have M min ≃ 1.41 x 10 -9 M ⊙ yr -1 , which is subdominant compared to the Eddington accretion of baryons. The spike density profile ρ0(r) within the self-gravitating regime cannot be fitted well by a single-power law but a double-power one. Despite that, we can fit ρ 0 (r) piecewise and find that ρ 0 (r) ∝ r -1.20 near the sound horizon, ρ 0 (r) ∝ r -1.00 towards the Bondi radius and ρ 0 (r) ∝ r -1.08 for the region in between. This contrasts with more cuspy ρ 0 (r) ∝ r -1.75 for dark matter with Coulomb-like self-interaction.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Hint of a new scalar interaction in LHCb data?

Abstract We explain recent LHCb measurements of the lepton universality ratios, $$R_{D^{(*)}}^{\tau /\ell }\equiv \frac{\mathcal {B}(\bar{B} \rightarrow D^{(*)+} \tau ^- \bar{\nu }_\tau )}{\mathcal {B}(\bar{B} \rightarrow D^{(*)+}\ell ^- \bar{\nu }_\ell )}$$ R D ( ∗ ) τ / ℓ ≡ B ( B ¯ → D ( ∗ ) + τ - ν ¯ τ ) B ( B ¯ → D ( ∗ ) + ℓ - ν ¯ ℓ ) and $${R(\Lambda _c^+)}^{\tau /\ell } \equiv \frac{\mathcal {B}(\Lambda _b \rightarrow \Lambda _c^+ \tau ^- \bar{\nu }_{\tau })}{\mathcal {B}(\Lambda _b \rightarrow \Lambda _c^+ \ell ^- \bar{\nu }_{\ell })}$$ R ( Λ c + ) τ / ℓ ≡ B ( Λ b → Λ c + τ - ν ¯ τ ) B ( Λ b → Λ c + ℓ - ν ¯ ℓ ) with $$\ell =\mu $$ ℓ = μ , via new physics that affects $$R_D^{\tau /\ell }$$ R D τ / ℓ and $$R(\Lambda _c^+)^{\tau /\ell }$$ R ( Λ c + ) τ / ℓ but not $$R_{D^*}^{\tau /\ell }$$ R D ∗ τ / ℓ . The scalar operator in the effective theory for new physics is indicated. We find that the forward-backward asymmetry and $$\tau $$ τ polarization in $$\bar{B} \rightarrow D^+ \tau ^{-} \bar{\nu }_{\tau }$$ B ¯ → D + τ - ν ¯ τ and $$\Lambda _b \rightarrow \Lambda _c^+ \tau ^- \bar{\nu }_{\tau }$$ Λ b → Λ c + τ - ν ¯ τ decays are significantly affected by the scalar interaction. We construct a simple two Higgs doublet model as a realization of our scenario and consider lepton universality in semileptonic charm and top decays, radiative B decay, B -mixing, and $$Z \rightarrow b \bar{b}$$ Z → b b ¯ .

Physics↗

DOE OSTI2022

Influences of Environmental Relative Humidity and Horizontal Scale of Subcloud Ascent on Deep Convective Initiation

This study examines two factors impacting initiation of moist deep convection: free-tropospheric environmental relative humidity (φ E ) and horizontal scale of subcloud ascent (R sub ), the latter exerting a dominant control on cumulus cloud width. A simple theoretical model is used to formulate a “scale selection” hypothesis: that a minimum R sub is required for moist convection to go deep, and that this minimum scale decreases with increasing φ E . Specifically, the ratio of R 2 sub to saturation deficit (1 - φ E ) must exceed a certain threshold value that depends on cloud-layer environmental lapse rate. Idealized, large-eddy simulations of moist convection forced by horizontally varying surface fluxes show strong sensitivity of maximum cumulus height to both φ E and R sub consistent with the hypothesis. Increasing R sub by only 300–400 m can lead to a large increase (>5 km) in cloud height. A passive tracer analysis shows that the bulk fractional entrainment rate decreases rapidly with R sub but depends little on φ E . However, buoyancy dilution increases as either R sub or φ E decreases; buoyancy above the level of free convection is rapidly depleted in dry environments when R sub is small. While deep convective initiation occurs with an increase in relative humidity of the near environment from moistening by earlier convection, the importance of this moisture preconditioning is inconclusive as it is accompanied by an increase in R sub . Overall, it is concluded that small changes to R sub driven by external forcing or by convection itself could be a dominant regulator of deep convective initiation.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Bridging Scales in Black Hole Accretion and Feedback: Relativistic Jet Linking the Horizon to the Host Galaxy

Simulating black hole (BH) accretion and feedback from the BH horizon to galactic scales is extremely challenging, as it involves a vast range of scales. Recently, our multizone method has successfully achieved global dynamical steady states of hot accretion flows in 3D general relativistic magnetohydrodynamic simulations by tracking the bidirectional interaction between a nonspinning BH and its host galaxy. In this paper, we present technical improvements to the method and apply it to spin a * = 0.9 BHs, which power relativistic jets. We first test the new multizone setup with a smaller Bondi radius, R B ≈ 400 r g , where r g is the gravitational radius. The strongly magnetized accretion launches a relativistic jet with an intermediate feedback efficiency η ∼ 30%, in between that of a prograde (η ∼ 100%) and retrograde (η ∼ 10%) torus. Interestingly, both prograde and retrograde simulations also eventually converge to the same intermediate efficiency when evolved long enough, as accumulated magnetic fields remove gas rotation. We then extend strongly magnetized simulations to larger Bondi radii, R B ≈ 2 × 10 3 , 2 × 10 4 , 2 × 10 5 r g . We find that the BH accretion rate $\dot{M}$ is suppressed with respect to the Bondi rate as $\dot{M}_{\textrm{B}}$ as $\dot{M}/\dot{M}_{\textrm{B}} ∝ R_{\textrm{B}}^{-1/2}$. However, despite some variability, the time-averaged feedback efficiency remains at η ∼ 30%, independent of R B . This suggests that BH feedback efficiency in hot accretion flows is mainly governed by the BH spin (a * ) rather than by the galactic properties (R B ). From these first-principles simulations, we provide a feedback subgrid prescription for cosmological simulations: $\dot{E}_{\textrm{fb}} = 2$ x $10^{-3}[R_{\textrm{B}}/(2$ x $10^5 r_g)]^{-1/2}$ $\dot{M}_{\textrm{B}}c^2$ for BH spin a * = 0.9.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

A ``Tail" of Two Cities: Fatality-based Modeling of COVID-19 Evolution in New York City and Cook County, IL

I describe $SIR$ modeling of the COVID-19 pandemic in two U.S. urban environments, New York City (NYC) and Cook County, IL, from onset through the month of June, 2020. Since testing was not widespread early in the pandemic in the U.S., I do not use data on confirmed cases and rely solely on public fatality data to estimate model parameters. Fits to the first 20 days of data determine a degenerate combination of the basic reproduction number, $R_0$, and the mean time to removal from the infectious population, $\gamma^{-1}$, with $\gamma (R_0-1) = 0.25 (0.21)$ inverse days for NYC (Cook County). Equivalently, the initial doubling time was $t_d = 2.8 (3.4)$ days for NYC (Cook). The early fatality data suggest that both locations had infections in early February. I model the mitigation measures implemented in mid-March in both locations (distancing, quarantine, isolation, etc) via a time-dependent reproduction number $R_t$ that declines monotonically from $R_0$ to a sma ller asymptotic value, with a parameterized functional form. The timing (mid-March) and duration (several days) of the transitions in $R_t$ appear well determined by the data. However, the fatality data determine only a degenerate combination of the parameters $R_0$, the percentage reduction in social contact due to mitigation measures, $X$, and the infection fatality rate (IFR), $f$. With flat priors, based on simulations the NYC model parameters have 95.45\% credible intervals of $R_0 = 3.0-5.4$, $X= 80-99.9\%$ and $f=2-6\%$, with $5-13\%$ of the population asymptotically infected. A strong external prior indicating a lower value of $f$ or of $\gamma^{-1}$ would imply lower values of $R_0$ and $X$ and higher percentage infection of the population. For Cook County, the evolution was qualitatively different: after mitigation measures were implemented, the daily fatality counts reached a plateau for about a month before tailing off. This is consistent with an $SIR$ model that exhibits ``critical slowing-down", in which $R_t$ plateaus ! at a value just above unity. For Cook County, the 95.45\% credible intervals for the model parameters are much broader and shifted downward, $R_0=1.4-4.7$, $X=26-54\%$, and $f=0.1-0.6\%$ with $15-88\%$ of the population asymptotically infected. Despite the apparently lower efficacy of its social contact reduction measures, Cook County has had significantly fewer fatalities per population than NYC, $D_\infty/N=100$ vs. 270 per 100,000. In the model, this is attributed to the lower inferred IFR for Cook; an external prior pointing to similar values of the IFR for the two locations would instead chalk up the difference in $D/N$ to differences in the relative growth rate of the disease. I derive a model-dependent threshold, $X^{crit}_e$, for `safe' re-opening, that is, for easing of contact reduction that would not trigger a second wave; for NYC, the models predict that increasing social contact by more than 20\% from post-mitigation levels will lead to renewed spread, while for Cook County the threshold value is very uncertain, given the parameter degeneracies. The timing of 2nd-wave growth will depend on the amplitude of contact increase relative to $X^{crit}_e$ and on the asymptotic growth rate, and the impact in terms of fatalities will depend on the parameter $f$.

Frieman, Joshua↗

DOE OSTI2024

Measurements of $$ {K}_S^0 $$-$$ {K}_L^0 $$ asymmetries in the decays $$ {\Lambda}_c^{+}\to p{K}_{L,S}^0 $$, $$ p{K}_{L,S}^0{\pi}^{+}{\pi}^{-} $$ and $$ p{K}_{L,S}^0{\pi}^0 $$

Abstract Usinge + e − annihilation data sets corresponding to an integrated luminosity of 4.5 fb −1 , collected with the BESIII detector at center-of-mass energies between 4.600 and 4.699 GeV, we report the first measurements of the absolute branching fractions$$ \mathcal{B}\left({\Lambda}_c^{+}\to p{K}_L^0\right) $$ B Λ c + → p K L 0 = (1.67±0.06±0.04)%,$$ \mathcal{B}\left({\Lambda}_c^{+}\to p{K}_L^0{\pi}^{+}{\pi}^{-}\right) $$ B Λ c + → p K L 0 π + π − = (1.69±0.10±0.05)%, and$$ \mathcal{B}\left({\Lambda}_c^{+}\to p{K}_L^0{\pi}^0\right) $$ B Λ c + → p K L 0 π 0 = (2.02±0.13±0.05)%, where the first uncertainties are statistical and the second systematic. Combining with the known branching fractions of$$ {\Lambda}_c^{+}\to p{K}_S^0 $$ Λ c + → p K S 0 ,$$ {\Lambda}_c^{+}\to p{K}_S^0{\pi}^{+}{\pi}^{-} $$ Λ c + → p K S 0 π + π − , and$$ {\Lambda}_c^{+}\to p{K}_S^0{\pi}^0 $$ Λ c + → p K S 0 π 0 , we present the first measurements of the$$ {K}_S^0 $$ K S 0 -$$ {K}_L^0 $$ K L 0 asymmetries$$ R\left({\Lambda}_c^{+},{K}_{S,L}^0X\right)=\frac{\mathcal{B}\left({\Lambda}_c^{+}\to {K}_S^0X\right)-\mathcal{B}\left({\Lambda}_c^{+}\to {K}_L^0X\right)}{\mathcal{B}\left({\Lambda}_c^{+}\to {K}_S^0X\right)+\mathcal{B}\left({\Lambda}_c^{+}\to {K}_L^0X\right)} $$ R Λ c + K S , L 0 X = B Λ c + → K S 0 X − B Λ c + → K L 0 X B Λ c + → K S 0 X + B Λ c + → K L 0 X in charmed baryon decays:$$ R\left({\Lambda}_c^{+},p{K}_{S,L}^0\right)=-0.025\pm 0.031 $$ R Λ c + p K S , L 0 = − 0.025 ± 0.031 ,$$ R\left({\Lambda}_c^{+},p{K}_{S,L}^0{\pi}^{+}{\pi}^{-}\right)=-0.027\pm 0.048 $$ R Λ c + p K S , L 0 π + π − = − 0.027 ± 0.048 and$$ R\left({\Lambda}_c^{+},p{K}_{S,L}^0{\pi}^0\right)=-0.015\pm 0.046 $$ R Λ c + p K S , L 0 π 0 = − 0.015 ± 0.046 . No significant asymmetries with statistical significance are observed.

Physics↗

DOE OSTI2022

Evidence of galaxy assembly bias in SDSS DR7 galaxy samples from count statistics

We present observational constraints on the galaxy–halo connection, focusing particularly on galaxy assembly bias from a novel combination of counts-in-cylinders statistics, P(N CIC ), with the standard measurements of the projected two-point correlation function w p (r p ), and number density n gal of galaxies. We measure n gal , w p (r p ), and P(N CIC ) for volume-limited, luminosity-threshold samples of galaxies selected from SDSS DR7, and use them to constrain halo occupation distribution (HOD) models, including a model in which galaxy occupation depends upon a secondary halo property, namely halo concentration. We detect significant positive central assembly bias for the M r < -20.0 and M r <-19.5 samples. Central galaxies preferentially reside within haloes of high concentration at fixed mass. Positive central assembly bias is also favoured in the M r < -20.5 and M r < -19.0 samples. We find no evidence of central assembly bias in the M r < -21.0 sample. We observe only a marginal preference for negative satellite assembly bias in the M r < -20.0 and M r < -19.0 samples, and non-zero satellite assembly bias is not indicated in other samples. Our findings underscore the necessity of accounting for galaxy assembly bias when interpreting galaxy survey data, and demonstrate the potential of count statistics in extracting information from the spatial distribution of galaxies, which could be applied to both galaxy–halo connection studies and cosmological analyses.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Measurement of the nuclear modification factor of $b$-jets in 5.02 TeV Pb+Pb collisions with the ATLAS detector

This paper presents a measurement of b-jet production in Pb+Pb and pp collisions at $\sqrt{sNN}$ = 5.02 TeV with the ATLAS detector at the LHC. The measurement uses 260 pb –1 of pp collisions collected in 2017 and 1.4 nb –1 of Pb+Pb collisions collected in 2018. In both collision systems, jets are reconstructed via the anti-$k_t$ algorithm. The b-jets are identified from a sample of jets containing muons from the semileptonic decay of b-quarks using template fits of the muon momentum relative to the jet axis. In pp collisions, b-jets are reconstructed for radius parameters R = 0.2 and R = 0.4, and only R = 0.2 jets are used in Pb+Pb collisions. For comparison, inclusive R = 0.2 jets are also measured using 1.7 nb –1 of Pb+Pb collisions collected in 2018 and the same pp collision data as the b-jet measurement. The nuclear modification factor, $R_{AA}$, is calculated for both b-jets and inclusive jets with R = 0.2 over the transverse momentum range of 80–290 GeV. The nuclear modification factor for b-jets decreases from peripheral to central collisions. The ratio of the b-jet $R_{AA}$ to inclusive jet $R_{AA}$ is also presented and suggests that the $R_{AA}$ for b-jets is larger than that for inclusive jets in central Pb+Pb collisions. The measurements are compared with theoretical calculations and suggest a role for mass and colour-charge effects in partonic energy loss in heavy-ion collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Evidence of Mass Ordering of Charm and Bottom Quark Energy Loss in Au+Au Collisions at RHIC

Abstract Partons traversing the strongly interacting medium produced in heavy-ion collisions are expected to lose energy depending on their color charge and mass. We measure the nuclear modification factors for charm- and bottom-decay electrons, defined as the ratio of yields, divided by the number of binary nucleon–nucleon collisions, in $$\sqrt{s_{\textrm{NN}}}=200$$ s NN = 200 GeV Au+Au collisions to p + p collisions ( $$R_{\textrm{AA}}$$ R AA ), or in central to peripheral Au+Au collisions ( $$R_{\textrm{CP}}$$ R CP ). We find the bottom-decay electron $$R_{\textrm{AA}}$$ R AA and $$R_{\textrm{CP}}$$ R CP to be significantly higher than those of charm-decay electrons. Model calculations including mass-dependent parton energy loss in a strongly coupled medium are consistent with the measured data. These observations provide evidence of mass ordering of charm and bottom quark energy loss when traversing through the strongly coupled medium created in heavy-ion collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

Sterile neutrino dark matter in left-right theories

SU(2) L × SU(2) R gauge symmetry requires three right-handed neutrinos (N i ), one of which, N 1 , can be sufficiently stable to be dark matter. In the early universe, W R exchange with the Standard Model thermal bath keeps the right-handed neutrinos in thermal equilibrium at high temperatures. N 1 can make up all of dark matter if they freeze-out while relativistic and are mildly diluted by subsequent decays of a long-lived and heavier right-handed neutrino, N 2 . We systematically study this parameter space, constraining the symmetry breaking scale of SU(2) R and the mass of N 1 to a triangle in the (v R , M 1 ) plane, with v R = (10 6 - 3 × 10 12 ) GeV and M 1 = (2 keV–1 MeV). Much of this triangle can be probed by signals of warm dark matter, especially if leptogenesis from N 2 decay yields the observed baryon asymmetry. The minimal value of v R is increased to 10 8 GeV for doublet breaking of SU(2) R , and further to 10 9 GeV if leptogenesis occurs via N 2 decay, while the upper bound on M 1 is reduced to 100 keV. In addition, there is a component of hot N 1 dark matter resulting from the late decay of N 2 → N 1 ℓ+ℓ- that can be probed by future cosmic microwave background observations. Interestingly, the range of v R allows both precision gauge coupling unification and the Higgs Parity understanding of the vanishing of the Standard Model Higgs quartic at scale vR. Finally, we study freeze-in production of N 1 dark matter via the WR interaction, which allows a much wider range of (v R , M 1 ).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Soil organic carbon accrual due to more efficient microbial utilization of plant inputs at greater long-term soil moisture

High long-term soil moisture may either stimulate or inhibit soil organic carbon (SOC) losses through changes to mineral and chemical composition, and resultant organo-mineral interactions. Yet, the trade-off between mineralization and accrual of SOC under long-term variation in unsaturated soil moisture remains uncertain. We tested the underexplored relationships between long-term soil moisture and organo-mineral chemical composition and its implications for SOC persistence in an experimental field in New York, USA, with differences in long-term mean soil volumetric water content (0–0.15 m depth) ranging from 0.40 to 0.63 (v/v) during the growing season. Long-term soil moisture across 20 subplots on four fallow plots were positively correlated with SOC (R 2 = 0.23; P = 0.019, n = 20), mineral-associated organic matter (MAOM) content (g fraction/g soil) (R 2 = 0.44; P = 0.001; n = 20) and occluded particulate organic matter (oPOM) content (R 2 = 0.18; P = 0.033; n = 20). Higher long-term soil moisture was associated with a decrease in the relative content of sodium pyrophosphate extractable Fe (R 2 = 0.33; P < 0.005; n = 20), an increase in sodium dithionite extractable Fe (R 2 = 0.44; P < 0.001; n = 20), and an increase in SOC retention by non-crystalline Al pools (R 2 = 0.51; P = 0.0002 for sodium pyrophosphate extracts, R 2 = 0.41; P = 0.0014 for hydroxylamine hydrochloride extracts; n = 20 for both). Increasing long-term soil moisture was associated with a four-fold increase in microbial biomass C (per unit SOC) and lower metabolic quotient (R 2 = 0.56, P < 0.001). MAOM fractions of high-moisture soils had lower C:N (from C:N 9.5 to 9.0, R 2 = 0.27, P = 0.011, n = 20). Consistent with decreasing C:N, increasing decomposition with increasing moisture was reflected by a 15% and 10% greater proportion of oxidized carboxylic-C to aromatic-C and O-alkyl C, respectively, as measured with 13 C NMR, and a more pronounced FTIR signature of N-containing proteinaceous compounds in high-moisture MAOM fractions, indicative of microbial metabolites and transformation products. A partial least squares regression showed that SOC content increased with greater long-term moisture (P = 0.019), pyrophosphate-extractable Al (P = 0.0001), and exchangeable Ca (P = 0.013). Here our results show that higher long-term soil moisture resulted in SOC accrual by enhancing conversion of plant inputs into microbial biomass that interacts with reactive minerals.

58 GEOSCIENCES↗

DOE OSTI2021

Airborne hyperspectral imaging of nitrogen deficiency on crop traits and yield of maize by machine learning and radiative transfer modeling

Nitrogen is an essential nutrient that directly affects plant photosynthesis, crop yield, and biomass production for bioenergy crops, but excessive application of nitrogen fertilizers can cause environmental degradation. To achieve sustainable nitrogen fertilizer management for precision agriculture, there is an urgent need for nondestructive and high spatial resolution monitoring of crop nitrogen and its allocation to photosynthetic proteins as that changes over time. Here, we used visible to shortwave infrared (400–2400 nm) airborne hyperspectral imaging with high spatial (0.5 m) and spectral (3–5 nm) resolutions to accurately estimate critical crop traits, i.e., nitrogen, chlorophyll, and photosynthetic capacity (CO 2 -saturated photosynthesis rate, V max,27 ), at leaf and canopy scales, and to assess nitrogen deficiency on crop yield. We conducted three airborne campaigns over a maize (Zea mays L.) field during the growing season of 2019. Physically based soil-canopy Radiative Transfer Modeling (RTM) and data-driven approaches i.e. Partial-Least Squares Regression (PLSR) were used to retrieve crop traits from hyperspectral reflectance, with ground truth of leaf nitrogen, chlorophyll, V max,27 , Leaf Area Index (LAI), and harvested grain yield. To improve computational efficiency of RTMs, Random Forest (RF) was used to mimic RTM simulations to generate machine learning surrogate models RTM-RF. The results show that prior knowledge of soil background and leaf angle distribution can significantly reduce the ill-posed RTM retrieval. RTM-RF achieved a high accuracy to predict leaf chlorophyll content (R 2 = 0.73) and LAI (R 2 = 0.75). Meanwhile, PLSR exhibited better accuracy to predict leaf chlorophyll content (R 2 = 0.79), nitrogen concentration (R 2 = 0.83), nitrogen content (R 2 = 0.77), and V max,27 (R 2 = 0.69) but required measured traits for model training. We also found that canopy structure signals can enhance the use of spectral data to predict nitrogen related photosynthetic traits, as combining RTM-RF LAI and PLSR leaf traits well predicted canopy-level traits (leaf traits × LAI) including canopy chlorophyll (R 2 = 0.80), nitrogen (R 2 = 0.85) and V max,27 (R 2 = 0.82). Compared to leaf traits, we further found that canopy-level photosynthetic traits, particularly canopy V max,27 , have higher correlation with maize grain yield. This study highlights the potential for synergistic use of process-based and data-driven approaches of hyperspectral imaging to quantify crop traits that facilitate precision agricultural management to secure food and bioenergy production.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Revealing Local and Directional Aspects of Catalytic Active Sites by the Nuclear and Surface Electrostatic Potential

This work examines the prospects of using the electrostatic potential, V(r), as a descriptor in heterogeneous catalysis. In particular, the subatomic spatial resolution of the property allows for analysis of both directionality and confinement effects in surface adsorption. This feature of V(r) is used to identify adsorption sites, orientations, and energetics for metal surfaces, particles, and nanoclusters upon interactions with catalytically relevant intermediates. The use of V(r) in assessing the 3D nature of catalytic sites in low-temperature and electrocatalysis is highlighted, and future directions in catalysis design are discussed. Ultimately, we provide a critical analysis of the use of V(r) in the predictions of local adsorption susceptibilities, and we address its limitations. The link between V(r) and other established descriptors in catalysis are motivated via physical relations and theoretical derivations; close ties are established between V(r) and the d-band center (ε d ), as well as the surface site stability (BE M ). In conclusion, by comparing the performance of V(r) evaluated on isodensity contours, i.e., the surface electrostatic potential, to that of V(r) evaluated at the nucleus of an atom, we investigate the application space for a directional and an atom-localized version of the V(r) descriptor for catalyst design.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Synthesis and Characterization of Uranium Complexes Supported by Substituted Aryldimethylsilylanilide Ligands

Here, we report the synthesis and characterization of substituted aryldimethylsilyldiisopropylanilide ligands and their respective bisamido complexes of U(III), (3,5-R 2 -PhMe 2 SiNDipp) 2 UI(dioxane) x (1, R = H, x= 0; 2, R = Me, x = 0; 3, R = t Bu, x = 1). We found that the steric bulk of the 3,5-R 2 -Ph ring affects the hapticity of the U–arene interaction. In the solid-state, 1 is a U–(η 6 -arene) complex, while 2 is a bis(U–(η 1 -arene)) complex. Theoretically calculated bond orders at PBE0 and PBE0-D3 levels of theory support these hapticity assignments. The 3,5- t Bu 2 -Ph rings of 3 are too bulky to interact with U and solid-state metrical parameters initially suggested a U–(η 1 -arene) interaction with one of the Dipp rings. However, bond order calculations show that this interaction is even weaker than in the previously reported ((PhMe 2 Si) 2 N) 3 U complex, leading to the conclusion that 3 is best described as a U–(η o -arene) complex. Molecular orbital analyses in conjunction with electron localization methods reveal that the U–(η 6 -arene) bonding in 1 is primarily electrostatic in nature. Some charge transfer takes place from the arene π orbitals to the U 6d/5f hybrid orbitals in addition to subtle δ-back-bonding. In 2 and 3, both π and δ interactions are substantially weaker, in agreement with the differences in the U–arene coordination modes. Surprisingly, attempts to generate less sterically bulky (3,5-R 2 -PhMe 2 SiNPh) 2 UI complexes results in disproportionation to homoleptic tetraamido (3,5-R 2 -PhMe 2 SiNPh) 4 U(IV) (4, R = H; 5, R = Me) complexes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Stabilization of reactive rare earth alkyl complexes through mechanistic studies

Rare earth tris(alkyl) complexes such as M(CH 2 SiMe 3 ) 3 (sol) n are widely used as precursors for many compounds and as homogeneous catalysts for alkene polymerization and alkane functionalization. However, the thermal instability of those most conveniently made from the commercially available lithium salt of the neosilyl anion, LiCH 2 SiMe 3 , Li(r), restricts their utility. We present a new range of synthetically useful, more kinetically stable rare earth neosilyl solvates, derived from a full kinetic study of the various possible decomposition mechanisms of 7 known and 12 new solvated rare earth neosilyl complexes M(CH 2 SiMe 3 ) 3 (sol) n M = Sc(III), Y(III), Lu(III), Sm(III), and sol = THF; TMEDA; DMPE; diglyme ((CH 3 ) 2 (OCH 2 CH 2 ) 2 O, G 2 ), triglyme ((CH 3 ) 2 (OCH 2 CH 2 ) 3 O, G 3 ). Surprisingly, simply using higher-denticity donors to sterically disfavor neosilyl γ-H elimination is not effective. While Sc(r) 3 ((CH 3 ) 2 (OCH 2 CH 2 ) 2 O) has a half-life, t 1/2 , of 258.1 h, six times longer than for Sc(r) 3 (C 4 H 8 O) 2 (t 1/2 = 43 h), Lu(r) 3 ((CH 3 ) 2 (OCH 2 CH 2 ) 2 O) and Y(r) 3 ((CH 3 ) 2 (OCH 2 CH 2 ) 2 O) do not show the expected, analogous increased t 1/2 . This is because new decomposition pathways appear for poorly fitting donors. Finally, kinetic studies demonstrate the impact of small, and increasing amounts of LiCl on the kinetics of the reactivity of the smaller alkyls Y(r) 3 (THF) 2 and Lu(r) 3 (THF) 2 ; molecules used in hydrocarbon chemistry and catalysis for fifty years. A new route to pure Y(r) 3 (THF) 2 , which avoids the traditional use of Li(r), is presented.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

The metabolic origins of non-photorespiratory CO 2 release during photosynthesis: a metabolic flux analysis

Respiration in the light (R L ) releases CO 2 in photosynthesizing leaves and is a phenomenon that occurs independently from photorespiration. Since R L lowers net carbon fixation, understanding R L could help improve plant carbon-use efficiency and models of crop photosynthesis. Although R L was identified more than 75 years ago, its biochemical mechanisms remain unclear. To identify reactions contributing to R L , we mapped metabolic fluxes in photosynthesizing source leaves of the oilseed crop and model plant camelina (Camelina sativa). We performed a flux analysis using isotopic labeling patterns of central metabolites during 13 CO 2 labeling time course, gas exchange, and carbohydrate production rate experiments. To quantify the contributions of multiple potential CO 2 sources with statistical and biological confidence, we increased the number of metabolites measured and reduced biological and technical heterogeneity by using single mature source leaves and quickly quenching metabolism by directly injecting liquid N 2 ; we then compared the goodness-of-fit between these data and data from models with alternative metabolic network structures and constraints. Our analysis predicted that R L releases 5.2 μmol CO2 g –1 FW h –1 of CO 2 , which is relatively consistent with a value of 9.3 μmol CO 2 g –1 FW h –1 measured by CO 2 gas exchange. The results indicated that ≤10% of R L results from TCA cycle reactions, which are widely considered to dominate R L . Further analysis of the results indicated that oxidation of glucose-6-phosphate to pentose phosphate via 6-phosphogluconate (the G6P/OPP shunt) can account for >% of CO 2 released by R L .

59 BASIC BIOLOGICAL SCIENCES↗