International Journal of Photoenergy

International Journal of Photoenergy / 2020 / Article
Special Issue

New Environmentally Friendly Solar Energy and Hydrogen Storage Materials

View this Special Issue

Research Article | Open Access

Volume 2020 |Article ID 8852665 |

Liwu Jiang, Meiling Wu, Peng Shi, Chuanhui Zhang, "Elastic, Mechanical, and Phonon Behavior of Orpiment Arsenic Trisulfide under Pressure", International Journal of Photoenergy, vol. 2020, Article ID 8852665, 6 pages, 2020.

Elastic, Mechanical, and Phonon Behavior of Orpiment Arsenic Trisulfide under Pressure

Academic Editor: Jiangwei LIU
Received13 Aug 2020
Revised20 Aug 2020
Accepted08 Sep 2020
Published19 Sep 2020


Arsenic trisulfide (As2S3) has been found to be an excellent glass former at high temperature and pressure. However, there is still some scarcity for the elastic and phonon behavior of the orpiment phase. By using the Dreiding force field of the geometry optimization computations, we investigated the elastic constants, mechanical moduli, and the phonon dispersion of orpiment As2S3 under the pressure from 0 to 5 GPa. Some results of the elastic parameters of orpiment-As2S3 at 0 GPa are consistent with the experimental data. The phonon dispersions for orpiment As2S3 under pressure are also reasonable with previous calculations.

1. Introduction

The computational simulation of materials has been a successful and rapid tool to understand the complex physical problems of materials. In addition, predicting the good elastic and mechanical properties of materials is a necessary demand for recent solid-state physics and science. And these mechanical properties at high pressure and temperature are also important for the development of modern technologies [14].

Arsenic trisulfide (As2S3) is the most studied compound due to its properties like large glass-forming ability, excellent infrared radiation transmission, and structural and mechanical properties [5, 6]. Under normal conditions, the stable phase is an orpiment layered structural semiconductor material. In its lattice, each As atom is surrounded by three S atoms, whereas each S atom is surrounded by two As atoms. At normal conditions, no other crystalline metastable phases of As2S3 could be artificially synthesized for a long time. The structure properties of arsenic sulfides were studied by the density functional theory [7], tight binding method [8], and nonempirical localized approach method [9]. The density of valence states of crystalline arsenic sulfides was investigated by using X-ray photoemission spectroscopy [10]. The optical and transport properties and electronic structure of Ni-doped arsenic sulfides were also studied [11]. Only one new phase of As2S3 named anorpiment has been naturally discovered recently [12]. It has low symmetry with triclinic structure and is made of covalently bonded layers of As and S atoms with a similar coordination arranged in a different stacking sequence by van der Waals force.

The high-pressure technique implies the active use of pressure as the fundamental thermodynamic variable for innovative material design. M. Güler and E. Güler reported that the classical force field can describe the surface, elastic, and phonon properties of wurtzite ZnS accurately, and the calculated data were in good agreement with the related experimental results [4]. Another previous work [13] pointed out that bonding disorder can be expected under high pressure in As2S3 glass, changing structural, vibrational, and optical properties. For example, using high-precision compressibility measurement and in situ high-energy X-ray diffraction to high pressure, the enhanced refractive index of the recovered As2S3 glass is promising for optical applications with improved functionalities [14]. A nearly permanent red shift in optical absorption after decompression can be used in high-impact-force optical sensors. In order to expand its application in the industrial field, it is necessary and interesting to study the elastic constants, mechanical moduli, and the phonon dispersion of orpiment As2S3 under high pressure.

In this paper, we provided the details of theoretical calculations with the employed interatomic potential in Section 2. Then, we discussed our results with the available experimental results in Section 3. Finally, we presented the main findings in Conclusions.

2. Computational Methods

The most significant feature of material modeling is the choice of the interatomic potentials. A proper interatomic potential can describe the physical properties of materials sufficiently and accurately. The empirical potentials are the modelling tools for materials because they can yield reasonable results. These kinds of potentials can provide reasonable results of the structures, energies, and elastic and mechanical properties of alloys [15, 16] and many other materials [17, 18].

