Abstract:We present the ElasTool package, an automated toolkit for calculating the second-order elastic constants (SOECs) of any two- (2D) and three-dimensional (3D) crystal systems.ElasTool uses three kinds of strain-matrix sets such as the high-efficiency strain-matrix sets (OHESS), the universal linear-independent coupling strains (ULICS), and the all-single-element strain-matrix sets (ASESS) to calculate the SOECs automatically. ElasTool can efficiently compute both zero- and high-temperature elastic constants. We describe in detail the theoretical background and computational method of elastic constants, the package structure, the installation, and run, the input/output files, the controlling parameters, and two representative examples of how to use the ElasTool package. ElasTool is useful for either the exploration of materials' elastic properties or high-throughput new materials screening and design. ElasTool is freely available on GitHub: https://github.com/elastoolProgram summaryProgram Title: ElasToolCPC Library link to program files: https://doi.org/10.17632/ktvmxrdhpz.1Code Ocean capsule: https://codeocean.com/capsule/1893813Licensing provisions: GNU General Public License, version 3Programming language: Python 3External routines: NumPy [1], Spglib [2], ASE [3], Pandas [4]Nature of problem: The stress-strain method of elastic constants calculation depends on accurate stresses calculated with first-principles methods, such as the density functional theory (DFT). Compared to the energy-strain method, the stress-strain approach needs a smaller number of strain sets to solve the equation sets needed to deduce the elastic constants; it is also more straightforward to implement. However, accurate stresses take a lot of time to compute within DFT. Thus, a smaller number of strain sets and more efficient strain sets are urgently needed to improve the computational efficiency of elastic constants. An automated solution coupled with DFT is necessary for the exploration of materials' elastic properties and high-throughput new materials screening and design.Solution method: The solution to improve the computational efficiency of the stress-strain method is to decrease the number of strain-matrix sets and optimize the strain-matrix sets. We coupled our previously proposed high-efficiency strain-matrix sets (OHESS) with DFT and automated the processes of calculating the elastic tensor using the stress-strain method in the ElasTool package. ElasTool can also adopt the all-single-element strain-matrix sets (ASESS) and the universal linear-independent coupling strains (ULICS) approaches. It can deal with both zero- and high-temperature elastic constants of any crystal systems belonging to 2D or 3D. Having obtained the elastic moduli, ElasTool also gives other essential mechanical and elastic properties of materials such as Young's modulus, bulk modulus, elastic anisotropy, Debye temperature, and the sound velocities.Additional comments including restrictions and unusual features: Currently, this package interfaces with Vienna Ab initio Simulation Package (VASP) code as the stress tensors calculator. [5-7] Extension to other electronic structures is straightforward.References[1]https://numpy.org/[2]https://atztogo.github.io/spglib/[3]https://wiki.fysik.dtu.dk/ase/[4]https://pandas.pydata.org/[5]https://www.vasp.at/[6]G. Kresse, J. Furthmüller, Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set, Phys. Rev. B 54 (1996) 11169.[7]G. Kresse, D. Joubert, From ultrasoft pseudopotentials to the projector augmented-wave method, Phys. Rev. B 59 (1999) 1758.

Accurate stress calculations based on numerical atomic orbital bases: Implementation and benchmarks

Accurate Stress Calculations Based On Numerical Atomic Orbital Bases: And Benchmarks

Introduction to First-Principles Simulation Package Abacus Based on Systematically Improvable Atomic Orbitals

On Stress Calculations in Atomistic Simulations

All-electron periodic G 0 W 0 implementation with numerical atomic orbital basis functions: Algorithm and benchmarks

Numeric Atom-Centered-Orbital Basis Sets With Valence-Correlation Consistency From H To Ar

DFT+ within the framework of linear combination of numerical atomic orbitals

Strategy for constructing compact numerical atomic orbital basis sets by incorporating the gradients of reference wavefunctions

Investigation on Applicability of Various Stress Definitions in Atomistic Simulation

Implementation of the SCAN Exchange-Correlation Functional with Numerical Atomic Orbitals

How to Compute the Atomic Stress Objectively?

ElasTool: An automated toolkit for elastic constants calculation

The pseudoatomic orbital basis: electronic accuracy and soft-mode distortions in ABO$_3$ perovskites

Atomistic Simulation Study of the FCC and BCC Crystal-Melt Interface Stresses

Ab initio Calculations in Atoms, Molecules, and Solids, Treating Spin-Orbit Coupling and Electron Interaction on Equal Footing

Ab Initio Perturbation Calculations of Realistic Effective Interactions in the Hartree–Fock Basis

Implementation of the meta-GGA exchange-correlation functional in numerical atomic orbital basis: With systematic testing on SCAN, rSCAN, and r2SCAN functionals

The Precision Calculations of Na Polarizabilities and Na−affinities

Deformed natural orbitals for ab initio calculations

Adaptive Local Basis Set for Kohn-Sham Density Functional Theory in a Discontinuous Galerkin Framework II: Force, Vibration, and Molecular Dynamics Calculations

Advancing non-atom-centered basis methods for more accurate interaction energies: benchmarks and large-scale applications