Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distribution 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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Local structure study on magnetostrictive material Tb 1–x Dy x Fe 2

Tb 1–x Dy x Fe 2 system has attracted more research interest due to the large magnetostrictive effect. The crystal structures and physical properties have been well studied, but research studies on their local structures are still rare. As such, in this work, the local structure of Tb 1–x Dy x Fe 2 samples was studied using the pair distribution function and x-ray absorption spectroscopy techniques. The results demonstrate that the system owns the same local crystal symmetry with its average structure in the ferromagnetic phase, and the crystal lattice of the system is more ordered with increasing Dy content, indicating that the Dy-rich tetragonal phase is more stable than the Tb-rich rhombohedral phase. The different roles of metallic bonds in affecting the crystal lattice are presented. The weak Fe 1 –Fe 2 bonds influenced by the local environment such as local stress from randomly distributed nanodomains could originate the anomalies in the lattice, resulting in the more ordered and stable Dy-rich phase than the Tb-rich phase.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Influences of δB contribution and parallel inertial term of energetic particles on MHD-kinetic hybrid simulations: a case study of the 1/1 internal kink mode

The magnetohydrodynamic-kinetic (MHD-kinetic) hybrid model (Park et al 1992 Phys. Fluids B 4 2033–7) has been widely applied in studying energetic particles (EPs) problems in fusion plasmas for past decades. The pressure-coupling scheme or the current-coupling scheme is adopted in this model. However, two noteworthy issues arise in the model application: firstly, the coupled term introduced in the pressure-coupling scheme, (∇•P h ) ⟂ , is often simplified by ∇•P h , which is equivalent to neglecting the parallel inertial term of EPs; secondly, besides the $δf$ contribution caused by changing in the EP distribution function, the magnetic field perturbation (the $δB$ contribution) generated during development of the instabilities should also be considered, but it is often ignored in existing hybrid simulations. In this paper, we derive the analytical formulations under these two coupling schemes and then numerically study the representative case of the linear stability of the $m/n$ = $1/1$ internal kink mode (IKM) (Fu et al 2006 Phys. Plasmas 13 052517) by using the CLT-K code. Further, it is found that the approximated models can still yield reasonable results when EPs are isotopically distributed. But it fails completely in cases with anisotropic EP distributions. In addition, we further investigate the influence of EP's orbit width on the stability of IKM and verify the equivalence between pressure-coupling scheme and the current-coupling scheme.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Detailed study of quark-hadron duality in spin structure functions of the proton and neutron

Background: The response of hadrons, the bound states of the strong force (QCD), to external probes can be described in two different, complementary frameworks: as direct interactions with their fundamental constituents, quarks and gluons, or alternatively as elastic or inelastic coherent scattering that leaves the hadrons in their ground state or in one of their excited (resonance) states. The former picture emerges most clearly in hard processes with high momentum transfer, where the hadron response can be described by the perturbative expansion of QCD, while at lower energy and momentum transfers, the resonant excitations of the hadrons dominate the cross section. The overlap region between these two pictures, where both yield similar predictions, is referred to as quark-hadron duality and has been extensively studied in reactions involving unpolarized hadrons. Some limited information on this phenomenon also exists for polarized protons, deuterons, and 3 He nuclei, but not yet for neutrons. Purpose: In this paper, we present comprehensive and detailed results on the correspondence between the extrapolated deep inelastic structure function g 1 of both the proton and the neutron with the same quantity measured in the nucleon resonance region. Thanks to the fine binning and high precision of our data, and using a well-controlled perturbative QCD (pQCD) fit for the partonic prediction, we can make quantitative statements about the kinematic range of applicability of both local duality and global duality. Method: We use the most updated QCD global analysis results at high x from the Jefferson Lab Angular Momentum Collaboration to extrapolate the spin structure function g 1 into the nucleon resonance region and then integrate over various intervals in the scaling variable x. We compare the results with the large data set collected in the quark-hadron transition region by the CLAS Collaboration, including, for the first time, deconvoluted neutron data, integrated over the same intervals. We present this comparison as a function of the momentum transfer Q 2 . Results: We find that, depending on the integration interval and the minimum momentum transfer chosen, a clear transition to quark-hadron duality can be observed in both nucleon species. Furthermore, we show, for the first time, the approach to scaling behavior for g 1 measured in the resonance region at sufficiently high momentum transfer. Conclusions: Here, our results can be used to quantify the deviations from the applicability of pQCD for data taken at moderate energies and can help with extraction of quark distribution functions from such data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Measurement of W±-boson differential cross-sections in proton–proton collisions with low pile-up data at s=5.02 TeV and 13TeV with the ATLAS detector

