Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Classical 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 631 records · Page 35

SULI Intern Final Report: Computationally Investigating Hydrogen Thermo-Diffusion in Yttrium Hydride Using Multiscale Methods

The renaissance of nuclear energy has arrived, heralding an age of abundant inexpensive clean energy, and renewed space exploration. In nuclear-powered spacecraft and microreactors, safety and size are of utmost importance. Yttrium Hydride (YHx) is being researched for its utility as a neutron moderator in nuclear reactors; the hydrogen in YHx slows down neutrons, enabling a continuous nuclear reaction in the reactor. This has the benefit of allowing reactors to be more safe, compact, and efficient. The goal of this effort is to computationally predict the coefficient of temperature-dependent hydrogen diffusion within YHx, the Soret coefficient. This parameter is essential for determining the safe operating modes of YHx moderators. Zirconium Hydride (ZrHx) is used in the Training, Research, Isotopes, General Atomics (TRIGA) reactor, is the reference material for these calculations. In this work, nanoscale atomic modeling in the Vienna Ab initio Simulation Package (VASP) is combined with the mesoscale finite element phase-field module in the Multiphysics Object-Oriented Simulation Environment (MOOSE); this culminates in a new multiscale computational method to simulate Soret diffusion of hydrogen in YHx. This data is useful for predicting experimental outcomes. This workflow involves convergence testing followed by static, Nudged Elastic Band (NEB), Quasi-Harmonic Approximation (QHA), and Molecular Dynamics (MD) calculations - linked with phase field simulation. NEB simulates hydrogen migration, while QHA and MD predict temperature-dependent properties. The static calculations align with literature, and preliminary NEB and QHA calculations yield accurate results. Once the atomic calculations are complete, we will incorporate Electron Backscatter Diffraction (EBSD) images and VASP-generated parameters into the phase field module to simulate intra- and intergranular transport of hydrogen in ZrHx and YHx. Future research will extend our approach to fuel-moderator materials systems such as Uranium-Yttrium Hydride (U-YHx). This work contributes to the development of advanced nuclear energy solutions for space travel.

36 - MATERIALS SCIENCE↗

Synergistic Enhancement of Lead and Selenate Uptake at the Barite (001)–Water Interface

The interactions of heavy metals with minerals influence the mobility and bioavailability of toxic elements in natural aqueous environments. The sorption of heavy metals on covalently bonded minerals is generally well described by surface complexation models (SCMs). However, understanding sorption on sparingly soluble minerals is challenging because of the dynamically evolving chemistry of sorbent surfaces. The interpretation can be even more complicated when multiple metal ions compete for sorption. In the present study, we observed synergistically enhanced uptake of lead and selenate on the barite (001) surface through two sorption mechanisms: lattice incorporation that dominates at lower coverages and two-dimensional monolayer growth that dominates at higher coverages. Further, we also observed a systematic increase in the sorption affinity with increasing co-sorbed ion coverages, different from the assumption of invariant binding constants for individual adsorption processes in classical SCMs. Computational simulations showed thermodynamically favorable co-incorporation of lead and selenate by simultaneously substituting for barium and sulfate in neighboring sites, resulting in the formation of molecular clusters that locally match the net dimension of the substrate lattice. These results emphasize the importance of ion–ion interactions at mineral–water interfaces that control the fate and transport of contaminants in the environment.

54 ENVIRONMENTAL SCIENCES↗

Geometric transformation and three-dimensional hopping of Hopf solitons

Arising in many branches of physics, Hopf solitons are three-dimensional particle-like field distortions with nontrivial topology described by the Hopf map. Despite their recent discovery in colloids and liquid crystals, the requirement of applied fields or confinement for stability impedes their utility in technological applications. Here we demonstrate stable Hopf solitons in a liquid crystal material without these requirements as a result of enhanced stability by tuning anisotropy of parameters that describe energetic costs of different gradient components in the molecular alignment field. Nevertheless, electric fields allow for inter-transformation of Hopf solitons between different geometric embodiments, as well as for their three-dimensional hopping-like dynamics in response to electric pulses. Numerical modelling reproduces both the equilibrium structure and topology-preserving out-of-equilibrium evolution of the soliton during switching and motions. Our findings may enable myriads of solitonic condensed matter phases and active matter systems, as well as their technological applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Optical-Gating of Spin-based Quantum States (Final Report)

