Showing 1–42 of 42 results for author: Kvaal, S

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.

Comments: 20 pages, 10 figures

Gaussians for Electronic and Rovibrational Quantum Dynamics

Authors: Aleksander P. Woźniak, Ludwik Adamowicz, Thomas Bondo Pedersen, Simen Kvaal

Abstract: The assumptions underpinning the adiabatic Born-Oppenheimer (BO) approximation are broken for molecules interacting with attosecond laser pulses, which generate complicated coupled electronic-nuclear wavepackets that generally will have components of electronic and dissociation continua as well as bound-state contributions. The conceptually most straightforward way to overcome this challenge is to… ▽ More The assumptions underpinning the adiabatic Born-Oppenheimer (BO) approximation are broken for molecules interacting with attosecond laser pulses, which generate complicated coupled electronic-nuclear wavepackets that generally will have components of electronic and dissociation continua as well as bound-state contributions. The conceptually most straightforward way to overcome this challenge is to treat the electronic and nuclear degrees of freedom on equal quantum-mechanical footing by not invoking the BO approximation at all. Explicitly correlated Gaussian (ECG) basis functions have proved successful for non-BO calculations of stationary molecular states and energies, reproducing rovibrational absorption spectra with very high accuracy. In this paper, we present a proof-of-principle study of the ability of fully flexible ECGs (FFECGs) to capture the intricate electronic and rovibrational dynamics generated by short, high-intensity laser pulses. By fitting linear combinations of FFECGs to accurate wave function histories obtained on a large real-space grid for a regularized 2D model of the hydrogen atom and for the 2D Morse potential we demonstrate that FFECGs provide a very compact description of laser-driven electronic and rovibrational dynamics. △ Less

Submitted 12 April, 2024; v1 submitted 22 January, 2024; originally announced January 2024.

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.

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.

Accelerated coupled cluster calculations with Procrustes orbital interpolation

Authors: Simon Elias Schrader, Simen Kvaal

Abstract: The coupled cluster method is considered a gold standard in quantum chemistry, reliably giving energies that are exact within chemical accuracy (1.6 mHartree). However, even in the CCSD approximation, where the cluster operator is truncated to include only single and double excitations, the method scales as $O(N^6)$ in the number of electrons, and the cluster operator needs to be solved for iterat… ▽ More The coupled cluster method is considered a gold standard in quantum chemistry, reliably giving energies that are exact within chemical accuracy (1.6 mHartree). However, even in the CCSD approximation, where the cluster operator is truncated to include only single and double excitations, the method scales as $O(N^6)$ in the number of electrons, and the cluster operator needs to be solved for iteratively, increasing computation time. Inspired by eigenvector continuation, we present here an algorithm making use of Gaussian processes that provides an improved initial guess for the coupled cluster amplitudes. The cluster operator is written as a linear combination of sample cluster operators which are obtained at particular sample geometries. By reusing the cluster operators from previous calculations in that way, it is possible to obtain a start guess for the amplitudes that surpasses both MP2-guesses and "previous geometry"-guesses in terms of the number of necessary iterations. As this improved guess is very close to the exact cluster operator, it can be used directly to calculate the CCSD energy to chemical accuracy, giving approximate CCSD energies scaling as $O(N^5)$. △ Less

Submitted 9 March, 2023; v1 submitted 15 December, 2022; originally announced December 2022.

Comments: 12 pages, 7 figures

Moreau--Yosida regularization in DFT

Authors: Simen Kvaal

Abstract: Moreau-Yosida regularization is introduced into the framework of exact DFT. Moreau-Yosida regularization is a lossless operation on lower semicontinuous proper convex functions over separable Hilbert spaces, and when applied to the universal functional of exact DFT (appropriately restricted to a bounded domain), gives a reformulation of the ubiquitous $v$-representability problem and a rigorous an… ▽ More Moreau-Yosida regularization is introduced into the framework of exact DFT. Moreau-Yosida regularization is a lossless operation on lower semicontinuous proper convex functions over separable Hilbert spaces, and when applied to the universal functional of exact DFT (appropriately restricted to a bounded domain), gives a reformulation of the ubiquitous $v$-representability problem and a rigorous and illuminating derivation of Kohn-Sham theory. The chapter comprises a self-contained introduction to exact DFT, basic tools from convex analysis such as sub- and superdifferentiability and convex conjugation, as well as basic results on the Moreau-Yosida regularization. The regularization is then applied to exact DFT and Kohn-Sham theory, and a basic iteration scheme based in the Optimal Damping Algorithm is analyzed. In particular, its global convergence established. Some perspectives are offered near the end of the chapter. △ Less

Submitted 10 August, 2022; originally announced August 2022.

Comments: To be published in a contributed Springer volume

No need for a grid: Adaptive fully-flexible gaussians for the time-dependent Schrödinger equation

Authors: Simen Kvaal, Caroline Lasser, Thomas Bondo Pedersen, Ludwik Adamowicz

Abstract: Linear combinations of complex gaussian functions, where the linear and nonlinear parameters are allowed to vary, are shown to provide an extremely flexible and effective approach for solving the time-dependent Schrödinger equation in one spatial dimension. The use of flexible basis sets has been proven notoriously hard within the systematics of the Dirac--Frenkel variational principle. In this wo… ▽ More Linear combinations of complex gaussian functions, where the linear and nonlinear parameters are allowed to vary, are shown to provide an extremely flexible and effective approach for solving the time-dependent Schrödinger equation in one spatial dimension. The use of flexible basis sets has been proven notoriously hard within the systematics of the Dirac--Frenkel variational principle. In this work we present an alternative time-propagation scheme that de-emphasizes optimal parameter evolution but directly targets residual minimization via the method of Rothe's method, also called the method of vertical time layers. We test the scheme using a simple model system mimicking an atom subjected to an extreme laser pulse. Such a pulse produces complicated ionization dynamics of the system. The scheme is shown to perform very well on this model and notably does not rely on a computational grid. Only a handful of gaussian functions are needed to achieve an accuracy on par with a high-resolution, grid-based solver. This paves the way for accurate and affordable solution of the time-dependent Schrödinger equation for atoms and molecules within and beyond the Born--Oppenheimer approximation. △ Less

