PaperlayerView exact revision on arXiv

Source-linked AI summary

ElATools: A tool for analyzing anisotropic elastic properties of the 2D and 3D materials

Shahram Yalameha, Zahra Nourbakhsh, Daryoosh Vashaee

arXiv:2105.07279v3cond-mat.mtrl-sci

TL;DR

Second-order elastic constants have been measured for only a small fraction of crystalline materials, and this limited data availability constrains materials design. ElATools computes mechanical properties, visualizes elastic-property surfaces and projections, and supports phase and group velocity and power-flow visualization in 3D materials. In ZnAu2(CN)4, negative linear compressibility is evident along the z-direction, consistent with the prediction in Ref..

  • Problem

    Second-order elastic constants have been measured for only a small fraction of crystalline materials, and this limited data availability constrains materials design.

    Source: §1 · p. 2

  • Method

    ElATools computes mechanical properties, visualizes elastic-property surfaces and projections, and supports phase and group velocity and power-flow visualization in 3D materials.

    Source: §1 · pp. 3, 5

  • Results

    In ZnAu2(CN)4, negative linear compressibility is evident along the z-direction, consistent with the prediction in Ref..

    Source: §3.4 · p. 18

  • Takeaways & Limitations

    Analysis and visualization of elastic tensors can reveal negative linear compressibility, negative Poisson’s ratio, and highly anisotropic moduli relevant to nanoscale electromechanical devices.

    Source: §1 · p. 2

Abstract

from arXiv · show

We introduce a computational method and a user-friendly code with a terminal-based graphical user interface (GUI), named ElATools, developed to analyze mechanical and anisotropic elastic properties. ElATools enables facile analysis of the second-order elastic stiffness tensor of two-dimensional (2D) and three-dimensional (3D) crystal systems. It computes and displays the main mechanical properties including the bulk modulus, Young's modulus, shear modulus, hardness, p-wave modulus, universal anisotropy index, Chung-Buessem anisotropy index, log-Euclidean anisotropy parameter, Cauchy pressures, Poisson's ratio, and Pugh's ratio, using three averaging schemes of Voigt, Reuss, and Hill. It includes an online and offline database from the Materials Project with more than 13,000 elastic stiffness constants for 3D materials. The program supports output files of the well-known computational codes IRelast, IRelast2D, ElaStic, and AELAS. Four types of plotting and visualization tools are integrated to conveniently interface with GNUPLOT, XMGRACE, view3dscene, and plotly libraries, offering immediate post-processing of the results. It can efficiently identify anomalous mechanical properties, such as negative linear compressibility, negative Poisson's ratio, and highly-anisotropic elastic modulus in 2D and 3D materials, which are central properties to design and develop high-performance nanoscale electromechanical devices. Six case studies on selected material systems, namely, ZnAu$_2$ (CN)$_4$, CrB$_2$, $\delta$-phosphorene, Pd$_2$O$_6$Se$_2$ monolayer, and GaAs, and a hypothetical set of systems with cubic symmetry are presented to demonstrate the descriptive and predictive capabilities of ElATools.

1. INTRODUCTION

Second-order elastic constants are important across engineering and scientific fields, but they have been measured for only a small fraction of crystalline materials. ElATools brings elastic-tensor analysis, visualization, and related capabilities together while adding features for studying 2D and 3D materials.

1. INTRODUCTION

  • Sparse experimental measurements limit materials design, while elastic-tensor analysis and visualization can reveal negative linear compressibility, negative Poisson’s ratio, and highly anisotropic moduli.The paper identifies these anomalous properties as relevant to developing nanoscale electromechanical devices.Source: §1 · pp. 2–3
  • The paper motivates analysis of anisotropic elasticity by noting that elastic constants also contain information about acoustic-wave behavior.The introduction presents phase- and group-velocity visualization as an additional capability for 3D materials.Source: §1 · pp. 2, 5
  • ElATools combines elastic-tensor analysis features in one program and adds capabilities beyond existing packages such as ELAM, ELATE, and MechElastic.The introduction describes these packages and frames ElATools as addressing their shortcomings.Source: §1 · p. 3
  • ElATools computes mechanical properties and anisotropy indices using Voigt, Reuss, and Hill averaging, and assesses mechanical stability from elastic-tensor eigenvalues.The listed properties include Young’s and shear moduli, hardness, Cauchy pressure, Poisson’s ratio, and Pugh’s ratio.Source: §1 · pp. 3, 5
  • Its visualization features include elastic-property surfaces and projections on desired planes, 3D phase and group velocities, and 2D polar covers for selected properties.The listed 2D polar covers show Poisson’s ratio, shear modulus, and Young’s modulus.Source: §1 · pp. 5–6
  • ElATools offers an online and offline Materials Project database of more than 13,000 elastic tensors, multiple visualization output formats, and a terminal-based graphical interface.Output formats support tools including view3dscene, XMGRACE, GNUPLOT, and web browsers.Source: §1 · p. 6

