Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “accelerated molecular dynamics”

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

Mechanistic Insights into the Demethylation of Lignin‐Derived Structures Using Protic Ionic Liquids: A Density Functional Theory Study

Lignin valorization is restricted by the stability of its methoxy groups, creating a critical need for efficient demethylation strategies. Here, in this study, density functional theory (DFT) is employed to dissect the mechanistic pathways of demethylation in lignin model compounds, guaiacol and syringol, using protic ionic liquids (PILs) that act as both solvent and catalyst. Conductor‐like Screening Model for Real Solvents (COSMO‐RS) analysis identifies monoethanolammonium acetate ([MEOA][Ace]) as the most promising medium, attributed to its strong hydrogen bonding network and solvation ability. By integrating implicit and explicit solvation models, it is revealed that an acid‐catalyzed hydrolytic mechanism governs demethylation, with PILs stabilizing crucial transition states and intermediates. Complementary electronic structure evaluations, including highest occupied molecular orbital–lowest unoccupied molecular orbital (HOMO‐LUMO) gap analysis, charge distribution, and electrostatic potential mapping, demonstrate how PILs lower energetic barriers and enhance reactivity. To simulate realistic environments, this study is extended to lignin dimer complexes with varying water content, uncovering how water‐bridged solvation reverses demethylation preference from guaiacyl to syringyl units. This mechanistic shift aligns with experimental observations showing faster S‐unit reactivity in hydrated systems. Together, these findings provide atomic‐level insight into lignin demethylation dynamics and highlight how tuning acid concentration in PILs can accelerate kinetics, enabling a rational pathway toward next‐generation biomass conversion technologies.

COSMO-RS↗
DOE OSTI2021

Semi-Automated Creation of Density Functional Tight Binding Models through Leveraging Chebyshev Polynomial-Based Force Fields

Density functional tight binding (DFTB) is an attractive method for accelerated quantum simulations of condensed matter due to its enhanced computational efficiency over standard density functional theory (DFT) approaches. However, DFTB models can be challenging to determine for individual systems of interest, especially for metallic and interfacial systems where different bonding arrangements can lead to significant changes in electronic states. In this regard, we have created a rapid-screening approach for determining systematically improvable DFTB interaction potentials that can yield transferable models for a variety of conditions. Our method leverages a recent reactive molecular dynamics force field where many-body interactions are represented by linear combinations of Chebyshev polynomials. This allows for the efficient creation of multi-center representations with relative ease, requiring only a small investment in initial DFT calculations. Here, we have focused our workflow on TiH 2 as a model system and show that a relatively small training set based on unit-cell-sized calculations yields a model accurate for both bulk and surface properties. Our approach is easy to implement and can yield reliable DFTB models over a broad range of thermodynamic conditions, where physical and chemical properties can be difficult to interrogate directly and there is historically a significant reliance on theoretical approaches for interpretation and validation of experimental results.

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

Evaluation of 6-OxP-CD, an Oxime-based cyclodextrin as a viable medical countermeasure against nerve agent poisoning: Experimental and molecular dynamic simulation studies on its inclusion complexes with cyclosarin, soman and VX

