Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ab initio data”

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 1,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Progress in modeling hydrogen assisted ammonia oxidation with new experiments and a further reconciliation of the NH 3 + OH rate constant

Hydrogen-assisted oxidation of ammonia in a premixed, laminar flow tubular reactor under reducing conditions was investigated experimentally and through chemical kinetic modeling. Due to its impact on the competition for OH among ammonia and hydrogen, the rate constant for NH 3 + OH (R1) was determined through state-of-the-art theoretical kinetics calculations, employing composite energies that include the effects of higher order electronic excitations on the electronic energies along the variational reaction path and treating the limitations in the kinetics posed by the passage through a hydrogen-bonded complex. The resulting rate constant was in close agreement with the recent experimental value from Zaczek et al. (2025), settling a long-term dispute about the high-temperature value of k 1 and confirming within 20% the value previously used in modeling. The chemical kinetic model, with no other changes, captured well measured concentrations of NH 3 , H 2 , NO, and N 2 O from flow reactor oxidation of NH 3 /H 2 at slightly reducing conditions over a range of temperature (900-1350 K) and NH 3 /H 2 ratios (0.5-2.0). Comparison of the present results with reported data from a non-premixed setup indicates that for laminar flow tubular reactors, the reactor configuration may have implications for the observed H 2 consumption due to the possibility of preferential oxidation during mixing.

Ab initio theory↗
DOE OSTI2020

Experimental and Computational Studies on Superhard Material Rhenium Diboride under Ultrahigh Pressures

An emerging class of superhard materials for extreme environment applications are compounds formed by heavy transition metals with light elements. In this work, ultrahigh pressure experiments on transition metal rhenium diboride (ReB 2 ) were carried out in a diamond anvil cell under isothermal and non-hydrostatic compression. Two independent high-pressure experiments were carried out on ReB 2 for the first time up to a pressure of 241 GPa (volume compression V/V 0 = 0.731 ± 0.004), with platinum as an internal pressure standard in X-ray diffraction studies. The hexagonal phase of ReB 2 was stable under highest pressure, and the anisotropy between the a-axis and c-axis compression increases with pressure to 241 GPa. The measured equation of state (EOS) above the yield stress of ReB 2 is well represented by the bulk modulus K 0 = 364 GPa and its first pressure derivative K 0 ´ = 3.53. Corresponding density-functional-theory (DFT) simulations of the EOS and elastic constants agreed well with the experimental data. DFT results indicated that ReB 2 becomes more ductile with enhanced tendency towards metallic bonding under compression. The DFT results also showed strong crystal anisotropy up to the maximum pressure under study. The pressure-enhanced electron density distribution along the Re and B bond direction renders the material highly incompressible along the c-axis. Our study helps to establish the fundamental basis for anisotropic compression of ReB 2 under ultrahigh pressures

36 MATERIALS SCIENCE↗
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 OSTI2021

Inelastic rate coefficients for collisions of C 4 H - with H 2

Carbon-chain anions were recently detected in the interstellar medium. These very reactive species are used as tracers of the physical and chemical conditions in a variety of astrophysical environments. However, the local thermodynamic equilibrium conditions are generally not fulfilled in these environments. Therefore, collisional as well as radiative rates are needed to accurately model the observed emission lines. We determine in this work the state-to-state rate coefficients of C 4 H- in collision with both ortho- and para-H 2 . A new ab initio 4D potential energy surface was computed using explicitly correlated coupled-cluster procedures. This surface was then employed to determine rotational excitation and de-excitation cross-sections and rate coefficients for the first 21 rotational levels (up to rotational level j 1 = 20) using the close-coupling method, while the coupled-state approximation was used to extend the calculations up to j 1 = 30. State-to-state rate coefficients were obtained for the temperature range 2–100K. The differences between the ortho- and para-H 2 rate coefficients are found to be small.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1993

Chromospheric heating by acoustic shocks - A confrontation of GHRS observations of Alpha Tauri (K5 III) with ab initio calculations

