Showing 1–31 of 31 results for author: Levitt, A

Linear response and resonances in adiabatic time-dependent density functional theory

Authors: Mi-Song Dupuy, Eloïse Letournel, Antoine Levitt

Abstract: We consider the electrons of a molecule in the adiabatic time-dependent density functional theory approximation. We establish the well-posedness of the time evolution and its linear response close to a non-degenerate ground state, and prove the appearance of resonances at relevant frequencies. The main mathematical difficulty is due to the structure of the linearized equations, which are not compl… ▽ More We consider the electrons of a molecule in the adiabatic time-dependent density functional theory approximation. We establish the well-posedness of the time evolution and its linear response close to a non-degenerate ground state, and prove the appearance of resonances at relevant frequencies. The main mathematical difficulty is due to the structure of the linearized equations, which are not complex-linear. We bypass this difficulty by reformulating the linearized problem as a real Hamiltonian system, whose stability is ensured by the second-order optimality conditions on the energy. △ Less

Submitted 24 July, 2024; originally announced August 2024.

Some mathematical insights on Density Matrix Embedding Theory

Authors: Eric Cancès, Fabian M. Faulstich, Alfred Kirsch, Eloïse Letournel, Antoine Levitt

Abstract: This article provides the first mathematical analysis of the Density Matrix Embedding Theory (DMET) method. We prove that, under certain assumptions, (i) the exact ground-state density matrix is a fixed-point of the DMET map for non-interacting systems, (ii) there exists a unique physical solution in the weakly-interacting regime, and (iii) DMET is exact at first order in the coupling parameter. W… ▽ More This article provides the first mathematical analysis of the Density Matrix Embedding Theory (DMET) method. We prove that, under certain assumptions, (i) the exact ground-state density matrix is a fixed-point of the DMET map for non-interacting systems, (ii) there exists a unique physical solution in the weakly-interacting regime, and (iii) DMET is exact at first order in the coupling parameter. We provide numerical simulations to support our results and comment on the physical meaning of the assumptions under which they hold true. We show that the violation of these assumptions may yield multiple solutions of the DMET equations. We moreover introduce and discuss a specific N-representability problem inherent to DMET. △ Less

Submitted 2 October, 2023; v1 submitted 25 May, 2023; originally announced May 2023.

Comments: 40 pages, 11 figures

MSC Class: 35Q40; 39B42; 49M27

Computing photoionization spectra in Gaussian basis sets

Authors: Ivan Duchemin, Antoine Levitt

Abstract: We present a method to compute the photoionization spectra of atoms and molecules in linear response time-dependent density functional theory. The electronic orbital variations corresponding to ionized electrons are expanded on a basis set of delocalized functions obtained as the solution of the inhomogeneous Helmholtz equation with gaussian basis set functions as right-hand side. The resulting sc… ▽ More We present a method to compute the photoionization spectra of atoms and molecules in linear response time-dependent density functional theory. The electronic orbital variations corresponding to ionized electrons are expanded on a basis set of delocalized functions obtained as the solution of the inhomogeneous Helmholtz equation with gaussian basis set functions as right-hand side. The resulting scheme is able to reproduce photoionization spectra without any need for artificial regularization or localization. We demonstrate that it is able to produce accurate spectra for semilocal exchange-correlation functionals even using relatively small standard gaussian basis sets. △ Less

Submitted 23 May, 2023; originally announced May 2023.

A robust, open-source implementation of the locally optimal block preconditioned conjugate gradient for large eigenvalue problems in quantum chemistry

Authors: Tommaso Nottoli, Ivan Giannì, Antoine Levitt, Filippo Lipparini

Abstract: We present two open-source implementations of the Locally Optimal Block Preconditioned Conjugate Gradient (LOBPCG) algorithm to find a few eigenvalues and eigenvectors of large, possibly sparse matrices. We then test LOBPCG for various quantum chemistry problems, encompassing medium to large, dense to sparse, wellbehaved to ill-conditioned ones, where the standard method typically used is Davidson… ▽ More We present two open-source implementations of the Locally Optimal Block Preconditioned Conjugate Gradient (LOBPCG) algorithm to find a few eigenvalues and eigenvectors of large, possibly sparse matrices. We then test LOBPCG for various quantum chemistry problems, encompassing medium to large, dense to sparse, wellbehaved to ill-conditioned ones, where the standard method typically used is Davidson's diagonalization. Numerical tests show that, while Davidson's method remains the best choice for most applications in quantum chemistry, LOBPCG represents a competitive alternative, especially when memory is an issue, and can even outperform Davidson for ill-conditioned, non diagonally dominant problems. △ Less

Submitted 11 May, 2023; originally announced May 2023.