Submitted 7 March, 2023; v1 submitted 1 July, 2022; originally announced July 2022.

Comments: 8 pages, 6 figures

Laser-induced dynamic alignment of the HD molecule without the Born-Oppenheimer approximation

Authors: Ludwik Adamowicz, Simen Kvaal, Caroline Lasser, Thomas Bondo Pedersen

Abstract: Laser-induced molecular alignment is well understood within the framework of the Born-Oppenheimer (BO) approximation Without the BO approximation, however, the concept of molecular structure is lost, making alignment hard to define precisely. In this work, we demonstrate the emergence of alignment from the first-ever non-BO quantum dynamics simulations, using the HD molecule exposed to ultrashort… ▽ More Laser-induced molecular alignment is well understood within the framework of the Born-Oppenheimer (BO) approximation Without the BO approximation, however, the concept of molecular structure is lost, making alignment hard to define precisely. In this work, we demonstrate the emergence of alignment from the first-ever non-BO quantum dynamics simulations, using the HD molecule exposed to ultrashort laser pulses as a few-body test case We extract the degree of alignment from the non-BO wave function by means of an operator expressed in terms of pseudo-proton coordinates that mimics the BO-based definition of alignment The only essential approximation, in addition to the semiclassical electric-dipole approximation for the matter-field interaction, is the choice of time-independent explicitly correlated Gaussian basis functions. We use a variational, electric-field-dependent basis-set construction procedure, which allows us to keep the basis-set dimension low whilst capturing the main effects of electric polarization on the nuclear and electronic degrees of freedom. The basis-set construction procedure is validated by comparing with virtually exact grid-based simulations for two one-dimensional model systems: laser-driven electron dynamics in a soft attractive Coulomb potential and nuclear rovibrational dynamics in a Morse potential. △ Less

Submitted 15 September, 2022; v1 submitted 30 May, 2022; originally announced May 2022.

Comments: 8 pages, 6 figures

Three Lagrangians for the complete-active space coupled-cluster method

Authors: Simen Kvaal

Abstract: Three fully variational formulations of the complete-active space coupled-cluster (CASCC) method are derived. The formulations include the ability to approximate the model vectors by smooth manifolds, thereby opening up the possibility for overcoming the exponential wall of scaling for model spaces of CAS type. In particular, model vectors of matrix-product states are considered, and it is argued… ▽ More Three fully variational formulations of the complete-active space coupled-cluster (CASCC) method are derived. The formulations include the ability to approximate the model vectors by smooth manifolds, thereby opening up the possibility for overcoming the exponential wall of scaling for model spaces of CAS type. In particular, model vectors of matrix-product states are considered, and it is argued that the present variational formulation allows not only favorably-scaling multireference coupled-cluster calculations, but also systematic correction of tailored coupled-cluster calculation and of quantum chemical density-matrix renormalization group methods, which are fast and polynomial scaling, but lacks the ability to properly resolve dynamical correlation at chemical accuracy. The extension of the variational formulations to the time-domain is also discussed, with derivations of abstract evolution equations. △ Less

Submitted 10 June, 2023; v1 submitted 18 May, 2022; originally announced May 2022.

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.

Lower semi-continuity of universal functional in paramagnetic current-density functional theory

Authors: Simen Kvaal, Andre Laestadius, Erik I. Tellgren, Trygve U. Helgaker

Abstract: A cornerstone of current-density functional theory (CDFT) in its paramagnetic formulation is proven. After a brief outline of the mathematical structure of CDFT, the lower semi-continuity and expectation valuedness of the CDFT constrained-search functional is proven, meaning that there is always a minimizing density matrix in the CDFT constrained-search universal density functional. These results… ▽ More A cornerstone of current-density functional theory (CDFT) in its paramagnetic formulation is proven. After a brief outline of the mathematical structure of CDFT, the lower semi-continuity and expectation valuedness of the CDFT constrained-search functional is proven, meaning that there is always a minimizing density matrix in the CDFT constrained-search universal density functional. These results place the mathematical framework of CDFT on the same footing as that of standard DFT. △ Less

Submitted 10 November, 2020; originally announced November 2020.

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.

Comments: 58 pages, 14 figures

A state-specific multireference coupled-cluster method based on the bivariational principle

Authors: Tilmann Bodenstein, Simen Kvaal

Abstract: A state-specific multireference coupled-cluster method based on Arponen's bivariational principle is presented, the bivar-MRCC method. The method is based on singlereference theory, and therefore has a relatively straightforward formulation and modest computational complexity. The main difference from established methods is the bivariational formulation, in which independent parameterizations of t… ▽ More A state-specific multireference coupled-cluster method based on Arponen's bivariational principle is presented, the bivar-MRCC method. The method is based on singlereference theory, and therefore has a relatively straightforward formulation and modest computational complexity. The main difference from established methods is the bivariational formulation, in which independent parameterizations of the wavefunction (ket) and its complex conjugate (bra) are made. Importantly, this allows manifest multiplicative separability (exact in the extended bivar-MRECC version of the method, and approximate otherwise), while preserving polynomial scaling of the working equations. A feature of the bivariational principle is that the formal bra and ket references can be included as bivariational parameters, which eliminates much of the bias towards the formal reference. A pilot implementation is described, and extensive benchmark calculations on several standard problems are performed. The results from the bivar-MRCC method are comparable to established state-specific multireference methods. Considering the relative affordability of the bivar-MRCC method, it may become a practical tool for non-experts. △ Less

Submitted 31 March, 2020; v1 submitted 30 March, 2020; originally announced March 2020.

Guaranteed convergence for a class of coupled-cluster methods based on Arponen's extended theory

Authors: Simen Kvaal, Andre Laestadius, Tilmann Bodenstein