This project developed a generalizable strategy for optical control of spin-based quantum states through the Photoisomerization-Induced Spin–Charge Excited State (PISCES) process, in which a photochromic spirooxazine ligand bound to an electronically bistable metal complex enables reversible, room-temperature optical switching of spin states in both dilute matrices and the solid state. A multidisciplinary team integrated synthesis, spectroscopy, theory, and surface science to establish design principles for opto-spintronic quantum materials. Systematic structure–property studies across four spirooxazine series revealed that donor–π–acceptor interactions govern thermal isomerization, which proceeds via a rotational pathway. Pulsed EPR measurements demonstrated millisecond-scale coherence times in cobalt-dioxolene-spirooxazine complexes, establishing their promise as optically addressable molecular qubits; the approach was extended to Cu(II), Ni(II), and lanthanide complexes, including optical modulation of zero-field splittings and clock transitions. Two-dimensional IR spectroscopy resolved picosecond ground-state population exchange dynamics, while ab initio calculations elucidated photoisomerization mechanisms and parametrized Hamiltonians for spin relaxation and decoherence. Toward device integration, the team fabricated 2D polymer scaffolds for covalent qubit organization, developed tip-enhanced Raman imaging of molecular monolayers, and demonstrated optical gating of gate voltage in spirooxazine-graphene field-effect transistors. The project produced numerous publications, dissertations, and presentations, and trained over 30 researchers.

2D polymers↗

Statistical mechanics in biology: how ubiquitous are long-range correlations?

The purpose of this opening talk is to describe examples of recent progress in applying statistical mechanics to biological systems. We first briefly review several biological systems, and then focus on the fractal features characterized by the long-range correlations found recently in DNA sequences containing non-coding material. We discuss the evidence supporting the finding that for sequences containing only coding regions, there are no long-range correlations. We also discuss the recent finding that the exponent alpha characterizing the long-range correlations increases with evolution, and we discuss two related models, the insertion model and the insertion-deletion model, that may account for the presence of long-range correlations. Finally, we summarize the analysis of long-term data on human heartbeats (up to 10(4) heart beats) that supports the possibility that the successive increments in the cardiac beat-to-beat intervals of healthy subjects display scale-invariant, long-range "anti-correlations" (a tendency to beat faster is balanced by a tendency to beat slower later on). In contrast, for a group of subjects with severe heart disease, long-range correlations vanish. This finding suggests that the classical theory of homeostasis, according to which stable physiological processes seek to maintain "constancy," should be extended to account for this type of dynamical, far from equilibrium, behavior.

Non-NASA Center↗

Mesoscale Modeling of Hydrogels Under Frictional Shear Stress

Hydrogels are three-dimensional networks of hydrophilic polymers often used as a simplified model of hydrated biological materials, from cartilaginous joints to the ocular tear film. However, the lubrication mechanisms of hydrogels remain poorly understood, partly due to their complex polymeric structure, which creates blurred interfaces during sliding that are challenging to study experimentally. In this study, we employ dissipative particle dynamics (DPD) to investigate the frictional behavior of a polymeric hydrogel network sliding against a solid wall in an explicit viscous solvent. This computational approach enables us to model hydrodynamic interactions and mesoscale polymer dynamics, capturing key aspects of hydrogel friction. Our simulations reveal that hydrogel friction is governed by the interplay between polymer relaxation and viscous shear, characterized by the Weissenberg number (Wi). At low Wi, friction coefficient remain nearly constant, dominated by polymer relaxation. However, at higher Wi, friction is dominated by viscous drag within a near-wall solvent layer, leading to a linear increase in friction coefficient with Wi. Furthermore, our results demonstrate an inverse relationship between the friction coefficient and the applied normal load, consistent with experimental observations. This work provides new insights into the fundamental tribological properties of hydrogels, shedding light on the micromechanics of hydrogel friction. Improving our understanding of hydrogel structure and dynamics under friction advances our knowledge of the mechanisms regulating biological lubrication in health and disease.