High precision single-differential W±$$W^\pm $$-boson production cross-sections as a function of electron or muon transverse momentum pT$$p_\textrm{T}$$ or their pseudorapity η$$\eta $$, as well as double-differential cross-sections as functions of these variables, are measured in proton–proton collisions at centre-of-mass energies s=5.02$$\sqrt{s}=5.02$$ TeV and 13 TeV. The W-boson charge asymmetry as a function of lepton η$$\eta $$ is also measured. The data, collected in dedicated runs at reduced instantaneous luminosity with the ATLAS detector at the Large Hadron Collider, correspond to integrated luminosities of 255 pb-1$$^{-1}$$ at 5.02 TeV and 338 pb-1$$^{-1}$$ at 13 TeV. The measurements are in agreement with Standard-Model predictions calculated at next-to-next-to-leading-order in the strong coupling constant αs$$\alpha _s$$ including transverse-momentum resummation at next-to-next-to-leading logarithmic accuracy using several parton distribution functions. The impact of the measured differential cross-sections as a function of lepton η$$\eta $$ on the determination of these functions is studied using a profiling technique.

Aad, G↗
DOE OSTI2023

An Efficient Reactive Force Field without Explicit Coordination Dependence for Studying Caustic Aluminum Chemistry

Reactive force fields (RFFs) are an expedient approach to sample chemical reaction paths in complex systems, relative to density functional theory. However, there is continued need to improve efficiencies, specifically in systems that have slow transverse degrees of freedom, as in highly viscous and superconcentrated solutions. Here, we present an RFF that is differentiated from current models (e.g., ReaxFF) by omitting explicit dependence on the atom coordination and employing a small parameter set based on Lennard-Jones, Gaussian, and Stillinger–Weber potentials. In this study, the model was parametrized from AIMD simulation data and is used to model aluminate reactivity in sodium hydroxide solutions with extensive validation against experimental radial distribution functions, computed free energy profiles for oligomerization, and formation energies. In conclusion, the model enables simulation of early stage Al(OH) 3 nucleation which has significant relevance to industrial processing of aluminum and has a computational cost that is reduced by 1 order of magnitude relative to ReaxFF.

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

Amorphous zinc–molybdenum–sulfide chalcogel as a long-cycle, high-capacity electrode for lithium-ion batteries

The inherent limitations of intercalation-based electrodes in lithium-ion batteries have prompted the search for alternative materials with higher specific capacities and robust electrochemical stability. Sulfur-based electrodes, despite their high theoretical capacities (1672 mAh g −1 ), typically suffer from poor cycling performance. In this work, zinc molybdenum polysulfide (Zn x Mo 3 S 13 , 0.5 ≤ x), an amorphous semiconductor chalcogel, exhibits high specific capacity and excellent cycling stability. Synchrotron X-ray pair distribution function and extended X-ray absorption fine structure analyses reveal a short-range atomic structure comprising Mo–Mo, M–S (M = Mo, Zn), and S–S bonding motifs. The coordination environment of Mo and S closely resembles that of Mo 3 S 13 clusters, interconnected via S–S bridges and Zn 2+ cations. The Li/Zn x Mo 3 S 13 cell delivers an initial discharge capacity of 844 mAh g −1 at C/3, and retains 386.2 mAh g −1 after 1000 cycles with an average coulombic efficiency of 99.99%. The distribution of relaxation times analysis confirms the formation of a stable solid electrolyte interphase, which underpins the cell's long-term stability. In conclusion, this outstanding performance is attributed to the synergistic effects of the chalcogel's unique amorphous framework, semiconductive character, Zn-mediated polysulfide anchoring, and structural resilience, positioning Zn x Mo 3 S 13 chalcogel among the most durable pure metal sulfide cathodes reported for next-generation LIBs.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A model for $F_L$ structure function at low values of $Q^2$ and $x$ - revisited

