Researcher profile

Mi-Song Dupuy

Mi-Song Dupuy contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
6works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

6 published item(s)

preprint2021arXiv

Finite-size effects in response functions of molecular systems

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.

preprint2020arXiv

Inversion symmetry of singular values and a new orbital ordering method in tensor train approximations for quantum chemistry

The tensor train approximation of electronic wave functions lies at the core of the QC-DMRG (Quantum Chemistry Density Matrix Renormalization Group) method, a recent state-of-the-art method for numerically solving the $N$-electron Schrödinger equation. It is well known that the accuracy of TT approximations is governed by the tail of the associated singular values, which in turn strongly depends on the ordering of the one-body basis. Here we find that the singular values $s_1\ge s_2\ge ... \ge s_d$ of tensors representing ground states of noninteracting Hamiltonians possess a surprising inversion symmetry, $s_1s_d=s_2s_{d-1}=s_3s_{d-2}=...$, thus reducing the tail behaviour to a single hidden invariant, which moreover depends explicitly on the ordering of the basis. For correlated wavefunctions, we find that the tail is upper bounded by a suitable superposition of the invariants. Optimizing the invariants or their superposition thus provides a new ordering scheme for QC-DMRG. Numerical tests on simple examples, i.e. linear combinations of a few Slater determinants, show that the new scheme reduces the tail of the singular values by several orders of magnitudes over existing methods, including the widely used Fiedler order.

preprint2020arXiv

Variational projector-augmented wave method: a full-potential approach for electronic structure calculations in solid-state physics

In solid-state physics, energies of crystals are usually computed with a plane-wave discretization of Kohn-Sham equations. However the presence of Coulomb singularities requires the use of large plane-wave cut-offs to produce accurate numerical results. In this paper, an analysis of the plane-wave convergence of the eigenvalues of periodic linear Hamiltonians with Coulomb potentials using the variational projector-augmented wave (VPAW) method is presented. In the VPAW method, an invertible transformation is applied to the original eigenvalue problem, acting locally in balls centered at the singularities. In this setting, a generalized eigenvalue problem needs to be solved using plane-waves. We show that cusps of the eigenfunctions of the VPAW eigenvalue problem at the positions of the nuclei are significantly reduced. These eigenfunctions have however a higher-order derivative discontinuity at the spheres centered at the nuclei. By balancing both sources of error, we show that the VPAW method can drastically improve the plane-wave convergence of the eigenvalues with a minor additional computational cost. Numerical tests are provided confirming the efficiency of the method to treat Coulomb singularities.

preprint2019arXiv

Discretization error cancellation in the plane-wave approximation of periodic Hamiltonians with Coulomb singularities

In solid-state physics, energies of molecular systems are usually computed with a plane-wave discretization of Kohn-Sham equations. A priori estimates of plane-wave convergence for periodic Kohn-Sham calculations with pseudopotentials have been proved , however in most computations in practice, plane-wave cut-offs are not tight enough to target the desired accuracy. It is often advocated that the real quantity of interest is not the value of the energy but of energy differences for different configurations. The computed energy difference is believed to be much more accurate because of `discretization error cancellation', since the sources of numerical errors are essentially the same for different configurations. For periodic linear Hamiltonians with Coulomb potentials, error cancellation can be explained by the universality of the Kato cusp condition. Using weighted Sobolev spaces, Taylor-type expansions of the eigenfunctions are available yielding a precise characterization of this singularity. This then gives an explicit formula of the first order term of the decay of the Fourier coefficients of the eigenfunctions. It enables one to prove that errors on eigenvalue differences are reduced but converge at the same rate as the error on the eigenvalue.

preprint2017arXiv

Projector augmented-wave method: an analysis in a one-dimensional setting

In this article, a numerical analysis of the projector augmented-wave (PAW) method is presented, restricted to the case of dimension one with Dirac potentials modeling the nuclei in a periodic setting. The PAW method is widely used in electronic ab initio calculations, in conjunction with pseudopotentials. It consists in replacing the original electronic Hamiltonian $H$ by a pseudo-Hamiltonian $H^{PAW}$ via the PAW transformation acting in balls around each nuclei. Formally, the new eigenvalue problem has the same eigenvalues as $H$ and smoother eigenfunctions. In practice, the pseudo-Hamiltonian $H^{PAW}$ has to be truncated, introducing an error that is rarely analyzed. In this paper, error estimates on the lowest PAW eigenvalue are proved for the one-dimensional periodic Schrödinger operator with double Dirac potentials.

preprint2017arXiv

Variational projector augmented-wave method: theoretical analysis and preliminary numerical results

In Kohn-Sham electronic structure computations, wave functions have singularities at nuclear positions. Because of these singularities, plane-wave expansions give a poor approximation of the eigenfunctions. In conjunction with the use of pseudo-potentials, the PAW (projector augmented-wave) method circumvents this issue by replacing the original eigenvalue problem by a new one with the same eigenvalues but smoother eigenvectors. Here a slightly different method, called VPAW (variational PAW), is proposed and analyzed. This new method allows for a better convergence with respect to the number of plane-waves. Some numerical tests on an idealized case corroborate this efficiency.