Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular evolution”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

The development and evolution of biological AMS at Livermore: a perspective

Biological accelerator mass spectrometry (AMS) provides ultrasensitive carbon-14 isotopic analysis enabling a deeper understanding of human health concerns by enabling quantification of pharmacokinetics and other molecular endpoints directly in humans. It enables environmentally and human relevant studies of metabolic pathways through the use of very low concentrations of labeled metabolic substrates in cells and organisms. Here, we discuss why AMS is an important tool for the biosciences, the development and evolution of biological AMS at Livermore and discuss technical refinements that will improve the efficiency of operation for the measurement of ultra-trace levels of 14 C, which, long term, will enable greater ease of use and sample throughput.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Multi-scale signaling and tumor evolution in high-grade gliomas

Although genomic anomalies in glioblastoma (GBM) have been well studied for over a decade, its 5-year survival rate remains lower than 5%. We seek to expand the molecular landscape of high-grade glioma, composed of IDH-wildtype GBM and IDH-mutant grade 4 astrocytoma, by integrating proteomic, metabolomic, lipidomic, and post-translational modifications (PTMs) with genomic and transcriptomic measurements to uncover multi-scale regulatory interactions governing tumor development and evolution. Applying 14 proteogenomic and metabolomic platforms to 228 tumors (212 GBM and 16 grade 4 IDH-mutant astrocytoma), including 28 at recurrence, plus 18 normal brain samples and 14 brain metastases as comparators, reveals heterogeneous upstream alterations converging on common downstream events at the proteomic and metabolomic levels and changes in protein-protein interactions and glycosylation site occupancy at recurrence. Recurrent genetic alterations and phosphorylation events on PTPN11 map to important regulatory domains in three dimensions, suggesting a central role for PTPN11 signaling across high-grade gliomas.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Engineering Enantiocomplementary Protoglobins for Stereoconvergent Construction of N -Alkylated α-Aminoketones

The synthesis of enantiopure compounds from a mixture of E/Z alkenes represents a notable challenge in synthetic chemistry. While enzymes excel in achieving unparalleled selectivity, their inherent specificity often confines activity to a single stereoisomeric substrate, consequently restricting the overall efficiency of such transformations. Here, we demonstrate that protoglobin-derived hemoproteins can catalyze stereoconvergent intermolecular amination using simple N-alkyl hydroxylamines as nitrene precursors, a transformation which remains elusive in synthetic chemistry. These engineered enzymes process E/Z mixtures of silyl enol ethers, enabling the precise incorporation of N-alkyl amino moieties (−NHAlkyl) into diverse molecular structures (up to 79% yield and 95% ee). Two complementary protoglobin variants were engineered using directed evolution to enable enantiodivergent synthesis of both enantiomers of α-aminoketones. This enzymatic platform achieves stereoconvergent and enantiodivergent transformations, facilitating the conversion of simple chemicals into an array of valuable pharmaceutical compounds featuring aminoketone functionalities.

Alcohols↗
DOE OSTI2023

Dynamic Molecular Investigation of the Solid-Electrolyte Interphase of an Anode-Free Lithium Metal Battery Using In Situ Liquid SIMS and Cryo-TEM

Solid electrolyte interphase (SEI) has been widely perceived to play a critical role in the stable cycling of rechargeable batteries. However, associated with the fragile and air-sensitive nature of the SEI layer, delineation of the formation process and the nature of SEI remains a big challenge. Here, we use in situ liquid time-of-flight secondary ion mass spectroscopy (TOF-SIMS), cryo- transmission electron microscope (TEM) and density functional theory (DFT) calculation to delineate molecular process on the formation of SEI layer under the dynamic operating condition. We discover that the onset potential for SEI layer formation and the thickness of the SEI show dependence on the solvation shell structure. Using LiCoO 2 as a cathode and Cu film as an anode, the SEI is noticed to start to form at around 2.0 V and reach its final thickness (irreversible part, ~ 40-50 nm) at about 3.0 V in the 1 M LiPF 6 –EC/DMC electrolyte, while for the case of 1 M LiFSI–DME, the SEI starts to form at around 1.5 V and reaches its final thickness (~ 20 nm) at about 2.0 V. The in situ TOF-SIMS clearly indicates the outer SEI layer formation and dissipation upon charging and discharging, implying a continued evolution of electrolyte structure with extended cycling. In conclusion, the present work establishes a direct correlation between the molecular signature of SEI layer with solvation feature of electrolytes in lithium batteries, providing insights for tailoring SEI layer toward improved electrochemical properties of lithium batteries.

