Search NASA⌕ Search

SEARCH · Search NASA

Results for “Molecular statics”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Excited State Bond Homolysis of Vanadium(V) Photocatalysts for Alkoxy Radical Generation

Advancements in photocatalysis have transformed synthetic organic chemistry, using light as a powerful tool to drive selective chemical transformations. Recent approaches have focused on metal-halide ligand-to-metal charge transfer (LMCT) photoactivated bond homolysis reactions leveraged by earth-abundant elements to generate valuable synthons for radical-mediated cross-coupling reactions. Of recent utility, oxovanadium(V) LMCT photocatalysts exhibit selective alkoxy radical generation from aliphatic alcohols upon blue light (UVA) irradiation under mild conditions. The selective photochemical liberation of alkoxy radicals is valuable for applying late-stage fragmentation approaches in organic synthesis and depolymerization strategies for nonbiodegradable polymers. Steady-state and time-resolved spectroscopy were used to assign the electronic structure of three well-defined V(V) photocatalysts in their ground and excited states. We assign the excited state for this transformation at earth-abundant vanadium(V), interrogating the electronic structure using static UV–visible absorption, ultrafast transient absorption, and electron paramagnetic resonance spectroscopy coupled to computational approaches. Furthermore, these findings afford assignments of the short-lived excited state intermediates that dictate selective homolytic bond cleavage in metal alkoxides, illustrating the valuable insight gleaned from fundamental investigations of the molecular photochemistry responsible for light-escalated chemical transformations.

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

High-pressure phase transition in 3-D printed nanolamellar high-entropy alloy by imaging and simulation insights

Abstract We report on the high-resolution imaging and molecular dynamics simulations of a 3D-printed eutectic high-entropy alloy (EHEA) Ni 40 Co 20 Fe 10 Cr 10 Al 18 W 2 consisting of nanolamellar BCC and FCC phases. The direct lattice imaging of 3D-printed samples shows the Kurdjumov–Sachs (K–S) orientation relation {111} FCC parallel to {110} BCC planes in the dual-phase lamellae. Unlike traditional iron and steels, this alloy shows an irreversible BCC-to-FCC phase transformation under high pressures. The nanolamellar morphology is maintained after pressure cycling to 30 GPa, and nano-diffraction studies show both layers to be in the FCC phase. The chemical compositions of the dual-phase lamellae after pressure recovery remain unchanged, suggesting a diffusion-less BCC–FCC transformation in this EHEA. The lattice imaging of the pressure-recovered sample does not show any specific orientation relation between the two resulting FCC phases, indicating that many grain orientations are produced during the BCC–FCC phase transformation. Molecular dynamics simulations on phase transformation in a nanolamellar BCC/FCC in K–S orientation show that phase transformation from BCC to FCC is completed under high pressures, and the FCC phase is retained on decompression aided by the stable interfaces. Our work elucidates the irreversible phase transformation under static compression, providing an understanding of the orientation relationships in 3-D printed EHEA under high pressures.

3-D Printing↗
DOE OSTI2024

Quantum tensor network algorithms for evaluation of spectral functions on quantum computers

We investigate quantum algorithms derived from tensor networks to simulate the static and dynamic properties of quantum many-body systems. Using a sequentially prepared quantum circuit representation of a matrix product state (MPS) that we call a quantum tensor network (QTN), we demonstrate algorithms to prepare ground and excited states on a quantum computer and apply them to molecular nanomagnets (MNMs) as a paradigmatic example. In this setting, we develop two approaches for extracting the spectral correlation functions measured in neutron-scattering experiments: (a) a generalization of the SWAP test for computing wave function overlaps and, (b) a generalization of the notion of matrix product operators to the QTN setting which generates a linear combination of unitaries. The latter method is discussed in detail for translationally invariant spin-half systems, where it is shown to reduce the qubit resource requirements compared with the SWAP method and may be generalized to other systems. We demonstrate the versatility of our approaches by simulating spin-1/2 and spin-3/2 MNMs, with the latter being an experimentally relevant model of a Cr$^{3+}_{8}$ ring. Here, our approach has qubit requirements that are independent of the number of constituents of the many-body system and scale only logarithmically with the bond dimension of the MPS representation, making them appealing for implementation on near-term quantum hardware with mid-circuit measurement and reset.

Neutron scattering↗
DOE OSTI2020

Coordinate transformation methodology for simulating quasistatic elastoplastic solids

