Search NASA⌕ Search

SEARCH · Search NASA

Results for “experimental comparison”

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

A robust spectral element implementation of the $k - τ$ RANS model in Nek5000/NekRS

The $k - ω$ Reynolds Averaged Navier Stokes (RANS) model is one of the industry standard approaches for modeling of turbulent flows. It performs better than the $k - ϵ$ model for low Reynolds number flows and is also more suitable for boundary layers with adverse pressure gradients. Major drawback of the model, however, is that the asymptotic value of $ω$ at the walls is singular, necessitating the use of a contrived “sufficiently” large value for $ω$ as the boundary condition for its transport equation. Here, this invariably leads to the solution being sensitive to near wall grid spacing. While an acceptable solution for low order (finite volume) methods, the excessive near wall gradients lead to persistent numerical stability issues in high order codes. To alleviate the problem, specifically in the context of the high order spectral element code Nek5000, a regularized $k - ω$ approach was formulated in our prior work (Tomboulides et al., 2018). The formulation, however, relies on the use of wall distance and its gradients for modeling the closure terms and can pose problems for simulations in complex geometries. This work presents a novel implementation of the $k - τ$ RANS model in Nek5000, where $τ = 1/ω$, eliminating the need for regularization, owing to the asymptotically bounded behavior of the source terms in the $τ$ transport equation, and also eliminating dependence on wall distance. Robustness and stability of the $k - τ$ model is ensured through implicit treatment of the source terms and their careful numerical implementation and demonstrated through several cases aimed at verification and validation. Studies include both canonical and engineering relevant problems, viz., turbulent channel flow, pipe flow, backward facing step, flow over NACA 0012 airfoil and flow in a T-junction. Results from the $k - τ$ model are shown to be consistent with regularized $k - ω$ model and also with the $k - ω$ SST model in OpenFOAM (for select studies). Comparison with experimental data is also shown, where available, to bolster validation efforts for the $k - τ$ model implementation through prediction of key turbulent quantities of interest.

Nek5000↗
DOE OSTI2025

2D CFD of lean premixed hydrogen–air flame quenching under locomotive engine conditions

Hydrogen (H 2 ) is a promising fuel for reducing emissions in heavy-duty internal combustion engines (ICEs), but its low quenching distance increases the risk of flame propagation into narrow crevice regions, such as piston-liner gaps. This work uses detailed CFD simulations with finite-rate chemistry to investigate premixed H 2 –air flame quenching in a two-dimensional (2D) region consistent with the piston-liner gap of a diesel ICE. Model accuracy was assessed by comparison with experimentally measured quenching distances in an annular stepwise diverging tube (ASDT). A parametric study was conducted to assess the influence of crevice width (0.05–1.18 mm), pressure (50–150 bar), unburned gas temperature (431–573 K), and equivalence ratio ( Φ= 0.3–0.6). Results show that the critical Péclet number for flame survival is no greater than 3.25, consistent with prior literature, even under ultra-lean and high-pressure conditions (Φ ≤ 0.3, > 50 bar). Additionally, flames with Péclet numbers exceeding 6.15 exhibited front wrinkling, suggesting the onset of velocity-driven instabilities and enhanced flame robustness. These findings help define thresholds for flame quenching in confined geometries and support the safe design of H 2 fueled ICEs.

08 HYDROGEN↗
DOE OSTI2022

Analytic model of dislocation density evolution in fcc polycrystals accounting for dislocation generation, storage, and dynamic recovery mechanisms

Here, an analytic model of the evolution of dislocation density in fcc polycrystals is described. The evolution equations approximately account for most known dislocation storage, dynamic recovery, and dislocation generation mechanisms in fcc polycrystals. Specifically, the model incorporates network (forest) and grain boundary storage, mobile-network and mobile–mobile annihilation, screw–screw annihilation via athermal and thermal single cross-slip, generation by double cross-slip (Koehler mechanism, including dipole formation), Frank-Read sources, grain boundary nucleation, and mobile–immobile dislocation nucleation due to shock loading. Single cross-slip is assumed to proceed through the Friedel–Escaig (FE) mechanism; the corresponding activation energy is calculated using a modified FE model. The activation energy for double cross-slip is calculated for the first time by extending the FE model. The exact evolution equations are integro-differential equations, and as such are difficult to implement in a code; hence, the evolution equations are simplified by making several approximations. Preliminary results on copper are presented, including comparisons to experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Crystal plasticity modeling of strain-induced martensitic transformations to predict strain rate and temperature sensitive behavior of 304 L steels: Applications to tension, compression, torsion, and impact

