Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userRESCU: a Real Space Electronic Structure Method
397
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
In this work we present RESCU, a powerful MATLAB-based Kohn-Sham density functional theory (KS-DFT) solver. We demonstrate that RESCU can compute the electronic structure properties of systems comprising many thousands of atoms using modest computer resources, e.g. 16 to 256 cores. Its computational efficiency is achieved from exploiting four routes. First, we use numerical atomic orbital (NAO) techniques to efficiently generate a good quality initial subspace which is crucially required by Chebyshev filtering methods. Second, we exploit the fact that only a subspace spanning the occupied Kohn-Sham states is required, and solving accurately the KS equation using eigensolvers can generally be avoided. Third, by judiciously analyzing and optimizing various parts of the procedure in RESCU, we delay the $O(N^3)$ scaling to large $N$, and our tests show that RESCU scales consistently as $O(N^{2.3})$ from a few hundred atoms to more than 5,000 atoms when using a real space grid discretization. The scaling is better or comparable in a NAO basis up to the 14,000 atoms level. Fourth, we exploit various numerical algorithms and, in particular, we introduce a partial Rayleigh-Ritz algorithm to achieve efficiency gains for systems comprising more than 10,000 electrons. We demonstrate the power of RESCU in solving KS-DFT problems using many examples running on 16, 64 and/or 256 cores: a 5,832 Si atoms supercell; a 8,788 Al atoms supercell; a 5,324 Cu atoms supercell and a small DNA molecule submerged in 1,713 water molecules for a total 5,399 atoms. The KS-DFT is entirely converged in a few hours in all cases. Our results suggest that the RESCU method has reached a milestone of solving thousands of atoms by KS-DFT on a modest computer cluster.
rate research
Read More
A characteristic feature of the state-of-the-art of real-space methods in electronic structure calculations is the diversity of the techniques used in the discretization of the relevant partial differential equations. In this context, the main approaches include finite-difference methods, various types of finite-elements and wavelets. This paper reports on the results of several code development projects that approach problems related to the electronic structure using these three different discretization methods. We review the ideas behind these methods, give examples of their applications, and discuss their similarities and differences.
Electronic states are responsible for most material properties, including chemical bonds, electrical and thermal conductivity, as well as optical and magnetic properties. Experimentally, however, they remain mostly elusive. Here, we report the real-space mapping of selected transitions between p and d states on the {AA}ngstrom scale in bulk rutile (TiO2) using electron energy-loss spectrometry (EELS), revealing information on individual bonds between atoms. On the one hand, this enables the experimental verification of theoretical predictions about electronic states. On the other hand, it paves the way for directly investigating electronic states under conditions that are at the limit of the current capabilities of numerical simulations such as, e.g., the electronic states at defects, interfaces, and quantum dots.
Magnetic topological materials, in which the time-reversal symmetry is broken, host various exotic quantum phenomena, including the quantum anomalous Hall effect, axion insulator states, and Majorana fermions. The study of magnetic topological materials is at the forefront of condensed matter physics. Recently, a variety of magnetic topological materials have been reported, such as Mn$_3$Sn, Co$_3$Sn$_2$S$_2$, Fe$_3$Sn$_2$, and MnBi$_2$Te$_4$. Here, we report the observation of a topological electronic structure in an antiferromagnet, HoSbTe, a member of the ZrSiS family of materials, by angle-resolved photoemission spectroscopy measurements and first-principles calculations. We demonstrate that HoSbTe is a Dirac nodal line semimetal when spin-orbit coupling (SOC) is neglected. However, our theoretical calculations show that the strong SOC in HoSbTe fully gaps out the nodal lines and drives the system to a weak topological insulator state, with each layer being a two-dimensional topological insulator. Because of the strong SOC in HoSbTe, the gap is as large as hundreds of meV along specific directions, which is directly observed by our ARPES measurements. The existence of magnetic order and topological properties in HoSbTe makes it a promising material for realization of exotic quantum devices.
First-principles electronic structure calculations are very widely used thanks to the many successful software packages available. Their traditional coding paradigm is monolithic, i.e., regardless of how modular its internal structure may be, the code is built independently from others, from the compiler up, with the exception of linear-algebra and message-passing libraries. This model has been quite successful for decades. The rapid progress in methodology, however, has resulted in an ever increasing complexity of those programs, which implies a growing amount of replication in coding and in the recurrent re-engineering needed to adapt to evolving hardware architecture. The Electronic Structure Library (esl) was initiated by CECAM (European Centre for Atomic and Molecular Calculations) to catalyze a paradigm shift away from the monolithic model and promote modularization, with the ambition to extract common tasks from electronic structure programs and redesign them as free, open-source libraries. They include heavy-duty ones with a high degree of parallelisation, and potential for adaptation to novel hardware within them, thereby separating the sophisticated computer science aspects of performance optimization and re-engineering from the computational science done by scientists when implementing new ideas. It is a community effort, undertaken by developers of various successful codes, now facing the challenges arising in the new model. This modular paradigm will improve overall coding efficiency and enable specialists (computer scientists or computational scientists) to use their skills more effectively. It will lead to a more sustainable and dynamic evolution of software as well as lower barriers to entry for new developers.
We present a real-space adaptive-coordinate method, which combines the advantages of the finite-difference approach with the accuracy and flexibility of the adaptive coordinate method. The discretized Kohn-Sham equations are written in generalized curvilinear coordinates and solved self-consistently by means of an iterative approach. The Poisson equation is solved in real space using the Multigrid algorithm. We implemented the method on a massively parallel computer, and applied it to the calculation of the equilibrium geometry and harmonic vibrational frequencies of the CO_2, CO, N_2 and F_2 molecules, yielding excellent agreement with the results of accurate quantum chemistry and Local Density Functional calculations.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
