Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density modification”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Modulating Surface Anionic Redox Chemistry toward Highly Stable Li-Rich Cathodes with Negligible Oxygen Loss

Low initial Coulombic efficiency and severe capacity/voltage fading during cycling caused by serious irreversible oxygen release, especially in the initial cycle, and resultantly induced unstable electrode/electrolyte interfacial chemistry, largely prohibit the commercial application of high-capacity Li-rich layered oxide cathodes (LLOs). Here, in this work, a dual reductive gas interface cotreatment strategy is applied to regulate the lattice oxygen redox activity and reversibility with a multiple defective structure design including Li/O/TM (TM = transition metal) vacancies and the intrinsic TM doping as well as a full-surface protective layer, which can suppress the irreversible TM migration and then undesirable phase transformation, resisting the corrosion of electrolyte during cycling effectively. Importantly, the introduced reversible SO 3 2- /SO 4 2- redox couple that provides extra capacity compensation could alleviate the distortion of oxygen-central octahedral structure and structural collapse caused by immoderate oxygen oxidation. Thus, the lattice oxygen redox chemistry is optimized, with negligible oxygen loss during the initial cycle. And the designed AS-LLO cathode with greatly enhanced structure stability shows high-capacity retentions of 99.2% at 0.3C after 100 cycles and 82.4% even after 1000 cycles at 5C. This work provides a guideline for manipulating the oxygen redox chemistry to achieve long-lifespan Li-rich layered oxide cathodes for high-energy-density lithium batteries.

25 ENERGY STORAGE↗
NASA NTRS2011

Modification of Jupiter's Stratosphere Three Weeks After the 2009 Impact

Infrared spectroscopy sensitive to thermal emission from Jupiter's stratosphere reveals effects persisting 23 days after the impact of a body in late July 2009. Measurements obtained on 2009 August II UT at the impact latitude of 56 S (planetocentric), using the Goddard Heterodyne Instrument for Planetary Wind and Composition mounted on the NASA Infrared Telescope Facility, reveal increased ethane abundance and the effects of aerosol opacity. An interval of reduced thermal continuum emission at 11. 744 ~lm is measured ~60o-80 towards planetary east of the impact site, estimated to be at 3050 longitude (System Ill). Retrieved stratospheric ethane mole fraction in the near vicinity of the impact site is enhanced by up to -60% relative to quiescent regions at this latitude. Thermal continuum emission at the impact site, and somewhat west of it, is significantly enhanced in the same spectra that retrieve enhanced ethane mole fraction. Assuming that the enhanced continuum brightness near the impact site results from thermalized aerosol debris blocking contribution from the continuum formed in the upper troposphere and indicating the local temperature, then continuum emission by a haze layer can be approximated by an opaque surface inserted at the 45-60 mbar pressure level in the stratosphere in an unperturbed thermal profile, setting an upper limit on the pressure and therefore a lower limit on the altitude of the top of the impact debris at this time. The reduced continuum brightness east of the impact site can be modeled by an opaque surface near the cold tropopause, which is consistent with a lower altitude of ejecta/impactor-formed opacity or significantly lesser column density of opaque haze material. The physical extent of the observed region of reduced continuum implies a minimum average velocity of 21 m/s transporting material prograde (planetary east) from the impact.

Fast, Kelly E.↗
DOE OSTI2022

Pressure-Driven Changes in the Electronic Bonding Environment of GeO 2 Glass above Megabar Pressures

Noncrystalline oxides under pressure undergo gradual structural modifications, highlighted by the formation of a dense noncrystalline network topology. The nature of the densified networks and their electronic structures at high pressures may account for the mechanical hardening and the anomalous changes in electromagnetic properties. Despite its importance, direct probing of the electronic structures in amorphous oxides under compression above the Mbar pressure (>100 GPa) is currently lacking. Here, we report the observation of pressure-driven changes in electronic configurations and their delocalization around oxygen in glasses using inelastic X-ray scattering spectroscopy (IXS). In particular, the first O K-edge IXS spectra for compressed GeO 2 glass up to 148 GPa, the highest pressure ever reached in an experimental study of GeO 2 glass, reveal that the glass densification results from a progressive increase of oxygen proximity. While the triply coordinated oxygen [3] O is dominant below ~50 GPa, the IXS spectra resolve multiple edge features that are unique to topologically disordered [4] O upon densification above 55 GPa. Topological compaction in GeO 2 glass above 100 GPa results in pronounced electronic delocalization, revealing the contribution from Ge d-orbitals to oxide densification. Strong correlations between the glass density and the electronic configurations beyond the Mbar conditions highlight the electronic origins of densification of heavy-metal-bearing oxide glasses. Furthermore, current experimental breakthroughs shed light on the direct probing of the electronic density of states in high-Z oxides above 1 Mbar, offering prospects for studies on the pressure-driven changes in magnetism, superconductivity, and electronic transport properties in heavy-metal-bearing oxides under compression.

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