36 MATERIALS SCIENCE↗

Hydrodynamically Controlled Self‐Organization in Mixtures of Active and Passive Colloids

Abstract Active particles are known to exhibit collective behavior and induce structure in a variety of soft‐matter systems. However, many naturally occurring complex fluids are mixtures of active and passive components. The authors examine how activity induces organization in such multi‐component systems. Mixtures of passive colloids and colloidal micromotors are investigated and it is observed that even a small fraction of active particles induces reorganization of the passive components in an intriguing series of phenomena. Experimental observations are combined with large‐scale simulations that explicitly resolve the near‐ and far‐field effects of the hydrodynamic flow and simultaneously accurately treat the fluid–colloid interfaces. It is demonstrated that neither conventional molecular dynamics simulations nor the reduction of hydrodynamic effects to phoretic attractions can explain the observed phenomena, which originate from the flow field that is generated by the active colloids and subsequently modified by the aggregating passive units. These findings not only offer insight into the organization of biological or synthetic active–passive mixtures, but also open avenues to controlling the behavior of passive building blocks by means of small amounts of active particles.

36 MATERIALS SCIENCE↗

Fast and accurate calculation of EXAFS Debye-Waller factors in U⁢O2 using the dynamical matrix method

Theoretical modeling of bonding dynamics in metal oxides is required for predicting their thermal conductivity, catalytic activity, and mechanical properties. A primary challenge is the scarcity of experimental methods for validating theoretical predictions of these atomic-scale dynamics. This work presents a workflow that uses experimental extended x-ray absorption fine structure (EXAFS) data collected at high temperatures to validate an interatomic force field for uranium dioxide (UO2), an important model material. The validated force field is then used to drive computationally intensive molecular dynamics (MD) simulations and as input for the much faster dynamical matrix Debye-Waller (DMDW) method. The predicted values of the Debye-Waller factors from the DMDW calculations are in good agreement with those obtained from the MD simulations, with residual pair-specific differences attributable to quantum zero-point motion at low temperatures and lattice anharmonicity at high temperatures. We further show that theoretical EXAFS spectra constructed directly from DMDW-derived Debye-Waller factors reproduce the experimental data (at relatively low temperatures) with accuracy comparable to full MD-EXAFS, providing an additional validation of the choice of the potential. This study establishes a validated, rapid computational pathway for modeling bond dynamics, naturally incorporating quantum nuclear\\\\r\\\\nstatistics absent in classical simulations, which are essential for the mechanistic understanding of complex oxide materials.

58 GEOSCIENCES↗

Quantum phases of Rydberg atoms on a kagome lattice

Significance Programmable quantum simulators based on Rydberg atom arrays have recently emerged as versatile platforms for exploring exotic many-body phases and quantum dynamics of strongly correlated systems. In this work, we theoretically investigate the quantum phases that can be realized by arranging such Rydberg atoms on a kagome lattice. Along with an extensive analysis of the states which break lattice symmetries due to classical correlations, we identify an intriguing regime that constitutes a promising candidate for hosting a phase with long-range quantum entanglement and topological order. Our results provide a route to experimentally realizing and probing highly entangled quantum matter.

74 ATOMIC AND MOLECULAR PHYSICS↗

Berni Julian Alder, theoretical physicist and inventor of molecular dynamics, 1925–2020

Berni Julian Alder, one of the leading figures in the invention of molecular dynamics simulations used for a wide array of problems in physics and chemistry, died on September 7th, 2020. His career, spanning more than 65 years, transformed statistical mechanics, many body physics, the study of chemistry and the microscopic dynamics of fluids, by making atomistic computational simulation (in parallel with traditional theory and experiment) a new pathway to unexpected discoveries. Among his many honors, the CECAM prize, recognizing exceptional contributions to the simulation of the microscopic properties of matter is named for him. He was awarded the National Medal of Science by President Obama in 2008.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A window on the composition of the early solar nebula: 2014MU69, Pluto, and Phoebe