Comments: Theoretical Chemistry Accounts: Theory, Computation, and Modeling, In press

Photoionization and core resonances from range-separated time-dependent density-functional theory for open-shell states: Example of the lithium atom

Authors: Julien Toulouse, Karno Schwinn, Felipe Zapata, Antoine Levitt, Eric Cancès, Eleonora Luppi

Abstract: We consider the calculations of photoionization spectra and core resonances of open-shell systems using range-separated time-dependent density-functional theory. Specifically, we use the time-dependent range-separated hybrid (TDRSH) scheme, combining a long-range Hartree-Fock (HF) exchange potential and kernel with a short-range potential and kernel from a local density-functional approximation, a… ▽ More We consider the calculations of photoionization spectra and core resonances of open-shell systems using range-separated time-dependent density-functional theory. Specifically, we use the time-dependent range-separated hybrid (TDRSH) scheme, combining a long-range Hartree-Fock (HF) exchange potential and kernel with a short-range potential and kernel from a local density-functional approximation, and the time-dependent locally range-separated hybrid (TDLRSH) scheme, which uses a local range-separation parameter. To efficiently perform the calculations, we formulate a spin-unrestricted linear-response Sternheimer approach in a non-orthogonal B-spline basis set and using appropriate frequency-dependent boundary conditions. We illustrate this approach on the Li atom, which suggests that TDRSH and TDLRSH are adequate simple methods for estimating single-electron photoionization spectra of open-shell systems. △ Less

Submitted 14 November, 2022; originally announced November 2022.

Comments: arXiv admin note: substantial text overlap with arXiv:2203.05951

Ab initio Canonical Sampling based on Variational Inference

Authors: Alois Castellano, Francois Bottin, Johann Bouchet, Antoine Levitt, Gabriel Stoltz

Abstract: Finite temperature calculations, based on ab initio molecular dynamics (AIMD) simulations, are a powerful tool able to predict material properties that cannot be deduced from ground state calculations. However, the high computational cost of AIMD limits its applicability for large or complex systems. To circumvent this limitation we introduce a new method named Machine Learning Assisted Canonical… ▽ More Finite temperature calculations, based on ab initio molecular dynamics (AIMD) simulations, are a powerful tool able to predict material properties that cannot be deduced from ground state calculations. However, the high computational cost of AIMD limits its applicability for large or complex systems. To circumvent this limitation we introduce a new method named Machine Learning Assisted Canonical Sampling (MLACS), which accelerates the sampling of the Born--Oppenheimer potential surface in the canonical ensemble. Based on a self-consistent variational procedure, the method iteratively trains a Machine Learning Interatomic Potential to generate configurations that approximate the canonical distribution of positions associated with the ab initio potential energy. By proving the reliability of the method on anharmonic systems, we show that the method is able to reproduce the results of AIMD with an ab initio accuracy at a fraction of its computational cost. △ Less

Submitted 20 October, 2022; originally announced October 2022.

Journal ref: Phys. Rev. B 106, L161110 (2022)

Numerical stability and efficiency of response property calculations in density functional theory

Authors: Eric Cancès, Michael F. Herbst, Gaspard Kemlin, Antoine Levitt, Benjamin Stamm

Abstract: Response calculations in density functional theory aim at computing the change in ground-state density induced by an external perturbation. At finite temperature these are usually performed by computing variations of orbitals, which involve the iterative solution of potentially badly-conditioned linear systems, the Sternheimer equations. Since many sets of variations of orbitals yield the same var… ▽ More Response calculations in density functional theory aim at computing the change in ground-state density induced by an external perturbation. At finite temperature these are usually performed by computing variations of orbitals, which involve the iterative solution of potentially badly-conditioned linear systems, the Sternheimer equations. Since many sets of variations of orbitals yield the same variation of density matrix this involves a choice of gauge. Taking a numerical analysis point of view we present the various gauge choices proposed in the literature in a common framework and study their stability. Beyond existing methods we propose a new approach, based on a Schur complement using extra orbitals from the self-consistent-field calculations, to improve the stability and efficiency of the iterative solution of Sternheimer equations. We show the success of this strategy on nontrivial examples of practical interest, such as Heusler transition metal alloy compounds, where savings of around 40% in the number of required cost-determining Hamiltonian applications have been achieved. △ Less

Submitted 16 February, 2023; v1 submitted 10 October, 2022; originally announced October 2022.

Journal ref: Letters in Mathematical Physics, 2023, 113 (1)

A priori error analysis of linear and nonlinear periodic Schr{ö}dinger equations with analytic potentials

Authors: Eric Cancès, Gaspard Kemlin, Antoine Levitt

