Source author record

Thierry Deutsch

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

15works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

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

15 published item(s)

preprint2022arXiv

Exact solution of the many-body problem with a $\mathcal{O}\left(n^6\right)$ complexity

In this article, we define a new mathematical object, called a pair $D=\left(A,C\right)$ of anti-commutation matrices (ACMP) based on the anti-commutation relation $a^†_{i}a_{j} + a_{j}a^†_{i} = δ_{ij}$ applied to the scalar product between the many-body wavefunctions. This ACMP explicitly separates the different levels of correlation. The one-body correlations are defined by a ACMP $D^0=\left(A^0,C^0\right)$ and the two-body ones by a set of $n$ ACMPs $D^i=\left(A^i,C^i\right)$ where $n$ is the number of states. We show that we can have a compact and exact parametrization with $n^4$ parameters of the two-body reduced density matrix (\TRDM) of any pure or mixed $N$-body state to determine the ground state energy with a $\mathcal{O}\left(n^6\right)$ complexity.

preprint2016arXiv

Challenges in Large Scale Quantum Mechanical Calculations

During the past decades, quantum mechanical methods have undergone an amazing transition from pioneering investigations of experts into a wide range of practical applications, made by a vast community of researchers. First principles calculations of systems containing up to a few hundred atoms have become a standard in many branches of science. The sizes of the systems which can be simulated have increased even further during recent years, and quantum-mechanical calculations of systems up to many thousands of atoms are nowadays possible. This opens up new appealing possibilities, in particular for interdisciplinary work, bridging together communities of different needs and sensibilities. In this review we will present the current status of this topic, and will also give an outlook on the vast multitude of applications, challenges and opportunities stimulated by electronic structure calculations, making this field an important working tool and bringing together researchers of many different domains.

preprint2016arXiv

Revisiting the domain model for lithium intercalated graphite

In this letter, we study the stability of the domain model for lithium intercalated graphite in stages III and II by means of Density Functional Theory and Kinetic Lattice Monte Carlo simulations. We find that the domain model is either thermodynamically or kinetically stable when compared to the standard model in stages III and II. The existence of domains in the intercalation sequence is well supported by recent high resolution transmission electron microscope observations in lithiated graphite. Moreover, we predict that such domain staging sequences leads to a wide range of diffusivity as reported in experiments.

preprint2015arXiv

Accurate and efficient linear scaling DFT calculations with universal applicability

Density Functional Theory calculations traditionally suffer from an inherent cubic scaling with respect to the size of the system, making big calculations extremely expensive. This cubic scaling can be avoided by the use of so-called linear scaling algorithms, which have been developed during the last few decades. In this way it becomes possible to perform ab-initio calculations for several tens of thousands of atoms or even more within a reasonable time frame. However, even though the use of linear scaling algorithms is physically well justified, their implementation often introduces some small errors. Consequently most implementations offering such a linear complexity either yield only a limited accuracy or, if one wants to go beyond this restriction, require a tedious fine tuning of many parameters. In our linear scaling approach within the BigDFT package, we were able to overcome this restriction. Using an ansatz based on localized support functions expressed in an underlying Daubechies wavelet basis -- which offers ideal properties for accurate linear scaling calculations -- we obtain an amazingly high accuracy and a universal applicability while still keeping the possibility of simulating large systems with only a moderate demand of computing resources.

preprint2015arXiv

Fragment Approach to Constrained Density Functional Theory Calculations using Daubechies Wavelets

In a recent paper we presented a linear scaling Kohn-Sham density functional theory (DFT) code based on Daubechies wavelets, where a minimal set of localized support functions is optimized in situ and therefore adapted to the chemical properties of the molecular system. Thanks to the systematically controllable accuracy of the underlying basis set, this approach is able to provide an optimal contracted basis for a given system: accuracies for ground state energies and atomic forces are of the same quality as an uncontracted, cubic scaling approach. This basis set offers, by construction, a natural subset where the density matrix of the system can be projected. In this paper we demonstrate the flexibility of this minimal basis formalism in providing a basis set that can be reused as-is, i.e. without reoptimization, for charge-constrained DFT calculations within a fragment approach. Support functions, represented in the underlying wavelet grid, of the template fragments are roto-translated with high numerical precision to the required positions and used as projectors for the charge weight function. We demonstrate the interest of this approach to express highly precise and efficient calculations for preparing diabatic states and for the computational setup of systems in complex environments.

