Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userStochastic Coupled Cluster Theory: Efficient Sampling of the Coupled Cluster Expansion
114
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We consider the sampling of the coupled cluster expansion within stochastic coupled cluster theory. Observing the limitations of previous approaches due to the inherently non-linear behaviour of a coupled cluster wavefunction representation we propose new approaches based upon an intuitive, well-defined condition for sampling weights and on sampling the expansion in cluster operators of different excitation levels. We term these modifications even and truncated selection respectively. Utilising both approaches demonstrates dramatically improved calculation stability as well as reduced computational and memory costs. These modifications are particularly effective at higher truncation levels owing to the large number of terms within the cluster expansion that can be neglected, as demonstrated by the reduction of the number of terms to be sampled at the level of CCSDT by 77% and at CCSDTQ56 by 98%.
rate research
Read More
We describe a modification of the stochastic coupled cluster algorithm that allows the use of multiple reference determinants. By considering the secondary references as excitations of the primary reference and using them to change the acceptance criteria for selection and spawning, we obtain a simple form of stochastic multireference coupled cluster which preserves the appealing aspects of the single reference approach. The method is able to successfully describe strongly correlated molecular systems using few references and low cluster truncation levels, showing promise as a tool to tackle strong correlation in more general systems. Moreover, it allows simple and comprehensive control of the included references and excitors thereof, and this flexibility can be taken advantage of to gain insight into some of the inner workings of established electronic structure methods.
A variant of coupled-cluster theory is described here, wherein the degrees of freedom are fluctuations of fragments between internally correlated states. The effects of intra-fragment correlation on the inter-fragment interaction are pre-computed and permanently folded into an effective Hamiltonian, thus avoiding redundant evaluations of local relaxations associated with coupled fluctuations. A companion article shows that a low-scaling step may be used to cast the electronic Hamiltonians of real systems into the form required. Two proof-of-principle demonstrations are presented here for non-covalent interactions. One uses harmonic oscillators, for which accuracy and algorithm structure can be carefully controlled in comparisons. The other uses small electronic systems (Be atoms) to demonstrate compelling accuracy and efficiency, also when inter-fragment electron exchange and charge transfer must be handled. Since the cost of the global calculation does not depend directly on the correlation models used for the fragments, this should provide a way to incorporate difficult electronic structure problems into large systems. This framework opens a promising path for building tunable, systematically improvable methods to capture properties of systems interacting with a large number of other systems. The extension to excited states is also straightforward.
We present a detailed discussion of our novel diagrammatic coupled cluster Monte Carlo (diagCCMC) [Scott et al. J. Phys. Chem. Lett. 2019, 10, 925]. The diagCCMC algorithm performs an imaginary-time propagation of the similarity-transformed coupled cluster Schrodinger equation. Imaginary-time updates are computed by stochastic sampling of the coupled cluster vector function: each term is evaluated as a randomly realised 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 realisation 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, independently of truncation level in coupled cluster theory. The algorithm only requires an asymptotic memory costs scaling linearly, as demonstrated previously. These scaling reductions require no ad hoc modifications to the approach.
We present a coupled cluster and linear response theory to compute properties of many-electron systems at non-zero temperatures. For this purpose, we make use of the thermofield dynamics, which allows for a compact wavefunction representation of the thermal density matrix, and extend our recently developed framework [J. Chem. Phys. 150, 154109 (2019)] to parameterize the so-called thermal state using an exponential ansatz with cluster operators that create thermal quasiparticle excitations on a mean-field reference. As benchmark examples, we apply this method to both model (one-dimensional Hubbard and Pairing) as well as ab-initio (atomic Beryllium and molecular Hydrogen) systems, while comparing with exact results.
A block-correlated coupled cluster (BCCC) method based on the generalized valence bond (GVB) wave function (GVB-BCCC in short) is proposed and implemented at the ab initio level, which represents an attractive multireference electronic structure method for strongly correlated systems. The GVB-BCCC method is demonstrated to provide accurate descriptions for multiple bond breaking in small molecules, although the GVB reference function is qualitatively wrong for the studied processes. For a challenging prototype of strongly correlated systems, tridecane with all 12 single C-C bonds at various distances, our calculations have shown that the GVB-BCCC2b method can provide highly comparable results as the density matrix renormalization group method for potential energy surfaces along simultaneous dissociation of all C-C bonds.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