Characterizing Density and Spatial Distribution of Trap States in Ta 3 N 5 Thin Films for Rational Defect Passivation

Tantalum nitride (Ta 3 N 5 ) has gained significant attention as a potential photoanode material, yet it has been challenged by material quality issues. Defect-induced trap states are detrimental to the performance of any semiconductor material. Beyond influencing the performance of Ta 3 N 5 films, defects can also accelerate the degradation in water during desired electrochemical applications. Defect passivation has provided an enormous boost to the development of many semiconductor materials but is currently in its infancy for Ta 3 N 5 . This is in part due to a lack of experimental understanding regarding the spatial and energetic distribution of trap states throughout Ta 3 N 5 thin films. Here, we employ drive-level capacitance profiling (DLCP) to experimentally resolve the spatial and energetic distribution of trap states throughout Ta 3 N 5 thin films. The density of deeper energetic traps is found to reach ∼2.5 to 6 × 10 22 cm –3 at the interfaces of neat Ta 3 N 5 thin films, over an order of magnitude greater than the bulk. In addition to the spatial profile of deep trap states, we report neat Ta 3 N 5 thin films to be highly n-type in nature, owning a free carrier density of ∼9.74 × 10 17 cm –3 . This information, coupled with the present understanding of native oxide layers on Ta 3 N 5 , has facilitated the rational design of a targeted passivation strategy that simultaneously provides a means for catalyst immobilization. Loading catalyst via silatrane moieties suppresses the density of defects at the surface of Ta 3 N 5 thin films by two orders of magnitude, while also reducing the free carrier density of films by over one order of magnitude, effectively dedoping the films to ∼2.40 × 10 16 cm –3 . The surface passivation of Ta 3 N 5 films translates to suppressed defect-induced trapping and recombination of photoexcited carriers, as determined through absorption, photoluminescence, and transient photovoltage. Here, this illustrates how developing a deeper understanding of the distribution and influence of defects in Ta 3 N 5 thin films has the potential to guide future works and ultimately accelerate the integration and development of high-performance Ta 3 N 5 thin film devices.

defect passivation↗
DOE OSTI2020

Novel PBX formulations containing thermally-expandable microspheres for on-demand control of explosive behavior

In this work, we present the formulation and analysis of inert Plastic-Bonded eXplosive (PBX)-surrogates loaded with thermally expandable microspheres (TEMs). TEMs consist of a thermoplastic acrylonitrile shell (10-50 µm diameter) encapsulating an inert low boiling hydrocarbon. Upon heating, the TEMs expand as the shell softens while the hydrocarbon gasifies, increasing the internal pressure and expanding the particle by as much as 120 vol%. We hypothesize that TEM expansion within a PBX will introduce tunable changes in local density and porosity, and ultimately the shock sensitivity. This paper focuses on microstructural analysis of surrogate PBX-TEM formulations, both before and after thermally- induced TEM expansion. The formulations were analyzed by scanning electron microscopy (SEM), ultra-small angle neutron scattering (USANS), and X-ray Computed Tomography (x-ray CT) under various thermal conditions. Evidence suggests that TEM expansion within a PBX is a viable method to impart microstructural changes that are predicted to have a measurable, on-demand effect on the initiation sensitivity. Future work includes development of PBX-TEM formulations with HMX and shock initiation experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Supported Electrophilic Organoruthenium Catalyst for the Hydrosilylation of Olefins