A reanalysis of the model for the longitudinal structure function F_L (x,Q^2) F L ( x , Q 2 ) at low x x and low Q^2 Q 2 was undertaken, in view of the advent of the EIC. The model includes the kinematic constraint F_L\sim Q^4 F L ∼ Q 4 as Q^2\rightarrow Q 2 → 0. It is based on the photon-gluon fusion mechanism suitably extrapolated to the region of low Q^2 Q 2 . Revised model was critically updated, extended to the EIC kinematic region, and e.g. contains new parameterisations of the parton distribution functions.

Badelek, Barbara↗
DOE OSTI2022

Determination of Molecular Structure and Dynamics of Molten Salts by Advanced Neutron and X-ray Scattering Measurements and Computer Modeling

The design and development of fully functional Molten-Salt Reactors (MSR) require detailed knowledge of the molten salt properties in order to understand and predict the salt’s behavior. Fundamental properties of interest include molecular structure, speciation, and dynamics (such as diffusion coefficients) of salt components and dissolved corrosion and fission products. Computer modeling is necessary to predict changes in physical and chemical properties due to irradiation, burning of dissolved fuel, and corrosion. The modeling requires experimental data, and advanced neutron and x-ray scattering and spectroscopy provide the most reliable and direct determination of the structure (Pair-Distribution Functions, PDF), and dynamics of ions in the melt. This project dealt with both fluoride and chloride salts. The PDFs have been measured by a combination of neutron and x-ray diffraction. We utilized the techniques of isotope substitutions, a very powerful tool available for neutron-scattering, to extract the details of the liquid structure. Although similar measurements have been done before, modern advanced neutron and x-ray-scattering techniques allow collecting the data at much higher resolution and in a wider range of temperatures. Importantly, we were among the first to study fluoride salts by neutron scattering. The importance of impurities and their effects on salt properties have become apparent recently and so new methods of salt purification were developed. We took advantage of these developments to produce reliable data, which have been used for computer simulations of both clean salts and those with added fission and corrosion products most relevant for MSRs. Ab initio molecular dynamics simulations have been performed to understand the multi-component liquid solution, in particular solubility of impurities and thermodynamic interactions in relation to the ionic-cluster structure of the fluid. We applied machine learning to regress from the simulation and experimental data in order to develop a fast-acting model that can handle molten salt with an arbitrary (≥ 10) number of chemical elements and be able to predict chemical potential as a function of composition and temperature. This project resulted in a number of experimental and computer-simulation publications, a patent application, and numerous conference presentations (American Physical Society, American Chemical Society, and The Electrochemical Society among others). Multiple students and postdocs participated and collaborated on aspects of this project. This project seeded new collaborations between MIT and other institutions, such as the University of Massachusetts Lowell, the University of Illinois Urbana-Champaign, the University of California Berkeley, and Oak Ridge and Los Alamos National Labs. As such, this project has had a broad and lasting impact beyond its original scientific scope.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Temperature-Dependent Optical and Structural Properties of Chiral Two-Dimensional Hybrid Lead-Iodide Perovskites

Layered hybrid halide perovskites gain chirality via the incorporation of chiral organic cations in the spacer layer. These so-called chiral two-dimensional (2D) halide perovskites have attracted considerable interest recently for chiral optoelectronic, spintronic, and ferroelectric applications. However, the effect of temperature on these materials, especially how the structure changes with temperature and its impact on the chiral optoelectronic properties, remains an open question. Here, we study the effect of temperature change on chiral optoelectronic and structural properties through temperature-dependent circularly polarized photoluminescence (CPPL) microscopy and synchrotron powder X-ray diffraction as well as pair distribution function analysis to elucidate the intrinsic chiral optoelectronic and structural variations for R-, S-, and racemic-methylbenzylamine lead iodide. Here our results show that the temperature-induced band gap changes indicate a strong electron–phonon coupling compared to the thermal-induced lattice expansion in chiral 2D perovskites. From powder diffraction measurements, a monotonic lattice expansion is observed on heating with no structural phase transitions over 90–360 K detected for the three samples studied herein. Locally, a strongly anisotropic and even negative expansion in some components of the instantaneous Pb–I pair-distance distribution is suggestive of coupling between dynamical intralayer distortions and lattice expansion. This work provides insights into the fundamental understanding of the temperature effect on chiral optoelectronic and structural properties of chiral 2D perovskites, which can lead to new chiral materials design strategies for future chiral optoelectronic applications.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Simulating energetic ions and enhanced fusion rates from ion-cyclotron resonance heating with a full-wave/Fokker–Planck model