The ability of the cyclodextrin-oxime construct 6-OxP-CD to bind and degrade the nerve agents Cyclosarin (GF), Soman (GD) and S -[2-[Di(propan-2-yl)amino]ethyl] O -ethyl methylphosphonothioate (VX) has been studied using 31 P-nuclear magnetic resonance (NMR) under physiological conditions. While 6-OxP-CD was found to degrade GF instantaneously under these conditions, it was found to form an inclusion complex with GD and significantly improve its degradation (t 1/2 ~ 2 hrs) relative over background (t 1/2 ~ 22 hrs). Consequently, effective formation of the 6-OxP-CD:GD inclusion complex results in the immediate neutralization of GD and thus preventing it from inhibiting its biological target. In contrast, NMR experiments did not find evidence for an inclusion complex between 6-OxP-CD and VX, and the agent’s degradation profile was identical to that of background degradation (t 1/2 ~ 24 hrs). As a complement to this experimental work, molecular dynamics (MD) simulations coupled with Molecular Mechanics-Generalized Born Surface Area (MM-GBSA) calculations have been applied to the study of inclusion complexes between 6-OxP-CD and the three nerve agents. These studies provide data that informs the understanding of the different degradative interactions exhibited by 6-OxP-CD with each nerve agent as it is introduced in the CD cavity in two different orientations (up and down). For its complex with GF, it was found that the oxime in 6-OxP-CD lies in very close proximity (P GF …O Oxime ~ 4–5 Å) to the phosphorus center of GF in the ‘down GF ’ orientation for most of the simulation accurately describing the ability of 6-OxP-CD to degrade this nerve agent rapidly and efficiently. Further computational studies involving the center of masses (COMs) for both components (GF and 6-OxP-CD) also provided some insight on the nature of this inclusion complex. Distances between the COMs (ΔCOM) lie closer in space in the ‘down GF ’ orientation than in the ‘up GF ’ orientation; a correlation that seems to hold true not only for GF but also for its congener, GD. In the case of GD, calculations for the ‘down GD ’ orientation showed that the oxime functional group in 6-OxP-CD although lying in close proximity (P GD …O Oxime ~ 4–5 Å) to the phosphorus center of the nerve agent for most of the simulation, adopts another stable conformation that increase this distance to ~ 12–14 Å, thus explaining the ability of 6-OxP-CD to bind and degrade GD but with less efficiency as observed experimentally (t 1/2 ~ 4 hr. vs. immediate). Lastly, studies on the VX:6-OxP-CD system demonstrated that VX does not form a stable inclusion complex with the oxime-bearing cyclodextrin and as such does not interact in a way that is conducive to an accelerated degradation scenario. Collectively, these studies serve as a basic platform from which the development of new cyclodextrin scaffolds based on 6-OxP-CD can be designed in the development of medical countermeasures against these highly toxic chemical warfare agents.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Repetitive proteins that undergo large conformational changes evade structural prediction algorithms

Protein structure prediction algorithms, such as AlphaFold, have accelerated protein design and advanced the understanding of the relationship between amino acid sequence and protein structure. However, these algorithms are limited in their ability to predict the structures of conformationally dynamic, intrinsically disordered, and stimuli-responsive proteins. To evaluate sequence-to-structure predictions of such challenging proteins, we explored a class of conformationally dynamic, repeats-in-toxin (RTX) proteins. RTX proteins adopt intrinsically disordered conformations in the absence of calcium and undergo reversible folding into β-roll structures upon binding to calcium. RTX proteins are characterized by tandem repeats of the sequence GGXGXDXUX, in which X can be any amino acid and U is an aliphatic amino acid. We designed RTX sequence variants with global substitutions of nonconserved amino acids, tandem repeats of consensus sequences GGAGXDTLY, and tandem repeats of scrambled sequences GGAGXDTYL. AlphaFold2 and AlphaFold3 predicted that all of these RTX variants adopt β-roll structures, characteristic of wild-type RTX bound to calcium. However, modeling the predicted structures with molecular dynamics simulations and characterizing the protein variants with circular dichroism spectroscopy, small-angle x-ray scattering, and x-ray crystallography revealed that variants adopt diverse, sequence-dependent structures in the absence and presence of calcium. To better design proteins for applications in biotechnology and sustainability, it is critical to build predictive tools that consider intrinsically disordered protein states and validate these tools with multi-mode, multi-scale experimental data.