2. THEORETICAL BACKGROUND

The background establishes how elastic tensors encode stress–strain response and how tensor transformations and wave equations yield directional mechanical and acoustic properties. It then describes averaging schemes and anisotropy measures for 2D and 3D materials.

2.1. Hooke’s law and elastic tensor of crystals

  • In the elastic regime, strain is proportional to stress; tensor symmetries reduce 81 stiffness or compliance components to 21 independent coordinates in the least-symmetric case.Further crystal-symmetry reduction gives 21, 15, 9, 7, 5, 5, and 3 independent constants for triclinic through cubic systems, respectively.Source: §2.1 · pp. 6–7

2.2. Transformation law, Christoffel equation, and representation surfaces of elastic properties

  • Tensor transformation laws and direction vectors support calculating orientation-dependent Young’s modulus, linear compressibility, bulk modulus, shear modulus, and Poisson’s ratio.The representation surfaces describe these elastic properties as functions of direction.Source: §2.2 · pp. 8–9
  • For 2D materials, directional Young’s and shear moduli and Poisson’s ratio are derived from compliance constants, with formulas adapted to crystal symmetry.The formulas apply to hexagonal, square, and rectangular systems; oblique systems have six independent elastic constants.Source: §2.2 · p. 10
  • Solving the Christoffel eigenvalue problem yields three wave frequencies and polarization vectors for each wave-vector direction, from which sound velocities can be calculated.The reduced formulation gives one primary and two secondary velocities, corresponding to longitudinal and transverse polarizations.Source: §2.2 · p. 11
  • Group velocity describes wave-packet travel and acoustic-energy transport, and its direction can differ from the phase-velocity direction.The angle between their directions is the power flow angle.Source: §2.2 · p. 12

2.3. Mechanical properties and elastic anisotropy of 2D and 3D materials

  • Voigt and Reuss bounds provide averages of single-crystal elastic constants, while their arithmetic mean is the Voigt-Reuss-Hill estimate for polycrystal behavior.The section gives formulas for the bulk and shear moduli under both bounds.Source: §2.3 · p. 12
  • The universal anisotropy index is zero for isotropic crystals, while a Pugh ratio below 0.57 indicates ductile behavior and higher values indicate brittleness.Cauchy pressure is also presented as a measure of ductility and brittleness, with its sign associated with bonding character.Source: §2.3 · pp. 12–13
  • For 2D systems, the relations among Young’s modulus, shear modulus, Poisson’s ratio, and elastic stiffness constants are derived.Young’s modulus is defined along an axis, Poisson’s ratio relates transverse response to tensile strain, and Gxy is the shear modulus in the xy-plane.Source: §2.3 · p. 13

3. SOFTWARE DESCRIPTION AND FEATURES

ElATools processes elastic stiffness tensors through stability checks and mechanical-property calculations, then saves data for visualization. Its post-processing and case studies demonstrate orientation-dependent properties and anisotropy across 2D and 3D materials.

3.1. Workflow and structure of ElATools

  • The workflow selects a 2D or 3D branch, reads elastic constants from supported files or a Materials Project database of more than 13,000 3D tensors, and checks stability before calculation.The stability check uses positive eigenvalues; tensors that fail it stop with a mechanical-instability error.Source: §3.1 · pp. 13–14
  • The calculation branches produce orientation-dependent and spatial properties, averaged polycrystalline moduli, anisotropy measures, and optional elastic-wave velocities before saving results for post-processing.Four visualization programs generate 2D or 3D files for GNUPLOT, XMGRACE, VRML-capable browsers, and web browsers.Source: §3.1 · pp. 14–15