Reproducing fast-ion enhanced fusion rates from ion-cyclotron resonance heating (ICRH) in tokamaks requires the self-consistent coupling of a full-wave solver and a Fokker–Planck solver, which evolves multiple simultaneously resonant ion species. We introduce a new self-consistent model that iterates the TORIC full-wave solver with the CQL3D Fokker–Planck solver using the integrated plasma simulator (IPS). This model evolves the bounce-averaged ion distribution functions in both parallel and perpendicular velocity-space with a quasilinear radio frequency (RF) diffusion operator valid in the ion finite Larmor radius (FLR) limit and the RF electric fields with the resultant non-Maxwellian FLR dielectric tensor. This produces non-Maxwellian ICRH simulations that are fully self-consistent, fast, and interoperable with integrated modeling frameworks, such as TRANSP/GACODE/IPS-FASTRAN. We demonstrate our model's capabilities by validating it against experimental data in Alcator C-Mod. We then perform the first RF heating simulations of SPARC using self-consistent non-Maxwellian ion distributions to investigate the potential to enhance fusion rates using ion cyclotron resonance heating generated fast ions.

Physics↗
DOE OSTI2024

Tracking copper nanofiller evolution in polysiloxane during processing into SiOC ceramic

Polymer-derived ceramics (PDCs) remain at the forefront of research for a variety of applications including ultra-high-temperature ceramics, energy storage and functional coatings. Despite their wide use, questions remain about the complex structural transition from polymer to ceramic and how local structure influences the final microstructure and resulting properties. This is further complicated when nanofillers are introduced to tailor structural and functional properties, as nanoparticle surfaces can interact with the matrix and influence the resulting structure. The inclusion of crystalline nanofiller produces a mixed crystalline–amorphous composite, which poses characterization challenges. With this study, we aim to address these challenges with a local-scale structural study that probes changes in a polysiloxane matrix with incorporated copper nanofiller. Composites were processed at three unique temperatures to capture mixing, pyrolysis and initial crystallization stages for the pre-ceramic polymer. We observed the evolution of the nanofiller with electron microscopy and applied synchrotron X-ray diffraction with differential pair distribution function (d-PDF) analysis to monitor changes in the matrix's local structure and interactions with the nanofiller. The application of the d-PDF to PDC materials is novel and informs future studies to understand interfacial interactions between nanofiller and matrix throughout PDC processing.

Chemistry↗
DOE OSTI2022

Quantifying radial migration in the Milky Way: inefficient over short time-scales but essential to the very outer disc beyond ∼15 kpc

ABSTRACT Stellar radial migration plays an important role in reshaping a galaxy’s structure and the radial distribution of stellar population properties. In this work, we revisit reported observational evidence for radial migration and quantify its strength using the age–[Fe/H] distribution of stars across the Milky Way with APOGEE data. We find a broken age–[Fe/H] relation in the Galactic disc at r > 6 kpc, with a more pronounced break at larger radii. To quantify the strength of radial migration, we assume stars born at each radius have a unique age and metallicity, and then decompose the metallicity distribution function (MDF) of mono-age young populations into different Gaussian components that originated from various birth radii at rbirth < 13 kpc. We find that, at ages of 2 and 3 Gyr, roughly half the stars were formed within 1 kpc of their present radius, and very few stars (<5 per cent) were formed more than 4 kpc away from their present radius. These results suggest limited short-distance radial migration and inefficient long-distance migration in the Milky Way during the last 3 Gyr. In the very outer disc beyond 15 kpc, the observed age–[Fe/H] distribution is consistent with the prediction of pure radial migration from smaller radii, suggesting a migration origin of the very outer disc. We also estimate intrinsic metallicity gradients at ages of 2 and 3 Gyr of −0.061 and −0.063 dex kpc−1, respectively.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Confronting lattice parton distributions with global QCD analysis

We present the first Monte Carlo based global QCD analysis of spin-averaged and spin-dependent parton distribution functions (PDFs) that includes nucleon isovector matrix elements in coordinate space from lattice QCD. We investigate the degree of universality of the extracted PDFs when the lattice and experimental data are treated under the same conditions within the Bayesian likelihood analysis. For the unpolarized sector, we find rather weak constraints from the current lattice data on the phenomenological PDFs, and difficulties in describing the lattice matrix elements at large spatial distances. In contrast, for the polarized PDFs we find good agreement between experiment and lattice data, with the latter providing significant constraints on the spin-dependent isovector quark and antiquark distributions.