A series of supported electrophilic organoruthenium complexes has been synthesized via surface organometallic chemistry (SOMC) techniques and applied to the selective hydrosilylation of olefins. The air-sensitive 16e - complex Cp*RuMes(PCy 3 ) (1) (Cp* = pentamethylcyclopentadienyl, Mes = mesityl) was synthesized by the treatment of Cp*RuCl(PCy 3 ) with mesityl Grignard MesMgBr. This species was chemisorbed onto sulfated zirconia SO 4 /ZrO 2 , but the resulting material was inactive toward cyclohexene hydrosilylation with phenylsilane. Instead, Cp*RuMes(PCy 3 ) was treated with phenylsilane (PhSiH 3 ) to provide a ruthenium disilyl hydride complex Cp*RuH(SiH 2 Ph) 2 (PCy 3 ) (3), which was fully characterized by NMR spectroscopy and single-crystal X-ray diffraction. Grafting this species onto SO 4 /ZrO 2 resulted in the formation of phenylsilane along with the surface electrophilic species [Cp*RuH(R)(X-SiHPh)(PCy 3 )] (R = H, O 3 S-O or O 3 Zr-O; 4a, 4b, X = O 3 S-O, and O 3 Zr-O, respectively) as the major species. Material 4 was characterized via a combination of spectroscopic techniques including dynamic nuclear polarization (DNP)-enhanced solid-state NMR spectroscopy, diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS), X-ray absorption spectroscopy (XAS), and density function theory (DFT) calculations. Further, capping the remaining acid sites on 4 with Me 3 Si-SiMe 3 provides 5, which significantly reduces side reactions, such as olefin isomerization and silane redistribution. Catalyst 5 is a highly robust and selective hydrosilylation catalyst and can be recycled up to 5 times without significant diminishment of activity. Exclusive anti-Markovnikov regiochemistry, cis-addition selectivity, and the inactivity of secondary and tertiary silanes provide support for the proposed Glaser-Tilley mechanism involving cationic ruthenium silylene species analogous to homogeneous systems.

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

Ultrafast Formation of Jahn–Teller Polarons Revealed by State-Selective Excitation in Correlated Spinel Co 3 O 4

Jahn–Teller polarons are quasiparticles that stem from symmetry breaking and strong local electron–phonon coupling. They originate from an excess charge carrier being dressed by a local lattice distortion, caused by the Jahn–Teller effect, and they critically impact electrical, structural, and magnetic properties in transition metal oxides. The observation of the microscopic steps involved in their formation is essential for enabling control over material properties through the targeted activation of local, site-specific modifications with light pulses. While Jahn–Teller distortion associated with polaron formation was predicted to contribute significantly to changes in electronic band gap and optical properties in Co 3 O 4 , its experimental observation remains elusive, requiring signatures of local symmetry reduction. In this work, we demonstrate Jahn–Teller polaron formation in spinel Co 3 O 4 . By exciting electronic transitions at 3.10 eV and 1.55 eV, we target either the Oh Cobalt(III) or the Td Cobalt(II) ions, and drive the subsequent coherent responses of the system through two different pathways. For the former, we demonstrate that ligand-to-metal charge transfer leads to Jahn–Teller polaron formation, which is linked to the deformation potential and magnetoelastic coupling. For the latter, we identify the coherent excitation of a T 2g phonon mode launched by on-site d-d electronic transitions. Key to our observations is the ability to target site-specific electronic excitations in spinel Co 3 O 4 using ultrafast optical pulses, while monitoring the ensuing low-energy collective modes through the coherent time-domain response of the material. Our approach, which combines the comparative analysis of experimental fingerprints with the support from density functional theory calculations, is broadly applicable to systems in which Jahn–Teller polaron physics has been theoretically predicted but remains experimentally unverified, and it underscores the potential of electronic-state targeting as a route to selectively excite and probe quasiparticle dynamics in solids.

Lattices↗
DOE OSTI2025

Single-Crystallization of O3-Type Layered Oxide Cathode for Na-Ion Battery

The development of high-energy-density Na-ion batteries places significant demands on single-crystal layered oxide cathodes, especially for further high-voltage, solid-state battery scenarios. In the O3-type structure, due to the original sluggish Na ion diffusion kinetics (approximately 1 order of magnitude lower than that of Li-ion), and further hindrance against diffusion kinetics caused by single-crystal architecture, these inherent defects lead to the decline in the electrochemical performance. Herein, we demonstrated that the single crystallization of O3-type NaNi 1/3 Fe 1/3 Mn 1/3 O 2 cathode (d 50 = 5.04 μm) aggravates surface-to-bulk phase inhomogeneity distribution, which is attributed to the uneven Na ions extraction. Moreover, the Na-depletion of the surface/shell region not only aggravates Na ion diffusion resistance but also leads to a higher valence state of transition-metal elements (e.g., Ni/Fe) near the surface of the single-crystal particle, which further compromises the cathode-electrolyte interface stability. Furthermore, not limited to revealing the challenges, tuning the particle size and moderating quasi-single-crystal strategies have been proven to effectively mitigate the negative uneven distributions of Na ions, phases, and valence/oxidative states, resulting in efficient modification for single crystallization of Na-layered oxide cathodes.

