Optical activity, the ability of certain media to rotate the plane of light polarization in the absence of a magnetic field, is a remarkable manifestation of polarization-dependent interaction between light and matter . It stems from the different strengths of coupling to right-handed and left-handed circularly polarized light leading to circular birefringence and circular dichroism . From the pioneering experiment by J.C. Bose in 1898 with left- and right-twisted jute elements , optical rotation is commonly associated with chirality. Optically active media attract much attention because they allow one to manipulate the polarization state of light. This is particularly important for superresolution imaging and the development of broadband optical components, such as filters and sensors with high polarization suppression ratio , . Moreover, photon polarization is currently considered as the degree of freedom to code and carry quantum information .
Optical activity, also referred to as gyrotropy, can be of either intrinsic or extrinsic nature. The former occurs in macroscopically homogeneous systems with chiral molecule components or crystals of the gyrotropic symmetry classes . A canonical example of optically active biological media is the tartaric acid, which was discovered by L. Pasteur and explained by the predominance of one of the two possible stereoisomers of the acid molecules . Tellurium is an example of gyrotropic crystals . This elemental semiconductor can exist in two enantiomorphic forms with the atoms, constituting the crystal lattice, bound into left-handed or right-handed chains. An example of gyrotropic but not chiral bulk crystals is silver gallium sulphide . Microscopically, optical activity originates from the spatial dispersion of susceptibility . It is increased in the vicinity of exciton resonances, as was observed in bulk wurtzite crystals  and II–VI quantum well structures . The term “extrinsic optical activity” coined recently is primarily used for metamaterials, metasurfaces, and photonic crystals composed of the arrays of chiral or achiral elements such as gammadions , , , G-shape nanostructures , and split-ring resonators  or chiral stacks of metasurfaces , . The efficiency of polarization conversion can be enhanced by plasmonic resonances in metal nanostructures , .
Apart from the above two distinct types of gyrotropic structures, we highlight the opportunity of an intermediate case of artificial materials with chiral stacking of atomic layers. Previously, they would be hypothetical. The emerging technology of van der Waals structures made of two-dimensional (2D) crystals of atomic thickness ,  enables the formation of such materials. Recently, it has been experimentally demonstrated that a pair of graphene layers stacked with a twist exhibits optical activity . The optical rotation angle normalized to the sample thickness was found to be several orders of magnitude larger than that in natural materials. The 2D crystals beyond the graphene, such as transition metal dichalcogenides (TMDC) MoS2 and WS2, are of particular interest for application in optics. TMDC layers have optical gaps and demonstrate strong light-matter coupling in the spectral range of exciton transitions , , , . Pronounced exciton resonances have been observed in the reflectance spectra of TMDC layers and thin films , , . The properties of twisted stacks including interlayer coupling strength are also being widely studied , , , , , . At the same time, the optical polarization effects related to chirality slipped away from the focus of previous research.
In this paper, we study the optical properties of chiral stacks of 2D semiconductors. We show that properly arranged stacks (see Figure 1) exhibit pronounced optical activity while individual layers do not. We develop a microscopic theory of the optical activity and the accompanying phenomenon of circular dichroism. It is also predicted that the photoluminescence (PL) of chiral stacks is circularly polarized.
The paper is organized as follows. In Section 2, we develop a theoretical approach to describe the exciton states in the chiral stacks of TMDC layers and the polarization optical properties of such stacks. We start by considering chiral bilayers and then generalize the model to multilayer stacks. In Section 3, we describe the spectral behavior of light transmission, optical activity, and circular dichroism for the chiral bilayers, multilayers, and thick TMDC stacks. In Section 4, we discuss the circularly polarized PL of chiral stacks. Section 5 summarizes the results of the paper.
The spatial symmetry of a single TMDC layer is described by the D3h point group. The layers are optically inactive and demonstrate isotropic linear response at the normal incident of radiation. We consider the optical response associated with bright direct excitons with the spin projections s=±1 along the normal to the layer that are excited by the σ± circularly polarized light, respectively. The excitons are formed by electrons and holes located at the K and K′ valleys of the 2D Brillouin zone. The exciton states with the spin projections s=±1 at zero in-plane wave vector are degenerate in energy.
Arranging two identical TMDC layers into a twisted stack reduces the point-group symmetry of the system to D3. The exciton states in such a bilayer are described by an effective Hamiltonian Hns,n′s′, where n, n′=1, 2 stand for the monolayer index and s, s′=±1 stand for the spin index. We note that in multilayer systems, there are also excitations consisting of electrons and holes localized in different layers , . Those interlayer excitons are far separated in energy from the intralayer excitons because of the reduced Coulomb interaction and are not considered here.
The D3 point group of the bilayer imposes restrictions on the form of the effective exciton Hamiltonian. The presence of the three-fold rotation axis along the bilayer normal eliminates the matrix elements of the Hamiltonian with s≠s′, so that
with real parameters ωx, J, and ϕ. Here, ωx is the exciton frequency in an isolated layer, J describes the interlayer hopping of excitons due to, e.g., tunneling or Förster excitation transfer, and ϕ is the phase acquired at the interlayer hopping due to the chirality of the bilayer. The Planck constant is set to unity. The dependence of the phase ϕ on the parameters of particular TMDC bilayers can be obtained from atomistic calculations, which are beyond the scope of the present paper. Instead, we focus below on the general optical properties of chiral bilayers and multilayers. We only note that the symmetry analysis yields that the dependence of the phase ϕ and the interlayer hopping parameter J on the twist angle is 2π/3-periodic and the optical activity vanishes in achiral stacks with the twist angles 0, π/3, 2π/3, etc. It was found in experiments and confirmed by atomistic calculations that in energetically favourable TMDC bilayers, the interlayer coupling varies with the twist angle , , .
The model above can be readily generalized to the chiral stack of N identical layers. Assuming the hopping of excitons between the nearest-neighbor layers only, we shall use the exciton Hamiltonian with the matrix elements
with n, n′=1, …, N.
Now we consider the interaction of excitons with the electromagnetic field. Because of high symmetry of an individual layer, the strength of exciton-photon coupling in a layer is described by a single parameter Γ0, which determines the exciton radiative decay rate. For the incident light propagating along the z axis (Figure 1), the exciton polarizations in the layers
The amplitude of a circularly polarized electromagnetic wave transmitted through the stack is given by
By solving Eq. (3), we calculate
Consider now the transmission of linearly polarized light through the chiral stack. The spectrum of transmission is determined by
Because of the difference in the transmission coefficients t(+) and t(−), the linearly polarized light, when transmitted through the stack, rotates its polarization plane (optical activity) and acquires ellipticity (circular dichroism). The corresponding rotation angle α and ellipticity angle ϵ are defined by
At |t(+)−t(−)|≪|t(±)|, the rotation and ellipticity angles are small, and Eq. (6) can be rewritten in the compact form
Below, we use these relations to analyze optical activity and circular dichroism in bilayers and multilayer stacks.
3 Results and discussion
3.1 Chiral bilayer
In a bilayer, the interlayer hopping of excitons leads to the formation of the symmetric and antisymmetric exciton modes denoted by the numbers m=1 and m=2, respectively . Each of them is two-fold degenerate in the spin index s. Diagonalization of the Hamiltonian (1) yields the eigenfrequencies of the modes
and the eigenfunctions
Figure 2A shows the transmission spectra of linearly polarized light for the chiral bilayers. Red and blue curves correspond to the cases of low and high nonradiative exciton decay rate Γ. The frequencies of the symmetric and antisymmetric exciton modes are indicated by vertical lines. The transmission spectra feature a single dip at the frequency of the symmetric mode with a barely noticeable trace of the antisymmetric mode.
At J≫Γ0, the symmetric and antisymmetric modes are well separated in frequency and do not interact with each other. It follows from Eq. (3) that the interaction of an individual exciton mode (s, m) with the circularly polarized light propagating along the z axis (→) or in the opposite direction (←) is characterized by the rates
The total radiative decay rate of excitons in the mode (s, m) is
Figure 2B and C show the frequency dependence of the rotation angle α(ω) and the ellipticity angle ϵ(ω), respectively. In contrast to the transmission spectra (Figure 2A), the spectra of optical rotation (Figure 2B) and circular dichroism (Figure 2C) comprise two resonances of the same strength at the frequencies of both eigen exciton modes. Polarization conversion at the frequency of the dark mode is as effective as that at the frequency of the bright mode, whereas the bilayer transparency at the dark mode frequency is very high.
This striking result can be explained as follows. The efficiency of polarization conversion at the exciton mode m is determined by the difference of the strengths of the exciton coupling to the right-handed and left-handed circularly polarized light propagating in the same direction:
These values for the bright and dark exciton modes in the chiral bilayer have the form
which follows from Eq. (11). As a result, the bright and dark exciton modes are revealed in the spectra of optical rotation and circular dichroism as the resonances of equal strengths and opposite signs.
Analytical solution of Eq. (3) shows that the rotation and ellipticity angles in the relevant case of qd≪1 are given by
The angles are proportional to the chiral hopping parameter Jsinϕ, which decays exponentially with the increase of the interlayer distance. Therefore, the interlayer distance should not exceed a few TMDC lattice constants. The widths of the resonance features in α(ω) and ϵ(ω) are determined by the nonradiative decay rate Γ. In high-quality structures with Γ≪Γ0, the resonances are very sharp and the rotation and ellipticity angles reach |α|max, |ϵ|max~(Γ0/Γ)qd|sinϕ|. In the case of Γ≫J, the resonances overlap and partly compensate each other leading to the decrease of the rotation angle and the ellipticity angle by the factor J/Γ.
3.2 Multilayer stacks
The stack consisting of N layers supports N spatial exciton modes. Diagonalization of the Hamiltonian (2) yields the eigenfrequencies of the modes
and the eigenfunctions Φ(s,m) with the elements
where the index m enumerates the spatial modes and the index s stands for the spin projection. Neglecting the phase factor, the functions Φ(s,m) are either even or odd with respect to the center of the stack. The eigen exciton modes in the stack of N=5 layers are sketched in the inset of Figure 3A.
The modes are coupled to the electromagnetic field differently. The rates of exciton-photon interaction for the modes (s, m) calculated after Eq. (10) have the form
The radiative decay rates of excitons in the modes with even m is much lower:
In the case of noninteracting modes, ω(m+1)−ω(m)≫ γ(s,m)+γ(s,m+1), the amplitude transmission coefficient of circularly polarized radiation through the stack is the sum of resonant contributions stemming from individual exciton modes
Using the definition Eq. (7), we obtain the expression for the rotation and ellipticity angles:
Figure 3 shows the spectra of transmission, optical activity, and circular dichroism of the chiral stacks of N=5 layers. Red and blue curves correspond to the stacks with low and high nonradiative exciton decay rates Γ. The transmission spectrum at low Γ (red curve in Figure 3A) features strong dips at the frequencies of the bright exciton modes (with odd m) and weak dips at the frequencies of the dark modes (with even m). The strengths of the resonances of both kinds decrease with the increase of the mode index m. At high Γ, the resonances associated with individual exciton modes are widen, and the fine structure of the transmission spectrum is not resolved.
The spectral dependences of the rotation angle (Figure 3B) and the ellipticity angle (Figure 3C) consist of five resonant contributions at the frequencies of eigen exciton modes ω(m). The resonances are well resolved in the case of low nonradiative exciton decay rate Γ. The strengths of the resonances are determined by
3.3 Thick stacks
Figure 4 shows the spectra of transmission, optical activity, and circular dichroism for the chiral stack of N=10 layers. In stacks with large N, the frequencies of eigen exciton modes, indicated by gray vertical lines in Figure 4, fill the miniband from ωx−2J to ωx+2J. At large enough N, the frequency separation between the neighboring modes Δω~4J/N becomes smaller than the broadening Γ, and the individual resonances are overlapped. The resulting angles of optical rotation and ellipticity are determined by the exciton modes with the frequencies close to the miniband bottom ωb=ωx−2J; see blue curves in Figure 4C and D.
To describe the optical properties of a thick stack, we replace the set of Eq. (3) for the exciton polarizations in individual layers
Equation (23) together with the Maxwell equation
The solution of Eqs. (23) and (25) in the bulk of the stack has the form P(s), E(s) ∝ exp(iQz) with the dispersion ω(Q) given by
where Q is the exciton-polariton wave vector.
Neglecting the mass term in the exciton dispersion and considering the Q-linear spin-dependent term as a small correction, we obtain
Figure 5 sketches the dispersions of exciton-polaritons with the spin projections s=±1 along the wave vector (solid red and blue curves, respectively) given by Eq. (26) for a lossless structure. The exciton-polariton dispersion is formed as a result of the avoided crossing of the light dispersion (dashed black line) and the exciton dispersion (dashed red and blue curves). The exciton spin-orbit interaction described by the last term in the Hamiltonian (24) leads to the spin-orbit splitting of the polariton branches. The splitting is particularly strong for the lower polariton branch at high wave vectors, i.e. at the frequencies close to ωb, where the original exciton splitting is strong and the exciton contribution to the polariton state is high. The spin-orbit splitting of polariton branches leads to a difference in the transmission coefficients for the right-handed and left-handed radiation.
To calculate the transmission coefficients, we solve Eqs. (23) and (25) with the polariton dispersion (27) and the boundary conditions of the continuity of the functions E(s) and dE(s)/dz at the front and back surfaces of the stack. Because the major contribution to the rotation and ellipticity angles is proportional to the stack thickness, small possible spin-dependent corrections to the boundary conditions can be neglected. The calculation yields
where L=(N−1)d is the stack thickness.
Finally, for the rotation and ellipticity angles in thick chiral stacks, we obtain
The rotation and ellipticity angles grow linearly with the stack thickness. The most pronounced conversion of the light polarization occurs for the light frequencies close to ωb, where the polariton spin-orbit splitting is strongest; see Figure 5. The dependences α(ω) and ϵ(ω) calculated after Eq. (31) for the stack of 10 layers are shown in Figure 4B and C, respectively, by green dashed curves. One can see that the analytical dependences agree well with the results of exact numerical calculations (blue curves). With the further increase of the number of layers N or the nonradiative decay rate of excitons Γ, the agreement becomes even better. An estimation after Eq. (31) yields α/L~102 rad/μm and α/L~0.3 rad/μm for stacks with low nonradiative decay rate Γ=0.1 meV and high nonradiative decay rate Γ=2 meV, respectively, and the other parameters Γ0=0.3 meV , J=2 meV, qd=0.01, and ϕ=0.1.
4 Polarized PL of chiral stacks
Finally, we discuss the PL of chiral stacks and show that the PL is circularly polarized. To calculate the PL spectrum, we apply the approach developed in Ref.  for multiple quantum well structures. This approach suggests that the PL is caused by random sources of excitons in the layers. The sources of exciton polarization in the layers are supposed to be identical and incoherent.
In this model, the exciton polarizations in the layers
The exciton generations in the layers are independent; therefore, the correlation function has the form
The amplitude of the electric field of the radiation emitted by the stack in the +z direction is given by Eq. (4) with
Following the procedure described in the previous section, we calculate the exciton polarizations [Eq. (32)] in the layers and then the PL spectrum [Eq. (33)]. Particularly, in the case of noninteracting exciton modes in the multilayer structure, the spectral density of the PL with the certain circular polarization s assumes the form
Figure 6A–C show the total PL spectra
Figure 6D–F show the spectral dependence of the PL circular polarization degree defined by
The PL is circularly polarized because of the chiral stacking of the layers. In structures with low nonradiative decay rate Γ (red curves), the degree of circular polarization at the frequencies of dark exciton modes is much higher than that at the frequencies of bright modes. The reason is that Pc of the radiation emitted by the m-th exciton mode is determined by the ratio
To summarize, we have developed a microscopic theory of optical activity and circular dichroism in chiral stacks of 2D crystals, such as the layers of transition metal dichalcogenides, in the spectral range of exciton transitions. The theory takes into account spin-dependent transport of excitons between the layers of the stack and coupling of excitons to an electromagnetic field. We have shown that the frequency dependence of the optical rotation angle and the ellipticity angle exhibits complex behavior reflecting the structure of eigen exciton modes in chiral stacks. In stacks with low nonradiative decay rate of excitons, the spectra of optical activity and circular dichroism comprise sharp resonances associated with individual exciton modes. In stacks with high nonradiative decay rate, the individual resonances overlap and partially compensate each other. The spectra of optical activity and circular dichroism in thick chiral stacks is well described by the developed analytical theory of exciton-polaritons with the effective spin-orbit coupling. We have also calculated the spectra of exciton PL of chiral stacks in the conditions of nonresonant pumping and shown that the PL is circularly polarized. While the consideration was focused on chiral stacks of identical layers, we expect similar polarization optical effects to occur in chiral van der Waals heterostructures in the spectral range of interlayer excitons.
We acknowledge fruitful discussions with A. N. Poddubny, E. L. Ivchenko, and L. E. Golub. This work was supported by the Government of the Russian Federation (contract # 14.W03.31.0011 at the Ioffe Institute of RAS). A. V. P. also acknowledges the support from the RF president grant SP-2912.2016.5 and the “Basis” Foundation.
Goldstein DH. Polarized light, 3rd ed. Boca Raton, FL, Taylor and Francis Group, 2011.
Kaminsky W. Experimental and phenomenological aspects of circular birefringence and related properties in transparent crystals. Rep Prog Phys 2000;63:1575–640.
Bose JC. On the rotation of plane of polarisation of electric waves by a twisted structure. Proc R Soc London 1898;63:146–52.
Wang J, Shen Z, Wu W. Broadband and high-efficiency circular polarizer based on planar-helix chiral metamaterials. Appl Phys Lett 2017;111:113503.
Chadha AS, Zhao D, Zhou W. Comparative study of metallic and dielectric helix photonic metamaterial. Opt Mater Express 2014;4:2460–7.
Kizel VA, Krasilov YI, Burkov VI. Experimental studies of gyrotropy of crystals. Sov Phys Usp 1975;18:745–73.
Mason SF. Molecular optical activity and the chiral discriminations. Cambridge, UK, Cambridge University Press, 1982.
Agranovich VM, Ginzburg V. Crystal optics with spatial dispersion, and excitons (Springer Series in Solid-State Sciences). Berlin, Springer, 2013.
Ivchenko EL, Sel’kin AV. Natural optical activity in semiconductors with wurtzite structure. Sov Phys JETP 1979;49:933–42.
Rogacheva AV, Fedotov VA, Schwanecke AS, Zheludev NI. Giant gyrotropy due to electromagnetic-field coupling in a bilayered chiral structure. Phys Rev Lett 2006;97:177401.
Lobanov SV, Weiss T, Gippius NA, et al. Polarization control of quantum dot emission by chiral photonic crystal slabs. Opt Lett 2015;40:1528–31.
Demenev AA, Kulakovskii VD, Schneider C, et al. Circularly polarized lasing in chiral modulated semiconductor microcavity with GaAs quantum wells. Appl Phys Lett 2016;109:171106.
Mamonov EA, Kolmychek IA, Vandendriessche S, et al. Anisotropy versus circular dichroism in second harmonic generation from fourfold symmetric arrays of G-shaped nanostructures. Phys Rev B 2014;89:121113.
Plum E, Liu X-X, Fedotov VA, Chen Y, Tsai DP, Zheludev NI. Metamaterials: optical activity without chirality. Phys Rev Lett 2009;102:113902.
Ibbotson LA, Demetriadou A, Croxall S, Hess O, Baumberg JJ. Optical nano-woodpiles: large-area metallic photonic crystals and metamaterials. Sci Rep 2015:5:8313.
Minovich AE, Miroshnichenko AE, Bykov AY, Murzina TV, Neshev DN, Kivshar YS. Functional and nonlinear optical metasurfaces. Laser Photonics Rev 2015;9:195–213.
Collins JT, Kuppe C, Hooper DC, Sibilia C, Centini M, Valev VK. Chirality and chiroptical effects in metal nanostructures: fundamentals and current trends. Adv Opt Mater 2017;5:1700182.
Kim C-J, Sánchez-Castillo A, Ziegler Z, Ogawa Y, Noguez C, Park J. Chiral atomically thin films. Nat Nanotechnol 2016;11:520.
Mak KF, Lee C, Hone J, Shan J, Heinz TF. Atomically thin MoS2: a new direct-gap semiconductor. Phys Rev Lett 2010;105:136805.
Mak KF, Shan J. Photonics and optoelectronics of 2D semiconductor transition metal dichalcogenides. Nat Photon 2016;10:216–26.
Cadiz F, Courtade E, Robert C, et al. Excitonic linewidth approaching the homogeneous limit in MoS2-based van der waals heterostructures. Phys Rev X 2017;7:021026.
Jun X, Mervin Z, Yuan W, Xiang Z. Excitons in atomically thin 2D semiconductors and their applications. Nanophotonics 2017;6:1309.
Arora A, Koperski M, Nogajewski K, Marcus J, Faugeras C, Potemski M. Excitonic resonances in thin films of WSe2: from monolayer to bulk material. Nanoscale 2015;7:10421–9.
Koperski M, Molas MR, Arora A, et al. Optical properties of atomically thin transition metal dichalcogenides: observations and puzzles. Nanophotonics 2017;6:1289.
Molas MR, Nogajewski K, Slobodeniuk AO, Binder J, Bartos M, Potemski M. The optical response of monolayer, few-layer and bulk tungsten disulfide. Nanoscale 2017;9:13128–41.
Liu K, Zhang L, Cao T, et al. Evolution of interlayer coupling in twisted molybdenum disulfide bilayers. Nat Commun 2014;5:4966.
Huang S, Ling X, Liang L, et al. Probing the interlayer coupling of twisted bilayer MoS2 using photoluminescence spectroscopy. Nano Lett 2014;14:5500–8.
van der Zande AM, Kunstmann J, Chernikov A, et al. Tailoring the electronic structure in bilayer molybdenum disulfide via interlayer twist. Nano Lett 2014;14:3869–75.
Yan J, Xia J, Wang X, et al. Stacking-dependent interlayer coupling in trilayer MoS2 with broken inversion symmetry. Nano Lett 2015;15:8155–61.
Plechinger G, Mooshammer F, Castellanos-Gomez A, Steele GA, Schüller C, Korn T. Optical spectroscopy of interlayer coupling in artificially stacked MoS2 layers. 2D Materials 2015;2:034016.
Xia J, Yan J, Shen ZX. Transition metal dichalcogenides: structural, optical and electronic property tuning via thickness and stacking. FlatChem 2017;4:1–19.
Arora A, Drüppel M, Schmidt R, et al. Interlayer excitons in a bulk van der Waals semiconductor. Nat Commun 2017;8:639.
Nagler P, Plechinger G, Ballottin MV, et al. Interlayer exciton dynamics in a dichalcogenide monolayer heterostructure. 2D Mater 2017;4:025112.
Ivchenko EL. Optical spectroscopy of semiconductor nanostructures. Harrow, UK, Alpha Science International, 2005.
Ivchenko EL, Nesvizhskii AI, Jorda S. Bragg reflection of light from quantum-well structures. Phys Solid State 1994;36:1156.
Poshakinskiy AV, Poddubny AN, Tarasenko SA. Reflection of short polarized optical pulses from periodic and aperiodic multiple quantum well structures. Phys Rev B 2012;86:205304.
Kazanov DR, Poshakinskiy AV, Shubina TV. Resonant photonic crystals based on van der waals heterostructures for effective light pulse retardation. Superlattices Microstruct 2017;112:639–43.
The eigen exciton modes are strictly symmetric or antisymmetric only at ϕ=0.
Robert C, Lagarde D, Cadiz F, et al. Exciton radiative lifetime in transition metal dichalcogenide monolayers. Phys Rev B 2016;93:205423.
Shahnazaryan V, Morina S, Tarasenko SA, Shelykh IA. Spin currents of exciton polaritons in microcavities with (110)-oriented quantum wells. Phys Rev B 2015;92;155305.
Voronov MM, Ivchenko EL, Erementchouk MV, Deych LI, Lisyansky AA. Photoluminescence spectroscopy of one-dimensional resonant photonic crystals. J Lumin 2007;125:112–7.