The initial chemical composition of any solar nebula will depend upon the degree to which 1) organic and ice components form on dust grains, 2) organic and molecular species form in the gas phase, 3) organics and ices are exchanged between the gas and solid state, and 4) the precursor and newly formed (more complex) materials survive and are modified in the developing planetary system. Infrared and radio observations of star-forming regions reveal that complex chemistry occurs on icy grains, sometimes before stars even form. Additional processing, through the protosolar disk and within the solar nebula further modifies most, but probably not all, of the initial materials. In fact, the modern Solar System still carries a fraction of its interstellar inheritance [Alexander et al., 2017]. Here we focus on three examples of small bodies in our Solar System, each containing chemical and dynamical clues to its origin and evolution: the small, cold-classical Kuiper Belt object (KBO) 2014MU69 , Pluto and Saturn’s moon, Phoebe. The New Horizons flyby of 2014 MU69 has given the first view of an unaltered body composed of material originally in the solar nebula at ~45 AU. The spectrum reveals methanol ice (not commonly found), a possible detection of water ice (Stern et al. 2019). Pluto’s internal and surface inventory of volatiles and complex organics, together with active geological processes including cryo-volcanism, indicate a surprising level of activity on a body in the outermost region of the Solar System, and the fluid that emerges from subsurface reservoirs may contain material inherited from the solar nebula (Cruikshank et al. 2019). Meanwhile, Saturn's captured moon, Phoebe, carries high D/H in H2O [Clark et al. 2018)], and complex organics (Cruikshank et al. 2008), both consistent with its formation in, and inheritance from, the outer solar nebula. Together, these objects provide windows on the origin and evolution of our Solar System and constraints to be considered in future chemical and physical models of PPDs. The Spitzer Space telescope discovered the ring of dust around Saturn that emanated from the bombardment of Phoebe in the recent past, dislodging primitive Solar System material from deep within that satellite. In addition, Spitzer has provided detailed studies of the chemistry of star forming regions where methanol and other ices are readily observed. As a result of these and other observations, protoplanetary disk models are beginning to combine both dynamical and chemical complexities into their codes, resulting in new insights for exoplanetary systems.

Yvonne J. Pendleton↗

High-pressure phase transitions in a laser directed energy deposited Fe-33Cu Alloy

Here, additively manufactured Fe-Cu alloys contain both equilibrium face-centered cubic (FCC) and metastable body-centered cubic (BCC) crystal structure Cu precipitates depending on their size. However, the stability of these nanoscale precipitates under extreme conditions such as high pressures has not been reported. This study investigates the phase transformations and microstructural stability of laser directed energy deposition (DED-LB) made Fe67Cu33 alloy (nominal composition in at.%) under high static pressure deformation using in-situ synchrotron X-ray diffraction under pressure, postmortem high-resolution scanning transmission electron microscopy (HR-STEM), and molecular dynamics (MD) simulations. In-situ XRD results reveal a reversible phase transformation from BCC to the hexagonal close-packed (HCP) structure in the Fe grains at an onset pressure of 16.4 GPa, significantly higher than reported in the literature for pure Fe. Although no high-pressure phase transition was observed in the FCC Cu grains through XRD, HR-STEM analysis uncovers a phase transition to the HCP structure in nanoscale metastable BCC Cu precipitates within the BCC Fe matrix. After decompression, the Fe matrix reverted back to the BCC structure with periodic lath martensite, while regions of the nanoscale BCC Cu precipitates retained the metastable HCP structure. MD simulations support the BCC → HCP transition in the embedded nanoscale coherent Cu precipitates, consistent with the classical Burgers mechanism. Thus, by leveraging in-situ XRD observation, postmortem high-resolution S/TEM microscopy, and MD simulations, this study offers profound insights into the distinctive phase transformations induced by high pressure in DED-LB Fe67Cu33 alloy, which is distinguished by its hierarchical microstructure.

36 MATERIALS SCIENCE↗

Tensor Network Path Integral Study of Dynamics in B850 LH2 Ring with Atomistically Derived Vibrations