Abstract: This paper is concerned with the numerical analysis of linear and nonlinear Schr{ö}dinger equations with analytic potentials. While the regularity of the potential (and the source term when there is one) automatically conveys to the solution in the linear cases, this is no longer true in general in the nonlinear case. We also study the rate of convergence of the planewave (Fourier) discretization… ▽ More This paper is concerned with the numerical analysis of linear and nonlinear Schr{ö}dinger equations with analytic potentials. While the regularity of the potential (and the source term when there is one) automatically conveys to the solution in the linear cases, this is no longer true in general in the nonlinear case. We also study the rate of convergence of the planewave (Fourier) discretization method for computing numerical approximations of the solution. △ Less

Submitted 21 November, 2023; v1 submitted 10 June, 2022; originally announced June 2022.

Efficient extraction of resonant states in systems with defects

Authors: Ivan Duchemin, Luigi Genovese, Eloïse Letournel, Antoine Levitt, Simon Ruget

Abstract: We introduce a new numerical method to compute resonances induced by localized defects in crystals. This method solves an integral equation in the defect region to compute analytic continuations of resolvents. Such an approach enables one to express the resonance in terms of a "resonance source", a function that is strictly localized within the defect region. The kernel of the integral equation, t… ▽ More We introduce a new numerical method to compute resonances induced by localized defects in crystals. This method solves an integral equation in the defect region to compute analytic continuations of resolvents. Such an approach enables one to express the resonance in terms of a "resonance source", a function that is strictly localized within the defect region. The kernel of the integral equation, to be applied on such a source term, is the Green function of the perfect crystal, which we show can be computed efficiently by a complex deformation of the Brillouin zone, named Brillouin Complex Deformation (BCD), thereby extending to reciprocal space the concept of complex coordinate transformations. △ Less

Submitted 16 September, 2022; v1 submitted 23 March, 2022; originally announced April 2022.

Photoionization and core resonances from range-separated density-functional theory: General formalism and example of the beryllium atom

Authors: Karno Schwinn, Felipe Zapata, Antoine Levitt, Eric Cancès, Eleonora Luppi, Julien Toulouse

Abstract: We explore the merits of linear-response range-separated time-dependent density-functional theory (TDDFT) for the calculation of photoionization spectra. We consider two variants of range-separated TDDFT, namely the timedependent range-separated hybrid (TDRSH) scheme which uses a global range-separation parameter and the timedependent locally range-separated hybrid (TDLRSH) which uses a local rang… ▽ More We explore the merits of linear-response range-separated time-dependent density-functional theory (TDDFT) for the calculation of photoionization spectra. We consider two variants of range-separated TDDFT, namely the timedependent range-separated hybrid (TDRSH) scheme which uses a global range-separation parameter and the timedependent locally range-separated hybrid (TDLRSH) which uses a local range-separation parameter, and compare with standard time-dependent local-density approximation (TDLDA) and time-dependent Hartree-Fock (TDHF). We show how to calculate photoionization spectra with these methods using the Sternheimer approach formulated in a nonorthogonal B-spline basis set with appropriate frequency-dependent boundary conditions. We illustrate these methods on the photoionization spectrum of the Be atom, focusing in particular on the core resonances. Both the TDRSH and TDLRSH photoionization spectra are found to constitute a large improvement over the TDLDA photoionization spectrum and a more modest improvement over the TDHF photoionization spectrum. △ Less

Submitted 11 March, 2022; originally announced March 2022.

Practical error bounds for properties in plane-wave electronic structure calculations

Authors: Eric Cancès, Geneviève Dusson, Gaspard Kemlin, Antoine Levitt

Abstract: We propose accurate computable error bounds for quantities of interest in plane-wave electronic structure calculations, in particular ground-state density matrices and energies, and interatomic forces. These bounds are based on an estimation of the error in terms of the residual of the solved equations, which is then efficiently approximated with computable terms. After providing coarse bounds bas… ▽ More We propose accurate computable error bounds for quantities of interest in plane-wave electronic structure calculations, in particular ground-state density matrices and energies, and interatomic forces. These bounds are based on an estimation of the error in terms of the residual of the solved equations, which is then efficiently approximated with computable terms. After providing coarse bounds based on an analysis of the inverse Jacobian, we improve on these bounds by solving a linear problem in a small dimension that involves a Schur complement. We numerically show how accurate these bounds are on a few representative materials, namely silicon, gallium arsenide and titanium dioxide. △ Less

Submitted 21 June, 2022; v1 submitted 2 November, 2021; originally announced November 2021.

A robust and efficient line search for self-consistent field iterations

Authors: Michael F. Herbst, Antoine Levitt