This paper advances crystallographically-based Olson-Cohen (direct γ → α’) and deformation mechanism (indirect γ→ε→α’) phase transformation models for predicting strain-induced austenite to martensite transformation. Here, the advanced transformation models enable predictions of not only strain-path sensitive, but also of strain-rate and temperature sensitive deformation of polycrystalline stainless steels (SSs). The deformation of constituent grains in SSs is modeled as a combination of anisotropic elasticity, crystallographic slip, and phase transformation, while the hardening is based on the evolution of dislocation density and explicit shifts in phase fractions. Such grain-scale deformation is implemented within the meso-scale elasto-plastic self-consistent (EPSC) homogenization model, which is coupled with the implicit finite element (FE) method to provide a constitutive response at each FE integration point for solving boundary value problems at the macro-scale. Parameters pertaining to the hardening and transformation models within FEEPSC are calibrated and validated on a suite of data including flow curves and phase fractions for monotonic compression, tension, and torsion as a function of strain-rate and temperature for wrought and additively manufactured (AM) SS304L. To illustrate the potential and accuracy of the integrated multi-level FE-EPSC simulation framework, geometry, mechanical response, phase fractions, and texture evolution are simulated during gas-gun impact deformation of a cylinder and quasi-static tension of a notched specimen made of AM SS304L. Details of the simulation framework, comparison between experimental and simulation results, and insights from the results are presented and discussed.

304L steels↗
DOE OSTI2025

An extension of the localized artificial diffusivity method for immiscible and high density ratio flows

The localized artificial diffusivity (LAD) method is widely regarded as the preferred multi-material regularization scheme for the compact finite difference method, because it is conservative, easy to implement, and generally robust for a wide range of multi-material problems. However, traditional LAD methods face significant challenges when applied to flows with large density ratios and when maintaining thermodynamic equilibrium across material interfaces. These limitations arise from the formulation of the artificial diffusivity flux and the reliance on enthalpy diffusion for interface regularization. Additionally, traditional LAD methods struggle to ensure stability under large density ratio conditions, fail to maintain a finite interface thickness, and are therefore unsuitable for modeling immiscible interfaces. Here, in this work, we discuss the origins of these issues in traditional LAD methods and propose modifications which enable the simulation of large density ratio and immiscible flows. The proposed method targets the artificial diffusion fluxes at gradients and ringing in the volume fraction, rather than the mass fraction in traditional methods, to consistently regularize large density ratio interfaces. Furthermore, the proposed method introduces an artificial bulk density diffusion term to enforce equilibrium conditions across interfaces. To address the challenge of modeling immiscible flows, a conservative diffuse interface term is incorporated into the formulation to ensure a finite interface thickness. Specific consideration is taken in the design of the method to ensure that these crucial properties are maintained for N -material flows. The effectiveness of the proposed method is demonstrated through a series of canonical test cases, and its accuracy is validated by comparison with experimental data on micro-bubble collapse in water. These results highlight the method’s robustness and its ability to overcome the limitations of traditional LAD approaches.

Artificial diffusivity↗
DOE OSTI2023

Failure analysis of nuclear transient-tested UN tristructural isotropic fuel particles in a 3D printed SiC matrix

Fully ceramic microencapsulated fuel elements containing UN tristructural isotropic (TRISO) fuel particles within a 3D printed SiC matrix were subjected to transient testing with varying energy depositions. Detailed post-irradiation examinations were performed, including leaching in hot HNO3 and post-leaching X-ray computed tomography, to quantify the percentage of failed TRISO particles and crack propagation within the particles and surrounding fuel matrix. In parallel, detailed finite element analyses were performed for comparison with experimental findings and to better evaluate transient failure modes. The lowest transient energy deposition—which still exceeded bounding values for high-temperature gas-cooled reactor applications—resulted in no detectable TRISO particle failures or matrix cracking, which was consistent with the simulations. Simulations of the higher-energy transients for which significant TRISO particle failure was expected were generally able to reproduce the transient temperatures and matrix cracking. Thus, the TRISO particle failures were explained based on the effects of local SiC matrix thickness and porosity. Results generally confirmed the high strength of the additively manufactured SiC matrix but also affirmed the need for a modified UN TRISO architecture to prevent SiC matrix cracks from propagating through TRISO layers. This unique failure mode has not historically been considered for TRISO fuels contained in weaker graphite matrices.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Rovibronic molecular line list for the $N_2(C^3Π_u–B^3Π_g)$ second positive system

