Research Article | Open Access
Xiaowen Li, Xiaobin Qiang, Zhenhao Gong, Yubo Zhang, Penglai Gong, Lang Chen, "Tunable Negative Poisson’s Ratio in Van der Waals Superlattice", Research, vol. 2021, Article ID 1904839, 11 pages, 2021. https://doi.org/10.34133/2021/1904839
Tunable Negative Poisson’s Ratio in Van der Waals Superlattice
Negative Poisson’s ratio (NPR) materials are functional and mechanical metamaterials that shrink (expand) longitudinally after being compressed (stretched) laterally. By using first-principles calculations, we found that Poisson’s ratio can be tuned from near zero to negative by different stacking modes in van der Waals (vdW) graphene/hexagonal boron nitride (G/-BN) superlattice. We attribute the NPR effect to the interaction of orbitals between the interfacial layers. Furthermore, a parameter calculated by analyzing the electronic band structure, namely, distance-dependent hopping integral, is used to describe the intensity of this interaction. We believe that this mechanism is not only applicable to G/-BN superlattice but can also explain and predict the NPR effect in other vdW layered superlattices. Therefore, the NPR phenomenon, which was relatively rare in 3D and 2D materials, can be realized in the vdW superlattices by different stacking orders. The combinations of tunable NPRs with the excellent electrical/optical properties of 2D vdW superlattices will pave a novel avenue to a wide range of multifunctional applications.
Negative Poisson’s ratio (NPR) material shrinks laterally when axially compressed or laterally expands when subjected to axial stretching. Compared to positive Poisson’s ratio (PPR) materials, NPR material has higher indentation resistance , larger impact resistance, more superior sound absorption performance, and better crack propagation resistance . In addition, when subjected to out-of-plane bending moments, the NPR material will exhibit a dome shape rather than the PPR material tending to saddle shape. These excellent properties indicate that the NPR materials have broad application prospects in the automotive, aerospace, marine, and other industrial fields .
The NPR phenomena have been found in many natural and artificial materials, such as cubic metals , porous polymers , honeycombs , perovskite [7, 8], silicon oxides [7, 8], ceramic aerogels , reentrant crystal structures [10–17], and liquid crystal elastomer . Recently, the NPR effect has also been found in some two-dimensional (2D) materials , such as black phosphorus [20, 21], δ-phosphorene , borophene [23, 24], graphene , h-BN , 1 T-type transition metal dichalcogenides , group-IV monochalcogenides , Be5C2 , silicon dioxide [30, 31], FeB6 , B4N , and Ag2S . In addition, there are hundreds of thousands of materials in the inorganic crystal structure database (ICSD); however, the number of NPR materials reported in the study is a few hundred, which is relatively small. Among them, artificial materials and structures often have a very high NPR, while crystal materials have a relatively small NPR. Especially in 2D materials, NPR is smaller. In 2014, NPR was first found in 2D materials, namely, black phosphorus, with a value of -0.027 . After that, the NPR of the 2D materials found in the later study mostly remains near -0.1. Based on our theoretical calculations, unlike pure 2D materials, the NPR in the vdW heterostructure can only be maintained around -0.1, which is due to the expansion amplitude of orbital under in-plane strain.
Moreover, although many in-depth studies have examined the possible existence of NPR effects in 3D and 2D materials, few studies have reported the NPR effect in vdW superlattice. With the development of manufacturing technology, graphene-based superlattices showed enhanced stability in experiments . Therefore, the designability and diversity of vdW superlattices provide a broad prospect for designing multifunctional NPR materials, such as NPR electrodes and molecular sieves . In previous studies, the most NPR phenomena can be attributed to its unique wrinkle or re-entrant structures. In addition to these special geometric reasons, we want to study the fundamental physical mechanisms that form NPR effect.
In our previous study, we reported near-zero Poisson’s ratio (ZPR) phenomena in G/h-BN and multilayer h-BN . Interestingly, in this study, using first-principles calculations, we found that Poisson’s ratios of G/h-BN superlattices are -0.109, -0.111, and 0.023 in different stacking modes. The dichotomy between NPR and PPR effects exhibited in the G/h-BN superlattice, which can be explained by a special electronic structure at the interfacial layer. Although G/h-BN is a kind of simple vdW heterostructure, it is convenient to make theoretical analysis and calculation clearly. It may open a beginning for the theoretical study of NPR effect in different stacking modes of vdW materials. In addition, we calculated the out-of-plane stiffness of the G/h-BN superlattice with different stacking modes. These modes with NPR also have out-of-plane negative shear modulus (NSM), i.e., when shear strain is applied to NSM materials, as the shear strain increases, the corresponding shear stress tends to decrease.
Ultimately, the NPR phenomenon, which was relatively rare in 3D and 2D materials, can be realized in the vdW superlattices by different stacking orders as designed. Furthermore, studying on how to change the PPR material into the NPR material not only has important practical engineering application value but also theoretical value for in-depth study on other possible related interesting physical properties, such as negative pressure electricity, negative stiffness, and negative thermal expansion.
The unit cell of the G/h-BN superlattice is composed of graphene unit cell (2 C atoms) and -BN unit cell (1 B atom and 1 N atom) in the - plane. The lattice constant of the pristine monolayer graphene and -BN are 2.465 Å and 2.509 Å, respectively. Then, the of G/h-BN superlattice calculated by first-principles calculations is 2.485 Å, and the lattice mismatch between graphene and -BN is less than 1%.
The interlayer binding energies () and equilibrium distances () of all stacking modes of G/-BN superlattices obtained by density functional theory (DFT) approach are given in Table S1. Here, we investigated three highly symmetric stacking modes of G/-BN superlattices: N atom sublattice on hexagonal C atom ring (stacking mode A), B atom sublattice on C atom ring center (stacking mode B), and N atom sublattice on C atom ring center (stacking mode C) . follows the order of , while follows the order of .
By analyzing the strain energy, the elastic constants of G/h-BN superlattices were derived from the linear fitting of the energy-strain relationship (Table S2). For hexagonal crystal, the in-plane mechanical properties of G/-BN superlattice are isotropic (, , ) . Young’s modulus is determined by elastic constants (see method section). Notably, the in-plane of the 2D material is the product of the of the corresponding 3D material and the effective thickness , and we took the as the effective thickness for each layer of 2D material.
Table 1 shows that the out-of-plane of the stacking mode A, B, and C is 44.9, 45.6, and 49.0 GPa, respectively. The smaller the of the stacking mode, the greater the corresponding . In addition, we calculated the in-plane of the G/-BN superlattice, which is almost equal to the sum of of the monolayer graphene and -BN. Therefore, this result explains the reason for the stability enhancement of carbon-based superlattices observed in experiments . However, the difference in of the superlattice is mainly due to the different in .
For monolayer materials, the unit of Young’s and shear modulus is Nm-1. For G/-BN superlattices, the unit of Young’s and shear modulus is GPa.
2.2. Poisson’s Ratio
We compared G/-BN superlattices in different stacking modes under different uniaxial strains along the direction () (Figure 1). For stacking modes A and B, the is auxetic for , but the same phenomenon was not found in stacking mode C. Interestingly, stacking modes A and B not only have the NPR () effect but they also have negative shear modulus (NSM) () in the out-of-plane direction. The shear force decreases with the increase of shear deformation, which is the NSM effect.
In order to study the anisotropy of Poisson’s ratio of these materials, orientation-dependent Poisson’s ratio was calculated (Figure 2). We found that the stacking modes A, B, and C have the ZPR () effect at , 36.8°, and 18.7°, respectively. Therefore, stacking modes A and B exhibit a NPR effect in a larger crystal orientation angle range than stacking mode C.
To our knowledge, it is very difficult to measure Poisson’ s ratio of several layers of two-dimensional (2D) material with the existing experimental method. Because for these ultra-thin films, when the in-plane strain is applied, the out-of-plane deformation is very small and difficult to observe. However, for multilayer 2D materials, X-ray diffraction can be easily used to measure Poisson’ s ratio when the thickness is close to 10 nm . The NPR effect is generated in the interfacial layer. Therefore, both multi and single-layer vdW materials can exhibit the same NPR. It is relatively easy to measure Poisson’ s ratio for the multilayer vdW materials with a certain thickness.
3.1. Interlayer Binding Energy
Assuming that the interaction between the two layers of the superlattice is additive, the binding potential can be expressed as the cumulative interaction of atoms between different layers . The binding energy of two atoms combined by vdW forces can be expressed by the Lennard-Jones potential: Here, represents the distance between the two atoms. The and are fitting constants. The first term represents the vdW attraction, and the second term represents Pauli’s repulsion . Therefore, the interlayer potential of the vdW superlattice can be expressed as where and are the mass densities of two layers of vdW superlattice, respectively. The distance is obtained from the geometric relationship between the coordinate and . For G/-BN superlattice, . In Figure 2, the fitting curves of ~ below the horizontal coordinate axis represent vdW attraction, and the curves of ~ above the horizontal coordinate axis represent Pauli repulsion. Therefore, the low-order vdW term plays a major role in the large , while high-order Pauli’s repulsion term plays a major role in the small . According to the first-principles calculation, Eq. (2) can well describe the of the vdW superlattice.
Figures 2(d) and 2(e) show that when , the Pauli repulsion energy increases significantly, while the vdW attraction energy has a negligible change. Therefore, with the decrease of , of the G/-BN superlattice expands, while the lowest point of the energy curve moves forward along the positive direction of the axis, resulting in a negative of -0.109 and -0.111. For the stacking mode C, vdW and Pauli repulsion show negligible changes under (see Figure 2(f)). The binding energy increases slightly; so, the G/h-BN superlattice exhibits the ZPR effect. According to the first-principles calculation, Poisson’s ratio of the material has a relationship with . If the material can exhibit the NPR effect, the value of under tensile strain () is greater than the initial value of (Table S4). However, we calculated that the increase of is not a necessary condition for the NPR effect, i.e., plays a major role in the NPR effect.
3.2. Relationship between Orbitals and NPR
The Bloch wave function of the orbital electron in a periodic lattice under the tight binding (TB) approximation can be expressed as For G/-BN superlattice, C atoms in graphene and N atoms in h-BN have orbitals. When the atom A is used as the origin of coordinates (Figures 3(a) and 3(b)), let the in-plane strain be the perturbation . The wave function of orbitals of the atom A is . Meanwhile, the wave function of orbitals of the atom B is . When the A and B atoms in the lattice are bonded, the wave function of the bonded pz orbital can be expressed by the linear combination of atomic orbitals as The constant C is a normalization constant, which should satisfy the normalization condition . The density-weighted length of the electrons in the out-of-plane direction can be expressed as Here, is the length of the electrons with momentum in the out-of-plane direction. Therefore, the length () of electrons with all momentum should be the integral of in the first Brillouin zone (BZ). Finally, the charge density-weighted length of the orbital in the out-of-plane direction can be obtained: where and , which is the length of the isolated orbital (Figure 3(c)). Therefore, according to Eq. (6), we got the analytical solution of the relationship between and (Figure 3(d)). The calculation details can be found in the Supporting Information. Meanwhile, partially differentiate to the in-plane perturbation is . Therefore, orbitals will extend out-of-plane under in-plane tensile strain.
Since the length of the orbital has auxetic effect under in-plane strain, we quantitatively studied the charge density distribution along the out-of-plane direction of G/-BN superlattice by using first-principles calculations. In the out-of-plane direction, the charge density at coordinate can be expressed as Here, is the charge density at the coordinate (, , ) with the energy of , and is the Fermi level of the system. In order to quantify the change of the charge density in the out-of-plane direction under stress, we calculated the weighted length of the electron density in the out-of-plane direction according to the following formula: Notably, the charge of the graphene in the out-of-plane direction is mainly contributed by the orbital. The of each layer in the G/-BN superlattice under in-plane strain and is shown in Table S4. In each layer of the G/-BN superlattice, the is elongated. When an in-plane strain was applied, the bond angle ∠NBN increased from 120 to 122.36, resulting in the charge localization (Fig. S2).
Quantitatively, we found that the value of of the monolayer -BN and graphene in G/-BN superlattices increased by 1.8 ~ 1.9% and 2.3 ~ 2.4%, as the in-plane tensile strain increases by 8% (Table S5), explaining the NPR effect in stacking modes A and B along the out-of-plane direction. This is consistent with the analytical solution obtained by TB approximation (Figure 3(d)). For stacking mode C, the N atom sublattice is on the C atom ring center. In -BN, the N atom has a fully filled orbital, while the B atom has an empty orbital. The orbitals of G/-BN superlattices hardly overlap (Figure 4(c)); so, the change of the orbitals has little effect on the Pauli repulsion between the interfacial layers, resulting in no significant NPR effect.
3.3. Relationship between Electronic Band Structures and NPR
Figure 5 shows the DFT and TB-based band structure of G/-BN superlattices in different stacking modes. To further understand the first-principles calculation results, we adopted the TB model to describe the electrons in G/-BN superlattices with different stacking modes near the Fermi level. During the TB calculation, a unit cell contains two C1 and C2 carbon atoms at different positions and one N atom. Since the electronic states of the three bands around the Fermi level are completely contributed by the orbitals of C1, C2, and N atoms, only the orbitals of C1, C2, and N atoms are included in the TB model. The Hamiltonian matrix can be written as where subscripts 1, 2, and 3 represent C1, C2, and N atoms, respectively. Because the interlayer distance is longer than the C-C bond length, the nearest-neighbor interaction between C and N atoms and the next nearest-neighbor interaction between C and C atoms are considered (detailed Hamiltonian matrix elements can refer to the Supporting Information). The distance-dependent hopping integral is determined by the formula Here, represents the interfacial layer equilibrium distance, and is the hopping integral between orbitals at . is the distance between the th and th atoms, and is the decay constant for the integral . For the G/-BN superlattices, the values of and can refer to the Supporting Information Table S6.
The distance-dependent hopping integral () describes the intensity of the interaction between the orbitals of C and N atoms. Therefore, is a power-exponential function of the interlayer spacing and is proportional to the NPR (Table 2). Furthermore, the higher the value of , the greater the value of the corresponding NPR. Note that after Taylor expansion of the , the quadratic term is the previous research results [46, 47].
Consequently, the vdW superlattice can exhibit an NPR effect only if they have orbitals in the out-of-plane direction, and the orbitals overlap between the interfacial layers. Meanwhile, the NPR effect in all vdW materials can be explained by the same physical mechanism given in this section. For example, for lattice-matched materials, a previous study showed that AA-stacked -BN (a N atom on a N atom in another layer) can exhibit an NPR effect  (Table 3). In addition, for lattice-mismatched vdW materials, the orbitals between the interfacial layers overlap; so, these materials should exhibit an NPR effect. For example, it has been observed that WS2/WSe2 heterostructure expands abnormally under engineering tensile strain . Therefore, according to this physical mechanism, the NPR phenomenon should exist in a large number of vdW materials, which was considered as a rare phenomenon in bulk and monolayer 2D materials.
Moreover, in the experiment, the isolated atomic layers can also be reassembled into the designed heterostructure layer by layer in a precisely selected order . Therefore, for the same kind of investigated material, it can also be switched in different stacking modes through experimental methods. For example, the G/-BN superlattices may also be tuned among stacking modes A, B, and C. Similarly, we can change the material without strong interlayer orbital interaction into a material with orbitals strongly overlapping between the interfacial layers, thus exhibiting an NPR effect.
In conclusion, we studied Poisson’s ratios and the binding energies of G/-BN superlattices in different stacking modes by using the first-principles method. We found that the stacking mode C has a ZPR effect at the interfacial layer, while the stacking modes A and B show NPR effects. The NPR effect is mainly due to the interaction of the orbitals between the interfacial layers. Furthermore, the distance-dependent hopping integral () calculated by analyzing that the electronic band structure can be used to describe the intensity of this interaction. The is a power-exponential function of the interlayer spacing and is proportional to the NPR. Moreover, we calculated their Young’s and shear modulus and found that the stacking modes A and B also have NSM effect in the out-of-plane direction. These materials with negative index coexistence will provide broad prospects for multifunctional and multipurpose materials. Finally, we expect that the theory can be verified by experiments and provide a solid foundation for the large-scale searching and predicting NPR materials in the future.
Based on density of functional theory, all first-principles calculations were implemented by the planewave projector augmented wave (PAW) method in Vienna ab initio simulation package (VASP) code . The exchange correlation functional adopted the generalized gradient approximation (GGA) of the Perdew−Burke−Ernzerhof (PBE) functional . In order to test the robustness of our results, the vdW-corrected functionals proposed by Grimme DFT + D2 , DFT + D3 , many-body dispersion (MBD) , and vdW-corrected functional optB88-vdW  methods were used in first-principles calculations. In this paper, the calculation results of functional optB88-vdW are given because of its good agreement with the experimental results, and the results obtained by different vdW-corrected methods are only slightly different in numerical value (the detailed results are in the Supporting Information).
The G/-BN superlattice was calculated by using Monkhorst–Pack K-point mesh. The energy cut-off value is 500 eV, and the structures were completely relaxed until their atomic Hellmann–Feynman forces were less than 0.005 eV/Å. The convergence criterion of energy in the self-consistency process is 10-6 eV. We also calculated electronic band structures for G/-BN superlattices by using the HSE06 hybrid functional .
To quantitatively characterize the mechanical properties of the interface, the interlayer binding energy () between the monolayer graphene and -BN is as follows: where , , and are the energies of the G/-BN superlattice, graphene, and h-BN, respectively. represents the in-plane area of the superlattice.
The elastic constant is defined by expanding the internal energy into Taylor series in elastic strain at constant entropy. The expansion coefficient in the Taylor series is the elastic constant : where and are the initial mass density and the Lagrangian strains of the material . In this work, we use contracted notations (1) as tensor indices. In addition, we define Lagrangian strains and . The compliance coefficients are defined as Young’s modulus for the material is computed by Poisson’s ratio is defined as where is the strain in the direction of uniaxial loading (in the -direction), and is the resulting strain in the transverse direction (the -direction). In our calculations, we applied different uniaxial strains to the lattice. This strained structure was then completely relaxed to evaluate the magnitude of the strain in the out-of-plane direction. The detailed calculation process of the relationship between θ and is provided in the Supporting Information.
All data needed to evaluate the conclusions in the paper are presented in the paper and supplementary materials. And additional data are available from the corresponding authors upon reasonable request.
Conflicts of Interest
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
This work is supported by the Center for Computational Science and Engineering of Southern University of Science and Technology, NSFC (National Natural Science Foundation of China) grant number 51972160, the Science and Technology Research Items of Shenzhen grant numbers JCYJ20170412153325679, JCYJ20180504165650580, JCYJ20190809142603695, and JCYJ20170817110302672, and High-level Special Funding (Nos. G02206303 and G02206403).
1: interlayer binding energies and equilibrium distances results-Supplementary Table S1 2: elastic constants results-Supplementary Table S2 3: volume modulus and shear modulus results-Supplementary Table S3 4: interlayer binding energy-Supplementary Table S4 5: relationship between orbitals and NPR-Supplementary Figure S1 6: charge density-Supplementary Figure S2, Table S5 7: charge density. Table 8: the relationship between Poisson's ratio and stiffness. The relationship between θ and 9. The relationship between electronic band structures and NPR-Supplementary Table S6. (Supplementary Materials)
- A. Alderson and K. L. Alderson, “Auxetic materials,” Proceedings of the Institution of Mechanical Engineers, Part G: Journal of Aerospace Engineering, vol. 221, no. 4, pp. 565–575, 2007.
- W. Yang, Z.-M. Li, W. Shi, B. H. Xie, and M. B. Yang, “Review on auxetic materials,” Journal of Materials Science, vol. 39, no. 10, pp. 3269–3279, 2004.
- C. Huang and L. Chen, “Negative Poisson's ratio in modern functional materials,” Advanced Materials, vol. 28, no. 37, pp. 8079–8096, 2016.
- F. Milstein and K. Huang, “Existence of a negative Poisson ratio in fcc crystals,” Physical Review B, vol. 19, no. 4, pp. 2030–2033, 1979.
- B. D. Caddock and K. E. Evans, “Microporous materials with negative Poisson's ratios. I. Microstructure and mechanical properties,” Journal of Physics D Applied Physics, vol. 22, no. 12, pp. 1877–1882, 1989.
- R. Lakes, “Foam structures with a negative Poisson's ratio,” Science, vol. 235, no. 4792, pp. 1038–1040, 1987.
- C. W. Huang, W. Ren, V. C. Nguyen et al., “Abnormal Poisson's ratio and linear compressibility in perovskite materials,” Advanced Materials, vol. 24, no. 30, pp. 4170–4174, 2012.
- A. Yeganeh-Haeri, D. J. Weidner, and J. B. Parise, “Elasticity of agr-cristobalite: a silicon dioxide with a negative Poisson's ratio,” Science, vol. 257, no. 5070, pp. 650–652, 1992.
- X. Xu, Q. Zhang, M. Hao et al., “Double-negative-index ceramic aerogels for thermal superinsulation,” Science, vol. 363, no. 6428, pp. 723–727, 2019.
- S. Babaee, J. Shim, J. C. Weaver, E. R. Chen, N. Patel, and K. Bertoldi, “3D soft metamaterials with negative Poisson's ratio,” Advanced Materials, vol. 25, no. 36, pp. 5044–5049, 2013.
- Z. Y. Wei, Z. V. Guo, L. Dudte, H. Y. Liang, and L. Mahadevan, “Geometric mechanics of periodic pleated origami,” Physical Review Letters, vol. 110, no. 21, article 215501, 2013.
- J. L. Silverberg, A. A. Evans, L. McLeod et al., “Using origami design principles to fold reprogrammable mechanical metamaterials,” Science, vol. 345, no. 6197, pp. 647–650, 2014.
- A. U. Ortiz, A. Boutin, A. H. Fuchs, and F.-X. Coudert, “Anisotropic elastic properties of flexible metal-organic frameworks: how soft are soft porous crystals?” Physical Review Letters, vol. 109, no. 19, article 195502, 2012.
- L. J. Hall, V. R. Coluci, D. S. GalvaO et al., “Sign change of Poisson's ratio for carbon nanotube sheets,” Science, vol. 320, no. 5875, pp. 504–507, 2008.
- F. Song, J. Zhou, X. Xu, Y. Xu, and Y. Bai, “Effect of a negative Poisson ratio in the tension of ceramics,” Physical Review Letters, vol. 100, no. 24, article 245502, 2008.
- J. J. Williams, C. W. Smith, K. E. Evans, Z. A. D. Lethbridge, and R. I. Walton, “Off-Axis elastic properties and the effect of extraframework species on structural flexibility of the NAT-type zeolites: simulations of structure and elastic properties,” Chemistry of Materials, vol. 19, no. 10, pp. 2423–2434, 2007.
- M.-S. Pham, C. Liu, I. Todd, and J. Lertthanasarn, “Damage-tolerant architected materials inspired by crystal microstructure,” Nature, vol. 565, no. 7739, pp. 305–311, 2019.
- D. Mistry, S. D. Connell, S. L. Mickthwaite, P. B. Morgan, J. H. Clamp, and H. F. Gleeson, “Coincident molecular auxeticity and negative order parameter in a liquid crystal elastomer,” Nature Communications, vol. 9, no. 1, article 5095, 2018.
- R. Peng, Y. Ma, Q. Wu, B. Huang, and Y. Dai, “Two-dimensional materials with intrinsic auxeticity: progress and perspectives,” Nanoscale, vol. 11, no. 24, pp. 11413–11428, 2019.
- J. W. Jiang and H. S. Park, “Negative poisson's ratio in single-layer black phosphorus,” Nature Communications, vol. 5, article 4727, 2014.
- Y. Du, J. Maassen, W. Wu, Z. Luo, X. Xu, and P. D. Ye, “Auxetic black phosphorus: a 2D material with negative Poisson's ratio,” Nano Letters, vol. 16, no. 10, pp. 6701–6708, 2016.
- H. Wang, X. Li, P. Li, and J. Yang, “δ-Phosphorene: a two dimensional material with a highly negative Poisson's ratio,” Nanoscale, vol. 9, no. 2, pp. 850–855, 2017.
- L. Kou, Y. Ma, C. Tang, Z. Sun, A. Du, and C. Chen, “Auxetic and ferroelastic borophane: a novel 2D material with negative Possion's ratio and switchable dirac transport channels,” Nano Letters, vol. 16, no. 12, pp. 7910–7914, 2016.
- A. J. Mannix, X. F. Zhou, B. Kiraly et al., “Synthesis of borophenes: anisotropic, two-dimensional boron polymorphs,” Science, vol. 350, no. 6267, pp. 1513–1516, 2015.
- Q. Zhang, X. Xu, D. Lin et al., “Hyperbolically patterned 3D graphene metamaterial with negative Poisson's ratio and superelasticity,” Advanced Materials, vol. 28, no. 11, pp. 2229–2237, 2016.
- S. Woo, H. C. Park, and Y.-W. Son, “Poisson's ratio in layered two-dimensional crystals,” Physical Review B, vol. 93, no. 7, article 075420, 2016.
- L. Yu, Q. Yan, and A. Ruzsinszky, “Negative Poisson's ratio in 1T-type crystalline two-dimensional transition metal dichalcogenides,” Nature Communications, vol. 8, article 15224, 2017.
- X. Kong, J. Deng, L. Li et al., “Tunable auxetic properties in group-IV monochalcogenide monolayers,” Physical Review B, vol. 98, no. 18, article 184104, 2018.
- Y. Wang, F. Li, Y. Li, and Z. Chen, “Semi-metallic Be5C2 monolayer global minimum with quasi-planar pentacoordinate carbons and negative Poisson's ratio,” Nature Communications, vol. 7, article 11488, 2016.
- Z. Gao, X. Dong, N. Li, and J. Ren, “Novel two-dimensional silicon dioxide with in-plane negative Poisson's ratio,” Nano Letters, vol. 17, no. 2, pp. 772–777, 2017.
- V. O. Özçelik, S. Cahangirov, and S. Ciraci, “Stable single-layer honeycomblike structure of silica,” Physical Review Letters, vol. 112, no. 24, article 246803, 2014.
- J. Li, Y. Wei, X. Fan et al., “Global minimum of two-dimensional FeB6and an oxidization induced negative Poisson's ratio: a new stable allotrope,” Journal of Materials Chemistry C, vol. 4, no. 40, pp. 9613–9621, 2016.
- B. Wang, Q. Wu, Y. Zhang, L. Ma, and J. Wang, “Auxetic B4N monolayer: a promising 2D material with in-plane negative Poisson’s ratio and large anisotropic mechanics,” ACS Applied Materials & Interfaces, vol. 11, no. 36, pp. 33231–33237, 2019.
- R. Peng, Y. Ma, Z. He, B. Huang, L. Kou, and Y. Dai, “Single-layer Ag2S: a two-dimensional bidirectional auxetic semiconductor,” Nano Letters, vol. 19, no. 2, pp. 1227–1233, 2019.
- C.-M. Park and H.-J. Sohn, “Black phosphorus and its composite for lithium rechargeable batteries,” Advanced Materials, vol. 19, no. 18, pp. 2465–2468, 2007.
- K. Gopinadhan, S. Hu, A. Esfandiar et al., “Complete steric exclusion of ions and proton transport through confined monolayer water,” Science, vol. 363, no. 6423, pp. 145–148, 2019.
- X. Li, C. Huang, S. Hu et al., “Negative and near-zero Poisson's ratios in 2D graphene/MoS2and graphene/h-BN heterostructures,” Journal of Materials Chemistry C, vol. 8, no. 12, pp. 4021–4029, 2020.
- R. Quhe, J. Zheng, G. Luo et al., “Tunable and sizable band gap of single-layer graphene sandwiched between hexagonal boron nitride,” NPG Asia Materials, vol. 4, article e6, 2012.
- C. Lee, X. Wei, J. W. Kysar, and J. Hone, “Measurement of the elastic properties and intrinsic strength of monolayer graphene,” Science, vol. 321, no. 5887, pp. 385–388, 2008.
- K. Kim, W. R. L. Lambrecht, and B. Segall, “Elastic constants and related properties of tetrahedrally bonded BN, AlN, GaN, and InN,” Physical Review B, vol. 53, no. 24, pp. 16310–16326, 1996.
- K. Brugger, “Pure modes for elastic waves in crystals,” Journal of Applied Physics, vol. 36, no. 3, pp. 759–768, 1965.
- S. Hu, A. Alsubaie, Y. Wang et al., “Poisson's ratio of BiFeO3 thin films: X-ray reciprocal space mapping under variable uniaxial strain,” physica status solidi (a), vol. 214, no. 1, article 1600356, 2017.
- H. C. Hamaker, “The London--van der Waals attraction between spherical particles,” Physica, vol. 4, no. 10, pp. 1058–1072, 1937.
- O. V. Gritsenko, P. R. T. Schipper, and E. J. Baerends, “Effect of Pauli repulsion on the molecular exchange-correlation Kohn-sham potential: a comparative calculation of Ne2 and N2,” Physical Review A, vol. 57, no. 5, pp. 3450–3457, 1998.
- Z. Gong, X. Shi, J. Li et al., “Theoretical prediction of low-energy Stone-Wales graphene with an intrinsic type-III Dirac cone,” Physical Review B, vol. 101, no. 15, article 155427, 2020.
- S. Froyen and W. A. Harrison, “Elementary prediction of linear combination of atomic orbitals matrix elements,” Physical Review B, vol. 20, no. 6, pp. 2420–2422, 1979.
- L. Hu, J. Zhao, and J. Yang, “Nano-scale displacement sensing based on van der Waals interactions,” Nanoscale, vol. 7, no. 19, pp. 8962–8967, 2015.
- S. Xie, L. Tu, Y. Han et al., “Coherent, atomically thin transition-metal dichalcogenide superlattices with engineered strain,” Science, vol. 359, no. 6380, pp. 1131–1136, 2018.
- A. K. Geim and I. V. Grigorieva, “Van der Waals heterostructures,” Nature, vol. 499, p. 419, 2013.
- G. Kresse and J. Furthmüller, “Efficient iterative schemes forab initiototal-energy calculations using a plane-wave basis set,” Physical Review B, vol. 54, no. 16, pp. 11169–11186, 1996.
- J. P. Perdew, J. A. Chevary, S. H. Vosko et al., “Atoms, molecules, solids, and surfaces: applications of the generalized gradient approximation for exchange and correlation,” Physical Review B, vol. 46, no. 11, pp. 6671–6687, 1992.
- S. Grimme, “Semiempirical GGA-type density functional constructed with a long-range dispersion correction,” Journal of Computational Chemistry, vol. 27, no. 15, pp. 1787–1799, 2006.
- S. Grimme, J. Antony, S. Ehrlich, and H. Krieg, “A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu,” The Journal of Chemical Physics, vol. 132, no. 15, article 154104, 2010.
- T. Bučko, S. Lebègue, T. Gould, and J. G. Ángyán, “Many-body dispersion corrections for periodic systems: an efficient reciprocal space implementation,” Journal of Physics: Condensed Matter, vol. 28, no. 4, article 045201, 2016.
- T. Thonhauser, V. R. Cooper, S. Li, A. Puzder, P. Hyldgaard, and D. C. Langreth, “Van der Waals density functional: self-consistent potential and the nature of the van der Waals bond,” Physical Review B, vol. 76, no. 12, article 125112, 2007.
- J. Heyd, G. E. Scuseria, and M. Ernzerhof, “Hybrid functionals based on a screened coulomb potential,” The Journal of Chemical Physics, vol. 118, no. 18, pp. 8207–8215, 2003.
- R. N. Thurston and K. Brugger, “Third-order elastic constants and the velocity of small amplitude elastic waves in homogeneously stressed media,” Physical Review, vol. 133, no. 6A, pp. A1604–A1610, 1964.
- J. Zhao, J. M. Winey, and Y. M. Gupta, “First-principles calculations of second- and third-order elastic constants for single crystals of arbitrary symmetry,” Physical Review B, vol. 75, no. 9, article 094105, 2007.
Copyright © 2021 Xiaowen Li et al. Exclusive Licensee Science and Technology Review Publishing House. Distributed under a Creative Commons Attribution License (CC BY 4.0).