Enhancing the Accuracy and Efficiency of Sample-Based Quantum Diagonalization with Phaseless Auxiliary-Field Quantum Monte Carlo
Supplemental Information for 10.1021/acs.jctc.5c01407
SEARCH · Engineering Papers
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.
Supplemental Information for 10.1021/acs.jctc.5c01407
Supplemental Information for 10.1021/acs.jctc.5c01407
Benchmark scalar-relativistic delta-coupled-cluster calculations of hetero-site double core ionization energies of small molecules containing second-row elements are reported. The present study has focused on the high-spin triplet components of two-site double core-ionized states, which are single reference in character and consistent with the use of standard coupledcluster methods. Contributions to computed double core ionization energies from electron-correlation and basis-set effects as well as corrections to the core-valence separation approximation have been analyzed. Furthermore, on the basis of systematic convergence of computational results with respect to these effects, delta-coupled-cluster calculations have been shown to be capable of providing accurate double core ionization energies with remaining errors estimated to be below 0.3 eV, and thus they are recommended for use to facilitate experimental studies of two-site double core-ionized states that are involved in X-ray pump/X-ray probe studies of electronic and molecular dynamics following inner shell ionization or excitation.
Not provided.
Understanding the equilibrium conditions at the metal oxide/aqueous interface is a key component toward visualizing the structure of water in confined environments and differentiating the catalytic activity of transition-metal oxides. While ambient pressure X-ray photoelectron spectroscopy (AP-XPS) has been the primary technique to investigate the formation of a hydration layer on many surfaces, results over the extended relative humidity (RH) range accessible experimentally have not been compared quantitatively to theoretical predictions. With the use of first-principles theoretical methods and accumulated knowledge of AP-XPS spectral analysis, we do so here for a model surface, TiO 2 -terminated undoped SrTiO 3 (100) (STO). The measured distribution of OH and H 2 O coverages from vacuum up to the first hydration layer is in good agreement with a static density functional theory (DFT) configuration involving partial dissociation of H 2 O per Ti-atom mediated by H-bonding. Furthermore, ab initio molecular dynamics (AIMD) simulations at 300 K for select coverages (1/4, 1/2, and 1 ML) test the role of fluctuations and entropy in the competition between adsorption and dissociation with coverage. This comparison between theory and experiment for OH and H 2 O coverages on STO provides a foundation for a more quantitative assessment of the first hydration layer and associated competition between adsorption, dissociation, and H-bonding on transition-metal oxide surfaces.
Not Available
Understanding molecular adsorption on surfaces underpins many problems in chemistry and materials science. Accurately and efficiently describing the adsorption has been a challenging task for first-principles methods as the process can involve both short-range chemical bond formations and long-range physical interactions, e.g., van der Waals (vdW) interaction. Density functional theory presents an appealing choice for modeling adsorption reactions, although calculations with many exchange-correlation density functional approximations struggle to accurately describe both chemical and physical molecular adsorptions. Here, we propose an efficient density functional approximation that is accurate for both chemical and physical adsorption by concurrently optimizing its semilocal component and the long-range vdW correction against the prototypical adsorption CO/Pt(111) and Ar 2 binding energy curve. The resulting function opens the door to accurate and efficient modeling of general molecular adsorption.
In the original manuscript, the value of parameter c was not explicitly provided. We wish to clarify that c was refitted to the hydrogen atom exchange energy, using the updated value of κ. With the updated κ = 0.2501, the optimized value of the parameter is c = 0.391596228. These values supersede any implicit assumptions about c in the original manuscript, and readers are encouraged to use these parameters for reproducing the results presented therein.
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.
Not Available
We present a transferable force field for water in proton-exchanged, alkali (Li, Na, K, Rb, and Cs) metal-exchanged, and alkaline-earth (Mg, Ca, Sr, and Ba) metal-exchanged zeolites. The fitting methodology is based on adsorbate–adsorbent interaction energies obtained from periodic density functional theory calculations and corrected using the coupled-cluster method applied to small model clusters. To ensure an accurate prediction of both adsorption and diffusion properties of water, sets of configurations that sample both adsorption sites and intracrystalline hopping transition states were used in the fitting. The quality of the force field is assessed for a wide range of zeolites with different topologies and chemical compositions, demonstrating good agreement between theoretical predictions and experimental measurements of water adsorption and diffusion.
Abstract We have developed a new procedure for combining lists of substorm onset times from multiple sources. We apply this procedure to observational data and to magnetohydrodynamic (MHD) model output from 1–31 January 2005. We show that this procedure is capable of rejecting false positive identifications and filling data gaps that appear in individual lists. The resulting combined onset lists produce a waiting time distribution that is comparable to previously published results, and superposed epoch analyses of the solar wind driving conditions and magnetospheric response during the resulting onset times are also comparable to previous results. Comparison of the substorm onset list from the MHD model to that obtained from observational data reveals that the MHD model reproduces many of the characteristic features of the observed substorms, in terms of solar wind driving, magnetospheric response, and waiting time distribution. Heidke skill scores show that the MHD model has statistically significant skill in predicting substorm onset times.
Simultaneous calculation of atmospheric processes is faster than calculating processes one at a time. This type of parallelism is beneficial or perhaps even necessary to provide good performance on modern supercomputers, which achieve faster performance through increased processor count rather than improved clock speed. The scalability of the Energy Exascale Earth System Model (E3SM) Atmosphere Model (EAM) is limited by the fluid dynamics which scales up to the number of mesh cells in the global mesh. In contrast, the suite of physics parameterizations in EAM is scalable up to the total number of physics columns, which is an order of magnitude greater than the number of mesh cells. A proposed solution to unlocking the greater potential performance from the physics suite is to solve the physics and dynamics in parallel. This work represents a first attempt at parallel splitting of the grid-scale fluid dynamics model and the subgrid-scale physics parameterizations in a global atmosphere model. We will demonstrate that switching to parallel physics-dynamics coupling extends the scalability of the EAM to up to 3 times the previous peak scalability limit and is up to 20% faster than the sequentially split coupling at the highest core counts and the same time step. Decadal simulations of both coupling approaches show very little impact to the model climate. This improved performance does not come without drawbacks, however. Parallel splitting requires a shorter time step and other modifications which largely offset performance gains. A mass fixer is required for conservation. Techniques for mitigating these issues are also discussed.
Predictions of hydrologic variables across the entire water cycle have significant value for water resources management as well as downstream applications such as ecosystem and water quality modeling. Recently, purely data-driven deep learning models like long short-term memory (LSTM) showed seemingly insurmountable performance in modeling rainfall runoff and other geoscientific variables, yet they cannot predict untrained physical variables and remain challenging to interpret. Here, we show that differentiable, learnable, process-based models (called δ models here) can approach the performance level of LSTM for the intensively observed variable (streamflow) with regionalized parameterization. We use a simple hydrologic model HBV as the backbone and use embedded neural networks, which can only be trained in a differentiable programming framework, to parameterize, enhance, or replace the process-based model's modules. Without using an ensemble or post-processor, δ models can obtain a median Nash-Sutcliffe efficiency of 0.732 for 671 basins across the USA for the Daymet forcing data set, compared to 0.748 from a state-of-the-art LSTM model with the same setup. For another forcing data set, the difference is even smaller: 0.715 versus 0.722. Meanwhile, the resulting learnable process-based models can output a full set of untrained variables, for example, soil and groundwater storage, snowpack, evapotranspiration, and baseflow, and can later be constrained by their observations. Both simulated evapotranspiration and fraction of discharge from baseflow agreed decently with alternative estimates. The general framework can work with models with various process complexity and opens up the path for learning physics from big data.