Here, a line list for the N second positive system, $B^3Π_g—C^3Π_u$, has been compiled using the PGOPHER spectral simulation software. The line list extends the number of vibrational states of the $B^3Π_g$ up to v=29 and a maximum rotational state of J=150 for simulation temperatures up to 7000 K. New electronic–vibrational transition moments were calculated using refined potential energy curves and a transition dipole moment with the DUO software. Comparisons to experimental data and the SPECAIR software have been used to validate the new line list. The results are available in ASCII ExoMol .state and .trans files and as a PGOPHER input file for use in spectral analysis.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Ab initio Molecular Dynamics Assessment of Thermodynamic and Transport Properties in (K,Li)Cl and (K, Na)Cl Molten Salt Mixtures

Molten salt mixtures are integral part of highly important technological applications such as nuclear reactors. However, due to inherent difficulties associated with experiment at high temperatures and the intrinsic complexity of liquid-phase multi-component systems, understanding their properties at a molecular level remains a challenge. Here, we report on an ab initio molecular dynamics investigation on structural, electronic, transport, and thermal properties of two common molten salt mixtures, (K, Li)Cl and (K,Na)Cl, at five different compositions and three temperatures. Most of the properties were found to depend on both composition and temperature. While properties, like atomic charges, show additive behaviors, other properties, such as electrical conductivity, show considerable deviations from additivity. We shall show that the mixing of the molten salt mixtures is mainly driven by entropy, and that the KCl and LiCl mix better than KCl and NaCl. Our computational results are in general consistent with available experimental data. Comparison with available theoretical data is also provided.

Nguyen, Manh Thuong↗
DOE OSTI2020

Longitudinal dependence of B and D meson nuclear modifications in heavy-ion collisions at RHIC and the LHC

It is widely acknowledged that heavy flavor probes are sensitive to the properties of the quark-gluon plasma and are often considered an important tool for the plasma tomography studies. Forward rapidity observables can provide further insight on the dynamics of the medium due to the interplay between the medium size and the differences in the production spectra of heavy quark probes. In this proceedings we present the nuclear modification factor R AA 's for B and D mesons, as well as heavy flavor leptons, in the rapidity range -4.0 < y < 4.0 obtained from relativistic Langevin equation with gluon radiation coupled with a (3+1)-dimensional viscous hydrodynamics medium background. Finally, we present comparison with experimental data at mid-rapidity as well as predictions for different rapidity ranges.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Laser Doppler vibrometry for piezoelectric coefficient ($d_{33}$) measurements in irradiated aluminum nitride

Sensors used for experiments in advanced reactors must survive in harsh environments. Few material systems can be used to construct sensors viable for extreme conditions. Aluminum nitride (AlN) is one such material because it has high thermal stability and radiation resistance and sustains good piezoelectric and dielectric properties at high temperatures. In this work, the piezoelectric coefficient $d_{33}$ of the AlN single-crystal with thermal and irradiation damage was investigated using an indirect method with a laser Doppler vibrometer (LDV). Surface electrodes were deposited on the AlN samples, and the vibration response of the samples to an applied voltage was monitored using the LDV as a function of the excitation frequency. The $d_{33}$ estimation was based on the excitation voltage and the thickness-mode displacement extracted from the LDV measurements. Further, six AlN substrates were irradiated with 8 MeV Al 2+ at three fluences (10 15 , 10 16 , and 10 17 ions/cm 2 ) and two temperatures (300 °C and 500 °C). The $d_{33}$ for the six irradiated samples and one pristine sample were measured, and the measurement uncertainty was estimated based on five repeated tests. All samples were also measured by a commercial piezometer for comparison. The experimental results demonstrate that the piezoelectric coefficients obtained by LDV were about 0.8–1.16 pm/V lower than those obtained by the piezometer. With the compensation of the clamping effect, the corrected LDV values are similar to the piezometer results. Both show similar trends in all samples, which validates the feasibility of the proposed method for $d_{33}$ measurement. Based on the LDV results, the irradiated samples show a 12%–22% decrease in $d_{33}$ compared with the pristine samples. The samples irradiated under the same fluence at a higher temperature (500 °C) demonstrated a lower $d_{33}$ than those at 300 °C. The effect of the retro-reflective tape, sample temperature, and sample size on the $d_{33}$ measurement were also studied.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Bimetallic-Derived Catalytic Structures for CO 2 -Assisted Ethane Activation

