Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userComplex Ground-State and Excitation Energies in Coupled-Cluster Theory
65
0
0.0
(
0
)
Authors
Simon Thomas
Ask ChatGPT about the research
No Arabic abstract
Since in coupled-cluster (CC) theory ground-state and excitation energies are eigenvalues of a non-Hermitian matrix, these energies can in principle take on complex values. In this paper we discuss the appearance of complex energy values in CC calculations from a mathematical perspective. We analyze the behaviour of the eigenvalues of Hermitian matrices that are perturbed (in a non-Hermitian manner) by a real parameter. Based on these results we show that for CC calculations with real-valued Hamiltonian matrices the ground-state energy generally takes a real value. Furthermore, we show that in the case of real-valued Hamiltonian matrices complex excitation energies only occur in the context of conical intersections. In such a case, unphysical consequences are encountered such as a wrong dimension of the intersection seam, large numerical deviations from full configuration-interaction (FCI) results, and the square-root-like behaviour of the potential surfaces near the conical intersection. In the case of CC calculations with complex-valued Hamiltonian matrix elements, it turns out that complex energy values are to be expected for ground and excited states when no symmetry is present. We confirm the occurrence of complex energies by sample calculations using a six-state model and by CC calculations for the H2O molecule in a strong magnetic field. We furthermore show that symmetry can prevent the occurrence of complex energy values. Lastly, we demonstrate that in most cases the real part of the complex energy values provides a very good approximation to the FCI energy.
rate research
Read More
The generator-coordinate method is a flexible and powerful reformulation of the variational principle. Here we show that by introducing a generator coordinate in the Kohn-Sham equation of density-functional theory, excitation energies can be obtained from ground-state density functionals. As a viability test, the method is applied to ground-state energies and various types of excited-state energies of atoms and ions from the He and the Li isoelectronic series. Results are compared to a variety of alternative DFT-based approaches to excited states, in particular time-dependent density-functional theory with exact and approximate potentials.
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.
The time-dependent density functional theory (TDDFT) has been broadly used to investigate the excited-state properties of various molecular systems. However, the current TDDFT heavily relies on outcomes from the corresponding ground-state density functional theory (DFT) calculations which may be prone to errors due to the lack of proper treatment in the non-dynamical correlation effects. Recently, thermally-assisted-occupation density functional theory (TAO-DFT) [J.-D. Chai, textit{J. Chem. Phys.} textbf{136}, 154104 (2012)], a DFT with fractional orbital occupations, was proposed, explicitly incorporating the non-dynamical correlation effects in the ground-state calculations with low computational complexity. In this work, we develop time-dependent (TD) TAO-DFT, which is a time-dependent, linear-response theory for excited states within the framework of TAO-DFT. With tests on the excited states of H$_{2}$, the first triplet excited state ($1^3Sigma_u^+$) was described well, with non-imaginary excitation energies. TDTAO-DFT also yields zero singlet-triplet gap in the dissociation limit, for the ground singlet ($1^1Sigma_g^+$) and the first triplet state ($1^3Sigma_u^+$). In addition, as compared to traditional TDDFT, the overall excited-state potential energy surfaces obtained from TDTAO-DFT are generally improved and better agree with results from the equation-of-motion coupled-cluster singles and doubles (EOM-CCSD).
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%.
We present detailed spectral calculations for small Lieb lattices having up to $N=4$ number of cells, in the regime of half-filling, an instance of particular relevance for the nano-magnetism of discrete systems such as quantum dot arrays, due to the degenerate levels at mid-spectrum. While for the Hubbard interaction model -and even number of sites- the ground state spin is given by the Lieb theorem, the inclusion of long range interaction -or odd number of sites- make the spin state not a priori known, which justifies our approach. We calculate also the excitation energies, which are of experimental importance, and find significant variation induced by the interaction potential. One obtains insights on the mechanisms involved that impose as ground state the Lieb state with lower spin rather than the Hund one with maximum spin for the degenerate levels, showing this in the first and second order of the interaction potential for the smaller lattices. The analytical results concorde with the numerical ones, which are performed by exact diagonalization calculations or by a combined mean-field and configuration interaction method. While the Lieb state is always lower in energy than the Hund state, for strong long-range interaction, when possible, another minimal spin state is imposed as ground state.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