Abstract: We propose a novel adaptive damping algorithm for the self-consistent field (SCF) iterations of Kohn-Sham density-functional theory, using a backtracking line search to automatically adjust the damping in each SCF step. This line search is based on a theoretically sound, accurate and inexpensive model for the energy as a function of the damping parameter. In contrast to usual damped SCF schemes, t… ▽ More We propose a novel adaptive damping algorithm for the self-consistent field (SCF) iterations of Kohn-Sham density-functional theory, using a backtracking line search to automatically adjust the damping in each SCF step. This line search is based on a theoretically sound, accurate and inexpensive model for the energy as a function of the damping parameter. In contrast to usual damped SCF schemes, the resulting algorithm is fully automatic and does not require the user to select a damping. We successfully apply it to a wide range of challenging systems, including elongated supercells, surfaces and transition-metal alloys. △ Less

Submitted 1 March, 2022; v1 submitted 28 September, 2021; originally announced September 2021.

Finite-size effects in response functions of molecular systems

Authors: Mi-Song Dupuy, Antoine Levitt

Abstract: We consider an electron in a localized potential submitted to a weak external, timedependent field. In the linear response regime, the response function can be computed using Kubo's formula. In this paper, we consider the numerical approximation of the response function by means of a truncation to a finite region of space. This is necessarily a singular approximation because of the discreteness of… ▽ More We consider an electron in a localized potential submitted to a weak external, timedependent field. In the linear response regime, the response function can be computed using Kubo's formula. In this paper, we consider the numerical approximation of the response function by means of a truncation to a finite region of space. This is necessarily a singular approximation because of the discreteness of the spectrum of the truncated Hamiltonian, and in practice a regularization (smoothing) has to be used. Our results provide error estimates for the response function past the ionization threshold with respect to both the smoothing parameter and the size of the computational domain. △ Less

Submitted 19 February, 2021; originally announced February 2021.

Black-box inhomogeneous preconditioning for self-consistent field iterations in density functional theory

Authors: Michael F. Herbst, Antoine Levitt

Abstract: We propose a new preconditioner based on the local density of states for computing the self-consistent problem in Kohn-Sham density functional theory. This preconditioner is inexpensive and able to cure the long-range charge sloshing known to hamper convergence in large, inhomogeneous systems such as clusters and surfaces. It is based on a parameter-free and physically motivated approximation to t… ▽ More We propose a new preconditioner based on the local density of states for computing the self-consistent problem in Kohn-Sham density functional theory. This preconditioner is inexpensive and able to cure the long-range charge sloshing known to hamper convergence in large, inhomogeneous systems such as clusters and surfaces. It is based on a parameter-free and physically motivated approximation to the independent-particle susceptibility operator, appropriate for both metals and insulators. It can be extended to semiconductors by using the macroscopic electronic dielectric constant as a parameter in the model. We test our preconditioner successfully on inhomogeneous systems containing metals, insulators, semiconductors and vacuum. △ Less

Submitted 17 November, 2020; v1 submitted 3 September, 2020; originally announced September 2020.

A posteriori error estimation for the non-self-consistent Kohn-Sham equations

Authors: Michael F. Herbst, Antoine Levitt, Eric Cancès

Abstract: We address the problem of bounding rigorously the errors in the numerical solution of the Kohn-Sham equations due to (i) the finiteness of the basis set, (ii) the convergence thresholds in iterative procedures, (iii) the propagation of rounding errors in floating-point arithmetic. In this contribution, we compute fully-guaranteed bounds on the solution of the non-self-consistent equations in the p… ▽ More We address the problem of bounding rigorously the errors in the numerical solution of the Kohn-Sham equations due to (i) the finiteness of the basis set, (ii) the convergence thresholds in iterative procedures, (iii) the propagation of rounding errors in floating-point arithmetic. In this contribution, we compute fully-guaranteed bounds on the solution of the non-self-consistent equations in the pseudopotential approximation in a plane-wave basis set. We demonstrate our methodology by providing band structure diagrams of silicon annotated with error bars indicating the combined error. △ Less

Submitted 15 June, 2020; v1 submitted 28 April, 2020; originally announced April 2020.

Convergence analysis of direct minimization and self-consistent iterations

Authors: Eric Cancès, Gaspard Kemlin, Antoine Levitt