In this work, all the theoretical calculations were carried out with the General Utility Lattice Program (GULP) molecular dynamics code [19]. This versatile code allows the concerned structures to be optimized at constant pressure (all internal and cell variables are included) or at constant volume (unit cell remains frozen). The Dreiding force field [20] was chosen to be the empirical potentials for simulations. It is widely used in research describing the properties of organic and inorganic materials [21, 22]. To avoid the constraints, constant pressure optimization was applied to the geometry of o-As2S3 cells with the Newton-Raphson method based on the Hessian matrix calculated from the second derivatives. In Figure 1, the cell geometry of o-As2S3 was assigned as , , , and with space group PNMA. In the geometry optimization calculations, the Hessian matrix was recursively updated with the BFGS [23] algorithm. The phonon and associated properties of o-As2S3 were also analyzed as a function of pressure within the ranges between 0 GPa and 5 GPa in steps of 1 GPa. We captured the phonon density of states and dispersions for a material after specifying a shrinking factor with phonon computations. Furthermore, phonons are described by calculating their values at points in a reciprocal space within the first Brillouin zone of the given crystal. To achieve the Brillouin zone integration and obtain the PDOS, we used a standard and reliable scheme developed by Monkhorst and Pack [24] with -point mesh.

3. Results and Discussion

As people know the volume of the crystal materials decreases, the density increases straightly under pressure. The density curve of o-As2S3 under pressure is shown in Figure 2. Under zero temperature, o-As2S3 has the lowest density value of 3.38 g/cm3 at 0 GPa and has the highest density value of 3.80 g/cm3 at 5 GPa. These results are consistent with the experimental value of 3.49 g/cm3 [25]. So our prediction crystal structure of o-As2S3 under pressure is accurate.

Then, we continued to compute the elastic and mechanical properties of o-As2S3 under pressure, such as the elastic constants (), bulk modulus (), shear modulus () and Young moduli (), and elastic wave velocities ( and ). The elastic constants were obtained by the second derivatives of the energy density with respect to strain. And the mechanical parameters were obtained by the linear combination of elastic constants using Voigt, Reuss, and Hill (VRH) approximations [26, 27].

Elastic constants deliver clear perceptions about the mechanical and other associated properties of materials. Though the elastic constants we obtained belong to the monocrystalline, these results also have important guiding significance for polycrystalline materials [28, 29]. For orthorhombic crystals, seven well-known elastic constants exist, which are specified as , , , , , , and [30]. Figure 3 shows the results of , , , , , , and constants for o-As2S3 under pressures between 0 GPa and 5 GPa. The calculated elastic constants are in the range of . The elastic constants and represent the longitudinal elastic character, and they increase under pressure for the reason of elastic wave propagation occurs easily under pressure. The values of and are similar to each other at 0 GPa, and the difference between them exists as the pressure increases. The elastic constant represents the shear elastic response to retarded wave propagation, and it also increases with the increasing pressure. In Table 1, our results for , , and of o-As2S3 at 0 GPa are consistent with the experimental results of McNeil et al. [31], but and are larger than the data in Ref. [31]. The reason is due to the anisotropy of o-As2S3 and the measuring technology of Brillouin scattering and the size of samples. McNeil’s measurement sample was made a natural crystal from Macedonia, and it was in the form of a thin plate roughly 5 mm square and 1 mm thick with only two good-quality surfaces (both perpendicular to the [010] axis). However, we used an orthorhombic crystal cell to do the simulation. So the anisotropy of o-As2S3 would influence the results of elastic parameters between us and the McNeil’s.

Parameter (GPa) (GPa) (GPa) (GPa) (GPa) (GPa) (GPa)

EXP. [30]99.620.627.121.91.1

Parameter (GPa) (GPa) (GPa) (km/s) (km/s)

According to stability, the mechanical stability criterion for orthorhombic crystals is ,, , and [30]. Our calculated results of elastic constants of o-As2S3 obey the mechanical stability criterion, which consequently indicates that the o-As2S3 is mechanically stable in its ground state.

The bulk modulus represents an essential elastic constant connected to the bonding strength and is used as a primary parameter for the calculation of a material’s hardness. The shear modulus represents the resistance to reversible deformations upon shear stress, while Young’s modulus is defined as the ratio of the tensile stress to the corresponding tensile strain. Figure 4 shows the , , and of o-As2S3 under pressure, and they all increase with the increasing pressure. However, the increase of and is obviously stronger than that of . This result points that o-As2S3 is an isotropic material, so the external pressure produces a uniform volume compression effect and much less effect on the shear deformation of the material. Usually, brittle materials display a considerable resistance to the deformation before fracture, whereas ductile materials would be deformed easily. Ductility is a mechanical property of the material that represents the ability to plastically deform before it undergoes the fracture. Brittleness refers to the nature of fracture failure of materials without obvious plastic deformation when the external force reaches a certain limit. Ductility and brittleness are important for the production of the desired materials, so we continued to study the ductile and brittle behaviors of o-As2S3 under pressure. The Pugh ratio evaluation () [32], C for ductile and brittle performances of materials, was considered in this section. A low value is associated with ductility, whereas a higher one corresponds to the brittle nature. The critical value which separates ductile and brittle materials is around 0.5. Our results show that values decrease from 0.42 () to 0.32 () at zero temperature for o-As2S3 as inserted in Figure 4. Therefore, o-As2S3 is a kind of a ductile material for the entire pressure range.