25 ENERGY STORAGE↗
DOE OSTI2020

Active Learning A Neural Network Model For Gold Clusters & Bulk From Sparse First Principles Training Data

Small metal clusters are of fundamental scientific interest and of tremendous significance in catalysis. These nanoscale clusters display diverse geometries and structural motifs depending on the cluster size; a knowledge of this size-dependent structural motifs and their dynamical evolution has been of longstanding interest. Given the high computational cost of first-principles calculations, molecular modeling and atomistic simulations such as molecular dynamics (MD) has proven to be an important complementary tool to aid this understanding. Classical MD typically employ predefined functional forms which limits their ability to capture such complex size-dependent structural and dynamical transformation. Neural Network (NN) based potentials represent flexible alternatives and in principle, well-trained NN potentials can provide high level of flexibility, transferability and accuracy on-par with the reference model used for training. A major challenge, however, is that NN models are interpolative and requires large quantities (similar to 10 4 or greater) of training data to ensure that the model adequately samples the energy landscape both near and far-from-equilibrium. A highly desirable goal is minimize the number of training data, especially if the underlying reference model is first-principles based and hence expensive. In this work, we introduce an active learning (AL) scheme that trains a NN model on-the-fly with minimal amount of first-principles based training data. Our AL workflow is initiated with a sparse training dataset (similar to 1 to 5 data points) and is updated on-the-fly via a Nested Ensemble Monte Carlo scheme that iteratively queries the energy landscape in regions of failure and updates the training pool to improve the network performance. Using a representative system of gold clusters, we demonstrate that our AL workflow can train a NN with similar to 500 total reference calculations. Using an extensive DFT test set of similar to 1100 configurations, we show that our AL-NN is able to accurately predict both the DFT energies and the forces for clusters of a myriad of different sizes. Our NN predictions are within 30 meV/atom and 40 meV/angstrom of the reference DFT calculations. Moreover, our AL-NN model also adequately captures the various size-dependent structural and dynamical properties of gold clusters in excellent agreement with DFT calculations and available experiments. We finally show that our AL-NN model also captures bulk properties reasonably well, even though they were not included in the training data.

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

Phase evolution in two-phase alloys during severe plastic deformation

Herein, phase evolution in FCC metals with strongly interacting alloy components during severe plastic deformation is investigated using molecular dynamics simulations. Specifically, we study the alloy microstructure in steady state, nucleation and growth of precipitates in supersaturated alloys, and the decomposition of precipitates in undersaturated alloys. The results are compared to a modified effective temperature model, providing a physical understanding for the atomic processes underlying the model and a perspective on its strengths and weaknesses. Key observations in this work are nucleation and growth of precipitates during SPD at a temperature of 100 K; Gibbs-Thomson-like behavior relating solubility to precipitate size under steady-state shearing; a direct relationship between the effective temperature and the shear modulus; and the importance of cluster agglomeration during precipitate growth. The study also reveals that the mechanisms of forced chemical mixing depends on precipitate size, adding complications for effective temperature models describing inhomogeneous systems. The simulations are shown to provide good semiquantitative agreement with experimental findings reported in the literature.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Experimental and theoretical investigation of the crystalline surface, film, and interface properties of antiperovskite Mn 3 GaN grown by molecular beam epitaxy on MgO(001)

Here, we present a study of the epitaxial growth, characterization, and theoretical modeling of thin film antiperovskite Mn 3 GaN, an antiferromagnetic material with kagome structure which is grown on MgO (001) substrates using N-plasma-assisted molecular beam epitaxy. Reflection high energy electron diffraction is used to assess the in-plane evolution of the film structure during growth, and the surface is investigated in-situ using scanning tunneling microscopy and Auger electron spectroscopy. These results are combined with precision measurements done ex-situ determining the film lattice constants using a combination of x-ray diffraction with reciprocal space mapping and scanning transmission electron microscopy. Overall, a uniform, homogeneous film with an atomically smooth vacuum surface and atomically sharp substrate interface is found having very small in-plane tensile strain and mild out-of-plane compressive strain. First-principles theoretical calculations are applied in order to ascertain the lowest energy models for both the Mn 3 GaN surface and the Mn 3 GaN/MgO film/substrate interface. Models including MnGa versus MnN surface layers and MnGa versus MnN interfacial layers are considered as functions of both the Mn and Ga chemical potentials. The predictions are discussed in comparison to the experimental results. The overall findings suggest that Mn 3 GaN on MgO(001) is a viable epitaxial film which can be further explored in connection with antiferromagnetic spintronics.

