ABSTRACT
The multireference coupled-cluster Monte Carlo (MR-CCMC) algorithm is a determinant-based quantum Monte Carlo (QMC) algorithm that is conceptually similar to Full Configuration Interaction QMC (FCIQMC). It has been shown to offer a balanced treatment of both static and dynamic correlation while retaining polynomial scaling, although application to large systems with significant strong correlation remained impractical. In this paper, we document recent algorithmic advances that enable rapid convergence and a more black-box approach to the multireference problem. These include a logarithmically scaling metric-tree-based excitation acceptance algorithm to search for determinants connected to the reference space at the desired excitation level and a symmetry-screening procedure for the reference space. We show that, for moderately sized reference spaces, the new search algorithm brings about an approximately 8-fold acceleration of one MR-CCMC iteration, while the symmetry screening procedure reduces the number of active reference space determinants with essentially no loss of accuracy. We also introduce a stochastic implementation of an approximate wall projector, which is the infinite imaginary time limit of the exponential projector, using a truncated expansion of the wall function in Chebyshev polynomials. Notably, this wall-Chebyshev projector can be used to accelerate any projector-based QMC algorithm. We show that it requires significantly fewer applications of the Hamiltonian to achieve the same statistical convergence. We benchmark these acceleration methods on the beryllium and carbon dimers, using initiator FCIQMC and MR-CCMC with basis sets up to cc-pVQZ quality.
ABSTRACT
The recent developments of quantum computing present novel potential pathways for quantum chemistry as the scaling of the computational power of quantum computers could be harnessed to naturally encode and solve electronic structure problems. Theoretically exact quantum algorithms for chemistry have been proposed (e.g., quantum phase estimation), but the limited capabilities of current noisy intermediate-scale quantum devices motivated the development of less demanding hybrid algorithms. In this context, the variational quantum eigensolver (VQE) algorithm was successfully introduced as an effective method to compute the ground-state energies of small molecules. This study investigates the folded spectrum (FS) method as an extension of the VQE algorithm for the computation of molecular excited states. It provides the possibility of directly computing excited states around a selected target energy using the same quantum circuit as for the ground-state calculation. Inspired by the variance-based methods from the quantum Monte Carlo literature, the FS method minimizes the energy variance, thus, in principle, requiring a computationally expensive squared Hamiltonian to be applied. We alleviate this potentially poor scaling by employing a Pauli grouping procedure to identify sets of commuting Pauli strings that can be evaluated simultaneously. This allows for a significant reduction in the computational cost. We applied the FS-VQE method to small molecules (H2, LiH), obtaining all electronic excited states with chemical accuracy on ideal quantum simulators. Furthermore, we explore the application of quantum error mitigation techniques, demonstrating improved energy accuracy on noisy simulators compared with simulations without mitigation.
ABSTRACT
We recently introduced an efficient methodology to perform density-corrected Hartree-Fock density functional theory [DC(HF)-DFT] calculations and an extension to it we called "corrected" HF DFT [C(HF)-DFT] [Graf and Thom, J. Chem. Theory Comput. 19 5427-5438 (2023)]. In this work, we take a further step and combine C(HF)-DFT, augmented with a straightforward orbital energy correction, with the random phase approximation (RPA). We refer to the resulting methodology as corrected HF RPA [C(HF)-RPA]. We evaluate the proposed methodology across various RPA methods: direct RPA (dRPA), RPA with an approximate exchange kernel, and RPA with second-order screened exchange. C(HF)-dRPA demonstrates very promising performance; for RPA with exchange methods, on the other hand, we often find over-corrections.
ABSTRACT
The crucial step in density-corrected Hartree-Fock density functional theory (DC(HF)-DFT) is to decide whether the density produced by the density functional for a specific calculation is erroneous and, hence, should be replaced by, in this case, the HF density. We introduce an indicator, based on the difference in noninteracting kinetic energies between DFT and HF calculations, to determine when the HF density is the better option. Our kinetic energy indicator directly compares the self-consistent density of the analyzed functional with the HF density, is size-intensive, reliable, and most importantly highly efficient. Moreover, we present a procedure that makes best use of the computed quantities necessary for DC(HF)-DFT by additionally evaluating a related hybrid functional and, in that way, not only "corrects" the density but also the functional itself; we call that procedure corrected Hartree-Fock density functional theory (C(HF)-DFT).
ABSTRACT
We investigate the possibility of using a transcorrelated (TC) Hamiltonian to describe electron correlation. A method to obtain TC wavefunctions was developed based on the mathematical framework of the bi-variational principle. This involves the construction of an effective TC Hamiltonian matrix, which can be solved in a self-consistent manner. This was optimized using a method we call second-order-moment minimization and demonstrate that it is possible to obtain highly accurate energies for some closed-shell atoms and helium-like ions. The effects of certain correlator terms on the description of electron-electron and electron-nuclear cusps were also examined graphically, and some TC wavefunctions were compared against near-exact Hylleraas wavefunctions.
ABSTRACT
The development of multireference coupled cluster (MRCC) techniques has remained an open area of study in electronic structure theory for decades due to the inherent complexity of expressing a multiconfigurational wavefunction in the fundamentally single-reference coupled cluster framework. The recently developed multireference-coupled cluster Monte Carlo (mrCCMC) technique uses the formal simplicity of the Monte Carlo approach to Hilbert space quantum chemistry to avoid some of the complexities of conventional MRCC, but there is room for improvement in terms of accuracy and, particularly, computational cost. In this paper, we explore the potential of incorporating ideas from conventional MRCC-namely, the treatment of the strongly correlated space in a configuration interaction formalism-to the mrCCMC framework, leading to a series of methods with increasing relaxation of the reference space in the presence of external amplitudes. These techniques offer new balances of stability and cost against accuracy, as well as a means to better explore and better understand the structure of solutions to the mrCCMC equations.
ABSTRACT
Electrochemical carbon dioxide capture recently emerged as a promising alternative approach to conventional energy-intensive carbon-capture methods. A common electrochemical capture approach is to employ redox-active molecules such as quinones. Upon electrochemical reduction, quinones become activated for the capture of CO2 through a chemical reaction. A key disadvantage of this method is the possibility of side-reactions with oxygen, which is present in almost all gas mixtures of interest for carbon capture. This issue can potentially be mitigated by fine-tuning redox potentials through the introduction of electron-withdrawing groups on the quinone ring. In this article, we investigate the thermodynamics of the electron transfer and chemical steps of CO2 capture in different quinone derivatives with a range of substituents. By combining density functional theory calculations and cyclic voltammetry experiments, we support a previously described trade-off between the redox potential and the strength of CO2 capture. We show that redox potentials can readily be tuned to more positive values to impart stability to oxygen, but significant decreases in CO2 binding free energies are observed as a consequence. Our calculations support this effect for a large series of anthraquinones and benzoquinones. Different trade-off relationships were observed for the two classes of molecules. These trade-offs must be taken into consideration in the design of improved redox-active molecules for electrochemical CO2 capture.
ABSTRACT
In quantum Monte Carlo (QMC) methods, energy estimators are calculated as (functions of) statistical averages of quantities sampled during a calculation. Associated statistical errors of these averages are often estimated. This error estimation is not straightforward and there are several choices of the error estimation methods. We evaluate the performance of three methods (the Straatsma method, an autoregressive model, and a blocking analysis based on von Neumann's ratio test for randomness) for the energy time series given by three QMC methods [diffusion Monte Carlo, full configuration interaction Quantum Monte Carlo (FCIQMC), and coupled cluster Monte Carlo (CCMC)]. From these analyses, we describe a hybrid analysis method which provides reliable error estimates for a series of various lengths of FCIQMC and CCMC's time series. Equally important is the estimation of the appropriate start point of the equilibrated phase. We establish that a simple mean squared error rule method as described by White [K. P. White, Jr., Simulation 69(6), 323 (1997)10.1177/003754979706900601] can provide reasonable estimations.
ABSTRACT
Current nonorthogonal configuration interaction (NOCI) methods often use a set of self-consistent field (SCF) states selected based on chemical intuition. However, it may be challenging to track these SCF states across a dissociation profile and the NOCI states recovered may be spin contaminated. In this Article, we propose a method of applying spin rotation on symmetry broken unrestricted Hartree-Fock (sb-UHF) states to generate a basis for NOCI. The dissociation of ethene was examined by localizing spin rotation on each resulting carbene fragment. We show that this gives a size-consistent description of its dissociation and results in spin-pure states at all geometries. The dissociation was also studied with different orbitals, namely, canonical UHF and absolutely localized molecular orbitals (ALMO). Furthermore, we demonstrate that the method can be used to restore spin symmetry of symmetry broken SCF wave functions for molecules of various sizes, marking an improvement over existing NOCI methods.
ABSTRACT
Unitary coupled cluster (UCC), originally developed as a variational alternative to the popular traditional coupled cluster method, has seen a resurgence as a functional form for use on quantum computers. However, the number of excitors present in the Ansatz often presents a barrier to implementation on quantum computers. Given the natural sparsity of wavefunctions obtained from quantum Monte Carlo methods, we consider here a stochastic solution to the UCC problem. Using the coupled cluster Monte Carlo framework, we develop cluster selection schemes that capture the structure of the UCC wavefunction, as well as its Trotterized approximation, and use these to solve the corresponding projected equations. Due to the fast convergence of the equations with order in the cluster expansion, this approach scales polynomially with the size of the system. Unlike traditional UCC implementations, our approach naturally produces a non-variational estimator for the energy in the form of the projected energy. For unitary coupled cluster singles and doubles (UCCSD) in small systems, we find that this agrees well with the expectation value of the energy and, in the case of two electrons, with full configuration interaction results. For the larger N2 system, the two estimators diverge, with the projected energy approaching the coupled cluster result, while the expectation value is close to results from traditional UCCSD.
ABSTRACT
Self-consistent-field (SCF) approximations formulated using Hartree-Fock (HF) or Kohn-Sham density-functional theory (KS-DFT) have the potential to yield multiple solutions. However, the formal relationship between multiple solutions identified using HF or KS-DFT remains generally unknown. We investigate the connection between multiple SCF solutions for HF or KS-DFT by introducing a parameterized functional that scales between the two representations. Using the hydrogen molecule and a model of electron transfer, we continuously map multiple solutions from the HF potential to a KS-DFT description. We discover that multiple solutions can coalesce and vanish as the functional changes, forming a direct analogy with the disappearance of real HF solutions along a change in molecular structure. To overcome this disappearance of solutions, we develop a complex-analytic extension of DFT-the "holomorphic DFT" approach-that allows every SCF stationary state to be analytically continued across all molecular structures and exchange-correlation functionals.
ABSTRACT
We present a detailed discussion of our novel diagrammatic coupled cluster Monte Carlo (diagCCMC) [Scott et al. J. Phys. Chem. Lett. 10, 925 (2019)]. The diagCCMC algorithm performs an imaginary-time propagation of the similarity-transformed coupled cluster Schrödinger equation. Imaginary-time updates are computed by the stochastic sampling of the coupled cluster vector function: each term is evaluated as a randomly realized diagram in the connected expansion of the similarity-transformed Hamiltonian. We highlight similarities and differences between deterministic and stochastic linked coupled cluster theory when the latter is re-expressed as a sampling of the diagrammatic expansion and discuss details of our implementation that allow for a walker-less realization of the stochastic sampling. Finally, we demonstrate that in the presence of locality, our algorithm can obtain a fixed errorbar per electron while only requiring an asymptotic computational effort that scales quartically with system size, independent of the truncation level in coupled cluster theory. The algorithm only requires an asymptotic memory cost scaling linearly, as demonstrated previously. These scaling reductions require no ad hoc modifications to the approach.
ABSTRACT
Nonorthogonal multireference methods can predict statically correlated adiabatic energies while providing chemical insight through the combination of diabatic reference states. However, reaching quantitative accuracy using nonorthogonal multireference expansions remains a significant challenge. In this work, we present the first rigorous perturbative correction to nonorthogonal configuration interaction, allowing the remaining dynamic correlation to be reliably computed. Our second-order "NOCI-PT2" theory exploits a zeroth-order generalized Fock Hamiltonian and builds the first-order interacting space using single and double excitations from each reference determinant. This approach therefore defines the rigorous nonorthogonal extension to conventional multireference perturbation theories. We find that NOCI-PT2 can quantitatively predict multireference potential energy surfaces and provides state-specific ground and excited states for adiabatic avoided crossings. Furthermore, we introduce an explicit imaginary-shift formalism requiring shift values that are an order of magnitude smaller than those used in conventional multireference perturbation theories.
ABSTRACT
Correction for 'Ionic-caged heterometallic bismuth-platinum complex exhibiting electrocatalytic CO2 reduction' by Takefumi Yoshida et al., Dalton Trans., 2020, 49, 2652-2660.
ABSTRACT
An air-stable heterometallic Bi-Pt complex with the formula [BiPt(SAc)5]n (1; SAc = thioacetate) was synthesized. The crystal structure, natural bond orbital (NBO) and local orbital locator (LOL) analyses, localized orbital bonding analysis (LOBA), and X-ray absorption fine structure (XAFS) measurements were used to confirm the existence of Bi-Pt bonding and an ionic cage of O atoms surrounding the Bi ion. From the cyclic voltammetry (CV) and controlled potential electrolysis (CPE) experiments, 1 in tetrahydrofuran reduced CO2 to CO, with a faradaic efficiency (FE) of 92% and a turnover frequency (TOF) of 8 s-1 after 30 min of CPE at -0.79 V vs. NHE. The proposed mechanism includes an energetically favored pathway via the ionic cage, which is supported by the results of DFT calculations and reflectance infrared spectroelectrochemistry data.
ABSTRACT
Acetato-bridged palladium-lanthanide tetranuclear heterometallic complexes of the form [Pd2 Ln2 (H2 O)2 (CH3 COO)10 ]â 2 CH3 COOH [Ln2 =Ce2 (1), Pr2 (2), Nd2 (3), Sm2 (4), Tb2 (5), Dy2 (6), Dy0.2 Y1.8 (6''), Ho2 (7), Er2 (8), Er0.24 Y1.7 (8''), Tm2 (9), Yb2 (10), Y2 (11)] were synthesised and characterised by experimental and theoretical techniques. All complexes containing Kramers lanthanide ions [Ln3+ =Ce (1), Nd (3), Sm (4), Dy (6), DyY (6''), Er (8), ErY (8''), Yb (10)] showed field-induced slow magnetic relaxation, characteristic of single-molecule magnetism and purely of molecular origin. In contrast, all non-Kramers lanthanide ions [Ln3+ =Pr (2), Tb (5), Ho (7), Tm (9), Y3+ (11) is diamagnetic and non-lanthanide] did not show any slow magnetic relaxation. The variation in the electronic structure and accompanying consequences across the complexes representing all Kramers and non-Kramers lanthanide ions were investigated. The origin of the magnetic properties and the extent to which the axial donor-acceptor interaction involving the lanthanide ions and an electron-deficient d z 2 orbital of palladium affects the observed magnetic and electronic properties across the lanthanide series are presented. Unique consistent electronic and magnetic properties of isostructural complexes spanning the lanthanide series with properties dependent on whether the ions are Kramers or non-Kramers are reported.
ABSTRACT
The convergence of full configuration interaction quantum Monte Carlo (FCIQMC) is accelerated using a quasi-Newton propagation (QN) which can also be applied to coupled cluster Monte Carlo (CCMC). The computational scaling of this optimized propagation is [Formula: see text], keeping the additional computational cost to a bare minimum. Its effects are investigated deterministically and stochastically on a model system, the uniform electron gas, with Hilbert space size up to 1040, and shown to accelerate convergence of the instantaneous projected energy by over an order of magnitude in the FCIQMC test case. Its capabilities are then demonstrated with FCIQMC on an archetypical quantum chemistry problem, the chromium dimer, in an all-electron basis set with Hilbert space size of about 1022 yielding highly accurate FCI energies.
ABSTRACT
The first three-dimensional (3D) conductive single-ion magnet (SIM), (TTF)2 [Co(pdms)2 ] (TTF=tetrathiafulvalene and H2 pdms=1,2-bis(methanesulfonamido)benzene), was electrochemically synthesised and investigated structurally, physically, and theoretically. The similar oxidation potentials of neutral TTF and the molecular precursor [HNEt3 ]2 [M(pdms)2 ] (M=Co, Zn) allow for multiple charge transfers (CTs) between the SIM donor [M(pdms)2 ]n- and the TTF.+ acceptor, as well as an intradonor CT from the pdms ligand to Co ion upon electrocrystallisation. Usually TTF functions as a donor, whereas in our system TTF is both a donor and an accepter because of the similar oxidation potentials. Furthermore, the [M(pdms)2 ]n- donor and TTF.+ acceptor are not segregated but strongly interact with each other, contrary to reported layered donor-acceptor electrical conductors. The strong intermolecular and intramolecular interactions, combined with CT, allow for relatively high electrical conductivity even down to very low temperatures. Furthermore, SIM behaviour with slow magnetic relaxation and opening of hysteresis loops was observed. (TTF)2 [Co(pdms)2 ] (2-Co) is an excellent building block for preparing new conductive SIMs.