Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Molecular machine”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

2D nanoconfinement distorts the solvation structure of hydroxide but not of hydronium

Understanding ion-specific behavior in nanoconfined water is essential for controlling charge transport and selectivity in two-dimensional membranes. Motivated by recent experiments revealing anomalous dielectric and transport behavior of water confined between hBN sheets, we use machine-learning-accelerated first-principles molecular dynamics to investigate the interfacial propensities of hydronium and hydroxide ions under similar confinement. We find that hydronium remains interfacial across all confinement regimes, whereas hydroxide shifts toward the interior as the environment becomes more bulk-like. This contrasting behavior reflects the combined influence of hydrogen bonding, interfacial water layering, and the polarization of hBN, which collectively stabilize hydronium at the surface while making hydroxide slightly more favorable within the structured interior. These findings expose an asymmetry in ion-surface coupling and establish a microscopic origin for hydronium’s enhanced interfacial affinity. The results provide mechanistic insight into ion partitioning in two-dimensional channels and highlight the collective structuring of confined water as a key determinant of interfacial ion thermodynamics.

organic↗
DOE OSTI2021

Computational optimal transport for molecular spectra: The fully discrete case

In this study, the use of computational optimal transport is investigated as a tool for comparing two molecular spectra. Unlike other techniques for comparing molecular spectra in a pattern-recognition framework, transport distances simultaneously encode information about line positions and intensities. In addition, it is shown that transport distances are a useful alternative to Euclidean distances as Euclidean distances are based on line-by-line comparisons, while transport distances reflect broader features of molecular spectra and adequately compare spectra with different resolutions. This paper includes a tutorial on the use of optimal transport and investigates several well-chosen examples to illustrate the utility of computational optimal transport for comparing molecular spectra.

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

Computational optimal transport for molecular spectra: The semi-discrete case

Comparing a discrete molecular spectrum to a continuous molecular spectrum in a quantitative manner is a challenging problem, for example, when attempting to fit a theoretical stick spectrum to a continuous spectrum. In this paper, the use of computational optimal transport is investigated for such a problem. In the optimal transport literature, the comparison of a discrete and a continuous spectrum is referred to as semi-discrete optimal transport and is a situation where a metric such as least-squares may be difficult to define except under special conditions. The merits of an optimal transport approach for this problem are investigated using the transport distance defined for the semi-discrete case. A tutorial on semi-discrete optimal transport for molecular spectra is included in this paper, and several well-chosen synthetic spectra are investigated to demonstrate the utility of computational optimal transport for the semi-discrete case. Among several types of investigations, we include calculations showing how the frequency resolution of the continuous spectrum affects the transport distance between a discrete and a continuous spectrum. We also use the transport distance to measure the distance between a continuous experimental electronic absorption spectrum of SO 2 and a theoretical stick spectrum for the same system. The comparison of the theoretical and experimental SO 2 spectra also allows us to suggest a theoretical value for the band origin that is closer to the observed band origin than previous theoretical values.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

Extended Rice–Thomson analysis and atomistic simulations revealing grain boundary effects on fracture in refractory high-entropy alloys

Significance This work serves to extend the fundamental ductile vs. brittle fracture theory, specifically the Rice–Thomson criterion, by introducing a grain boundary ahead of an initiating crack which propagates at an oblique angle to impinge the boundary. Atomistic fracture simulations on two refractory complex concentrated alloys, the brittle NbMoTaW and the ductile Nb 45 Ta 25 Ti 15 Hf 15 , demonstrate qualitative correspondence with the extended Rice–Thomson criterion and experimental observations. Abstract Understanding how grain boundaries mediate fracture remains a critical challenge in designing ductile, high-performance refractory alloys. Here, we extend the Rice–Thomson criterion to account for the angle between cracks and the impinging grain boundaries (GBs), capturing the competition between intergranular fracture and dislocation-mediated plasticity. Using machine learning interatomic potentials, we performed molecular statics simulations to probe fracture mechanisms in nanocrystalline NbMoTaW and Nb 45 Ta 25 Ti 15 Hf 15 , each with two different grain sizes, revealing trends consistent with experimental observations and the extended Rice model. Comparison with averaged R-curves for bulk samples demonstrates that GBs enhance ductility in Nb 45 Ta 25 Ti 15 Hf 15 in both grain sizes investigated. In contrast, GBs only locally improve fracture resistance in NbMoTaW when cracks are temporarily pinned at GBs inclined at high angles from the crack, but generally promote brittle intergranular fracture. These contrasting behaviors are attributed to differences in GB cohesion, reflecting clear alloying trends that align with ab-initio calculations and trends observed experimentally. Our results bridge classical fracture theory, atomistic simulations, and experimental observations, providing a comprehensive understanding of the fracture mechanisms in nanocrystalline refractory complex concentrated alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Generalized representative structures for atomistic systems