Molecular dynamics simulations frequently employ periodic boundary conditions where the positions of the periodic images are manipulated in order to apply deformation to the material sample. For example, Lees-Edwards conditions use moving periodic images to apply simple shear. Here, we examine the problem of precisely comparing this type of simulation to continuum solid mechanics. We employ a hypo-elastoplastic mechanical model, and develop a projection method to enforce quasi-static equilibrium. We introduce a simulation framework that uses a fixed Cartesian computational grid on a reference domain, and imposes deformation via a time-dependent coordinate transformation to the physical domain. As a test case for our method, we consider the evolution of shear bands in a bulk metallic glass using the shear transformation zone theory of amorphous plasticity. We examine the growth of shear bands in simple shear and pure shear conditions as a function of the initial preparation of the bulk metallic glass.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

A Modified Network Alteration Model to Predict Quasi-Static Behavior of the Cross-Linked Polymers During Hydrolytic Aging

Abstract The network alteration theory is utilized to capture the effects of hydrolytic aging on the mechanical integrity of the polymer matrix. Hydrolytic degradation generally considered as a result of two simultaneous chemical phenomena; (i) reduction of the cross-links, and (ii) active chain loss. Assuming the alteration of the chain density is identical to the change of the number average molecular weight along the aging trajectory, the strain energy of polymer matrix can be rewritten as a function of deformation, deformation history, storage time, and aging temperature. The proposed model includes only five physically inspired material parameters. Thus, while it is computationally efficient, it shows good agreement with own experimental data, which performed on various range of accelerated aging temperatures and times. With respect to its computational efficiency, simplicity, accuracy, and interpret-ability, the model is a good choice for advanced implementations in FE programs.

Bahrololoumi, Amir↗
DOE OSTI2020

Multireference Electron Correlation Methods: Journeys along Potential Energy Surfaces

Multireference electron correlation methods describe static and dynamical electron correlation in a balanced way and, therefore, can yield accurate and predictive results even when single-reference methods or multiconfigurational self-consistent field theory fails. One of their most prominent applications in quantum chemistry is the exploration of potential energy surfaces. This includes the optimization of molecular geometries, such as equilibrium geometries and conical intersections and on-the-fly photodynamics simulations, both of which depend heavily on the ability of the method to properly explore the potential energy surface. Because such applications require nuclear gradients and derivative couplings, the availability of analytical nuclear gradients greatly enhances the scope of quantum chemical methods. This review focuses on the developments and advances made in the past two decades. A detailed account of the analytical nuclear gradient and derivative coupling theories is presented. Emphasis is given to the software infrastructure that allows one to make use of these methods. Notable applications of multireference electron correlation methods to chemistry, including geometry optimizations and on-the-fly dynamics, are summarized at the end followed by a discussion of future prospects.

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

Shock compression of crystalline TeO 2 to the high-pressure fluid regime: Insights from ab initio molecular dynamics simulations

The shock response of fully-dense and porous crystalline tellurium dioxide (TeO 2 ⁠) to the high-pressure and high-temperature fluid regime was investigated within the framework of density functional theory with Mermin’s generalization to finite temperatures. The principal and porous shock Hugoniot curves were predicted from canonical ab initio molecular dynamics (AIMD) simulations, with the phase space sampled along isotherms up to 80 000 K, for densities ranging from ρ = 3 to 17 g/cm 3 . The polymorphs investigated are α-TeO 2 paratellurite (⁠P4 1 2 1 2), TeO 2 cotunnite (⁠Pnma⁠), and TeO 2 post-cotunnite (⁠P2 1 /m⁠). Based on the discontinuity found in the calculated U s – u p slope of TeO 2 post-cotunnite at a shock velocity of U s ≃ 8.35km/s and a particle velocity of u p ≃ 3.64 km/s, the shock melting temperature and pressure are predicted to be ≃ 6500 K and ≃ 170 GPa. Results from the AIMD simulations are in line with the static compression data of TeO 2 paratellurite and cotunnite, and with the recent shock Hugoniot data for single-crystal α- TeO 2 for pressures up to 85 GPa, obtained using the inclined-mirror method and the velocity interferometer system for any reflector combined with powder gun and two-stage light-gas gun.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Trigonal Bipyramidal V 3+ Complex as an Optically Addressable Molecular Qubit Candidate

