Search NASA⌕ Search

SEARCH · Search NASA

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

Elucidating molecular scale interactions underlying the freezing behavior of salt solutions in silica nanopores

Elucidating the influence of nanoscale confinement on the freezing behavior of salt solutions is of fundamental interest to environmental security and materials science and engineering. Specific structural information such as the coordination environment of ions in confined aqueous medium, effects on the extended hexagonal network of water molecules and their mutual non-bonding interactions in confinement are sparse in the literature along with the change in the dynamical characteristics of water in the presence of ions. To address these knowledge gaps, the current study is focused on investigating the influence of reduced dimensionality arising from nanoscale confinement on the structural evolution of salt solutions on freezing and the associated fluid–surface interactions. In this regard, operando wide-angle X-ray scattering (WAXS) measurements and classical molecular dynamics (MD) simulations are conducted with water and 0.5 M CaCl 2 , MgCl 2 and KCl solutions confined in 4 nm sized SBA-15 silica pores upon cooling from 300 K to 200 K. The freezing point of the salt solutions is depressed to 235 K which is about 10 K lower than that of confined water. The translational dynamics of confined salt solutions indicates shifts in the fragile-to-strong dynamical crossover to a lower temperature compared to pure water. The strong electrostatic attraction between the cations and the surrounding water molecules contributes to the freezing point depression in confined salt solutions. These insights unlock the molecular-scale basis and mechanisms underlying the freezing behavior of confined salt solutions.

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

Controlling reversible phase transitions in rare-earth nickelates for novel memory devices

Resistive switching in correlated complex oxides is lucrative for emerging applications in neuromorphic computing, and densely scaled non-volatile memory. Electrical conductance of such complex oxides can be controllable switched across multiple orders of magnitude by either (a) electroforming a conduction channel (e.g., in tungsten oxide), or (b) inducing Mott-Hubbard transition (e.g., in rare-earth nickelates)– both via controlled migration of defects (such as oxygen vacancies) under applied bias. Nevertheless, the promise of such defect-driven electronic transitions are far from realized due to a lack of fundamental understanding of the atomic-scale processes that underlie migration and spatiotemporal evolution of oxygen vacancies over nano-to-mesoscopic length/timescales under applied electric field. In this project, we employ a synergistic integration of density functional theory (DFT) calculations, ab initio/classical molecular dynamics (AIMD/CMD) simulations, machine learning (ML), precision synthesis, and multi-modal X-ray imaging experiments to address this knowledge gap. Such an integrated approach offers to elucidate the correlations between subtle structural distortion and oxidation states; treat localized charge carriers; describe defect/ion transport in the presence of electric field; and, in turn, greatly advance the current understanding of microstructural evolution in complex oxides under applied bias. The fundamental knowledge gained from this work will enable precise control over hierarchical defect structures and unravel new routes to manipulate resistance states in complex oxides. This, in turn, will accelerate design of novel devices with desired set of neural functionalities, and high-speed densely-scaled resistive random access memory technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Modulating Cu electrode microenvironments with MOF coatings: insights from molecular dynamics and electrochemical experiments of CO reduction

Metal-organic frameworks (MOFs) present a compelling strategy for tuning electrochemical interfaces by reshaping interfacial solvent structure. In this study, we examine how MOF coatings influence the microenvironment at copper electrodes during the CO electroreduction reaction (CORR) using a combined approach of molecular dynamics (MD) simulations and electrochemical experiments. Two MOFs, NU-901 and ZIF-8, are selected to explore the impact of pore size and channel hydrophobicity on electrochemical activity and interfacial concentration in acetonitrile (ACN) and dimethyl sulfoxide (DMSO) electrolytes. Electrochemical measurements reveal that MOF@Cu electrodes exhibit lower Faradaic efficiencies for CO hydrogenation products (ethylene and methane) compared to bare copper but have dramatic impacts on the interfacial microenvironment. NU-901, with its larger pores and strong interactions with DMSO, traps DMSO molecules and enhances CO coordination in DMSO but suppresses CORR selectivity in favor of the hydrogen evolution reaction (HER). ZIF-8, with smaller pores and hydrophobic channels, limits the interfacial water concentration, and, in ACN, promotes CO coordination. The simulations provide insights into how MOFs can act as physical modulators of reactant delivery and interfacial structure to control electrochemical microenvironments. This work highlights the value of molecular dynamics in uncovering how structural features of MOFs influence interfacial phenomena, even when catalytic performance is not directly improved.