36 MATERIALS SCIENCE↗
NASA NTRS1995

Analysis of the progressive failure of brittle matrix composites

This report investigates two of the most common modes of localized failures, namely, periodic fiber-bridged matrix cracks and transverse matrix cracks. A modification of Daniels' bundle theory is combined with Weibull's weakest link theory to model the statistical distribution of the periodic matrix cracking strength for an individual layer. Results of the model predictions are compared with experimental data from the open literature. Extensions to the model are made to account for possible imperfections within the layer (i.e., nonuniform fiber lengths, irregular crack spacing, and degraded in-situ fiber properties), and the results of these studies are presented. A generalized shear-lag analysis is derived which is capable of modeling the development of transverse matrix cracks in material systems having a general multilayer configuration and under states of full in-plane load. A method for computing the effective elastic properties for the damaged layer at the global level is detailed based upon the solution for the effects of the damage at the local level. This methodology is general in nature and is therefore also applicable to (0(sub m)/90(sub n))(sub s) systems. The characteristic stress-strain response for more general cases is shown to be qualitatively correct (experimental data is not available for a quantitative evaluation), and the damage evolution is recorded in terms of the matrix crack density as a function of the applied strain. Probabilistic effects are introduced to account for the statistical nature of the material strengths, thus allowing cumulative distribution curves for the probability of failure to be generated for each of the example laminates. Additionally, Oh and Finney's classic work on fracture location in brittle materials is extended and combined with the shear-lag analysis. The result is an analytical form for predicting the probability density function for the location of the next transverse crack occurrence within a crack bounded region. The results of this study verified qualitatively the validity of assuming a uniform crack spacing (as was done in the shear-lag model).

Thomas, David J.↗
DOE OSTI2025

First-Principles Insights into the Thermocatalytic Cracking of Ammonia-Hydrogen Blends on Fe(110). 2. Kinetics

Ammonia (NH 3 ) is an energy-rich molecule that is routinely synthesized from nitrogen (N 2 ) and hydrogen (H 2 ). NH 3 ’s more favorable physical properties compared to H 2 suggests it may offer a way to more conveniently store, transport, and, when needed, extract H 2 via thermal decomposition. However, the high kinetic barrier and endoergicity to decompose to H 2 and N 2 require high temperatures. The standard reaction free energy indicates nearly 100% thermodynamic conversion to the diatomic molecules only at ~673 K and higher. However, even at these temperatures, a catalyst, e.g., iron (Fe), is needed for favorable kinetic conversion. Here, in this study, we explore via density functional theory the kinetics of NH 3 decomposition on the most stable facet of body-centered cubic Fe, namely, (110), under typical high-temperature and finite-pressure operando conditions. We predict coverage-dependent energetics of elementary surface reactions, often neglected in atomic-scale modeling. From these models, we find the recombinative desorption of adsorbed N as N 2 is rate-determining at 573.15–773.15 K and even at an extreme case of 1173.15 K. From microkinetic modeling, we find that the steady-state turnover frequencies (TOFs) for N 2 and H 2 generation rates (r$_{H_2}$) depend exponentially on temperature. The catalyst achieves a steady-state TOF of 36.4 s –1 and an r$_{H_2}$ of 0.107 μmol cm –2 s –1 for a feed of 1.8 bar NH 3 with 0.2 bar H 2 at 1173.15 K. However, at 773.15 K, with the same feed composition and velocity, the steady-state TOF and r$_{H_2}$ decrease to 0.14 s –1 and 4.10 × 10 –4 μmol cm –2 s –1 , respectively, as the process is significantly hindered by slow N 2 desorption. Although at first glance counterintuitive, our simulations suggest that surface modifications that reduce Fe’s reactivity toward NH x species should enhance its overall NH 3 decomposition activity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRSDate not supplied

The BioSensor Instrument Beyond LEIA: a Versatile Platform for Lunar Biology