Sato, Nobuo↗
DOE OSTI2020

The chemical properties of the Milky Way’s on-bar and off-bar regions: evidence for inhomogeneous star formation history in the bulge

ABSTRACT Numerous studies of integrated starlight, stellar counts, and kinematics have confirmed that the Milky Way is a barred galaxy. However, far fewer studies have investigated the bar’s stellar population properties, which carry valuable independent information regarding the bar’s formation history. Here, we conduct a detailed analysis of chemical abundance distributions ([Fe/H] and [Mg/Fe]) in the on-bar and off-bar regions to study the azimuthal variation of star formation history (SFH) in the inner Galaxy. We find that the on-bar and off-bar stars at Galactocentric radii 3 kpc < rGC < 5 kpc have remarkably consistent [Fe/H] and [Mg/Fe] distribution functions and [Mg/Fe]–[Fe/H] relation, suggesting a common SFH shared by the long bar and the disc. In contrast, the bar and disc at smaller radii (2 kpc < rGC < 3 kpc) show noticeable differences, with relatively more very metal-rich ($\rm [Fe/H] \sim 0.4$) stars but fewer solar abundance stars in the bar. Given the three-phase star formation history proposed for the inner Galaxy in Lian et al., these differences could be explained by the off-bar disc having experienced either a faster early quenching process or recent metal-poor gas accretion. Vertical variations of the abundance distributions at small rGC suggest a wider vertical distribution of low-α stars in the bar, which may serve as chemical evidence for vertical heating through the bar buckling process. The lack of such vertical variations outside the bulge may then suggest a lack of vertical heating in the long bar.

Lian, Jianhui↗
DOE OSTI2021

Modelling the stellar halo with RR-Lyrae stars

ABSTRACT A seven-parameter distribution function (DF) is fitted to $20\, 000$ RR-Lyrae stars for which only astrometric data are available. The observational data are predicted by the DF in conjunction with the gravitational potential of a self-consistent model Galaxy defined by DFs for the dark halo, the bulge, and a four-component disc. Tests of the technique developed to deal with missing line-of-sight velocities show that adding such velocities tightens constraints on the DF only slightly. The recovered model of the RR-Lyrae population confirms that the population is flattened and has a strongly radially biased velocity distribution. At large radii, its density profile tends to ρ ∼ r−4.5 but no power law provides a good fit inside the solar sphere. The model is shown to provide an excellent fit to the data for stars brighter than r = 16.5 but at certain longitudes it predicts too few faint stars at Galactocentric radii $\sim 20\, \mathrm{kpc}$, possibly signalling that the halo is not axisymmetric. The DF is used to predict the velocity distribution of BHB stars for which space velocities are available. The z components are predicted successfully but too much anisotropy in the vRvϕ plane is expected.

Li, Chengdong↗
DOE OSTI2023

Spatiotemporal Studies of Soluble Inorganic Nanostructures with X‐rays and Neutrons

Abstract This Review addresses the use of X‐ray and neutron scattering as well as X‐ray absorption to describe how inorganic nanostructured materials assemble, evolve, and function in solution. We first provide an overview of techniques and instrumentation (both large user facilities and benchtop). We review recent studies of soluble inorganic nanostructure assembly, covering the disciplines of materials synthesis, processes in nature, nuclear materials, and the widely applicable fundamental processes of hydrophobic interactions and ion pairing. Reviewed studies cover size regimes and length scales ranging from sub‐Ångström (coordination chemistry and ion pairing) to several nanometers (molecular clusters, i.e. polyoxometalates, polyoxocations, and metal‐organic polyhedra), to the mesoscale (supramolecular assembly processes). Reviewed studies predominantly exploit 1) SAXS/WAXS/SANS (small‐ and wide‐angle X‐ray or neutron scattering), 2) PDF (pair‐distribution function analysis of X‐ray total scattering), and 3) XANES and EXAFS (X‐ray absorption near‐edge structure and extended X‐ray absorption fine structure, respectively). While the scattering techniques provide structural information, X‐ray absorption yields the oxidation state in addition to the local coordination. Our goal for this Review is to provide information and inspiration for the inorganic/materials science communities that may benefit from elucidating the role of solution speciation in natural and synthetic processes.