Synthetic chemistry enables a bottom-up approach to quantum information science, where atoms can be deterministically positioned in a quantum bit or qubit. Two key requirements to realize quantum technologies are qubit initialization and read-out. By imbuing molecular spins with optical initialization and readout mechanisms, analogous to solid-state defects, molecules could be integrated into existing quantum infrastructure. To mimic the electronic structure of optically addressable defect sites, we designed the spin-triplet, V 3+ complex, (C 6 F 5 ) 3 trenVCN t Bu (1).We measured the static spin properties as well as the spin coherence time of1demonstrating coherent control of this spin qubit with a 240 GHz electron paramagnetic resonance spectrometer powered by a free electron laser. We found that1exhibited narrow, near-infrared photoluminescence (PL) from a spin-singlet excited state. Using variable magnetic field PL spectroscopy, we resolved emission into each of the ground-state spin sublevels, a crucial component for spin-selective optical initialization and readout. This work demonstrates that trigonally symmetric, heteroleptic V 3+ complexes are candidates for optical spin addressability.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2021

Evidence of a liquid–liquid transition in a glass-forming ionic liquid

A liquid–liquid transition (LLT) is a transformation from one liquid to another through a first-order transition. The LLT is fundamental to the understanding of the liquid state and has been reported in a few materials such as silicon, phosphorus, triphenyl phosphite, and water. Furthermore, it has been suggested that the unique properties of materials such as water, which is critical for life on the planet, are linked to the existence of the LLT. However, the experimental evidence for the existence of an LLT in many molecular liquids remains controversial, due to the prevalence and high propensity of the materials to crystallize. In this work, we show evidence of an LLT in a glass-forming trihexyltetradecylphosphonium borohydride ionic liquid that shows no tendency to crystallize under normal laboratory conditions. We observe a step-like increase in the static dielectric permittivity at the transition. Furthermore, the sizes of nonpolar local domains and ion-coordination numbers deduced from wide-angle X-ray scattering also change abruptly at the LLT. We independently corroborate these changes in local organization using Raman spectroscopy. The experimental access to the evolution of local order and structural dynamics across a liquid–liquid transition opens up unprecedented possibilities to understand the nature of the liquid state.

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

Sweep-tracing algorithm: in silico slip crystallography and tension-compression asymmetry in BCC metals

Abstract Direct Molecular Dynamics (MD) simulations are being increasingly employed to model dislocation-mediated crystal plasticity with atomic resolution. Thanks to the dislocation extraction algorithm (DXA), dislocation lines can be now accurately detected and positioned in space and their Burgers vector unambiguously identified in silico, while the simulation is being performed. However, DXA extracts static snapshots of dislocation configurations that by themselves present no information on dislocation motion. Referred to as a sweep-tracing algorithm (STA), here we introduce a practical computational method to observe dislocation motion and to accurately quantify its important characteristics such as preferential slip planes (slip crystallography). STA reconnects pairs of successive snapshots extracted by DXA and computes elementary slip facets thus precisely tracing the motion of dislocation segments from one snapshot to the next. As a testbed for our new method, we apply STA to the analysis of dislocation motion in large-scale MD simulations of single crystal plasticity in BCC metals. We observe that, when the crystal is subjected to uniaxial deformation along its [001] axis, dislocation slip predominantly occurs on the {112} maximum resolved shear stress plane under tension, while in compression slip is non-crystallographic (pencil) resulting in asymmetric mechanical response. The marked contrast in the observed slip crystallography is attributed to the twinning/anti-twinning asymmetry of shears in the {112} planes relatively favoring dislocation motion in the twinning sense while hindering dislocations from moving in the anti-twinning directions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

AIMD‐Based Protocols for Modeling Exciplex Fluorescence Spectra and Inter‐System Crossing in Photocatalytic Chromophores

ABSTRACT This study introduces a computational protocol for modeling the emission spectra of exciplexes using excited‐state ab initio molecular dynamics (AIMD) simulations. The protocol is applied to a model exciplex formed by oligo‐p‐phenylenes (OPPs) and triethylamine (TEA), which is of interest in the context of photocatalytic reduction of . AIMD facilitates efficient sampling of the conformational space of OPP3 and OPP4 exciplexes with TEA, offering a dynamic alternative to previously employed static methods. The AIMD‐based protocol successfully reproduces experimental emission spectra for OPP‐TEA exciplexes, agreeing with previous computational and experimental findings. The results show that AIMD simulations provide an efficient means of sampling the conformational space of these exciplexes, requiring less user input and, in some instances, fewer computational resources than multiple excited‐state optimizations initiated from user‐specified initial structures. The study also evaluates the yield of intersystem crossing (ISC) using AIMD and Landau‐Zener probability. The results suggest that ISC is a minor decay channel for OPP3 and OPP4. This work provides new insights into the structural flexibility and emission characteristics of OPP‐TEA photoredox catalyst systems, potentially contributing to improved design strategies for organic chromophores in reduction applications.