preprint2015arXiv

Identification of Resonant States via the Generalized Virial Theorem

The numerical extraction of resonant states of open quantum systems is usually a difficult problem. Regularization techniques, such as the mapping to complex coordinates or the addition of Complex Absorbing Potentials are typically employed, as they render resonant wavefunctions localized and therefore normalizable. Physically relevant metastable states have energies that do not depend on the chosen regularization method. Their identification therefore involves cumbersome comparisons between multiple regularised calculations, often performed graphically, which require fine-tuning and specific intuition to avoid approximated, if not wrong, results. In this Letter, we define an operator that explicitly measures such invariance, valid for any arbitrary mapping of spatial coordinates. Resonant states of the system can eventually simply be identified evaluating the expectation value of this operator. Our method eases the extraction of resonant states even for numerical potentials that are difficult to scale to complex coordinates, and avoids the need for ad hoc complex absorbing potentials. We provide explicit evidence of our findings discussing one-dimensional case-studies,also in the presence of external electric fields.

preprint2015arXiv

Multipole-Preserving Quadratures for Discretization of Functions in Real-Space Electronic Structure Calculations

Discretizing an analytic function on a uniform real-space grid is often done via a straightforward collocation method. This is ubiquitous in all areas of computational physics and quantum chemistry. An example in Density Functional Theory (DFT) is given by the external potential or the pseudo-potential describing the interaction between ions and electrons. The accuracy of the collocation method used is therefore very important for the reliability of subsequent treatments like self-consistent field solutions of the electronic structure problems. By construction, the collocation method introduces numerical artifacts typical of real-space treatments, like the so-called egg-box error, that may spoil the numerical stability of the description when the real-space grid is too coarse. As the external potential is an input of the problem, even a highly precise computational treatment cannot cope this inconvenience. We present in this paper a new quadrature scheme that is able to exactly preserve the moments of a given analytic function even for large grid spacings, while reconciling with the traditional collocation method when the grid spacing is small enough. In the context of real-space electronic structure calculations, we show that this method improves considerably the stability of the results for large grid spacings, opening the path towards reliable low-accuracy DFT calculations with reduced number of degrees of freedom.

preprint2014arXiv

Daubechies Wavelets for Linear Scaling Density Functional Theory

We demonstrate that Daubechies wavelets can be used to construct a minimal set of optimized localized contracted basis functions in which the Kohn-Sham orbitals can be represented with an arbitrarily high, controllable precision. Ground state energies and the forces acting on the ions can be calculated in this basis with the same accuracy as if they were calculated directly in a Daubechies wavelets basis, provided that the amplitude of these contracted basis functions is sufficiently small on the surface of the localization region, which is guaranteed by the optimization procedure described in this work. This approach reduces the computational costs of DFT calculations, and can be combined with sparse matrix algebra to obtain linear scaling with respect to the number of electrons in the system. Calculations on systems of 10,000 atoms or more thus become feasible in a systematic basis set with moderate computational resources. Further computational savings can be achieved by exploiting the similarity of the contracted basis functions for closely related environments, e.g. in geometry optimizations or combined calculations of neutral and charged systems.

preprint2013arXiv

Accurate Complex Scaling of Three Dimensional Numerical Potentials

The complex scaling method, which consists in continuing spatial coordinates into the complex plane, is a well-established method that allows to compute resonant eigenfunctions of the time-independent Schroedinger operator. Whenever it is desirable to apply the complex scaling to investigate resonances in physical systems defined on numerical discrete grids, the most direct approach relies on the application of a similarity transformation to the original, unscaled Hamiltonian. We show that such an approach can be conveniently implemented in the Daubechies wavelet basis set, featuring a very promising level of generality, high accuracy, and no need for artificial convergence parameters. Complex scaling of three dimensional numerical potentials can be efficiently and accurately performed. By carrying out an illustrative resonant state computation in the case of a one-dimensional model potential, we then show that our wavelet-based approach may disclose new exciting opportunities in the field of computational non-Hermitian quantum mechanics.

preprint2013arXiv

Norm-conserving pseudopotentials with chemical accuracy compared to all-electron calculations