Abstract: This article is concerned with the numerical solution of subspace optimization problems, consisting of minimizing a smooth functional over the set of orthogonal projectors of fixed rank. Such problems are encountered in particular in electronic structure calculation (Hartree-Fock and Kohn-Sham Density Functional Theory -DFT- models). We compare from a numerical analysis perspective two simple repr… ▽ More This article is concerned with the numerical solution of subspace optimization problems, consisting of minimizing a smooth functional over the set of orthogonal projectors of fixed rank. Such problems are encountered in particular in electronic structure calculation (Hartree-Fock and Kohn-Sham Density Functional Theory -DFT- models). We compare from a numerical analysis perspective two simple representatives, the damped self-consistent field (SCF) iterations and the gradient descent algorithm, of the two classes of methods competing in the field: SCF and direct minimization methods. We derive asymptotic rates of convergence for these algorithms and analyze their dependence on the spectral gap and other properties of the problem. Our theoretical results are complemented by numerical simulations on a variety of examples, from toy models with tunable parameters to realistic Kohn-Sham computations. We also provide an example of chaotic behavior of the simple SCF iterations for a nonquadratic functional. △ Less

Submitted 27 October, 2020; v1 submitted 20 April, 2020; originally announced April 2020.

Hypocoercivity with Schur complements

Authors: E. Bernard, M. Fathi, A. Levitt, G. Stoltz