Giudetti, Goran [Department of Chemistry Universit↗
DOE OSTI2023

Simulating dielectric spectra: A demonstration of the direct electric field method and a new model for the nonlinear dielectric response

Here we demonstrate a method to compute the dielectric spectra of fluids in molecular dynamics (MD) by directly applying electric fields to the simulation. We obtain spectra from MD simulations with low magnitude electric fields (≈0.01 V/Å) in agreement with spectra from the fluctuation–dissipation method for water and acetonitrile. We examine this method’s trade-off between noise at low field magnitudes and the nonlinearity of the response at higher field magnitudes. We then apply the Booth equation to describe the nonlinear response of both fluids at low frequency (0.1 GHz) and high field magnitude (up to 0.5 V/Å). We develop a model of the frequency-dependent nonlinear response by combining the Booth description of the static nonlinear dielectric response of fluids with the frequency-dependent linear dielectric response of the Debye model. We find good agreement between our model and the MD simulations of the nonlinear dielectric response for both acetonitrile and water.

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

A spall and diffraction study of nanosecond pressure release across the iron ε-α phase boundary

The extreme response of polycrystalline iron at high pressures and high strain rates is revealed by means of high-power laser pulses. The compression portion of the pulse coupled with x-ray diffraction identifies the expected body-centered cubic (α) to hexagonal close packed (ε) displacive transformation. Upon release, observation shows that the complete reverse transformation takes approximately 8 ns and that the structure returns to its initial microstructural configuration, in a reversible transformation path. This is in good agreement with molecular dynamics (MD) simulations which predict an inverse dependence between transformation time and strain rate. The grain size is reduced from μm to nm range during compression and begins increasing back to the original grain size on decompression. The kinetics of the transition is dictated by heterogenous nucleation as it follows the Johnson-Mehl-Avrami-Kolmogorov equation with the appropriate time exponent of ~1. This is confirmed by MD simulations which also identify profuse twinning and dislocation generation. Here, the tensile pulse generated upon reflection at the free surface is captured by time-resolved free surface velocity measurements from which a peak tensile stress of 7 GPa is obtained, in stark contrast with its quasi-static value of ~200 MPa. At these strain rates, the strength of grain interiors, which is determined by twinning and slip exceeds the strength of the boundaries, and failure initiates preferentially in the latter.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Solar-driven selective conversion of millimolar dissolved carbon to fuels with molecular flux generation

Abstract The direct utilization of dissolved inorganic carbon in seawater for CO 2 conversion promises chemical production on-demand and with zero carbon footprint. Photoelectrochemical (PEC) CO 2 reduction (CO 2 R) devices promise the sustainable conversion of dissolved carbon in seawater to carbon products using sunlight as the only energy input. However, the diffusion-dominant transport mechanism and the near-zero concentration of CO 2 (aq) (CO 2 dissolved in aqueous solution) in static seawater has made it extremely challenging to achieve high solar-to-fuel (STF) efficiency and high carbon-product selectivity. Here, where CO 2 (aq) as a reactant generated in situ by acidification of HCO 3 - flows continuously from BiVO 4 photoanodes to Si photocathodes, enabling a single-step conversion of dissolved carbon into products. Our PEC device significantly increases the CO selectivity from 3% to 21%, which approaches the 30% theoretical limit according to multi-physics modeling. Meanwhile, the Si/BiVO 4 PEC CO 2 R device achieved a STF efficiency of 0.71%. Such flow engineering achieves flow-dependent selectivity, rate, and stability in simulated seawater, thus promising practical solar fuel production at scale.

Liu, Bin↗
DOE OSTI2021

Effects of interaction strength of associating groups on linear and star polymer dynamics

A small number of associating groups incorporated onto a polymer backbone have dramatic effects on the mobility and viscoelastic response of the macromolecules in melts. These associating groups assemble, driving the formation of clusters, whose lifetime affects the properties of the polymers. In this study, we probe the effects of the interaction strength on the structure and dynamics of two topologies, linear and star polymer melts, and further investigate blends of associative and non-associating polymers using molecular dynamics simulations. Polymer chains of approximately one entanglement length are described by a bead–spring model, and the associating groups are incorporated in the form of interacting beads with an interaction strength between them that is varied from 1 to 20 k B T. We find that, for all melts and blends, interaction of a few k B T between the associating groups drives cluster formation, where the size of the clusters increases with increasing interaction strength. These clusters act as physical crosslinkers, which slow the chain mobility. Blends of chains with and without associating groups macroscopically phase separate for interaction strength between the associating groups of a few k B T and above. For weakly interacting associating groups, the static structure function S(q) is well fit by functional form predicted by the random phase approximation where a clear deviation occurs as phase segregation takes place, providing a quantitative assessment of phase segregation.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Gas Transport in the Insensitive High Explosive PBX 9502

Small-scale laboratory experiments were performed to analyze gas transport mechanics in PBX 9502, an insensitive high explosive (IHE) composition. Two independent investigations are reported here. First, gas permeametry and dynamic pycnometry techniques were used to measure the molecular flow (Knudsen) coefficient and the internally interconnected void fraction of pristine and thermally damaged samples in two uniaxial pressing orientations. The permeability of PBX 9502 was found to be unmeasurably low (of the same order of magnitude as PTFE) with gas transport being diffusion-dominated. Secondly, a pressure vessel experiment was developed to measure quasi-static and dynamic gas generation as the explosive was heated to self-ignition (cookoff). The gas generation results and the permeability/diffusivity findings provide evidence that PBX 9502 remains impermeable until seconds prior to self-ignition. At ignition, internal void-pressure drives macro-scale cracking and the sample becomes uniformly incorporated in the following deconsolidating deflagration. These results are discussed within the context of previous observations of pressure-dependent cookoff behavior of PBX 9502 and provide a more complete description of thermal damage evolution in this explosive composition.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2024

Quantitative measurements of reactive oxygen species partitioning in electron transfer flavoenzyme magnetic field sensing

Biological magnetic field sensing that gives rise to physiological responses is of considerable importance in quantum biology. The radical pair mechanism (RPM) is a fundamental quantum process that can explain some of the observed biological magnetic effects. In magnetically sensitive radical pair (RP) reactions, coherent spin dynamics between singlet and triplet pairs are modulated by weak magnetic fields. The resulting singlet and triplet reaction products lead to distinct biological signaling channels and cellular outcomes. A prevalent RP in biology is between flavin semiquinone and superoxide (O 2 •–) in the biological activation of molecular oxygen. This RP can result in a partitioning of reactive oxygen species (ROS) products to form either O 2 •– or hydrogen peroxide (H 2 O 2 ). Here, we examine magnetic sensing of recombinant human electron transfer flavoenzyme (ETF) reoxidation by selectively measuring O 2 •– and H 2 O 2 product distributions. ROS partitioning was observed between two static magnetic fields at 20 nT and 50 μT, with a 13% decrease in H 2 O 2 singlet products and a 10% increase in O 2 •– triplet products relative to 50 µT. RPM product yields were calculated for a realistic flavin/superoxide RP across the range of static magnetic fields, in agreement with experimental results. For a triplet born RP, the RPM also predicts about three times more O 2 •– than H 2 O 2 , with experimental results exhibiting about four time more O 2 •– produced by ETF. The method presented here illustrates the potential of a novel magnetic flavoprotein biological sensor that is directly linked to mitochondria bioenergetics and can be used as a target to study cell physiology.

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

In situ X-ray diffraction of silicate liquids and glasses under dynamic and static compression to megabar pressures

Properties of liquid silicates under high-pressure and high-temperature conditions are critical for modeling the dynamics and solidification mechanisms of the magma ocean in the early Earth, as well as for constraining entrainment of melts in the mantle and in the present-day core–mantle boundary. In this paper, we present in situ structural measurements by X-ray diffraction of selected amorphous silicates compressed statically in diamond anvil cells (up to 157 GPa at room temperature) or dynamically by laser-generated shock compression (up to 130 GPa and 6,000 K along the MgSiO 3 glass Hugoniot). The X-ray diffraction patterns of silicate glasses and liquids reveal similar characteristics over a wide pressure and temperature range. Beyond the increase in Si coordination observed at 20 GPa, we find no evidence for major structural changes occurring in the silicate melts studied up to pressures and temperatures exceeding Earth’s core mantle boundary conditions. This result is supported by molecular dynamics calculations. Our findings reinforce the widely used assumption that the silicate glasses studies are appropriate structural analogs for understanding the atomic arrangement of silicate liquids at these high pressures.

58 GEOSCIENCES↗