Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userSolving the Bethe-Salpeter equation on massively parallel architectures
115
0
0.0
(
0
)
Added by
Edoardo di Napoli
Publication date
2020
fields
Informatics Engineering
and research's language is
English
Authors
Xiao Zhang
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
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.
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 molecular dynamics simulation code ls1 mardyn is presented. It is a highly scalable code, optimized for massively parallel execution on supercomputing architectures, and currently holds the world record for the largest molecular simulation with over four trillion particles. It enables the application of pair potentials to length and time scales which were previously out of scope for molecular dynamics simulation. With an efficient dynamic load balancing scheme, it delivers high scalability even for challenging heterogeneous configurations. Presently, multi-center rigid potential models based on Lennard-Jones sites, point charges and higher-order polarities are supported. Due to its modular design, ls1 mardyn can be extended to new physical models, methods, and algorithms, allowing future users to tailor it to suit their respective needs. Possible applications include scenarios with complex geometries, e.g. for fluids at interfaces, as well as non-equilibrium molecular dynamics simulation of heat and mass transfer.
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.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