Copper electrode microenvironments↗
DOE OSTI2025

An atomistic study connecting underlying dislocation behavior with superior mechanical properties of NiCoCr medium entropy alloy

NiCoCr-based medium-entropy alloy (MEA) with a simple face-centered cubic crystal phase exhibits excellent mechanical properties, often attributed to the synergy of multiple deformation mechanisms. However, the atomistic origin of their outstanding mechanical response, including microstructural evolution and dislocation behavior under varying strain-rates and orientation, remains unclear. In this work, we employ large-scale molecular dynamics (MD) simulations to investigate the changes in deformation mechanisms along three distinct orientations ([110], [111], [100]) under varying strain rates (1 ×10 8 /sec, 1 ×10 10 /sec, 1 ×10 12 /sec) in the NiCoCr MEA. The presence of the stair-rod and the Shockley partial dislocations under uniaxial tensile strain are found to play a key role in the formation of deformation twinning and ε-martensite, which positively correlates with strain-rate dependent dislocation analysis. These findings further establish the role of the dislocations in controlling the superior mechanical response and excellent fracture toughness of the NiCoCr MEA. Systematic transmission-electron microscopy tests performed on the [111]-oriented crystals, deformed at different strain levels, at room temperature provide clear evidence of both the extended stacking-fault and the stair rods, confirming the predicted microstructural features. Finally, this study offers key insights into the complex nucleation mechanisms of deformation twinning and ε-martensite, such as twinning – and transformation–induced plasticity (TWIP-TRIP), providing valuable guidelines for studying similar material classes.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Identifying impacts of contact tracing on HIV epidemiological inference from phylogenetic data

Abstract Robust sampling methods are foundational to inferences using phylogenies. Yet the impact of using contact tracing, a type of non-uniform sampling used in public health applications such as infectious disease outbreak investigations, has not been investigated in the molecular epidemiology field. To understand how contact tracing influences a recovered phylogeny, we developed a new simulation tool called SEEPS (Sequence Evolution and Epidemiological Process Simulator) that allows for the simulation of contact tracing and the resulting transmission tree, pathogen phylogeny, and corresponding virus genetic sequences. Importantly, SEEPS takes within-host evolution into account when generating pathogen phylogenies and sequences from transmission histories. Using SEEPS, we demonstrate that contact tracing can significantly impact the structure of the resulting tree, as described by popular tree statistics. Contact tracing generates phylogenies that are less balanced than the underlying transmission process, less representative of the larger epidemiological process, and affects the internal/external branch length ratios that characterize specific epidemiological scenarios. We also examined real data from a 2007–2008 Swedish HIV-1 outbreak and the broader 1998–2010 European HIV-1 epidemic to highlight the differences in contact tracing and expected phylogenies. Aided by SEEPS, we show that the data collection of the Swedish outbreak was strongly influenced by contact tracing even after downsampling, while the broader European Union epidemic showed little evidence of universal contact tracing, agreeing with the known epidemiological information about sampling and spread. Overall, our results highlight the importance of including possible non-uniform sampling schemes when examining phylogenetic trees. For that, SEEPS serves as a useful tool to evaluate such impacts, thereby facilitating better phylogenetic inferences of the characteristics of a disease outbreak. SEEPS is available at https://github.com/MolEvolEpid/SEEPS.

Virology↗
DOE OSTI2024

Photoinitiated Reactions of Molecules and Radicals in Molecular Beams