Abstract: A wide class of coupled-cluster methods is introduced, based on Arponen's extended coupled-cluster theory. This class of methods is formulated in terms of a coordinate transformation of the cluster operators. The mathematical framework for the error analysis of coupled-cluster methods based on Arponen's bivariational principle is presented, in which the concept of local strong monotonicity of the… ▽ More A wide class of coupled-cluster methods is introduced, based on Arponen's extended coupled-cluster theory. This class of methods is formulated in terms of a coordinate transformation of the cluster operators. The mathematical framework for the error analysis of coupled-cluster methods based on Arponen's bivariational principle is presented, in which the concept of local strong monotonicity of the flipped gradient of the energy is central. A general mathematical result is presented, describing sufficient conditions for coordinate transformations to preserve the local strong monotonicity. The result is applied to the presented class of methods, which include the standard and quadratic coupled-cluster methods, and also Arponen's canonical version of extended coupled-cluster theory. Some numerical experiments are presented, and the use of canonical coordinates for diagnostics is discussed. △ Less

Submitted 15 March, 2020; originally announced March 2020.

Comments: 11 pages

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.

Comments: 5 pages, 6 figures (supplemental material, 7 pages, 11 figures)

Journal ref: The Journal of Chemical Physics, 152, (2020), 071102

Kohn-Sham theory with paramagnetic currents: compatibility and functional differentiability

Authors: Andre Laestadius, Erik I. Tellgren, Markus Penz, Michael Ruggenthaler, Simen Kvaal, Trygve Helgaker

Abstract: Recent work has established Moreau-Yosida regularization as a mathematical tool to achieve rigorous functional differentiability in density-functional theory. In this article, we extend this tool to paramagnetic current-density-functional theory, the most common density-functional framework for magnetic field effects. The extension includes a well-defined Kohn-Sham iteration scheme with a partial… ▽ More Recent work has established Moreau-Yosida regularization as a mathematical tool to achieve rigorous functional differentiability in density-functional theory. In this article, we extend this tool to paramagnetic current-density-functional theory, the most common density-functional framework for magnetic field effects. The extension includes a well-defined Kohn-Sham iteration scheme with a partial convergence result. To this end, we rely on a formulation of Moreau-Yosida regularization for reflexive and strictly convex function spaces. The optimal $L^p$-characterization of the paramagnetic current density $L^1\cap L^{3/2}$ is derived from the $N$-representability conditions. A crucial prerequisite for the convex formulation of paramagnetic current-density-functional theory, termed compatibility between function spaces for the particle density and the current density, is pointed out and analyzed. Several results about compatible function spaces are given, including their recursive construction. The regularized, exact functionals are calculated numerically for a Kohn-Sham iteration on a quantum ring, illustrating their performance for different regularization parameters. △ Less

Submitted 21 May, 2019; v1 submitted 24 February, 2019; originally announced February 2019.

Journal ref: J. Chem. Theory Comput.2019, 15, 7, 4003-4020 (2019)

Symplectic integration and physical interpretation of time-dependent coupled-cluster theory

Authors: Thomas Bondo Pedersen, Simen Kvaal

Abstract: The formulation of the time-dependent Schrodinger equation in terms of coupled-cluster theory is outlined, with emphasis on the bivariational framework and its classical Hamiltonian structure. An indefinite inner product is introduced, inducing physical interpretation of coupled-cluster states in the form of transition probabilities, autocorrelation functions, and explicitly real values for observ… ▽ More The formulation of the time-dependent Schrodinger equation in terms of coupled-cluster theory is outlined, with emphasis on the bivariational framework and its classical Hamiltonian structure. An indefinite inner product is introduced, inducing physical interpretation of coupled-cluster states in the form of transition probabilities, autocorrelation functions, and explicitly real values for observables, solving interpretation issues which are present in time-dependent coupled-cluster theory and in ground-state calculations of molecular systems under influence of external magnetic fields. The problem of the numerical integration of the equations of motion is considered, and a critial evaluation of the standard fourth-order Runge--Kutta scheme and the symplectic Gauss integrator of variable order is given, including several illustrative numerical experiments. While the Gauss integrator is stable even for laser pulses well above the perturbation limit, our experiments indicate that a system-dependent upper limit exists for the external field strengths. Above this limit, time-dependent coupled-cluster calculations become very challenging numerically, even in the full configuration interaction limit. The source of these numerical instabilities is shown to be rapid increases of the amplitudes as ultrashort high-intensity laser pulses pump the system out of the ground state into states that are virtually orthogonal to the static Hartree-Fock reference determinant. △ Less

Submitted 14 April, 2019; v1 submitted 11 December, 2018; originally announced December 2018.

Comments: 14 pages, 13 figures

Journal ref: J. Chem. Phys. 150 (2019) 144106

Numerical and Theoretical Aspects of the DMRG-TCC Method Exemplified by the Nitrogen Dimer

Authors: Fabian M. Faulstich, Mihály Máté, Andre Laestadius, Mihály András Csirik, Libor Veis, Andrej Antalik, Jiří Brabec, Reinhold Schneider, Jiří Pittner, Simen Kvaal, Örs Legeza

Abstract: In this article, we investigate the numerical and theoretical aspects of the coupled-cluster method tailored by matrix-product states. We investigate chemical properties of the used method, such as energy size extensivity and the equivalence of linked and unlinked formulation. The existing mathematical analysis is here elaborated in a quantum chemical framework. In particular, we highlight the use… ▽ More In this article, we investigate the numerical and theoretical aspects of the coupled-cluster method tailored by matrix-product states. We investigate chemical properties of the used method, such as energy size extensivity and the equivalence of linked and unlinked formulation. The existing mathematical analysis is here elaborated in a quantum chemical framework. In particular, we highlight the use of a so-called CAS-ext gap describing the basis splitting between the complete active space and the external part. Moreover, the behavior of the energy error as a function of the optimal basis splitting is discussed. We show numerical investigations on the robustness with respect to the bond dimensions of the single orbital entropy and the mutual information, which are quantities that are used to choose the complete active space. Furthermore, we extend the mathematical analysis with a numerical study on the complete active space dependence of the error. △ Less

