Molecular Structure and Modeling of Water–Air and Ice–Air Interfaces Monitored by Sum-Frequency Generation
Not Available
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
SEARCH · Engineering Papers
Results for “molecular modelling”
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.
A data-driven peridynamic continuum model for upscaling molecular dynamics.
Abstract not provided.
D'Elia, Marta↗
Autonomous alloy composition optimization using molecular dynamics guided by a large language model
Here, we present an autonomous materials discovery framework that couples a large language model (LLM) with molecular dynamics (MD) simulations to optimize Fe–Cr–Mn alloy compositions for tensile strength. Starting from six distinct compositions, the LLM operated as an intelligent agent, iteratively proposing changes based on prior simulation results and constraints. Over 50 iterations per case, the LLM adaptively explored the composition space, identifying high-strength regions, not easily accessible by conventional methods. The highest strength, 18.7 GPa, was achieved with Fe 71 Cr 25 Mn 4 composition, identified from a Fe 75 Cr 20 Mn 5 starting point. The LLM autonomously adjusted its strategy in real time, demonstrating closed-loop decision-making using commodity hardware. This approach showcases the potential of LLMs as scientific co-pilots, capable of accelerating materials discovery and generalizable to other domains like biology and drug design.
Autonomy↗
RLMolLM: Reinforcement Learning-Enhanced Language Model Framework for Inverse Molecular Design
Inverse molecular design faces significant challenges due to vast chemical space and complex property requirements. While language models show promise for molecular generation, they struggle with validity, multi-property optimization, and structural constraints. This work presents RLMolLM, a reinforcement learning framework combining Proximal Policy Optimization (PPO) with genetic algorithms to address these limitations. Our approach optimizes multiple user-specified properties including quantitative estimates of drug-likeness (QED), synthetic accessibility (SA), and ADMET (absorption, distribution, metabolism, excretion, and toxicity) endpoints without requiring complete model retraining, while maintaining capability for scaffold-constrained generation where specific substructures must be preserved. We outperform state-of-the-art methods for molecular optimization, achieving best QED scores across GDB13, Moses, and Zinc datasets with up to 31% improvement over previous methods while maintaining excellent validity, uniqueness, and novelty metrics. For simultaneous multi-property optimization, our framework achieves substantial improvements in ADMET properties including 4.5-fold reduction in hERG toxicity and enhanced Caco-2 permeability compared to Moses dataset. Under structural constraints, the framework significantly improves molecular validity while preserving scaffolds and effectively optimizing properties. In conclusion, this versatile solution advances pharmaceutical and materials molecular design through effective integration of reinforcement learning and genetic algorithms with multi-property optimization and scaffold preservation.
Genetic algorithms↗
Resolving Dynamics in the Ensemble: Finding Paths through Intermediate States and Disordered Protein Structures
Proteins have been found to inhabit a diverse set of three-dimensional structures. The dynamics that govern protein interconversion between structures happen over a wide range of time scales-picoseconds to seconds. Our understanding of protein functions and dynamics is largely reliant upon our ability to elucidate physically populated structures. From an experimental structural characterization perspective, we are often limited to measuring the ensemble-averaged structure both in the steady-state and time-resolved regimes. Generating kinetic models and understanding protein structure–function relationships require atomistic knowledge of the populated states in the ensemble. Here, we present ensemble refinement methodologies that integrate time-resolved experimental signals with molecular dynamics models. We first discuss integration of experimental structural restraints to molecular models in disordered protein systems that adhere to the principle of maximum entropy for creating a complete set of ensemble structures. We then propose strategies to find kinetic pathways between the refined structures, using time-resolved inputs to guide molecular dynamics trajectories and the use of inference to generate tailored stimuli to prepare a desired ensemble of protein states.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Experimental investigation of abnormal transverse flow enhancement of α particles in heavy-ion collisions
The mass dependence of the transverse flow for Z = 1–5 fragments from the collisions of 40 Ar + 27 Al, 40 Ar + 48 Ti, and 40 Ar + 58 Ni at 47 MeV/nucleon is investigated experimentally in this article. The transverse flow values are determined using the in-plane components of the fragment transverse momenta, where three conventional methods, i.e., the kinetic flow tensor method, the transverse momentum analysis method, and the azimuthal correlation method, are applied to reconstruct the reaction plane in an event-by-event basis. It is demonstrated from the comparison of the present experimental mass dependent flow measurements and the model simulations using an improved antisymmetrized molecular dynamics model that the experimentally observed abnormal α transverse flow enhancement is closely related to the reaction plane reconstruction procedure in the flow extraction. We further investigate the physical existence of the abnormal α flow behavior using a two-particle azimuthal correlation method, which allows us to provide the relative flow magnitude information with an identification of fragment charge number without the knowledge of the reaction plane differing from the three conventional methods. It is found that the relative flow magnitudes deduced from the two-particle azimuthal correlation functions with an identification of Z, with the correction for the recoil effect imposed by the momentum conservation, show a monotonically increasing trend as a function of fragment charge number, with no exception of the α flow enhancement. These results, in addition to those from the improved antisymmetrized molecular dynamics model simulations, definitely provide experimental evidences for the inexistence of the abnormal α flow behavior in the heavy-ion collisions at the present incident energy region in nature.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Quantum Molecular Charge-Transfer Model for Multistep Auger–Meitner Decay Cascade Dynamics
The fragmentation of molecular cations following inner-shell decay processes in molecules containing heavy elements underpins the X-ray damage effects observed in X-ray scattering measurements of biological and chemical materials, as well as in medical applications involving Auger electron-emitting radionuclides. Traditionally, these processes are modeled using simulations that describe the electronic structure at an atomic level, thereby omitting molecular bonding effects. This work addresses the gap by introducing a novel approach that couples Auger–Meitner decay to nuclear dynamics across multiple decay steps, by developing a decay spawning dynamics algorithm and applying it to potential energy surfaces characterized with ab initio molecular dynamics simulations. We showcase the approach on a model decay cascade following K-shell ionization of IBr and subsequent Kβ fluorescence decay. We examine two competing channels that undergo two decay steps, resulting in ion pairs with a total 3+ charge state. This approach provides a continuous description of the electron transfer dynamics occurring during the multistep decay cascade and molecular fragmentation, revealing the combined inner-shell decay and charge transfer time scale to be approximately 75 fs. In conclusion, our computed kinetic energies of ion fragments show good agreement with experimental data.
Ab initio molecular dynamics↗
Short-wavelength-sensitive 2 (Sws2) visual photopigment models combined with atomistic molecular simulations to predict spectral peaks of absorbance
For many species, vision is one of the most important sensory modalities for mediating essential tasks that include navigation, predation and foraging, predator avoidance, and numerous social behaviors. The vertebrate visual process begins when photons of the light interact with rod and cone photoreceptors that are present in the neural retina. Vertebrate visual photopigments are housed within these photoreceptor cells and are sensitive to a wide range of wavelengths that peak within the light spectrum, the latter of which is a function of the type of chromophore used and how it interacts with specific amino acid residues found within the opsin protein sequence. Minor differences in the amino acid sequences of the opsins are known to lead to large differences in the spectral peak of absorbance (i.e. the λ max value). In our prior studies, we developed a new approach that combined homology modeling and molecular dynamics simulations to gather structural information associated with chromophore conformation, then used it to generate statistical models for the accurate prediction of λ max values for photopigments derived from Rh1 and Rh2 amino acid sequences. In the present study, we test our novel approach to predict the λ max of phylogenetically distant Sws2 cone opsins. To build a model that can predict the λ max using our approach presented in our prior studies, we selected a spectrally-diverse set of 11 teleost Sws2 photopigments for which both amino acid sequence information and experimentally measured λ max values are known. The final first-order regression model, consisting of three terms associated with chromophore conformation, was sufficient to predict the λ max of Sws2 photopigments with high accuracy. This study further highlights the breadth of our approach in reliably predicting λ max values of Sws2 cone photopigments, evolutionary-more distant from template bovine RH1, and provided mechanistic insights into the role of known spectral tuning sites.
59 BASIC BIOLOGICAL SCIENCES↗
The impact of the atmospheric turbulence-development tendency on new particle formation: a common finding on three continents
Abstract A new mechanism of new particle formation (NPF) is investigated using comprehensive measurements of aerosol physicochemical quantities and meteorological variables made in three continents, including Beijing, China; the Southern Great Plains site in the USA; and SMEAR II Station in Hyytiälä, Finland. Despite the considerably different emissions of chemical species among the sites, a common relationship was found between the characteristics of NPF and the stability intensity. The stability parameter (ζ = Z/L, where Z is the height above ground and L is the Monin–Obukhov length) is found to play an important role; it drops significantly before NPF as the atmosphere becomes more unstable, which may serve as an indicator of nucleation bursts. As the atmosphere becomes unstable, the NPF duration is closely related to the tendency for turbulence development, which influences the evolution of the condensation sink. Presumably, the unstable atmosphere may dilute pre-existing particles, effectively reducing the condensation sink, especially at coarse mode to foster nucleation. This new mechanism is confirmed by model simulations using a molecular dynamic model that mimics the impact of turbulence development on nucleation by inducing and intensifying homogeneous nucleation events.
54 ENVIRONMENTAL SCIENCES↗
Atomic-Level Features for Kinetic Monte Carlo Models of Complex Chemistry from Molecular Dynamics Simulations
The high computational cost of evaluating atomic interactions recently motivated the development of computationally inexpensive kinetic models, which can be parametrized from MD simulations of complex chemistry of thousands of species or other processes and accelerate the prediction of the chemical evolution by up to four order of magnitude. Such models go beyond the commonly employed potential energy surface fitting methods in that they are aimed purely at describing kinetic effects. So far, such kinetic models utilize molecular descriptions of reactions and have been constrained to only reproduce molecules previously observed in MD simulations. Therefore, these descriptions fail to predict the reactivity of unobserved molecules, for example in the case of large molecules or solids. In this work, we propose a new approach for the extraction of reaction mechanisms and reaction rates from MD simulations, namely the use of atomic-level features. Using the complex chemical network of hydrocarbon pyrolysis as example, it is demonstrated that kinetic models built using atomic features are able to explore chemical reaction pathways never observed in the MD simulations used to parametrize them, a critical feature to describe rare events. Atomic-level features are shown to construct reaction mechanisms and estimate reaction rates of unknown molecular species from elementary atomic events. Through comparisons of the model ability to extrapolate to longer simulation timescales and different chemical compositions than the ones used for parameterization, it is demonstrated that kinetic models employing atomic features retain the same level of accuracy and transferability as the use of features based on molecular species, while being more compact and parametrized with less data. We also find that atomic features can better describe the formation of large molecules enabling the simultaneous description of small molecules and condensed phases.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Advancing quantitative models of soil microbiology, ecology, and biochemistry
Rapid advances in soil biology are increasingly reflected in mathematical models. From molecular to global scales, models contribute to fundamental understanding and prediction of critical soil processes, such as carbon sequestration. Recent models explicitly incorporate microbial control over soil carbon and nutrient cycling, an approach that offers greater biological realism but raises new challenges in model parameterization and stability. Although they have proliferated at community and ecosystem scales, microbial-explicit models have yet to be incorporated into most Earth system models. Applying soil models across scales requires integration with data to support parameterization and validation using approaches such as Bayesian data assimilation. In conclusion, to catalyze future progress, we recommend breaking down barriers between modeling and empirical disciplines while broadening access to the computational infrastructure that supports soil biological modeling.
Xie, Wally↗
Multinucleon transfer in 48 Ti + 48 Ti collisions at 11.5 MeV/nucleon
In this work, experimental data of projectile fragment distributions from a preliminary experimental study of the reaction of 11.5 MeV/nucleon 48 Ti on 48 Ti analyzed with the MARS recoil separator at the Cyclotron Institute of Texas A&M University are presented. Production cross sections, momentum distributions and excitation-energy distributions are extracted and compared with calculations with the Deep Inelastic Transfer (DIT) model and the Constrained Molecular Dynamics model (CoMD). Despite the limited extent of the present dataset, an overall agreement of the models with the data is found pointing at the prevailing nucleon-exchange character of the collisions at this energy.
Souliotis, Georgios↗
KCl-UCl 3 molten salts investigated by Ab Initio Molecular Dynamics (AIMD) simulations: A comparative study with three dispersion models
Ab Initio Molecular Dynamics (AIMD) simulations are performed on molten KCl-UCl 3 salt mixtures to determine energies, heat capacities, and densities. The density-dependent energy correction (DFT-dDsC), Grimme et al.’s DFT-D3, and Langreth & Lundqvist (vdW-cx) models are used for dispersion forces and combined with the Perdew-Burke-Ernzerhof (PBE) exchange-correlation potential with a Hubbard U param eter for the 5$f$ electrons of uranium. After validating predictions for the end-member systems to literature data, KCl-UCl 3 mixtures are studied at select temperatures. Densities and energies both deviate from ideal solution behavior, with the maximum deviation occurring around 36% UCl 3 for mixing energies and slightly lower (29% UCl 3 ) for densities. Compared to the NaCl-UCl 3 system, which was previously investigated using the same simulation methodologies, the KCl-UCl 3 density and mixing energy deviations from ideal solution behavior are larger by almost a factor of two. No deviation from ideal solution behavior for heat capacity was observed. The AIMD predictions for mixing energies and densities agree qualitatively with experimental data, though the spread in data obtained from the various dispersion force models utilized, measurements, and empirical estimates makes strong conclusions difficult. The dependence of thermodynamic and thermophysical properties on composition is correlated with the local chemistry of the solution phase, in particular, the tendency of UCl 3 to form network structures.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Geometry-complete diffusion for 3D molecule generation and optimization
Abstract Generative deep learning methods have recently been proposed for generating 3D molecules using equivariant graph neural networks (GNNs) within a denoising diffusion framework. However, such methods are unable to learn important geometric properties of 3D molecules, as they adopt molecule-agnostic and non-geometric GNNs as their 3D graph denoising networks, which notably hinders their ability to generate valid large 3D molecules. In this work, we address these gaps by introducing the Geometry-Complete Diffusion Model (GCDM) for 3D molecule generation, which outperforms existing 3D molecular diffusion models by significant margins across conditional and unconditional settings for the QM9 dataset and the larger GEOM-Drugs dataset, respectively. Importantly, we demonstrate that GCDM’s generative denoising process enables the model to generate a significant proportion of valid and energetically-stable large molecules at the scale of GEOM-Drugs, whereas previous methods fail to do so with the features they learn. Additionally, we show that extensions of GCDM can not only effectively design 3D molecules for specific protein pockets but can be repurposed to consistently optimize the geometry and chemical composition of existing 3D molecules for molecular stability and property specificity, demonstrating new versatility of molecular diffusion models. Code and data are freely available on GitHub .
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Shadow Molecular Dynamics and Atomic Cluster Expansions for Flexible Charge Models
Here, a shadow molecular dynamics scheme for flexible charge models is presented where the shadow Born–Oppenheimer potential is derived from a coarse-grained approximation of range-separated density functional theory. The interatomic potential, including the atomic electronegativities and the charge-independent short-range part of the potential and force terms, is modeled by the linear atomic cluster expansion (ACE), which provides a computationally efficient alternative to many machine learning methods. The shadow molecular dynamics scheme is based on extended Lagrangian (XL) Born–Oppenheimer molecular dynamics (BOMD). XL-BOMD provides stable dynamics while avoiding the costly computational overhead associated with solving an all-to-all system of equations, which normally is required to determine the relaxed electronic ground state prior to each force evaluation. To demonstrate the proposed shadow molecular dynamics scheme for flexible charge models using atomic cluster expansion, we emulate the dynamics generated from self-consistent charge density functional tight-binding (SCC-DFTB) theory using a second-order charge equilibration (QEq) model. The charge-independent potentials and electronegativities of the QEq model are trained for a supercell of uranium oxide (UO 2 ) and a molecular system of liquid water. The combined ACE+XL-QEq molecular dynamics simulations are stable over a wide range of temperatures both for the oxide and for the molecular systems and provide a precise sampling of the Born–Oppenheimer potential energy surfaces. Accurate ground Coulomb energies are produced by the ACE-based electronegativity model during an NVE simulation of UO 2 , predicted to be within 1 meV of those from SCC-DFTB on average during comparable simulations.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Stochastic symplectic reduced-order modeling for model-form uncertainty quantification in molecular dynamics simulations in various statistical ensembles
Here, this work focuses on the representation of model-form uncertainties in molecular dynamics simulations in various statistical ensembles. In prior contributions, the modeling of such uncertainties was formalized and applied to quantify the impact of, and the error generated by, pair-potential selection in the microcanonical ensemble (NVE). In this work, we extend this formulation and present a linear-subspace reduced-order model for the canonical (NVT) and isobaric (NPT) ensembles. The symplectic reduced-order basis is randomized on the tangent space of the Stiefel manifold to provide topological relationships and capture model-form uncertainty. Using the Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS), we assess the relevance of these stochastic reduced-order atomistic models on canonical problems involving a Lennard-Jones fluid and an argon crystal melt.
42 ENGINEERING↗
Attention-based generative models for de novo molecular design
An implementation of attention within the variational autoencoder framework for continuous representation of molecules. The addition of attention significantly increases model performance for complex tasks such as exploration of novel chemistries.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
AIMSim : An accessible cheminformatics platform for similarity operations on chemicals datasets
The recent advances in deep learning, generative modeling, and statistical learning have ushered in a renewed interest in traditional cheminformatics tools and methods. Quantifying molecular similarity is essential in molecular generative modeling, exploratory molecular synthesis campaigns, and drug-discovery applications to assess how new molecules differ from existing ones. Further, most tools target advanced users and lack general implementations accessible to the larger community. In this work, we introduce Artificial Intelligence Molecular Similarity (AIMSim), an accessible cheminformatics platform for performing similarity operations on collections of molecules called molecular datasets. AIMSim provides a unified platform to perform similarity-based tasks on molecular datasets, such as diversity quantification, outlier and novelty analysis, clustering, dimensionality reduction, and inter-molecular comparisons. AIMSim implements all major binary similarity metrics and molecular fingerprints and is provided as a Python package that includes support for command-line use as well as a Graphical User Interface for code-free utilization with fully interactive plots.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