Introduction: The BioSensor is a deep-space-compatible automated microfluidic culturing instrument. While originally designed to measure the effects of deep space radiation on yeast growth for the BioSentinel mission, it has the potential to host a diverse range of life science experiments with single- and/or multi-celled organisms and can be adapted to interface with a diverse range of platforms in both crewed and uncrewed settings. It is therefore a leading candidate for hosting life sciences experimentation associated with a lunar surface habitat. BioSensor design: The function of the BioSensor is to monitor the growth and metabolic activity of samples in batch fluidic culture, without the need for crew involvement. The current configuration houses organisms in 16 wells within 16 microfluidic cards, accommodating a total of 256 samples, replicates, and controls. Each well has an optical system including three LEDs and a photodiode detector to measure absorbance at three wavelengths, enabling measurement of optical density, color change in dyes such as alamarBlue, and bioproduction of pigmented compounds. Organisms are loaded into fluidic wells and air-dried for storage during integration, launch, and transit, then activated by the introduction of culture medium from storage bags via manifolds that fill one card at a time. Temperature is controlled by individual card heaters, and timing of all activities (fluidics fills, optical measurements, temperatures) is directed by an experiment script. The self-contained BioSensor payload is roughly 4U in volume; with electrical/mechanical/thermal interface, e.g. for operation on ISS or a lunar lander, as well as a linear energy-transfer (LET) charged-particle radiation spectrometer, the entire system is closer to 6U. BioSentinel and LEIA: Flown on the ISS and in a deep-space free flyer for BioSentinel, the BioSensor has been modified for use in the LEIA mission, including improvements to reduce the sensitivity to lengthy launch delays. LEIA will monitor yeast growth in the radiation and reduced-gravity environment of the lunar surface no earlier than 2026, on a CLPS lander [4]. Changes include accommodating additional culture media and an additional LED color for a new biological assay (bioproduction of carotenoids-- dietary antioxidants), as well as modifications to the housing to allow late-load biology changeout and improved isolation between fluidics and electronics. Future prospects: Future work with the BioSensor, beyond LEIA, will include expanding the range and diversity of organisms and assays accommodated. Preliminary work has demonstrated the growth of Arabidopsis seedlings in BioSensor fluidic cards, including optical measurements of growth rate over time. Minor modifications could allow measurement of phenotypes related to photosynthetic capacity in both plants and cyanobacteria. The experimental capabilities of the BioSensor could be dramatically increased by introducing the capability for fluorescence measurements, and/or the design of novel biological assays using luminescence. The BioSensor can also be adapted for new platforms and experiment settings; in addition to free-flyer, ISS, and CLPS lander, a preliminary design concept has been developed for crewed deployment directly to the lunar surface. The instrument could be accommodated inside a lunar habitat, where its automated operation would make it an excellent candidate for experiments from fundamental investigations into the response of organisms to lunar surface conditions to applied-science purposes such as screening engineered strains of various organisms for bioproduction capability.

J A Lee↗
NASA NTRS2008

Evaluation of Cloud Physical Properties of ECMWF Analysis and Re-Analysis (ERA-40 and ERA Interim) against CERES Tropical Deep Convective Cloud Object Observations

This study presents an approach that converts the vertical profiles of grid-averaged cloud properties from large-scale models to probability density functions (pdfs) of subgrid-cell cloud physical properties measured at satellite footprints. Cloud physical and radiative properties, rather than just cloud and precipitation occurrences, of assimilated cloud systems by the European Center for Medium-range Weather Forecasts (ECMWF) operational analysis (EOA) and ECMWF Re-Analyses (ERA-40 and ERA Interim) are validated against those obtained from Earth Observing System satellite cloud object data for January-August 1998 and March 2000 periods. These properties include ice water path (IWP), cloud-top height and temperature, cloud optical depth and solar and infrared radiative fluxes. Each cloud object, a contiguous region with similar cloud physical properties, is temporally and spatially matched with EOA and ERA-40 data. Results indicate that most pdfs of EOA and ERA-40 cloud physical and radiative properties agree with those of satellite observations of the tropical deep convective cloud-object type for the January-August 1998 period. There are, however, significant discrepancies in selected ranges of the cloud property pdfs such as the upper range of EOA cloud top height. A major discrepancy is that the dependence of the pdfs on the cloud object size for both EOA and ERA-40 is not as strong as in the observations. Modifications to the cloud parameterization in ECMWF that occurred in October 1999 eliminate the clouds near the tropopause but shift power of the pdf to lower cloud-top heights and greatly reduce the ranges of IWP and cloud optical depth pdfs. These features persist in ERA-40 due to the use of the same cloud parameterizations. The downgrade of data assimilation technique and the lack of snow water content information in ERA-40, not the coarser horizontal grid resolution, are also responsible for the disagreements with observed pdfs of cloud physical properties although the detection rates of cloud object occurrence are improved for small size categories. A possible improvement to the convective parameterization is to introduce a stronger dependence of updraft penetration heights with grid-cell dynamics. These conclusions will be rechecked using the ERA Interim data, due to recent changes in the ECMWF convective parameterization (Bechtold et al. 2004, 2008). Results from the ERA Interim will be presented at the meeting.