The UV photochemistry of organic molecules is a fundamental process that governs reactions in the atmosphere, synthetic chemistry, processing of organic aerosols, and biological damage in living tissues. In many organic molecules, the ensuing evolution involves pathways that are in competition, giving rise to different products, isomerization, coupling to other electronic states, and secondary reactive collisions. Because photodissociation is usually fast (picoseconds to microseconds), these processes are far from equilibrium and are controlled by kinetic competition and dynamical forces. The work accomplished was focused primarily on the photochemistry of alpha-keto carboxylic acids, a group of acids produced from natural sources, which are implicated in aerosol formation and biological processes. In the atmosphere, they are destroyed mainly by solar radiation. Studying their photochemistry has been surprisingly difficult because of the complexity of their excited electronic states, and the effect of collisions and secondary reactions. The second project investigated the lifetime of excited electronic states of pyrazine and picoline and had both experimental and theoretical aspects. The work focused on a model molecule, pyruvic acid (PA), because it was hypothesized that it should be a good source of the unstable carbene, methylhydroxycarbene (MHC). PA has an internal hydrogen bond that controls the evolution of its decomposition. The decomposition was studied following excitation to two of its lowest excited states reached by laser irradiation at 351 and 193 nm. The photodissociation dynamics in the absence and presence of collisions was monitored and compared. Understanding the UV photochemistry requires the use of complementary experimental approaches. Two methods were enlisted, which together generated a comprehensive and detailed set of results: (i) The time-sliced velocity map imaging (SVMI) instrument at USC was exploited to determine kinetic energy release of fragments, fast dissociation timescales, and internal state distributions of fragments for which Resonance Enhanced Multiphoton Ionization (REMPI) schemes exist; and (ii) The multiplexed photoionization mass spectrometer (MPIMS) setup developed at the Sandia Combustion Research Facility was used for product discovery, achieved by exploiting tunable narrowband VUV radiation at the Advanced Light Source (ALS), and to follow in real time their subsequent unimolecular and bimolecular reactions. The experimental conditions ranged from collisionless molecular beams to study nascent products to flow reactors of variable pressures to study subsequent bimolecular reactions of the products. The goals stated above were successfully accomplished. By exciting PA to its lowest excited state, MHC was identified as the only primary product and its isomerization to vinylalcohol and acetaldehyde was directly observed in real time. Moreover, we reported the first bimolecular reaction of MHC, which was with acetaldehyde, and identified its reaction product. This paves the way for the study of other reactions of this important carbene intermediate. When excitation was carried out at 193 nm, which imparted much higher energy to PA, many more products were directly observed in real time, some deriving from three-body dissociation. Quantitative branching ratios and secondary reactions of radical products were also determined and analyzed. The studies of pyrazine and picoline focused on the second excited state of these molecules and showed (via ionization studies) that their excited states were very short lived (<100 fs) because of efficient couplings to lower electronic states via vibronic coupling. Theoretical work on modeling the absorption spectrum and decay mechanisms of the excited states are in progress in collaboration with theoreticians.

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

Molecular dynamics study of hydrogen Cottrell atmosphere in aluminum: Influence of solute-solute interactions in the dislocation core

Cottrell atmospheres form when solute atoms segregate to dislocations, and is one of the most basic processes by which solutes alter mechanical properties of materials. In the case of hydrogen, Cottrell atmospheres are believed to contribute to hydrogen embrittlement. Here, in this study, we use direct molecular dynamics simulations to study the formation of hydrogen Cottrell atmospheres around edge dislocations in aluminum. Using coarse-graining techniques, we resolve the temporal and spatial evolution of the hydrogen concentration field. By comparing our results with theories of Cottrell atmosphere formation and thermodynamics, we show that a repulsive, concentration-dependent H-H interaction reduces the peak atmosphere concentration at the dislocation core. Furthermore, this interaction differs significantly from the H-H interaction that is experienced in the bulk away from a dislocation.

Aluminum↗
DOE OSTI2025

A Route to Design Novel Functional Peptides by Applying a Denoising Diffusional Model to mRNA Display Libraries

In vitro directed evolution techniques, such as mRNA display, enable peptide ligand discovery and optimization. However, physical libraries that rely on a genetic code can only search a small fraction of sequence space due to inherent biases in the genetic code and experimental limitations. To address this challenge, denoising diffusion implicit models (DDIMs) are applied to generate novel peptide ligands against B‐cell lymphoma extra‐large (Bcl‐x L ), a key cancer target. Starting with high‐throughput sequencing data from previous selections, a DDIM is trained to produce novel sequences with high affinity binding. Experimental validation confirms that most generated sequences are functionally equivalent to the original library members for Bcl‐x L binding and demonstrated comparable binding kinetics and affinity relative to the wildtype and nearest original neighbors. Importantly, this approach generated rare sequences not easily accessible via mutation and directed evolution. These results indicate that DDIMs can complement and expand directed evolution data, efficiently exploring underrepresented regions of sequence space. This approach provides a broadly applicable framework for accelerating ligand discovery and optimizing molecular properties across diverse targets.