We compare ab initio calculations of semiforbidden C II line profiles near 2325 A with recently published observations of the inactive red giant Alpha Tau (K5 III) obtained using the GHRS on board the Hubble Space Telescope. Our one-dimensional, time-dependent calculations assume that the chromosphere is heated by stochastic acoustic shocks generated by photospheric convection. We calculate various models using results from traditional (mixing length) convection zone calculations as input to hydrodynamical models. The semiforbidden C II line profiles and ratios provide sensitive diagnostics of chromospheric velocity fields, electron densities, and temperatures. We identify major differences between observed and computed line profiles which are related to basic gas dynamics and which are probably not due to technical modeling restrictions. If the GHRS observations are representative of chromospheric conditions at all epochs, then one (or more) of our model assumptions must be incorrect. Several possibilities are examined. We predict time variability of semiforbidden C II lines for comparison with observations. Based upon data from the IUE archives, we argue that photospheric motions associated with supergranulation or global pulsation modes are unimportant in heating the chromosphere of Alpha Tau.

Judge, P. G.↗
DOE OSTI2016

The interaction of reactants, intermediates and products with Cu ions in Cu-SSZ-13 NH3 SCR catalysts: An energetic and ab initio X-ray absorption modeling study

Copper-exchanged SSZ-13 (Cu-SSZ-13) provides excellent catalytic activity and hydrothermal stability in the selective catalytic reduction (SCR) of NOx by using NH3 as a reductant. In this contribution, we compare the most likely positions for Cu with a single charge compensating Al atom (ZCu), including the effect of the adsorption of reactants, intermediates and products of a NH3 SCR reaction by using first-principles calculations based on density functional theory. We find that the 6-membered ring (6MR) site is the most energetically favorable while the 8-membered ring (8MR) sites are less favorable with energy differences about 0.5 eV with respect to the 6MR site. Upon molecular adsorption, the energy differences between Cu in the 8MR and 6MR sites decreases and almost disappears, which means that the thermodynamic stability of ZCu with molecularly adsorbed species is nearly identical for Cu in different sites. Considering more complicated scenarios of NO (or CO) adsorption, we find that the co-adsorption of 2 NO (or 2 CO) molecules weakens their interaction with Cu, as well as their co-adsorption with OH and H2O. The X-ray absorption near edge structure (XANES) of Cu in Cu-SSZ-13 under different conditions was also modeled from first principles. We find that the pre-edge feature in the K-edge XANES of Cu in clean Cu-SSZ-13 with Cu in a 8MR site is induced by a 1s to 4s transition while there is no pre-edge peak when Cu is in a 6MR site. Molecular adsorption onto Cu in a 6MR or a 8MR site results in the same K-edge XANES. When NO is adsorbed onto Cu in SSZ-13, we find that the origin of the small pre-edge peak in the K-edge XANES lies in a 1s to 3d electron transition. On the other hand, when CO or N2 is adsorbed onto an isolated Cu atom in SSZ-13, a transition from the 1s state to the hybridized state of 3d and 4p results in a pre-edge peak in the K-edge of the XANES spectrum. It is also found that the edge position and intensity of the maximum peak of Cu K-edge XANES depend on the Cu locations and the adsorbed molecules. By studying the Cu K-edge XANES of ZCu with adsorption of species of different oxidizing power, we conclude that a higher oxidation state of Cu will result in a higher energy edge position in the Cu K-edge XANES. These results are expected to greatly aid our understanding of Cu-SSZ-13 SCR catalysts. This work was supported by institutional funds provided to JSM from the Voiland School of Chemical Engineering and Bioengineering. Financial support was also provided by the National Science Foundation GOALI program under contract Nos. CBET-1258717, CBET-1258715 and CBET-1258690. We thank Prof. Fabio Ribeiro, Mr. Atish Parekh, Prof. W. F. Schneider, Mr. Christopher Paolucci, Mr. Trunjoyo Anggara and Prof. Jeff Miller for stimulating discussions on the modeling of the XANES spectrum and Ms. Alyssa Hensley for her comments on the manuscript. J.Sz. acknowledges the financial support of this work by the U.S. Department of Energy (DOE), Office of Energy Efficiency and Renewable Energy, Vehicle Technologies Program. A portion of the computer time for the computational work was performed using EMSL, a national scientific user facility sponsored by the Department of Energy’s Office of Biological and Environmental Research and located at PNNL. PNNL is a multi-program national laboratory operated for the US DOE by Battelle.

Zhang, Renqin↗
DOE OSTI2016