Chang, Marina P. [Stanford Univ., CA (United State↗
DOE OSTI2021

Robust Molecular Predictive Methods for Novel Polymer Discovery and Applications

Polymeric materials are ubiquitous in modern society and they play an instrumental role in almost all industries, undoubtedly including the energy and environment sectors. Increased demand of energy and awareness to sustainability both necessitates the development of novel polymers with enhanced properties. Unfortunately, their structural and behavioral complexity render such discovery challenging and impeded. To address this problem, scientists are developing various computational modeling techniques and leveraging their power to depict the relationship between structural characteristics of polymers and their properties (such as rheological behaviors), and use such prediction to guide the design and syntheses of novel polymeric materials with enhanced performances. Unfortunately, predicting the relationships between polymer structure and composition with rheological properties via atomistic modeling is still a major challenge because of the extended time and length scales involved. Studying dynamic shear viscosity and linear viscoelasticity using molecular models requires capabilities that have been elusive, including representation of large molecular weight chains with an effective internal scale capable of describing entanglement, shear-rates that are in the s-1 scale with accurate quantitative stresses, and chemically-realistic combinations of both homogeneous and heterogeneous systems. Motivated by these unmet challenges, the overall technical objective of this DOE-STTR Phase II project is to develop robust molecular predictive methods for advanced polymer discovery and applications and especially for designing and demonstrating the “smart” polymer-based waterflooding enhanced oil recovery (EOR) process. In particular, we apply state-of-the-art molecular modeling methods developed by our academic partner, Materials Stimulation Center (MSC) at California Institute of Technology (Caltech), to facilitate and accelerate the experimental discovery processes. During the Phase I of this project, we had focused on development and demonstration of the molecular modeling methods to describe rheological properties of non-Newtonian polymer fluids, and to improve our fundamental understandings of shear-thickening mechanism and kinetics. In Phase II, we further apply the theoretical models to guide our experimental programs to improve our design of smart rheology modifier (SRM) polymers and their optimization for EOR. Specifically, we have three objectives in the Phase II study: (1) to further improve out computational modeling methods, coupling with the advanced machine learning algorithms; (2) to develop cost-effective and efficient SRM-flooding process suitable for EOR applications under typical reservoir conditions; and (3) to further explore the application of our molecular predictive models for innovative material discovery in other industrial applications. The recent development of our multiscale predictive framework allows the successful prediction of rheological properties from the chemical structure for polymers of experimentally relevant molecular weights, and provides an in-silico machine learning engine for screening novel compositions and structures with optimized non-Newtonian response, required for both shear-thinning and shear-thickening applications. Our framework provides: (1) procedures and tools for systematic coarsening from atomistic models and reverse mapping of coarse-grain models to atomistic, (2) unique ab initio methods to characterize the atomistic origin of colloidal and interfacial interactions and phenomena, (3) systematic structure and composition builders based on practical descriptors that drive rheological changes in polymer melts and diluted polymer mixtures, (4) a rheological properties engine capable of predicting viscosity in the zero-shear limit and under realistic dynamic conditions (for shear-rates commensurate with experiments) for large heterogeneous systems, (5) coarse-grain force fields with improved non-bond descriptions based on accurate quantum mechanics, (6) an in-silico screening machine learning engine that feeds from the systematic model builders to cover the descriptors search space, computes the rheological properties from converged trajectories spanning sub-milliseconds and ranks them for each structure/composition using an automated viscosity-vs-shear rate fitness function that can be tuned for shear-thickening, shear-thinning and other rheological responses.

02 PETROLEUM↗
NASA NTRS1985

Radiative processes

Solar radiation and the processes that control its deposition in the Earth atmosphere are considered. The published data obtained since 1978 define a reference solar spectral irradiance for use in atmospheric chemical and dynamical studies, while long term satellite measurements are now providing information on variations in the output of the Sun over a range of time scales. As concerns absorption of solar radiation in the atmosphere, new cross section data for molecular oxygen and ozone are now available. Line-by-line calculations used to predict infrared flux divergences, both as regards assumptions made in radiative transfer calculations and in the spectroscopic parameters used as inputs are examined. Also examined are the influence of radiative processes on planetary scale wave activity, photochemical acceleration of radiative damping, and the breakdown of local thermodynamic equilibrium at mesospheric altitudes.

Frederick, J. E.↗
DOE OSTI2022

Development and application of ReaxFF methodology for understanding the chemical dynamics of metal carbonates in aqueous solutions

A new ReaxFF reactive force field has been developed for metal carbonate systems including Na + , Ca 2+ , and Mg 2+ cations and the CO 3 2– anion. This force field is fully transferable with previous ReaxFF water and water/electrolyte descriptions. The Me–O–C (Me = metal) three-body valence angle parameters and Me–C non-reactive parameters of the force field have been optimized against quantum mechanical calculations including equations of state, heats of formation, heats of reaction, angle distortions and vibrational frequencies. The new metal carbonate force field has been validated using molecular dynamics simulations to study the solvation and reactivity of metal and carbonate ions in water at 300 K and 700 K. The coordination radius and self-diffusion coefficient show good consistency with existing experimental and simulation results. The angular distribution analysis explains the structural preference of carbonate ions to form carbonates and bicarbonates, where Na + predominantly forms carbonates due to weaker angular strain, while Ca 2+ and Mg 2+ prefer to form bicarbonate monodentate in nature. Residence time distribution analyses on different systems reveal the role of ions in accelerating and decelerating the dynamics of water and carbonate ions under different thermodynamic conditions. The formation and dissolution of bicarbonates and carbonates in solution were explored on the basis of the protonation capability in different systems. Here, the nucleation phenomenon of metal carbonates at ambient and supercritical conditions is explained from the perspective of cluster formation over time: Ca 2+ ions can form prenucleation clusters at ambient temperature but show saturation with increasing temperature, whereas Na + and Mg 2+ ions show a rapid increase in cluster size and amount upon increasing time and temperature.

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

TurboRVB: A many-body toolkit for ab initio electronic simulations by quantum Monte Carlo

TurboRVB is a computational package for ab initio Quantum Monte Carlo (QMC) simulations of both molecular and bulk electronic systems. The code implements two types of well established QMC algorithms: Variational Monte Carlo (VMC) and diffusion Monte Carlo in its robust and efficient lattice regularized variant. A key feature of the code is the possibility of using strongly correlated many-body wave functions (WFs), capable of describing several materials with very high accuracy, even when standard mean-field approaches [e.g., density functional theory (DFT)] fail. The electronic WF is obtained by applying a Jastrow factor, which takes into account dynamical correlations, to the most general mean-field ground state, written either as an antisymmetrized geminal power with spin-singlet pairing or as a Pfaffian, including both singlet and triplet correlations. This WF can be viewed as an efficient implementation of the so-called resonating valence bond (RVB) Ansatz, first proposed by Pauling and Anderson in quantum chemistry [L. Pauling, The Nature of the Chemical Bond (Cornell University Press, 1960)] and condensed matter physics [P.W. Anderson, Mat. Res. Bull 8, 153 (1973)], respectively. The RVB Ansatz implemented in TurboRVB has a large variational freedom, including the Jastrow correlated Slater determinant as its simplest, but nontrivial case. Moreover, it has the remarkable advantage of remaining with an affordable computational cost, proportional to the one spent for the evaluation of a single Slater determinant. Therefore, its application to large systems is computationally feasible. The WF is expanded in a localized basis set. Several basis set functions are implemented, such as Gaussian, Slater, and mixed types, with no restriction on the choice of their contraction. The code implements the adjoint algorithmic differentiation that enables a very efficient evaluation of energy derivatives, comprising the ionic forces. Thus, one can perform structural optimizations and molecular dynamics in the canonical NVT ensemble at the VMC level. For the electronic part, a full WF optimization (Jastrow and antisymmetric parts together) is made possible, thanks to state-of-the-art stochastic algorithms for energy minimization. In the optimization procedure, the first guess can be obtained at the mean-field level by a built-in DFT driver. The code was efficiently parallelized by using a hybrid MPI-OpenMP protocol, which is also an ideal environment for exploiting the computational power of modern Graphics Processing Unit accelerators.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2000

HST Studies of the Chromospheres, Wind, and Mass-Loss Rates of Cool Giant and Supergiant Stars

UV spectra of K-M giant and supergiant stars and of carbon stars have been acquired with the Goddard High Resolution Spectrograph (GHRS) on the Hubble Space Telescope (HST). These spectra have been used to measure chromospheric flow and turbulent velocities, study the acceleration of their stellar winds, acquire constraints on their outer atmospheric structure, and enable estimates of their mass-loss rates. Results from our observations of the giant stars Gamma Dra (K5 III hybrid), Alpha Tau (K5 III), Gamma Cru (M3.4 III), Mu Gem (M3 IIIab), and 30 Her (MG III), the supergiants Alpha Ori (M2 Iab) and Lambda Vel (K5 Ib), and the carbon stars TX Psc (NO; C6,2) and TW Hor (NO; C7,2) will be summarized and compared. The high resolution and wavelength accuracy of these data have allowed the direct measurement of the acceleration of the stellar winds in the chromospheres of several of these stars (from initial velocities of 3-9 km/s to upper velocities of 15-25 km/s) and of the chromospheric macroturbulence (-25-35 km/s). The high signal-to-noise and large dynamic range of these spectra have allowed the detection and identification of numerous new emission features, including weak C IV emission indicative of hot transition-region plasma in the non-coronal giant Alpha Tau, many new fluorescent lines of Fe II, and the first detection of fluorescent molecular hydrogen emission and of Ca II recombination lines in the UV spectrum of a giant star. The UV spectrum of two carbon stars have been studied with unprecedented resolution and reveal extraordinarily complicated Mg II lines nearly smothered by circumstellar absorptions. Finally, comparison of synthetic UV emission line profiles computed with the Lamers et al. (1987) Sobolev with Exact Integration (SEI) code with observations of chromospheric emission lines overlain with wind absorption features provides estimates of the mass-loss rates for four of these stars.

Carpenter, Kenneth G.↗
DOE OSTI2023

Spin–vibronic coherence drives singlet–triplet conversion

Design-specific control over the transitions between excited electronic states with different spin multiplicities is of the utmost importance in molecular and materials chemistry. Previous studies have indicated that the combination of spin-orbit and vibronic effects, collectively termed the spin-vibronic effect, can accelerate quantum-mechanically forbidden transitions at non-adiabatic crossings. However, it has been difficult to identify precise experimental manifestations of the spin-vibronic mechanism. Here we present coherence spectroscopy experiments that reveal the interplay between the spin, electronic and vibrational degrees of freedom that drive efficient singlet-triplet conversion in four structurally related dinuclear Pt(II) metal-metal-to-ligand charge-transfer (MMLCT) complexes. Photoexcitation activates the formation of a Pt-Pt bond, launching a stretching vibrational wavepacket. The molecular-structure-dependent decoherence and recoherence dynamics of this wavepacket resolve the spin-vibronic mechanism. Further, we find that vectorial motion along the Pt-Pt stretching coordinates tunes the singlet and intermediate-state energy gap irreversibly towards the conical intersection and subsequently drives formation of the lowest stable triplet state in a ratcheting fashion. This work demonstrates the viability of using vibronic coherences as probes to clarify the interplay among spin, electronic and nuclear dynamics in spin-conversion processes, and this could inspire new modular designs to tailor the properties of excited states.

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

Efficient Treatment of Large Active Spaces through Multi-GPU Parallel Implementation of Direct Configuration Interaction

In this study, we have extended our graphical processing unit (GPU)-accelerated direct configuration interaction program to multiple devices, reducing iteration times for configuration spaces of 165 million determinants to only 3 s using NVIDIA P100 GPUs. Similar improvements in the one- and two-particle reduced density matrix formation allow for fast analytical energy gradients and electronic properties. Our parallel algorithm enables the calculation of arbitrarily large configuration spaces (limited only by available system memory), with iteration times of 13 min for an active space of 18 electrons in 18 orbitals (2.4 billion determinants) using six consumer grade NVIDIA 1080Ti GPUs. These advances enable routine molecular dynamics simulations, geometry optimizations, and absorption spectrum calculations for molecules with large configuration spaces, a task that has heretofore required massive computational effort. In this work, we demonstrate the utility of our program by generating the absorption spectrum for diphenyl acetylene at the floating occupation molecular orbital complete active space configuration interaction level of theory. Lastly, several active spaces were investigated to assess the dependence of spectral features on orbital space dimension.

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

Cabana: A Performance Portable Library for Particle-Based Simulations

Particle-based simulations are ubiquitous throughout many fields of computational science and engineering, spanning the atomistic level with molecular dynamics (MD), to mesoscale particle-in-cell (PIC) simulations for solid mechanics, device-scale modeling with PIC methods for plasma physics, and massive N-body cosmology simulations of galaxy structures, with many other methods in between (Hockney & Eastwood, 1989). While these methods use particles to represent significantly different entities with completely different physical models, many low-level details are shared including performant algorithms for short- and/or long-range particle interactions, multi-node particle communication patterns, and other data management tasks such as particle sorting and neighbor list construction. Cabana is a performance portable library for particle-based simulations, developed as part of the Co-Design Center for Particle Applications (CoPA) within the Exascale Computing Project (ECP) (Alexander et al., 2020). The CoPA project and its full development scope, including ECP partner applications, algorithm development, and similar software libraries for quantum MD, is described in (Mniszewski et al., 2021). Cabana uses the Kokkos library for on-node parallelism (Edwards et al., 2014; Trott et al., 2022), enabling simulation on multi-core CPU and GPU architectures, and MPI for GPU-aware, multi-node communication. Cabana provides particle simulation capabilities on almost all current Kokkos backends, including serial execution, OpenMP (including OpenMP-Target for GPUs), CUDA (NVIDIA GPUs), HIP (AMD GPUs), and SYCL (Intel GPUs), providing a clear path for the coming generation of accelerator-based exascale hardware. Cabana builds on Kokkos by providing new particle data structures and particle algorithms resulting in a similar execution policy-based, node-level programming model that is intended to be used in addition to the core Kokkos library within an application. Cabana is designed as an application and physics agnostic, but particle-specific toolkit which can either be used to generate a new application, or to be used as needed in existing applications at various levels of invasiveness including through interfaces that wrap user memory in existing data structures.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Machine learning dielectric screening for the simulation of excited state properties of molecules and materials

Accurate and efficient calculations of absorption spectra of molecules and materials are essential for the understanding and rational design of broad classes of systems. Solving the Bethe–Salpeter equation (BSE) for electron–hole pairs usually yields accurate predictions of absorption spectra, but it is computationally expensive, especially if thermal averages of spectra computed for multiple configurations are required. We present a method based on machine learning to evaluate a key quantity entering the definition of absorption spectra: the dielectric screening. We show that our approach yields a model for the screening that is transferable between multiple configurations sampled during first principles molecular dynamics simulations; hence it leads to a substantial improvement in the efficiency of calculations of finite temperature spectra. We obtained computational gains of one to two orders of magnitude for systems with 50 to 500 atoms, including liquids, solids, nanostructures, and solid/liquid interfaces. Importantly, the models of dielectric screening derived here may be used not only in the solution of the BSE but also in developing functionals for time-dependent density functional theory (TDDFT) calculations of homogeneous and heterogeneous systems. Overall, our work provides a strategy to combine machine learning with electronic structure calculations to accelerate first principles simulations of excited-state properties.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Active Learning Framework

Machine learning (ML) of interatomic potentials show great promise to accelerate scientific simulation, e.g., by emulating expensive computations at a high accuracy but much reduced computational cost. Training datasets are calculated from computationally expensive ab initio quantum mechanics methods, density functional theory (DFT). Trained on this data, an ML model can be very successful in predicting energy and forces for new atomic configurations. A critical factor is the quality and diversity of the training dataset. Thus, a highly automated approach to dataset construction based on active learning framework is designed suitable for material physics. The active learning scheme begins with fully randomized atomic configurations. Then, many Molecular Dynamics (MD) trajectories are simulated using current ML potentials, where each MD trajectory is initialized to a random disordered configuration. The temperature is varied in order to diversify the sampled configuration during these simulations. The variance of predictions for eight neural networks within an ensemble is analyzed to determine whether the model is operating as expected. This helps in determining whether collecting more data would be helpful to the model by checking the ensemble variance is greater than the threshold. In this case, the MD trajectory is terminated and the final atomic configuration is placed on a queue (SQL database) for DFT calculations and added to training dataset. Periodically, ML model is retrained to the updated training model. This Active Learning loop is iterated until the cost of MD simulations becomes prohibitively expensive. The MD simulations will hopefully be sufficiently robust to support nucleation after many active learning iterations. In this sense, active learning scheme must automatically discover the important low energy and nonequilibrium physics.

Nebgen, Benjamin↗
DOE OSTI2022

Dynamic, Adaptive, Systems and Materials: Complex, Simple and Emergent Behaviors

This program has been funded by DoE/BES for twenty years. It has moved into and out of various subjects as it has developed, but it has retained its focus on complexity and complex systems. The project has evolved in the following way: Self-Assembly and Biomimetic Self-Assembly: All self-assembling systems depend upon a minimum of two types of interaction: a repulsion and an attraction. For the familiar molecular systems, attractive interactions are typically hydrogen bonds and electrostatic interactions. Repulsive interactions include steric effects, hydrophobic effects (in biological systems), and charge-charge repulsion. We have expanded this repertoire to include surface interactions, magnetic interactions, and others. I list these systems in the order in which we have explored them: i) A key emphasis in current work is in understanding how the movement of ions in a magnetic field (the Lorentz effect) interacts with catalytic systems. We have demonstrated that an acceleration in rate of reduction of CO 2 to CO can be accomplished by applying an external magnetic field. This acceleration is largely due to the application of the Lorentz effect on mass transport at the catalyst’s surface. ii) We have also extensively explored the influence of electrostatics, as exhibited in self-assembling systems, by tribocharging. iii) Another key system involves surface tension effects; examples include interactions between heavy particles floating at a liquid-air interface, and interacting by changes in surface area; interactions of bubbles and bubble rafts, behaviors of bubble trains in microfluidic networks, and behaviors of microorganisms in constraining environments. iv) This work has intentionally de-emphasized biological systems; but it does include some work on protein-ligand interactions and interactions among microorganisms. v) We have also explored applications of some of these effects, these explorations include bubble rafts as diffraction gratings, exploration of the structures that can be obtained by tribocharging and uses of these structures in exploring nucleation and melting of crystals. vi) Although not a major focus of this work, several other topics have emerged and offer opportunities for future work. These include the behavior of bubble trains and bubble rafts in microfluidic systems. A particularly interesting example is the formation of bubble trains that repeat in the alteration of large and small bubbles according to rules we do not presently understand, but are uniquely large-period oscillating systems. These systems offer a new route into understanding the instabilities of the type represented by oscillations. vii) We have also begun exploratory projects on magnetic levitation (especially to determine molecular density), and information storage (in molecules). Magnetic Levitation: Self-assembling and biomimetic systems require both attraction and repulsion. We have used electrostatics (tribocharging), interfacial free-energies (surface tension and related forces) and others. Potential uses include reconfigurable diffraction gratings and liquid lenses; exploration of mechanisms in tribocharging; tunneling in EGaIn junctions; and bubble trains (especially in micro-fluidic systems). Examples of systems representing these topics is included in the following papers: Complexity: Disks rotating at a water-air interface; Benard-Marangoni effects; Vortex-Crystals from spinning magnetic disks (Marangoni effects); EGaIn Electrode to study quantum tunneling; Self-Assembly of electrostatically-charged metallic spheres (electrets); Dynamically reconfigurable lens; Using computational designs of ligands for enzymes; Electrostatic self-assembly by tribocharging; Monodisperse bubble trains in microchannel systems; Inverted dripping faucet; Flames; Printing of micro-organisms to regenerate the “ink” of printing device; Using micro-organisms to move loads (“microoxen”); Motion of bacterial swarms near surfaces; Making monodisperse particles in microfluidic systems; Coding/decoding of information stored in droplet trains in microfluidic networks; Magnetic levitation; and Information storage. i) Tribocharging. The change in focus of this work on electrets from the fundamentals of charging to applications of these materials in studying self-assembly using electrostatic interactions. ii) Bubbles in Microchannels. The realization that systems of bubbles in microchannels represented a major opportunity to study complexity in a very tractable system, and the development of a semi-quantitative theory of this subject. iii) Flames. The growth of “flames” remains an exploratory subject for the research, although their currently relatively little active work involving it ongoing. iv) Systems with Microorganisms. The removal of work in biological systems from this project. Based on work supported in this program, we now have a significant project on the development of microfluidic tools for studying C. elegans (a nematode), but this work was not appropriate for a program focused on complexity, and we developed separate support for it. (It is, however, an example of successful seeding of a new area by BES.) The work on electrets has gone through a period in which a part of the program was the subject of a MURI; the focus of this work was to develop materials that did not charge electrostatically on friction or contact. The MURI is now over, and the work on dynamic self-assembly (supported by BES) is the major focus. “Flames” has also enjoyed synergistic support, in terms of a project supported by DARPA on flame suppression (in the absence of extinguishing agents, using acoustic and electrostatic interactions). This work was helpful in understanding some of the basics of flames, but is entirely distinct from the BES focus in complexity. A growing interest is in the Lorentz effect. The Lorentz effect is the force exerted on charged particles (electrons, ions, charged molecules) when they move through a perpendicular magnetic field. The Lorentz effect is almost ubiquitous in modern technology: examples of applications include electric motors, dynamos, cathode ray tubes, many batteries, and most systems that control electrical currents with magnetic forces. We have begun to explore the Lorentz effect in electrochemical systems and heterogeneous catalytic systems involving charged organic species and inorganic ions. This work is still at an early stage, but initial studies that Lorentz effects can be large when ions move through magnetic fields, or magnetic fields move in the presence of ions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Thermo-mechanical deterioration and molecular degradation of 3D-printed methacrylate-based polymer in various chemical environments

