Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userTwo-level Chebyshev filter based complementary subspace method: pushing the envelope of large-scale electronic structure calculations
100
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We describe a novel iterative strategy for Kohn-Sham density functional theory calculations aimed at large systems (> 1000 electrons), applicable to metals and insulators alike. In lieu of explicit diagonalization of the Kohn-Sham Hamiltonian on every self-consistent field (SCF) iteration, we employ a two-level Chebyshev polynomial filter based complementary subspace strategy to: 1) compute a set of vectors that span the occupied subspace of the Hamiltonian; 2) reduce subspace diagonalization to just partially occupied states; and 3) obtain those states in an efficient, scalable manner via an inner Chebyshev-filter iteration. By reducing the necessary computation to just partially occupied states, and obtaining these through an inner Chebyshev iteration, our approach reduces the cost of large metallic calculations significantly, while eliminating subspace diagonalization for insulating systems altogether. We describe the implementation of the method within the framework of the Discontinuous Galerkin (DG) electronic structure method and show that this results in a computational scheme that can effectively tackle bulk and nano systems containing tens of thousands of electrons, with chemical accuracy, within a few minutes or less of wall clock time per SCF iteration on large-scale computing platforms. We anticipate that our method will be instrumental in pushing the envelope of large-scale ab initio molecular dynamics. As a demonstration of this, we simulate a bulk silicon system containing 8,000 atoms at finite temperature, and obtain an average SCF step wall time of 51 seconds on 34,560 processors; thus allowing us to carry out 1.0 ps of ab initio molecular dynamics in approximately 28 hours (of wall time).
rate research
Read More
The GW method is a many-body electronic structure technique capable of generating accurate quasiparticle properties for realistic systems spanning physics, chemistry, and materials science. Despite its power, GW is not routinely applied to large complex assemblies due to its large computational overhead and quartic scaling with particle number. Here, the GW equations are recast, exactly, as Fourier-Laplace time integrals over complex time propagators. The propagators are then shredded via energy partitioning and the time integrals approximated in a controlled manner using generalized Gaussian quadrature(s) while discrete variable methods are employed to represent the required propagators in real-space. The resulting cubic scaling GW method has a sufficiently small prefactor to outperform standard quartic scaling methods on small systems ($gtrapprox$ 10 atoms) and also represents a substantial improvement over other cubic methods tested for all system sizes studied. The approach can be applied to any theoretical framework containing large sums of terms with energy differences in the denominator.
We propose an adaptive planewave method for eigenvalue problems in electronic structure calculations. The method combines a priori convergence rates and accurate a posteriori error estimates into an effective way of updating the energy cut-off for planewave discretizations, for both linear and nonlinear eigenvalue problems. The method is error controllable for linear eigenvalue problems in the sense that for a given required accuracy, an energy cut-off for which the solution matches the target accuracy can be reached efficiently. Further, the method is particularly promising for nonlinear eigenvalue problems in electronic structure calculations as it shall reduce the cost of early iterations in self-consistent algorithms. We present some numerical experiments for both linear and nonlinear eigenvalue problems. In particular, we provide electronic structure calculations for some insulator and metallic systems simulated with Kohn--Sham density functional theory (DFT) and the projector augmented wave (PAW) method, illustrating the efficiency and potential of the algorithm.
Motivated by the recently proposed parallel orbital-updating approach in real space method, we propose a parallel orbital-updating based plane-wave basis method for electronic structure calculations, for solving the corresponding eigenvalue problems. In addition, we propose two new modified parallel orbital-updating methods. Compared to the traditional plane-wave methods, our methods allow for two-level parallelization, which is particularly interesting for large scale parallelization. Numerical experiments show that these new methods are more reliable and efficient for large scale calculations on modern supercomputers
A parameterized tight-binding (TB) model based on the first-principles GW calculations is developed for single layer tin diselenide (SnSe$_2$) and used to study its electronic and optical properties under external magnetic field. The truncated model is derived from six maximally localized wannier orbitals on Se site, which accurately describes the quasi-particle electronic states of single layer SnSe$_2$ in a wide energy range. The quasi-particle electronic states are dominated by the hoppings between nearest wannier orbitals ($t_1$-$t_6$). Our numerical calculation shows that, due to the electron-hole asymmetry, two sets of Landau Level spectrum are obtained when a perpendicular magnetic field is applied. The Landau Level spectrum follows linear dependence on the level index and magnetic field, exhibiting properties of two-dimensional electron gas in traditional semiconductors. The optical conductivity calculation shows that the optical gap is very close to the GW value, and can be tuned by external magnetic field. Our proposed TB model can be used for further exploring the electronic, optical, and transport properties of SnSe$_2$, especially in the presence of external magnetic fields.
Using large-scale DFT calculations, we have investigated the structural and electronic properties of both armchair and zigzag graphdiyne nanotubes as a function of size. To provide insight in these properties, we present new detailed calculations of the structural relaxation energy, effective electron/hole mass, and size-scaling of the bandgap as a function of size and chirality using accurate screened-exchange DFT calculations. These calculations provide a systematic evaluation of the structural and electronic properties of the largest graphdiyne nanotubes to date - up to 1,296 atoms and 23,328 basis functions. Our calculations find that zigzag graphdiyne nanotubes (GDNTs) are structurally more stable compared to armchair GDNTs of the same size. Furthermore, these large-scale calculations allow us to present simple analytical formulae to guide future experimental efforts for estimating the fundamental bandgaps of these unique nanotubes as a function of chirality and diameter. While the bandgaps for both the armchair and zigzag GDNTs can be tuned as a function of size, the conductivity in each of these two different chiralities is markedly different. Zigzag GDNTs have wider valence and conduction bands and are expected to have a higher electron- and hole-mobility than their armchair counterparts.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