Submitted 20 September, 2018; originally announced September 2018.

Generalized Kohn-Sham iteration on Banach spaces

Authors: Andre Laestadius, Markus Penz, Erik I. Tellgren, Michael Ruggenthaler, Simen Kvaal, Trygve Helgaker

Abstract: A detailed account of the Kohn-Sham algorithm from quantum chemistry, formulated rigorously in the very general setting of convex analysis on Banach spaces, is given here. Starting from a Levy-Lieb-type functional, its convex and lower semi-continuous extension is regularized to obtain differentiability. This extra layer allows to rigorously introduce, in contrast to the common unregularized appro… ▽ More A detailed account of the Kohn-Sham algorithm from quantum chemistry, formulated rigorously in the very general setting of convex analysis on Banach spaces, is given here. Starting from a Levy-Lieb-type functional, its convex and lower semi-continuous extension is regularized to obtain differentiability. This extra layer allows to rigorously introduce, in contrast to the common unregularized approach, a well-defined Kohn-Sham iteration scheme. Convergence in a weak sense is then proven. This generalized formulation is applicable to a wide range of different density-functional theories and possibly even to models outside of quantum mechanics. △ Less

Submitted 17 October, 2018; v1 submitted 23 April, 2018; originally announced April 2018.

Journal ref: J. Chem. Phys. 149, 164103 (2018)

Analysis of The Tailored Coupled-Cluster Method in Quantum Chemistry

Authors: Fabian M. Faulstich, Andre Laestadius, Örs Legeza, Reinhold Schneider, Simen Kvaal

Abstract: In quantum chemistry, one of the most important challenges is the static correlation problem when solving the electronic Schrödinger equation for molecules in the Born--Oppenheimer approximation. In this article, we analyze the tailored coupled-cluster method (TCC), one particular and promising method for treating molecular electronic-structure problems with static correlation. The TCC method comb… ▽ More In quantum chemistry, one of the most important challenges is the static correlation problem when solving the electronic Schrödinger equation for molecules in the Born--Oppenheimer approximation. In this article, we analyze the tailored coupled-cluster method (TCC), one particular and promising method for treating molecular electronic-structure problems with static correlation. The TCC method combines the single-reference coupled-cluster (CC) approach with an approximate reference calculation in a subspace [complete active space (CAS)] of the considered Hilbert space that covers the static correlation. A one-particle spectral gap assumption is introduced, separating the CAS from the remaining Hilbert space. This replaces the nonexisting or nearly nonexisting gap between the highest occupied molecular orbital and the lowest unoccupied molecular orbital usually encountered in standard single-reference quantum chemistry. The analysis covers, in particular, CC methods tailored by tensor-network states (TNS-TCC methods). The problem is formulated in a nonlinear functional analysis framework, and, under certain conditions such as the aforementioned gap, local uniqueness and existence are proved using Zarantonello's lemma. From the Aubin--Nitsche-duality method, a quadratic error bound valid for TNS-TCC methods is derived, e.g., for linear-tensor-network TCC schemes using the density matrix renormalization group method. △ Less

Submitted 13 October, 2019; v1 submitted 15 February, 2018; originally announced February 2018.

Uniform magnetic fields in density-functional theory

Authors: Erik I. Tellgren, Andre Laestadius, Trygve Helgaker, Simen Kvaal, Andrew M. Teale

Abstract: We construct a density-functional formalism adapted to uniform external magnetic fields that is intermediate between conventional Density Functional Theory and Current-Density Functional Theory (CDFT). In the intermediate theory, which we term LDFT, the basic variables are the density, the canonical momentum, and the paramagnetic contribution to the magnetic moment. Both a constrained-search formu… ▽ More We construct a density-functional formalism adapted to uniform external magnetic fields that is intermediate between conventional Density Functional Theory and Current-Density Functional Theory (CDFT). In the intermediate theory, which we term LDFT, the basic variables are the density, the canonical momentum, and the paramagnetic contribution to the magnetic moment. Both a constrained-search formulation and a convex formulation in terms of Legendre--Fenchel transformations are constructed. Many theoretical issues in CDFT find simplified analogues in LDFT. We prove results concerning $N$-representability, Hohenberg--Kohn-like mappings, existence of minimizers in the constrained-search expression, and a restricted analogue to gauge invariance. The issue of additivity of the energy over non-interacting subsystems, which is qualitatively different in LDFT and CDFT, is also discussed. △ Less

Submitted 29 September, 2017; originally announced September 2017.

Analysis of the Extended Coupled-Cluster Method in Quantum Chemistry

Authors: Andre Laestadius, Simen Kvaal

Abstract: The mathematical foundation of the so-called extended coupled-cluster method for the solution of the many-fermion Schrödinger equation is here developed. We prove an existence and uniqueness result, both in the full infinite-dimensional amplitude space as well as for discretized versions of it. The extended coupled-cluster method is formulated as a critical point of an energy function using a gene… ▽ More The mathematical foundation of the so-called extended coupled-cluster method for the solution of the many-fermion Schrödinger equation is here developed. We prove an existence and uniqueness result, both in the full infinite-dimensional amplitude space as well as for discretized versions of it. The extended coupled-cluster method is formulated as a critical point of an energy function using a generalization of the Rayleigh-Ritz principle: the bivariational principle. This gives a quadratic bound for the energy error in the discretized case. The existence and uniqueness results are proved using a type of monotonicity property for the flipped gradient of the energy function. Comparisons to the analysis of the standard coupled-cluster method is made, and it is argued that the bivariational principle is a useful tool, both for studying coupled-cluster type methods, and for developing new computational schemes in general. △ Less

Submitted 14 February, 2017; originally announced February 2017.

Ground-state densities from the Rayleigh--Ritz variation principle and from density-functional theory

Authors: Simen Kvaal, Trygve Helgaker

