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.

[1]  A. Levitt Screening in the Finite-Temperature Reduced Hartree–Fock Model , 2018, Archive for Rational Mechanics and Analysis.

[2]  S. Teufel Non-equilibrium Almost-Stationary States and Linear Response for Gapped Quantum Systems , 2017, Communications in Mathematical Physics.

[3]  P. Stefanov,et al.  Approximating Resonances with the Complex Absorbing Potential Method , 2004, math-ph/0409020.

[4]  M. Radosz The principles of limit absorption and limit amplitude for periodic operators , 2010 .

[5]  D. Gontier,et al.  Supercell calculations in the reduced Hartree-Fock model for crystals with local defects , 2015, 1512.08636.

[6]  D. Fujiwara A construction of the fundamental solution for the Schrödinger equation , 1979 .

[7]  E. Cancès,et al.  Coherent Electronic Transport in Periodic Crystals , 2020, Annales Henri Poincaré.

[8]  A. Klein,et al.  Linear response theory for magnetic Schrödinger operators in disordered media , 2004, math-ph/0408028.

[9]  G. Stoltz,et al.  A mathematical formulation of the random phase approximation for crystals , 2011, 1109.2416.

[10]  Shmuel Agmon,et al.  Spectral properties of Schrödinger operators and scattering theory , 1975 .

[11]  Shmuel Agmon,et al.  Analyticity properties in scattering and spectral theory for Schrödinger operators with long-range radial potentials , 1992 .

[12]  Luigi Genovese,et al.  Locality and computational reliability of linear response calculations for molecular systems , 2018, Physical Review Materials.

[13]  Kristian Kirsch,et al.  Methods Of Modern Mathematical Physics , 2016 .

[14]  Virginie Ehrlacher,et al.  Numerical quadrature in the Brillouin zone for periodic Schrödinger operators , 2020, Numerische Mathematik.

[15]  Yvon Maday,et al.  Numerical analysis of the planewave discretization of some orbital-free and Kohn-Sham models , 2010, 1003.1612.

[16]  Yvon Maday,et al.  Non-consistent approximations of self-adjoint eigenproblems: application to the supercell method , 2012, Numerische Mathematik.

[17]  M. Fraas,et al.  The Adiabatic Theorem and Linear Response Theory for Extended Quantum Systems , 2017, 1705.02838.

[18]  J. G. Muga,et al.  Complex absorbing potentials , 2004 .

[19]  E. Prodan Quantum Transport in Disordered Systems Under Magnetic Fields: A Study Based on Operator Algebras , 2012, 1204.6490.

[20]  Xavier Antoine,et al.  A friendly review of absorbing boundary conditions and perfectly matched layers for classical and relativistic quantum waves equations , 2017 .