Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Quantum chemical 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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Energy Frontier Research Centers: Center for the Computational Design of Functional Layered Materials (CCDM) August 1, 2014 - July 31, 2018; Center for Complex Materials from First Principles (CCM) August 1, 2018 - July 31, 2021 (Final Report)

The mission of the DOE Energy Frontier Research Centers CCDM (2014-2018) and CCM (2018-2021) was to theoretically develop, computationally apply, and experimentally validate electronic structure methods for all materials, with a focus on the complex materials, especially layered and two-dimensional materials, strongly-correlated materials, and liquid water. This was achieved by over 200 published journal articles authored by about 17 senior investigators from physics and chemistry and from theory, computation, and experiment, plus their collaborators. In particular, the Centers confirmed the predictive power of the SCAN (strongly constrained and appropriately normed) density functional, which was constructed to satisfy 17 known exact constraints and several appropriate norms. Without being fitted to real bonded systems, and at a modest computational cost, SCAN correctly predicted covalent, ionic, metallic, hydrogen, and van der Waals bonds in many challenging materials. SCAN gave an improved description of defects in semiconductors, surface properties of metals, seven phases of ice, liquid water, liquid and supercooled silicon, subtle structural distortions in ferroelectrics, formation energies and structural predictions for solids, and critical pressures for structural phase transitions. Perhaps most remarkably, SCAN correctly described some strongly-correlated materials that were previously believed to be beyond the reach of density-functional approximations. SCAN is the only density functional that correctly predicts the band gap closing under chemical doping of the cuprate high-temperature superconducting materials. SCAN also predicts a landscape of competing stripe and magnetic phases in the cuprates. For some materials with some codes, SCAN has convergence problems that are greatly reduced by the CCM-developed r 2 SCAN, without loss of accuracy or rigor. SCAN and r 2 SCAN still make some self-interaction error, which is greatly reduced by the CCDM/ CCM-developed local orbital scaling correction (LOSC). These Centers further proved that the fundamental energy gaps of a solid from an orbital energy difference and from total energy differences are the same for a large class of generalized Kohn-Sham (GKS) functionals, including SCAN and standard hybrid functionals, and that symmetry breaking arises when a dynamic density fluctuation drops to zero frequency. The Centers identified new mechanisms for catalysis in layered materials with ions intercalated between the layers, investigated charge density waves both in model systems and in real layered materials, studied changes of band gap with the number of layers, and explored topological ultrathin films, bent nanoribbons, and defects.

08 HYDROGEN↗
DOE OSTI2023

Optically Enhanced Solid-State 1 H NMR Spectroscopy

Low sensitivity is the primary limitation to extending nuclear magnetic resonance (NMR) techniques to more advanced chemical and structural studies. Photochemically induced dynamic nuclear polarization (photo-CIDNP) is an NMR hyperpolarization technique where light is used to excite a suitable donor–acceptor system, creating a spin-correlated radical pair whose evolution drives nuclear hyperpolarization. Systems that exhibit photo-CIDNP in solids are not common, and this effect has, up to now, only been observed for 13 C and 15 N nuclei. However, the low gyromagnetic ratio and natural abundance of these nuclei trap the local hyperpolarization in the vicinity of the chromophore and limit the utility for bulk hyperpolarization. Here, we report the first example of optically enhanced solid-state 1 H NMR spectroscopy in the high-field regime. This is achieved via photo-CIDNP of a donor–chromophore–acceptor molecule in a frozen solution at 0.3 T and 85 K, where spontaneous spin diffusion among the abundant strongly coupled 1 H nuclei relays polarization through the whole sample, yielding a 16-fold bulk 1 H signal enhancement under continuous laser irradiation at 450 nm. These findings enable a new strategy for hyperpolarized NMR beyond the current limits of conventional microwave-driven DNP.

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

Zero-Field NMR and Millitesla-SLIC Spectra for >200 Molecules from Density Functional Theory and Spin Dynamics