Density functional theory↗
DOE OSTI2025

Protocol for engineering poly(ethylene terephthalate) hydrolases via directed evolution using a high-throughput screening assay

Poly(ethylene terephthalate) (PET) hydrolases, which depolymerize PET to its monomers, have gained attention for their potential to facilitate bio-industrial recycling of this waste plastic. Here, we present a protocol for screening large, random mutagenesis enzyme libraries simultaneously for enhanced activity, solubility, and stability. We outline steps for library construction, screening using plate-based split GFP and model substrate assays, and determination of enzyme thermostability. We then detail procedures for validation assays on PET substrates and characterization of final variants.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Miscibility of rock and ice in the interiors of water worlds

Abstract Super-Earths and sub-Neptunes are the most common planet types in our galaxy. A subset of these planets is predicted to be water worlds, bodies that are rich in water and poor in hydrogen gas. The interior structures of water worlds have been assumed to consist of water surrounding a rocky mantle and iron core. In small planets, water and rock form distinct layers with limited incorporation of water into silicate phases, but these materials may interact differently during the growth and evolution of water worlds due to greater interior pressures and temperatures. Here, we use density functional molecular dynamics (DFT-MD) simulations to study the miscibility and interactions of enstatite (MgSiO 3 ), a major end-member silicate phase, and water (H 2 O) at extreme conditions in water world interiors. We explore pressures ranging from 30 to 120 GPa and temperatures from 500 to 8000 K. Our results demonstrate that enstatite and water are miscible in all proportions if the temperature exceeds the melting point of MgSiO 3 . Furthermore, we performed smoothed particle hydrodynamics simulations to demonstrate that the conditions necessary for rock-water miscibility are reached during giant impacts between water-rich bodies of 0.7–4.7 Earth masses. Our simulations lead to water worlds that include a mixed layer of rock and water.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Investigating shock-induced chemical reactions in Ni+Al multilayers: A continuum-based mesoscale approach with Arrhenius kinetics and artificial thermal conduction

A new continuum-based mesoscale modeling approach for shock-induced chemical reactions (SICRs) in Ni+Al multilayers is demonstrated in Sandia’s shock physics hydrocode, CTH. The approach utilizes Arrhenius-type kinetics and artificial thermal conduction. Our work builds upon previous efforts to parameterize equations of state for Ni x Al y [J. Appl. Phys. 137, 075102 (2025)], as well as simulations of inert shocks in realistic 2D microstructures [J. Appl. Phys. 137, 225301 (2025)]. To calibrate the reaction kinetics, pairs of the reaction coordinate, R′, vs time are extracted from the molecular dynamics (MD) literature. Here, the MD-informed kinetics are used to simulate the dynamic evolution of pressure and temperature in 2D mesoscale simulations. Overall, the MD-informed kinetics obtained for planar interfaces are too slow, as initial reaction is not observed on a nanosecond time scale. Even with quasi-periodic shock focusing leading to the formation of so-called “hot-spots,” the hot spots are unable to grow and coalesce using the fitted Arrhenius rate constants. However, by increasing the rate constants by two orders of magnitude, SICRs are observed at a shock pressure near 30 GPa, which is supported by experiments. Consequently, these mesoscale simulations suggest that unresolved shear-based mechanical mixing might possibly account for the discrepancies in kinetic rates, with shock-generated intense perturbations, interfacial vortical flows, and elevated temperatures serving as favorable reaction conditions. Future work will calibrate a shear-dependent reaction rate from the MD simulations with realistic interfaces that are wavy, diffuse, and disordered.