Abstract: The relationship between the densities of ground-state wave functions (i.e., the minimizers of the Rayleigh--Ritz (RR) variation principle) and the ground-state densities in density-functional theory (i.e., the minimizers of the Hohenberg--Kohn (HK) variation principle) is studied within the framework of convex conjugation, in a generic setting covering molecular systems, solid-state systems, and… ▽ More The relationship between the densities of ground-state wave functions (i.e., the minimizers of the Rayleigh--Ritz (RR) variation principle) and the ground-state densities in density-functional theory (i.e., the minimizers of the Hohenberg--Kohn (HK) variation principle) is studied within the framework of convex conjugation, in a generic setting covering molecular systems, solid-state systems, and more. Having introduced admissible density functionals as functionals that produce the exact ground ground-state energy for a given external potential by minimizing over densities in the HK variation principle, necessary sufficient conditions on such functionals are established to ensure that the RR ground-state densities and the HK ground-state densities are identical. We apply the results to molecular systems in the BO-approximation. For any given potential $v \in L^{3/2}(\mathbb{R}^3) + L^{\infty}(\mathbb{R}^3)$, we establish a one-to-one correspondence between the mixed ground-state densities of the RR variation principle and the mixed ground-state densities of the HK variation principle when the Lieb density-matrix constrained-search universal density functional is taken as the admissible functional. A similar one-to-one correspondence is established between the pure ground-state densities of the RR variation principle and the pure ground-state densities obtained using the HK variation principle with the Levy--Lieb pure-state constrained-search functional. In other words, all physical ground-state densities (pure or mixed) are recovered with these functionals and no false densities (i.e., minimizing densities that are not physical) exist. The importance of topology (i.e., choice of Banach space of densities and potentials) is emphasized and illustrated. The relevance of these results for current-density-functional theory is examined. △ Less

Submitted 25 June, 2015; originally announced June 2015.

Differentiable but exact formulation of density-functional theory

Authors: Simen Kvaal, Ulf Ekström, Andrew M. Teale, Trygve Helgaker

Abstract: The universal density functional $F$ of density-functional theory is a complicated and ill-behaved function of the density-in particular, $F$ is not differentiable, making many formal manipulations more complicated. Whilst $F$ has been well characterized in terms of convex analysis as forming a conjugate pair $(E,F)$ with the ground-state energy $E$ via the Hohenberg-Kohn and Lieb variation princi… ▽ More The universal density functional $F$ of density-functional theory is a complicated and ill-behaved function of the density-in particular, $F$ is not differentiable, making many formal manipulations more complicated. Whilst $F$ has been well characterized in terms of convex analysis as forming a conjugate pair $(E,F)$ with the ground-state energy $E$ via the Hohenberg-Kohn and Lieb variation principles, $F$ is nondifferentiable and subdifferentiable only on a small (but dense) set of its domain. In this article, we apply a tool from convex analysis, Moreau-Yosida regularization, to construct, for any $ε>0$, pairs of conjugate functionals $({}^ε\!E,{}^ε\!F)$ that converge to $(E,F)$ pointwise everywhere as $ε\rightarrow 0^+$, and such that ${}^ε\!F$ is (Fréchet) differentiable. For technical reasons, we limit our attention to molecular electronic systems in a finite but large box. It is noteworthy that no information is lost in the Moreau-Yosida regularization: the physical ground-state energy $E(v)$ is exactly recoverable from the regularized ground-state energy ${}^ε\!E(v)$ in a simple way. All concepts and results pertaining to the original $(E,F)$ pair have direct counterparts in results for $({}^ε\! E, {}^ε\!F)$. The Moreau-Yosida regularization therefore allows for an exact, differentiable formulation of density-functional theory. In particular, taking advantage of the differentiability of ${}^ε\!F$, a rigorous formulation of Kohn-Sham theory is presented that does not suffer from the noninteracting representability problem in standard Kohn-Sham theory. △ Less

Submitted 12 February, 2014; v1 submitted 13 December, 2013; originally announced December 2013.

Fermion $N$-representability for prescribed density and paramagnetic current density

Authors: Erik Tellgren, Simen Kvaal, Trygve Helgaker

Abstract: The $N$-representability problem is the problem of determining whether or not there exists $N$-particle states with some prescribed property. Here we report an affirmative solution to the fermion $N$-representability problem when both the density and paramagnetic current density are prescribed. This problem arises in current-density functional theory and is a generalization of the well-studied cor… ▽ More The $N$-representability problem is the problem of determining whether or not there exists $N$-particle states with some prescribed property. Here we report an affirmative solution to the fermion $N$-representability problem when both the density and paramagnetic current density are prescribed. This problem arises in current-density functional theory and is a generalization of the well-studied corresponding problem (only the density prescribed) in density functional theory. Given any density and paramagnetic current density satisfying a minimal regularity condition (essentially that a von Weizäcker-like the canonical kinetic energy density is locally integrable), we prove that there exist a corresponding $N$-particle state. We prove this by constructing an explicit one-particle reduced density matrix in the form of a position-space kernel, i.e.\ a function of two continuous position variables. In order to make minimal assumptions, we also address mathematical subtleties regarding the diagonal of, and how to rigorously extract paramagnetic current densities from, one-particle reduced density matrices in kernel form. △ Less

Submitted 28 November, 2013; v1 submitted 4 October, 2013; originally announced October 2013.

An inverse iteration method for eigenvalue problems with eigenvector nonlinearities

Authors: Elias Jarlebring, Simen Kvaal, Wim Michiels