A new method is presented to generate atomic structures that reproduce the essential characteristics of arbitrary material systems, phases, or ensembles. Previous methods allow one to reproduce the essential characteristics (e.g. the chemical disorder) of a large random alloy within a small crystal structure. The ability to generate small representations of random alloys, along with the restriction to crystal systems, results from using the fixed-lattice cluster correlations to describe structural characteristics. A more general description of the structural characteristics of atomic systems is obtained using complete sets of atomic environment descriptors. These are used within for generating representative atomic structures without restriction to fixed lattices. A general data-driven approach is provided here utilizing the atomic cluster expansion (ACE) basis. The N-body ACE descriptors are a complete set of atomic environment descriptors that span both chemical and spatial degrees of freedom and are used within for describing atomic structures. The generalized representative structure (GRS) method presented within generates small atomic structures that reproduce ACE descriptor distributions corresponding to arbitrary structural and chemical complexity. It is shown that systematically improvable representations of crystalline systems on fixed parent lattices, amorphous materials, liquids, and ensembles of atomic structures may be produced efficiently through optimization algorithms. With the GRS method, we highlight reduced representations of atomistic machine-learning training datasets that contain similar amounts of information and small 40–72 atom representations of liquid phases. The ability to use GRS methodology as a driver for informed novel structure generation is also demonstrated. The advantages over other data-driven methods and state-of-the-art methods restricted to high-symmetry systems are highlighted.

atomic cluster expansion↗
DOE OSTI2020

Molecular-Level Exploration of the Structure-Function Relations Underlying Interfacial Charge Transfer in the Subphthalocyanine/C 60 Organic Photovoltaic System

The arrangement of organic molecules at the donor-acceptor interface in an organic photovoltaic (OPV) cell can have a strong effect on the generation of charge carriers and thereby cell performance. In this paper, we report the molecular-level exploration of the ensemble of interfacial donor-acceptor pair geometries and the charge-transfer (CT) rates to which they give rise. Our approach combines molecular-dynamics simulations, electronic structure calculations, machine learning, and rate theory. This approach is applied to the boron subphthalocyanine chloride (donor) and C 60 (acceptor) OPV system. In this work we find that the interface is dominated by a previously unreported donor-acceptor pair edge geometry, which contributes significantly to device performance in a manner that depends on the initial conditions. Quantitative relations between the morphology and CT rates are established, which can be used to advance the design of more efficient OPV devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Nonperturbative phonon scatterings and the two-channel thermal transport in Tl 3 VSe 4

We review the role of nonperturbative phonon scattering in strongly anharmonic materials having ultralow lattice thermal conductivity with unusual temperature dependence. We take Tl 3 VSe 4 as an example and investigate its lattice dynamics using perturbation theory (PT) up to the fourth order and molecular dynamics (MD) with a machine-learning potential. We find distinct differences of phonon linewidth between PT and MD in the whole Brillouin zone. The comparison between the theoretical phonon linewidths and experiments suggests that PT severely underestimates the phonon scatterings, even when the fourth-order anharmonicity is included. Moreover, we extend our calculations to higher temperatures and evaluate the two-channel thermal conductivity based on the unified theory developed by Simoncelli et al. [Nat. Phys. 15, 809 (2019)]. We find a crucial coherence contribution to the total thermal conductivity at high temperatures. Our results pave the path for future studies of phonon properties and lattice thermal conductivities of strongly anharmonic crystals beyond the conventional PT realm.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Correlated dynamic disorder, octahedral tilts, and acoustic phonon softening in CsSnBr 3 and CsPbBr 3

