Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Atomic ensemble”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Stiffnesses by (TtN) ensemble molecular dynamics

Calculation of elastic constants requires initializing a lattice of atoms, an h matrix defining the computational cell boundaries, and the parameters defining the interatomic potentials, then numerically integrating the equations of motion for the atoms and the computational cell boundaries, calculating strains and stresses at each timestep, and updating at each timestep the average quantities required for calculation of fluctuations. Enforcements of the DYNAMO FORTRAN code included addition of documentation, improved memory management and data flow, and enabling use of several interatomic potentials. Documentation included headers for each subroutine, which document the parameters input to the routine, algorithms employed, and output.

Tibbits, Patrick↗
DOE OSTI2020

Ensembles of Metastable States Govern Heterogeneous Catalysis on Dynamic Interfaces

Heterogeneous catalysis is at the heart of chemical industry. Being able to tune and design efficient catalysts for processes of interest is of an utmost importance, and for this, the molecular-level understanding of heterogeneous catalysts is the first step, and indeed a prime focus of the modern catalysis research. For a long time, a single most thermodynamically stable structure of the catalytic interface attained in reaction conditions had been envisioned as the reactive phase. However, some catalytic interfaces continue to undergo structural dynamics in the steady state, triggered by high temperatures, pressures, and binding and changing reagents. Among particularly dynamic interfaces are such widely-used catalysts as crystalline and amorphous surfaced supporting (sub-)nano metallic clusters. Recently, it became clear that this dynamic fluxionality causes the supported clusters to populate many distinct structural and stoichiometric states in catalytic conditions. Hence, the catalytic interface should be viewed as an evolving statistical ensemble of many (not one) structures. As such, every member in the ensemble contributes to the properties of the catalyst differently, and in proportion to its probability to be populated. This new notion flips the established paradigm and calls for new theory, modeling approaches, operando measurements, and updated design strategies. The statistical ensemble nature of surface-supported sub-nano cluster catalysts can be exemplified by oxide-supported and adsorbate-covered Pt, Pd, Cu, CuPd clusters, catalytic toward oxidative and non-oxidative dehydrogenation. They have access to a variety of 3D and quasi-2D shapes. The compositions of their thermal ensembles are dependent on the cluster size, leading to size-specific catalytic activities and the famous “every atom counts” phenomenon. The support and adsorbates affect catalyst structures, and state of the reacting species causes the ensemble to change in every reaction intermediate. The most stable member of the ensemble dominates the thermodynamic properties of the corresponding intermediate, whereas the kinetics can be determined by more active but less populated metastable catalyst states, and that suggests that many earlier studies might have overlooked the actual active sites. Both effects depend on the relative time-scales of catalyst restructuring and reaction dynamics. The catalyst may routinely operate off-equilibrium. Ensemble phenomena lead to surprising exceptions from established rules of catalysis, such as scaling relations, and Arrhenius behavior. Catalyst deactivation is also an ensemble property, and its extent of mitigation can be predicted through the new paradigm. These findings were enabled by advances in theory, such as global optimization and subsequent utilization of multiple local minima, and pathways sampling, as well as operando catalyst characterization. The fact that the per-site and per-species resolution is needed for the description and predicting of catalyst properties gives theory the central role in catalysis research, as most experiments provide ensemble-average information and cannot detect the crucial minority species that may be responsible for the catalytic activity.

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

Development of a wide bandwidth heterodyne dispersion interferometer for electron density measurement of atmospheric pressure plasmas

One of the challenges of electron density measurements of an atmospheric pressure plasma (APP) with a laser interferometer is the significant and unwanted phase shift caused by changes in the neutral gas density. These unwanted phase shifts can be mitigated and plasma density measured using an interferometer architecture called a dispersion interferometer (DI). A DI is composed of two nonlinear orientation patterned GaAs crystals for frequency doubling and measures the phase shift induced by a plasma in the interference signal between two second-harmonic beams. Measurement of plasma dynamics or a short plasma pulse in less than a millisecond is enabled with heterodyne detection in a DI using an acousto-optic cell with a frequency of 40 MHz. This heterodyne DI (HDI) is targeted to measure APPs in an electron density range of 10 20 –10 24 m −3 . Finally, the HDI achieves a line-integrated density resolution of 2 × 10 15 m −2 (a phase resolution of 0.005°) with a time constant of 1 μs using ensemble averaging techniques.

