Search NASA⌕ Search

SEARCH · Search NASA

Results for “Wave function”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Majorana Edge Modes in Isolated Wires

Topological superconductors are believed to host exotic quasiparticle excitations known as Majorana zero modes (MZMs), with much of the evidence based on BCS mean-field theory. The direct application of mean-field arguments is tenuous in finite, isolated systems relevant in some experiments. Here, we develop a new correlation-based method for identifying MZMs in interacting, number-conserving systems. Using the density matrix renormalization group, we study fermion number-conserving models with long-range interactions, which under periodic boundary conditions exhibit robust topological and nontopological superconductivity, tuned by the strength of interaction [Ortiz et al., Phys. Rev. Lett. 113, 267002 (2014)]. We find evidence that, on the topological side, Majorana edge modes appear in open chains, manifesting as the vanishing of the energy splitting between odd- and even-parity ground states with increasing system size. Additionally, off-diagonal two-point correlation functions show nonlocal, parity-dependent edge effects. These correlations reveal the spatial structure of Majorana modes in the many-body wave function. We show that the correlation diagnostic applies broadly, including to short-range interacting models, where topological superconductivity is more fragile due to the absence of a bulk excitation gap.

Thomas-Markarian, Jaden↗
DOE OSTI2025

Datasets for Custom-trained Machine-learning Interatomic Potentials: Nitric Acid Aqueous Solution

