Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userOptimization Algorithm for the Generation of ONCV Pseudopotentials
147
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We present an optimization algorithm to construct pseudopotentials and use it to generate a set of Optimized Norm-Conserving Vanderbilt (ONCV) pseudopotentials for elements up to Z=83 (Bi) (excluding Lanthanides). We introduce a quality function that assesses the agreement of a pseudopotential calculation with all-electron FLAPW results, and the necessary plane-wave energy cutoff. This quality function allows us to use a Nelder-Mead optimization algorithm on a training set of materials to optimize the input parameters of the pseudopotential construction for most of the periodic table. We control the accuracy of the resulting pseudopotentials on a test set of materials independent of the training set. We find that the automatically constructed pseudopotentials provide a good agreement with the all-electron results obtained using the FLEUR code with a plane-wave energy cutoff of approximately 60 Ry.
rate research
Read More
We show that efficient norm-conserving pseudopotentials for electronic structure calculations can be obtained from a polynomial Ansatz for the potential. Our pseudopotential is a polynomial of degree ten in the radial variable and fulfills the same smoothness conditions imposed by the Troullier-Martins method [Phys. Rev. B 43, 1993 (1991)] where pseudopotentials are represented by a polynomial of degree twenty-two. We compare our method to the Troullier-Martins approach in electronic structure calculations for diamond and iron in the bcc structure and find that the two methods perform equally well in calculations of the total energy. However, first and second derivatives of the total energy with respect to atomic coordinates converge significantly faster with the plane wave cutoff if the standard Troullier-Martins potentials are replaced by the pseudopotentials introduced here.
Strong multiple scattering of the probe in scanning transmission electron microscopy (STEM) means image simulations are usually required for quantitative interpretation and analysis of elemental maps produced by electron energy-loss spectroscopy (EELS). These simulations require a full quantum-mechanical treatment of multiple scattering of the electron beam, both before and after a core-level inelastic transition. Current algorithms scale quadratically and can take up to a week to calculate on desktop machines even for simple crystal unit cells and do not scale well to the nano-scale heterogeneous systems that are often of interest to materials science researchers. We introduce an algorithm with linear scaling that typically results in an order of magnitude reduction in compute time for these calculations without introducing additional error and discuss approximations that further improve computational scaling for larger scale objects with modest penalties in calculation error. We demonstrate these speed-ups by calculating the atomic resolution STEM-EELS map using the L-edge transition of Fe, for of a nanoparticle 80 AA in diameter in 16 hours, a calculation that would have taken at least 80 days using a conventional multislice approach.
We report exact expressions for atomic forces in the diffusion Monte Carlo (DMC) method when using nonlocal pseudopotentials. We present approximate schemes for estimating these expressions in both mixed and pure DMC calculations, including the pseudopotential Pulay term which has not previously been calculated and the Pulay nodal term which has not been calculated for real systems in pure DMC simulations. Harmonic vibrational frequencies and equilibrium bond lengths are derived from the DMC forces and compared with those obtained from DMC potential energy curves. Results for four small molecules show that the equilibrium bond lengths obtained from our best force and energy calculations differ by less than 0.002 Angstrom.
We extend the recently developed converse NMR approach [T. Thonhauser, D. Ceresoli, A. Mostofi, N. Marzari, R. Resta, and D. Vanderbilt, J. Chem. Phys. textbf{131}, 101101 (2009)] such that it can be used in conjunction with norm-conserving, non-local pseudopotentials. This extension permits the efficient ab-initio calculation of NMR chemical shifts for elements other than hydrogen within the convenience of a plane-wave pseudopotential approach. We have tested our approach on several finite and periodic systems, finding very good agreement with established methods and experimental results.
Fully-nonlocal two-projector norm-conserving pseudopotentials are shown to be compatible with a systematic approach to the optimization of convergence with the size of the plane-wave basis. A new formulation of the optimization is developed, including the ability to apply it to positive-energy atomic scattering states, and to enforce greater continuity in the pseudopotential. The generalization of norm-conservation to multiple projectors is reviewed and recast for the present purposes. Comparisons among the results of all-electron and one- and two-projector norm-conserving pseudopotential calculations of lattice constants and bulk moduli are made for a group of solids chosen to represent a variety of types of bonding and a sampling of the periodic table.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