By adding a non-linear core correction to the well established Dual Space Gaussian type pseudopotentials for the chemical elements up to the third period, we construct improved pseudopotentials for the Perdew Burke Ernzerhof (PBE) functional and demonstrate that they exhibit excellent accuracy. Our benchmarks for the G2-1 test set show average atomization energy errors of only half a kcal/mol. The pseudopotentials also remain highly reliable for high pressure phases of crystalline solids. When supplemented by empirical dispersion corrections the average error in the interaction energy between molecules is also about half a kcal/mol. The accuracy that can be obtained by these pseudopotentials in combination with a systematic basis set is well superior to the accuracy that can be obtained by commonly used medium size Gaussian basis sets in all-electron calculations.

preprint2011arXiv

Real-Time Dynamics and Conical Intersections

This book chapter is a general introduction to photo-dynamics. An outlined structure to the theoretical studies on classical/quantum/mixed dynamics were described from the literature point of view. Such methods get constrained for the practical applications due to parametrizations/computational cost. As an alternate TDDFT-based dynamics studies are promising. we have shown that the TDDFT calculations could be as accurate as high quality monte-carlo calculations and the reliability of the TDDFT calculations for locating the conical intersection through the simple ring-opening reaction of oxirane. Pros and cons of TDDFT for photo-dynamics calculations have been presented for the development of this fast evolving field.

preprint2011arXiv

Wavelet-Based Linear-Response Time-Dependent Density-Functional Theory

Linear-response time-dependent (TD) density-functional theory (DFT) has been implemented in the pseudopotential wavelet-based electronic structure program BigDFT and results are compared against those obtained with the all-electron Gaussian-type orbital program deMon2k for the calculation of electronic absorption spectra of N2 using the TD local density approximation (LDA). The two programs give comparable excitation energies and absorption spectra once suitably extensive basis sets are used. Convergence of LDA density orbitals and orbital energies to the basis-set limit is significantly faster for BigDFT than for deMon2k. However the number of virtual orbitals used in TD-DFT calculations is a parameter in BigDFT, while all virtual orbitals are included in TD-DFT calculations in deMon2k. As a reality check, we report the x-ray crystal structure and the measured and calculated absorption spectrum (excitation energies and oscillator strengths) of the small organic molecule N-cyclohexyl-2-(4-methoxyphenyl)imidazo[1,2-a]pyridin-3-amine.

preprint2011arXiv

Wavelets for Density-Functional Theory and Post-Density-Functional-Theory Calculations

We give a fairly comprehensive review of wavelets and of their application to density-functional theory (DFT) and to our recent application of a wavelet-based version of linear-response time-dependent DFT (LR-TD-DFT). Our intended audience is quantum chemists and theoretical solid-state and chemical physicists. Wavelets are a Fourier-transform-like approach which developed primarily in the latter half of the last century and which was rapidly adapted by engineers in the 1990s because of its advantages compared to standard Fourier transform techniques for multiresolution problems with complicated boundary conditions. High performance computing wavelet codes now also exist for DFT applications in quantum chemistry and solid-state physics, notably the BigDFT code described in this chapter. After briefly describing the basic equations of DFT and LR-TD-DFT, we discuss how they are solved in BigDFT and present new results on the small test molecule carbon monoxide to show how BigDFT results compare against those obtained with the quantum chemistry gaussian-type orbital (GTO) based code deMon2k. In general, the two programs give essentially the same orbital energies, but the wavelet basis of BigDFT converges to the basis set limit much more rapidly than does the GTO basis set of deMon2k. Wavelet-based LR-TD-DFT is still in its infancy, but our calculations confirm the feasibility of implementing LR-TD-DFT in a wavelet-based code.

preprint2010arXiv

Structural Metastability of Endohedral Silicon Fullerenes

Endohedrally doped Si20 fullerenes appear as appealing building blocks for nanoscale materials. We investigate their structural stability with an unbiased and systematic global geometry optimization method within density-functional theory. For a wide range of metal doping atoms, it was sufficient to explore the Born Oppenheimer surface for only a moderate number of local minima to find structures that clearly differ from the initial endohedral cages, but are considerably more favorable in terms of energy. Previously proposed structures are thus all metastable.

preprint2008arXiv

Daubechies wavelets as a basis set for density functional pseudopotential calculations

Daubechies wavelets are a powerful systematic basis set for electronic structure calculations because they are orthogonal and localized both in real and Fourier space. We describe in detail how this basis set can be used to obtain a highly efficient and accurate method for density functional electronic structure calculations. An implementation of this method is available in the ABINIT free software package. This code shows high systematic convergence properties, very good performances and an excellent efficiency for parallel calculations.