Metal halide perovskites (MHPs) have emerged as highly promising materials for optoelectronic applications, with all-inorganic MHPs presenting enhanced stability compared to their hybrid counterparts. Here, in this study, we investigate the atomic dynamics and structural fluctuations in single crystals of CsSnBr⁢ 3 and CsPbBr⁢ 3 through systematic inelastic neutron scattering (INS) measurements as a function of temperature. Our experiments are compared with first-principle simulations, augmented with large-scale molecular dynamics modeling, based on machine-learned neural network potentials. Through both INS and simulations, we find quasi-elastic diffuse rods in reciprocal space in both compounds, originating from fluctuating planar domains featuring correlated tilts of Br octahedron. The diffuse rods exhibit a slow, overdamped dynamic process, modulated across 𝑸 space, reflecting the strong lattice anharmonicity of the inorganic framework. We do not find evidence for dynamic off-centering of the Sn 2+ ions besides phonon vibrations at the center of the Br octahedron. These results offer valuable insights into the unusual anharmonic atomic dynamics and intricate correlated structural distortions in MHPs, which will be critical for rationalizing and further tailoring their thermal and optoelectronic properties.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Crystal structure of vancomycin bound to the resistance determinant D -alanine- D -serine

Vancomycin is a glycopeptide antibiotic that for decades has been a mainstay of treatment for persistent bacterial infections. However, the spread of antibiotic resistance threatens its continued utility. In particular, vancomycin-resistant enterococci (VRE) have become a pressing clinical challenge. Vancomycin acts by binding and sequestering the intermediate Lipid II in cell-wall biosynthesis, specifically recognizing a D-alanine-D-alanine dipeptide motif within the Lipid II molecule. VRE achieve resistance by remodeling this motif to either D-alanine-D-lactate or D-alanine-D-serine; the former substitution essentially abolishes recognition by vancomycin of Lipid II, whereas the latter reduces the affinity of the antibiotic by roughly one order of magnitude. The complex of vancomycin bound to D-alanine-D-serine has been crystallized, and its 1.20 Å X-ray crystal structure is presented here. This structure reveals that the D-alanine-D-serine ligand is bound in essentially the same position and same pose as the native D-alanine-D-alanine ligand. The serine-containing ligand appears to be slightly too large to be comfortably accommodated in this way, suggesting one possible contribution to the reduced binding affinity. In addition, two flexible hydroxyl groups – one from the serine side chain of the ligand, and the other from a glucose sugar on the antibiotic – are locked into single conformations in the complex, which is likely to contribute an unfavorable entropic component to the recognition of the serine-containing ligand.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Rational Design of High-Performing Electrodes in Energy Storage Devices

Rational design of interfaces with proper physico-chemical properties is necessary to optimize the performance of electrochemical devices, which requires fundamental understanding of the interfaces. In this project, we leveraged quantum mechanical and molecular dynamics simulations as well as machine learning (ML) technique to correlate the relationships between structure and properties including ion transport and electronic structures of cathode host materials in lithium-sulfur (Li-S) and carbon anode in sodium (Na) batteries for improved electrochemical performance.

25 ENERGY STORAGE↗
DOE OSTI2023

Electrochemical Characterization of Molten Salt Chemistry During Atmospheric Ingressions

This report describes the salt chemistry and sedimentation studies that were completed in FY23 as part of the Pyrochemical Fuel Cycles – ANL project. The primary objective of these activities was to employ electrochemical methods in order to quantify and gain insights into the reaction mechanisms associated with the interaction of O2 and moisture impurities with the molten chloride salts used for the pyrochemical processing of nuclear materials. If not monitored and controlled, the ingression of these atmospheric impurities can lead to changes in the salt redox conditions and initiate the formation of oxide particles, potentially resulting in operational and safeguards challenges. By elucidating the mechanisms underlying the formation of oxide particles, we can effectively monitor the conditions of the salt and design systems to mitigate the influence of atmospheric O2 and H2O. To achieve these goals, we designed a test system that allowed for precisely controlled gas ingressions into the molten salt vessel. This system included on-line monitoring provided by electrochemical probes that enable near real-time measurements of the salt conditions throughout the course of the experiment. Using this system, we systematically varied experimental conditions, including O2/Ar flow rates and O2 concentrations, to comprehensively understand their impact on solid particle formation. We also used a particle size and shape analyzer to characterize the particles that were generated. Tests in FY23 concentrated on the formation of CeO2 particles, but we also prepared a separate apparatus targeting UO2 particles for use in FY24. We additionally conducted extensive modeling activities in support of this work. This included multiphysics simulations of the gas ingressions along with the development of a machine learning workflow to enable accelerated molecular dynamics modeling of the molten salt chemistry of LiCl-KCl-UCl3. The combination of experimental and modeling tools has allowed us to begin to get a more complete understanding of the reactions that occur when atmospheric ingressions occur in pyroprocessing systems.

