Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userSolving the Bethe-Salpeter equation with exponential convergence
95
0
0.0
(
0
)
Added by
Markus Wallerberger
Publication date
2020
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
The Bethe-Salpeter equation plays a crucial role in understanding the physics of correlated fermions, relating to optical excitations in solids as well as resonances in high-energy physics. Yet, it is notoriously difficult to control numerically, typically requiring an effort that scales polynomially with energy scales and accuracy. This puts many interesting systems out of computational reach. Using the intermediate representation and sparse modelling for two-particle objects on the Matsubara axis, we develop an algorithm that solves the Bethe-Salpeter equation in $O(L^8)$ time with $O(L^4)$ memory, where $L$ grows only logarithmically with inverse temperature, bandwidth, and desired accuracy, This opens the door for computations in hitherto inaccessible regimes. We benchmark the method on the Hubbard atom and on the multi-orbital weak-coupling limit, where we observe the expected exponential convergence to the analytical results. We then showcase the method for a realistic impurity problem.
rate research
Read More
We present a method to compute optical spectra and exciton binding energies of molecules and solids based on the solution of the Bethe-Salpeter equation (BSE) and the calculation of the screened Coulomb interaction in finite field. The method does not require the explicit evaluation of dielectric matrices nor of virtual electronic states, and can be easily applied without resorting to the random phase approximation. In addition it utilizes localized orbitals obtained from Bloch states using bisection techniques, thus greatly reducing the complexity of the calculation and enabling the efficient use of hybrid functionals to obtain single particle wavefunctions. We report exciton binding energies of several molecules and absorption spectra of condensed systems of unprecedented size, including water and ice samples with hundreds of atoms.
The last ten years have witnessed fast spreading of massively parallel computing clusters, from leading supercomputing facilities down to the average university computing center. Many companies in the private sector have undergone a similar evolution. In this scenario, the seamless integration of software and middleware libraries is a key ingredient to ensure portability of scientific codes and guarantees them an extended lifetime. In this work, we describe the integration of the ChASE library, a modern parallel eigensolver, into an existing legacy code for the first-principles computation of optical properties of materials via solution of the Bethe-Salpeter equation for the optical polarization function. Our numerical tests show that, as a result of integrating ChASE and parallelizing the reading routine, the code experiences a remarkable speedup and greatly improved scaling behavior on both multi- and many-core architectures. We demonstrate that such a modernized BSE code will, by fully exploiting parallel computing architectures and file systems, enable domain scientists to accurately study complex material systems that were not accessible before.
We present a method to directly solving the Bethe-Salpeter equation in Minkowski space, both for bound and scattering states. It is based on a proper treatment of the singularities which appear in the kernel, propagators and Bethe-Salpeter amplitude itself. The off-mass shell scattering amplitude for spinless particles interacting by a one boson exchange is computed for the first time.
The method of solving the Bethe-Salpeter equation in Minkowski space, which we developed previously for spinless particles, is extended to a system of two fermions. The method is based on the Nakanishi integral representation of the amplitude and on projecting the equation on the light-front plane. The singularities in the projected two-fermion kernel are regularized without modifying the original Bethe-Salpeter amplitudes. The numerical solutions for the J=0 bound state with the scalar, pseudoscalar and massless vector exchange kernels are found. The stability of the scalar and positronium states without vertex form factor is discussed. Binding energies are in close agreement with the Euclidean results. Corresponding amplitudes in Minkowski space are obtained.
The off-mass shell scattering amplitude, satisfying the Bethe-Salpeter equation for spinless particles in Minkowski space with the ladder kernel, is computed for the first time.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