Xu, Kuan-Man↗
NASA NTRS2005

Evaluation of a New Cloud Droplet Activation Parameterization wtih in Situ Data from CRYSTAL-FACE and CSTRIPE

The accuracy of the 2003 prognostic, physically based aerosol activation parameterization of A. Nenes and J. H. Seinfeld (NS) with modification introduced by C. Fountoukis and A. Nenes in 2005 (modified NS) is evaluated against extensive microphysical data sets collected on board the Center for Interdisciplinary Remotely Piloted Aircraft Studies (CIRPAS) Twin Otter aircraft for cumuliform and stratiform clouds of marine and continental origin. The cumuliform cloud data were collected during NASA's Cirrus Regional Study of Tropical Anvils and Cirrus Layers-Florida Area Cirrus Experiment (CRYSTAL-FACE, Key West, Florida, July 2002), while the stratiform cloud data were gathered during Coastal Stratocumulus Imposed Perturbation Experiment (CSTRIPE, Monterey, California, July 2003). In situ data sets of aerosol size distribution, chemical composition, and updraft velocities are used as input for the NS parameterization, and the evaluation is carried out by comparing predicted cloud droplet number concentrations (CDNC) with observations. This is the first known study in which a prognostic cloud droplet activation parameterization has been evaluated against a wide range of observations. On average, predicted droplet concentration in adiabatic regions is within -20% of observations at the base of cumuliform clouds and -30% of observations at different altitudes throughout the stratiform clouds, all within experimental uncertainty. Furthermore, CDNC is well parameterized using either a single mean updraft velocity w or by weighting droplet nucleation rates with a Gaussian probability density function of w. This study suggests that for nonprecipitating warm clouds of variable microphysics, aerosol composition, and size distribution the modified NS parameterization can accurately predict cloud droplet activation and can be successfully implemented for describing the aerosol activation process in global climate models.

Meskhidze, Nicholas↗
NASA NTRS1998

The Postshock Chemical Lifetimes of Outflow Tracers and a Possible New Mechanism to Produce Water Ice Mantles

We have used a coupled time-dependent chemical and dynamical model to investigate the lifetime of the chemical legacy in the wake of C-type shocks. We concentrate this study on the chemistry of H2O and O2, two molecules which are predicted to have abundances that are significantly affected in shock-heated gas. Two models are presented: (1) a three-stage model of preshock, shocked, and postshock gas; and (2) a Monte Carlo cloud simulation where we explore the effects of stochastic shock activity on molecular gas over a cloud lifetime. For both models we separately examine the pure gas-phase chemistry as well as the chemistry including the interactions of molecules with grain surfaces. In agreement with previous studies, we find that shock velocities in excess of 10 km/s are required to convert all of the oxygen not locked in CO into H2O before the gas has an opportunity to cool. For pure gas phase models the lifetime of the high water abundances, or "H2O legacy," in the postshock gas is approximately (4-7) x 10(exp 5) yr, independent of the gas density. A density dependence for the lifetime of H2O is found in gas-grain models as the water molecules deplete onto grains at the depletion timescale. Through the Monte Carlo cloud simulation we demonstrate that the time-average abundance of H2O, the weighted average of the amount of time gas spends in preshock, shock, and postshock stages, is a sensitive function of the frequency of shocks. Thus we predict that the abundance of H2O, and to a lesser extent O2, can be used to trace the history of shock activity in molecular gas. We use previous large-scale surveys of molecular outflows to constrain the frequency of 10 km/s shocks in regions with varying star formation properties and discuss the observations required to test these results. We discuss the postshock lifetimes for other possible outflow tracers (e.g., SiO and CH3OH) and show that the differences between the lifetimes for various tracers can produce potentially observable chemical variations between younger and older outflows. For gas-grain models we find that the abundance of water-ice on grain surfaces can be quite large and is comparable to that observed in molecular clouds. This offers a possible alternative method to create water mantles without resorting to grain surface chemistry: gas heating and chemical modification due to a C-type shock and subsequent depletion of the gas-phase species onto grain mantles.