The advancement of additive manufacturing (AM) has accelerated the development of stereolithographic (SLA) photo-curable resins, particularly methacrylate-based polymers, due to the ability to their high-resolution, robust mechanical properties, and suitability. However, their long-term performance in chemical environments remains poorly understood. This study investigates the extent and mechanisms of degradation on an SLA-printed methacrylate-based polymer subjected to various chemicals, including polar and non-polar solvents, as well as strongly acidic aqueous solutions over a 12-week accelerated aging period. A comprehensive analytical approach incorporating swelling kinetics, surface morphology, tensile and dynamic mechanical analysis (DMA), Fourier-transform infrared (FTIR) spectroscopy, and mass spectrometry (MS) was employed to characterize chemical absorption, structural integrity, and leached products. Results reveal that degradation severity is governed by both the polarity and reactivity of the chemical environment. Notably, exposure to 6 mol L -1 HNO₃ induced the most severe deterioration, with over threefold higher swelling compared to other media, and significant reductions in tensile strength, tensile modulus, and glass transition temperature (T g ). In contrast, specimens aged in non-polar solvents (xylene and dodecane) exhibited negligible chemical interaction and retained mechanical performance. FTIR and MS analyses identified acid-catalyzed hydrolysis of ester groups as prominent degradation pathways in acidic media, while diffusion-controlled plasticization prevailed in polar solvents. Furthermore, this study provides valuable insights into the chemical stability of SLA-printed polymers and develops predictive degradation profiles that are crucial for designing durable polymer systems for advanced industrial use.