Qi, Pearl [Mork Family Department of Chemical Engi↗
DOE OSTI2026

Adaptive_Evolution_2026

Thermoanaerobacterium saccharolyticum, an anaerobic and thermophilic bacterium capable of metabolizing sugar monomers and soluble oligomers into ethanol, has been proposed for use in consolidated bioprocessing in coculture with compatible cellulolytic bacteria such asClostridium thermocellum. Although the mixed acid fermentation of both of these strains has been engineered to produce ethanol as the major fermentation product, the maximum titer produced thus far byC. thermocellum, about 3% (w/v) ethanol, is half that produced byT. saccharolyticum. There is thus motivation to understand the mechanistic basis of the robust ethanol pathway inT. saccharolyticumso that key features can be recapitulated inC. thermocellumand other organisms.Previously, we characterized theindividualrole of the main genes responsible for electron transfer in the ethanol production ofT. saccharolyticum. However, the consequences of thecombinedloss of function of all these genes have not been investigated, nor has the way in which fermentative metabolism adapts to such constraints. In this work, we combined knockouts of ferredoxin nicotinamide oxidoreductase (fnor) genes (nfnAandnfnB) and hydrogenase genes (hydAandhfsD) and studied their effects on fermentation and grow. We showed that these genetic modifications together impair growth and decrease electron transfer from reduced ferredoxin, thereby redirecting flux from the pyruvate ferredoxin oxidoreductase enzyme to the pyruvate formate lyase enzyme. We also performed adaptive evolution of these mutants to rescue their growth, and most notably, we observed a single nucleotide variation in the alcohol dehydrogenaseadhAgene. Through molecular dynamics simulations and enzymatic assays, we determined that this point mutation causes a structural change that impairs the AdhA specificity for the NADPH cofactor and increases NADH-linked activity to restore redox balance. These findings consolidate our understanding of the functioning of electron transfer pathways in this organism.

alcohol dehydrogenase↗
DOE OSTI2026

Enhanced Catalytic Dechlorination of Polyvinyl Chloride (PVC) and H2 Production Enabled by Synergistic Gaδ+/Ga0 Active Sites in Liquid Metal Particles

Polyvinyl chloride (PVC) is ubiquitous yet challenging to recycle due to its tendency to thermally decompose above 250 °C, releasing toxic, corrosive chlorinated compounds, and its inability to melt. Here, we report a catalytic strategy for PVC upcycling at 160 °C using gallium liquid metal particles (Ga-LMP) featuring a dynamic Ga-GaOOH core–shell architecture. These catalysts enable concurrent dechlorination and hydrogen evolution, yielding up to 7% H2 (based on initial hydrogen atoms in PVC) along with a highly dechlorinated (>95%) carbonaceous solid and aqueous HCl. Mechanistic investigations combining X-ray photoelectron spectroscopy, infrared spectroscopy, solid-state NMR, inelastic neutron scattering, and ab initio molecular dynamics reveal a synergistic interplay between Gaδ+ sites in the GaOOH shell and metallic Ga0 in the core. Cationic Ga initiates C–Cl bond activation and HCl formation, while progressive reduction of the shell exposes Ga0 sites that promote C–H activation and H2 evolution. Control experiments with a Ga salt and bulk Ga liquid metal confirmed that neither oxidation state alone can achieve both transformations efficiently. This work establishes a dynamic dual-site paradigm for liquid metal catalysis, in which the in situ evolution and coexistence of oxidized and metallic species enable sequential and cooperative bond activation pathways. These findings provide a general design principle for novel liquid metal catalysts that target challenging polymer transformations under mild conditions.

Zingg, Benjamin [ORNL] (ORCID:0009000914530153)↗
DOE OSTI2024

The effect of low temperature on poly(3-methyl- N -vinylcaprolactam)- b -poly( N -vinylpyrrolidone) diblock copolymer nanovesicles assembled from all-aqueous media

Nanosized polymeric vesicles (polymersomes) self-assembled from double hydrophilic copolymers of poly(3-methyl-N-vinylcaprolactam) n -b-poly(N-vinylpyrrolidone) m (PMVC n -b-PVPON m ) using all aqueous media are a promising platform for biomedical applications, because of their superior stability over liposomes in vivo and high loading capacity. Herein, we explored the temperature-sensitive behavior of PMVC 58 -b-PVPON 65 vesicles using transmission electron microscopy (TEM), dynamic light scattering (DLS), atomic force microscopy (AFM), and small-angle neutron scattering (SANS) in response to lowering the solution temperature from 37 to 25, 20, 14 and 4 °C. The copolymer vesicles with an average size of 350 nm at 37 °C were assembled from the diblock copolymer dissolved in aqueous solution at 4 °C. We show that while the polymersome's size gradually decreases upon the temperature decrease from 37 to 4 °C, the average shell thickness increases from 17 nm to 25 nm, respectively. SANS study revealed that the PMVC 58 -b-PVPON 65 vesicle undergoes a gradual structure evolution from a dense-shell vesicle at 37–25 °C to a highly-hydrated shell vesicle at 20–14 °C to molecular chain aggregates at 4 °C. From SANS contrast matching study, this vesicle behavior is found to be driven by the gradual rehydration of PMVC block at 37–14 °C. The shell hydration at 20–14 °C also correlated with the 4.4-fold decrease in the relative fluorescence intensity from vesicle-encapsulated fluorescent dye, indicating ~80% of the dye release within 12 hours after the vesicle exposure to 14 °C. No significant (<5%) dye release was observed for the vesicle solutions at 37–20 °C, indicating excellent cargo retention inside the vesicles. Our study provides new fundamental insights on temperature-sensitive polymer vesicles and demonstrates that the copolymer assembly into polymersomes can be achieved by decreasing a copolymer aqueous solution temperature below 14 °C followed by solution exposure to ≥20 °C. This type of all-aqueous assembly, instead of nanoprecipitation from organic solvents or solvent exchange, can be highly desirable for encapsulating a wide range of biological molecules, including proteins, peptides, and nucleic acids, into stable polymer vesicles without a need for organic solvents for dissolution of the copolymers that are amphiphilic at physiologically relevant temperatures of 20–37 °C.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Electrochemical Oxidation of Seawater for Renewable Energy Generation and Storage in Coastal Areas

In this work, we demonstrated how the Clean Sustainable Electrochemical Treatment process can be used to process seawater while simultaneously generating molecular hydrogen. The specific goal of this project is to (1) demonstrate selective removal of chemical oxygen demand over the oxygen and chlorine evolution reactions, (2) optimize electrode composition to improve electrochemical performance, and (3) demonstrate long-term (>100 h) stability experiments.

08 HYDROGEN↗
DOE OSTI2025

Coupling of Charge Regulation and Geometry in Soft Ionizable Molecular Assemblies

The size, shape, and charge of structures, such as proteins and amphiphile assemblies, respond in an interconnected manner to solution ionic conditions. Here, we analyze assemblies of an amphiphile (C 16 K 2 ), with two ionizable amino acids [lysine (K)] coupled to a 16-carbon alkyl tail, via small-angle X-ray scattering (SAXS), nonlinear Poisson–Boltzmann theory (nl-PB), and hybrid Monte Carlo-molecular dynamics (MC-MD) simulations. SAXS revealed structural transitions from spherical micelles to cylindrical micelles to bilayers with increasing pH. By combining SAXS-determined structural information and nl-PB, we derived the molecular degree of ionization as a function of pH. The back-calculated titration curves matched the experimental data over an extended pH range, without adjustable parameters. Similarly, the SAXS data on the evolution of spherical micelle structure with ionic strength were combined with nl-PB and MC-MD to derive the bare and effective charges. MC-MD, which considered finite ion sizes, showed that bare and effective charges saturate quickly with increasing salt concentration. Furthermore, the calculated effective charges closely matched results from Zeta-potential measurements. The presented approach has advantages over customary methods for charge regulation, such as the Henderson–Hasselbalch (HH) or Hill models, where molecular ionization/deionization in assemblies is described by effective pKs that are distinct from the pK for isolated molecules. However, these models lack a physical explanation for these pK shifts. By contrast, our approach of combining structural details with an electrostatic model and simulations provides a more intuitive understanding of structure-charge coupling and a framework for understanding charge regulation in many synthetic and biological systems.

ionization↗
DOE OSTI2025

Unusual Ferromagnetic Band Evolution and High Curie Temperature in Monolayer 1T‐CrTe2 on Bilayer Graphene

2D van der Waals ferromagnets hold immense promise for spintronic applications due to their controllability and versatility. Despite their significance, the realization and in-depth characterization of ferromagnetic materials in atomically thin single layers, close to the true 2D limit, has been scarce. Here, a successful synthesis of monolayer (ML) 1T-CrTe2 is reported on a bilayer graphene (BLG) substrate via molecular beam epitaxy. Using angle-resolved photoemission spectroscopy and magneto-optical Kerr effect measurements, that the ferromagnetic transition is observed at the Curie temperature (TC) of 150 K in ML 1T-CrTe2 on BLG, accompanied by unconventional temperature-dependent band evolutions. The spectroscopic analysis and first-principle calculations reveal that the ferromagnetism may arise from Goodenough-Kanamori super-exchange and double-exchange interactions, enhanced by the lattice distortion and the electron doping from the BLG substrate. These findings provide pivotal insight into the fundamental understanding of mechanisms governing 2D ferromagnetism and offer a pathway for engineering higher TC in 2D materials for future spintronic devices.

Park, Kyoungree↗
DOE OSTI2024

High-temperature 205 Tl decay clarifies 205 Pb dating in early Solar System

Radioactive nuclei with lifetimes on the order of millions of years can reveal the formation history of the Sun and active nucleosynthesis occurring at the time and place of its birth. Among such nuclei whose decay signatures are found in the oldest meteorites, 205 Pb is a powerful example, as it is produced exclusively by slow neutron captures (the s process), with most being synthesized in asymptotic giant branch (AGB) stars. However, making accurate abundance predictions for 205 Pb has so far been impossible because the weak decay rates of 205 Pb and 205 Tl are very uncertain at stellar temperatures. To constrain these decay rates, we measured for the first time the bound-state β - decay of fully ionized 205 Tl 81+ , an exotic decay mode that only occurs in highly charged ions. The measured half-life is 4.7 times longer than the previous theoretical estimate and our 10% experimental uncertainty has eliminated the main nuclear-physics limitation. With new, experimentally backed decay rates, we used AGB stellar models to calculate 205 Pb yields. Propagating those yields with basic galactic chemical evolution (GCE) and comparing with the 205 Pb/ 204 Pb ratio from meteorites, we determined the isolation time of solar material inside its parent molecular cloud. We find positive isolation times that are consistent with the other s-process short-lived radioactive nuclei found in the early Solar System. Our results reaffirm the site of the Sun’s birth as a long-lived, giant molecular cloud and support the use of the 205 Pb– 205 Tl decay system as a chronometer in the early Solar System.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Active Learning Surrogates for Integrating Electron Microscopy and Computational Insights from Simulations in Autonomous Experiments

Artificial Intelligence (AI) combined with simulations and experiments has great potential to accelerate scientific discovery across technology and pharmaceuticals. However, the gap between simulations and experiments is challenging due to disparities in time and scale, making it difficult to estimate properties like energy and electronic states from experiments, and to provide feedback based on theoretical insights.Our research addresses the challenge by developing unique deep kernel based surrogate models that learns from microscopic images, mapping structural features to energy differences from defect formation. We start with full-training using simulated images to determine optimal settings, establishing a baseline for active learning. Using these settings from the baseline, active learning is trained, and predicts structures along simulation trajectories based on uncertainty and energetic stability, thus reducing data requirements, simulation time and computational costs. The results demonstrate that the model achieves a low average error margin of approximately 0.03 meV, indicating good performance. To enhance feature extraction and reconstruction capabilities, we developed an autoencoder-decoder as additional surrogate to create latent space to capture essential features, enabling precise comparisons between simulations and experiments. The results from this model achieved a reconstruction loss of around 0.2 and accurately reconstructed molecular structures.Overall, this work advances the steering of experiments through computational simulations by employing a surrogate models that actively predicts the trajectories of structural evolution, achieving time-to-solution comparable to experimental measurements.

Saranathan, Gayathri [Hewlett-Packard]↗
DOE OSTI2025

Proton radiation effects in indium oxide using cascade molecular dynamics simulations

Metal oxide (MO) semiconductors, characterized by their wide band gaps and notable charge transport properties, are promising candidates for electronic applications in extreme environments, including near-Earth space. However, atomistic simulations of radiation–matter interactions in MOs remain challenging due to the limitations of existing interatomic potentials, which often fail to capture both the short-range repulsive forces essential for radiation damage modeling and the long-range electrostatic effects governing defect evolution. In this work, we develop a customized interatomic potential tailored for radiation damage simulations in indium oxide (In 2 O 3 ) as a model system, a representative MO material. Our potential integrates the Ziegler-Biersack-Littmark potential to accurately describe short-range interactions with Buckingham and Coulombic potentials to account for long-range forces. We perform molecular dynamics simulations of low-energy proton irradiation using this custom potential. We employ the primary knock-on atom (PKA) cascade method to study atomic displacements and primary defect formation. Simulations were conducted for 1 keV proton irradiation in four randomly chosen directions, and PKA-driven defect analyses at 5, 10, and 15 keV to examine the effects of direction and energy level on damage generation. Our results provide insight into the impact of irradiation direction and energy level on the cascade evolution and defect formation mechanisms.

Atomistic simulations↗