It is known that the longitudinal () and shear elastic () waves may arise in solids at low temperatures due to vibrational excitations originating from the acoustic modes [33]. And the pressure behavior of the and of o-As2S3 pressure at is shown in Figure 5. The indicates the longitudinal elastic wave velocity, and the indicates the shear wave velocity. We can see that has a significant increment compared to , and this is the most common case for materials. The reason is that the bulk modulus and Young’s modulus compared to the shear modulus have the same trend.

The ground state phonon dispersions of o-As2S3 were analyzed by using Dreiding potential quantitatively and compared in Reference [34] with theoretical results. We calculated the partial ground state phonon density of states (PDOS) of o-As2S3 to explain the contribution of As and S elements to the total phonon density of state (TDOS) of the material and showed them in Figure 6. The PDOS and TDOS appear with some separated regions corresponding to the longitudinal acoustic mode, transverse acoustic mode, longitudinal optic, and transverse optic mode of o-As2S3. Besides, the contribution of the As element to acoustic phonon modes is the same as that of S element, but in the high-frequency area, the case is different, and the contribution of the S element is dominant to optical modes. There is a big gap of 200 cm-1 (between the frequencies 300 cm-1 and 500 cm-1) originating from the mass differences of As and S elements of o-As2S3.

In addition, Figure 7 shows the phonon dispersion of o-As2S3 along the chosen -A path (the same as with original Ref [29]) in a reciprocal space for pressures 0 GPa and 5 GPa. It is seen that the pressure at 5 GPa shifts the phonon dispersion curves to higher frequency values. The reason is the atoms of o-As2S3 move towards each other and have to locate in the steeper potential wells under high pressures. In Figure 8, the corresponding PDOS curves of five pressures (0 GPa, 1 GPa, 2 GPa, 3 GPa, 4 GPa, and 5 GPa) are presented. The increasing pressure increases the PDOS peaks of o-As2S3 under the high pressures, and the gap between acoustic and optical modes also shifts to higher frequencies from 350 cm-1 to 500 cm-1.

In summary, the calculated parameters of o-As2S3 in this work are consistent with experiments and those of some published theoretical data. And some of the elastic constants and phonon properties of o-As2S3 at 0 GPa demonstrate fair accordance with the experimental and theoretical data. We hope that our results add value to the forthcoming researches about o-As2S3 under pressure.

4. Conclusions

In this paper, we applied the Dreiding interatomic potential for the geometry optimization on the structural properties of o-As2S3 and studied both the ground state and pressure-dependent elastic, mechanical, and phonon dispersions of o-As2S3. The present results for the ground state of o-As2S3 were former experiments for the elastic constants, and the phonon characteristics were consistent with the published theoretical data. Moreover, we presented the effect of pressure on o-As2S3 and obtained some reasonable results. The mechanical properties of bulk modulus, shear modulus, Young’s modulus, and longitudinal elastic character were increments under pressure. The o-As2S3 exhibits ductile character under high pressure. At last, the increasing pressure increases the PDOS peaks of o-As2S3 under the high pressures, and the gap between acoustic and optical modes also shifts to higher frequencies from 350 cm-1 to 500 cm-1.

Data Availability

The research data used to support the findings of this study are included in the article.

Conflicts of Interest

We declare that we have no conflict of interest about the manuscript.