Abstract: Consider a symmetric matrix $A(v)\in\RR^{n\times n}$ depending on a vector $v\in\RR^n$ and satisfying the property $A(αv)=A(v)$ for any $α\in\RR\backslash{0}$. We will here study the problem of finding $(λ,v)\in\RR\times \RR^n\backslash\{0\}$ such that $(λ,v)$ is an eigenpair of the matrix $A(v)$ and we propose a generalization of inverse iteration for eigenvalue problems with this type of eigenve… ▽ More Consider a symmetric matrix $A(v)\in\RR^{n\times n}$ depending on a vector $v\in\RR^n$ and satisfying the property $A(αv)=A(v)$ for any $α\in\RR\backslash{0}$. We will here study the problem of finding $(λ,v)\in\RR\times \RR^n\backslash\{0\}$ such that $(λ,v)$ is an eigenpair of the matrix $A(v)$ and we propose a generalization of inverse iteration for eigenvalue problems with this type of eigenvector nonlinearity. The convergence of the proposed method is studied and several convergence properties are shown to be analogous to inverse iteration for standard eigenvalue problems, including local convergence properties. The algorithm is also shown to be equivalent to a particular discretization of an associated ordinary differential equation, if the shift is chosen in a particular way. The algorithm is adapted to a variant of the Schrödinger equation known as the Gross-Pitaevskii equation. We use numerical simulations toillustrate the convergence properties, as well as the efficiency of the algorithm and the adaption. △ Less

Submitted 3 December, 2012; originally announced December 2012.

The choice of basic variables in current-density functional theory

Authors: Erik I. Tellgren, Simen Kvaal, Espen Sagvolden, Ulf Ekström, Andrew M. Teale, Trygve Helgaker

Abstract: The selection of basic variables in current-density functional theory and formal properties of the resulting formulations are critically examined. Focus is placed on the extent to which the Hohenberg--Kohn theorem, constrained-search approach and Lieb's formulation (in terms of convex and concave conjugation) of standard density-functional theory can be generalized to provide foundations for curre… ▽ More The selection of basic variables in current-density functional theory and formal properties of the resulting formulations are critically examined. Focus is placed on the extent to which the Hohenberg--Kohn theorem, constrained-search approach and Lieb's formulation (in terms of convex and concave conjugation) of standard density-functional theory can be generalized to provide foundations for current-density functional theory. For the well-known case with the gauge-dependent paramagnetic current density as a basic variable, we find that the resulting total energy functional is not concave. It is shown that a simple redefinition of the scalar potential restores concavity and enables the application of convex analysis and convex/concave conjugation. As a result, the solution sets arising in potential-optimization problems can be given a simple characterization. We also review attempts to establish theories with the physical current density as a basic variable. Despite the appealing physical motivation behind this choice of basic variables, we find that the mathematical foundations of the theories proposed to date are unsatisfactory. Moreover, the analogy to standard density-functional theory is substantially weaker as neither the constrained-search approach nor the convex analysis framework carry over to a theory making use of the physical current density. △ Less

Submitted 8 October, 2012; originally announced October 2012.

Ab initio quantum dynamics using coupled-cluster

Authors: Simen Kvaal

Abstract: The curse of dimensionality (COD) limits the current state-of-the-art {\it ab initio} propagation methods for non-relativistic quantum mechanics to relatively few particles. For stationary structure calculations, the coupled-cluster (CC) method overcomes the COD in the sense that the method scales polynomially with the number of particles while still being size-consistent and extensive. We general… ▽ More The curse of dimensionality (COD) limits the current state-of-the-art {\it ab initio} propagation methods for non-relativistic quantum mechanics to relatively few particles. For stationary structure calculations, the coupled-cluster (CC) method overcomes the COD in the sense that the method scales polynomially with the number of particles while still being size-consistent and extensive. We generalize the CC method to the time domain while allowing the single-particle functions to vary in an adaptive fashion as well, thereby creating a highly flexible, polynomially scaling approximation to the time-dependent Schrödinger equation. The method inherits size-consistency and extensivity from the CC method. The method is dubbed orbital-adaptive time-dependent coupled-cluster (OATDCC), and is a hierarchy of approximations to the now standard multi-configurational time-dependent Hartree method for fermions. A numerical experiment is also given. △ Less

Submitted 3 February, 2012; v1 submitted 26 January, 2012; originally announced January 2012.

Comments: 5 figures

Multiconfigurational time-dependent Hartree method for describing particle loss due to absorbing boundary conditions

Authors: Simen Kvaal

Abstract: Absorbing boundary conditions in the form of a complex absorbing potential are routinely introduced in the Schrödinger equation to limit the computational domain or to study reactive scattering events using the multi-configurational time-dependent Hartree method (MCTDH). However, it is known that a pure wave-function description does not allow the modeling and propagation of the remnants of a syst… ▽ More Absorbing boundary conditions in the form of a complex absorbing potential are routinely introduced in the Schrödinger equation to limit the computational domain or to study reactive scattering events using the multi-configurational time-dependent Hartree method (MCTDH). However, it is known that a pure wave-function description does not allow the modeling and propagation of the remnants of a system of which some parts are removed by the absorbing boundary. It was recently shown [S. Selstø and S. Kvaal, J. Phys. B: At. Mol. Opt. Phys. {\bfseries 43} (2010), 065004] that a master equation of Lindblad form was necessary for such a description. We formulate a multiconfigurational time-dependent Hartree method for this master equation, usable for any quantum system composed of any mixture of species. The formulation is a strict generalization of pure-state propagation using standard MCTDH. We demonstrate the formulation with a numerical experiment. △ Less

Submitted 7 March, 2011; v1 submitted 18 February, 2011; originally announced February 2011.

Comments: 5 figures, 13 pages

Ab initio computation of circular quantum dots

Authors: M. Pedersen Lohne, G. Hagen, M. Hjorth-Jensen, S. Kvaal, F. Pederiva

Abstract: We perform coupled-cluster and diffusion Monte Carlo calculations of the energies of circular quantum dots up to 20 electrons. The coupled-cluster calculations include triples corrections and a renormalized Coulomb interaction defined for a given number of low-lying oscillator shells. Using such a renormalized Coulomb interaction brings the coupled-cluster calculations with triples correlations in… ▽ More We perform coupled-cluster and diffusion Monte Carlo calculations of the energies of circular quantum dots up to 20 electrons. The coupled-cluster calculations include triples corrections and a renormalized Coulomb interaction defined for a given number of low-lying oscillator shells. Using such a renormalized Coulomb interaction brings the coupled-cluster calculations with triples correlations in excellent agreement with the diffusion Monte Carlo calculations. This opens up perspectives for doing ab initio calculations for much larger systems of electrons. △ Less