This dataset was generated using an iterative active learning strategy with the ArcaNN software package (https://github.com/arcann-chem/arcann_training) to train machine-learning interatomic potentials (MLIPs) for aqueous nitric acid. Each active-learning cycle consisted of three stages: (1) training, (2) exploration, and (3) labeling. The initial training set comprised approximately 800 randomly selected configurations from a previous study by Lewis et al. (https://doi.org/10.1021/jp205510q), which investigated nitric acid solutions at 2, 3, 4, and 5 mol/L. For all configurations, single-point calculations of atomic forces and total energies were performed at the quantum density functional theory BLYP-D2 and PBE-D3 levels of theory using the CP2K Quickstep module. Valence electrons were treated explicitly, while core electrons on all atoms were represented by norm-conserving Goedecker–Teter–Hutter (GTH) pseudopotentials. Long-range dispersion interactions were accounted for using Grimme dispersion corrections. Wave functions were expanded in a mixed Gaussian-and-plane-wave scheme using TZV2P-MOLOPT basis sets for all elements and an 800 Ry auxiliary plane-wave cutoff for the electron density. Self-consistent field convergence was accelerated using orbital transformation and Direct Inversion in the Iterative Subspace, with a convergence threshold of 10^{-6}. All single-point calculations were carried out in periodic orthorhombic cells whose dimensions match those of the molecular configurations sampled from earlier trajectories. The CELL_REF keyword in CP2K was used to define a fixed reference cell, ensuring consistency in the reference data used for MLIP training, particularly when cell fluctuations are present in NpT simulations. The resulting high-fidelity energies and forces constitute the ground-truth labels used to train the MLIPs contained in this dataset.

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗
DOE OSTI2024

Tetraquarks made of sufficiently unequal-mass heavy quarks are bound in QCD

Tetraquarks, bound states composed of two quarks and two antiquarks, have been the subject of intense study but are challenging to understand from first principles. We apply variational and Green’s function Monte Carlo methods to compute tetraquark ground-state energies in potential nonrelativistic QCD using a wide range of color and spatial wave functions. We find no evidence for bound tetraquarks composed of equal-mass quarks and antiquarks. Conversely, we find clear evidence for the existence of bound tetraquarks for sufficiently unequal quark/antiquark mass ratios at all overall mass scales where our effective theory results are applicable. We predict the critical mass ratios for bound state formation and study tetraquark bound states’ spatial and color structure at leading order and next-to-leading order in potential nonrelativistic QCD. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Custom-trained Machine-learning Interatomic Potentials: ZnCl2 Aqueous Solution

This dataset was generated using an iterative active-learning strategy implemented in the ArcaNN software package (https://github.com/arcann-chem/arcann_training) to train machine-learning interatomic potentials for aqueous ZnCl2 solutions. Each active-learning cycle consisted of three stages: training, exploration, and labeling. The initial training set combined configurations generated in this work from enhanced-sampling ab initio molecular dynamics simulations with configurations from a previously reported neural-network-potential study of aqueous ZnCl2. The enhanced-sampling ab initio molecular dynamics simulations involved Zn–Cl separation and the chloride coordination number around Zn²? as collective variables. These configurations served as the seed dataset. Subsequent active-learning cycles expanded the training set by identifying and labeling configurations that were poorly represented by the current models, thereby improving coverage of ion-association states and changes in local coordination and charge-state environments relevant to the solution free-energy landscape. For all selected configurations, single-point calculations of the total energies and atomic forces were performed within density functional theory using the CP2K Quickstep module. Reference calculations employed the revPBE-D3 and r2SCAN exchange-correlation functionals. Motivated by recent work on aqueous Zn²?, the main revPBE calculations omitted D3 dispersion contributions involving Zn²?, while retaining the D3 correction for water and chloride. For comparison, fully dispersion-corrected revPBE-D3 reference calculations were also performed, with D3 applied to all species, including Zn²?. Valence electrons were treated explicitly, while core electrons were represented using norm-conserving Goedecker–Teter–Hutter pseudopotentials. The wave functions were expanded using the mixed Gaussian-and-plane-wave scheme with TZV2P-MOLOPT basis sets for all elements and a 600 Ry auxiliary plane-wave cutoff for the electron density. Self-consistent-field convergence was accelerated using the orbital-transformation and Direct Inversion in the Iterative Subspace algorithms, with a convergence threshold of 10?6. All single-point calculations were performed in periodic orthorhombic cells. The CELL_REF keyword in CP2K was used to define a fixed reference cell with a box length of 25 Å. This treatment ensured a consistent reference for configurations extracted from NpT trajectories with fluctuating cell dimensions. The resulting DFT energies and atomic forces constitute the ground-truth labels used to train the MLIPs. The resulting MLIP was trained for aqueous ZnCl2 solutions spanning concentrations from 0 to 30 molal and a broad pH range, from strongly acidic to strongly basic conditions. Representative examples of configurations included in the MLIP training dataset are provided below. These include 1) Representative configurations from the dataset labeled at the revPBE-D3 level, with D3 dispersion interactions involving Zn2+ excluded (revPBE-wo-D3). 2) Representative configurations from the dataset labeled at the fully dispersion-corrected revPBE-D3 level, with D3 interactions applied to all species, including Zn2+ (revPBE-D3). 3) Representative configurations from the dataset labeled at the r2SCAN level of theory (r2SCAN).

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗
DOE OSTI2025

Real-Space Quantification of Exciton Localization in Acene Crystals Using Wannier Function Decomposition

Here, we introduce the Wannier function decomposition of excitons (WFDX) method to quantify exciton localization in solids within the ab initio Bethe-Salpeter equation framework. By decomposing each Bloch exciton wave function into products of single-particle electron and hole maximally localized Wannier functions, this real-space approach provides well-defined orbital- and spatial-resolved measures of both Frenkel and charge-transfer excitons at low computational cost. We apply WFDX to excitons in acene crystals, quantifying how the number of rings, the exciton spin state, and the center-of-mass momentum affect spatial localization. Additionally, we show how this real-space representation reflects structural nonsymmorphic symmetries that are hidden in standard reciprocal-space descriptions. We demonstrate how the WFDX framework can be used to efficiently interpolate exciton expansion coefficients in reciprocal-space and outline how it may facilitate evaluation of observables involving position operators, highlighting its potential as a general tool for both analyzing and computing excitonic properties in solids.

Tao, Zui [University of California, Berkeley, CA (↗
DOE OSTI2024

Out-of-time-ordered-correlators for the pure inverted quartic oscillator: classical chaos meets quantum stability

Out-of-time-ordered-correlators (OTOCs) have been suggested as a means to diagnose chaotic behavior in quantum mechanical systems. Recently, it was found that OTOCs display exponential growth for the inverted quantum harmonic oscillator, mirroring the fact that this system is classically and quantum mechanically unstable. In this work, I study OTOCs for the inverted anharmonic (pure quartic) oscillator in quantum mechanics, finding only oscillatory behavior despite the classically unstable nature of the system. For higher temperature, OTOCs seem to exhibit saturation consistent with a value of –2< x 2 > T < p 2 > T at late times. I provide analytic evidence from the spectral zeta-function and the WKB method as well as direct numerical solutions of the Schrödinger equation that the inverted quartic oscillator possesses a real and positive energy eigenspectrum, and normalizable wave-functions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Unraveling the magnetic ground state in the alkali-metal lanthanide oxide Na 2 ⁢Pr⁡O 3

Here, a comprehensive set of muon spin spectroscopy and neutron scattering measurements supported by ab initio and model Hamiltonian simulations have been used to investigate the magnetic ground state of Na 2 ⁢PrO 3 . μSR reveals a Néel antiferromagnetic order below T N ~ 4.9K, with a small static magnetic moment m static ≤ 0.22 μ B /Pr collinearly aligned along the c axis. Inelastic neutron measurements reveal the full spectrum of crystal field excitations and confirm that the Pr 4+ ground-state wave function deviates significantly from the Γ 7 limit that is relevant to the Kitaev model. Single- and two-magnon excitations are observed in the ordered state below T N =4.6K and are well described by nonlinear spin wave theory from the Néel state using a magnetic Hamiltonian with Heisenberg exchange J=1 meV and symmetric anisotropic exchange Γ/J=0.1, corresponding to an XY model. Intense two magnon excitations are accounted for by g-factor anisotropy g z /g ± = 1.29. A fluctuating moment δm 2 = 0.57 ⁢(22) ⁢μ$^2_B$/Pr extracted from the energy and momentum integrated inelastic neutron signal is reduced from expectations for a local J = 1/2 moment with average g factor g avg ≈ 1.1. Together, the results demonstrate that the small moment in Na 2 ⁢PrO 3 arises from crystal field and covalency effects and the material does not exhibit significant quantum fluctuations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Anisotropic excitonic magnetism from discrete C 4 symmetry in CeRhIn 5

Anisotropy in strongly correlated materials is a central parameter in determining the electronic ground state and is tuned through the local crystalline electric field. This is notably the case in the CeCo x Rh 1 − x In 5 system where the ground-state wave function can provide the basis for antiferromagnetism and/or unconventional superconductivity. We develop a methodology to understand the local magnetic anisotropy and experimentally investigate with neutron spectroscopy applied to antiferromagnetic ( T N = 3.8 K ) CeRhIn 5 , which is isostructural to d -wave superconducting ( T c = 2.3 K ) CeCoIn 5 . Through diagonalizing the local crystal field Hamiltonian with discrete tetragonal C 4 point group symmetry and coupling these states with the random phase approximation, we find two distinct modes polarized along the crystallographic c and a − b planes, agreeing with experiment. The anisotropy and bandwidth, underlying the energy scale of these modes, are tuneable with a magnetic field which we use experimentally to separate in energy single and multiparticle excitations thereby demonstrating the instability of excitations polarized within the crystallographic a − b plane in CeRhIn 5 . We compare this approach to a S eff = 1 2 parametrizations and argue for the need to extend conventional SU(2) theories of magnetic excitations to utilize the multilevel nature of the underlying crystal-field basis states constrained by the local point-group C 4 symmetry. Published by the American Physical Society 2024

Brener, D. J. (ORCID:0000000333400177)↗
DOE OSTI2025

Hydrogen Bond Benchmark: Focal‐Point Analysis and Assessment of DFT Functionals

We performed a hierarchical, convergent ab initio benchmark study and systematically analyzed the performance of density functional approximations for describing hydrogen bonds in small neutral, cationic, and anionic complexes, as well as in larger systems involving amide, urea, deltamide, and squaramide moieties. Focal point analyses (FPA), extrapolating to the ab initio limit, were carried out using correlated wave function methods up to CCSDT(Q) for the small complexes and CCSD(T) for the larger systems, together with correlation-consistent Gaussian basis sets up to the complete basis set limit. Optimized geometries and vibrational frequencies were obtained at the CCSD(T) level. The resulting FPA hydrogen-bond energies converge within a few tenths of a kcal mol −1 . These reference data were used to evaluate 60 density functionals (including 12 dispersion-corrected), spanning the local-density approximation (LDA), generalized gradient approximations (GGAs), meta-GGAs, hybrids, meta-hybrids, double-hybrids, and range-separated hybrids. Overall, the meta-hybrid M06-2X provides the best performance for both hydrogen bond energies and geometries, while the dispersion-corrected GGAs BLYP-D3(BJ) and BLYP-D4 also yield accurate hydrogen-bond data and can serve as cost-effective options for studying large and complex systems.

coupled cluster theory↗
DOE OSTI2025

First principles study of the Fermi surface topology of CeCu 2 ⁢Si 2

Since the discovery of heavy-fermion superconductivity in CeCu 2 ⁢Si 2 , the material has attracted great interest, particularly with regard to the nature of the superconducting pairing and its mechanism. Consequently, it is essential to better understand the electronic Fermi surface topology and its role in strong antiferromagnetic fluctuations. The standard density functional theory method is insufficient to model the interplay of strong on-site Coulomb repulsion in localized 4⁢𝑓 electrons and their hybridization with itinerant ligand-orbital electrons. We have performed electronic ground-state calculations on CeCu 2 ⁢Si 2 using the Gutzwiller wave function approximation. The Gutzwiller approximation captures the quasiparticle band renormalization from the strong on-site Coulomb repulsion. We have performed an analysis of this effect on the electronic structure and the Fermi surface topology by varying the interaction strength and taking into account the crystal-field splitting. Using the de Haas-van Alphen effect, the extremal Fermi surface cross-sectional areas were calculated to quantify the effects of quasiparticle mass renormalization on the Fermi surface. Our results confirm the presence of two Fermi surface sheets corresponding to the heavy (488⁢𝑚 𝑒 ) and light (4.35⁢𝑚 𝑒 ) quasiparticles when the crystal-field splitting is accounted for on equal footing with the electronic correlations. This method gives the best agreement with experimental measurements as well as the renormalized band method.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Nonadiabatic dynamics of photoexcited thiopyridone isomers: An interplay between El-Sayed’s conditions and energy gap law

We investigated the non-adiabatic dynamics of photoexcited thiopyridone systems across their ortho-, meta-, and para-isomeric forms. The relaxation pathways of the three isomers in both gas phase and solvent environments are mapped using surface hopping dynamics based on time-dependent density functional theory. Our analysis highlights the influence of isomeric structures on photophysical behavior, offering insights into design principles to control photochemical phenomena. The simulations suggest a systematic reduction in the rate of intersystem crossing (ISC) from ortho- to meta- to para-isomer. Comparisons with multiconfigurational wave function methods in the gas phase further demonstrate how electronic structure influences the predicted dynamical pathways. The simulated dynamics demonstrates that the spin–orbit coupling strength alone does not determine the rate of ISC, as both state energetics and underlying electronic and structural features play decisive roles. These aspects explain the much slower ISC in the para-isomer, as well as the non-negligible role of the El-Sayed forbidden pathway in the computed ISC dynamics.

Complete-active space self-consistent field↗
DOE OSTI2025

Diverse electronic topography in a distorted kagome metal LaTi 3 ⁢Bi 4

Recent reports on a family of kagome metals of the form 𝐿⁢𝑛⁢Ti 3 ⁢Bi 4 (𝐿⁢𝑛 = Lanthanide) have stoked interest due to the combination of highly anisotropic magnetism and a rich electronic structure. The electronic structure near the Fermi level is proposed to exhibit Dirac points and van Hove singularities (VHSs). Here, in this manuscript, we use angle-resolved photoemission spectroscopy measurements in combination with density functional theory calculations to investigate the electronic structure of an interesting kagome metal LaTi 3 ⁢Bi 4 . Our results reveal multiple VHSs with one VHS located in the vicinity of the Fermi level. We clearly observe two flat bands, which originate from the destructive interference of wave functions within the Ti kagome motif. These flat bands and VHSs originate from Ti 𝑑 orbitals and are very responsive to the polarization of the incident beam. We notice a significant anisotropy in the electronic structure, resulting from the breaking of sixfold rotational symmetry in this material. Our findings demonstrate this member of Ti based kagome materials as a promising platform to explore novel emerging phenomena in the wider 𝐿⁢𝑛⁢Ti 3⁢ Bi 4 (𝐿⁢𝑛 = lanthanide) family.

electronic structure↗
DOE OSTI2025

Crystal field splittings and magnetic ground state of the square-lattice antiferromagnets YbBi 2 ⁢ClO 4 and YbBi 2 ⁢IO 4 with J eff = $\frac{1}{2}$

Here, we report on the crystal field level splitting and magnetic ground state of the J eff = $\frac{1}{2}$ square lattice antiferromagnets YbBi 2 ⁢ClO 4 and YbBi 2 ⁢IO 4 using powder inelastic neutron scattering (INS) and neutron diffraction measurements. Both compounds exhibit a well-isolated Γ 7 doublet ground state under a tetragonal crystal field environment, confirming a robust J eff = $\frac{1}{2}$ picture with slight XY-type anisotropic character in the g-tensor. Notably, the ground state wave functions closely resemble the Γ 7 doublet expected in the perfect cubic limit, consistent with the nearly cubic ligand configuration of eight O 2- ions surrounding Yb 3+ . Below T N = 0.21 K, YbBi 2 ⁢IO 4 exhibits a stripe long-range magnetic order characterized by an ordering wave vector q m = (1/2, 0, 0) or its symmetry-equivalent (0, 1/2, 0), with magnetic moments aligned along q m . The ordered moment is approximately 79% of the classical prediction, significantly larger than expected from the isotropic J 1 -J 2 model, suggesting the possible involvement of exchange anisotropy in explaining this observation. We show that symmetry-allowed XXZ and bond-dependent anisotropic exchange terms in a square lattice can play a critical role in stabilizing the stripe order and suppressing the moment reduction as observed. These findings establish YbBi 2 ⁢ClO 4 and YbBi 2 ⁢IO 4 as unique platforms for exploring rich J eff = $\frac{1}{2}$ magnetism from two less investigated perspectives: (i) on a square lattice and (ii) within a (nearly) cubic ligand environment.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Projector‐Based Quantum Embedding Study of Iron Complexes

Projection‐based embedding theory (PBET) is used to calculate and assess the challenging spin‐crossover energies for a selection of small Fe‐containing systems by embedding the metal center into the frozen potential of the ligands. MP2, CCSD, and CCSD(T) are embedded in potentials from the SCAN and r 2 SCAN functionals and compared with the canonical values for the constituent methods and previously reported reference values. Considering the PBET calculations as a correction for the underlying DFT, the embedding calculations are able to provided improvement for most cases. In some cases, the PBET methods are able to compensate for limitations in the wave function methods and produce results similar to more rigorous calculations from the literature. For the systems with spin‐crossover energies near zero, the current methodology fails to provide consistent improvement. In conclusion, the isolated recalculation of the electronic structure around the metal center when embedded into a DFT treatment of the ligand field shows promise as a pragmatic and lower cost treatment compared to the canonical treatment of the whole system of the difficult class of spin‐crossover complexes.

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

Transverse-momentum-dependent pion structures from lattice QCD: Collins-Soper kernel, soft factor, TMDWF, and TMDPDF

We present the first lattice quantum chromodynamics (QCD) calculation of the pion valence-quark transverse-momentum-dependent parton distribution function (TMDPDF) within the framework of large-momentum effective theory (LaMET). Using correlators fixed in the Coulomb gauge (CG), we computed the quasi-TMD beam function for a pion with a mass of 300 MeV, a fine lattice spacing of 𝑎 =0.06 fm, and multiple large momenta up to 3 GeV. The intrinsic soft functions in the CG approach are extracted from form factors with large momentum transfer, and as a byproduct, we also obtain the corresponding Collins-Soper (CS) kernel. Our determinations of both the soft function and the CS kernel agree with perturbation theory at small transverse separations (𝑏 ⊥ ) between the quarks. At larger 𝑏 ⊥ , the CS kernel remains consistent with recent results obtained using both CG and gauge-invariant TMD correlators in the literature. By combining next-to-leading logarithmic factorization of the quasi-TMD beam function and the soft function, we obtain an 𝑥-dependent pion valence-quark TMDPDF for transverse separations 𝑏 ⊥ ≳1 fm. Interestingly, we find that the 𝑏 ⊥ dependence of the phenomenological parametrizations of TMDPDF for moderate values of 𝑥 are in reasonable agreement with our QCD determinations. In addition, we present results for the transverse-momentum-dependent wave function for a heavier pion with 670 MeV mass.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

From Oxo to Oxyl to Biradical: Systematic Multireference Calculations of Methane Activation at MOF Nodes

Methane C–H activation at transition-metal sites often involves electronic structures that challenge conventional single-reference electronic structure descriptions. Although Kohn–Sham density functional theory (DFT) is widely used to study catalytic trends, its reliability for reactions involving strongly correlated species remains uncertain. Here we present a systematic multireference investigation of methane activation at metal–organic framework (MOF) node catalysts across the 3d transition-metal series. We introduce an automated workflow for active space selection to enable consistent application of multireference methods, including multiconfiguration pair-density functional theory and n-electron valence state perturbation theory, to these catalytic systems. These calculations show substantial static correlation in the C–H activation reaction step and predict activation barriers that differ from DFT by 30–70 kJ mol–1, with DFT often qualitatively disagreeing in barrier height trends across transition metals. Analysis of multireference wave functions shows that reactivity is governed by the electronic structure of the M–O moiety along a continuum from metal–oxo to oxyl radical and O biradical character. Increased oxygen-centered spin density and weakened M–O bonding are identified as descriptors of catalytic activity which correlate with lower activation barriers.

Wardzala, Jacob↗
DOE OSTI2025

Connectivity-Dependent Exciton–Phonon Coupling in Cesium Bismuth Halide Quantum Dots

Metal halide octahedra form the fundamental functional building blocks of metal halide perovskites, dictating their structures, optical properties, electronic structures, and dynamics. Here, in this study, we show that the connectivity of bismuth halide octahedra in Cs 3 Bi 2 Br 9 and Cs 3 Bi 2 I 9 quantum dots (QDs) changes with different halide elements. We use first-principles calculations to reveal the key role of the connectivity of bismuth halide octahedra on the wave function symmetry, Huang-Rhys factor, and exciton-phonon interaction strength. Following QD synthesis via a ligand-mediated transport method, the effect of connectivity is verified with transient absorption spectroscopy, where we contrast Cs 3 Bi 2 Br 9 and Cs 3 Bi 2 I 9 QD exciton dynamics. In photoexcited Cs 3 Bi 2 I 9 QDs, phonons related to the vibrational motions of face-sharing [BiI 6 ] 3- bioctahedra couple strongly to the electronic state and drive rapid carrier relaxation. Equivalent signals are not observed for photoexcited Cs 3 Bi 2 Br 9 QDs, implying a lack of phonon involvement in band-edge absorption and subsequent exciton relaxation. Our findings suggest that structural engineering can effectively tune the exciton-phonon coupling and therefore influence exciton relaxation and recombination in perovskite nanomaterials.

TDDFT↗
DOE OSTI2025

Hydrogen and deuterium tunneling in niobium

We use density functional methods to identify the atomic configurations of H and D atoms trapped by O impurities embedded in bulk Nb. The O atoms are located at the octahedral position in the Nb body-centered cubic (BCC) lattice, and H (D) atoms tunnel between two degenerate tetrahedral sites separated by a mirror plane. Using nudged elastic band (NEB) methods, we calculate the double-well potential for O-H and O-D impurities and the wave functions and tunnel splittings for H and D atoms. Our results agree with those obtained from analysis of heat capacity and neutron scattering measurements on Nb with low concentrations of O-H and O-D.

density functional theory↗