No Arabic abstract
The MechElastic Python package evaluates the mechanical and elastic properties of bulk and 2D materials using the elastic coefficient matrix ($C_{ij}$) obtained from any ab-initio density-functional theory (DFT) code. The current version of this package reads the output of VASP, ABINIT, and Quantum Espresso codes (but it can be easily generalized to any other DFT code) and performs the appropriate post-processing of elastic constants as per the requirement of the user. This program can also detect the input structures crystal symmetry and test the mechanical stability of all crystal classes using the Born-Huang criteria. Various useful material-specific properties such as elastic moduli, longitudinal and transverse elastic wave velocities, Debye temperature, elastic anisotropy, 2D layer modulus, hardness, Pughs ratio, Cauchys pressure, Kleinman parameter, and Lames coefficients, can be estimated using this program. Another existing feature of this program is to employ the ELATE package [J. Phys.: Condens. Matter 28, 275201 (2016)] and plot the spatial variation of several elastic properties such as Poissons ratio, linear compressibility, shear modulus, and Youngs modulus in three dimensions. Further, the MechElastic package can plot the equation of state (EOS) curves for energy and pressure for a variety of EOS models such as Murnaghan, Birch, Birch-Murnaghan, and Vinet, by reading the inputted energy/pressure versus volume data obtained via numerical calculations or experiments. This package is particularly useful for the high-throughput analysis of elastic and mechanical properties of materials.
The dynamical, dielectric and elastic properties of GeTe, a ferroelectric material in its low temperature rhombohedral phase, have been investigated using first-principles density functional theory. We report the electronic energy bands, phonon dispersion curves, electronic and low frequency dielectric tensors, infra-red reflectivity, Born effective charges, elastic and piezoelectric tensors and compare them with the existing theoretical and experimental results, as well as with similar quantities available for other ferroelectric materials, when appropriate.
Anisotropic materials, with orientation-dependent properties, have attracted more and more attention due to their compelling tunable and flexible performance in electronic and optomechanical devices. So far, two-dimensional (2D) black phosphorus shows the largest known anisotropic behavior, which is highly desired for synaptic and neuromorphic devices, multifunctional directional memories, and even polarization-sensitive photodetector, whereas it is unstable at ambient conditions. Recently, 2D few-layered As2S3 with superior chemical stability was successfully exfoliated in experiments. However, the electronic and mechanical properties of monolayer and bilayer As2S3 is still lacking. Here, we report the large anisotropic electronic and mechanical properties of As2S3 systems through first-principles calculations and general angle-dependent Hookes law. Monolayer and bilayer As2S3 exhibit anisotropic factors of Youngs modulus of 3.15 and 3.32, respectively, which are larger than the black phosphorous with experimentally confirmed and an anisotropic factor of 2. This study provides an effective route to flexible orientation-dependent nanoelectronics, nanomechanics, and offers implications in promoting related experimental investigations.
Recently, a new class of carbon allotrope called protomene was proposed. This new structure is composed of sp2 and sp3 carbon-bonds. Topologically, protomene can be considered as an sp3 carbon structure (~80% of this bond type) doped by sp2 carbons. First-principles simulations have shown that protomene presents an electronic bandgap of ~3.4 eV. However, up to now, its mechanical properties have not been investigated. In this work, we have investigated protomene mechanical behavior under tensile strain through fully atomistic reactive molecular dynamics simulations using the ReaxFF force field, as available in the LAMMPS code. At room temperature, our results show that the protomene is very stable and the obtained ultimate strength and ultimate stress indicates an anisotropic behavior. The highest ultimate strength was obtained for the x-direction, with a value of ~110 GPa. As for the ultimate strain, the highest one was for the z-direction (~25% of strain) before protomene mechanical fracture.
Carbon nitride-based nanostructures have attracted special attention (from theory and experiments) due to their remarkable electromechanical properties. In this work we have investigated the mechanical properties of some graphene-like carbon nitride membranes through fully atomistic reactive molecular dynamics simulations. We have analyzed three different structures of these CN families, the so-called graphene-based g-CN, triazine-based g-C3N4 and heptazine-based g-C3N4. The stretching dynamics of these membranes was studied for deformations along their two main axes and at three different temperatures: 10K, 300K and 600K. We show that g-CN membranes have the lowest ultimate fracture strain value, followed by heptazine-based and triazine-based ones, respectively. This behavior can be explained in terms of their differences in terms of density values, topologies and types of chemical bonds. The dependency of the fracture patterns on the stretching directions is also discussed.
Molecular dynamics simulations play an increasingly important role in the rational design of (nano)-materials and in the study of biomacromolecules. However, generating input files and realistic starting coordinates for these simulations is a major bottleneck, especially for high throughput protocols and for complex multi-component systems. To eliminate this bottleneck, we present the polyply software suite that leverages 1) a multi-scale graph matching algorithm designed to generate parameters quickly and for arbitrarily complex polymeric topologies, and 2) a generic multi-scale random walk protocol capable of setting up complex systems efficiently and independent of the target force-field or model resolution. We benchmark quality and performance of the approach by creating melt simulations of six different polymers using two force-fields with different resolution. We further demonstrate the power of our approach by setting up a multi lamellar microphase-separated block copolymer system for next generation batteries, and by generating a liquid-liquid phase separated polyethylene oxide-dextran system inside a lipid vesicle, featuring both branching and molecular weight distribution of the dextran component.