Submitted 24 September, 2010; originally announced September 2010.

Comments: 17 pages, 4 figures, submitted to Physical Review B

Journal ref: Phys.Rev.B84:115302,2011

Computing singularities of perturbation series

Authors: Simen Kvaal, Elias Jarlebring, Wim Michiels

Abstract: Many properties of current \emph{ab initio} approaches to the quantum many-body problem, both perturbational or otherwise, are related to the singularity structure of Rayleigh--Schrödinger perturbation theory. A numerical procedure is presented that in principle computes the complete set of singularities, including the dominant singularity which limits the radius of convergence. The method approxi… ▽ More Many properties of current \emph{ab initio} approaches to the quantum many-body problem, both perturbational or otherwise, are related to the singularity structure of Rayleigh--Schrödinger perturbation theory. A numerical procedure is presented that in principle computes the complete set of singularities, including the dominant singularity which limits the radius of convergence. The method approximates the singularities as eigenvalues of a certain generalized eigenvalue equation which is solved using iterative techniques. It relies on computation of the action of the perturbed Hamiltonian on a vector, and does not rely on the terms in the perturbation series. Some illustrative model problems are studied, including a Helium-like model with $δ$-function interactions for which Møller--Plesset perturbation theory is considered and the radius of convergence found. △ Less

Submitted 15 September, 2010; originally announced September 2010.

Comments: 11 figures, submitted

Adaptive time-dependent coupled cluster method for multiparticle wavepacket dynamics

Authors: Simen Kvaal

Abstract: A dynamical formulation of coupled cluster theory is derived using a variational principle. By allowing time-dependent single-particle functions, a high degree of adaptivity is introduced, allowing complex systems to be simulated with high accuracy. Equations of motion are derived which are shown to be suitable for computer implementation. The method, called adaptive time-dependent coupled cluster… ▽ More A dynamical formulation of coupled cluster theory is derived using a variational principle. By allowing time-dependent single-particle functions, a high degree of adaptivity is introduced, allowing complex systems to be simulated with high accuracy. Equations of motion are derived which are shown to be suitable for computer implementation. The method, called adaptive time-dependent coupled cluster, is a strict generalization of the formulation used in standard coupled cluster response theory, and also represents a systematic hierarchy of size-consistent approximations, including standard time-dependent Hartree--Fock as a trivial case. △ Less

Submitted 24 November, 2014; v1 submitted 23 August, 2010; originally announced August 2010.

Comments: Paper has been temporarily withdrawn due to developments in theory. The updated version is at arXiv:1201.5548

The role of final state correlation in double ionization of helium: a master equation approach

Authors: Sølve Selstø, Tore Birkeland, Simen Kvaal, Raymond Nepstad, Morten Førre

Abstract: The process of nonsequential two-photon double ionization of helium is studied by two complementary numerical approaches. First, the time-dependent Schr{ö}dinger equation is solved and the final wave function is analyzed in terms of projection onto eigenstates of the uncorrelated Hamiltonian, i.e., with no electron-electron interaction included in the final states. Then, the double ionization prob… ▽ More The process of nonsequential two-photon double ionization of helium is studied by two complementary numerical approaches. First, the time-dependent Schr{ö}dinger equation is solved and the final wave function is analyzed in terms of projection onto eigenstates of the uncorrelated Hamiltonian, i.e., with no electron-electron interaction included in the final states. Then, the double ionization probability is found by means of a recently developed approach in which the concept of absorbing boundaries has been generalized to apply to systems consisting of more than one particle. This generalization is achieved through the Lindblad equation. A model of reduced dimensionality, which describes the process at a qualitative level, has been used. The agreement between the methods provides a strong indication that procedures using projections onto uncorrelated continuum states are adequate when extracting total cross sections for the direct double ionization process. △ Less

Submitted 10 September, 2010; v1 submitted 10 June, 2010; originally announced June 2010.

Comments: 8 pages, 4 figures

Many-body interactions and nuclear structure

Authors: M. Hjorth-Jensen, D. J. Dean, G. Hagen, S. Kvaal

Abstract: This article presents several challenges to nuclear many-body theory and our understanding of the stability of nuclear matte r. In order to achieve this, we present five different cases, starting with an idealized toy model. These cases expose problems that need to be understood in order to match recent advances in nuclear theory with current experimental programs in low-energy nuclear physics. I… ▽ More This article presents several challenges to nuclear many-body theory and our understanding of the stability of nuclear matte r. In order to achieve this, we present five different cases, starting with an idealized toy model. These cases expose problems that need to be understood in order to match recent advances in nuclear theory with current experimental programs in low-energy nuclear physics. In particular, we focus on our current understanding, or lack thereof, of many-body forces, and how they evolve as functions of the number of particles . We provide examples of discrepancies between theory and experiment and outline some selected perspectives for future research directions. △ Less

Submitted 7 March, 2010; originally announced March 2010.

Comments: 26 pages and 6 figures, IOP latex style. Subimtted to Journal of Physics G, Focus issue on open problems in Nuclear structure theory

Journal ref: J.Phys.G37:064035,2010

Manifestly gauge invariant discretizations of the Schrödinger equation

Authors: Tore Gunnar Halvorsen, Simen Kvaal

Abstract: Grid-based discretizations of the time dependent Schrödinger equation coupled to an external magnetic field are converted to manifest gauge invariant discretizations. This is done using generalizations of ideas used in classical lattice gauge theory, and the process defined is applicable to a large class of discretized differential operators. In particular, popular discretizations such as pseudo… ▽ More Grid-based discretizations of the time dependent Schrödinger equation coupled to an external magnetic field are converted to manifest gauge invariant discretizations. This is done using generalizations of ideas used in classical lattice gauge theory, and the process defined is applicable to a large class of discretized differential operators. In particular, popular discretizations such as pseudospectral discretizations using the fast Fourier transform can be transformed to gauge invariant schemes. Also generic gauge invariant versions of generic time integration methods are considered, enabling completely gauge invariant calculations of the time dependent Schrödinger equation. Numerical examples illuminating the differences between a gauge invariant discretization and conventional discretization procedures are also presented. △ Less