atmospheric pressure plasma↗
DOE OSTI2021

Persistent Homology Metrics Reveal Quantum Fluctuations and Reactive Atoms in Path Integral Dynamics

Nuclear quantum effects (NQEs) are known to impact a number of features associated with chemical reactivity and physicochemical properties, particularly for light atoms and at low temperatures. In the imaginary time path integral formalism, each atom is mapped onto a “ring polymer” whose spread is related to the quantum mechanical uncertainty in the particle’s position, i.e., its thermal wavelength. A number of metrics have previously been used to investigate and characterize this spread and explain effects arising from quantum delocalization, zero-point energy, and tunneling. Many of these shape metrics consider just the instantaneous structure of the ring polymers. However, given the significant interest in methods such as centroid molecular dynamics and ring polymer molecular dynamics that link the molecular dynamics of these ring polymers to real time properties, there exists significant opportunity to exploit metrics that also allow for the study of the fluctuations of the atom delocalization in time. Here we consider the ring polymer delocalization from the perspective of computational topology, specifically persistent homology, which describes the 3-dimensional arrangement of point cloud data, (i.e. atomic positions). We employ the Betti sequence probability distribution to define the ensemble of shapes adopted by the ring polymer. The Wasserstein distances of Betti sequences adjacent in time are used to characterize fluctuations in shape, where the Fourier transform and associated principal components provides added information differentiating atoms with different NQEs based on their dynamic properties. We demonstrate this methodology on two representative systems, a glassy system consisting of two atom types with dramatically different de Broglie thermal wavelengths, and ab initio molecular dynamics simulation of an aqueous 4 M HCl solution where the H-atoms are differentiated based on their participation in proton transfer reactions.

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

Spin Squeezing as a Probe of Emergent Quantum Orders

Nuclear magnetic resonance (NMR) experiments can reveal local properties in materials, but are often limited by the low signal-to-noise ratio. Spin squeezed states have an improved resolution below the Heisenberg limit in one of the spin components, and have been extensively used to improve the sensitivity of atomic clocks, for example [1]. Interacting and entangled spin ensembles with non-linear coupling are a natural candidate for implementing squeezing. Furthermore, we propose measurement of the spin-squeezing parameter that itself can act as a local probe of emergent orders in quantum materials. In particular, we demonstrate how to investigate an anisotropic electric field gradient via its coupling to the nuclear quadrupole moment. While squeezed spin states are pure, the squeezing parameter can be estimated for both pure and mixed states. We evaluate the range of fields and temperatures for which a thermal-equilibrium state is sufficient to improve the resolution in an NMR experiment and probe relevant parameters of the quadrupole Hamiltonian, including its anisotropy.

36 MATERIALS SCIENCE↗
DOE OSTI2021

SINGLE: Atomic-resolution structure identification of nanocrystals by graphene liquid cell EM

Analysis of the three-dimensional (3D) structures of nanocrystals with solution-phase transmission electron microscopy is beginning to reveal their unique physiochemical properties. We developed a “one-particle Brownian 3D reconstruction method” based on imaging of ensembles of colloidal nanocrystals using graphene liquid cell electron microscopy. Projection images of differently rotated nanocrystals are acquired using a direct electron detector with high temporal (<2.5 ms) resolution and analyzed to obtain an ensemble of 3D reconstructions. Here, we introduce computational methods required for successful atomic-resolution 3D reconstruction: (i) tracking of the individual particles throughout the time series, (ii) subtraction of the interfering background of the graphene liquid cell, (iii) identification and rejection of low-quality images, and (iv) tailored strategies for 2D/3D alignment and averaging that differ from those used in biological cryo–electron microscopy. Our developments are made available through the open-source software package SINGLE.

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