Thermal relaxation of molecular oxygen in collisions with nitrogen atoms

Investigation of O{sub 2}–N collisions is performed by means of the quasi-classical trajectory method on the two lowest ab initio potential energy surfaces at temperatures relevant to hypersonic flows. A complete set of bound–bound and bound–free transition rates is obtained for each precollisional rovibrational state. Special attention is paid to the vibrational and rotational relaxations of oxygen as a result of chemically non-reactive interaction with nitrogen atoms. The vibrational relaxation of oxygen partially occurs via the formation of an intermediate NO{sub 2} complex. The efficient energy randomization results in rapid vibrational relaxation at low temperatures, compared to other molecular systems with a purely repulsive potential. The vibrational relaxation time, computed by means of master equation studies, is nearly an order of magnitude lower than the relaxation time in N{sub 2}–O collisions. The rotational nonequilibrium starts to play a significant effect at translational temperatures above 8000 K. The present work provides convenient relations for the vibrational and rotational relaxation times as well as for the quasi-steady dissociation rate coefficient and thus fills a gap in data due to a lack of experimental measurements for this system.

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

First-principles study of the infrared spectrum in liquid water from a systematically improved description of H-bond network

An accurate ab initio theory of the H-bond structure of liquid water requires a high-level exchange correlation approximation from density functional theory. Based on the liquid structures modeled by ab initio molecular dynamics by using maximally localized Wannier functions as a basis, we study the infrared spectrum of water within the canonical ensemble. In particular, we employ both the Perdew-Burke-Ernzerhof (PBE) functional within the generalized gradient approximation (GGA) and the state-of-the-art meta-GGA level approximation provided by the strongly constrained and appropriately normed (SCAN) functional. In this work, we demonstrate that the SCAN functional improves not only the water structure but also the theoretical infrared spectrum of water. Our analyses show that the improvement in the stretching and bending bands can be mainly attributed to better descriptions of directional H bonding and the covalency at the inter- and intramolecular levels, respectively. On the other hand, better agreements in libration and hindered translation bands are due to the improved dynamics of the H-bond network enabled by a less structured liquid in the experimental direction. The spectrum predicted by SCAN shows much better agreement with experimental data than the conventionally widely adopted PBE functional at the GGA level.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2014

Intermolecular potential energy surface and thermophysical properties of the CH{sub 4}–N{sub 2} system

A five-dimensional potential energy surface (PES) for the interaction of a rigid methane molecule with a rigid nitrogen molecule was determined from quantum-chemical ab initio calculations. The counterpoise-corrected supermolecular approach at the CCSD(T) level of theory was utilized to compute a total of 743 points on the PES. The interaction energies were calculated using basis sets of up to quadruple-zeta quality with bond functions and were extrapolated to the complete basis set limit. An analytical site-site potential function with nine sites for methane and five sites for nitrogen was fitted to the interaction energies. The PES was validated by calculating the cross second virial coefficient as well as the shear viscosity and binary diffusion coefficient in the dilute-gas limit for CH{sub 4}–N{sub 2} mixtures. An improved PES was obtained by adjusting a single parameter of the analytical potential function in such a way that quantitative agreement with the most accurate experimental values of the cross second virial coefficient was achieved. The transport property values obtained with the adjusted PES are in good agreement with the best experimental data.

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

Equilibrium carbon isotope fractionation factors of hydrocarbons: Semi-empirical force-field method

Here, we have calculated the reduced partition function ratios for carbon isotopes (β-factor) of 67 hydrocarbons (alkanes, alkenes, alkynes, cycloalkanes, and aromatics), including their 267 single-substituted isotopomers. The calculations were performed using the harmonic oscillator – rigid rotator model and the Urey / Bigeleisen-Mayer approach. Normal frequencies of molecular vibrations of the isotopologues were calculated from the molecular structures data and valence field force constants stored in the Light-handled Elucidation of Vibrations (LEV) database elaborated by Dr. Gribov and his colleagues in Russia. The LEV database is constructed from experimental spectroscopic and structural data by solving inverse problems of molecular vibrations. The LEV is internally consistent and the most comprehensive database available to date. The β-factor were calculated in the temperature range of 200–800 K with a 10 K step. Our calculations predict that the β-factors increase with an increasing number of C atoms within the same groups of hydrocarbons (e.g., alkanes). Our calculations also show a general descending order of 13 C enrichments among the different groups of hydrocarbons: cycloalkanes, aromatics, alkenes (double bonds) and isoalkanes, alkanes, alkynes (triple bonds). Position-specific, intramolecular isotope effects within hydrocarbons are determined by the β-factors of C in the different functional groups in the order: quaternary (C), tertiary (methine, CH), secondary (methylene, CH 2 ), primary (methyl, CH 3 ) and double bond (C=C), saturated bond (C - C), triple bond (HC≡). Our calculations on bulk and position-specific carbon isotope β-factors of the hydrocarbons, which are generally consistent with very limited ab initio calculations in the literature, are internally consistent and the most comprehensive to date for future applications to position-specific isotope geochemistry of hydrocarbons.

