BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Simen Kvaal

Simen Kvaal appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
physics.chem-phquant-phphysics.comp-phmath-phmath.MPmath.NAnucl-thcond-mat.mtrl-sciNumerical Analysisphysics.atom-phcond-mat.mes-hallcond-mat.quant-gascond-mat.str-el

Catalog footprint

What is connected

19works
13topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

19 published item(s)

preprint2026arXiv

Reduced Density Matrix Functional Theory And A Reduced Formulation Of Density Functional Theory

A mathematical framework for reduced density matrix functional theory (RDMFT) is proposed. The work is inspired by and generalizes the work by E.H.~Lieb [E.H. Lieb, Int. J. Quant. Chem. 24(1983), pp.243--277] on density-functional theory (DFT). We introduce a Banach space for density matrices with finite kinetic energy. The dual space is a rich class of single-particle potentials, i.e., Hermitian forms. The ground state energy of an $N$-fermion system with external forces given by any such Hermitian form is expressed as the Legendre--Fenchel transform of a convex and lower semicontinuous ``universal'' reduced density matrix functional. The formalism is employed to provide a mathematical framework for density-functional theory (DFT). The main tool here is a rigorous definition of diagonals of reduced density matrices. The result is a refinement of Lieb's results on DFT applicable to a wide variety of models.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Linear and nonlinear optical properties from TDOMP2 theory

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Moreau--Yosida regularization in DFT

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Numerical stability of time-dependent coupled-cluster methods for many-electron dynamics in intense laser pulses

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Adaptive time-dependent coupled cluster method for multiparticle wavepacket dynamics

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Differentiable but exact formulation of density-functional theory

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Ab initio quantum dynamics using coupled-cluster

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

An inverse iteration method for eigenvalue problems with eigenvector nonlinearities

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The choice of basic variables in current-density functional theory

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Absorbing boundary conditions for dynamical many-body quantum systems

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Computing singularities of perturbation series

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Harmonic oscillator eigenfunction expansions, quantum dots, and effective interactions

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Manifestly gauge invariant discretizations of the Schrödinger equation

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.

0 citations0 reviewsNo code linkedOpen access