NMR is usually performed at magnetic fields of 1 T and above to obtain sufficient sensitivity and spectral dispersion to identify chemicals based on chemical shifts and J couplings. At lower fields, the advent of hyperpolarization technologies and sensitive detectors can address sensitivity concerns. However, it remains disputed whether spectral signatures at zero and ultra-low fields are sufficient for chemical identification. Here, we report an all–electron DFT-based batch calculation of J-coupling constants, which are used to generate J coupling NMR spectra at zero field and 6.5 mT for over 200 small molecules. In the developed computational tool chain, we first used the all-electron FHI-aims code to calculate the molecular J couplings and chemical shifts. We then fed the calculated NMR parameters into the NMR simulation package SPINACH to simulate both heteronuclear J coupling spectra at zero-field, and homonuclear J coupling spectra as spin-lock induced crossing (SLIC) spectra at ultra-low field (6.5 mT). The resulting spectra demonstrate that zero and ultra-low field NMR spectra can represent unique identifiers of chemical structure for small molecules.

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

Importance of Electron Correlation on the Geometry and Electronic Structure of [2Fe–2S] Systems: A Benchmark Study of the [Fe 2 S 2 (SCH 3 ) 4 ] 2–,3–,4– , [Fe 2 S 2 (SCys) 4 ] 2– , [Fe 2 S 2 (S- p -tol) 4 ] 2– , and [Fe 2 S 2 (S- o -xyl) 4 ] 2– Complexes

Iron–sulfur clusters are crucial for biological electron transport and catalysis. Obtaining accurate geometries, energetics, manifolds of their excited electronic states, and reduction energies is important to understand their role in these processes. Here, using a [2Fe–2S] model complex with Fe II and Fe III oxidation states, which leads to different charges, i.e., [Fe 2 S 2 (SMe) 4 ] 2–,3–,4– , we benchmarked a variety of computational methodologies ranging from density functional theory (DFT) to post-Hartree–Fock methods, including complete active space self-consistent field (CASSCF), multireference configuration interaction, the second-order N-electron valence state perturbation theory (NEVPT2), and the linearized integrand approximation of adiabatic connection (AC0) approaches. Additionally, we studied three experimentally well-characterized complexes, [Fe 2 S 2 (SCys) 4 ] 2– , [Fe 2 S 2 (S-o-tol) 4 ] 2– , and [Fe 2 S 2 (S-o-xyl) 4 ] 2– , via DFT methods. We conclude that the dynamic electron correlation is important for accurately predicting the geometry of these complexes. Broken symmetry (BS) DFT correctly predicts experimental geometries of low-spin multiplicity, while CASSCF does not. However, BS-DFT significantly overestimates the difference between the low- and high-spin electronic states for a given oxidation state. At the same time, CASSCF underestimates it but provides relative energies closer to the reference NEVPT2 results. Finally, AC0 provides energetics of NEVPT2 quality with the additional advantage of being able to use large CASSCF sizes. NEVPT2 gives the best estimates of the Fe III /Fe III → Fe II /Fe III (4.27 eV) and Fe II /F III → Fe II /F II (7.72 eV) reduction energies. The results provide insight into the electronic structure of these complexes and assist in the understanding of their physical properties.

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

Magnon-Magnon Interaction Induced by Dynamic Coupling in a Hybrid Magnonic Crystal