The recently introduced multisite tensor network path integral (MS-TNPI) allows simulation of extended quantum systems coupled to dissipative media. We use MS-TNPI to simulate the exciton transport and the absorption spectrum of a B850 bacteriochlorophyll (BChl) ring. The MS-TNPI network is extended to account for the ring topology of the B850 system. Accurate molecular-dynamics-based description of the molecular vibrations and the protein scaffold is incorporated through the framework of Feynman–Vernon influence functional. To relate the present work with the excitonic picture, an exploration of the absorption spectrum is done by simulating it using approximate and topologically consistent transition dipole moment vectors. Comparison of these numerically exact MS-TNPI absorption spectra are shown with second-order cumulant approximations. Finally, the effect of temperature on both the exact and the approximate spectra is also explored.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Minimization of Disorder as a Key Design Principle for Natural Sizes of Light Harvesting 2 Complexes

The light harvesting 2 (LH2) complex of purple bacteria has excellent energy conversion efficiency. Clarifying the design principle behind such efficiency at the atomistic level is crucial for understanding its structure–function relationship and can be utilized for the design of artificial light harvesting systems. To this end, we conducted comprehensive computational investigation of the dynamical and statistical nature of electronic excited states of pigment molecules in a natural LH2 complex with 9-fold symmetry and its two non-natural in silico analogues with 6- and 12-fold symmetries. To ensure reliable and efficient all-atomistic molecular dynamics simulations, we combined a well established interpolation approach for the construction of the potential energy surface with a neural network machine learning approach. Outcomes of these calculations clarify that non-natural forms of LH2-type complexes have significantly larger quasistatic disorder than those for the natural one. In addition, non-natural systems have more disruptions of the hydrogen bonding, underscoring its crucial role for reducing the disorder. On the other hand, local environmental dynamics are relatively insensitive to the structural changes although there is moderate enhancement in the anharmonic or interatomic components for the synthetic ones. These findings based on all-atomistic simulations provide direct computational evidence that the structure and sizes of natural LH2 complexes are designed to minimize the energetic disorder. We analyze quantitative implications of these for the energy transferring capability of the LH2 complex.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

New principles of self‐organization created through the interplay of DNA condensates, microtubules, and motors

Bioinspired design—which holds great promise for a new generation of materials that are robust to defects, scalable under green manufacture, environmentally responsive, and programmably reconfigurable—requires mastery over molecular self-organization. Yet, from its specific mechanisms to most general architectures, the principles governing self-organization remain poorly understood and not even fully enumerated. For living systems, one obvious architectural principle is the modular reuse of a few simple molecular components in myriad combinations to achieve more complex phenomena. For example, the mechanical tasks of a cell are driven by the nonequilibrium dynamics of cytoskeletal filaments and molecular motors—the same filaments and motors, reprogrammed by a variety of modulators, perform tasks ranging from cell movement to division. Similarly, many compartmentalization tasks are performed by liquid-like condensates of simple components, which act as membraneless organelles to localize particular molecules in space and time (e.g. for gene regulation or RNA processing). In a few cases, condensates combine and interact with the cytoskeleton to create still more complex phenomena, e.g. the nucleation of microtubule asters from the centrosome (a protein condensate) to form the mitotic spindle during cell division. Very little is known about the fundamental mechanisms of such filament-plus-condensate phenomena. Despite few examples, the landscape of behaviors that can be achieved through the combination of condensates, filaments, and motors appears vast. However, exploration has been hindered by a lack of systems that have sufficiently programmable and dynamically tunable interactions between component condensates, filaments, and motors. We proposed to combine programmable DNA condensates, filamentous microtubules, and light-controlled motors into self-organizing systems whose principles go beyond those that have been observed in nature. In one limit, our systems will use microtubules and motors to create the molecular analog of a network of roads, which will organize droplets of DNA condensates capable of carrying molecular cargo. DNA condensates coupled to motors will flow from one microtubule aster hub to another, with their direction and timing controlled by DNA circuits. In another limit, microtubules will swim through bulk DNA condensates and exhibit strong interactions with boundaries between different types of condensates. Microtubule swimmers will reflect, get trapped, or refract at boundaries, under a mechanical analog of the classical optical index of refraction. DNA condensates having different mechanical indexes of refraction will be used to construct the analog of optical lenses, so that microtubule swimmers can be manipulated like light—collimated, diffracted, focused, and sorted based on properties analogous to wavelength. These two limits define two new architectures, within which multiple new mechanistic principles for self-organization will be discovered and explored. To explore these architectures, the motor-based coupling between DNA condensates and filaments will be controlled in time and space through the use of opto-proteins that create reversible links between DNA condensates and motors upon illumination. For each principle of interest, patterns of light will create virtual experiments by defining patterns of activity where DNA condensates walk along filaments, or filaments swim through condensates, and patterns of inactivity which will serve either as controls, or as boundary conditions vital to create the desired phenomena. This research serves the goals of Basic Energy Sciences Biomolecular Material Program by elucidating the principles by which the emergent, nonequilibrium behavior of collections of DNA condensates, motors, and microtubules can be programmed by environmental light patterns to create complex motion and materials transport. Because DNA condensates can be readily coupled to virtually any high performance nanomaterial, from carbon nanotubes, to metal nanoparticles, to light harvesting systems, this work provides a path to the construction, self-maintenance and reconfiguration of materials relevant to the Department of Energy.