Guo, Jicheng↗
DOE OSTI2022

Data-Driven and Multiscale Modeling of DNA-Templated Dye Aggregates

Dye aggregates are of interest for excitonic applications, including biomedical imaging, organic photovoltaics, and quantum information systems. Dyes with large transition dipole moments (μ) are necessary to optimize coupling within dye aggregates. Extinction coefficients (ε) can be used to determine the μ of dyes, and so dyes with a large ε (>150,000 M –1 cm –1 ) should be engineered or identified. However, dye properties leading to a large ε are not fully understood, and low-throughput methods of dye screening, such as experimental measurements or density functional theory (DFT) calculations, can be time-consuming. In order to screen large datasets of molecules for desirable properties (i.e., large ε and μ), a computational workflow was established using machine learning (ML), DFT, time-dependent (TD-) DFT, and molecular dynamics (MD). ML models were developed through training and validation on a dataset of 8802 dyes using structural features. A Classifier was developed with an accuracy of 97% and a Regressor was constructed with an R2 of above 0.9, comparing between experiment and ML prediction. Using the Regressor, the ε values of over 18,000 dyes were predicted. The top 100 dyes were further screened using DFT and TD-DFT to identify 15 dyes with a μ relative to a reference dye, pentamethine indocyanine dye Cy5. Two benchmark MD simulations were performed on Cy5 and Cy5.5 dimers, and it was found that MD could accurately capture experimental results. The results of this study exhibit that our computational workflow for identifying dyes with a large μ for excitonic applications is effective and can be used as a tool to develop new dyes for excitonic applications.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Auto-deconvolution and molecular networking of gas chromatography–mass spectrometry data

We engineered a machine learning approach, MSHub, to enable auto-deconvolution of gas chromatography–mass spectrometry (GC–MS) data. We then designed workflows to enable the community to store, process, share, annotate, compare and perform molecular networking of GC–MS data within the Global Natural Product Social (GNPS) Molecular Networking analysis platform. MSHub/GNPS performs auto-deconvolution of compound fragmentation patterns via unsupervised non-negative matrix factorization and quantifies the reproducibility of fragmentation patterns across samples.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Role of Local Structure in the Enhanced Dynamics of Deformed Glasses

External stress can accelerate molecular mobility of amorphous solids by several orders of magnitude. The changes in mobility are commonly interpreted through the Eyring model, which invokes an empirical activation volume. Here, we analyze constant-stress molecular dynamics simulations and propose a structure-dependent Eyring model, connecting activation volume to a machine-learned field, softness. We show that stress has a heterogeneous effect on the mobility that depends on local structure through softness. The barrier impeding relaxation reduces more for well-packed particles, which explains the narrower distribution of relaxation time observed under stress.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Machine learning enabled quantification of the hydrogen bonds inside the polyelectrolyte brush layer probed using all-atom molecular dynamics simulations

The configuration of densely grafted charged polyelectrolyte (PE) brushes is strongly dictated by the properties and behavior of the counterions that screen the PE brush charges and the solvent molecules (typically water) that solvate the brush molecules and these screening counterions. Only recently, efforts have been made to study the PE brushes atomistically, thereby shedding light on the properties of brush-supported ions and water molecules. However, even for such efforts, there are limitations associated with using a generic definition to estimate certain properties of water and ions inside the brush layer. For example, water–water hydrogen bonds (HBs) will behave differently for locations outside and inside the brush layer, given the fact that the densely closely grafted PE brush molecules create a soft nanoconfinement where the water connectivity becomes highly disrupted: therefore, using the same definition to quantify the HBs inside and outside the brush layer will be unwise. In this paper, we address this limitation by employing an unsupervised machine learning (ML) approach to predict the water–water hydrogen bonding inside a cationic PE brush layer modeled using all-atom molecular dynamics (MD) simulations. Here, the ML method, which relies on a clustering approach and uses the equilibrium coordinates of the water molecules (obtained from the all-atom MD simulations) as the input, is capable of identifying the structural modification of water–water HBs (revealed through appropriate clustering of the data) inside the PE brush layer induced soft nanoconfinement. Such capabilities would not have been possible by using a generic definition of the HBs. Our calculations lead to four key findings: (1) the clusters formed inside and outside the brush layer are structurally similar; (2) the margin of the cluster is shorter inside the PE brush layer confirming the possible disruption of the HBs inside the PE brush layer; (3) the average “hydrogen–acceptor-oxygen–donor-oxygen” angle that defines the HB is reduced for the HBs formed inside the brush layer; (4) the use of the generic definition (definition usable for characterizing the HBs in brush-free bulk) leads to an overprediction of the number of HBs formed inside the PE brush layer.

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