We report a combined experimental and numerical investigation of spin-wave dynamics in a hybrid magnonic crystal consisting of a CoFeB artificial spin ice (ASI) of stadium-shaped nanoelements patterned atop a continuous NiFe film separated by a 5 nm Al 2 O 3 spacer. Using Brillouin light scattering spectroscopy, we probe the frequency dependence of thermal spin waves as functions of applied magnetic field and wavevector, revealing the decisive role of interlayer dipolar coupling in the magnetization dynamics. Micromagnetic simulations complement the experiments, showing a strong interplay between ASI edge modes and backward volume modes in the NiFe film. The contrast in saturation magnetization between CoFeB and NiFe enhances this coupling, leading to a pronounced hybridization manifested as a triplet of peaks in the BLS spectra predicted by simulations and observed experimentally. This magnon−magnon coupling persists over a wide magnetic field range, shaping both the spin-wave dispersion at fixed fields and the full frequency-field response throughout the magnetic hysteresis loop. Our findings establish how ASI geometry can selectively enhance specific spin-wave wavelengths in the underlying film, thereby boosting their amplitude and identifying them as preferential channels for spin wave transmission and manipulation.

Brillouin light scattering spectroscopy↗
DOE OSTI2025

Observation of the axion quasiparticle in 2D MnBi 2 Te 4

The axion is a hypothetical fundamental particle that is conjectured to correspond to the coherent oscillation of the θ field in quantum chromodynamics. Its existence would solve multiple fundamental questions, including the strong CP problem of quantum chromodynamics and dark matter, but the axion has never been detected. Electrodynamics of condensed-matter systems can also give rise to a similar θ, so far studied as a static, quantized value to characterize the topology of materials. Coherent oscillation of θ in condensed matter has been proposed to lead to physics directly analogous to the high-energy axion particle—the dynamical axion quasiparticle (DAQ). Here we report the observation of the DAQ in MnBi 2 Te 4 . By combining a two-dimensional electronic device with ultrafast pump–probe optics, we observe a coherent oscillation of θ at about 44 gigahertz, which is uniquely induced by its out-of-phase antiferromagnetic magnon. This represents direct evidence for the presence of the DAQ, which in two-dimensional MnBi 2 Te 4 is found to arise from the magnon-induced coherent modulation of the Berry curvature. The DAQ also has implications in light–matter interaction and coherent antiferromagnetic spintronics, as it might lead to axion polaritons and electric control of ultrafast spin polarization. Finally, the DAQ could be used to detect axion particles. We estimate the detection frequency range and sensitivity in the millielectronvolt regime, which has so far been poorly explored.

topological matter↗
DOE OSTI2021

Predicted Reaction Mechanisms, Product Speciation, Kinetics, and Detonation Properties of the Insensitive Explosive 2,6-Diamino-3,5-dinitropyrazine-1-oxide (LLM-105)

2,6-Diamino-3,5-dinitropyrazine-1-oxide (LLM-105) is a relatively new and promising insensitive high-explosive (IHE) material that remains only partially characterized. IHEs are of interest for a range of applications and from a fundamental science standpoint, as the root causes behind insensitivity are poorly understood. In this work, we adopt a multitheory approach based on reactive molecular dynamic simulations performed with density functional theory, density functional tight-binding, and reactive force fields to characterize the reaction pathways, product speciation, reaction kinetics, and detonation performance of LLM-105. We compare and contrast these predictions to 1,3,5-triamino-2,4,6-trinitrobenzene (TATB), a prototypical IHE, and 1,3,5,7-tetranitro-1,3,5,7-tetrazoctane (HMX), a more sensitive and higher performance material. The combination of different predictive models allows access to processes operative on progressively longer timescales while providing benchmarks for assessing uncertainties in the predictions. We find that the early reaction pathways of LLM-105 decomposition are extremely similar to TATB; they involve intra- and intermolecular hydrogen transfer. Additionally, the detonation performance of LLM-105 falls between that of TATB and HMX. We find agreement between predictive models for first-step reaction pathways but significant differences in final product formations. Predictions of detonation performance result in a wide range of values, and one-step kinetic parameters show the similar reaction rates at high temperatures for three out of four models considered.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Identifying the Molecular Edge Termination of Exfoliated Hexagonal Boron Nitride Nanosheets with Solid-State NMR Spectroscopy and Plane-Wave DFT Calculations