Yin, Jia‐Fu↗
DOE OSTI2020

Measurement of single-diffractive dijet production in proton-proton collisions at $\sqrt{s} =$ 8 TeV with the CMS and TOTEM experiments

Measurements are presented of the single-diffractive dijet cross section and the diffractive cross section as a function of the proton fractional momentum loss $\xi $ and the four-momentum transfer squared t. Both processes ${\text{ p }{}{}} {\text{ p }{}{}} \rightarrow {\text{ p }{}{}} {\text{ X }} $ and ${\text{ p }{}{}} {\text{ p }{}{}} \rightarrow {\text{ X }} {\text{ p }{}{}} $, i.e. with the proton scattering to either side of the interaction point, are measured, where ${\text{ X }} $ includes at least two jets; the results of the two processes are averaged. The analyses are based on data collected simultaneously with the CMS and TOTEM detectors at the LHC in proton–proton collisions at $\sqrt{s} = 8\,\text {Te}\text {V} $ during a dedicated run with $\beta ^{*} = 90\,\text {m} $ at low instantaneous luminosity and correspond to an integrated luminosity of $37.5{\,\text {nb}^{-1}} $. The single-diffractive dijet cross section $\sigma ^{{\text{ p }{}{}} {\text{ X }}}_{\mathrm {jj}}$, in the kinematic region $\xi < 0.1$, $0.03< |t | < 1\,\text {Ge}\text {V} ^2$, with at least two jets with transverse momentum $p_{\mathrm {T}} > 40\,\text {Ge}\text {V} $, and pseudorapidity $|\eta | < 4.4$, is $21.7 \pm 0.9\,\text {(stat)} \,^{+3.0}_{-3.3}\,\text {(syst)} \pm 0.9\,\text {(lumi)} \,\text {nb} $. The ratio of the single-diffractive to inclusive dijet yields, normalised per unit of $\xi $, is presented as a function of x, the longitudinal momentum fraction of the proton carried by the struck parton. The ratio in the kinematic region defined above, for x values in the range $-2.9 \le \log _{10} x \le -1.6$, is $R = (\sigma ^{{\text{ p }{}{}} {\text{ X }}}_{\mathrm {jj}}/\Delta \xi )/\sigma _{\mathrm {jj}} = 0.025 \pm 0.001\,\text {(stat)} \pm 0.003\,\text {(syst)} $, where $\sigma ^{{\text{ p }{}{}} {\text{ X }}}_{\mathrm {jj}}$ and $\sigma _{\mathrm {jj}}$ are the single-diffractive and inclusive dijet cross sections, respectively. The results are compared with predictions from models of diffractive and nondiffractive interactions. Monte Carlo predictions based on the HERA diffractive parton distribution functions agree well with the data when corrected for the effect of soft rescattering between the spectator partons.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Minimal Experimental Bias on the Hydrogen Bond Greatly Improves Ab Initio Molecular Dynamics Simulations of Water

Experiment Directed Simulations (EDS) is a method within a class of techniques seeking to improve molecular simulations by minimally biasing the system Hamiltonian to reproduce certain experimental observables. In a previous application of EDS to ab initio molecular dynamics (AIMD) simulation based on electronic density functional theory (DFT), the AIMD simulations of water were biased to reproduce its experimentally derived solvation structure. In particular, by solely biasing the O-O pair correlation functions, other structural and dynamical properties that were not biased were improved. In this work, the hypothesis is tested that directly biasing the OH pair correlation (and hence the H-O∙∙∙H hydrogen bonding), will provide an even better improvement of DFT-based water properties in AIMD simulations. The logic behind this hypothesis is that for most electronic DFT descriptions of water the hydrogen bonding is known to be deficient due to anomalous charge transfer and over polarization in the DFT. Using recent advances to the EDS learning algorithm, we thus train a minimal bias on AIMD water that reproduces the O-H radial distribution function derived from the highly accurate MB-pol model of water. Finally, it is then confirmed that biasing the O-H pair correlation alone can lead to improved AIMD water properties, with structural and dynamical properties in even closer to experiment than the previous EDS-AIMD model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