Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “favorability”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Computational Predictions of the Hydrolysis of 2,4,6-Trinitrotoluene (TNT) and 2,4-Dinitroanisole (DNAN)

Hydrolysis is a common transformation reaction that can affect the environmental fate of many organic compounds. In this study, three proposed mechanisms of alkaline hydrolysis of 2,4,6-trinitrotoluene (TNT) and 2,4-dinitroaniline (DNAN) were investigated with plane-wave density functional theory (DFT) combined with ab initio and classical molecular dynamics (AIMD/MM) free energy simulations, Gaussian basis set DFT calculations, and correlated molecular orbital theory calculations. Most of the computations in this study were carried out using the Arrows web based tools. For each mechanism: Meisenheimer complex formation, nucleophilic aromatic substitution, and proton abstraction reaction energies and activation barriers were calculated for the reaction at each relevant site. For TNT, it was found that the most kinetically favorable first hydrolysis steps involve Meisenheimer complex formation by attachment of OH$^{-}$ at the C1 and C3 arene carbons and proton abstraction from the methyl group. The nucleophilic aromatic substitution reactions at the C2 and C4 arene carbons were found to be thermodynamically favorable. However, the calculated activation barriers were slightly lower than previous studies, but still found to be $\Delta$G$^{\ddagger} \approx 18$ kcal/mol using PBEo AIMD/MM free energy simulations, suggesting the reactions are not kinetically significant. For DNAN, the barriers of nucleophilic aromatic substitution were even greater ($\Delta$G$^{\ddagger} > 29$ kcal/mol PBEo AIMD/MM). Further, the most favorable hydrolysis reaction for DNAN was found to be a two-step process in which the hydroxyl first attacks the C1 carbon to form a Meisenheimer complex at the C1 arene carbon C1-(OCH$_3$)OH$^{-}$, subsequently, the methoxy anion (-OCH$_3$) at the C1 arene carbon dissociates and the proton shuttles from the C1-OH to the dissociated methoxy group, resulting in methanol and an aryloxy anion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Quantum Sieving for Isotopic Separations of Gases Using Porous Materials─30 Years of Progress

Here, this paper reviews theoretical and experimental efforts to establish microporous materials that can be used to separate isotopologues of molecular gases such as H 2 and D 2 . Emphasis is placed on use of simplified models to highlight the quantum phenomena that make these separations possible. In equilibrium adsorption, differences in zero point energy in the adsorbed states of molecules favor binding of heaver isotopologues (e.g. D 2 relative to H 2 ). Experimental data showing this effect was reported as early as 1933, but the theoretical work of Beenakkers et al. in 1995 [Chem. Phys. Lett. 232 (1995) 379-382] spurred modern efforts to develop materials with strong so-called quantum sieving. In porous materials where the transition state for molecular diffusion is more strongly confined than the energy minima associated with equilibrium adsorption, differences in zero point energies between these two sites can lead to isotopic differences in molecular diffusivities. This effect favors diffusion of heavier species (e.g. D 2 ) relative to lighter species (e.g. H 2 ). This so-called kinetic quantum sieving has been observed experimentally in porous carbons and in porous organic cage materials. We show that quantum tunneling, which favors hopping of lighter species across energy barriers, diminishes the strength of kinetic quantum sieving but that it appears to make only a small contribution to the net molecular diffusivities in many porous materials of practical interest.

