Even More Efficient Quantum Computations of Chemistry Through Tensor Hypercontraction
Not Available
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
SEARCH · Engineering Papers
Results for “quantum mechanical calculations”
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.
An under-approximation of entropy for elemental multiconfigurational ground state electronic structures
A combinatorial approach has been applied to the allowable permutations of quantum electronic configurations under the constraints of Hund's rule for established ground state configurations toward an under-approximation of electronic structure entropy. Combined with a previously reported over-approximation, the approximations are used in conjunction in an attempt to bracket the upper and lower entropy limits for multiconfigurational ground state electronic structure entropy and compared to known standard molar entropies for the elements. This formality has been used for the application of a classical statistical mechanics methodology to be applied to the discrete sets of quantum mechanical states of Pu in order to calculate orbital occupancies in Pu's multiconfigurational ground state. Without consideration of the relative energies of various possible electronic configurations contributing to the multiconfigurational ground state, the calculations are performed under a general energy degeneracy assumption weighted to the number of permutations for specific configurations. The number of configurations assumed to significantly contribute is gradually constrained in order to approach a low-order approximation of orbital occupancies in Pu that are then compared to experimental and other calculated results from the literature.
Beaux, II, Miles F. (ORCID:000000032192626X)↗
Predicting Hugoniot equation of state in erythritol with ab initio and reactive molecular dynamics
Erythritol has been proposed as an inert surrogate for developing theoretical and computational models to study aging in energetic materials. In this work, we present a comparison of mechanical and shock properties of erythritol computed using the ReaxFF reactive force field and from ab initio calculations employing density functional theory (DFT). We screened eight different ReaxFF parameterizations, of which the CHO parameters developed for hydrocarbon oxidation provide the most accurate predictions of mechanical properties and the crystal structure of erythritol. Further validation of the applicability of this ReaxFF parameterization for modeling erythritol is demonstrated by comparing predictions of the elastic constants, crystal structure, vibrational density of states, and Hugoniot curves against DFT calculations. The ReaxFF predictions are in close agreement with the DFT simulations for the elastic constants and shock Hugoniot when the crystal is loaded along its c axis but show as much as 30% disagreement in the elastic constants in the ab plane and 12% difference in shock pressures when shocked along the a or b crystal axes. Last, we compare thermomechanical properties predicted from classical molecular dynamics with those calculated using the quasi-harmonic approximation and show that quantum mechanical effects produce large discrepancies in the computed values of heat capacity and thermal expansion coefficients compared with classical assumptions. Combining classical molecular dynamics predictions of mechanical behavior with phonon-based calculations of thermal behaviors, we show that predicted shock-induced temperatures for pressures up to 6.5 GPa do not exceed the pressure-dependent melting point of erythritol.
Hu, Jing↗
Ab initio calculation of infrared intensities for hydrogen peroxide
Results of an ab initio SCF quantum mechanical study are used to derive estimates for the infrared intensities of the fundamental vibrations of hydrogen peroxide. Atomic polar tensors (APTs) were calculated on the basis of a 4-31G basis set, and used to derive absolute intensities for the vibrational transitions. Comparison of the APTs calculated for H2O2 with those previously obtained for H2O and CH3OH, and of the absolute intensities derived from the H2O2 APTs with those derived from APTs transferred from H2O and CH3OH, reveals the sets of values to differ by no more than a factor of two, supporting the validity of the theoretical calculation. Values of the infrared intensities obtained correspond to A1 = 14.5 km/mol, A2 = 0.91 km/mol, A3 = 0.058 km/mol, A4 = 123 km/mol, A5 = 46.2 km/mol, and A6 = 101 km/mol. Charge, charge flux and overlap contributions to the dipole moment derivatives are also computed.
Rogers, J. D.↗
An extended trajectory-mechanics approach for calculating two-phase flow paths
A technique originating in quantum dynamics is used to derive a trajectory-based, semi-analytical solution for two-phase flow. The partial differential equation governing the evolution of the aqueous phase is equivalent to a family of ordinary differential equations defined along a path through the porous medium. The trajectories may be found by solving the differential equations directly or by post-processing the output of a numerical solution to the full set of governing equations. The trajectories, which differ from conventional streamlines, are found to bend downward in response to gravitational forces. The curvature is more pronounced as the dip of the porous layer containing the flow increases. Subtle changes in the relative permeability curve can lead to significant variations in the trajectories. The ordinary differential equation for the trajectory provides an expression for the travel time along the path. The expression produces a semi-analytical approximation to the model parameter sensitivities, the partial derivatives of the travel times with respect to changes in the permeability model. The semi-analytical trajectory-based sensitivities generally agree with those computed using a numerical reservoir simulator and a perturbation approach. The sensitivities are useful in tomographic imaging algorithms designed to estimate the spatial variation in permeability within a porous medium using multiphase observations.
58 GEOSCIENCES↗
Ab Initio Potential Energy Surfaces and the Calculation of Accurate Vibrational Frequencies
Due to advances in quantum mechanical methods over the last few years, it is now possible to determine ab initio potential energy surfaces in which fundamental vibrational frequencies are accurate to within plus or minus 8 cm(exp -1) on average, and molecular bond distances are accurate to within plus or minus 0.001-0.003 Angstroms, depending on the nature of the bond. That is, the potential energy surfaces have not been scaled or empirically adjusted in any way, showing that theoretical methods have progressed to the point of being useful in analyzing spectra that are not from a tightly controlled laboratory environment, such as vibrational spectra from the interstellar medium. Some recent examples demonstrating this accuracy will be presented and discussed. These include the HNO, CH4, C2H4, and ClCN molecules. The HNO molecule is interesting due to the very large H-N anharmonicity, while ClCN has a very large Fermi resonance. The ab initio studies for the CH4 and C2H4 molecules present the first accurate full quartic force fields of any kind (i.e., whether theoretical or empirical) for a five-atom and six-atom system, respectively.
Lee, Timothy J.↗
Benchmarking an Embedded Adaptive Sampling Configuration Interaction Method for Surface Reactions: H 2 Desorption from and CH 4 Dissociation on Cu(111)
Embedded (emb-) correlated wavefunction (CW) theory enables accurate assessments of both ground- and excited-state reaction mechanisms involved in heterogeneous catalysis. Embedded multireference second-order perturbation theory (emb-MRPT2) based on reference wavefunctions generated via embedded complete active space self-consistent field (emb-CASSCF) theory is currently state-of-the-art. However, the factorial scaling of CASSCF limits the size of active space and the complexity of systems that can be studied. In this work, we assess the efficacy of an alternative CW method, adaptive sampling configuration interaction (ASCI)–which enables large active spaces to be used–for studying surface reactions. We couple ASCI with density functional embedding theory (DFET) and benchmark its performance for two reactions: H 2 desorption from and CH 4 dissociation on the Cu(111) surface. Unlike embedded complete active space second-order perturbation theory (emb-CASPT2) that accurately reproduces a measured H 2 desorption barrier, embedded ASCI, using a very large active space (though one that still comprises a small portion of the full set of orbitals) fails to do so. Adding an extra correlation term from embedded Møller–Plesset second-order perturbation theory (emb-MP2) improves the desorption barrier and endothermicity predictions. Thus, the inaccuracy of embedded ASCI comes from the missing dynamic correlation from the many other electrons and orbitals not included in the active space. For CH 4 dissociation, again embedded ASCI overestimates the dissociation barrier compared to emb-CASPT2 predictions. Adding dynamic correlation from emb-MP2 helps correct the barrier. However, this composite approach suffers from double counting of correlation within embedded ASCI followed by emb-MP2 calculations. We therefore conclude that the state-of-the-art emb-MRPT2 based on reference wavefunctions generated via emb-CASSCF remains the method of choice for studying surface reactions. emb-ASCI is useful when large active spaces beyond the limit of emb-CASSCF are essential, such as to study complex surface reactions with significant multiconfigurational character (static correlation) but weak dynamic correlation.
08 HYDROGEN↗
Quantum chemical calculation of the equilibrium structures of small metal atom clusters
The application of ab initio quantum mechanical approaches in the study of metal atom clusters requires simplifying techniques that do not compromise the reliability of the calculations. Various aspects of the implementation of the effective core potential (ECP) technique for the removal of the metal atom core electrons from the calculation were examined. The ECP molecular integral formulae were modified to bring out the shell characteristics as a first step towards fulfilling the increasing need to speed up the computation of the ECP integrals. Work on the relationships among the derivatives of the molecular integrals that extends some of the techniques pioneered by Komornicki for the calculation of the gradients of the electronic energy was completed and a formulation of the ECP approach that quite naturally unifies the various state-of-the-art "shape- and Hamiltonian-consistent" techniques was discovered.
Kahn, L. R.↗
Experimental and theoretical studies of the He(2+)-He system - Differential cross sections for direct, single-, and double-charge-transfer scattering at keV energies
Measurements and calculations of differential cross sections for direct scattering, single-charge transfer, and double-charge transfer in collisions of 1.5-, 2.0-, 6.0-, and 10.0-keV (He-3)2+ with an He-4 target are reported. The measurements cover laboratory scattering angles below 1.5 deg with an angular resolution of about 0.03 deg. A quantum-mechanical molecular-state representation is employed in the calculations; in the case of single-charge transfer a two-state close-coupling calculation is carried out taking into account electron-translation effects. The theoretical calculations agree well with the experimental results for direct scattering and double-charge transfer. The present calculation identifies the origins of oscillatory structures observed in the differential cross sections.
Gao, R. S.↗
Quantum mechanical study of elastic scattering and rotational excitation of CO by electrons
Coupling calculations of differential, integral, and momentum transfer cross sections for pure elastic scattering and rotational excitation of CO by electron impact are reported. The calculations are based on a static charge distribution that has correct dipole and quadrupole moments, has cusps at the nuclei, and is augmented by an SCF treatment of charge polarization and a local approximation for exchange. The rotationally summed cross sections, with no adjustable parameters in the scattering calculation, are in reasonably good agreement with the experimental cross sections but are somewhat larger at small scattering angles.
Onda, K.↗
Spin Splitting Energy of Transition Metals: A New, More Affordable Wave Function Benchmark Method and Its Use to Test Density Functional Theory
Accurately predicting the spin splitting energy of chemical species is important for understanding their reactivity and magnetic properties, but it is very challenging, especially for molecules containing transition metals. One impediment to progress is the scarcity of accurate benchmark data. Here we report a set of calculations designed to yield reliable benchmarks for simple transition-metal complexes that can be used to test density functional methods that are affordable for large systems of more practical interest. Various wave function methods are tested against experiment for Fe 2+ , Fe 3+ , and Co 3+ , including CASSCF, CASPT2, CASPT3, MRCISD, MRCISD+Q, ACPF, AQCC, CCSD(T), and CASPT2/CCSD(T) and also a new method called CASPT2.5, which is performed by taking the average of the CASPT2 and CASPT3 energies. Furthermore, we find that MRCISD+Q, ACPF, and AQCC require smaller active spaces for good accuracy than are required by CASPT2 and CASPT3, and this aspect may be important for calculations on larger molecules; here we find that CASPT2.5 extrapolated to a complete basis set is the most suitable method—in terms of computational cost and in terms of accuracy on monatomic systems—and therefore we chose this method for molecular benchmarks. Then Kohn–Sham density functional calculations with 60 exchange-correlation functionals are tested for FeF 2 , FeCl 2 , and CoF 2 . We find that MN15-L, M06-SX, and revM06 have very good agreement with CASPT2.5 benchmarks in terms of both the spin splitting energy and the optimized geometry for each spin state. In addition, we recommend def2-TZVP as the most suitable basis set to perform density functional calculations for molecular spin splitting energies; extra polarization functions in the basis set do not help to increase the accuracy of the spin splitting energy in KS calculations.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Localized Active Space Pair-Density Functional Theory
Accurate quantum chemical methods for the prediction of spin-state energy gaps for strongly correlated systems are computationally expensive and scale poorly with the size of the system. This makes calculations for many experimentally interesting molecules impractical even with abundant computational resources. Previous work has shown that the localized active space (LAS) self-consistent field (SCF) method can be an efficient way to obtain multiconfiguration SCF wave functions of comparable quality to the corresponding complete active space (CAS) ones. To obtain quantitative results, a post-SCF method is needed to estimate the complete correlation energy. One such method is multiconfiguration pair-density functional theory (PDFT), which calculates the energy based on the density and on-top pair density obtained from a multiconfiguration wave function. In this work, we introduce localized-active-space PDFT, which uses a LAS wave function for subsequent PDFT calculations. The method is tested by computing spin-state energies and gaps in conjugated organic molecules and a bimetallic compound and comparing to the corresponding CAS-PDFT values.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Parallel Implementation of Nonadditive Gaussian Process Potentials for Monte Carlo Simulations
A strategy is presented to implement Gaussian process potentials in molecular simulations through parallel programming. Attention is focused on the three-body nonadditive energy, though all algorithms extend straightforwardly to the additive energy. The method to distribute pairs and triplets between processes is general to all potentials. Results are presented for a simulation box of argon, including full box and atom displacement calculations, which are relevant to Monte Carlo simulation. Data on speed-up are presented for up to 120 processes across four nodes. A 4-fold speed-up is observed over five processes, extending to 20-fold over 40 processes and 30-fold over 120 processes.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Calculated Ionization Potentials of MO 3 and MO 2 for M = U, Mo, W, and Nd
Here, ionization potentials (IPs) for MO 3 and MO 2 for M = U, Mo, W, and Nd have been predicted using the Feller–Peterson–Dixon (FPD) approach at the coupled cluster CCSD(T)/complete basis set level including additional corrections. The additional corrections are mostly small, with spin-orbit effects contributing less than 0.05 eV, except for NdO 2 where the correction lowers the IP by 0.26 eV. The IPs for UO 3 and UO 2 are calculated to be 9.59 and 6.09 eV, respectively. The calculated IPs for MoO 3 and WO 3 are very similar, 11.13 and 11.11 eV, respectively, and MoO 2 and WO 2 are 8.51 and 8.79 eV, respectively. MoO 2 has a triplet ground state, whereas WO 2 has a singlet ground state. The calculated IP for NdO 2 is 7.90 eV. NdO 3 does not achieve a high +VI formal oxidation state on the lanthanide and has an IP of 7.80 eV. These calculated IPs are expected to have error bars of ±0.04 eV.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Improving Predictions of Spin-Crossover Complex Properties through DFT Calculations with a Local Hybrid Functional
We conducted a study on the performance of the local hybrid exchange-correlation functional PBE0r for a set of 95 experimentally-characterized iron spin crossover (SCO) complexes. The PBE0r functional is a variant of PBE0 where the exchange correction is restricted to on-site terms formulated within the basis of local orbitals. We determine the free parameters of the PBE0r functional against experimental data and other hybrid functionals. With a Hartree-Fock (HF) exchange factor of 4%, the PBE0r functional accurately reproduces the electronic and free energy trends predicted in prior DFT studies for these 95 complexes using the B3LYP functional. Larger values of HF exchange stabilize high-spin states. The PBE0r-predicted bond lengths tend to exceed the experimental bond lengths, and bond lengths are less sensitive to HF exchange. The predicted SCO transition temperatures T 1/2 from PBE0r correlate moderately with the experimental transition temperatures, showing a slight improvement compared to the previous modB3LYP-predicted T 1/2 . Furthermore, this study suggests the PBE0r functional as computationally cost-effective and offers the possibility of simulating larger complexes with accuracy comparable to other global hybrid functionals, provided the HF exchange parameter is carefully optimized.
25 ENERGY STORAGE↗
Application of First-Principles-Based Artificial Neural Network Potentials to Multiscale-Shock Dynamics Simulations on Solid Materials
The use of artificial neural network (ANN) potentials trained with first-principles calculations has emerged as a promising approach for molecular dynamics (MD) simulations encompassing large space and time scales while retaining first-principles accuracy. To date, however, the application of ANN-MD has been limited to near-equilibrium processes. Here we combine first-principles-trained ANN-MD with multiscale shock theory (MSST) to successfully describe far-from-equilibrium shock phenomena. Overall, our ANN-MSST-MD approach describes shock-wave propagation in solids with first-principles accuracy but a 5000 times shorter computing time. Accordingly, ANN-MD-MSST was able to resolve fine, long-time elastic deformation at low shock speed, which was impossible with first-principles MD because of the high computational cost. This work thus lays a foundation of ANN-MD simulation to study a wide range of far-from-equilibrium processes.
36 MATERIALS SCIENCE↗
Exploiting a Shortcoming of Coupled-Cluster Theory: The Extent of Non-Hermiticity as a Diagnostic Indicator of Computational Accuracy
The fundamental non-Hermitian nature of the forms of the coupled-cluster (CC) theory widely used in quantum chemistry has usually been viewed as a negative, but the present paper shows how this can be used to an advantage. Specifically, the non-symmetric nature of the reduced one-particle density matrix (in the molecular orbital basis) is advocated as a diagnostic indicator of computational quality. In the limit of the full coupled-cluster theory [which is equivalent to full configuration interaction (FCI)], the electronic wave function and correlation energy are exact within a given one-particle basis set, and the symmetric character of the exact density matrix is recovered. The extent of the density matrix asymmetry is shown to provide a measure of “how difficult the problem is” (like the well-known T 1 diagnostic), but its variation with the level of theory also gives information about “how well this particular method works”, irrespective of the difficulty of the problem at hand. The proposed diagnostic is described and applied to a select group of small molecules, and an example of its overall utility for the practicing quantum chemist is illustrated through its application to the beryllium dimer (Be 2 ). Future application of this idea to excited states, open-shell systems, and symmetry-breaking problems and an extension of the method to the two-particle density are then proposed.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Tuning the Magnetic Properties of CrI 3 Using Ni Thin Film Deposition for Applications in Spintronic Devices
Chromium triiodide (CrI 3 ), a van der Waals magnet, has recently been shown to host Ising ferromagnetism down to the monolayer limit. It is a potentially important material in 2D magnet-based applications, such as magnetic sensors and spintronic devices. Prior studies have revealed the coexistence of two different types of interlayer magnetic coupling, the antiferromagnetic (AFM) coupling near the surface and ferromagnetic (FM) coupling in the deep bulk layers, in pristine CrI 3 crystals below the Curie temperature. In this study, we used Ni thin film deposition to tune the surface magnetic states in bulk CrI 3 . A nanometer thickness (4 nm thick) of Ni was deposited on the surface of CrI 3 using electron-beam evaporation to form a Ni/CrI 3 heterostructure. The deposited Ni thin layer forms nanoclusters that completely cover the CrI 3 surface. Magnetic states of CrI 3 before and after Ni deposition are probed by ultralow-frequency magneto-Raman spectroscopy. Instead of seeing three spin wave branches in Raman scattering below 2 T, as in pristine bulk CrI 3 , in Ni/CrI 3 , we observe only a single spin wave branch, which is softened compared to that in pristine CrI 3 and displays a Zeeman shift under an out-of-plane magnetic field up to 7 T. This observation reveals that the AFM layers on pristine CrI 3 crystal surfaces transit into FM layers, so the entire CrI 3 crystal is in the FM state after Ni deposition. First-principles calculations show that Ni atoms tend to diffuse into the CrI 3 lattice, and the FM interlayer coupling has a much lower energy than AFM coupling in the presence of Ni atoms. As a result, our studies show that the magnetic state of CrI 3 can be modified through deposition of a thin metal layer on the surface, offering a route for controlling the spin degree of freedom in van der Waals magnets.
36 MATERIALS SCIENCE↗