Department of Photonics Engineering, Technical University of Denmark, Kgs. Lyngby, 2800, Denmark
College of Science and Mathematics, The University of West Georgia, 1601 Maple St, Carrollton, GA 30117. USA
PSE Division, King Abdullah University of Science and Technology, Thuwal 23955-6900, Kingdom of Saudi Arabia
PSE Division, King Abdullah University of Science and Technology, Thuwal 23955-6900, Kingdom of Saudi Arabia
^{1}Department of Photonics Engineering, Technical University of Denmark, Kgs. Lyngby, 2800, Denmark
^{2}College of Science and Mathematics, The University of West Georgia, 1601 Maple St, Carrollton, GA 30117. USA
^{3}PSE Division, King Abdullah University of Science and Technology, Thuwal 23955-6900, Kingdom of Saudi Arabia
This article is part of the Thematic Series "Silicene, germanene and other group IV 2D materials".
Guest Editor: P. Vogt
Beilstein J. Nanotechnol. 2017, 8, 1345–1356. https://doi.org/10.3762/bjnano.8.136
Received 16 Feb 2017,
Accepted 14 Jun 2017,
Published 30 Jun 2017
A general three-dimensional continuum model of phonons in two-dimensional materials is developed. Our first-principles derivation includes full consideration of the lattice anisotropy and flexural modes perpendicular to the layers and can thus be applied to any two-dimensional material. In this paper, we use the model to not only compare the phonon spectra among the group-IV materials but also to study whether these phonons differ from those of a compound material such as molybdenum disulfide. The origin of quadratic modes is clarified. Mode coupling for both graphene and silicene is obtained, contrary to previous works. Our model allows us to predict the existence of confined optical phonon modes for the group-IV materials but not for molybdenum disulfide. A comparison of the long-wavelength modes to density-functional results is included.
Keywords: graphene; molybdenum disulfide; phonon; silicene; two-dimensional materials
Phonon spectra in two-dimensional (2D) nanomaterials have almost exclusively been computed using density-functional theory (DFT) based codes. One of the earliest applications to group-IV elemental 2D materials was for the important prediction of the stability of silicene and germanene [1]. These are complex calculations and prone to qualitative errors due to the various approximations such as convergence criteria and use of approximate functionals [2]. Even the stability or not of a given structure could be incorrectly inferred. For example, borophene and indiene have been predicted to be unstable in one paper [3], though other calculations (and in the case of borophene, even experiments [4]) have obtained opposite results [5,6]. In addition to obtaining a spectrum, it is often also useful to be able to predict and/or interpret properties of the spectrum based upon either microscopic or symmetric arguments. An excellent example is the prediction of a Dirac cone for silicene on the basis of symmetry [7] when DFT calculations either failed to recognize it [8] or were unable to explain it [1].
An alternative model of lattice vibrations is a classical continuum model, which is expected to reproduce most accurately phonons with wavelengths longer than lattice separations, i.e., near k = 0. One of the earliest such models applied to an ionic crystal slab is that of Fuchs and Kliewer [9], from which they deduced that transverse optical (TO) and longitudinal optical (LO) modes have different frequencies at k ≈ 0 as well as the existence of surface optical (SO) modes with an exponential dependence away from the slab. Slightly different models have been introduced by a number of authors for graphene [10-13]. One commonality is to treat the sheet as strictly two-dimensional. Additionally, instead of deriving the phonon dispersion relations from first principles, they all assumed the known results that there are in-plane acoustic modes with linear dispersions and out-of-plane transverse acoustic modes with quadratic dispersions (the latter being consistent with the elastic theory of thin plates) to construct either a Lagrangian [10] or equations of motion [11]. Goupalov also considered optical phonons but simply parameterized the dispersion relations to match the experimental data. In all of the above, the out-of-plane vibrations were assumed decoupled from the in-plane ones.
In this paper, a continuum theory of acoustic and optical phonons in 2D nanomaterials is derived from first principles, contrary to earlier approaches, by starting with the elastic and electric equations, and taking into account the full crystalline symmetry and piezoelectric couplings when allowed by symmetry. We apply the theory to obtain the phonons in group-IV elemental 2D materials. Given that there are two fundamental structures for the free-standing sheets (the flat hexagonal structure of graphene and the buckled hexagonal structure of the other elements – we will refer to silicene as the prototypical example of the latter), we will consider both of them. It should be recognized that, to date, silicene has been grown on substrates in different reconstruction state. The reconstruction is an atomic-scale distinction that is not describable by the current continuum model. Substrate effects on the phonons can be considered in an extended model that would need to be solved numerically. This can be studied in a future publication as it is beyond the analytical solutions sought for in the current manuscript. Furthermore, we have included a study of a well-known compound 2D material (molybdenum disulfide MoS_{2}) in order to further understand the properties derived for the elemental materials.
The 2D materials will be treated as 3D thin-plate materials in the following, well aware that the out-of-plane dimension contains one or a few atomic layers. In general, this would allow one to study multilayers though we will only focus on monolayers in this paper. Nonetheless, it is important to keep the third coordinate in the analysis to reveal the true phonon dispersions as observed experimentally and in DFT calculations.
The general 3D elastic equations are given by the equation of motion [14]
where T_{ik} is the stress tensor, ρ is the mass density, and u_{i} is the displacement. Equation 1 contains all the physics of the problem. In order to simplify and then solve it, it is necessary to specify the crystal symmetry of the vibrating system.
The three problems considered in this article are graphene, silicene and MoS_{2}. The Bravais lattice symmetries of the single-layer graphene (D_{6h} ≡ 6/mmm) and MoS_{2} () structures belong to the hexagonal system, while silicene belongs to the trigonal system (point group D_{3}_{d}). Graphene and silicene are non-piezoelectric materials because of the inversion symmetry of the unit cell, while MoS_{2} is piezoelectric because its unit cell exhibits inversion asymmetry.
For both graphene and MoS_{2}, the general form of the stiffness tensor for hexagonal structures is
and the stress–strain relations T_{I} = c_{IJ}S_{J} for graphene become
where T_{I} and S_{J} denote stress and strain, respectively. Here we have used Voigt notation for tensors. The latter equations are different for MoS_{2} due to the presence of piezoelecticity.
Combining Equation 1 and Equation 3, we obtain for graphene
where ω is the vibrational (angular) frequency. The displacements u_{x} and u_{y} are in-plane displacements, and u_{z} is the out-of-plane displacement.
Expressing the latter set of equations in the displacements, we get
A solution of the combined system (Equation 7–Equation 9) with appropriate boundary conditions allows for the determination of the displacements u_{x}, u_{y}, u_{z}. Note that the displacements in the three directions are coupled. This is different from the 2D model assumed before, which led to a decoupling of the out-of-plane vibrations. The coupling is a consequence of the finite thickness of the sheet with no mirror symmetry imposed. Thus, our model is sufficiently general to apply to multilayers. Earlier DFT calculations [15] had argued that there is no coupling between the in-plane and out-of-plane modes due to the mirror symmetry of graphene, leading to fewer scattering channels and, therefore, a higher thermal conductivity compared to, e.g., silicene. Our model reveals that such mode coupling, even when present, would occur for large k_{z} values (due to the small thicknesses) and, hence, would be unlikely to have a significant impact.
The phonons are normal mode solutions to the equations of motion. To proceed from the graphene elastic equations, we make the following plane wave ansatz
where k_{x} and k_{y} are wave numbers. Inserting Equation 10–Equation 12 into Equation 7–Equation 9 yields the following matrix expression in the unknown functions f_{x}, f_{y}, and f_{z}:
where , and
A semi-analytic solution can be easily obtained for the case when k_{y} = 0. In this case, C = D_{2} = 0 and the u_{y} displacement decouples from the displacements u_{x} and u_{z}. It follows that f_{y} obeys the wave equation
The solution to this differential equation can be found immediately by imposing the vacuum boundary condition
i.e.,
where 2h is the graphene layer thickness and −h,h define the graphene layer boundaries.
We obtain
where n = 0,1,2,3,…
For the coupled system f_{x}–f_{z}, the determinantal system obtained from Equation 13 leads to the solutions
where
and
Let us introduce the notation {α_{1},α_{2},α_{3},α_{4}} ≡ {α_{1},−α_{1},α_{2},−α_{2}} and
Then we have
where β_{i} are unknown coefficients.
A 4 × 4 matrix equation in β_{i} is finally obtained by invoking the mechanical boundary conditions (graphene in vacuum)
Observe that
Solving the determinantal equation for the latter 4 × 4 matrix equation as a function of k_{x} specifies a discrete set of (band) eigenfrequencies ω_{i}(k_{x}) and the corresponding eigenmodes f_{x}, f_{z} where i denotes the band index. A numerical solution reveals the out-of-plane mode to be quadratic in nature.
Surface optical phonons can be derived by solving for the electrostatic potential via the Maxwell–Poisson equation for the displacement field D_{i}.
For graphene,
Using the fact that the dielectric function is given by
we find
The Maxwell–Poisson equation in the absence of free charges reads
and finally becomes
The solution in the graphene slab becomes
where
This gives for the electric field components
Continuity in E_{x} at the slab interfaces with vacuum at z = ±h requires
and, similarly, continuity in the normal electric displacement requires
Solving the secular equation in leads to
which is equivalent to
We note in passing that the latter equation agrees with the result of Licari and Evrard [16] for interface optical phonon modes in cubic crystals where symmetry forces ε_{xx} = ε_{zz} ≡ ε. In the cubic case, however, confined optical phonon modes also exist at the LO phonon frequency since ε(ω_{LO}) = 0.
Confined optical phonons are also obtained by starting with the Maxwell–Poisson equation (Equation 31). The transverse electric field and normal electric displacement are
neglecting the DC component (static polarization) to the electric displacement. Confinement of the optical phonon modes implies . Assuming [type I]
where are constants, the Maxwell–Poisson equation and the boundary conditions are fulfilled if ε_{zz}(ω) = 0. Hence confined optical phonons exist in graphene.
We refer to silicene as the canonical example of the other group-IV materials even though the following derivation and qualitative properties apply to all of them since they all have the same symmetry of the buckled hexagonal structure, leading to the trigonal D_{3}_{d} point group.
The general form of the stiffness tensor for D_{3}_{d} trigonal structures is
and the stress–strain relations become
The only difference compared to the corresponding equations for graphene are the terms containing c_{14}.
The elastic equations for silicene then read
or, in terms of the displacements,
Again, the differences compared to graphene are the c_{14} terms. The electric equation for trigonal structures is exactly the same as for hexagonal structures since the permittivity matrix is the same. Hence, the general method for finding surface optical phonon modes repeats the description of graphene. Qualitatively, there is therefore no difference among the phonons for all the group-IV elemental 2D materials.
While the primary focus of this paper is on the properties of phonons of the group-IV elemental materials, it is important to know if there are properties of the phonons that are due to these being elemental. Hence, we will now consider the phonons for MoS_{2} as a prototypical compound 2D material because of its extensively studied properties. As already mentioned, the lack of inversion symmetry leads to the new phenomenon of piezoelectricity.
The stiffness tensor for MoS_{2} is the same as for graphene, Equation 2, because of the same hexagonal symmetry. However, the stress–strain relations are different because of the presence of piezoelectricity. Specifically, there are additional contributions to the stress–strain constitutive relations:
by use of the piezoelectric e-tensor:
and the equation
In the above, E_{j} is the internal electric field.
With the above result for the stress tensor including piezoelectric contributions, we find the elastic equations for MoS_{2}
This system of equations is only complete if we also solve for the Maxwell–Poisson equation. For MoS_{2}, the displacement field D_{i} is
with P_{sp} being the spontaneous polarization contribution. The dielectric function is the same as for graphene (Equation 28) and we find
where P_{sp,x}, P_{sp,y} and P_{sp,z} denote the spontaneous polarization components along the x-, x- and z-direction, respectively. The Maxwell–Poisson equation (Equation 30) becomes
because the spontaneous polarization P_{sp} is constant in space. A solution of the combined system Equation 51–Equation 53 and Equation 56 with appropriate boundary conditions allows for the determination of the electric field (or electric potential ) and the displacements u_{x}, u_{y} and u_{z}.
In the case of MoS_{2}, the general solutions are more complicated. One can still make the plane wave ansatz,
Now, inserting Equation 57–Equation 60 in Equation 51–Equation 53 and Equation 56 yields the 4 × 4 matrix expression in the unknown functions f_{x}, f_{y}, f_{z} and
where
For the case, k_{y} = 0, we have C = D_{2} = E_{2} = 0 and u_{y} decouples from u_{x}, u_{z} and . The wave equation for the u_{y} mode is the same as for graphene. To determine the general solution for the other displacement components, we once again solve the determinantal equation for the 3 × 3 sub-matrix in the components f_{x}, f_{z} and . The result is
where
and six roots (three pairs of opposite signs) exist. It also follows that
Hence the general solution is, using a notation similar to the case of graphene,
where β_{i} are unknown coefficients.
The Maxwell–Poisson equation in vacuum () has the general solution
where are unknown constants. An 8 × 8 matrix equation in β_{i}, is finally obtained by invoking the four mechanical boundary conditions at z = h,−h
and the four electric boundary conditions at z = h,−h,
The mechanical boundary conditions are the same as for graphene,
but the electric boundary conditions are (written out)
where ε_{0} is the vacuum permittivity. Solving the determinantal equation for the 8 × 8 matrix equation as a function of k_{x} specifies a discrete set of (band) eigenfrequencies ω_{i}(k_{x}) and the corresponding eigenmodes f_{x}, f_{z} and , where i denotes the band index.
An important result concerning confined optical phonons can be obtained without solving the determinantal equations. From the elastic equations (Equation 51 and Equation 53) and the electric equation (Equation 56) we have, assuming k_{y} = 0,
Inspection of the above set of equations reveals that confined phonon solutions can be sought in the form [type I]
or in the form [type II]
apart from a multiplying factor exp(ik_{x}x − iωt). The above choice reflects an expectation that four roots are real (two pairs of opposite-signed roots) and two imaginary (one pair of opposite-signed roots).
From Equation 75–Equation 77 follows
where i = 1,2,3. Note that a maximum of two of the latter relations can be independent since the values of δ_{i} are found by setting the system determinant to zero.
A general confined optical phonon mode can now be written as a Fourier series expansion [type I]
since this construction implies E_{x}(z = ±h) = 0 and using continuity in the transverse electric field at the interfaces. Equation 75–Equation 77 must still apply, however, in particular combining Equation 81–Equation 82 demands
For each term m, this latter condition, in general, would lead to a different ω_{m}. However, a normal mode is characterized by a unique frequency ω. Hence, for confined optical phonon modes, only one term in m is allowed in the general Fourier series expansions above. Further, imposing continuity in the normal electric displacement component at z = ±h by use of the third equation in Equation 55 gives
and unless accidentally ε_{zz}(ω_{m}) = 0 (treated separately in the next paragraph). For ε_{zz}(ω_{m}) ≠ 0, it follows from Equation 80 that A_{x,m} = 0 and Equation 81 yields A_{z,m} = 0. The conclusion is that confined optical phonons in the general case cannot exist in MoS_{2}! Note that for graphene, since it is a non-piezoelectric material, confined optical phonon modes do exist. The difference between the two materials lies in Equation 77, specifically in the term containing the piezoelectric coefficient.
If ε_{zz} = 0, continuity of the normal electric displacement component requires, as before, and we return to the situation treated in the previous subsection about confined optical phonons. We need to consider two possible cases: (a) ε_{xx}(ω) = 0 when ε_{zz}(ω) = 0 and (b) ε_{xx}(ω) ≠ 0 when ε_{zz}(ω) = 0. In case (a), we immediately obtain from continuity in the normal electric displacement. Further, Equation 77 yields u_{x} = 0 everywhere as a function of z. Equation 83–Equation 85 must hold from continuity in the transverse electric field. Finally, Equation 75 gives
and Equation 76 requires [type I]
and the condition
which is not fulfilled, unless accidental degeneracy applies, simultaneous with ε_{xx}(ω) = ε_{zz}(ω) = 0. In case (b), Equation 86 can be used and, unless accidental degeneracy applies, the obtained ω_{m} values do not fulfill ε_{zz}(ω) = 0. In conclusion, confined optical phonon modes do not exist in the case ε_{zz} = 0.
In order to illustrate the validity of the phonon dispersion relations obtained from our model, we compare them to DFT calculations. The continuum theory will require as input elasticity constants, piezoelectric coefficients, and dielectric functions.
We first give the standard phonon dispersion relation as obtained from DFT calculations (Figure 1). They are obtained from first principles calculations using the Vienna ab initio simulation package (VASP) [17] with a kinetic energy cut-off of 500 eV in the expansion of the electronic wave functions. Four C and six Mo and S valence electrons are considered. The generalized gradient approximation of the exchange–correlation potential in the Perdew–Burke–Ernzerhof flavor is employed [18]. Brillouin-zone integrations are performed using the tetrahedron method with Blöchl corrections [19]. We construct cells consisting of the monolayer and an approximately 15 Å thick vacuum region along the c direction. Structure optimizations are performed on Γ-centered 24 × 24 × 1 k-meshes. A direct method based on 4 × 4 × 1 supercells is used for obtaining the phonon dispersions within the harmonic approximation [20]. Forces are evaluated on 3 × 3 × 1 k-meshes, including long-range dipole contributions to the dynamical matrix following the method of [21]. Born effective charges and dielectric tensors are obtained within perturbation theory [22] and elastic constants by the homogeneous deformation method [23]. The elastic stiffness tensor and frequency-dependent dielectric tensor (independent particle approximation [24]) are calculated on Γ-centered 36 × 36 × 1 k-meshes.
The following parameters are found using DFT calculations as explained above: c_{11}d = 345 Pa·m, c_{12}d = 73 Pa·m, c_{13}d = 0.00387 Pa·m, c_{33}d = 0.531 Pa·m, c_{44}d = 0.0535 Pa·m, c_{66}d = 136 Pa·m, d = 3.4·10^{−8} m, and ρ_{2D} = 7.61·10^{−7} kg/m^{2}. For the permittivity data we used the DC values ε_{xx} = 4.4ε_{0} and ε_{zz} = 1.3ε_{0}.
In Figure 2, we show the frequency vs wavenumber (ω–k_{x}) dispersion in the vicinity of the Γ point (k_{y} = 0). Evidently, one mode shows a parabolic dispersion and one mode is linear. There are two higher-order modes originating from the boundary conditions along the z coordinate. The dispersion curves associated with the u_{y} vibrations decoupled from the u_{x} − u_{z} vibrations are not shown in the plots. We emphasize that for graphene optical and acoustic phonon modes decouple and are computed separately as described above. In Figure 3, the dispersion curves for optical phonon modes near the Γ point are shown.
For single-layer MoS_{2}, we use the following parameters: c_{11}d = 140 Pa·m, c_{12}d = 33 Pa·m, c_{13}d=-0.013 Pa·m, c_{33}d = 0.078 Pa·m, c_{44}d = −1.07 Pa·m, c_{66}d = 53.7 Pa·m, d = 6.2·10^{−10} m, ρ_{2D} = 3.1 ·10^{−6} kg/m^{2}, e_{x}_{1} = 0.5 C/m^{2}. These elasticity parameters were again computed using VASP. The frequency-dependent permittivity values are taken from [25] and the piezoelectric constant e_{x}_{1} is from Ref. [26]. In the case of MoS_{2}, all optical and acoustic phonon modes are computed by solving the combined set of elastic and electric equations and associated slab boundary conditions as described above. We note in passing that both confined and surface acousto-optical phonon modes are found using the present formalism. In Figure 4, we show the frequency vs wavenumber (ω–k_{x}) dispersion in the vicinity of the Γ point (k_{y} = 0).
A three-dimensional, first-principles, continuum elasticity theory was developed for two-dimensional materials. Piezoelectric materials required the simultaneous consideration of the electrostatics equations. Application to graphene, silicene and MoS_{2} revealed a number of interesting results. The out-of-plane vibrations were coupled to the in-plane motion for graphene (and for silicene as well), contrary to previous results using infinitely thin sheets. Acoustic modes with linear and quadratic dispersions were obtained, in agreement with experimental results and other models. We predict the existence of confined optical modes in all of the elemental group-IV materials and the nonexistence of confined optical modes in MoS_{2}. Our model will be applied to other 2D materials as well as to multilayers in future work. Additionally, the results of this model can be combined with the solution to the electron problem to compute electron–phonon scattering rates. The latter would be useful for understanding a number of physical properties as well as for applications.
1. | Cahangirov, S.; Topsakal, M.; Aktürk, E.; Şahin, H.; Ciraci, S. Phys. Rev. Lett. 2009, 102, 236804. doi:10.1103/PhysRevLett.102.236804 |
5. | Peng, B.; Zhang, H.; Shao, H.; Xu, Y.; Zhang, R.; Zhu, H. J. Mater. Chem. C 2016, 4, 3592–3598. doi:10.1039/C6TC00115G |
6. | Singh, D.; Gupta, S. K.; Lukačevic, I.; Sonvane, Y. RSC Adv. 2016, 6, 8006–8014. doi:10.1039/C5RA25773E |
16. | Licari, J. J.; Evrard, R. Phys. Rev. B 1977, 15, 2254–2264. doi:10.1103/PhysRevB.15.2254 |
4. | Mannix, A. J.; Zhou, X.-F.; Kiraly, B.; Wood, J. D.; Alducin, D.; Myers, B. D.; Liu, X.; Fisher, B. L.; Santiago, U.; Guest, J. R.; Yacaman, M. J.; Ponce, A.; Oganov, A. R.; Hersam, M. C.; Guisinger, N. P. Science 2015, 350, 1513–1516. doi:10.1126/science.aad1080 |
17. | Kresse, G.; Joubert, D. Phys. Rev. B 1999, 59, 1758–1775. doi:10.1103/PhysRevB.59.1758 |
3. | Kamal, C.; Chakrabarti, A.; Ezawa, M. New J. Phys. 2015, 17, 083014. doi:10.1088/1367-2630/17/8/083014 |
14. | Landau, L. D.; Lifshitz, E. M. Theory of Elasticity, 2nd ed.; Course of theoretical physics, Vol. 7; Pergamon Press: Oxford, United Kingdom, 1970. |
2. | Roome, N. J. Electronic and Phonon Properties of 2D Layered Materials. Ph.D. Thesis, University of Surrey, United Kingdom, 2015. |
15. | Xie, H.; Hu, M.; Bao, H. Appl. Phys. Lett. 2014, 104, 131906. doi:10.1063/1.4870586 |
9. | Fuchs, R.; Kliewer, K. L. Phys. Rev. 1965, 140, A2076–A2088. doi:10.1103/PhysRev.140.A2076 |
10. | Suzuura, H.; Ando, T. Phys. Rev. B 2002, 65, 235412. doi:10.1103/PhysRevB.65.235412 |
1. | Cahangirov, S.; Topsakal, M.; Aktürk, E.; Şahin, H.; Ciraci, S. Phys. Rev. Lett. 2009, 102, 236804. doi:10.1103/PhysRevLett.102.236804 |
11. | Goupalov, S. V. Phys. Rev. B 2005, 71, 085420. doi:10.1103/PhysRevB.71.085420 |
8. | Takeda, K.; Shiraishi, K. Phys. Rev. B 1994, 50, 14916–14922. doi:10.1103/PhysRevB.50.14916 |
7. | Guzmán-Verri, G. G.; Lew Yan Voon, L. C. Phys. Rev. B 2007, 76, 075131. doi:10.1103/PhysRevB.76.075131 |
10. | Suzuura, H.; Ando, T. Phys. Rev. B 2002, 65, 235412. doi:10.1103/PhysRevB.65.235412 |
11. | Goupalov, S. V. Phys. Rev. B 2005, 71, 085420. doi:10.1103/PhysRevB.71.085420 |
12. | Qian, J.; Allen, M. J.; Yang, Y.; Dutta, M.; Stroscio, M. A. Superlattices Microstruct. 2009, 46, 881–888. doi:10.1016/j.spmi.2009.09.001 |
13. | Droth, M.; Burkard, G. Phys. Rev. B 2011, 84, 155404. doi:10.1103/PhysRevB.84.155404 |
20. | Alfè, D. Comput. Phys. Commun. 2009, 180, 2622–2633. doi:10.1016/j.cpc.2009.03.010 |
18. | Perdew, J. P.; Burke, K.; Ernzerhof, M. Phys. Rev. Lett. 1996, 77, 3865–3868. doi:10.1103/PhysRevLett.77.3865 |
19. | Blöchl, P. E.; Jepsen, O.; Andersen, O. K. Phys. Rev. B 1994, 49, 16223–16233. doi:10.1103/PhysRevB.49.16223 |
25. | Mukherjee, B.; Tseng, F.; Gunlycke, D.; Amara, K. K.; Eda, G.; Simsek, E. Opt. Mater. Express 2015, 5, 447–455. doi:10.1364/OME.5.000447 |
26. | Zhu, H.; Wang, Y.; Xiao, J.; Liu, M.; Xiong, S.; Wong, Z. J.; Ye, Z.; Ye, Y.; Yin, X.; Zhang, X. Nat. Nanotechnol. 2015, 10, 151–155. doi:10.1038/nnano.2014.309 |
23. | Le Page, Y.; Saxe, P. Phys. Rev. B 2002, 65, 104104. doi:10.1103/PhysRevB.65.104104 |
24. | Gajdoš, M.; Hummer, K.; Kresse, G.; Furthmüller, J.; Bechstedt, F. Phys. Rev. B 2006, 73, 045112. doi:10.1103/PhysRevB.73.045112 |
21. | Cochran, W.; Cowley, R. A. J. Phys. Chem. Solids 1962, 23, 447–450. doi:10.1016/0022-3697(62)90084-7 |
22. | Baroni, S.; Giannozzi, P.; Testa, A. Phys. Rev. Lett. 1987, 58, 1861–1864. doi:10.1103/PhysRevLett.58.1861 |
© 2017 Willatzen et al.; licensee Beilstein-Institut.
This is an Open Access article under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
The license is subject to the Beilstein Journal of Nanotechnology terms and conditions: (http://www.beilstein-journals.org/bjnano)