Hexagonal boron nitride nanosheets (h-BNNS), the isoelectronic analog to graphene, have received interest over the past decade due to their high thermal oxidative resistance, high bandgap, catalytic activity, and low cost. The functional groups that terminate boron and nitrogen zigzag and/or armchair edges directly affect their chemical, physical, and electronic properties. However, an understanding of the molecular edge termination present in h-BNNS is lacking. Here, high-resolution magic-angle spinning (MAS) solid-state NMR (SSNMR) spectroscopy, and plane-wave density-functional theory (DFT) calculations are used to determine the molecular edge termination in exfoliated h-BNNS. 1 H → 11 B cross-polarization MAS (CPMAS) SSNMR spectra of h-BNNS revealed multiple hydroxyl/oxygen coordinated boron edge sites that were not detectable in direct excitation experiments. A dynamic nuclear polarization (DNP)-enhanced 1 H → 15 N CPMAS spectrum of h-BNNS displayed four distinct 15 N resonances while a 2D 1 H{ 14 N} dipolar-HMQC spectrum acquired with fast MAS revealed three distinct 14 N environments. Plane-wave DFT calculations were used to construct model edge structures and predict the corresponding 11 B, 14 N and 15 N SSNMR spectra. Comparison of the experimental and predicted SSNMR spectra confirms that zigzag and armchair edges with both amine and boron hydroxide/oxide termination are present. The detailed characterization of h-BNNS molecular edge termination will prove useful for many material science applications. The techniques outlined here should also be applicable to understand the molecular edge terminations in other 2D materials

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

Determining the Conformation of Supported Complexes Using an 17 O TEDOR-like NMR Experiment

Dynamic nuclear polarization surface-enhanced nuclear magnetic resonance (NMR) spectroscopy has enabled the determination of the three-dimensional configuration of surface sites, in particular supported metal complexes of relevance to single-site heterogeneous catalysis. These approaches have chiefly leveraged the application of NMR double-resonance experiments that either reveal the complex conformation via point-to-point intramolecular distances between spin-labeled atoms or the complex-surface orientation via distances between the spins and the surface plane. Either method typically requires expensive isotope labeling and each reports on different structural features. The application of an experiment that simultaneously reveals both types of distances with chemical resolution would be ideal. Here, in this article, we describe an 17 O{ 1 H} pseudo-3D correlation experiment that achieves this goal. Specifically, Si–O–Si and Si–O–M oxygens are well-resolved by 17 O NMR; therefore, distances can be simultaneously measured radially, between Si– 17 O–M and the 1 H’s of the ligands, and vertically to the Si– 17 O–Si linkages of the silica support. We demonstrate the experiment using supported yttrium and zirconium complexes. Good agreement is obtained when comparing the experimental results to theoretical predictions from density functional theory calculations, highlighting the reliability of this relatively simple experiment.

alkyls↗
DOE OSTI2025

Magnetic semiconductors from functionalized Cr 2 ⁢C MXenes

Here, we report an ab initio investigation of functionalized and 3⁢𝑑-electron-doped Cr 2 ⁢C MXenes. Upon functionalization, the Cr 2⁢ C becomes chemically, dynamically, and mechanically stable, and it exhibits magnetic semiconducting behavior. Cr 2 ⁢CF 2 stands out as a wide band gap semiconductor, possessing superexchange interaction mediated by F atoms within the layer; however, the applied strain transforms it from an indirect to a direct band gap semiconductor. Strong spin-phonon coupling found in Cr 2 ⁢CH 2 is supported by the distorted Cr spin density due to the hydrogen environment. Two magnon branches, associated with two sublattice spins, are found in the ferromagnetic Cr 2 ⁢CO 2 and antiferromagnetic Cr 2 ⁢CF 2 . Depending on the types of 3⁢𝑑-electron dopants and functionalization, Cr 2 ⁢C MXenes (except for Cr 2 ⁢CO 2 ) change from an indirect band gap magnetic semiconductor to different states of electronic and magnetic matter, including an exotic direct band gap magnetic bipolar semiconductor. In addition, we reveal a band inversion between the two highest valence bands in the Fe-doped Cr 2 ⁢CCl 2 .