Kittell, David E. [Sandia National Laboratories (S↗
DOE OSTI2023

The structure of fly Teneurin-m reveals an asymmetric self-assembly that allows expansion into zippers

Teneurins are conserved cell adhesion molecules essential for embryogenesis and neural development in animals. Key to teneurin function is the ability of its extracellular region to form homophilic interactions in cis and/or in trans. However, our molecular understanding of teneurin homophilic interaction remains largely incomplete. Here, we showed that an extracellular fragment of Teneurin-m, the major teneurin homolog in flies, behaves as a homodimer in solution. The structure of Teneurin-m revealed that the transthyretin-related domain from one protomer and the β-propeller domain from the other mediates Teneurin-m self-association, which is abolished by point mutation of conserved residues. Strikingly, this architecture generates an asymmetric oligomerization interface that enables expansion of Teneurin-m into long zipper arrays reminiscent of protocadherins. An alternatively spliced site that exists only in vertebrates and regulates homophilic interaction in mammalian teneurins overlaps with the fly Teneurin-m self-association interface. Our work provides a molecular understanding of teneurin homophilic interaction and sheds light on its role in teneurin function throughout evolution.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Off-Equilibrium Reactivity of Boron-Enriched Metal Diboride Surfaces in Electroreduction Conditions

Boron-based materials, featuring B-dependent reactivity and diverse phases, are emerging as promising catalyst systems. However, the catalytic mechanism on many borides remains poorly understood due to complex surface reconstructions under reaction conditions. Here, we investigate the MoB 2 surface in conditions of hydrogen evolution reaction in acidic media, using grand canonical global optimization, grand canonical density functional theory, ab initio molecular dynamics, free energy surface sampling, and an analytical model for electrochemical barrier evaluation. We propose a boron-enrichment strategy to tune the surface reactivity of the hexagonal face of MoB 2 . We reveal the dynamic nature of the B-enriched surface under H coverage and kinetic trapping of the system in the metastable regime with an extensive examination of the deactivation pathways. The metastable center B site on B-enriched surfaces, featuring buckled-up configuration and a usual relaxation effect, is found to be highly active toward HER via the Volmer–Heyrovsky mechanism. In conclusion, this work demonstrates how off-equilibrium behaviors can arise from the interplay between adsorbate coverage and surface reconstruction on a seemingly simple surface, and we present a theoretical framework and computational workflows to address these behaviors, along with other realistic complexities, in kinetics simulations.

Adsorption↗
DOE OSTI2023

Limited impact of jet-induced feedback in the multi-phase nuclear interstellar medium of 4C12.50

4C12.50 (IRAS 13451+1232) at z = 0.122 is an ultraluminous infared radio galaxy that has often been proposed as a prime candidate for the link between ultraluminous infared galaxies and young radio galaxies. It is also an interesting target to investigate whether and how radio-induced feedback affects the evolution of galaxies in the early phases of radio activity. We study, in detail for the first time, the hot (≥1500 K) molecular gas in 4C12.50. The potential impact of the radio jet on this gas phase, as well as on the star formation activity, are investigated. We also study the ionised (including coronal) gas as traced by the near-infrared lines. Using near-infrared long slit spectroscopy obtained with EMIR on GTC and X-shooter on VLT, we analyse the emission line spectrum of the ionised, coronal, and, specially, the hot molecular gas in the western nucleus hosting the compact radio jet. Based on high spatial resolution ALMA CO(2–1) data, we also revise the location of 4C12.50 in the Kennicutt-Schmidt diagram in order to investigate whether star formation is suppressed. 4C12.50 hosts (2.1 ± 0.4) × 10 4 M ⊙ of hot molecular gas. An unusually high rotational temperature T rot = 3020 ± 160 K is inferred. The molecular gas mass obeys a power-law temperature distribution, dM H2 /dT ∝ T –5 , from T ~300 K and up to ~3000 K. Both results support the idea that shocks (probably induced by the radio jet) contribute to the heating and excitation of the hot molecular gas. A molecular outflow is not detected. The coupling of the outflowing ionised and neutral outflows with the hot molecular gas is poor. Contrary to other studies, we claim that there is no evidence for star formation suppression in this object. If radio-induced feedback can regulate the star formation activity in galaxies, 4C12.50 is a promising candidate to reveal this phenomenon in action. However, we find no solid evidence for a current or past impact of this mechanism on the evolution of this system, neither by clearing out the dusty central cocoon efficiently, nor by suppressing the star formation activity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Dislocation loop evolution in Kr-irradiated ThO 2

Here, the early stage of microstructural evolution of ThO 2 , under krypton irradiation at 600, 800, and 1000°C, was investigated using in situ transmission electron microscopy (TEM). Dislocation loops grew faster, whereas their number density decreased with increasing irradiation temperature. Loop density was found to decrease with ion dose. Interstitial dislocation loops, including Frank loops with Burgers vector of $\textit{a}/3\langle111\rangle$ and perfect loops with Burgers vector of $\textit{a}/2\langle110\rangle$, were determined by traditional TEM and atomic resolution–scanning TEM techniques. Atomistic and mesoscale level modeling are performed to interpret experimental observations. The migration energy barriers of defects in ThO 2 were calculated using density-functional theory. The energetics of different dislocation loop types were studied using molecular dynamics simulations. Loop density and diameter were analyzed using a kinetic rate theory model that considers stoichiometric loop evolution. This analysis reveals that loop growth is governed by the mobility of cation interstitials, whereas loop nucleation is determined by the mobility of anion defects. Lastly, a rate theory model was used to extract the diffusion coefficients of thorium interstitials, oxygen interstitials, and vacancies.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Structurally Driven Environmental Degradation of Friction in MoS2 Films

Abstract We report an investigation of the friction mechanisms of MoS 2 thin films under changing environments and contact conditions using a variety of computational and experimental techniques. Molecular dynamics simulations were used to study the effects of water and molecular oxygen on friction and bonding of MoS 2 lamellae during initial sliding. Characterization via photoelectron emission microscopy (PEEM) and Kelvin probe force microscopy (KPFM) were used to determine work function changes in shear modified material within the top few nanometers of MoS 2 wear scars. The work function was shown to change with contact conditions and environment, and shown by density functional theory (DFT) calculations and literature reports to be correlated with lamellae size and thickness of the basally oriented surface layer. Results from nanoscale simulations and macroscale experiments suggest that the evolution of the friction behavior of MoS 2 is linked primarily to the formation or inhibition of a basally oriented, molecularly thin surface film with long-range order.

Curry, John F.↗
DOE OSTI2024

Enhanced mobility of dislocation network nodes and its effect on dislocation multiplication and strain hardening

Understanding plastic deformation of crystals in terms of the fundamental physics of dislocations has remained a grand challenge in materials science for decades. To overcome this, the Discrete Dislocation Dynamics (DDD) method has been developed, but its lack of atomistic resolution leaves open the possibility that certain key mechanisms may be overlooked. Here, by comparing large-scale Molecular Dynamics (MD) with DDD simulations performed under identical conditions we uncover significant discrepancies in the predicted strength and microstructure evolution in BCC crystals under high-strain rate conditions. These are traced to unexpected behaviors of dislocation network nodes forming at dislocation intersections, that can move in ways not previously anticipated as revealed by MD. Once these newfound freedoms of nodal motion are incorporated, DDD simulations begin to closely match plastic evolution observed in MD. This additional mechanism of motion whereby non-screw dislocations can change their glide plane profoundly affects fundamental processes of dislocation multiplication, recovery and storage that define strength of metals.

36 MATERIALS SCIENCE↗
DOE OSTI2021

In Silico Discovery of Multistep Chemistry Initiated by a Conical Intersection: The Challenging Case of Donor–Acceptor Stenhouse Adducts

Detailed mechanistic understanding of multistep chemical reactions triggered by internal conversion via a conical intersection is a challenging task that emphasizes limitations in theoretical and experimental techniques. We present a discovery-based, hypothesis-free computational approach based on first-principles molecular dynamics to discover and refine the switching mechanism of donor–acceptor Stenhouse adducts (DASAs). We simulate the photochemical experiment in silico, following the “hot” ground state dynamics for 10 ps after photoexcitation. Using state-of-the-art graphical processing units-enabled electronic structure calculations we performed in total ~2 ns of nonadiabatic ab initio molecular dynamics discovering (a) critical intermediates that are involved in the open-to-closed transformation, (b) several competing pathways which lower the overall switching yield, and (c) key elements for future design strategies. Finally, our dynamics describe the natural evolution of both the nuclear and electronic degrees of freedom that govern the interconversion between DASA ground-state intermediates, exposing significant elements for future design strategies of molecular switches.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