Sholl, David S. [Oak Ridge National Laboratory (OR↗
DOE OSTI2024

N and OH-Immobilized Cu 3 Clusters In Situ Reconstructed from Single-Metal Sites for Efficient CO 2 Electromethanation in Bicontinuous Mesochannels

Cu-based catalysts hold promise for electrifying CO 2 to produce methane, an extensively used fuel. However, the activity and selectivity remain insufficient due to the lack of catalyst design principles to steer complex CO 2 reduction pathways. Herein, we develop a concept to design carbon-supported Cu catalysts by regulating Cu active sites’ atomic-scale structures and engineering the carbon support’s mesoscale architecture. This aims to provide a favorable local reaction microenvironment for a selective CO 2 reduction pathway to methane. In situ X-ray absorption and Raman spectroscopy analyses reveal the dynamic reconstruction of nitrogen and hydroxyl-immobilized Cu 3 (N,OH-Cu 3 ) clusters derived from atomically dispersed Cu–N 3 sites under realistic CO 2 reduction conditions. The N,OH-Cu 3 sites possess moderate *CO adsorption affinity and a low barrier for *CO hydrogenation, enabling intrinsically selective CO 2 -to-CH 4 reduction compared to the C–C coupling with a high energy barrier. Importantly, a block copolymer-derived carbon fiber support with interconnected mesopores is constructed. The unique long-range mesochannels offer an H 2 O-deficient microenvironment and prolong the transport path for the CO intermediate, which could suppress the hydrogen evolution reaction and favor deep CO 2 reduction toward methane formation. Thus, the newly developed catalyst consisting of in situ constructed N,OH-Cu 3 active sites embedded into bicontinuous carbon mesochannels achieved an unprecedented Faradaic efficiency of 74.2% for the CO 2 reduction to methane at an industry-level current density of 300 mA cm –2 . In conclusion, this work explores effective concepts for steering desirable reaction pathways in complex interfacial catalytic systems via modulating active site structures at the atomic level and engineering pore architectures of supports on the mesoscale to create favorable microenvironments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Divergent Bimetallic Mechanisms in Copper(II)-Mediated C–C, N–N, and O–O Oxidative Coupling Reactions

Copper-catalyzed aerobic oxidative coupling of diaryl imines provides a route for conversion of ammonia to hydrazine. Here the present study uses experimental and density functional theory computational methods to investigate the mechanism of N–N bond formation, and the data support a mechanism involving bimolecular coupling of Cu-coordinated iminyl radicals. Computational analysis is extended to Cu II -mediated C–C, N–N, and O–O coupling reactions involved in the formation of cyanogen (NC–CN) from HCN, 1,3-butadiyne from ethyne (i.e., Glaser coupling), hydrazine from ammonia, and hydrogen peroxide from water. The results reveal two different mechanistic pathways. Heteroatom ligands with an uncoordinated lone pair (iminyl, NH 2 , OH) undergo charge transfer to Cu II , generating ligand-centered radicals that undergo facile bimolecular radical–radical coupling. Ligands lacking a lone pair (CN and CCH) form bridged binuclear diamond-core structures that undergo C–C coupling. This mechanistic bifurcation is rationalized by analysis of spin densities in key intermediates and transition states, as well as multiconfigurational calculations. Radical–radical coupling is especially favorable for N–N coupling owing to energetically favorable charge transfer in the intermediate and thermodynamically favorable product formation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Alkane Dehydrogenation and H/D Exchange by a Cationic Pincer-Ir(III) Hydride: Cooperative C–H Addition and β -H Elimination Modes Induce Anomalous Selectivity

We report that the cationic iridium complex ( iPr PCP)IrH + catalyzes the transfer-dehydrogenation of alkanes to give alkenes and hydrogen isotope exchange (HIE) of alkanes and arenes. Contrary to established selectivity trends found for C–H activation by transition metal complexes, strained cycloalkanes, including cyclopentane, cycloheptane, and cyclooctane, undergo C–H addition much more readily than n-alkanes, which in turn are much more reactive than cyclohexane. Aromatic C–H bonds also undergo H/D exchange much less rapidly than those of the strained cycloalkanes, but much more favorably than cyclohexane. The order of reactivity toward dehydrogenation correlates qualitatively with the reaction thermodynamics, but the magnitude is much greater than can be explained by thermodynamics. Accordingly, the cycloalkenes corresponding to the strained cycloalkanes undergo hydrogenation much more readily than cyclohexene, despite the less favorable thermodynamics of such hydrogenations. Here, computational (DFT) studies allow rationalization of the origin of reactivity and the unusual selectivity. Specifically, the initial C–H addition is strongly assisted by β-agostic interactions, which are particularly favorable for the strained cycloalkanes. Subsequent to α-C–H addition, the H atom of the β-agostic C–H bond is transferred directly to the hydride ligand of ( iPr PCP)IrH + to give a dihydrogen ligand.

02 PETROLEUM↗
DOE OSTI2021

Crucial roles of eastward propagating environments in the summer MCS initiation over the U.S. Great Plains

Mesoscale convective systems (MCSs) frequently occur over the U.S. Great Plains during summer. An analysis using self-organizing map is conducted to identify four types of summer MCS initiation environments during 2004-2017. The first two types feature favorable large-scale environments at both upper and low levels, while Type-3 has favorable lower-level and surface conditions but unfavorable upper-level circulation, and Type-4 features the most unfavorable large-scale environments for MCS initiation. Despite the unfavorable large-scale environment, the convection-centered environments in Type-4 are favorable for MCS initiation and similar to the first two types, suggesting a role of sub-synoptic disturbances as an MCS precursor. All four types of the MCS initiation delineate a clear eastward propagating feature in many fields, such as upper-level potential vorticity/geopotential height, surface pressure and surface equivalent potential temperature, upstream up to 25°-longitude west of and ~36 hours before the MCS initiation. The propagating environments and local, non-propagating low-level moisture are found to be important in MCS initiation at the foothill of the Rocky Mountains, but over the central Great Plains, it is the coupling of dynamical and moisture anomalies associated with propagating waves that results in the MCS initiation. By tracking MCSs and mid-tropospheric perturbations (MPs), a type of sub-synoptic disturbances with Rocky Mountains origin, ~30% of MPs is associated with MCS initiation, mostly in Type-4. Although MPs are related to a small fraction of MCS initiation, MCSs that are associated with MPs tend to produce more rainfall in a larger area with a stronger convective intensity, suggesting MPs to be a source of intense MCSs in summer.

Song, Fengfei↗
DOE OSTI2023

Extension of Large Fire Emissions From Summer to Autumn and Its Drivers in the Western US

Abstract Burned areas in the western US have increased ten‐fold since 1980s, which are attributable to multiple factors, including increasing heat, changing precipitation patterns, and extended drought. To better understand how these factors contribute to large fire emissions (gridded monthly fire emissions >95th percentile of all the fire emissions in the western US; 0.009 Gg/month), we build a machine learning model to predict fire emissions (PM 2.5 ) over the western US at 0.25° resolution, interpreted using explainable artificial intelligence (XAI). From the predictor contributions derived from XAI, we conduct k‐means clustering analysis to identify four clusters of predictor variables representing different drivers of large fire emissions. The four clusters feature the contributions of fuel load (Cluster 1) and different levels of dryness (Cluster 2–4), controlled by fuel moisture, drought condition, and fire‐favorable large‐scale meteorological patterns featuring high temperature, high pressure, and low relative humidity. In the past two decades, large fire emissions peak in summer. However, large fire emissions increased significantly in September and October in 2010–2020 relative to 2000–2009, extending the peak large fire emissions from summer to autumn. The larger enhancements of large fire emissions during autumn compared to summer are contributed by decreased fuel moisture, along with more frequent concurrent fire‐favorable large‐scale meteorological patterns and drought. These results highlight fuel drying as a common driver supported by multiple drivers, such as warmer temperature and more frequent synoptic patterns favorable for fires, in increasing the autumn risk of large fire emissions across the western US.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Towards cost-competitive middle distillate fuels from ethanol within a market-flexible C2 platform-based biorefinery concept

Ethanol to middle distillates (ETMD) is a promising pathway to produce sustainable liquid fuels to decarbonize the hard-to-electrify transportation sectors due to (1) the abundant sugar/starch and lignocellulosic biomass, (2) the existing deployment scale of fuel ethanol production (similar to 29 billion gallons per year globally), and (3) emerging opportunities in C2+ alcohol synthesis from CO2. Here we report a conceptual market-responsive biorefinery centered around a new ETMD pathway based on one-step ethanol to butene-rich olefins (ETO) over a Cu-Zn-Y/Beta catalyst. Specifically, this ethanol conversion pathway comprises one-step ETO, oligomerization, and hydrotreating. This ETO is distinct from that in the conventional ethanol-to-jet process which is based on two-step ethanol to ethylene and ethylene oligomerization to butenes. Butene-rich olefins can be shifted to butadiene-rich products by simply changing the reaction atmosphere from hydrogen to inert gas over the same ETO catalyst. Leveraging the experimental results, baseline techno-economic analysis (TEA) and sensitivity analysis indicate that the ethanol conversion cost is $0.60 per gallon gasoline equivalent (GGE), with opportunities for further cost reduction via improving the liquid hydrocarbon yield and space velocities, and process optimization on balancing dewatering of ethanol feed prior to the ETO step. The minimum fuel selling price (MFSP) of liquid hydrocarbons derived from corn starch ethanol with butadiene as coproduct is $1.64 per GGE, in the range that is cost competitive with petroleum kerosene-type jet fuel. Projected MFSP for cellulosic ethanol (corn stover) derived hydrocarbons is below $3.00 per GGE and co-production of butadiene further reduces the MFSP to $1.70 per GGE. The Well-to-Wake life-cycle analysis indicates that 85% greenhouse gas emission reduction can be achieved when using corn stover compared to petroleum reference and the associated carbon credits will provide significant economic incentives to favor the cellulosic ethanol-derived hydrocarbon fuels. This study demonstrates a low-cost pathway to middle distillate fuels leveraging existing ethanol infrastructure, where catalysis innovation drives the reduction of process complexity and flexible coproduction of a value-added chemical product.

Zhang, Junyan↗
DOE OSTI2026

In situ catalyst activation and regeneration enable energy-efficient high-current CO 2 reduction to ethanol-rich C 2+ mixtures

Electrochemical conversion of dissolved CO 2 in bicarbonate electrolytes, i.e., bicarbonate electrolysis, offers distinct advantages over gas diffusion electrode systems by enabling direct utilization of the CO 2 capture electrolyte while bypassing the energy-intensive CO 2 release step. However, bicarbonate electrolysis faces challenges such as CO 2 mass-transfer limitation, local pH-driven CO 2 depletion, and high cathodic potentials. The higher potential often causes catalyst surface reorganization, leading to a gradual loss of active sites and variations in selectivity during CO 2 reduction. Here, we report a directed, in situ activation and regeneration method that allows precatalysts to equilibrate under dynamic (pulsed) electrolysis conditions. We demonstrate in situ activation of a scalable Cu 2 O/Cu mesh that, under short-width (t = 4 s) pulsed electrolysis, provides stable mixed oxidation states of Cu, favoring the formation of an ethanol-rich crude mixture. The pulsed electrolysis waveform, consisting of six distinct segments, is tuned to form Cu + oxides, which are then reduced to generate local alkaline conditions favoring C–C coupling. This synergistic effect results in FEs of 73% for C2+ products and 39% for ethanol at an applied current density of −150 mA cm −2 and a cathodic potential of −1.45 V (vs. RHE). The overall half-cell energy efficiency is ∼30% for C 2+ products. The in situ Raman experiments confirm the role of pCO 2 R in dynamically regenerating Cu+-containing surface species during pulsed operation, thereby steering selectivity towards C 2+ products. A comprehensive multiscale, multiphysics model is developed to investigate the dynamic behavior of copper surface species (Cu, Cu + , and Cu 2+ ) and local microenvironmental conditions during the pCO 2 R. The results reveal that the coexistence of different copper oxidation states, especially the Cu+ intermediate, is critical in steering selectivity towards multicarbon (C 2+ ) products. The dynamic modulation of surface redox states via tailored pulsing strategies favors C–C coupling pathways by inducing localized alkaline conditions and stabilizing reactive intermediates. This work establishes a predictive modeling platform that links pulse waveform design with mechanistic insights into catalyst state evolution and product selectivity. Overall, this study provides valuable insights into the synergistic effect of in situ activation of pre-catalysts and pulsed electrolysis for higher selectivity towards C 2+ products.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Experimental study of the edge radial electric field in different drift configurations and its role in the access to H-mode at ASDEX Upgrade

The formation of the equilibrium radial electric field (Er) has been studied experimentally at ASDEX Upgrade (AUG) in L-modes of “favorable” (ion ∇ B-drift toward primary X-point) and “unfavorable” (ion ∇ B-drift away from primary X-point) drift configurations, in view of its impact on H-mode access, which changes with drift configurations. Edge electron and ion kinetic profiles and impurity velocity and mean-field Er profiles across the separatrix are investigated, employing new and improved measurement techniques. The experimental results are compared to local neoclassical theory as well as to a simple 1D scrape-off layer (SOL) model. It is found that in L-modes of matched heating power and plasma density, the upstream SOL Er and the main ion pressure gradient in the plasma edge are the same for either drift configurations, whereas the Er well in the confined plasma is shallower in unfavorable compared to the favorable drift configuration. The contributions of toroidal and poloidal main ion flows to Er, which are inferred from local neoclassical theory and the experiment, cannot account for these observed differences. Furthermore, it is found that in the L-mode, the intrinsic toroidal edge rotation decreases with increasing collisionality and it is co-current in the banana-plateau regime for all different drift configurations at AUG. This gives rise to a possible interaction of parallel Pfirsch–Schlüter flows in the SOL with the confined plasma. Thus, the different H-mode power threshold for the two drift configurations cannot be explained in the same way at AUG as suggested by LaBombard et al. [Phys. Plasmas 12, 056111 (2005)] for Alcator C-Mod. Finally, comparisons of Er profiles in favorable and unfavorable drift configurations at the respective confinement transitions show that also the Er gradients are all different, which indirectly indicates a different type or strength of the characteristic edge turbulence in the two drift configurations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Early-time measure in eternal inflation

In a situation like eternal inflation, where our data is replicated at infinitely-many other space-time events, it is necessary to make a prior assumption about our location to extract predictions. The principle of mediocrity entails that we live at asymptotic late times, when the occupational probabilities of vacua has settled to a near-equilibrium distribution. In this paper we further develop the idea that we instead exist during the approach to equilibrium, much earlier than the exponentially-long mixing time. In this case we are most likely to reside in vacua that are easily accessed dynamically. Using first-passage statistics, we prove that vacua that maximize their space-time volume at early times have: 1. maximal ever-hitting probability; 2. minimal mean first-passage time; and 3. minimal decay rate. These requirements are succinctly captured by an early-time measure. The idea that we live at early times is a predictive guiding principle, with many phenomenological implications. First, our vacuum should lie deep in a funneled region, akin to folding energy landscapes of proteins. Second, optimal landscape regions are characterized by relatively short-lived vacua, with lifetime of order the de Sitter Page time. For our vacuum, this lifetime is ~ 10 130 years, which is consistent with the Standard Model estimate due to Higgs metastability. Third, the measure favors vacua with small, positive vacuum energy. This can address the cosmological constant problem, provided there are sufficiently many vacua in the entire ensemble of funnels. As a concrete example, we study the Bousso-Polchinski lattice of flux vacua, and find that the early-time measure favors lattices with the fewest number of flux dimensions. This favors compactifications with a large hierarchy between the lightest modulus and all other Kähler and complex structure moduli.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Neutrino mass and mass ordering: no conclusive evidence for normal ordering

The extraction of the neutrino mass ordering is one of the major challenges in particle physics and cosmology, not only for its implications for a fundamental theory of mass generation in nature, but also for its decisive role in the scale of future neutrinoless double beta decay experimental searches. It has been recently claimed that current oscillation, beta decay and cosmological limits on the different observables describing the neutrino mass parameter space provide robust decisive Bayesian evidence in favor of the normal ordering of the neutrino mass spectrum. We further investigate these strong claims using a rich and wide phenomenology, with different sampling techniques of the neutrino parameter space. Contrary to the findings of Jimenez et al., no decisive evidence for the normal mass ordering is found. Neutrino mass ordering analyses must rely on priors and parameterizations that are ordering-agnostic: robust results should be regarded as those in which the preference for the normal neutrino mass ordering is driven exclusively by the data, while we find a difference of up to a factor of 33 in the Bayes factors among the different priors and parameterizations exploited here. An ordering-agnostic prior would be represented by the case of parameterizations sampling over the two mass splittings and a mass scale, or those sampling over the individual neutrino masses via normal prior distributions only. In this regard, we show that the current significance in favor of the normal mass ordering should be taken as 2.7σ (i.e. moderate evidence), mostly driven by neutrino oscillation data. Let us stress that, while current data favor NO only mildly, we do not exclude the possibility that this may change in the future. Eventually, upcoming oscillation and cosmological data may (or may not) lead to a more significant exclusion of IO.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

The role of B T -dependent flows on W accumulation at the edge of the confined plasma

Abstract Near-separatrix impurity accumulation between the crown and the outer midplane of tokamaks is a common feature in results from codes such as SOLPS-ITER and DIVIMP; however, experimental evidence of accumulation has only recently been obtained and is reported here. The codes find that the poloidal distribution of impurity ions in the scrape-off layer (SOL) depends primarily on toroidal field ( B T )-dependent parallel flow patterns of the background plasma and the parallel ion temperature gradient (∇ ‖ T ion ) force. Experimentally, Mach probes used in L-mode plasmas with favorable (for H-mode access) B T measure fast ( M ∼ 0.3–0.5) inner-target-directed (ITD) background plasma flows at the crown of single-null discharges. This study reports a set of DIVIMP simulations for two similar H-mode discharges from the DIII-D W metal rings campaign differing primarily in B T -direction to assess the effect that fast ITD flows have on the distribution of W ions in the SOL. It is found that for imposed ITD flows of M = 0.3, W ions that otherwise accumulate due to the ∇ ‖ T ion -force are largely flushed out. It is also found that doubling the radial diffusion coefficient from 0.3 to 0.6 m 2 s −1 prevents accumulation due to rapid cross-field transport into the far-SOL, where background plasma flows drain W ions to the divertors. Far-SOL W distributions from DIVIMP are then used to specify input to the impurity transport code 3DLIM, which is used to interpretively model collector probe (CP) deposition patterns measured in the ‘wall-SOL’. It is demonstrated that the deposition patterns are consistent with the DIVIMP predictions of near-SOL accumulation for the unfavorable- B T direction, and little/no accumulation for the favorable- B T direction. The wall-SOL CPs have thus provided the first experimental evidence, albeit indirect, of near-SOL W accumulation—finding it occurs for the unfavorable- B T direction only. For the favorable- B T direction, fast flows can largely prevent accumulation from occurring.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Role of electromagnetic gauge-field fluctuations in the selection between chiral and nematic superconductivity

Motivated by the observation of nematic superconductivity in several systems, we revisit the problem of the leading pairing instability of two-component unconventional superconductors on the triangular lattice—such as (p x , p y )-wave and (d x 2 –y 2 , d xy )-wave superconductors. Such a system has two possible superconducting states: the chiral state (e.g., p + ip or d + id), which breaks time-reversal symmetry, and the nematic state (e.g., p + p or d + d), which breaks the threefold rotational symmetry of the lattice. Weak-coupling calculations generally favor the chiral over the nematic superconducting state, raising the question of what mechanism can stabilize the latter. Here, we show that the electromagnetic field fluctuations can play a crucial role in selecting between these two states. Specifically, we derive and analyze the effective free energy for the two-component superconducting order parameter after integrating out the gauge-field fluctuations, which is formally justified if the spatial order parameter fluctuations can be neglected. A nonanalytic cubic term arises, as in the case of a conventional s -wave superconductor. However, unlike the latter, the cubic term depends on the relative phase and on the relative amplitudes between the two order parameter components, in such a way that it generally favors the nematic state. This result is a direct consequence of the fact that the stiffness of the superconducting order parameter is not isotropic. Competition with the quartic term, which favors the chiral state, leads to a renormalized phase diagram in which the nematic state displaces the chiral state over a wide region in the parameter space. Here, we analyze the stability of the fluctuation-induced nematic phase, generalize our results to tetragonal lattices, and discuss their applicability to candidate nematic superconductors, including twisted bilayer graphene.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Nonstoichiometric Salt Intercalation as a Means to Stabilize Alkali Doping of 2D Materials

Although doping with alkali atoms is a powerful technique for introducing charge carriers into physical systems, the resulting charge-transfer systems are generally not air stable. Here, in this paper, we describe computationally a strategy towards increasing the stability of alkali-doped materials that employs stoichiometrically unbalanced salt crystals with excess cations (which could be deposited during, e.g., in situ gating) to achieve doping levels similar to those attained by pure alkali metal doping. The crystalline interior of the salt crystal acts as a template to stabilize the excess dopant atoms against oxidation and deintercalation, which otherwise would be highly favorable. We characterize this doping method for graphene, NbSe 2 , and Bi 2 Se 3 and its effect on direct-to-indirect band gap transitions, 2D superconductivity, and thermoelectric performance. Salt intercalation should be generally applicable to systems which can accommodate this “ionic crystal” doping (and particularly favorable when geometrical packing constraints favor nonstoichiometry).

2-dimensional systems↗
DOE OSTI2023

Effects of the surface termination and the oxygen vacancy position on LaNiO 3 ultra-thin films: First-principles study

In this study, while ultra-thin layers of the LaNiO 3 film exhibit a remarkable metal-insulator transition as the film thickness becomes smaller than a few unit cell (u.c.), the formation of possible oxygen vacancies and their effects on the correlated electronic structure have been rarely studied using first principles. Here, we investigate the effects of the surface termination and the oxygen vacancy position on the electronic properties and vacancy energetics of LaNiO 3 ultra-thin films under the compressive strain using density functional theory plus U (DFT + U). We find that oxygen vacancies can be easily formed in the Ni layers with the NiO 2 terminated surface (0.5 u.c. and 1.5 u.c. thickness) compared to the structures with the LaO terminated surface and the in-plane vacancy is energetically favored than the out-of-plane vacancy. When two vacancy sites are allowed, the Ni square plane geometry is energetically more stable in most cases as two oxygen vacancies tend to stay near a Ni ion. Strong anisotropy between the in-plane and out-of-plane vacancy formation as well as the layer and orbital dependent electronic structure occur due to strain, surface termination, charge reconstruction, and quantum confinement effects. The in-plane vacancy of the NiO 2 terminated structure is favored since the released charge due to the oxygen vacancy can be easily accommodated in the d x 2 -y 2 orbital, which is less occupied than the d z 2 orbital. Remarkably, the oxygen vacancy structure containing the Ni square-plane geometry becomes an insulating state in DFT + U with a sizable band gap of 1.2 eV because the large crystal field splitting between d z 2 and d x 2 -y 2 orbitals in the square-plane favors an insulating state and the Mott insulating state is induced in other Ni sites due to strong electronic correlations. In the thin-film structure without oxygen vacancies, the strong correlation effect in DFT + U drives a pseudogap ground state at the Fermi energy, similarly as the experimental photo-emission spectra; however, the variation of the DFT + U electronic structure depending on the surface termination becomes weaker compared to those obtained in DFT.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Thermodynamics and kinetics of H adsorption and intercalation for graphene on 6 H -SiC(0001) from first-principles calculations

Previous experimental observations for H intercalation under graphene on SiC surfaces motivate the clarification of configuration stabilities and kinetic processes related to intercalation. From first-principles density-functional-theory calculations, we analyze H adsorption and intercalation for graphene on a 6H-SiC(0001) surface, where the system includes two single-atom-thick graphene layers: the top-layer graphene (TLG) and the underling buffer-layer graphene (BLG) above the terminal Si layer. Our chemical potential analysis shows that in the low-H coverage regime (described by a single H atom within a sufficiently large supercell), intercalation into the gallery between TLG and BLG or into the gallery underneath BLG is more favorable thermodynamically than adsorption on top of TLG. However, intercalation into the gallery between TLG and BLG is most favorable. We obtain energy barriers of about 1.3 and 2.3 eV for a H atom diffusing on and under TLG, respectively. From an additional analysis of the energy landscape in the vicinity of a step on the TLG, we assess how readily one guest H atom on the TLG terrace can directly penetrate the TLG into the gallery between TLG and BLG versus crossing a TLG step to access the gallery. Finally, we also perform density functional theory calculations for higher H coverages revealing a shift in favorability to intercalation of H underneath BLG and characterizing the variation with H coverage in interlayer spacings.

36 MATERIALS SCIENCE↗
DOE OSTI2025

When to vaccinate for seasonal influenza: check the peak forecast

Background Seasonal influenza infects 5-20% of people every year in the United States, resulting in hospitalizations, deaths, and adverse economic impacts. To mitigate these impacts, influenza vaccines are developed and distributed annually; however, growing evidence suggests that vaccine effectiveness (VE) wanes over the course of a flu season. Delaying influenza vaccination for older adults has attracted attention as a potential public health strategy. However, given the uncertainties in seasonal peak, vaccine effectiveness, and waning rates, postponing vaccination could also lead to increased morbidity, motivating an evaluation of a range of potential scenarios. The aim of this study was to investigate favorable age group-specific vaccination schedules that could lead to the greatest disease burden reduction. Methods We systematically investigated a broad range of vaccination start times for five age groups under six combinations of initial effectiveness and waning rates, based on influenza cases and vaccine uptake data from 10 influenza seasons. We defined the most favorable vaccination schedule as the one that resulted in the greatest reduction in disease burden. Results In scenarios with fast waning, all age groups benefit from delaying vaccination regardless of initial VE and peak timing. In scenarios with slower waning, results are mixed. For the ≥65 group, high initial VE and slow waning suggests that in early-peaking seasons, early vaccination most effectively reduces disease burden, while in late-peaking seasons delaying vaccination is most effective. For the ≥65 group in medium and low initial VE, and slow waning scenarios, delaying vaccination appears to prevent the greatest number of cases, regardless of whether the season peaks early or late. Conclusion The most favorable vaccination schedule is sensitive to changes in initial VE, waning rate, and peak timing. Given estimates of these quantities from statistical and immunological models and observations, our methods can inform vaccination recommendations in order to most effectively reduce the annual disease burden caused by seasonal influenza. Specifically, accurate peak timing forecasts for the upcoming season have the potential to guide decisions on when to vaccinate.

59 BASIC BIOLOGICAL SCIENCES↗