58 GEOSCIENCES↗
DOE OSTI2013

From Vibrational Spectroscopy to Force Fields and Structures of Saccharides: New Computational Algorithms and Applications

This work was undertaken with the main objective to investigate basic reactions that take place in relatively simple saccharides (mono-saccharides and cellobiose - the building block of cellulose) , in isolation and in cluster with few water molecules or with (gas-phase) clusters of few waters and ionic compounds (salt, isolated ions like H + or OH - ). Within the context of this work, different potentials were investigated; among them, were the PM3 semi empirical potential, DFT/BLYP and a new hybrid potential constructed from MP2 for the harmonic part and from adjusted Hartree-Fock anharmonic interactions (VSCF-PT2). These potentials were evaluated by comparison with experimental data from published sources and from several collaborating groups. The findings show excellent agreement between experiments and predictions with the hybrid VSCF-PT2 potential and very good agreement with predictions obtained from dynamics with dispersion corrected DFT/BLYP potential. Investigation of hydration of cellobiose, was another topic of interest. Guided by a hydration motif demonstrated by our experimental collaborators (team of Prof J.P. Simons), we demonstrated large energetic and structural differences between the two species of cellobiose: cis and trans. The later, which is dominant in solid and liquid phases, is higher in energy in the gas-phase and compared to pure water, it does not disturb as much the network of H bonds. In contrast, the cis species exhibits asymmetric hydration in cluster with up to 25 waters, indicating that it has surfactant properties. Another highlight of this research effort was the successful first time spectrometric and spectroscopic study of a gas-phase protonated sugar derivative (alpha-D-Galactopyranoside) and its interpretation by Ab Initio molecular dynamics (AIMD) simulations. The findings demonstrate the formation of a motif in which a proton bridges between two Oxygen atoms (belonging to OH groups) at the sugar; The vibrational bands involving the shared proton were strongly shifted to lower frequencies ( by about ~ 500 cm -1 for the symmetric mode, in this case). A similar motif was also observed recently by us in protonated cellobiose, indicating that this might be a common mechanism for interaction of a proton with sugars, perhaps similar to the proton wires observed in proteins. The simulations with protonated sugars also shed light on different mechanisms of interaction of a sugar with a proton, including formation of a carboxonium ion, mutarotation events, ring puckering and in the disaccharide cellobiose, the breaking of the glycosidic bond (in both forms of cis and trans). One final highlight to note in this summary, is the finding that in the monosaccharide beta-D-Galactose a hydroxyl ion abstracts a proton (forming water) in a barrierless process at room temperature, but the water remains bound to the sugar backbone, though it migrates around it; actual degradation occurs at ~500 K when water leaves the sugar. However, the study also shows that the water abstraction reaction can be reversed in the presence of 2 additional water molecules complexed with the sugar.

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

Muon capture in nuclei: An ab initio approach based on Green's function Monte Carlo methods

An ab initio Green's function Monte Carlo (GFMC) method is introduced for calculating total rates of muon weak capture in light nuclei with mass number $A ≤ 12$. As a first application of the method, we perform a calculation of the rate in $^3$$\textbf{H}$ and $^4$$\textbf{He}$ in a dynamical framework based on realistic two- and three-nucleon interactions and realistic nuclear charge-changing weak currents. The currents include one- and two-body terms induced by $π$- and $ρ$-meson exchange, and $N$-to-$Δ$ excitation, and are constrained to reproduce the empirical value of the Gamow-Teller matrix element in tritium. We investigate the sensitivity of theoretical predictions to current parametrizations of the nucleon axial and induced pseudoscalar form factors as well as to two-body contributions in the weak currents. The large uncertainties in the measured $^4$$\textbf{He}$ rates obtained from bubble-chamber experiments (carried out over 50 years ago) prevent us from drawing any definite conclusions. No data exist for $^3$$\textbf{H}$, but results are compared to those of a recent Faddeev calculation as a validation of the present GFMC method.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