In recent years, the simultaneous upgrading of CO 2 and ethane has emerged as a promising approach for generating valuable gaseous (CO, H 2 , and ethylene) and liquid (aromatics and C 3 oxygenates) chemicals from the greenhouse gas CO 2 and large-reserved shale gas. The key challenges for controlling product selectivity lie in the selective C-H and C-C bond cleavage of ethane with the assistance of CO 2 . Bimetallic-derived catalysts likely undergo alloying or oxygen-induced segregation under reaction conditions, thus providing diverse types of interfacial sites, e.g., metal/support (M/M'O x ) interface and metal oxide/metal (M'O x /M) inverse interface, that are beneficial for selective CO 2 -assisted ethane upgrading. The alloying extent can be initially predicted by cohesive energy and atomic radius (or Wigner-Seitz radius), while the preference for segregation to form the on-top suboxide can be approximated using the work function, electronegativity, and binding strength of adsorbed oxygen. Furthermore, bimetallic-derived catalysts are typically supported on high surface area oxides. Modifying the reducibility and acidity/basicity of the oxide supports and introducing surface defects facilitate CO 2 activation and oxygen supplies for ethane activation. Using in situ synchrotron characterization and density functional theory (DFT) calculations, we found that the electronic properties of oxygen species influence the selective cleavage of C-H/C-C bonds in ethane, with electron-deficient oxygen over the metal (or alloy) surface promoting nonselective bond scission to produce syngas and electron-enriched oxygen over the metal oxide/metal interface enhancing selective C-H scission to yield ethylene. Here, we further demonstrate that the preferred structures of the catalyst surfaces, either alloy surfaces or metal oxide/metal inverse interfaces, can be controlled through the appropriate choice of metal combinations and their atomic ratios. Through a comprehensive comparison of experimental results and DFT calculations, the selectivity of C-C/C-H bond scission is correlated with the thermodynamically favorable bimetallic-derived structures (i.e., alloy surfaces or metal oxide/metal inverse interfaces) under reaction conditions over a wide range of bimetallic catalysts. These findings not only offer structural and mechanistic insights into bimetallic-derived catalysts but also provide design principles for selective catalysts for CO 2 -assisted activation of ethane and other light alkanes. This Account concludes by discussing challenges and opportunities in designing advanced bimetallic-derived catalysts, incorporating new reaction chemistries for other products, employing precise synthesis strategies for well-defined structures with optimized site densities, and leveraging time/spatial/energy-resolved in situ spectroscopy/scattering/microscopy techniques for comprehensive structural analysis. The research methodologies established here are helpful for the investigation of dynamic alloy and interfacial structures and should inspire more efforts toward the simultaneous upgrading of CO 2 and shale gas.

03 NATURAL GAS↗
DOE OSTI2021

Computational Spectroscopy of the Cr–Cr Bond in Coordination Complexes

In this work, we report the accurate computational vibrational analysis of the Cr–Cr bond in dichromium complexes using second-order multireference complete active space methods (CASPT2), allowing direct comparison with experimental spectroscopic data both to facilitate interpreting the low-energy region of the spectra and to provide insights into the nature of the bonds themselves. Recent technological development by the authors has realized such computation for the first time. Accurate simulation of the vibrational structure of these compounds has been hampered by their notorious multiconfigurational electronic structure that yields bond distances that do not correlate with bond order. Some measured Cr–Cr vibrational stretching modes, ν(Cr 2 ), have suggested weaker bonding, even for so-called ultrashort Cr–Cr bonds, while others are in line with the bond distance. Here, we optimize geometries and compute ν(Cr 2 ) with CASPT2 for three well-characterized complexes, Cr 2 (O 2 CCH 3 ) 4 (H 2 O) 2 , Cr 2 (mhp) 4 , and Cr 2 (dmp) 4 . We obtain CASPT2 harmonic ν(Cr 2 ) modes in good agreement with experiment at 282 cm –1 for Cr 2 (mhp) 4 and 353 cm –1 for Cr 2 (dmp) 4 , compute 50 Cr and 54 Cr isotope shifts, and demonstrate that the use of the so-called IPEA shift leads to improved Cr–Cr distances. Additionally, normal mode sampling was used to estimate anharmonicity along ν(Cr 2 ), leading to an anharmonic mode of 272 cm –1 for Cr 2 (mhp) 4 and 333 cm –1 for Cr 2 (dmp) 4 .

36 MATERIALS SCIENCE↗
DOE OSTI2025