36 MATERIALS SCIENCE↗
NASA NTRS2000

Applications in Bioastronautics and Bioinformatics: Early Radiation Cataracts Detected by Noninvasive, Quantitative, and Remote Means

Human exploration of Mars is a key goal in NASA's exploration planning in the next 20 years. Maintaining crew health and good vision is certainly an important aspect of achieving a successful mission. Continuous radiation exposure is a risk factor for radiation-induced cataracts in astronauts because radiation exposure in space travel has the potential of accelerating the aging process (ref. 1). A patented compact device (ref. 2) based on the technique of dynamic light scattering (DLS) was designed for monitoring an astronaut's ocular health during long-duration space travel. This capability of early diagnosis, unmatched by any other clinical technique in use today, may enable prompt initiation of preventive/curative therapy. An Internet web-based system integrating photon correlation data and controlling the hardware to monitor cataract development in vivo at a remote site in real time (teleophthalmology) is currently being developed. The new technology detects cataracts very early (at the molecular level). Cataract studies onboard the International Space Station will be helpful in quantifying any adverse effect of radiation to ocular health. The normal lens in a human eye, situated behind the cornea, is a transparent tissue. It contains 35 wt % protein and 65 wt % water. Aging, disease (e.g., diabetes), smoking, dehydration, malnutrition, and exposure to ultraviolet light and ionizing radiation can cause agglomeration of the lens proteins. Protein aggregation can take place anywhere in the lens, causing lens opacity. The aggregation and opacification could produce nuclear (central portion of the lens) or cortical (peripheral) cataracts. Nuclear and posterior subcapsular (the membrane's capsule surrounds the whole lens) cataracts, being on the visual optical axis of the eye, cause visual impairment that can finally lead to blindness. The lens proteins, in their native state, are small in size. As a cataract develops, this size grows from a few nanometers (transparent) to several micrometers (cloudy). Ansari and Datiles have shown that DLS can detect cataracts at least two to three orders of magnitude earlier noninvasively and quantitatively than the best imaging (Scheimpflug) techniques in clinical use today (ref. 3).

Ansari, Rafat R.↗
DOE OSTI2024

Direct observation of ultrafast cluster dynamics in supercritical carbon dioxide using X-ray Photon Correlation Spectroscopy

Supercritical fluids exhibit distinct thermodynamic and transport properties, making them of particular interest for a wide range of scientific and engineering applications. These anomalous properties emerge from structural heterogeneities due to the formation of molecular clusters at conditions above the critical point. While the static behavior of these clusters and their effects on the thermodynamic response functions have been recognized, the relation between the ultrafast cluster dynamics and transport properties remains elusive. By measuring the intermediate scattering function in carbon dioxide at conditions near the critical point with X-ray photon correlation spectroscopy, we directly capture the cross-over dynamics between 4 and 13 picoseconds, revealing the transition between ballistic and diffusive motion. Complementary analysis using large-scale molecular dynamics simulations reveals that this behavior arises from collisions between unbound molecules and clusters. This study provides direct evidence of the ultrafast momentum exchange between clusters, which has significant impact on transport properties, solvation processes, and reaction kinetics in supercritical fluids.

carbon capture and storage↗