The Electronic Structure of Zirconium and Hafnium Monochalcogenides

High-level ab initio CCSD(T) and spin–orbit icMRCI+Q calculations were used to predict potential energy curves (PECs) for the lowest-lying states of ZrO, ZrS, HfO, and HfS. The prediction of the ground state is basis set dependent at the icMRCI+Q level for ZrO and ZrS due to the small singlet–triplet splitting between the lowest 1 Σ + and 3 Δ states. CCSD(T) with a spin orbit correction predicted the 1 Σ + ground state in agreement with experiment. New all-electron basis sets were developed for Hf to improve the results over those predicted by use of effective core potentials (ECPs) that subsume the 4f electrons into the definition of the core. The use of the new DK-4f basis sets rather than ECPs became more important for HfO and HfS where there is a lack of a good core–valence separation. icMRCI+Q, CCSD(T), and DFT calculations for the spectroscopic parameters of ZrO, ZrS, HfO, and HfS were benchmarked with available experimental data. Bond dissociation energies (BDEs) of these four systems were calculated at the Feller–Peterson–Dixon (FPD) level to be 762.1 (ZrO), 543.5 (ZrS), 803.8 (HfO), and 575.1 kJ/mol (HfS), in excellent agreement with experiment. The HfS BDE was remeasured using the R3PI method, providing an updated experimental measurement of D 0 (HfS) = 5.978 ± 0.002 eV = 576.8 ± 0.2 kJ/mol. This experimental value, combined with experimental measurements of the ionization energies of Hf and HfS, gives the cationic BDE of D 0 (Hf + -S) = 5.124 ± 0.002 eV = 494.4 ± 0.2 kJ/mol.

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

Calculation of ion–ion mutual neutralization rate constants using Landau–Zener theory coupled with trajectory simulations for Ar + –Cl − , Br − , I −

In this computational study, we self-consistently calculate the rate constants of mutual neutralization reactions by incorporating the electron transfer probability, using Landau–Zener state transition theory with inputs derived from ab initio quantum chemistry calculations, into classical trajectory simulations. Electronic structure calculations are done using correlation consistent basis sets with multi-reference configuration interaction to map all the molecular electronic states below the ion-dissociation limit as a function of the distance between the reacting species. Our electronic structure calculations have been significantly improved from our previous work through improved selection of molecular electronic configurations maintaining a fine grid of 1a 0 over a wide range of bond lengths and accurate treatment of spin–orbit couplings. Non-adiabatic coupling matrix elements are calculated with the three-point central difference method near each avoided crossing to estimate the exact crossing point R x and coupling parameter H if , which are inputs to the multi-channel Landau–Zener theory to calculate the electron transition probability. Our approach is applied to estimate the mutual neutralization rate constants for the following ion pairs: Ar + –Cl − , Ar + –Br − , Ar + –I − at ∼133 Pa. Furthermore, our predictions are compared against the experimental data reported. It is seen that the improvement in the electronic structure calculation results in excellent agreement between the simulation results and the available experimental data to within a factor of ∼2 or ∼±50%.

Complete-active space self-consistent field↗
NASA NTRSDate not supplied

Quantum Hardware-Enabled Molecular Dynamics via Transfer Learning

