Unitary Coupled Cluster: Seizing the Quantum Moment
Not Available
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
SEARCH · Search NASA
Results for “coupled cluster method”
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.
Characterizing Signatures of Geothermal Exploration Data with Machine Learning Techniques: An Application to the Nevada Play Fairway Analysis
We are introducing machine learning methods to the play fairway analysis to generate geothermal potential maps to support the evaluation of geothermal resource potential and the exploration for undiscovered blind geothermal systems in the Nevada Great Basin region. Our project aims to identify new ways to combine the play fairway data and empirically organize relationships between feature weights and labels in an improved workflow. As a means of doing this, we introduce machine learning methods to evaluate the influence of certain geological and geophysical features/feature sets in predicting geothermal favorability. This report highlights promising approaches based on supervised and unsupervised learning methods. First, we demonstrate a filter method applied to supervised classification modeling. The supervised filter method is based on permutation analysis to evaluate every possible feature combination/drop out scenario and rank feature influence based on the performance variance of supervised classification models. Additionally, we present an unsupervised factor analysis based on principal component analysis coupled with a semi-supervised kmeans clustering algorithm. This analysis allows us to identify the optimal number of groups/clusters for training sites and structural settings to identify feature patterns including correlation, variance, and latent and dominant feature relationships. The results from these methods offer a promising avenue for identifying favorable sources of predictive information to identify the locations of blind geothermal systems and furthering our understanding of complex geothermal feature and label relationships in the Great Basin region and beyond.
15 GEOTHERMAL ENERGY↗
Electronic structure of strongly correlated systems: recent developments in multiconfiguration pair-density functional theory and multiconfiguration nonclassical-energy functional theory
Strong electron correlation plays an important role in transition-metal and heavy-metal chemistry, magnetic molecules, bond breaking, biradicals, excited states, and many functional materials, but it provides a significant challenge for modern electronic structure theory. The treatment of strongly correlated systems usually requires a multireference method to adequately describe spin densities and near-degeneracy correlation. However, quantitative computation of dynamic correlation with multireference wave functions is often difficult or impractical. Multiconfiguration pair-density functional theory (MC-PDFT) provides a way to blend multiconfiguration wave function theory and density functional theory to quantitatively treat both near-degeneracy correlation and dynamic correlation in strongly correlated systems; it is more affordable than multireference perturbation theory, multireference configuration interaction, or multireference coupled cluster theory and more accurate for many properties than Kohn–Sham density functional theory. This perspective article provides a brief introduction to strongly correlated systems and previously reviewed progress on MC-PDFT followed by a discussion of several recent developments and applications of MC-PDFT and related methods, including localized-active-space MC-PDFT, generalized active-space MC-PDFT, density-matrix-renormalization-group MC-PDFT, hybrid MC-PDFT, multistate MC-PDFT, spin–orbit coupling, analytic gradients, and dipole moments. We also review the more recently introduced multiconfiguration nonclassical-energy functional theory (MC-NEFT), which is like MC-PDFT but allows for other ingredients in the nonclassical-energy functional. We discuss two new kinds of MC-NEFT methods, namely multiconfiguration density coherence functional theory and machine-learned functionals.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
The benzene radical anion: A computationally demanding prototype for aromatic anions
The benzene radical anion is studied with ab initio coupled-cluster theory in large basis sets. Unlike the usual assumption, we find that, at the level of theory investigated, the minimum energy geometry is non-planar with tetrahedral distortion at two opposite carbon atoms. The anion is well known for its instability to auto-ionization which poses computational challenges to determine its properties. Despite the importance of the benzene radical anion, the considerable attention it has received in the literature so far has failed to address the details of its structure and shape-resonance character at a high level of theory. Here, we examine the dynamic Jahn-Teller effect and its impact on the anion potential energy surface. We find that a minimum energy geometry of C{sub 2} symmetry is located below one D{sub 2h} stationary point on a C{sub 2h} pseudo-rotation surface. The applicability of standard wave function methods to an unbound anion is assessed with the stabilization method. The isotropic hyperfine splitting constants (A{sub iso}) are computed and compared to data obtained from experimental electron spin resonance experiments. Satisfactory agreement with experiment is obtained with coupled-cluster theory and large basis sets such as cc-pCVQZ.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
CASM — A software package for first-principles based study of multicomponent crystalline solids
CASM is a software package that enables first-principles based studies of crystalline materials. It has been designed to treat coupled chemical, mechanical, vibrational, and magnetic degrees of freedom to determine ground state and finite temperature properties of crystals. The symmetry of the underlying parent crystal structure is used to enumerate perturbations of the parent crystal structure and generate derivative structures which can be input to first-principles calculations in order to explore the ground state energy landscape. CASM algorithmically constructs cluster expansions that fully couple discrete and continuous degrees of freedom, and generates highly efficient code to evaluate the cluster expansion basis functions. Widely used machine learning methods are integrated for fitting expansion coefficients to first-principle calculations. The fully parameterized cluster expansions can be combined with (kinetic) Monte Carlo methods to calculate finite temperature thermodynamic and kinetic properties. CASM Alloy Manager identifies distinct parent crystal structures in an alloy system, creating individual projects for each, and integrating the results. The integrated infrastructure facilitates the linkage between first-principles statistical mechanics predictions with higher length scale computational methods, such as phase field simulations. In conclusion, CASM projects are designed to be easy to share in repositories, to re-use, and to extend to include additional chemical species or types of degrees of freedom.
36 MATERIALS SCIENCE↗
Coupled cluster calculations on TiO2 nanoclusters
The excitation energies of the four lowest-lying singlet excited states of the TiO2 Ti2O4 and Ti3O6 clusters are calculated by a variety of different Equation-of-Motion Coupled Cluster (EOM-CC) approaches in order to obtain benchmark values for the optical excitations of titanium dioxide clusters. More specifically we investigate what is the effect of the inclusion of triple excitations "triples" in the (EOM-)CC scheme on the calculated excited states of those clusters. While for the monomer and dimer the inclusion of triples causes only a rigid shift in the excitation energies, in the case of the trimer the crossing of the interested states is observed. Coupled cluster approaches where triples are treated perturbatively were found to offer no advantage over EOM-CCSD, whereas the active-space methods (EOM-CCSDt(II/I)) proved to yield results very close to the full EOM-CCSDT, but at a much lower computational cost.
calculations↗
Vibrational quasi-degenerate perturbation theory with optimized coordinates: Applications to ethylene and trans-1,3-butadiene
A perturbative extension to optimized coordinate vibrational self-consistent field (oc-VSCF) is proposed based on the quasi-degenerate perturbation theory (QDPT). A scheme to construct the degenerate space (P space) is developed, which incorporates degenerate configurations and alleviates the divergence of perturbative expansion due to localized coordinates in oc-VSCF (e.g., local O–H stretching modes of water). An efficient configuration selection scheme is also implemented, which screens out the Hamiltonian matrix element between the P space configuration (p) and the complementary Q space configuration (q) based on a difference in their quantum numbers (λ{sub pq} = ∑{sub s}|p{sub s} − q{sub s}|). It is demonstrated that the second-order vibrational QDPT based on optimized coordinates (oc-VQDPT2) smoothly converges with respect to the order of the mode coupling, and outperforms the conventional one based on normal coordinates. Furthermore, an improved, fast algorithm is developed for optimizing the coordinates. First, the minimization of the VSCF energy is conducted in a restricted parameter space, in which only a portion of pairs of coordinates is selectively transformed. A rational index is devised for this purpose, which identifies the important coordinate pairs to mix from others that may remain unchanged based on the magnitude of harmonic coupling induced by the transformation. Second, a cubic force field (CFF) is employed in place of a quartic force field, which bypasses intensive procedures that arise due to the presence of the fourth-order force constants. It is found that oc-VSCF based on CFF together with the pair selection scheme yields the coordinates similar in character to the conventional ones such that the final vibrational energy is affected very little while gaining an order of magnitude acceleration. The proposed method is applied to ethylene and trans-1,3-butadiene. An accurate, multi-resolution potential, which combines the MP2 and coupled-cluster with singles, doubles, and perturbative triples level of electronic structure theory, is generated and employed in the oc-VQDPT2 calculation to obtain the fundamental tones as well as selected overtones/combination tones coupled to the fundamentals through the Fermi resonance. The calculated frequencies of ethylene and trans-1,3-butadiene are found to be in excellent agreement with the experimental values with a mean absolute error of 8 and 9 cm{sup −1}, respectively.
74 ATOMIC AND MOLECULAR PHYSICS↗
Melting behaviour of gold-platinum nanoalloy clusters by molecular dynamics simulations
The melting behavior of bimetallic gold-platinum nanoclusters is studied by applying Brownian-type isothermal molecular dynamics (MD) simulation, a program modified from the cubic coupling scheme (CCS). The process begins with the ground-state structures obtained from global minimum search algorithm and proceeds with the investigation of the effect of temperature on the thermal properties of gold-platinum nanoalloy clusters. N-body Gupta potential has been employed in order to account for the interactions between gold and platinum atoms. The ground states of the nanoalloy clusters, which are core-shell segregated, are heated until they become thermally segregated. The detailed melting mechanism of the nanoalloy clusters is studied via this approach to provide insight into the thermal stability of the nanoalloy clusters.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Multi-scale modeling of fuel fragmentation and microstructural evolution
An overview of NEAMS program activities to model fuel fragmentation and microstructual evolution to a group of industry attendees at a meeting hosted by EPRI. Topics include atomistic methods to determine grain boundary strength, phase-field fracture modeling to determine critical pressure at grain boundaries, and coupled phase-field-cluster dynamics to model microstructural evolution.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
A Coupled-Cluster Study of the Molecular Structure, Vibrational Spectrum, and Relative Energies of the XCN and XNC (X=F, Cl) Isomers
The XCN and XNC (X=F, Cl) isomers have been investigated using the CCSD and CCSD(T) methods in conjunction with a TZ2P basis set. Equilibrium geometries, dipole moments, harmonic frequencies, IR intensities and relative energies have been evaluated. The CCSD(T) geometries and vibrational frequencies for the XCN isomers are in good agreement with the available experimental data. The CCSD(T) results for FCN and FNC are in good agreement with the CEPA calculations of Botshwina et al., with the exception of the energy difference, which the CEPA method underestimates by about 1.2 kcal/mol. FCN and CICN are shown to be lower in energy than the FNC and ClNC isomers by 69.511.0 and 42.711.0 kcal/mol (0 K), respectively.
Lee, Timothy J.↗
Character of intermolecular interaction in pyridine-argon complex: Ab initio potential energy surface, internal dynamics, and interrelations between SAPT energy components
The pyridine-Ar (PAr) van der Waals (vdW) complex is studied using a high level ab initio method. Its structure, binding energy, and intermolecular vibrational states are determined from the analytical potential energy surface constructed from interaction energy (IE) values computed at the coupled cluster level of theory with single, double, and perturbatively included triple excitations with the augmented correlation consistent polarized valence double-ζ (aug-cc-pVDZ) basis set complemented by midbond functions. The structure of the complex at its global minimum with Ar at a distance of 3.509 Å from the pyridine plane and shifted by 0.218 Å from the center of mass towards nitrogen agrees well with the corresponding equilibrium structure derived previously from the rotational spectrum of PAr. The PAr binding energy D{sub e} of 392 cm{sup −1} is close to that of 387 cm{sup −1} calculated earlier at the same ab initio level for the prototypical benzene-Ar (BAr) complex. However, under an extension of the basis set, D{sub e} for PAr becomes slightly lower than D{sub e} for BAr. The ab initio vdW vibrational energy levels allow us to estimate the reliability of the methods for the determination of the vdW fundamentals from the rotational spectra. To disclose the character of the intermolecular interaction in PAr, the symmetry-adapted perturbation theory (SAPT) is employed for the analysis of different physical contributions to IE. It is found that SAPT components of IE can be approximately expressed in the binding region by only two of them: the exchange repulsion and dispersion energy. The total induction effect is negligible. The interrelations between various SAPT components found for PAr are fulfilled for a few other complexes involving aromatic molecules and Ar or Ne, which indicates that they are valid for all rare gas (Rg) atoms and aromatics.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Testing chameleon gravity with the Coma cluster
We propose a novel method to test the gravitational interactions in the outskirts of galaxy clusters. When gravity is modified, this is typically accompanied by the introduction of an additional scalar degree of freedom, which mediates an attractive fifth force. The presence of an extra gravitational coupling, however, is tightly constrained by local measurements. In chameleon modifications of gravity, local tests can be evaded by employing a screening mechanism that suppresses the fifth force in dense environments. While the chameleon field may be screened in the interior of the cluster, its outer region can still be affected by the extra force, introducing a deviation between the hydrostatic and lensing mass of the cluster. Thus, the chameleon modification can be tested by combining the gas and lensing measurements of the cluster. We demonstrate the operability of our method with the Coma cluster, for which both a lensing measurement and gas observations from the X-ray surface brightness, the X-ray temperature, and the Sunyaev-Zel'dovich effect are available. Using the joint observational data set, we perform a Markov chain Monte Carlo analysis of the parameter space describing the different profiles in both the Newtonian and chameleon scenarios. We report competitive constraints on the chameleon field amplitude and its coupling strength to matter. In the case of f(R) gravity, corresponding to a specific choice of the coupling, we find an upper bound on the background field amplitude of |f{sub R0}| < 6 × 10{sup −5}, which is currently the tightest constraint on cosmological scales.
79 ASTRONOMY AND ASTROPHYSICS↗
On the difference between variational and unitary coupled cluster theories
There have been assertions in the literature that the variational and unitary forms of coupled cluster theory lead to the same energy functional. Numerical evidence from previous authors was inconsistent with this claim, yet the small energy differences found between the two methods and the relatively large number of variational parameters precluded an unequivocal conclusion. Using the Lipkin Hamiltonian, we here present conclusive numerical evidence that the two theories yield different energies. The ambiguities arising from the size of the cluster parameter space are absent in the Lipkin model, particularly when truncating to double excitations. We show that in the symmetry adapted basis under strong correlation, the differences between the variational and unitary models are large, whereas they yield quite similar energies in the weakly correlated regime previously explored. We also provide a qualitative argument rationalizing why these two models cannot be the same. Additionally, we study a generalized non-unitary and non-hermitian variant that contains excitation, de-excitation, and mixed operators with different amplitudes and show that it works best when compared to the traditional, variational, unitary, and extended forms of coupled cluster doubles theories
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
A Clustering-based biased Monte Carlo Approach to Protein Titration Curve Prediction
We develop and implement a novel approach to computing the ensemble averages in systems characterized by pair-wise interactions between the entities. Methods involving full enumeration of the configuration space result in exponential complexity. Sampling methods such as Markov Chain Monte Carlo (MCMC) algorithms have been proposed to tackle the exponential complexity of these problems. In certain scenarios where significant energetic coupling exists between the entities, the accuracy of the such algorithms can be diminished. We propose a strategy to improve the accuracy of the MCMC runs by taking advantage of the cluster structure in the interaction energy matrix. We propose two different schemes for performing the biased MCMC runs on the partitioned systems and show that they are valid MCMC schemes. We then apply these algorithms to the problem of computing the protonation fractions and hence the titration curves of titratable protein residues that constitute a given protein. We leverage both synthesized and real-world systems and show the improved performance of our biased MCMC methods when compared to the regular MCMC method.
Visweswara Sathanur, Arun↗
Size-extensivity-corrected multireference configuration interaction schemes to accurately predict bond dissociation energies of oxygenated hydrocarbons
Oxygenated hydrocarbons play important roles in combustion science as renewable fuels and additives, but many details about their combustion chemistry remain poorly understood. Although many methods exist for computing accurate electronic energies of molecules at equilibrium geometries, a consistent description of entire combustion reaction potential energy surfaces (PESs) requires multireference correlated wavefunction theories. Here we use bond dissociation energies (BDEs) as a foundational metric to benchmark methods based on multireference configuration interaction (MRCI) for several classes of oxygenated compounds (alcohols, aldehydes, carboxylic acids, and methyl esters). We compare results from multireference singles and doubles configuration interaction to those utilizing a posteriori and a priori size-extensivity corrections, benchmarked against experiment and coupled cluster theory. We demonstrate that size-extensivity corrections are necessary for chemically accurate BDE predictions even in relatively small molecules and furnish examples of unphysical BDE predictions resulting from using too-small orbital active spaces. We also outline the specific challenges in using MRCI methods for carbonyl-containing compounds. The resulting complete basis set extrapolated, size-extensivity-corrected MRCI scheme produces BDEs generally accurate to within 1 kcal/mol, laying the foundation for this scheme's use on larger molecules and for more complex regions of combustion PESs.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Monte Carlo configuration interaction with perturbation corrections for dissociation energies of first row diatomic molecules: C{sub 2}, N{sub 2}, O{sub 2}, CO, and NO
Dissociation energies for the diatomic molecules C{sub 2}, N{sub 2}, O{sub 2}, CO, and NO are estimated using the Monte Carlo configuration interaction (MCCI) and augmented by a second order perturbation theory correction. The calculations are performed using the correlation consistent polarized valence “triple zeta” atomic orbital basis and resulting dissociation energies are compared to coupled cluster calculations including up to triple excitations (CCSDT) and Full Configuration Interaction Quantum Monte Carlo (FCIQMC) estimates. It is found that the MCCI method readily describes the correct behavior for dissociation for the diatomics even when capturing only a relatively small fraction (∼80%) of the correlation energy. At this level only a small number of configurations, typically O(10{sup 3}) from a FCI space of dimension O(10{sup 14}), are required to describe dissociation. Including the perturbation correction to the MCCI estimates, the difference in dissociation energies with respect to CCSDT ranges between 1.2 and 3.1 kcal/mol, and the difference when comparing to FCIQMC estimates narrows to between 0.5 and 1.9 kcal/mol. Discussions on MCCI's ability to recover static and dynamic correlations and on the form of correlations in the electronic configuration space are presented.
74 ATOMIC AND MOLECULAR PHYSICS↗
Simple eigenvalue-self-consistent $\overline{Δ}GW_0$
We show that a rigid scissors-like GW self-consistency approach, labeled here $\overline{Δ}GW_0$, can be trivially implemented at zero additional cost for large scale one-shot $G_0W_0$ calculations. The method significantly improves one-shot $G_0W_0$ and for large systems is very accurate. $\overline{Δ}GW_0$ is similar in spirit to ev$GW_0$ where the self-consistency is only applied on the eigenvalues entering Green's function, while both W and the eigenvectors of Green's function are held fixed. $\overline{Δ}GW_0$ further assumes that the shift of the eigenvalues is rigid scissors-like so that all occupied states are shifted by the same amount and analogously for all the unoccupied states. We show that this results in a trivial modification of the time-dependent $G_0W_0$ self-energy, enabling an a posteriori self-consistency cycle. The method is applicable for our recent stochastic-$GW$ approach, thereby enabling self-consistent calculations for giant systems with thousands of electrons. The accuracy of $\overline{Δ}GW_0$ increases with the system size. For molecules, it is up to 0.4-0.5 eV away from coupled-cluster single double triple (CCSD(T)), but for tetracene and hexacene, it matches the ionization energies from both CCSD(T) and ev$GW_0$ to better than 0.05 eV. For solids, as exemplified here by periodic supercells of semiconductors and insulators with 6192 valence electrons, the method matches ev$GW_0$ quite well and both methods are in good agreement with the experiment.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Emergent properties of nuclei from ab initio coupled-cluster calculations
Emergent properties such as nuclear saturation and deformation, and the effects on shell structure due to the proximity of the scattering continuum and particle decay channels are fascinating phenomena in atomic nuclei. In recent years, ab initio approaches to nuclei have taken the first steps towards tackling the computational challenge of describing these phenomena from Hamiltonians with microscopic degrees of freedom. Our endeavor is now possible due to ideas from effective field theories, novel optimization strategies for nuclear interactions, ab initio methods exhibiting a soft scaling with mass number, and ever-increasing computational power. We review some of the recent accomplishments. We also present new results. The recently optimized chiral interaction NNLO${}_{{\rm{sat}}}$ is shown to provide an accurate description of both charge radii and binding energies in selected light- and medium-mass nuclei up to 56Ni. We derive an efficient scheme for including continuum effects in coupled-cluster computations of nuclei based on chiral nucleon–nucleon and three-nucleon forces, and present new results for unbound states in the neutron-rich isotopes of oxygen and calcium. Finally, the coupling to the continuum impacts the energies of the ${J}^{\pi }=1/{2}^{-},3/{2}^{-},7/{2}^{-},3/{2}^{+}$ states in ${}^{\mathrm{17,23,25}}$O, and—contrary to naive shell-model expectations—the level ordering of the ${J}^{\pi }=3/{2}^{+},5/{2}^{+},9/{2}^{+}$ states in ${}^{\mathrm{53,55,61}}$Ca.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