Abstract: We propose an approach to obtaining explicit estimates on the resolvent of hypocoercive operators by using Schur complements, rather than from an exponential decay of the evolution semigroup combined with a time integral. We present applications to Langevin-like dynamics and Fokker--Planck equations, as well as the linear Boltzmann equation (which is also the generator of randomized Hybrid Monte C… ▽ More We propose an approach to obtaining explicit estimates on the resolvent of hypocoercive operators by using Schur complements, rather than from an exponential decay of the evolution semigroup combined with a time integral. We present applications to Langevin-like dynamics and Fokker--Planck equations, as well as the linear Boltzmann equation (which is also the generator of randomized Hybrid Monte Carlo in molecular dynamics). In particular, we make precise the dependence of the resolvent bounds on the parameters of the dynamics and on the dimension. We also highlight the relationship of our method with other hypocoercive approaches. △ Less

Submitted 10 September, 2021; v1 submitted 2 March, 2020; originally announced March 2020.

Coherent electronic transport in periodic crystals

Authors: Eric Cancès, Clotilde Fermanian Kammerer, Antoine Levitt, Sami Siraj-Dine

Abstract: We consider independent electrons in a periodic crystal in their ground state, and turn on a uniform electric field at some prescribed time. We rigorously define the current per unit volume and study its properties using both linear response and adiabatic theory. Our results provide a unified framework for various phenomena such as the quantification of Hall conductivity of insulators with broken… ▽ More We consider independent electrons in a periodic crystal in their ground state, and turn on a uniform electric field at some prescribed time. We rigorously define the current per unit volume and study its properties using both linear response and adiabatic theory. Our results provide a unified framework for various phenomena such as the quantification of Hall conductivity of insulators with broken time-reversibility, the ballistic regime of electrons in metals, Bloch oscillations in the long-time response of metals, and the static conductivity of graphene. We identify explicitly the regime in which each holds. △ Less

Submitted 9 November, 2020; v1 submitted 5 February, 2020; originally announced February 2020.

Numerical construction of Wannier functions through homotopy

Authors: David Gontier, Antoine Levitt, Sami Siraj-Dine

Abstract: We provide a mathematically proven, simple and efficient algorithm to build localised Wannier functions, with the only requirement that the system has vanishing Chern numbers. Our algorithm is able to build localised Wannier for topological insulators such as the Kane-Mele model. It is based on an explicit and constructive proof of homotopies for the unitary group. We provide numerical tests valid… ▽ More We provide a mathematically proven, simple and efficient algorithm to build localised Wannier functions, with the only requirement that the system has vanishing Chern numbers. Our algorithm is able to build localised Wannier for topological insulators such as the Kane-Mele model. It is based on an explicit and constructive proof of homotopies for the unitary group. We provide numerical tests validating the methods for several systems, including the Kane-Mele model. △ Less

Submitted 17 December, 2018; originally announced December 2018.

Screening in the finite-temperature reduced Hartree-Fock model

Authors: Antoine Levitt

Abstract: We prove the existence of solutions of the reduced Hartree-Fock equations at finite temperature for a periodic crystal with a small defect, and show total screening of the defect charge by the electrons. We also show the convergence of the damped self-consistent field iteration using Kerker preconditioning to remove charge sloshing. As a crucial step of the proof, we define and study the propertie… ▽ More We prove the existence of solutions of the reduced Hartree-Fock equations at finite temperature for a periodic crystal with a small defect, and show total screening of the defect charge by the electrons. We also show the convergence of the damped self-consistent field iteration using Kerker preconditioning to remove charge sloshing. As a crucial step of the proof, we define and study the properties of the dielectric operator. △ Less

Submitted 27 July, 2020; v1 submitted 8 October, 2018; originally announced October 2018.

Journal ref: Archive for Rational Mechanics and Analysis, Springer Verlag, In press

Numerical quadrature in the Brillouin zone for periodic Schrodinger operators

Authors: Eric Cances, Virginie Ehrlacher, David Gontier, Antoine Levitt, Damiano Lombardi

Abstract: As a consequence of Bloch's theorem, the numerical computation of the fermionic ground state density matrices and energies of periodic Schrodinger operators involves integrals over the Brillouin zone. These integrals are difficult to compute numerically in metals due to discontinuities in the integrand. We perform an error analysis of several widely-used quadrature rules and smearing methods for B… ▽ More As a consequence of Bloch's theorem, the numerical computation of the fermionic ground state density matrices and energies of periodic Schrodinger operators involves integrals over the Brillouin zone. These integrals are difficult to compute numerically in metals due to discontinuities in the integrand. We perform an error analysis of several widely-used quadrature rules and smearing methods for Brillouin zone integration. We precisely identify the assumptions implicit in these methods and rigorously prove error bounds. Numerical results for two-dimensional periodic systems are also provided. Our results shed light on the properties of these numerical schemes, and provide guidance as to the appropriate choice of numerical parameters. △ Less

Submitted 23 September, 2021; v1 submitted 18 May, 2018; originally announced May 2018.

Variational formulation for Wannier functions with entangled band structure

Authors: Anil Damle, Antoine Levitt, Lin Lin

Abstract: Wannier functions provide a localized representation of spectral subspaces of periodic Hamiltonians, and play an important role for interpreting and accelerating Hartree-Fock and Kohn-Sham density functional theory calculations in quantum physics and chemistry. For systems with isolated band structure, the existence of exponentially localized Wannier functions and numerical algorithms for finding… ▽ More Wannier functions provide a localized representation of spectral subspaces of periodic Hamiltonians, and play an important role for interpreting and accelerating Hartree-Fock and Kohn-Sham density functional theory calculations in quantum physics and chemistry. For systems with isolated band structure, the existence of exponentially localized Wannier functions and numerical algorithms for finding them are well studied. In contrast, for systems with entangled band structure, Wannier functions must be generalized to span a subspace larger than the spectral subspace of interest to achieve favorable spatial locality. In this setting, little is known about the theoretical properties of these Wannier functions, and few algorithms can find them robustly. We develop a variational formulation to compute these generalized maximally localized Wannier functions. When paired with an initial guess based on the selected columns of the density matrix (SCDM) method, our method can robustly find Wannier functions for systems with entangled band structure. We formulate the problem as a constrained nonlinear optimization problem, and show how the widely used disentanglement procedure can be interpreted as a splitting method to approximately solve this problem. We demonstrate the performance of our method using real materials including silicon, copper, and aluminum. To examine more precisely the localization properties of Wannier functions, we study the free electron gas in one and two dimensions, where we show that the maximally-localized Wannier functions only decay algebraically. We also explain using a one dimensional example how to modify them to obtain super-algebraic decay. △ Less

Submitted 25 January, 2018; originally announced January 2018.

Comments: 26 pages, 10 figures

MSC Class: 65Z05; 82D25; 65F30; 65K10

Localised Wannier functions in metallic systems

Authors: Horia Cornean, David Gontier, Antoine Levitt, Domenico Monaco

Abstract: The existence and construction of exponentially localised Wannier functions for insulators is a well-studied problem. In comparison, the case of metallic systems has been much less explored, even though localised Wannier functions constitute an important and widely used tool for the numerical band interpolation of metallic condensed matter systems. In this paper we prove that, under generic condit… ▽ More The existence and construction of exponentially localised Wannier functions for insulators is a well-studied problem. In comparison, the case of metallic systems has been much less explored, even though localised Wannier functions constitute an important and widely used tool for the numerical band interpolation of metallic condensed matter systems. In this paper we prove that, under generic conditions, $N$ energy bands of a metal can be exactly represented by $N+1$ Wannier functions decaying faster than any polynomial. We also show that, in general, the lack of a spectral gap does not allow for exponential decay. △ Less

Submitted 21 December, 2017; originally announced December 2017.

Comments: 18 pages, 4 figures

MSC Class: 35Q40; 81Q30; 81Q70

Efficient spectral computation of the stationary states of rotating Bose-Einstein condensates by the preconditioned nonlinear conjugate gradient method

Authors: Xavier Antoine, Antoine Levitt, Qinglin Tang

Abstract: We propose a preconditioned nonlinear conjugate gradient method coupled with a spectral spatial dis-cretization scheme for computing the ground states (GS) of rotating Bose-Einstein condensates (BEC), modeled by the Gross-Pitaevskii Equation (GPE). We first start by reviewing the classical gradient flow (also known as imaginary time (IMT)) method which considers the problem from the PDE standpoint… ▽ More We propose a preconditioned nonlinear conjugate gradient method coupled with a spectral spatial dis-cretization scheme for computing the ground states (GS) of rotating Bose-Einstein condensates (BEC), modeled by the Gross-Pitaevskii Equation (GPE). We first start by reviewing the classical gradient flow (also known as imaginary time (IMT)) method which considers the problem from the PDE standpoint, leading to numerically solve a dissipative equation. Based on this IMT equation, we analyze the forward Euler (FE), Crank-Nicolson (CN) and the classical backward Euler (BE) schemes for linear problems and recognize classical power iterations, allowing us to derive convergence rates. By considering the alternative point of view of minimization problems, we propose the preconditioned gradient (PG) and conjugate gradient (PCG) methods for the GS computation of the GPE. We investigate the choice of the preconditioner, which plays a key role in the acceleration of the convergence process. The performance of the new algorithms is tested in 1D, 2D and 3D. We conclude that the PCG method outperforms all the previous methods, most particularly for 2D and 3D fast rotating BECs, while being simple to implement. △ Less

Submitted 7 November, 2016; originally announced November 2016.

Convergence and Cycling in Walker-type Saddle Search Algorithms

Authors: Antoine Levitt, Christoph Ortner

Abstract: Algorithms for computing local minima of smooth objective functions enjoy a mature theory as well as robust and efficient implementations. By comparison, the theory and practice of saddle search is destitute. In this paper we present results for idealized versions of the dimer and gentlest ascent (GAD) saddle search algorithms that show-case the limitations of what is theoretically achievable with… ▽ More Algorithms for computing local minima of smooth objective functions enjoy a mature theory as well as robust and efficient implementations. By comparison, the theory and practice of saddle search is destitute. In this paper we present results for idealized versions of the dimer and gentlest ascent (GAD) saddle search algorithms that show-case the limitations of what is theoretically achievable within the current class of saddle search algorithms: (1) we present an improved estimate on the region of attraction of saddles; and (2) we construct quasi-periodic solutions which indicate that it is impossible to obtain globally convergent variants of dimer and GAD type algorithms. △ Less

Submitted 29 July, 2016; originally announced July 2016.

Robust determination of maximally-localized Wannier functions

Authors: Éric Cancès, Antoine Levitt, Gianluca Panati, Gabriel Stoltz

Abstract: We propose an algorithm to determine Maximally Localized Wannier Functions (MLWFs). This algorithm, based on recent theoretical developments, does not require any physical input such as initial guesses for the Wannier functions, unlike popular schemes based on the projection method. We discuss how the projection method can fail on fine grids when the initial guesses are too far from MLWFs. We demo… ▽ More We propose an algorithm to determine Maximally Localized Wannier Functions (MLWFs). This algorithm, based on recent theoretical developments, does not require any physical input such as initial guesses for the Wannier functions, unlike popular schemes based on the projection method. We discuss how the projection method can fail on fine grids when the initial guesses are too far from MLWFs. We demonstrate that our algorithm is able to find localized Wannier functions through tests on two-dimensional systems, simplified models of semiconductors, and realistic DFT systems by interfacing with the Wannier90 code. We also test our algorithm on the Haldane and Kane-Mele models to examine how it fails in the presence of topological obstructions. △ Less

Submitted 24 January, 2017; v1 submitted 23 May, 2016; originally announced May 2016.

Report number: Roma01.Math.MP

Journal ref: Phys. Rev. B 95, 075114 (2017)

Parallel eigensolvers in plane-wave Density Functional Theory

Authors: Antoine Levitt, Marc Torrent

Abstract: We consider the problem of parallelizing electronic structure computations in plane-wave Density Functional Theory. Because of the limited scalability of Fourier transforms, parallelism has to be found at the eigensolver level. We show how a recently proposed algorithm based on Chebyshev polynomials can scale into the tens of thousands of processors, outperforming block conjugate gradient algorith… ▽ More We consider the problem of parallelizing electronic structure computations in plane-wave Density Functional Theory. Because of the limited scalability of Fourier transforms, parallelism has to be found at the eigensolver level. We show how a recently proposed algorithm based on Chebyshev polynomials can scale into the tens of thousands of processors, outperforming block conjugate gradient algorithms for large computations. △ Less

Submitted 7 October, 2014; v1 submitted 17 June, 2014; originally announced June 2014.

Solutions of the multiconfiguration Dirac-Fock equations

Authors: Antoine Levitt

Abstract: The multiconfiguration Dirac-Fock (MCDF) model uses a linear combination of Slater determinants to approximate the electronic $N$-body wave function of a relativistic molecular system, resulting in a coupled system of nonlinear eigenvalue equations, the MCDF equations. In this paper, we prove the existence of solutions of these equations in the weakly relativistic regime. First, using a new variat… ▽ More The multiconfiguration Dirac-Fock (MCDF) model uses a linear combination of Slater determinants to approximate the electronic $N$-body wave function of a relativistic molecular system, resulting in a coupled system of nonlinear eigenvalue equations, the MCDF equations. In this paper, we prove the existence of solutions of these equations in the weakly relativistic regime. First, using a new variational principle as well as results of Lewin on the multiconfiguration nonrelativistic model, and Esteban and Séré on the single-configuration relativistic model, we prove the existence of critical points for the associated energy functional, under the constraint that the occupation numbers are not too small. Then, this constraint can be removed in the weakly relativistic regime, and we obtain non-constrained critical points, i.e. solutions of the multiconfiguration Dirac-Fock equations. △ Less

Submitted 22 July, 2014; v1 submitted 16 May, 2013; originally announced May 2013.

Comments: Reviews in Mathematical Physics (2014) Accepté pour publication

Best constants in Lieb-Thirring inequalities: a numerical investigation

Authors: Antoine Levitt

Abstract: We investigate numerically the optimal constants in Lieb-Thirring inequalities by studying the associated maximization problem. We use a monotonic fixed-point algorithm and a finite element discretization to obtain trial potentials which provide lower bounds on the optimal constants. We examine the one-dimensional and radial cases in detail. Our numerical results provide new lower bounds, insight… ▽ More We investigate numerically the optimal constants in Lieb-Thirring inequalities by studying the associated maximization problem. We use a monotonic fixed-point algorithm and a finite element discretization to obtain trial potentials which provide lower bounds on the optimal constants. We examine the one-dimensional and radial cases in detail. Our numerical results provide new lower bounds, insight into the behavior of the maximizers and confirm some existing conjectures. Based on our numerical results, we formulate a complete conjecture about the best constants for all possible values of the parameters. △ Less

Submitted 7 June, 2012; originally announced June 2012.

Convergence of gradient-based algorithms for the Hartree-Fock equations

Authors: Antoine Levitt

Abstract: The numerical solution of the Hartree-Fock equations is a central problem in quantum chemistry for which numerous algorithms exist. Attempts to justify these algorithms mathematically have been made, notably in by Cances and Le Bris in 2000, but, to our knowledge, no complete convergence proof has been published. In this paper, we prove the convergence of a natural gradient algorithm, using a grad… ▽ More The numerical solution of the Hartree-Fock equations is a central problem in quantum chemistry for which numerous algorithms exist. Attempts to justify these algorithms mathematically have been made, notably in by Cances and Le Bris in 2000, but, to our knowledge, no complete convergence proof has been published. In this paper, we prove the convergence of a natural gradient algorithm, using a gradient inequality for analytic functionals due to Lojasiewicz. Then, expanding upon the analysis of Cances and Le Bris, we prove convergence results for the Roothaan and Level-Shifting algorithms. In each case, our method of proof provides estimates on the convergence rate. We compare these with numerical results for the algorithms studied. △ Less

Submitted 1 February, 2012; v1 submitted 26 September, 2011; originally announced September 2011.

Forced motion of a probe particle near the colloidal glass transition

Authors: P. Habdas, D. Schaar, Andrew C. Levitt, Eric R. Weeks

Abstract: We use confocal microscopy to study the motion of a magnetic bead in a dense colloidal suspension, near the colloidal glass transition volume fraction $φ_g$. For dense liquid-like samples near $φ_g$, below a threshold force the magnetic bead exhibits only localized caged motion. Above this force, the bead is pulled with a fluctuating velocity. The relationship between force and velocity becomes… ▽ More We use confocal microscopy to study the motion of a magnetic bead in a dense colloidal suspension, near the colloidal glass transition volume fraction $φ_g$. For dense liquid-like samples near $φ_g$, below a threshold force the magnetic bead exhibits only localized caged motion. Above this force, the bead is pulled with a fluctuating velocity. The relationship between force and velocity becomes increasingly nonlinear as $φ_g$ is approached. The threshold force and nonlinear drag force vary strongly with the volume fraction, while the velocity fluctuations do not change near the transition. △ Less

Submitted 19 May, 2004; v1 submitted 29 August, 2003; originally announced August 2003.

Comments: 7 pages, 4 figures revised version, accepted for publication in Europhysics Letters

Journal ref: Europhys. Lett. 67, 477-483 (2004)