This research was supported by the National Science and Technology Major Project (Grant No. 2017-VI-0012-0084), the Fundamental Research Funds for the Central Universities (Grant No. FRF-GF-19-029B), and the National Natural Science Foundation of China (Grant No. 51471022).


  1. O. Boudrifa, A. Bouhemadou, N. Guechi, S. Bin-Omran, Y. Al-Douri, and R. Khenata, “First-principles prediction of the structural, elastic, thermodynamic, electronic and optical properties of Li4Sr3Ge2N6 quaternary nitride,” Journal of alloys and compounds, vol. 618, pp. 84–94, 2015. View at: Publisher Site | Google Scholar
  2. D. Duan, X. Huang, F. Tian et al., “Pressure-induced decomposition of solid hydrogen sulfide,” Physical Review B, vol. 91, no. 18, article 180502, 2015. View at: Publisher Site | Google Scholar
  3. Z. Zeng, C. S. Garoufalis, S. Baskoutas, and G. Bester, “Excitonic optical properties of wurtzite ZnS quantum dots under pressure,” The Journal of Chemical Physics, vol. 142, no. 11, article 114305, 2015. View at: Publisher Site | Google Scholar
  4. M. Güler and E. Güler, “Theoretical analysis of elastic, Mechanical and Phonon Properties of Wurtzite Zinc Sulfide under Pressure,” Crystals, vol. 7, no. 6, 2017. View at: Publisher Site | Google Scholar
  5. V. V. Brazhkin, Y. Katayama, M. V. Kondrin, A. G. Lyapin, and H. Saitoh, “Structural transformation yielding an unusual metallic state in liquid As2S3 under high pressure,” Physical Review B, vol. 82, no. 14, article 140202, 2010. View at: Publisher Site | Google Scholar
  6. G. V. Gibbs, A. F. Wallace, R. Zallen et al., “Bond paths and van der Waals interactions in orpiment, As2S3,” The Journal of Physical Chemistry. A, vol. 114, no. 23, pp. 6550–6557, 2010. View at: Publisher Site | Google Scholar
  7. P. Srivastava, H. Singh Mund, and Y. Sharma, “Investigation of electronic properties of crystalline arsenic chalcogenides: theory and experiment,” Physica B: Condensed Matter, vol. 406, no. 15-16, pp. 3083–3088, 2011. View at: Publisher Site | Google Scholar
  8. Y. Watanabe, H. Kawazoe, and M. Yamane, “Imperfections in amorphous chalcogenides. III. Interacting-lone-pair model for localized gap states based on a tight-binding energy-band calculation for As2S3,” Physical Review B, vol. 38, no. 8, pp. 5677–5683, 1988. View at: Publisher Site | Google Scholar
  9. D. W. Bullett, “Electronic structure of arsenic chalcogenides,” Physical Review B, vol. 14, no. 4, pp. 1683–1692, 1976. View at: Publisher Site | Google Scholar
  10. S. G. Bishop and N. J. Shevchik, “Densities of valence states of amorphous and crystalline As2S3, As2Se3, and As2Te3: X-ray photoemission and theory,” Physical Review B, vol. 12, no. 4, pp. 1567–1578, 1975. View at: Publisher Site | Google Scholar
  11. Y. Sharma and P. Srivastava, “Optical and transport properties and electronic structure of nickel doped arsenic chalcogenides,” Computational Materials Science, vol. 53, no. 1, pp. 451–459, 2012. View at: Publisher Site | Google Scholar
  12. A. R. Kampf, R. T. Downs, R. M. Housley, R. A. Jenkins, and J. Hyrsl, “Anorpiment, As2S3, the triclinic dimorph of orpiment,” Mineralogical Magazine, vol. 75, no. 6, pp. 2857–2867, 2011. View at: Publisher Site | Google Scholar
  13. C. Y. Yang, M. A. Paesler, and D. E. Sayers, “Chemical order in the glassy AsxS1−x system: an X-ray-absorption spectroscopy study,” Physical Review B, vol. 39, no. 14, pp. 10342–10351, 1989. View at: Publisher Site | Google Scholar
  14. E. Soignard, O. B. Tsiok, A. S. Tverjanovich et al., “Pressure-Driven Chemical Disorder in Glassy As2S3 up to 14.7 GPa, Postdensification Effects, and Applications in Materials Design,” The Journal of Physical Chemistry B, vol. 124, pp. 430–442, 2020. View at: Publisher Site | Google Scholar
  15. C. H. Zhang, S. Huang, J. Shen, and X. N. Chen, “Chen’s lattice inversion embedded-atom method for Ni Al alloy,” Chinese Physics B, vol. 21, no. 11, article 113401, 2012. View at: Publisher Site | Google Scholar
  16. C. H. Zhang, S. Huang, J. Shen, and X. N. Chen, “Structural and mechanical properties of Fe-Al compounds: an atomistic study by EAM simulation,” Intermetallics, vol. 52, pp. 86–91, 2014. View at: Publisher Site | Google Scholar
  17. L. C. He, S. S. Guo, J. C. Lei, Z. D. Sha, and Z. S. Liu, “The effect of stone-thrower-Wales defects on mechanical properties of graphene sheets – a molecular dynamics study,” Carbon, vol. 75, pp. 124–132, 2014. View at: Publisher Site | Google Scholar
  18. C. H. Zhang, Y. Wang, and D. B. Sun, “The molecular dynamics simulation on the mechanical properties of Ni glass with external pressure,” International Journal of Modern Physics B, vol. 31, no. 20, article 1750138, 2017. View at: Publisher Site | Google Scholar
  19. J. D. Gale and A. L. Rohl, “The General Utility Lattice Program (GULP),” Molecular Simulation, vol. 29, no. 5, pp. 291–341, 2003. View at: Publisher Site | Google Scholar
  20. S. L. Mayo, B. D. Olafson, and W. A. Goddard, “DREIDING: a generic force field for molecular simulations,” The Journal of Physical Chemistry, vol. 94, no. 26, pp. 8897–8909, 1990. View at: Publisher Site | Google Scholar
  21. R. Devanathan, A. Venkatnathan, and M. Dupuis, “Atomistic simulation of nafion membrane: I. Effect of hydration on membrane nanostructure,” The Journal of Physical Chemistry. B, vol. 111, no. 28, pp. 8069–8079, 2007. View at: Publisher Site | Google Scholar
  22. P. G. Boyd, S. M. Moosavi, M. Witman, and B. Smit, “Force-field prediction of materials properties in metal-organic frameworks,” Journal of Physical Chemistry Letters, vol. 8, no. 2, pp. 357–363, 2017. View at: Publisher Site | Google Scholar
  23. R. Fletcher, “A new approach to variable metric algorithms,” The Computer Journal, vol. 13, no. 3, pp. 317–322, 1970. View at: Publisher Site | Google Scholar
  24. H. J. Monkhorst and J. D. Pack, “Special points for Brillouin-zone integrations,” Physical Review B, vol. 13, no. 12, pp. 5188–5192, 1976. View at: Publisher Site | Google Scholar
  25. D. J. E. Mullen and W. Nowacki, “Refinement of the crystal structures of realgar, AsS and orpiment, As2S3,” Zeitschrift für Kristallographie - Crystalline Materials, vol. 136, no. 1-2, pp. 48–65, 1972. View at: Publisher Site | Google Scholar
  26. W. Voigt, Lehrbuch der Kristallphysik (mit Ausschluss der Kristalloptik), Vieweg+Teubner Verlag, 1928.
  27. A. Reuss, “Berechnung der Fließgrenze von Mischkristallen auf Grund der Plastizitätsbedingung für Einkristalle,” ZAMM‐Journal of Applied Mathematics and Mechanics/Zeitschrift für Angewandte Mathematik und Mechanik, vol. 9, no. 1, pp. 49–58, 1929. View at: Publisher Site | Google Scholar
  28. S. Hirsekorn, “Elastic properties of polycrystals: a review,” Textures and Microstructures, vol. 12, no. 1-3, pp. 1–14, 1990. View at: Publisher Site | Google Scholar
  29. R. L. Fleischer, “Substitutional solutes in AlRu-I. Effects of solute on moduli, lattice parameters and vacancy production,” Acta Metallurgica et Materialia, vol. 41, no. 3, pp. 863–869, 1993. View at: Publisher Site | Google Scholar
  30. L. E. McNeil and M. H. Grimsditch, “Elastic constants of As2S3,” Physical Review B, vol. 44, no. 9, pp. 4174–4177, 1991. View at: Publisher Site | Google Scholar
  31. S. F. Pugh, “XCII. Relations between the elastic moduli and the plastic properties of polycrystalline pure metals,” Philosophical Magazine, vol. 45, no. 367, pp. 823–843, 1954. View at: Publisher Site | Google Scholar
  32. B. G. Dick and A. W. Overhauser, “Theory of the dielectric constants of alkali halide crystals,” Physics Review, vol. 112, no. 1, pp. 90–103, 1958. View at: Publisher Site | Google Scholar
  33. S. Radescu, A. Mujica, P. R. Hernandez et al., “Study of the orpiment and anorpiment phases of As2S3under pressure,” Journal of Physics: Conference Series, vol. 950, article 042018, no. 4, 2017. View at: Publisher Site | Google Scholar
  34. R. Hill, “The elastic behaviour of a crystalline aggregate,” Proceedings of the Physical Society. Section A, vol. 65, no. 5, pp. 349–354, 1952. View at: Publisher Site | Google Scholar

Copyright © 2020 Liwu Jiang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly as possible. We will be providing unlimited waivers of publication charges for accepted research articles as well as case reports and case series related to COVID-19. Review articles are excluded from this waiver policy. Sign up here as a reviewer to help fast-track new submissions.