Bergin, Edwin A.↗
DOE OSTI2024

Co-Hydrotreatment of Pyrolytic Lignin and Waste Cooking Oil to Produce Hydrocarbons

The co-hydrotreatment of pyrolysis bio-oils and vegetable oils could facilitate the integration of bioderived fuels into existing fuel infrastructures without significant modifications. Co-hydrotreatment studies with different PL/WCO blend ratios (0, 10, 20, 30, and 40 wt %) were conducted over the NiMo/γ-Al 2 O 3 catalyst. The coke formation value for WCO was 0.7 wt % and ranged between 1.5 and 2.5 wt % with the increase of pyrolytic lignin in the blend. The data suggest that coke is formed from both the sugar- and lignin-derived oligomers since the coke yield reported in this study is comparable with the coke yield obtained by coprocessing the bio-oil fraction without light oxygenate compounds, based on previous studies. All blends are recommended for co-hydrotreatment based on coke yield. The resulting organic phase was distilled into hydrocarbons at <150 °C, 150 to 250 °C, and 250 to 350 °C. The organic product’s overall hydrocarbon distribution was strongly dependent on the pyrolytic lignin concentration. The yield of the distillation cut with a boiling point range between 250 and 350 °C was higher in the hydrotreated samples, with more than 20% during distillation. Two-dimensional GC×GC was used to determine carbon distribution on the combustibles. FTIR and UV fluorescence spectra showed that most feedstocks were converted to potential combustibles for transportation applications after hydrotreatment. The distillation cuts properties, such as density, viscosity, and surface tension, were reported and compared to transportation fuel properties. Hydrocarbon yields reveal the potential of cohydrotreatment to upgrade pyrolytic lignin into valuable products. Furthermore, these findings provide useful insights into utilizing pyrolytic lignin to produce promising combustibles suitable for transportation applications.

co-hydrotreatment↗
DOE OSTI2025

Tailoring the Physicochemical Properties of Nb Thin Films via Surface Engineering Methods

The modification of surface oxide layers formed on niobium (Nb) thin films via chemical mechanical planarization (CMP) and accelerated neutral atom beam (ANAB) processing provides a promising route toward tailoring their emergent properties and performance when used as superconducting qubits. Here, in this study, we show that CMP- and ANAB-formed Nb oxides are significantly thinner and smoother than the native oxide, as revealed by transmission electron microscopy (TEM) and atomic force microscopy. Scanning TEM and energy-dispersive X-ray spectroscopy along with X-ray photoelectron spectroscopy identified an oxidation gradient within the native and surface-engineered oxides. The topside layer is dominated by Nb 5+ (Nb 2 O 5 ), with various Nb suboxides present closer to the oxide/metal interface. Time-of-flight secondary ion mass spectrometry (ToF-SIMS) depth profiling confirmed the presence of an oxygen content gradient and demonstrated the enhanced resistance of the CMP- and ANAB-formed oxides to oxygen surface exchange and subsequent diffusion via 18 O 2 isotopic labeling experiments. ToF-SIMS also identified an interfacial layer containing trapped hydrogen (H)-containing species at the Nb oxide/metal interface. In situ ToF-SIMS and TEM revealed migration of the H/OH interfacial layer coinciding with decomposition of the surface oxide. Furthermore, our density functional theory calculations indicated that both H from moisture present in ambient air and bulk H in Nb films tend to segregate at the interface. These findings underscore the importance of understanding surface oxidation mechanisms, hydrogen incorporation, and their impact on the designed functionalities of Nb-based devices.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Geodyn Material Library: Pseudocap models for dry porous tocks