The ability to perform ab initio molecular dynamics simulations using potential energy surfaces provided by quantum computers would open the door to virtually exact dynamics for a variety of chemical and biochemical systems, with impacts on catalysis and biophysics. Nonetheless, performing molecular dynamics on surfaces produced by quantum hardware has been hampered by the noisy energies typically produced by quantum computers and challenges associated with computing gradients and scaling to large systems interest. A recent set of advances in machine learning, known as transfer learning, provides a new path forward for molecular dynamics simulations on quantum hardware. Transfer learning offers a workaround, where one first trains models on larger, less accurate classical datasets and then refines them on smaller, more accurate quantum datasets. We explore this approach by training machine learning models to predict a molecule's potential energy based on its geometric structure using Behler-Parrinello neural networks. When successfully trained, the model enables energy gradient predictions necessary for dynamic simulations. To reduce the quantum resources needed, the model is initially trained with data derived from classical density functional theory and subsequently refined with a smaller dataset obtained from a variational quantum eigensolver optimization of the unitary coupled cluster ansatz. We show that this approach significantly reduces the size of the needed quantum training dataset while capturing the high accuracies needed within quantum chemistry simulations. The success of this two-step training method opens more opportunities to apply machine learning models on quantum data, a significant stride towards efficient quantum-classical hybrid computational models.

quantum computing↗
DOE OSTI2021

Statistical analysis of correlations in the x-ray induced Coulomb explosion of iodopyridine

Coulomb explosion imaging is a promising experimental tool to study individual molecules. Here, in this work simulation data for the x-ray induced Coulomb explosion of 2-iodopyridine is analyzed and the involved fragmentation dynamics are described. It is found that particular final ion momenta show correlations that reflect a collision of two atoms during the explosion. Variations of the proton kinetic energies can be associated with variations in the charge build-up time. Covariances of the forces during the explosion can be utilized to simplify the description of the dynamics in reduced dimensionality using only four collective coordinates.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2017

Shock equation of state of LiH 6 to 1.1 TPa

Using laser-generated shock waves, we have measured pressure, density, and temperature of LiH on the principal Hugoniot between 260 and 1100 GPa (2.6–11 Mbar) and on a second-shock Hugoniot up to 1400 GPa to near fivefold compression, extending the maximum pressure reached in non-nuclear experiments by a factor of two. Here, we observe the onset of metal-like reflectivity consistent with temperature-induced ionization of the Li 2s electron, and no sign of additional changes in ionization up to the maximum pressure. Our measurements are in good agreement with gas gun, Z-machine, and underground test data and are accurately described by quantum molecular dynamics simulations. The results confirm the validity of equation of state models built on an average-atom description of the electron-thermal contribution to the free energy and a density-dependent Grüneisen parameter to describe shock response of LiH over this pressure range.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

The origins of intra- and inter-molecular vibrational couplings: A case study of H{sub 2}O-Ar on full and reduced-dimensional potential energy surface

The origin and strength of intra- and inter-molecular vibrational coupling is difficult to probe by direct experimental observations. However, explicitly including or not including some specific intramolecular vibrational modes to study intermolecular interaction provides a precise theoretical way to examine the effects of anharmonic coupling between modes. In this work, a full-dimension intra- and inter-molecular ab initio potential energy surface (PES) for H{sub 2}O–Ar, which explicitly incorporates interdependence on the intramolecular (Q{sub 1}, Q{sub 2}, Q{sub 3}) normal-mode coordinates of the H{sub 2}O monomer, has been calculated. In addition, four analytic vibrational-quantum-state-specific PESs are obtained by least-squares fitting vibrationally averaged interaction energies for the (v{sub 1}, v{sub 2}, v{sub 3}) = (0, 0, 0), (0, 0, 1), (1, 0, 0), (0, 1, 0) states of H{sub 2}O to the three-dimensional Morse/long-range potential function. Each vibrationally averaged PES fitted to 442 points has root-mean-square (rms) deviation smaller than 0.15 cm{sup −1}, and required only 58 parameters. With the 3D PESs of H{sub 2}O–Ar dimer system, we employed the combined radial discrete variable representation/angular finite basis representation method and Lanczos algorithm to calculate rovibrational energy levels. This showed that the resulting vibrationally averaged PESs provide good representations of the experimental infrared data, with rms discrepancies smaller than 0.02 cm{sup −1} for all three rotational branches of the asymmetric stretch fundamental transitions. The infrared band origin shifts associated with three fundamental bands of H{sub 2}O in H{sub 2}O–Ar complex are predicted for the first time and are found to be in good agreement with the (extrapolated) experimental values. Upon introduction of additional intramolecular degrees of freedom into the intermolecular potential energy surface, there is clear spectroscopic evidence of intra- and intermolecular vibrational couplings.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