Submitted 19 June, 2009; v1 submitted 26 May, 2009; originally announced May 2009.

Comments: 22 pages, 3 figures

Absorbing boundary conditions for dynamical many-body quantum systems

Authors: Sølve Selstø, Simen Kvaal

Abstract: In numerical studies of the dynamics of unbound quantum mechanical systems, absorbing boundary conditions are frequently applied. Although this certainly provides a useful tool in facilitating the description of the system, its applications to systems consisting of more than one particle is problematic. This is due to the fact that all information about the system is lost upon absorption of one… ▽ More In numerical studies of the dynamics of unbound quantum mechanical systems, absorbing boundary conditions are frequently applied. Although this certainly provides a useful tool in facilitating the description of the system, its applications to systems consisting of more than one particle is problematic. This is due to the fact that all information about the system is lost upon absorption of one particle; a formalism based solely on the Scrh{ö}dinger equation is not able to describe the remainder of the system as particles are lost. Here we demonstrate how the dynamics of a quantum system with a given number of identical fermions may be described in a manner which allows for particle loss. A consistent formalism which incorporates the evolution of sub-systems with a reduced number of particles is constructed through the Lindblad equation. Specifically, the transition from an $N$-particle system to an $(N-1)$-particle system due to a complex absorbing potential is achieved by relating the Lindblad operators to annihilation operators. The method allows for a straight forward interpretation of how many constituent particles have left the system after interaction. We illustrate the formalism using one-dimensional two-particle model problems. △ Less

Submitted 28 January, 2010; v1 submitted 14 April, 2009; originally announced April 2009.

Comments: 16 pages, 5 figures

Open source FCI code for quantum dots and effective interactions

Authors: Simen Kvaal

Abstract: We describe OpenFCI, an open source implementation of the full configuration-interaction method (FCI) for two-dimensional quantum dots with optional use of effective renormalized interactions. The code is written in C++ and is available under the Gnu General Public License. The code and core libraries are well documented and structured in a way such that customizations and generalizations to oth… ▽ More We describe OpenFCI, an open source implementation of the full configuration-interaction method (FCI) for two-dimensional quantum dots with optional use of effective renormalized interactions. The code is written in C++ and is available under the Gnu General Public License. The code and core libraries are well documented and structured in a way such that customizations and generalizations to other systems and numerical methods are easy tasks. As examples we provide a matrix element tabulation program and an implementation of a simple model from nuclear physics, in addition to the quantum dot application itself. △ Less

Submitted 15 October, 2008; originally announced October 2008.

Comments: 13 pages, 2 figures

Harmonic oscillator eigenfunction expansions, quantum dots, and effective interactions

Authors: Simen Kvaal

Abstract: We give a thorough analysis of the convergence properties of the configuration-interaction method as applied to parabolic quantum dots among other systems, including \emph{a priori} error estimates. The method converges slowly in general, and in order to overcome this, we propose to use an effective two-body interaction well-known from nuclear physics. Through numerical experiments we demonstrat… ▽ More We give a thorough analysis of the convergence properties of the configuration-interaction method as applied to parabolic quantum dots among other systems, including \emph{a priori} error estimates. The method converges slowly in general, and in order to overcome this, we propose to use an effective two-body interaction well-known from nuclear physics. Through numerical experiments we demonstrate a significant increase in accuracy of the configuration interaction method. △ Less

Submitted 15 May, 2009; v1 submitted 15 August, 2008; originally announced August 2008.

Comments: 3 figures

Geometry of effective Hamiltonians

Authors: Simen Kvaal

Abstract: We give a complete geometrical description of the effective Hamiltonians common in nuclear shell model calculations. By recasting the theory in a manifestly geometric form, we reinterpret and clarify several points. Some of these results are hitherto unknown or unpublished. In particular, commuting observables and symmetries are discussed in detail. Simple and explicit proofs are given, and nume… ▽ More We give a complete geometrical description of the effective Hamiltonians common in nuclear shell model calculations. By recasting the theory in a manifestly geometric form, we reinterpret and clarify several points. Some of these results are hitherto unknown or unpublished. In particular, commuting observables and symmetries are discussed in detail. Simple and explicit proofs are given, and numerical algorithms are proposed, that improve and stabilize common methods used today. △ Less

Submitted 13 August, 2008; originally announced August 2008.

Comments: 1 figure

Journal ref: Phys.Rev.C78:044330,2008

Effective interactions and large-scale diagonalization for quantum dots

Authors: Simen Kvaal, Morten Hjorth-Jensen, Halvor Moll Nilsen

Abstract: The widely used large-scale diagonalization method using harmonic oscillator basis functions (an instance of the Rayleigh-Ritz method, also called a spectral method, configuration-interaction method, or ``exact diagonalization'' method) is systematically analyzed using results for the convergence of Hermite function series. We apply this theory to a Hamiltonian for a one-dimensional model of a q… ▽ More The widely used large-scale diagonalization method using harmonic oscillator basis functions (an instance of the Rayleigh-Ritz method, also called a spectral method, configuration-interaction method, or ``exact diagonalization'' method) is systematically analyzed using results for the convergence of Hermite function series. We apply this theory to a Hamiltonian for a one-dimensional model of a quantum dot. The method is shown to converge slowly, and the non-smooth character of the interaction potential is identified as the main problem with the chosen basis, while on the other hand its important advantages are pointed out. An effective interaction obtained by a similarity transformation is proposed for improving the convergence of the diagonalization scheme, and numerical experiments are performed to demonstrate the improvement. Generalizations to more particles and dimensions are discussed. △ Less

Submitted 22 April, 2007; v1 submitted 19 April, 2007; originally announced April 2007.

Comments: 7 figures, submitted to Physical Review B Single reference error fixed

Journal ref: Phys. Rev. B 76 085421 (2007)