Equivariant Graph Attention Network - 3D Conformers & Feature Fusion

EGAN-3F (Equivariant Graph Attention Network - 3D Conformers & Feature Fusion) presents an innovative approach for predicting binding affinity between small molecules and protein targets, a fundamental task in drug discovery. Traditional structure-based methods often depend on protein-ligand complex structures obtained from crystallography or molecular docking. In contrast, ligand-only machine learning models using 1D or 2D representations such as SMILES have been developed to predict binding affinity without structural information about the target; however, their accuracy is often limited due to the lack of 3D ligand information. EGAN-3F addresses this limitation by integrating spatially aware graph learning with traditional descriptor-based features. We systematically investigate how combining 2D and 3D molecular representations enhances binding affinity prediction from SMILES strings. This approach underscores the importance of modeling conformational diversity and incorporating chemically meaningful descriptors to improve predictive accuracy. The key innovation of EGAN-3F lies in its ability to achieve robust ligand-based binding affinity predictions without requiring protein-ligand complex structures, effectively bridging the gap between purely structural and ligand-only modeling paradigms.

Shim, Heesung [Lawrence Livermore National Laborat↗
DOE OSTI2024

Machine learning assisted prediction of tungsten heavy alloy plasma facing component performance for fusion energy applications

Tungsten and tungsten heavy alloys (WHAs), known for their remarkably high hardness, durability, and corrosion resistance, play a critical role in the thriving development of nuclear fusion reactors in recent years. However, the exploration in tungsten alloys for the nuclear-related applications has been limited by the difficulty of manufacturing and the complexity of experiments to reproduce the environment of nuclear reaction. Therefore, this project aims to utilize nanoscale simulation methods such as density functional theory (DFT) and molecular dynamics (MD) with the help of machine learning techniques to not only understand the mechanisms of tungsten alloys but also allow us to computationally predict their mechanical behaviors under extreme environments. One critical problem of the application of WHAs in nuclear reactors is the surface melting. In the current design of the SPARC reactor, the WHA, W97Ni2.1Fe0.9 or W97NiFe, is chosen to be the first wall components to confine the plasma where the particles are fiercely moving and colliding into each other to create nuclear fusion reaction. This process will generate extremely high heat flux onto these WHA tiles, leaving high surface temperature that could possibly melt the surface of the WHA tiles, As illustrated in Fig. 1(a). a laser experiment previously done illustrates that a rough surface damage would be made after the surface melting where the matrix area mainly composed of nickel and iron as shown in Fig. 1(b), will first melt and then leave vacancies between these tungsten grains. Unfortunately, these kinds of roughness on the first-wall components could be deadly to the plasma inside a Tokmak reactor because the heat that is supposed to dissipate at a designed ratio through the tiles may in turn be excessively absorbed and accumulated on any uneven area of the surface, which will eventually make the whole nuclear reaction fail. In this project, we will introduce a machine learning potential, Allegro, based on DFT calculation and then build a MD model for W-Ni-Fe alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Microscopic insights into the solvation of polyethylene glycol chains in water: A machine learning potential approach

Polyethylene glycol (PEG) is a structurally simple, nontoxic, and water-soluble polymer widely utilized in medical and pharmaceutical applications. Notably, when a PEG chain is immersed in water, the surrounding water molecules play a key role in driving conformational changes of this macromolecule. In this study, we explore the solvation behavior of PEG under mechanical strain using molecular dynamics simulations, with an interatomic potential obtained from machine learning. Our focus is on the transition from the favored coil-like conformation to an extended one under external force. Through analyses of radial distribution functions, hydrogen bonding, and solvation dynamics, we uncover how mechanical stretching influences the local hydration environment. Furthermore, we disentangle the enthalpic and entropic contributions to the conformational stability of PEG in water. Surprisingly, our neural network potential model identifies dewetting of PEG C-atoms, and not water H-bonding with PEG O-atoms, as the main enthalpic driving force for the coiling of PEG in water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