3.4. Visualization and Post-processing

  • ZnAu2(CN)4 shows negative linear compressibility along z and a small Poisson’s ratio of -0.02 on the (110) plane.Its angular heat maps also make negative Poisson’s ratio visible.Source: §3.4 · p. 18
  • CrB2’s bulk and Young’s moduli are anisotropic, with (010)-plane projections showing their orientation dependence.The paper contrasts these shapes with the sphere, circle, and straight line expected for an isotropic system.Source: §3.4 · p. 19
  • GaAs phase- and group-velocity plots distinguish fast and slow secondary modes by phase velocity, while primary modes show the greatest anisotropy.The group velocity of the fast-secondary mode can be lower than that of the slow-secondary mode for some propagation directions.Source: §3.4 · p. 19
  • δ-phosphorene reaches a maximum negative Poisson’s ratio of -0.267 at 90 degrees, and ElATools also calculates its shear modulus.The reported Young’s-modulus orientation dependence agrees with the cited work.Source: §3.4 · p. 19
  • Pd2O6Se2 monolayer has a Poisson’s ratio of -0.49, visible in the polar heat maps of its elastic properties.The authors describe polar heat maps as useful for finding small negative Poisson’s ratios in 2D materials.Source: §3.4 · p. 19
  • For hypothetical cubic systems with coefficients spanning 0–100 GPa, C44 has a more critical role in negative Poisson’s ratio than C11 or C12.The minimum Poisson’s ratio changes from almost positive at C44 = 1 GPa to negative as C44 rises; it can reach -3 at 100 GPa, while the maximum can reach 4.Source: §3.4 · pp. 20, 22

4. SUMMARY AND OUTLOOK

ElATools analyzes second-order elastic tensors for 2D and 3D crystals and supplies elastic-property and anisotropy measures. Its appendices describe the quantities, estimation methods, and planned extensions.

4. SUMMARY AND OUTLOOK

  • ElATools is a Fortran90 code for analyzing second-order elastic tensors in 2D and 3D crystal systems.The authors also report that it can generate data for machine learning to detect and predict elastic and anisotropy properties.Source: §4 · p. 25
  • The authors plan to extend ElATools to analyze other tensorial properties, including piezoelectric and photoelastic tensors.Source: §4 · p. 25

4.1. Appendix A

  • The appendix defines elastic moduli, Poisson’s ratio, anisotropy indices, wave velocities, and related parameters for describing material properties.It notes that the universal anisotropy index is zero for isotropic materials and that deviation from zero measures elastic anisotropy.Source: §4.1 · pp. 26–28
  • The appendix describes Pugh’s ratio as a ductility-or-brittleness indicator, with a critical value of 1.75.Materials above 1.75 are described as ductile and those below it as brittle.Source: §4.1 · p. 26
  • The appendix gives Clarke’s model as the method currently used by ElaTools to calculate minimum thermal conductivity.Source: §4.1 · p. 26
  • For hardness prediction, ElATools uses semi-empirical correlations, with a selection guide based on crystal class and bandgap.The guide draws on the model proposed by Singh et al..Source: §4.1 · p. 28, Table II · p. 29

Type of material Cubic Hexagonal Orthorhombic Rhombohedral General

The appendices organize the elastic properties and anisotropy indices calculated for different material classes and document file and visualization outputs. They also report that infinite anisotropy of an elastic modulus corresponds to anomalous mechanical properties.

Type of material Cubic Hexagonal Orthorhombic Rhombohedral General

4.2. Appendix B

  • Appendix B lists input commands and files for the visualization executables and their generated output formats.The listed tools include dat2gnu.x, dat2agr.x, dat2wrl.x, and dat2html.x.Source: §4.2 · p. 30, Table IV · p. 30, Table V · p. 31

4.3. Appendix C

  • Appendix C lists main elastic properties and anisotropy indices for selected ZnAu2(CN)4, CrB2, GaAs, δ-phosphorene, and Pd2O6Se2 compounds.Source: §4.3 · p. 33, Table VIII · p. 34, Table IX · p. 34
  • When the anisotropy measure A_M is infinite, the material has anomalous mechanical properties.The authors note that marked anisotropy is often associated with negative Poisson’s ratio and negative linear compressibility.Source: §4.3 · p. 33

4.4. Appendix D

  • Appendix D lists default colors used to visualize elastic properties and points to supplementary information for color personalization.Source: §4.4 · p. 36, Table XIV · p. 36
Loading 2105.07279v3…