FLAMES─Fast, Low-Storage, Accurate, and Memory-Efficient Adaptive Sampling─Approach to Resolve Spatially Dependent Dynamics of Molecular Liquids

Many critical phenomena in soft matter occur at large length scales, necessitating the resolution of their structure and dynamics at low wavenumbers. However, resolving wavenumber-dependent dynamics computationally via molecular dynamics simulations presents significant challenges, as these phenomena span several orders of magnitude in both time and length scales, resulting in high computational costs and memory demands. Here, this work highlights the computational and memory challenges associated with analyzing molecular trajectories in reciprocal space and demonstrates a method to address them. We introduce FLAMESFast, Low-storage, Accurate, and Memory-Efficient adaptive Sampling, which is a direct method for calculation of structure factors, allowing us to select only the required number of wavevectors for binning. We also use wavenumber-dependent time steps to extract dynamics. Our FLAMES approach effectively mitigates computational and memory/storage bottlenecks. We demonstrate the method using simulations of a model system, liquid octane, at various temperatures. Comparisons with experimental data and real space computation show that the FLAMES technique achieves high accuracy in resolving temperature- and spatially dependent dynamics while being significantly more computationally efficient and requiring less memory and storage than methods based on a uniform wavevector grid and fixed temporal spacing.

Chen, Guang [Argonne National Laboratory (ANL), Ar↗
DOE OSTI2025

Multireference Equation-of-Motion Driven Similarity Renormalization Group: Theoretical Foundations and Applications to Ionized States

We present a formulation and implementation of an equation-of-motion (EOM) extension of the multireference driven similarity renormalization group (MR-DSRG) formalism for ionization potentials (IP-EOM-DSRG). The IP-EOM-DSRG formalism results in a Hermitian generalized eigenvalue problem, delivering accurate ionization potentials for strongly correlated systems. The EOM step scales as O(N 5 ) with the basis set size N, allowing for efficient calculation of spectroscopic properties, such as transition energies and intensities. The IP-EOM-DSRG formalism is combined with three truncation schemes of the parent MR-DSRG theory: an iterative nonperturbative method with up to two-body excitations [MR-LDSRG(2)] and second- and third-order perturbative approximations [DSRG-MRPT2/3]. We benchmark these variants by computing (1) the vertical valence ionization potentials of a series of small molecules at both equilibrium and stretched geometries; (2) the spectroscopic constants of several low-lying electronic states of the OH, CN, N 2 + , and CO + radicals; and (3) the binding curves of low-lying electronic states of the CN radical. A comparison with experimental data and theoretical results shows that all three IP-EOM-DSRG methods accurately reproduce the vertical ionization potentials and spectroscopic constants of these systems. Notably, the DSRG-MRPT3 and MR-LDSRG(2) versions outperform several state-of-the-art multireference methods of comparable or higher cost.

Hamiltonians↗
DOE OSTI2022

Inelastic Neutron Scattering from Thin Film Biaxially Oriented Polyethylene Terephthalate

Recent interest in emerging processes for polymer manufacturing and bio-based chemistries for direct chemical recycling/upcycling has motivated new research focused on a deeper understanding of atomic-scale polymer properties and how they influence macroscopic phenomena. Uncovering the fundamental properties of polymers that give rise to macroscopic behavior could enable new pathways for improved recyclability or utilization of alternative “greener” polymer analogues. In this study, the neutron vibrational spectrum was measured for a film of biaxially oriented polyethylene terephthalate (BoPET) using inelastic neutron scattering (INS), to investigate the relationship between the structure and dynamics of a widely used polymer. Compared to conventional spectroscopic techniques, the use of INS is advantageous for polymeric materials due to the absence of selection rules (i.e., all transitions are allowed), broad-band energy range, and considerable sensitivity to hydrogen modes. In order to distinguish the vibrational modes caused by trans and gauche rotational isomerism, the normal modes of vibration were calculated from a density functional theory-optimized structure of crystalline PET (cPET), representative of the all-trans state, and compared with INS from “highly crystalline” PET powder. Although in- and out-of-plane wagging of hydrogens on the ring structure exhibit significant contribution to both BoPET and cPET spectra, the wagging, rocking, and twisting modes of hydrogen on the ethylene glycol group are, in most cases, conformation-specific. These results were further rationalized by investigating the role of hyperconjugation in stabilizing both conformations using the natural bond order method. Here, through comparison of experimental and calculated INS results, this work provides the fundamental basis for discovering the role of structure and dynamics in shaping the macroscopic properties of PET and polymer analogues.

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

X-ray Spectroscopy Fingerprints of Pristine and Functionalized Graphene

In this work, we demonstrate how to identify and characterize the atomic structure of pristine and functionalized graphene materials from a combination of computational simulation of X-ray spectra, on the one hand, and computer-aided interpretation of experimental spectra, on the other. Despite the enormous scientific and industrial interest, the precise structure of these 2D materials remains under debate. As we show in this study, a wide range of model structures from pristine to heavily oxidized graphene can be studied and understood with the same approach. We move systematically from pristine to highly oxidized and defective computational models, and we compare the simulation results with experimental data. Comparison with experiments is valuable also the other way around; this method allows us to verify that the simulated models are close to the real samples, which in turn makes simulated structures amenable to several computational experiments. Our results provide ab initio semiquantitative information and a new platform for extended insight into the structure and chemical composition of graphene-based materials.

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

Dipolar Heteronuclear Correlation Solid-State NMR Experiments between Half-Integer Quadrupolar Nuclei: The Case of 11 B– 17 O

With 73% of all NMR-active nuclei being quadrupolar, there is great interest in the development of NMR experiments that can probe the proximity of quadrupolar spins. Here, pulse sequences for magic-angle spinning (MAS) 11 B– 17 O resonance-echo saturation-pulse double-resonance (RESPDOR) and dipolar heteronuclear multiple quantum correlation (D-HMQC) solid-state NMR experiments were investigated. In these pulse sequences, rotational-echo double-resonance (REDOR) recoupling was used with central transition (CT)-selective π-pulses applied to either the 11 B or 17 O spins to recouple 11 B– 17 O dipolar interactions. 11 B{ 17 O} RESPDOR experiments on 17 O-enriched boric acid and benzene diboronic acid showed that application of dipolar recoupling on the 11 B channel yielded more signal dephasing than when recoupling is applied on the 17 O channel; however, short effective 11 B transverse relaxation time constants (T 2 ') hinder the acquisition of dephasing curves out to long recoupling durations. Application of REDOR recoupling to 17 O spins was found to produce significant dephasing without compromising the 11 B T 2 '. Comparison of experimental 11 B{ 17 O} RESPDOR curves to those of numerical simulations enabled the 17 O isotopic enrichment to be estimated. 2D 11 B{ 17 O} D-HMQC spectra were recorded with either 11 B or 17 O REDOR recoupling under a variety of radio frequency field conditions. Lastly, 2D 11 B{ 17 O} and 23 Na{ 17 O} D-HMQC spectra of an 17 O-enriched sodium borate glass were acquired to demonstrate the practical application of these heteronuclear correlation experiments to probe structural connectivity between two quadrupolar spins. Importantly, the high-field 2D 11 B– 17 O D-HMQC NMR spectrum revealed two unique 17 O sites correlating to 4-coordinate BO 4 ( [4] B), which were attributed to the [3] B–O– [4] B and [4] B–O– [4] B bridging O atoms. Furthermore, the heteronuclear correlation experiments outlined here should be applicable to a variety of quadrupolar spin pairs.

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

Is Deprotonation of the Oxygen-Evolving Complex of Photosystem II during the S 1 → S 2 Transition Suppressed by Proton Quantum Delocalization?

We address the protonation state of the water-derived ligands in the oxygen-evolving complex (OEC) of photosystem II (PSII), prepared in the S 2 state of the Kok cycle. We perform quantum mechanics/molecular mechanics calculations of isotropic proton hyperfine coupling constants, with direct comparisons to experimental data from two-dimensional hyperfine sublevel correlation (HYSCORE) spectroscopy and extended X-ray absorption fine structure (EXAFS). We find a low-barrier hydrogen bond with significant delocalization of the proton shared by the water-derived ligand, W1, and the aspartic acid residue D1–D61 of the D1 polypeptide. The lowering of the zero-point energy of a shared proton due to quantum delocalization precludes its release to the lumen during the S 1 → S 2 transition. Retention of the proton facilitates the shuttling of a proton during the isomerization of the tetranuclear manganese–calcium–oxo (Mn 4 Ca–oxo) cluster, from the “open” to “closed” conformation, a step suggested to be necessary for oxygen evolution from previous studies. Our findings suggest that quantum-delocalized protons, stabilized by low-barrier hydrogen bonds in model catalytic systems, can facilitate the accumulation of multiple oxidizing equivalents at low overpotentials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