This report describes the second edition of the Pseudocap Strength models for porous rocks implemented in GEODYN material library. The first model was developed in 2007 and calibrated for concrete. Then, the model parameters were calibrated based on triaxial tests reported for limestones and sandstones of various porosities. In these models some key parameters were chosen as functions of the reference porosity,Φ. Since then multiple modifications were implemented in the model, therefore, it has been recalibrated for some common porous materials such as limestones, sandstones, alluvium, tuffs and granite. Two types of models are described in this repot. The first type (called Pseudocap Model or PM) is for rocks from a specific location. Parameters were calibrated for several specific geologic materials. The second type (called Generic Pseudocap Model or GPM) is useful for the sites where only basic information (rock type, porosity) is available. Generic models include built-in correlations between porosities and other mechanical properties observed for certain rock types. The models of both types were validated by comparing not only to quasi-static triaxial tests for these materials but also to shock Hugoniot data and spherical explosion data for some materials. All models were derived in the frame of isotropic plasticity. They are designed to be used in explicit finite element/difference codes. Tangent stiffness tensor is not provided but can be calculated numerically for the model to be used in implicit finite element codes. For an isotropic material the stress can be decomposed into volumetric and deviatoric parts. The volumetric part is modeled using an Equation of state (EOS) which calculates the pressure and the bulk sound speed as functions of the internal specific energy and density. Here a simple, Mie-Gruneisen EOS is presented, but tabulated EOS (LEOS) provided by the library can be used as well. The stress is limited by the yield surface which depends on three invariants of the stress tensor and specific internal energy. In addition, to capture the strain-rate dependence a simple multiplier is used for the yield surface which depends on the equivalent plastic strain rate. The failure surface (the ultimate yield, Y f , defined later) is chosen in the Hoek-Brown form, commonly used in rock mechanics. It includes measurable parameters such as Unconfined Compressive Strength (UCS) as well as scale parameters characterizing the quality of the rock such as GSI (Geologic Strength Index). Thus, even though the model is calibrated for small samples it offers a way to extrapolate the strength to the field scale using geological characterization of the rock mass. The model captures effects of brittle-ductile transition in rocks by introducing a cap multiplier to the yield function. The rate of dilatancy (bulking) is proportional to the slope of the yield surface affected by the cap. Therefore, it takes place only at low confinements when the pressure is less than the brittle-ductile transition pressure, P BD . On the contrary, the porous compaction takes place at pressures higher than P BD . The cap moves as the porosity is compacted or new porosity is generated due to dilatancy. The porous compaction is modeled using an evolution equation which includes deviatoric stress so that the onset of compaction corresponds to the cap surface. The model captures effects of shear-enhanced compaction which is an important for porous rocks. Section 2 describes the modeling framework and Section 3 presents the model calibration procedure. Section 4 compares experimental data for various rocks versus model predictions. The model parameters used for this comparison are given in Appendix. The files with material constants are available with the latest GEODYN material library distribution.

58 GEOSCIENCES↗
DOE OSTI2024

Computing equilibrium free energies through a nonequilibrium quench

Many methods to accelerate sampling of molecular configurations are based on the idea that temperature can be used to accelerate rare transitions. These methods typically compute equilibrium properties at a target temperature using reweighting or through Monte Carlo exchanges between replicas at higher temperatures. A recent paper [G. M. Rotskoff and E. Vanden-Eijnden, Phys. Rev. Lett. 122, 150602 (2019)] demonstrated that accurate equilibrium densities of states can also be computed through a nonequilibrium “quench” process, where sampling is performed at a higher temperature to encourage rapid mixing and then quenched to lower energy states with dissipative dynamics. Here, we provide an implementation of the quench dynamics in LAMMPS and evaluate a new formulation of nonequilibrium estimators for the computation of partition functions or free energy surfaces (FESs) of molecular systems. We show that the method is exact for a minimal model of N-independent harmonic springs and use these analytical results to develop heuristics for the amount of quenching required to obtain accurate sampling. We then test the quench approach on alanine dipeptide, where we show that it gives an FES that is accurate near the most stable configurations using the quench approach but disagrees with a reference umbrella sampling calculation in high FE regions. We then show that combining quenching with umbrella sampling allows the efficient calculation of the free energy in all regions. Moreover, by using this combined scheme, we obtain the FES across a range of temperatures at no additional cost, making it much more efficient than standard umbrella sampling if this information is required. Finally, we discuss how this approach can be extended to solute tempering and demonstrate that it is highly accurate for the case of solvated alanine dipeptide without any additional modifications.

Chemistry↗