Conformational ensembles reveal the origins of serine protease catalysis

Enzymes exist in ensembles of states that encode the energetics underlying their catalysis. Conformational ensembles built from 1231 structures of 17 serine proteases revealed atomic-level changes across their reaction states. By comparing the enzymatic and solution reaction, we identified molecular features that provide catalysis and quantified their energetic contributions to catalysis. Serine proteases precisely position their reactants in destabilized conformers, creating a downhill energetic gradient that selectively favors the motions required for reaction while limiting off-pathway conformational states. The same catalytic features have repeatedly evolved in proteases and additional enzymes across multiple distinct structural folds. In conclusion, our ensemble-function analyses revealed previously unknown catalytic features, provided quantitative models based on simple physical and chemical principles, and identified motifs recurrent in nature that may inspire enzyme design.

Du, Siyuan [Stanford Univ., CA (United States)] (O↗
DOE OSTI2023

Reduced dynamic complexity allows structure elucidation of an excited state of KRAS G13D

Localized dynamics of RAS, including regions distal to the nucleotide-binding site, is of high interest for elucidating the mechanisms by which RAS proteins interact with effectors and regulators and for designing inhibitors. Among several oncogenic mutants, methyl relaxation dispersion experiments reveal highly synchronized conformational dynamics in the active (GMPPNP-bound) KRAS G13D , which suggests an exchange between two conformational states in solution. Methyl and 31 P NMR spectra of active KRAS G13D in solution confirm a two-state ensemble interconverting on the millisecond timescale, with a major Pγ atom peak corresponding to the dominant State 1 conformation and a secondary peak indicating an intermediate state different from the known State 2 conformation recognized by RAS effectors. High-resolution crystal structures of active KRAS G13D and KRAS G13D -RAF1 RBD complex provide snapshots of the State 1 and 2 conformations, respectively. We use residual dipolar couplings to solve and cross-validate the structure of the intermediate state of active KRAS G13D , showing a conformation distinct from those of States 1 and 2 outside the known flexible switch regions. The dynamic coupling between the conformational exchange in the effector lobe and the breathing motion in the allosteric lobe is further validated by a secondary mutation in the allosteric lobe, which affects the conformational population equilibrium.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Observation of millihertz-level cooperative Lamb shifts in an optical atomic clock

Collective couplings of atomic dipoles to a shared electromagnetic environment produce a wide range of many-body phenomena. We report on the direct observation of resonant electric dipole-dipole interactions in a cubic array of atoms in the many-excitation limit. The interactions produce spatially dependent cooperative Lamb shifts when spectroscopically interrogating the millihertz-wide optical clock transition in strontium-87. We show that the ensemble-averaged shifts can be suppressed below the level of evaluated systematic uncertainties for optical atomic clocks. Additionally, we demonstrate that excitation of the atomic dipoles near a Bragg angle can enhance these effects by nearly an order of magnitude compared with nonresonant geometries. Our work demonstrates a platform for precise studies of the quantum many-body physics of spins with long-range interactions mediated by propagating photons.

Science & Technology - Other Topics↗
DOE OSTI2026

Ensemble Effects on Hydroxide Bond Dissociation Free Energies in Polyoxovanadate Clusters

Understanding structure-property relationships is foundational to numerous modern chemistries, such as proton-coupled electron transfer (PCET). However, an experimentally measured property is the result of the behavior from an ensemble of molecules. Neglecting ensemble effects, especially under complex chemical environments, may obfuscate these relationships and lead to discrepancies between theory and experiment. In this work, we demonstrate the impact of configurational entropy and local chemical environments on hydroxide bond dissociation free energies [BDFE- (O−H)] for a set of polyoxovanadate nanoclusters, at ambient conditions. The O−H bond strengths are investigated via density functional theory (DFT) coupled with statistical thermodynamic analysis and bilinear modeling, and compared with previous experimental results on the same systems, namely electrochemical solutions of: [V 6 O 13−x (OH) x (TRIOL R ) 2 ] −2 (x = 2, 4, 6; R = NO 2 , Me) and [V 6 O 11−x (OMe) 2 (OH) x (TRIOL NO 2 ) 2 ] −2 (x = 2, 4). Interestingly, we find that ensemble effects, even at room temperature, can account for a significant portion of the BDFE(O−H) trend with the degree of reduction via H atom binding, which cannot be fully captured by single-structure, static DFT calculations. Moreover, we find that the ensemble effects may be replicated statistically, requiring only enumeration of energetically accessible H-binding sites. With the ensemble effects resolved, we present a simple bilinear model to reconcile remaining biases between experiment and ensemble-informed theory, which corelate with clusterspecific electronic environment differences. The bilinear model achieves outstanding accuracy vs experiments with a root-mean squared error of 0.4 kcal/mol. Finally, based on the physicochemical characteristics of hydrogen interaction with polyoxometalates, we present a simple methodology that captures the BDFE(O−H) trend while dramatically reducing required DFT calculations by 98% and achieving accuracy within 1 kcal/mol. Overall, this work elucidates the roles and structural origins of configurational entropy and chemical effects on polyoxometalate hydroxide bond energies, with potential applicability to various atomically precise metal oxide systems. Importantly, it introduces models for rapid and highly accurate property calculations in connection with experiments.

Cluster chemistry↗
DOE OSTI2020

Finite-temperature many-body perturbation theory in the canonical ensemble

Benchmark data are presented for the zeroth- through third-order many-body perturbation corrections to the electronic Helmholtz energy, internal energy, and entropy in the canonical ensemble in a wide range of temperature. They are determined as numerical λ-derivatives of the respective quantities computed by thermal full configuration interaction with a perturbation-scaled Hamiltonian, $\hat{H} = \hat{H}_0 + λ\hat{V}$. Sum-over-states analytical formulas for up to the third-order corrections to these properties are also derived as analytical λ-derivatives. These formulas, which are verified by exact numerical agreement with the benchmark data, are given in terms of the Hirschfelder–Certain degenerate perturbation energies and should be valid for both degenerate and nondegenerate reference states at any temperature down to zero. Further, the results in the canonical ensemble are compared with the same in the grand canonical ensemble.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

KRAS4a and KRAS4b show distinct lipid-dependent regulation of RAS-RAF membrane dynamics

KRAS4a and KRAS4b are important regulators of signaling, and their interactions with the plasma membrane are dynamic and influenced by lipid composition. KRAS 4a and 4b have nearly identical globular domains but differ in their membrane-associated hyper variable region (HVR). The functional distinctions between these isoforms remain unclear, particularly with regards to their dependence on specific lipids and the membrane environment. Previous work showed that the membrane orientation of KRAS4b affects its ability to bind to RAF kinase RBDCRD and that the KRAS–RBDCRD complex adopts different poses on the membrane as well as influences the size and composition of the lipid environment. To model differences between KRAS 4a and 4b protein–lipid interactions, we extended the Multiscale Machine-Learned Modeling Infrastructure (MuMMI) to incorporate continuum simulations in the grand canonical ensemble, enabling sampling across macroscopic, coarse-grained, and all-atom resolutions. Using this framework, we systematically altered PIP2 concentrations, KRAS 4a versus 4b, and RAF RBDCRD complexation to assess impacts on membrane–protein interactions and dynamics. Our results reveal that reducing PIP2 shifts and broadens the membrane orientational preference of both KRAS 4b and 4a, with stronger effects on 4b HVR localization versus 4a. We demonstrate that with depletion of the strong negatively charged PIP2 lipid, the less charged phosphatidylserine replaces PIP2. Our findings highlight similarities and distinctions in the dynamics and lipid dependency of KRAS isoforms and suggest that ordering of the local lipid composition by HVRs is a shared property and key modulator of RAS-mediated signaling at the plasma membrane.

Biological and medical sciences↗
DOE OSTI2022

Enhanced Oxide Reduction by Hydrogen at Cuprous Oxide–Copper Interfaces near Ascending Step Edges

We report that understanding the dynamic processes involved in the interaction of hydrogen with oxides is of fundamental importance in catalysis. This paper probes the reduction of Cu 2 O-‘29’ surfaces by hydrogen at room temperature combining in situ ambient pressure scanning tunneling microscopy (STM), X-ray photoelectron spectroscopy (XPS), and density functional theory (DFT) calculations. Reduction of the atomic layer thin Cu 2 O film is observed to be preferentially initiated2 at step edges and terrace defects, where perfect Cu 2 O(111) terraces are found to be stable towards hydrogenation under the same conditions. After a long induction period, regions of partially reduced Cu 2 O-‘29’ and metallic Cu co-exist before the surface is fully reduced to Cu(111). The reduction rate strongly depends on the nature of nearby Cu step edges. We propose a mechanism for the reduction of Cu 2 O-‘29’ by hydrogen where free copper atoms from ascending metallic step edges facilitate the formation of active ensembles for H 2 dissociation and transfer H to the edges of Cu2O regions

08 HYDROGEN↗
NASA NTRS1988

High altitude atmospheric modeling

Five empirical models were compared with 13 data sets, including both atmospheric drag-based data and mass spectrometer data. The most recently published model, MSIS-86, was found to be the best model overall with an accuracy around 15 percent. The excellent overall agreement of the mass spectrometer-based MSIS models with the drag data, including both the older data from orbital decay and the newer accelerometer data, suggests that the absolute calibration of the (ensemble of) mass spectrometers and the assumed drag coefficient in the atomic oxygen regime are consistent to 5 percent. This study illustrates a number of reasons for the current accuracy limit such as calibration accuracy and unmodeled trends. Nevertheless, the largest variations in total density in the thermosphere are accounted for, to a very high degree, by existing models. The greatest potential for improvements is in areas where we still have insufficient data (like the lower thermosphere or exosphere), where there are disagreements in technique (such as the exosphere) which can be resolved, or wherever generally more accurate measurements become available.

Hedin, Alan E.↗
DOE OSTI2025

Spin squeezing in an ensemble of nitrogen–vacancy centres in diamond

Spin-squeezed states provide a seminal example of how the structure of quantum mechanical correlations can be controlled to produce metrologically useful entanglement. These squeezed states have been demonstrated in a wide variety of quantum systems ranging from atoms in optical cavities to trapped ion crystals. By contrast, despite their numerous advantages as practical sensors, spin ensembles in solid-state materials have yet to be controlled with sufficient precision to generate targeted entanglement such as spin squeezing. Here we report the experimental demonstration of spin squeezing in a solid-state spin system. Our experiments are performed on a strongly interacting ensemble of nitrogen–vacancy colour centres in diamond at room temperature, and squeezing (−0.50 ± 0.13 dB) below the noise of uncorrelated spins is generated by the native magnetic dipole–dipole interaction between nitrogen–vacancy centres. To generate and detect squeezing in a solid-state spin system, we overcome several challenges. First, we develop an approach, using interaction-enabled noise spectroscopy, to characterize the quantum projection noise in our system without directly resolving the spin probability distribution. Second, noting that the random positioning of spin defects severely limits the generation of spin squeezing, we implement a pair of strategies aimed at isolating the dynamics of a relatively ordered sub-ensemble of nitrogen–vacancy centres. Furthermore, our results open the door to entanglement-enhanced metrology using macroscopic ensembles of optically active spins in solids.

Atomic and molecular physics↗
DOE OSTI2021

Formation of a Ti–Cu(111) single atom alloy: Structure and CO binding

A single atom Ti–Cu(111) surface alloy can be generated by depositing small amounts of Ti onto Cu(111) at slightly elevated surface temperatures (~500 to 600 K). Here, scanning tunneling microscopy shows that small Ti-rich islands covered by a Cu single layer form preferentially on ascending step edges of Cu(111) during Ti deposition below about 400 K but that a Ti–Cu(111) alloy replaces these small islands during deposition between 500 and 600 K, producing an alloy in the brims of the steps. Larger partially Cu-covered Ti-containing islands also form on the Cu(111) terraces at temperatures between 300 and 700 K. After surface exposure to CO at low temperatures, reflection absorption infrared spectroscopy (RAIRS) reveals distinct C–O stretch bands at 2102 and 2050 cm –1 attributed to CO adsorbed on Cu-covered Ti-containing domains vs sites in the Ti–Cu(111) surface alloy. Calculations using density functional theory (DFT) suggest that the lower frequency C–O stretch band originates specifically from CO adsorbed on isolated Ti atoms in the Ti–Cu(111) surface alloy and predicts a higher C–O stretch frequency for CO adsorbed on Cu above subsurface Ti ensembles. DFT further predicts that CO preferentially adsorbs in flat-lying configurations on contiguous Ti surface structures with more than one Ti atom and thus that CO adsorbed on such structures should not be observed with RAIRS. The ability to generate a single atom Ti–Cu(111) alloy will provide future opportunities to investigate the surface chemistry promoted by a representative early transition metal dopant on a Cu(111) host surface.

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

RASPA3

RASPA3, a molecular simulation code for computing adsorption and diffusion in nanoporous materials and thermodynamic and transport properties of fluids. It implements force field based classical Monte Carlo/molecular dynamics in various ensembles. RASPA3 is rewritten from the ground up in C++23 with speed and code readability in mind. Transition-matrix Monte Carlo is added to compute the density of states and free energies. The Monte Carlo code for rigid molecules is based on quaternions, and the atomic positions needed in the energy evaluation are recreated from the center of mass position and quaternion orientation. The expanded ensemble methodology for fractional molecules, with a scaling parameter λ between 0 and 1, now also keeps track of analytic expressions of dU/dλ, allowing independent verification of the chemical potential using thermodynamic integration. The source code is freely available under the MIT license on GitHub.

Dubbeldam, David↗
DOE OSTI2023

Adaptive Pore Opening to Form Tailored Adsorption Sites in a Cooperatively Flexible Framework Enables Record Inverse Propane/Propylene Separation

A proposed low-energy alternative to the separation of alkanes from alkenes by energy-intensive cryogenic distillation is separation by porous adsorbents. Unfortunately, most adsorbents preferentially take up the desired, high-value major component alkene, requiring frequent regeneration. Adsorbents with inverse selectivity for the minor component alkane would enable the direct production of purified, reagent-grade alkene, greatly reducing global energy consumption. However, such materials are exceedingly rare, especially for propane/propylene separation. Here, we report that through adaptive and spontaneous pore size and shape adaptation to optimize an ensemble of weak noncovalent interactions, the structurally responsive metal-organic framework CdIF-13 (sod-Cd(benzimidazolate) 2 ) exhibits inverse selectivity for propane over propylene with record-setting separation performance under industrially relevant temperature, pressure, and mixture conditions. Powder synchrotron X-ray diffraction measurements combined with first-principles calculations yield atomic-scale insight and reveal the induced fit mechanism of adsorbate-specific pore adaptation and ensemble interactions between ligands and adsorbates. Dynamic column breakthrough measurements confirm that CdIF-13 displays selectivity under mixed-component conditions of varying ratios, with a record measured selectivity factor of α ≈ 3 at 95:5 propylene:propane at 298 K and 1 bar. When sequenced with a low-cost rigid adsorbent, we demonstrated the direct purification of propylene under ambient conditions. In conclusion, this combined atomic-level structural characterization and performance testing firmly establishes how cooperatively flexible materials can be capable of unprecedented separation factors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