2D materials↗
DOE OSTI2026

Schwinger-Keldysh effective action for hydrodynamics with approximate symmetries

We study hydrodynamic theories with approximate symmetries in the recently developed effective action approach on the Schwinger-Keldysh contour. We employ the method of spurious symmetry transformation for small explicit symmetry-breaking parameters to systematically constrain symmetry-breaking effects in the nonequilibrium effective action for hydrodynamics. We apply our method to the hydrodynamic theory of chiral symmetry in quantum chromodynamics at finite temperature and density and its explicit breaking by quark masses. We show that the spurious symmetry and the Kubo-Martin-Schwinger relation dictate that the Ward-Takahashi identity for the axial symmetry, i.e., the partial conservation of axial vector current (PCAC) relation, contains a relaxational term proportional to the axial chemical potential, whose kinetic coefficient is at least of the second order in the quark mass. In the phase where the chiral symmetry is spontaneously broken, and the pseudo-Nambu-Goldstone pions appear as hydrodynamic variables, this relaxation effect is subleading compared to the conventional pion mass term in the PCAC relation, which is of the first order in the quark mass. On the other hand, in the chiral symmetry restored phase, we show that our relaxation term, which is of the second order in the quark mass, becomes the leading contribution to the axial charge relaxation. Therefore, the leading axial charge relaxation mechanism is parametrically different in the quark mass across a chiral phase transition.

Goldstone bosons↗
DOE OSTI2026

Detecting the full photoemission cone from laser-based ARPES experiments by leveraging deflector technology

Angle-resolved photoemission spectroscopy (ARPES) provides a direct access to the electronic band structure of solid and molecular systems. The momentum range accessible by this technique depends directly on the photon energy used, and low-photon-energy sources are insufficient to photoemit electrons over the full Brillouin zone of most quantum materials. In addition, while electrons are emitted over a 2π solid angle, conventional hemispherical analyzers only collect a small subset of those electrons. A previous work [Gauthier et al., Rev. Sci. Instrum. 92, 123907 (2021)] demonstrated that electrons emitted over a larger field-of-view can be acquired in one fixed configuration by accelerating them toward the analyzer with a bias voltage. Here, in this study, we extend this work by leveraging the deflector technology of novel ARPES hemispherical analyzers. We demonstrate the ability to detect all 2π photoemitted electrons in a fixed configuration for various materials, such as gold, cuprates, and transition-metal dichalcogenides. This approach is especially advantageous for time-resolved ARPES, as electron dynamics over a large momentum range can be accessed with identical measurement conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Cavity molecular dynamics simulations of vibrational polariton-enhanced molecular nonlinear absorption

We report recent experiments have observed that the chemical and photophysical properties of molecules can be modified inside an optical Fabry– Pérot microcavity under collective vibrational strong coupling (VSC) conditions, and such modification is currently not well understood by theory. In an effort to understand the origin of such cavity-induced phenomena, some recent studies have focused on the effect of the cavity environment on the nonlinear optical response of the molecular subsystem. Here, we use a recently proposed protocol for classical cavity molecular dynamics simulations to numerically investigate the linear and the nonlinear response of liquid carbon dioxide under such VSC conditions following an optical pulse excitation. We find that applying a strong pulse of excitation to the lower hybrid light–matter state, i.e., the lower polariton (LP), can lead to an overall molecular nonlinear absorption that is enhanced by up to two orders of magnitude relative to the excitation outside the cavity. This polariton-enhanced multiphoton absorption also causes an ultrashort LP lifetime (0.2 ps) under strong illumination. Unlike usual polariton relaxation processes—whereby polaritonic energy transfers directly to the manifold of singly excited vibrational dark states—under the present mechanism, the LP transfers energy directly to the manifold of higher vibrationally excited dark states; these highly excited dark states subsequently relax to the manifold of singly excited states with a lifetime of tens of ps. Because the present mechanism is generic in nature, we expect these numerical predictions to be experimentally observed in different molecular systems and in cavities with different volumes.

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

