-
Real-Time Coupled Cluster Theory with Approximate Triples
Authors:
Zhe Wang,
Håkon Emil Kristiansen,
Thomas Bondo Pedersen,
T. Daniel Crawford
Abstract:
The formalism of real-time (RT) methods has been well-established during recent years, while no inclusion beyond the double excitation has been discussed. In this article, we introduce an implementation of real-time coupled cluster singles, doubles and approximate triples (CC3) method to explore the potential of a high excitation level. The CC3 method is well-known for its advantages in calculatin…
▽ More
The formalism of real-time (RT) methods has been well-established during recent years, while no inclusion beyond the double excitation has been discussed. In this article, we introduce an implementation of real-time coupled cluster singles, doubles and approximate triples (CC3) method to explore the potential of a high excitation level. The CC3 method is well-known for its advantages in calculating dynamic properties and combining with the response theory. It is a well-qualified candidate for handling the interaction between the system and the applied field, and therefore suitable for a RT implementation. The derivation and implementation are first demonstrated following applications on calculating frequency-dependent properties. Terms with triples are calculated and added upon the existing CCSD equations, giving the method a formally $N^{7}$ scaling. The Graphics Processing Unit (GPU) accelerated implementation is utilized to reduce the computational cost. It is been verified that the GPU implementation can speed up the calculation by up to a factor of 17 for water cluster test cases. Additionally, the single-precision arithmetic is used and compared to the conventional double-precision arithmetic. No significant difference is found in the polarizabilities and $G'$ tensor results, but a higher percentage error for the first hyperpolarizabilities is observed. Compared to the linear response (LR) CC3 results, the percentage errors of RT-CC3 polarizabilities and RT-CC3 first hyperpolarizabilities are under 0.1% and 1%, respectively for the $H_2O$/cc-pVDZ test case. Furthermore, a discussion on the calculation of polarizabilities is included, which compares RT-CC3 with RT-CCSD and time-dependent nonorthogonal orbital-optimized coupled cluster doubles (TDNOCCD), in order to examine the performance of RT-CC3 and the orbital-optimization effect using a group of ten-electron systems.
△ Less
Submitted 11 July, 2024;
originally announced July 2024.
-
Time evolution as an optimization problem: The hydrogen atom in strong laser fields in a basis of time-dependent Gaussian wave packets
Authors:
Simon Elias Schrader,
Håkon Emil Kristiansen,
Thomas Bondo Pedersen,
Simen Kvaal
Abstract:
Recent advances in attosecond science have made it increasingly important to develop stable, reliable and accurate algorithms and methods to model the time evolution of atoms and molecules in intense laser fields. A key process in attosecond science is high-harmonic generation, which is challenging to model with fixed Gaussian basis sets, as it produces high-energy electrons, with a resulting rapi…
▽ More
Recent advances in attosecond science have made it increasingly important to develop stable, reliable and accurate algorithms and methods to model the time evolution of atoms and molecules in intense laser fields. A key process in attosecond science is high-harmonic generation, which is challenging to model with fixed Gaussian basis sets, as it produces high-energy electrons, with a resulting rapidly varying and highly oscillatory wave function that extends over dozens of ångström. Recently, Rothe's method, where time evolution is rephrased as an optimization problem, has been applied to the one-dimensional Schrödinger equation. Here, we apply Rothe's method to the hydrogen wave function and demonstrate that complex-valued Gaussian wave packets with time-dependent width, center, and momentum parameters are able to reproduce spectra obtained from essentially exact grid calculations for high-harmonic generation with only 50-181 Gaussians for field strengths up to $5\times 10^{14}$W/cm$^2$. This paves the way for the inclusion of continuum contributions into real-time, time-dependent electronic-structure theory with Gaussian basis sets for strong fields, and eventually accurate simulations of the time evolution of molecules without the Born-Oppenheimer approximation.
△ Less
Submitted 11 April, 2024;
originally announced April 2024.
-
Coulomb interaction-driven entanglement of electrons on helium
Authors:
Niyaz R. Beysengulov,
Johannes Pollanen,
Øyvind S. Schøyen,
Stian D. Bilek,
Jonas B. Flaten,
Oskar Leinonen,
Håkon Emil Kristiansen,
Zachary J. Stewart,
Jared D. Weidman,
Angela K. Wilson,
Morten Hjorth-Jensen
Abstract:
The generation and evolution of entanglement in quantum many-body systems is an active area of research that spans multiple fields, from quantum information science to the simulation of quantum many-body systems encountered in condensed matter, subatomic physics, and quantum chemistry. Motivated by recent experiments exploring quantum information processing systems with electrons trapped above the…
▽ More
The generation and evolution of entanglement in quantum many-body systems is an active area of research that spans multiple fields, from quantum information science to the simulation of quantum many-body systems encountered in condensed matter, subatomic physics, and quantum chemistry. Motivated by recent experiments exploring quantum information processing systems with electrons trapped above the surface of cryogenic noble gas substrates, we theoretically investigate the generation of \emph{motional} entanglement between two electrons via their unscreened Coulomb interaction. The model system consists of two electrons confined in separate electrostatic traps which establish microwave frequency quantized states of their motion. We compute the motional energy spectra of the electrons, as well as their entanglement, by diagonalizing the model Hamiltonian with respect to a single-particle Hartree product basis. This computational procedure can in turn be employed for device design and guidance of experimental implementations. In particular, the theoretical tools developed here can be used for fine tuning and optimization of control parameters in future experiments with electrons trapped above the surface of superfluid helium or solid neon.
△ Less
Submitted 22 January, 2024; v1 submitted 7 October, 2023;
originally announced October 2023.
-
Magnetic Optical Rotation from Real-Time Simulations in Finite Magnetic Fields
Authors:
Benedicte Sverdrup Ofstad,
Meilani Wibowo-Teale,
Håkon Emil Kristiansen,
Einar Aurbakken,
Marios Petros Kitsaras,
Øyvind Sigmundson Schøyen,
Eirill Hauge,
Simen Kvaal,
Stella Stopkowicz,
Andrew M. Wibowo-Teale,
Thomas Bondo Pedersen
Abstract:
We present a numerical approach to magnetic optical rotation based on real-time time-dependent electronic-structure theory. Not relying on perturbation expansions in the magnetic-field strength, the formulation allows us to test the range of validity of the linear relation between the rotation angle per unit path length and the magnetic-field strength that was established empirically by Verdet 160…
▽ More
We present a numerical approach to magnetic optical rotation based on real-time time-dependent electronic-structure theory. Not relying on perturbation expansions in the magnetic-field strength, the formulation allows us to test the range of validity of the linear relation between the rotation angle per unit path length and the magnetic-field strength that was established empirically by Verdet 160 years ago. Results obtained from time-dependent coupled-cluster and time-dependent current density-functional theory are presented for the closed-shell molecules H2, HF, and CO in magnetic fields up to 55 kT at standard temperature and pressure conditions. We find that Verdet's linearity remains valid up to roughly 10-20 kT, above which significant deviations from linearity are observed. Among the three current density-functional approximations tested in this work, the current-dependent Tao-Perdew-Staroverov-Scuseria hybrid functional performs the best in comparison with time-dependent coupled-cluster singles and doubles results for the magnetic optical rotation.
△ Less
Submitted 11 August, 2023;
originally announced August 2023.
-
Transient spectroscopy from time-dependent electronic-structure theory without multipole expansions
Authors:
Einar Aurbakken,
Benedicte Sverdrup Ofstad,
Håkon Emil Kristiansen,
Øyvind Sigmundson Schøyen,
Simen Kvaal,
Lasse Kragh Sørensen,
Roland Lindh,
Thomas Bondo Pedersen
Abstract:
Based on the work done by an electromagnetic field on an atomic or molecular electronic system, a general gauge invariant formulation of transient absorption spectroscopy is presented within the semi-classical approximation. Avoiding multipole expansions, a computationally viable expression for the spectral response function is derived from the minimal-coupling Hamiltonian of an electronic system…
▽ More
Based on the work done by an electromagnetic field on an atomic or molecular electronic system, a general gauge invariant formulation of transient absorption spectroscopy is presented within the semi-classical approximation. Avoiding multipole expansions, a computationally viable expression for the spectral response function is derived from the minimal-coupling Hamiltonian of an electronic system interacting with one or more laser pulses described by a source-free, enveloped electromagnetic vector potential. With a fixed-basis expansion of the electronic wave function, the computational cost of simulations of laser-driven electron dynamics beyond the dipole approximation is the same as simulations adopting the dipole approximation. We illustrate the theory by time-dependent configuration interaction and coupled-cluster simulations of core-level absorption and circular dichroism spectra.
△ Less
Submitted 5 July, 2023;
originally announced July 2023.
-
Cost-Efficient High-Resolution Linear Absorption Spectra Through Extrapolating the Dipole Moment from Real-Time Time-Dependent Electronic-Structure Theory
Authors:
Eirill Hauge,
Håkon Emil Kristiansen,
Lukas Konecny,
Marius Kadek,
Michal Repisky,
Thomas Bondo Pedersen
Abstract:
We present a novel function fitting method for approximating the propagation of the time-dependent electric dipole moment from real-time electronic structure calculations. Real-time calculations of the electronic absorption spectrum require discrete Fourier transforms of the electric dipole moment. The spectral resolution is determined by the total propagation time, i.e. the trajectory length of t…
▽ More
We present a novel function fitting method for approximating the propagation of the time-dependent electric dipole moment from real-time electronic structure calculations. Real-time calculations of the electronic absorption spectrum require discrete Fourier transforms of the electric dipole moment. The spectral resolution is determined by the total propagation time, i.e. the trajectory length of the dipole moment, causing a high computational cost. Our developed method uses function fitting on shorter trajectories of the dipole moment, achieving arbitrary spectral resolution through extrapolation. Numerical testing shows that the fitting method can reproduce high-resolution spectra using short dipole trajectories. The method converges with as little as 100 a.u. dipole trajectories for some systems, though the difficulty converging increases with the spectral density. We also introduce an error estimate of the fit, reliably assessing its convergence and hence the quality of the approximated spectrum.
△ Less
Submitted 31 October, 2023; v1 submitted 4 July, 2023;
originally announced July 2023.
-
Adiabatic extraction of nonlinear optical properties from real-time time-dependent electronic-structure theory
Authors:
Benedicte Sverdrup Ofstad,
Håkon Emil Kristiansen,
Einar Aurbakken,
Øyvind Sigmundson Schøyen,
Simen Kvaal,
Thomas Bondo Pedersen
Abstract:
Real-time simulations of laser-driven electron dynamics contain information about molecular optical properties through all orders in response theory. These properties can be extracted by assuming convergence of the power series expansion of induced electric and magnetic multipole moments. However, the accuracy relative to analytical results from response theory quickly deteriorates for higher-orde…
▽ More
Real-time simulations of laser-driven electron dynamics contain information about molecular optical properties through all orders in response theory. These properties can be extracted by assuming convergence of the power series expansion of induced electric and magnetic multipole moments. However, the accuracy relative to analytical results from response theory quickly deteriorates for higher-order responses due to the presence of high-frequency oscillations in the induced multipole moment in the time domain. This problem has been ascribed to missing higher-order corrections. We here demonstrate that the deviations are caused by nonadiabatic effects arising from the finite-time ramping from zero to full strength of the external laser field. Three different approaches, two using a ramped wave and one using a pulsed wave, for extracting electrical properties from real-time time-dependent electronic-structure simulations are investigated. The standard linear ramp is compared to a quadratic ramp, which is found to yield highly accurate results for polarizabilities, and first and second hyperpolarizabilities, at roughly half the computational cost. Results for the third hyperpolarizability are presented along with a simple, computable measure of reliability.
△ Less
Submitted 6 February, 2023;
originally announced February 2023.
-
The $S$-diagnostic -- an a posteriori error assessment for single-reference coupled-cluster methods
Authors:
Fabian M. Faulstich,
Håkon E. Kristiansen,
Mihaly A. Csirik,
Simen Kvaal,
Thomas Bondo Pedersen,
Andre Laestadius
Abstract:
We propose a novel a posteriori error assessment for the single-reference coupled-cluster (SRCC) method called the $S$-diagnostic. We provide a derivation of the $S$-diagnostic that is rooted in the mathematical analysis of different SRCC variants. We numerically scrutinized the $S$-diagnostic, testing its performance for (1) geometry optimizations, (2) electronic correlation simulations of system…
▽ More
We propose a novel a posteriori error assessment for the single-reference coupled-cluster (SRCC) method called the $S$-diagnostic. We provide a derivation of the $S$-diagnostic that is rooted in the mathematical analysis of different SRCC variants. We numerically scrutinized the $S$-diagnostic, testing its performance for (1) geometry optimizations, (2) electronic correlation simulations of systems with varying numerical difficulty, and (3) the square-planar copper complexes [CuCl$_4$]$^{2-}$, [Cu(NH$_3$)$_4$]$^{2+}$, and [Cu(H$_2$O)$_4$]$^{2+}$. Throughout the numerical investigations, the $S$-diagnostic is compared to other SRCC diagnostic procedures, that is, the $T_1$, $D_1$, and $D_2$ diagnostics as well as different indices of multi-determinantal and multi-reference character in coupled-cluster theory. Our numerical investigations show that the $S$-diagnostic outperforms the $T_1$, $D_1$, and $D_2$ diagnostics and is comparable to the indices of multi-determinantal and multi-reference character in coupled-cluster theory in their individual fields of applicability. The experiments investigating the performance of the $S$-diagnostic for geometry optimizations using SRCC reveal that the $S$-diagnostic correlates well with different error measures at a high level of statistical relevance. The experiments investigating the performance of the $S$-diagnostic for electronic correlation simulations show that the $S$-diagnostic correctly predicts strong multi-reference regimes. The $S$-diagnostic moreover correctly detects the successful SRCC computations for [CuCl$_4$]$^{2-}$, [Cu(NH$_3$)$_4$]$^{2+}$, and [Cu(H$_2$O)$_4$]$^{2+}$, which have been known to be misdiagnosed by $T_1$ and $D_1$ diagnostics in the past. This shows that the $S$-diagnostic is a promising candidate for an a posteriori diagnostic for SRCC calculations.
△ Less
Submitted 26 January, 2023;
originally announced January 2023.
-
Linear and nonlinear optical properties from TDOMP2 theory
Authors:
Håkon Emil Kristiansen,
Benedicte Sverdrup Ofstad,
Eirill Hauge,
Einar Aurbakken,
Øyvind Sigmundson Schøyen,
Simen Kvaal,
Thomas Bondo Pedersen
Abstract:
In this work we present a derivation of the real-time time-dependent orbital-optimized Møller-Plesser TDOMP2 and its biorthogonal companion, time-dependent non-orthogonal OMP2 (TDNOMP2), theory starting from the time-dependent bivariational principle and a parametrization based on the exponential orbital-rotation operator formulation commonly used in time-independent molecular electronic-structure…
▽ More
In this work we present a derivation of the real-time time-dependent orbital-optimized Møller-Plesser TDOMP2 and its biorthogonal companion, time-dependent non-orthogonal OMP2 (TDNOMP2), theory starting from the time-dependent bivariational principle and a parametrization based on the exponential orbital-rotation operator formulation commonly used in time-independent molecular electronic-structure theory. We apply the TDOMP2 method to extract absorption spectra and frequency-dependent polarizabilities and first hyperpolarizabilities from real-time simulations, comparing the results with those obtained from conventional time-dependent coupled-cluster singles and doubles (TDCCSD) simulations and from its second-order approximation TDCC2. We also compare with results from CCSD and CC2 linear and quadratic response theory. Our results indicate that while TDOMP2 absorption spectra are of the same quality as TDCC2 spectra, frequency-dependent polarizabilities and hyperpolarizabilties from TDOMP2 simulations are significantly closer to TDCCSD results than those from TDCC2 simulations.
△ Less
Submitted 21 April, 2022; v1 submitted 27 December, 2021;
originally announced December 2021.
-
Evaluation of Gaussian integrals for the modeling of two-dimensional quantum systems
Authors:
Øyvind Sigmundson Schøyen,
Håkon Emil Kristiansen,
Alfred Alocias Mariadason
Abstract:
We have developed a McMurchie-Davidson-like recursion formula for efficient evaluation of the Coulomb attraction and interaction matrix elements between two-dimensional primitive Cartesian Gaussian type orbitals. We also present recurrence schemes for combined position and differential operator integrals, and three-center Gaussian integrals. The Cartesian Gaussian orbitals are isotropic in the exp…
▽ More
We have developed a McMurchie-Davidson-like recursion formula for efficient evaluation of the Coulomb attraction and interaction matrix elements between two-dimensional primitive Cartesian Gaussian type orbitals. We also present recurrence schemes for combined position and differential operator integrals, and three-center Gaussian integrals. The Cartesian Gaussian orbitals are isotropic in the exponent, but with arbitrary centers and angular momentum.
△ Less
Submitted 29 November, 2021; v1 submitted 18 November, 2021;
originally announced November 2021.
-
Interpretation of Coupled-Cluster Many-Electron Dynamics in Terms of Stationary States
Authors:
Thomas Bondo Pedersen,
Håkon Emil Kristiansen,
Tilmann Bodenstein,
Simen Kvaal,
Øyvind Sigmundson Schøyen
Abstract:
We demonstrate theoretically and numerically that laser-driven many-electron dynamics, as described by bivariational time-dependent coupled-cluster theory, may be analyzed in terms of stationary-state populations. Projectors heuristically defined from linear response theory and equation-of-motion coupled-cluster theory are proposed for the calculation of stationary-state populations during interac…
▽ More
We demonstrate theoretically and numerically that laser-driven many-electron dynamics, as described by bivariational time-dependent coupled-cluster theory, may be analyzed in terms of stationary-state populations. Projectors heuristically defined from linear response theory and equation-of-motion coupled-cluster theory are proposed for the calculation of stationary-state populations during interaction with laser pulses or other external forces, and conservation laws of the populations are discussed. Numerical tests of the proposed projectors, involving both linear and nonlinear optical processes for the He and Be atoms, and for the LiH, CH$^+$, and LiF molecules, show that the laser-driven evolution of the stationary-state populations at the coupled-cluster singles-and-doubles (CCSD) level is very close to that obtained by full configuration-interaction theory provided all stationary states actively participating in the dynamics are sufficiently well approximated. When double-excited states are important for the dynamics, the quality of the CCSD results deteriorate. Observing that populations computed from the linear-response projector may show spurious small-amplitude, high-frequency oscillations, the equation-of-motion projector emerges as the most promising approach to stationary-state populations.
△ Less
Submitted 19 December, 2020; v1 submitted 21 September, 2020;
originally announced September 2020.
-
Numerical stability of time-dependent coupled-cluster methods for many-electron dynamics in intense laser pulses
Authors:
Håkon Emil Kristiansen,
Øyvind Sigmundson Schøyen,
Simen Kvaal,
Thomas Bondo Pedersen
Abstract:
We investigate the numerical stability of time-dependent coupled-cluster theory for many-electron dynamics in intense laser pulses, comparing two coupled-cluster formulations with full configuration interaction theory. Our numerical experiments show that orbital-adaptive time-dependent coupled-cluster doubles (OATDCCD) theory offers significantly improved stability compared with the conventional H…
▽ More
We investigate the numerical stability of time-dependent coupled-cluster theory for many-electron dynamics in intense laser pulses, comparing two coupled-cluster formulations with full configuration interaction theory. Our numerical experiments show that orbital-adaptive time-dependent coupled-cluster doubles (OATDCCD) theory offers significantly improved stability compared with the conventional Hartree-Fock-based time-dependent coupled-cluster singles-and-doubles (TDCCSD) formulation. The improved stability stems from greatly reduced oscillations in the doubles amplitudes, which, in turn, can be traced to the dynamic biorthonormal reference determinants of OATDCCD theory. As long as these are good approximations to the Brueckner determinant, OATDCCD theory is numerically stable. We propose the reference weight as a diagnostic quantity to identify situations where the TDCCSD and OATDCCD theories become unstable.
△ Less
Submitted 24 February, 2020; v1 submitted 19 December, 2019;
originally announced December 2019.