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.

Deformcell: a Python Script to Simplify and Fasten Mechanical Properties Calculations of Molecular Crystals in VASP Package for Research and Teaching Purposes

Computational crystal plastiticy : from single crystal to homogenized polycrystals

MyElas: An automatized tool-kit for high-throughput calculation, post-processing and visualization of elasticity and related properties of solids

VASPKIT: A user-friendly interface facilitating high-throughput computing and analysis using VASP code

doped: Python toolkit for robust and repeatable charged defect supercell calculations

Modeling a Unit Cell: Crystallographic Refinement Procedure Using the Biomolecular MD Simulation Platform Amber

A toolbox for improving the workflow of NMR crystallography

X-ray Crystallography Module in MD Simulation Program Amber 2023. Refining the Models of Protein Crystals

3DStructGen: an interactive web-based 3D structure generation for non-periodic molecule and crystal

Nudged-Elastic Band Calculations of Polymorph Transitions and Solid-State Reactions in Molecular Crystals

FaVAD: A software workflow for characterisation and visualizing of defects in crystalline structures

qvasp: A flexible toolkit for VASP users in materials simulations

FlεX: a computer vision program to evaluate strain in flexible crystals

MolDy: molecular dynamics simulation made easy

A method to extract slip system dependent information for crystal plasticity models

Improved Description of the Structure of Molecular and Layered Crystals: Ab Initio DFT Calculations with van der Waals Corrections

VASPKIT: A Pre- and Post-Processing Program for VASP code

ElasTool: An automated toolkit for elastic constants calculation

MD-assisted Refinement of X-Ray Coordinates

"On‐the‐fly" Crystal : How to reliably and automatically characterize and construct potential energy surfaces

MechElastic: A Python Library for Analysis of Mechanical and Elastic Properties of Bulk and 2D Materials