Path integral Monte Carlo approach to the structural properties and collective excitations of liquid $$^3{\text {He}}$$ without fixed nodes

Abstract Due to its nature as a strongly correlated quantum liquid, ultracold helium is characterized by the nontrivial interplay of different physical effects. Bosonic $$^4{\text {He}}$$ 4 He exhibits superfluidity and Bose-Einstein condensation. Its physical properties have been accurately determined on the basis of ab initio path integral Monte Carlo (PIMC) simulations. In contrast, the corresponding theoretical description of fermionic $$^3{\text {He}}$$ 3 He is severely hampered by the notorious fermion sign problem, and previous PIMC results have been derived by introducing the uncontrolled fixed-node approximation. In this work, we present extensive new PIMC simulations of normal liquid $$^3{\text {He}}$$ 3 He without any nodal constraints. This allows us to to unambiguously quantify the impact of Fermi statistics and to study the effects of temperature on different physical properties like the static structure factor $$S({\mathbf {q}})$$ S ( q ) , the momentum distribution $$n({\mathbf {q}})$$ n ( q ) , and the static density response function $$\chi ({\mathbf {q}})$$ χ ( q ) . In addition, the dynamic structure factor $$S({\mathbf {q}},\omega )$$ S ( q , ω ) is rigorously reconstructed from imaginary-time PIMC data. From simulations of $$^3{\text {He}}$$ 3 He , we derived the familiar phonon–maxon–roton dispersion function that is well-known for $$^4{\text {He}}$$ 4 He and has been reported previously for two-dimensional $$^3{\text {He}}$$ 3 He films (Nature 483:576–579 (2012)). The comparison of our new results for both $$S({\mathbf {q}})$$ S ( q ) and $$S({\mathbf {q}},\omega )$$ S ( q , ω ) with neutron scattering measurements reveals an excellent agreement between theory and experiment.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Edge States Drive Exciton Dissociation in Ruddlesden–Popper Lead Halide Perovskite Thin Films

Efficient photovoltaic cells based on thin films of solution-processed 2D Ruddlesden–Popper hybrid perovskites (RPPs) represent an exciting breakthrough due to their enhanced tunability and chemical stability relative to those fabricated from 3D phases. However, reports of efficient charge separation and current collection are in apparent contradiction with the well-known enhancement of the exciton binding energy in multilayered halide perovskites, which should lower the device’s internal quantum efficiency and voltage. This controversy has led to various proposals for the electronic and physical structure of RPP thin films, including phase inhomogeneity as the driving force for exciton dissociation and transport. We address this apparent paradox in high-quality hot-cast RPP films by correlating ultrafast transient absorption spectroscopy with X-ray scattering measurements. We show that a hot-casting fabrication method produces highly phase pure n = 3 (BA) 2 (MA) n-1 Pb n I 3n+1 RPP structures. The high-phase purity and large grain sizes allow us to observe vertical transport of excitons via a diffusive process and allow us to determine that charge separation is primarily driven by dissociation at surface localized subgap electronic states. Furthermore, we analyze the differential absorption kinetics in films of varying thickness to directly determine that the excitonic diffusion constant is ~0.18 cm 2 s –1 . We propose that a surface localized structural distortion, observed using surface selective grazing incidence X-ray scattering measurements, is responsible for the creation of the surface localized defect states. We find that the density and spatial distribution of these defect states is a function of preparation conditions.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Supercomputer-Based Ensemble Docking Drug Discovery Pipeline with Application to Covid-19