60 APPLIED LIFE SCIENCES↗

Beyond Thermodynamics: Assessing the Dynamical Softness of Hydrated Ions from First Principles

Ion (de)hydration is a key rate-determining step in interfacial processes from corrosion to electrochemical energy storage. However, predicting the kinetics of ion (de)hydration remains challenging, prompting the use of static proxies such as hydration energy and valence. While useful for assessing thermodynamic preferences, such descriptors cannot fully capture the dynamical softness of the hydration shell that dictates kinetics. Accordingly, we use first-principles molecular dynamics to analyze hydration shell softness for a diverse set of metal cations. Additionally, three dynamic metrics are introduced to intuitively describe the bond rigidity, shape deformability, and exchange fluidity of the solvation shell. Together, these metrics capture the relevant physics in the static descriptors, while offering a far more complete and efficient representation for the overall propensity for (de)hydration. Application to the hydrated ion set demonstrates a weak connection between dynamical softness and hydration energy, confirming that dynamical descriptors of hydration are key for correctly describing ion transfer processes.

36 MATERIALS SCIENCE↗

Temperature and radiation effects on brittle versus ductile fracture behavior in miscible phase boundaries: insight from atomistic simulations

Temperature- and irradiation-assisted failure mechanisms in miscible phase boundaries are clarified via atomistic calculations. We first establish the temperature-dependent brittle-to-ductile transition in U–Zr miscible phase boundaries. Our results confirm that these boundaries are mostly brittle at low temperatures and ductile at elevated temperatures. We then investigate the changes induced by irradiation on the fracture mechanisms in such phase boundaries. The irradiation-induced defect accumulation follows a multi-stage process that starts with the accumulation of isolated small dislocation loops before transitioning to the saturation and growth of larger dislocation loops and end up with a reorganization into forest dislocations. The accumulation of loops is the primary feature to participate in the delineation between brittle and ductile interfacial fracture in irradiated phase boundaries. At low damage levels, radiation defect interactions with the crack tip are limited and U–Zr miscible boundaries fail through the emission of dislocations ahead of the crack tip followed by brittle cleavage in agreement with the classical Griffith’s criterion for crack stability. At higher damage levels, the failure mode transitions from brittle crack growth to ductile void growth. In this case, the microcrack tip is blunted by the high density of pre-existing, radiation-induced defects in the vicinity of the crack. This interaction leads to the development and growth of a cavity at the interface as opposed to interfacial cleavage.

36 MATERIALS SCIENCE↗

Possibility of reducing nonlinear distortions of a frequency response of a gas ring laser with a periodic alternating frequency bias

We report a theoretical study of the dependence of the beat frequency of counterpropagating waves on the rotation velocity in a gas ring laser with a periodic alternating frequency bias. Rectangular biases of two types are considered: with constant amplitude and with periodic amplitude modulation. We have found dynamic zones for locking-in the counterpropagating wave frequencies and determined deviations of the frequency response from the ideal one that arise when measuring frequency nonreciprocities, close in magnitude to the amplitude of an alternating frequency bias. It is shown that in the case of a rectangular bias with periodic amplitude modulation, nonlinear distortions of the frequency response can be reduced by an order of magnitude. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