In this work, we present a supercomputer-driven pipeline for in silico drug discovery using enhanced sampling molecular dynamics (MD) and ensemble docking. Ensemble docking makes use of MD results by docking compound databases into representative protein binding-site conformations, thus taking into account the dynamic properties of the binding sites. We also describe preliminary results obtained for 24 systems involving eight proteins of the proteome of SARS-CoV-2. The MD involves temperature replica exchange enhanced sampling, making use of massively parallel supercomputing to quickly sample the configurational space of protein drug targets. Using the Summit supercomputer at the Oak Ridge National Laboratory, more than 1 ms of enhanced sampling MD can be generated per day. We have ensemble docked repurposing databases to 10 configurations of each of the 24 SARS-CoV-2 systems using AutoDock Vina. Comparison to experiment demonstrates remarkably high hit rates for the top scoring tranches of compounds identified by our ensemble approach. We also demonstrate that, using Autodock-GPU on Summit, it is possible to perform exhaustive docking of one billion compounds in under 24 h. Finally, we discuss preliminary results and planned improvements to the pipeline, including the use of quantum mechanical (QM), machine learning, and artificial intelligence (AI) methods to cluster MD trajectories and rescore docking poses.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

From bulk to surface: Structure and dynamics of amorphous alumina from deep potential molecular dynamics

Understanding the atomic-scale structure and dynamics of amorphous oxide surfaces is essential for interpreting their chemical reactivity, mechanical stability, and interfacial behavior, yet direct experimental characterization remains challenging. We employ Deep Potential (DP) molecular dynamics to generate large-scale, ab initio -quality models of amorphous Al 2 O 3 bulk glasses and melt-quenched free surfaces, enabling a quantitative analysis of both structure and relaxation dynamics with statistical confidence inaccessible to direct ab initio simulation. The trained DP model reproduces experimental liquid and glass structure, captures the cooling-rate dependence of the bulk glass transition, and corrects systematic biases in the polyhedral populations predicted by widely used classical force fields. At the free surface, mass density recovers to bulk values over ~10 Å, while local coordination requires a slightly wider subsurface region to fully converge. The outermost layer is oxygen-enriched, exhibits altered polyhedral connectivity with contracted Al–O bonds, and hosts a broad population of under-coordinated motifs (notably AlO 3 and OAl 2 ) whose abundances are governed by glass stability. These under-coordinated surface motifs exhibit distinct vibrational signatures and occur as locally paired Lewis acid and Brønsted base sites consistent with bond-valence compensation, yet remain spatially dispersed rather than aggregating into extended clusters. Despite this pronounced structural heterogeneity, surface relaxation and the glass-transition temperature remain comparable to their bulk counterparts, suggesting that the disordered surface is kinetically stable once formed. Together, these results establish a molecular-level picture of amorphous alumina surfaces and demonstrate the capability of machine-learned potentials to resolve structure–property relationships in disordered oxide interfaces.

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

Anionic Lipids Regulate the Light-Harvesting Complex 1-Reaction Center Photocycle in Purple Bacteria

Photosynthetic purple bacteria can capture and convert sunlight with a remarkable, nearly 100% quantum efficiency. The light-harvesting complex 1-reaction center (LH1-RC) core complex is the membrane complex fundamentally responsible for solar energy conversion. LH1-RC has a highly conserved surrounding lipid composition known to favor anionic lipids for an unknown function. Here, in this work, we compared experimentally the rate of LH1-to-RC energy transfer in detergent, membrane nanodiscs with varying lipid compositions, purified membrane fragments, and live cells. The energy transfer rate indicated that RC turnover decreased in neutral lipids, yet was partially restored in anionic lipids, revealing an unexpected lipid dependence. In complementary molecular dynamics simulations, the anionic lipid cardiolipin showed electrostatic interactions with LH1-RC that may mediate quinone exchange, providing a mechanism for the observed lipid dependence. Overall, these results revealed that anionic lipids facilitate LH1-RC redox cycling, identifying a functional role for membrane composition in photosynthetic solar energy conversion.

bacteria↗