Symmetry-Based Design Rules for Second-Harmonic Generation in Stacked and Twisted MoS2 Bilayers

Caterina Cocchi; Sumanti Patra

REVIEW 4 major objections 3 minor 46 references

Symmetry-Based Design Rules for Second-Harmonic Generation in Stacked and Twisted MoS2 Bilayers

T0 review · 4 major / 3 minor · reviewed 2026-08-01 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read The twist angle of a MoS2 bilayer can be read directly from the rotation of its second-harmonic generation lobes.

desk verdict Solid symmetry analysis and one clean numerical observation, but the θ_tw = 2φ0 rule is single-angle empiricism wearing a universal-design-rule costume. read the letter →

arxiv 2607.19937 v1 pith:5WM2NO4I submitted 2026-07-22 cond-mat.mtrl-sci

Sumanti Patra , Caterina Cocchi This is my paper · ORCID

classification cond-mat.mtrl-sci

keywords second-harmonicgenerationMoS2bilayertwistanglepoint-groupsymmetryNeumann'sprincipleχ(2)tensortransitionmetaldichalcogenides

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper argues that the second-harmonic generation (SHG) response of a bilayer MoS2 is fully controlled by the stacking point group, and that twisting the layers activates a new in-plane susceptibility component that was silent in all untwisted stackings. That activation rigidly rotates the six-lobed SHG polar pattern by an angle equal to half the twist angle, and the rotation does not change with photon energy across the computed optical spectrum. If correct, SHG becomes a direct, wavelength-independent optical ruler for twist angles, alongside a set of design rules that connect stacking configuration to allowed and forbidden nonlinear tensor components. The argument combines group-theoretic tensor analysis with first-principles spectra for five configurations spanning four point groups.

What carries the argument

The central object is the second-order susceptibility tensor χ^(2) — the coefficient relating an induced polarization at twice the optical frequency to the product of two fundamental fields — constrained by Neumann's principle: every symmetry operation of the crystal point group must leave the tensor invariant. The decisive symmetry element is the vertical mirror plane σv; in all untwisted stackings it enforces χ_xxx = 0, while in the twisted C3 configuration it is absent, so χ_xxx becomes an independent component. The ratio χ_xxx/χ_yyy enters a trigonometric phase φ0 in the SHG polar pattern, and the calculations show this phase equals half the structural twist angle across the whole spectr

What would settle it

Measure the co-polarized SHG polar pattern of a MoS2 twisted bilayer with an independently known twist angle (e.g., by electron diffraction) at two or more fundamental wavelengths, and check whether the node/lobe rotation is exactly θ_tw/2 and identical at both wavelengths. Alternatively, repeat the first-principles calculation in a larger or symmetrized supercell that includes the C2′ axes: if the lobe shift changes or becomes frequency-dependent, the central claim is refuted.

Watch

Extended reading notes

Core claim

Using Neumann's principle, the paper reduces the 27-component χ^(2) tensor to the allowed independent components for the point groups D3h (monolayer/AA), D3d (2H), C3v (3R), and C3 (twisted). In the twisted C3 bilayer the in-plane components χ_xxx and χ_yyy are both independent and non-zero, which turns the co-polarized SHG intensity into sin²(3(θ−φ0)) with φ0 = (1/3) arctan(χ_xxx/χ_yyy). First-principles results for a 21.8° twisted bilayer give φ0 = 10.9°, exactly half the twist angle, at every photon energy considered; the paper calls this a geometric locking of the lobe shift to the twist.

Load-bearing premise

The load-bearing premise is that a 42-atom supercell with exact C3 symmetry faithfully represents the infinite twisted bilayer — an ideal macroscopic twisted homobilayer would also possess twofold in-plane rotation axes that force out-of-plane tensor components to zero — and that the independent-particle approximation preserves the frequency independence of the in-plane ratio.

Editorial extensions

If this is right

In twisted MoS2 bilayers, the SHG polar pattern keeps its six-lobed shape but rotates by half the twist angle, so measuring the lobe orientation gives the twist angle directly.
The rotation is independent of excitation wavelength, so the probe does not require tuning to specific resonances.
Out-of-plane tensor components activated by broken horizontal mirror symmetry can be detected with oblique-incidence p-polarized light, distinguishing 3R and twisted bilayers from monolayers and AA stacks.
AB(2H) stacking has an inversion center and produces identically zero electric-dipole SHG, so SHG can also distinguish stacking configurations.
Because the rules are set by point group alone, they should transfer to other 2H TMD bilayers and related van der Waals crystals.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

A direct test of the locking relation not reported in the paper: compute or measure φ0 at other commensurate twist angles; if φ0 = θ_tw/2 holds for 9.4°, 13.2°, or 38.2°, the relation is universal, and if it does not, the 21.8° result may be specific to that cell.
The paper's own caveat about the missing C2′ axes suggests the out-of-plane components in twisted bilayers may be an artefact of the small supercell; the in-plane ratio driving φ0 could still survive in a truly D3-symmetric macroscopic layer, but that is an open question the authors do not settle.
Since the independent-particle approximation omits excitons, a natural next step is to check whether many-body corrections preserve the frequency independence of χ_xxx/χ_yyy; excitonic renormalization could alter peak positions without necessarily changing the geometric ratio.
The symmetry logic should extend beyond MoS2 to other 2H TMDs and to local moiré domains, suggesting SHG lobe mapping as a way to image spatial variations of twist angle in heterostructures.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

The paper is worth a serious referee, but the central rule is overclaimed. The group-theory part is clean and correct: the tensor tables for D3h, D3d, C3v, and C3 are standard, and the DFT spectra confirm the symmetry-allowed components. That confirmation is useful as a reference, even if not surprising. The genuinely new observation is that in their 21.8° C3 supercell, φ0(ω) stays locked at θ_tw/2 across the whole computed spectrum. If that holds at other angles, it would be a genuinely handy calibration rule for SHG-based twist metrology. The oblique-incidence discussion for accessing out-of-plane components is also clear and useful.

The soft spots are the universality claim. The φ0 = θ_tw/2 relation rests on a single twist angle, a single stacking, and the independent-particle approximation. There are no error bars, no second-angle test, and no analytical derivation. More importantly, the paper itself notes that the ideal macroscopic twisted bilayer would have C2′ axes that are absent in the C3 supercell. The stress-test note pushes this correctly: those C2′ axes would force χ_yyy = 0 in D3 symmetry, which changes the in-plane tensor and the lobe shift entirely. The paper never addresses why the C3 supercell is the representative case for real twisted bilayers. That does not invalidate the DFT result for the specific C3 cell they studied, but it does mean the headline — a wavelength-independent optical route for twist-angle determination — is not established as a universal design rule.

The symmetry tables alone are textbook-neat, and the paper is clearly written. I would cite it for the C3 tensor structure and the frequency-flat φ0 observation, but I would not build anything on the 2φ0 relation without more angles and an explicit treatment of the C2′ problem. A serious referee should see this; the correct outcome is major revision, not rejection.

Referee Report

4 major / 3 minor

Summary. The paper develops a group-theoretic framework for the second-order susceptibility tensor of MoS2 in four point groups (D3h, D3d, C3v, C3), and validates the symmetry-allowed components with independent-particle DFT calculations and Wannier interpolation. For the C3 twisted bilayer at θ_tw = 21.8°, the authors find that the activation of χ_xxx rotates the SHG polar-lobe pattern by φ0 ≈ 10.9°, and they claim the universal relation θ_tw = 2φ0 that is frequency-independent across the optical spectrum, proposing SHG as a wavelength-independent twist-angle probe. The Neumann-principle reductions per point group are internally consistent, but the universal and frequency-independent claims are not adequately supported.

Significance. If the central relation θ_tw = 2φ0 were established, it would provide a simple, non-destructive optical route for twist-angle metrology in vdW heterostructures, with broad applicability to TMDs. The paper also offers a useful systematic classification of χ^(2) tensor structures across stacking symmetries and verifies these selection rules numerically, including the complete suppression in centrosymmetric 2H stacking and the activation of out-of-plane components when σ_h is broken. The data are shared on Zenodo, and the computational workflow is standard. However, the headline universal law rests on a single C3 supercell and ignores the in-plane consequences of the C2′ axes that the authors themselves attribute to the idealized twisted homobilayer. These issues undermine the generality of Eq. (30), although the underlying tensor tables and symmetry analysis remain valuable.

major comments (4)

[Sec. IV A, Eq. (10), Eq. (30)] The authors acknowledge that an idealized macroscopic twisted homobilayer has in-plane C2′ axes, but they apply the constraint only to out-of-plane components. For point group D3, a C2′ rotation about an in-plane x axis (R=diag(1,-1,-1)) also constrains the in-plane tensor: components with an odd number of y indices vanish, so the C3 relations in Eq. (10) reduce to χ_yyy = χ_yxx = χ_xxy = χ_xyx = 0 and χ_xxx = -χ_xyy = -χ_yxy = -χ_yyx ≠ 0. In a lab frame with the C2′ axis at β = θ_tw/2 from the bottom zigzag axis, this gives χ_xxx/χ_yyy = -cot(3β), and Eq. (22) yields φ0 = β - 30°, which for θ_tw = 21.8° is -19.1°, not +10.9°. The 42-atom supercell is explicitly C3, so the computed +10.9° and activated out-of-plane components are consequences of the reduced cell symmetry, not of the ideal twisted bilayer. Eq. (30) cannot be asserted as universal without a D3-preserving calculation or a c
[Eqs. (19)–(22)] The compact phase-shifted form Eq. (21) assumes that χ_xxx/χ_yyy is real. For complex tensor components, I_∥(θ) = |χ_yyy sin3θ - χ_xxx cos3θ|² contains the cross term -2 Re(χ_xxx χ_yyy*) sin3θ cos3θ and does not factor into (|χ_yyy|²+|χ_xxx|²) sin²(3(θ-φ0)) with a real φ0; Eq. (22) is then undefined. The manuscript neither states this reality condition nor provides the relative phase of the computed components. Since the six-lobed nodal pattern and the interpretation of φ0 as a rigid lobe shift depend on this factorization, this missing step is load-bearing.
[Sec. IV C, Fig. 4(d), Eq. (30)] The universal relation θ_tw = 2φ0 is validated at a single commensurate angle, 21.8°, from one C3 supercell. No second twist angle, no incommensurate cell, and no error bar or convergence estimate are reported. Deviations in Fig. 4(d) are described only qualitatively. The statement that the ratio is “solely dictated by the spatial projection of the atomic coordinates” is an assertion; without an analytical derivation or multiple test angles, it cannot support a material-agnostic law.
[Sec. III and Sec. IV C] The frequency-independence claim covers only the independent-particle DFT window shown in Figs. 2 and 4 (roughly 0.6–1.6 eV), not “the entire optical spectrum” as stated in the Abstract. Many-body effects are neglected; excitonic and self-energy corrections can renormalize the relative spectral weights of χ_xxx and χ_yyy and could make φ0 frequency-dependent. The claims should be re-scoped or supported by calculations beyond the independent-particle approximation.

minor comments (3)

[Sec. II C] Typo: “although the become accessible” should read “although they become accessible.”
[Appendix A 2] The text says C3 leaves “13 unique independent tensor components”; the later C3 tensor structure in Eq. (10) contains fewer independent parameters after permutation and C3 relations. Please clarify the counting convention.
[Fig. 4(d)] The vertical axis label “Peak shift ¢peak” is unclear; it should be φ0(ω) or a related notation. Also, the color-coding by normalized |χ_yyy| should be explained in the caption.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the symmetry analysis and the first-principles SHG calculations are independent, and no fitted parameter or self-citation is load-bearing.

full rationale

The derivation chain is self-contained. The tensor structures in Table II are derived from Neumann's principle applied to the point groups, so they are independent of the DFT results and are not fitted. The C3 lobe-shift formula, Eq. (22), is an algebraic identity obtained by projecting the polarization, and it is then evaluated using the DFT-computed independent tensor components χ_xxx and χ_yyy, which are never adjusted to enforce θ_tw/2. The central numerical result, φ0 = 10.9° for θ_tw = 21.8°, is therefore an output of the ab initio calculation rather than an input; it could in principle have taken a different value, and the paper even notes frequency-dependent fluctuations where the tensor magnitudes approach zero. There are no self-citations that carry the argument, no fitted parameters relabeled as predictions, and no ansatz imported from prior work by the same authors. The acknowledged limitations—the absence of C2′ axes in the 42-atom supercell and the use of the independent-particle approximation—are physical/correctness concerns about whether the supercell represents the ideal macroscopic twisted bilayer, not circularity: they do not make any equation equal to its input by construction.

Assumptions & free parameters 3 free parameters · 5 assumptions · 0 invented entities

No experimentally fitted parameters; the free parameters are numerical settings and the single tested twist angle. The load-bearing assumptions are the independent-particle approximation, the representativeness of the finite supercell, and the standard symmetry principles.

free parameters (3)

Lorentzian broadening η = 0.01 eV
Chosen for spectral resolution; not fitted to experiment, but affects the width of computed resonances and the noise floor of φ0(ω).
Twist angle θ_tw = 21.8°
Single commensurate angle used for the twisted supercell; the claim φ0=θ_tw/2 is only tested at this value.
Evaluation energies for polar patterns = 1.31 eV (1L), 1.28 eV (AA), 1.30 eV (tw)
Angular patterns evaluated at selected C-peak resonance energies from the computed spectra; these choices do not affect the symmetry of the pattern but set the amplitudes.

assumptions (5)

domain assumption Neumann's principle: χ^(2) is invariant under every operation of the crystal point group (Eq. 4).
Standard principle of crystal physics; not proved in the paper.
standard math Intrinsic permutation symmetry χ^(2)_ijk = χ^(2)_ikj for SHG (Appendix A1).
Follows from indistinguishability of the two photons.
domain assumption The independent-particle approximation (no excitonic effects) is sufficient to determine the tensor component ratio χ_xxx/χ_yyy and its frequency dependence.
All χ^(2) spectra are computed with the IP approximation (Sec. III). The frequency-independence claim of φ0 rests on this.
domain assumption The 42-atom supercell of the θ=21.8° twisted bilayer correctly represents the symmetry and the SHG tensor of the infinite twisted bilayer.
The paper itself notes the supercell lacks C2′ axes that an idealized macroscopic twisted bilayer would have (Sec. IV A), which would forbid χ_zxx and χ_xxz.
domain assumption Wannier interpolation (Wannier90/postw90) yields converged χ^(2) spectra on the dense 108×108×1 k-mesh.
Validation described only as band-structure match in Supplemental Material.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Symmetry-Based Design Rules for Second-Harmonic Generation in Stacked and Twisted MoS2 Bilayers." pith.science (2026). https://pith.science/paper/5WM2NO4I

@misc{pith2026260719937,
  author       = {Pith},
  title        = {Pith review of: Symmetry-Based Design Rules for Second-Harmonic Generation in Stacked and Twisted MoS2 Bilayers},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/5WM2NO4I}},
  note         = {Machine review of arXiv:2607.19937}
}

read the original abstract

Understanding how stacking controls the nonlinear optical response of two-dimensional materials is key to designing van der Waals heterostructures with tailored functionalities. Here, we establish a comprehensive symmetry-based framework mapping the structural configuration of MoS2 bilayers across four point groups (D3h, D3d, C3v, C3) to their second-order susceptibility tensor chi^(2). Using group-theory arguments benchmarked against first-principles calculations, we demonstrate how symmetry breaking controls the activation and suppression of individual tensor elements in these systems. We show that the emergence of the in-plane component chi_xxx in twisted configurations (C3 group) induces a rigid azimuthal rotation of the second-harmonic generation polar lobes, which remains frequency-independent across the entire optical spectrum, locking to half of the structural twist angle. Our findings establish a direct, wavelength-independent optical route for twist-angle determination and provide a clear roadmap for engineering nonlinear optical responses in two-dimensional materials.

Figures

Figures reproduced from arXiv: 2607.19937 by the authors.

**Figure 2.** FIG. 2. Frequency-dependent susceptibility spectra [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗

**Figure 3.** FIG. 3. Computed [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗

Figures from the paper (1 more)

**Figure 4.** Figure 4: FIG. 4. Co-polarized SHG polar patterns [PITH_FULL_IMAGE:figures/full_fig_p010_4.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

46 extracted references

[1]

(A2) This constraint reduces the 27 initial components to 3× 6 = 18 independent elements

Intrinsic Permutation Symmetry In the SHG experimental setting, two incoming pho- tons oscillate at the identical frequencyω, making the last two indices ofχ (2) ijk completely interchangeable: χ(2) ijk =χ (2) ikj .(A1) For any fixed indexi∈ {x, y, z}, the remainingjkpair possesses 3×3 = 9 permutations, only 6 of which are linearly independent underj↔kexc...
[2]

Reduction byC 3 Rotation The threefold rotation operationC 3 around the out- of-planez-axis represents a counterclockwise rotation by 120◦ about thez-axis, defined by the standard transfor- mation matrix: R(C3) =   − 1 2 − √ 3 2 0√ 3 2 − 1 2 0 0 0 1   .(A3) Neumann’s principle requires that the tensor must re- main invariant under any valid symmetry o...
[3]

Horizontal Mirror Planeσ h The horizontal mirror planeσ h reflects spatial coordi- nates across the basal plane, mappingx→x,y→y, z→ −z, with matrixR(σ h) = diag(1,1,−1). Under this transformation, eachzindex (n z) contributes a factor of −1: χ(2) ijk σh − − →(−1)nz χ(2) ijk .(A12) Applying Neumann’s principle then forcesχ (2) ijk = 0 for any combination w...
[4]

Here, eachyindex (n y) contributes a factor of−1: χ(2) ijk σv − →(−1)ny χ(2) ijk .(A15) Components with an odd value ofn y vanish

V ertical Mirror Planeσv The vertical mirror planeσ v reflects coordinates across a vertical plane perpendicular to the layer, mappingx→ x,y→ −y,z→z, expressed asR(σ v) = diag(1,−1,1). Here, eachyindex (n y) contributes a factor of−1: χ(2) ijk σv − →(−1)ny χ(2) ijk .(A15) Components with an odd value ofn y vanish. Taking as an example the in-plane compone...
[5]

Inversion SymmetryI A macroscopic spatial inversion operation maps all co- ordinates through the origin,r→ −r, yielding a trans- formation matrixR(I) = diag(−1,−1,−1). For a third- rank polar tensor, this operation introduces an odd trans- formation parity across all indices simultaneously: χ(2) ijk = 0∀i, j, k.(A18) Enforcing Neumann’s principle requires...
[6]

Y. Li, Y. Rao, K. F. Mak, Y. You, S. Wang, C. R. Dean, and T. F. Heinz, Probing symmetry properties of few- layer MoS2 and h-BN by optical second-harmonic gener- ation, Nano Lett.13, 3329 (2013)

2013
[7]

Autere, H

A. Autere, H. Jussila, Y. Dai, Y. Wang, H. Lipsanen, and Z. Sun, Nonlinear optics with 2D layered materials, Adv. Mater.30, 1705963 (2018)

2018
[8]

Klimmer, O

S. Klimmer, O. Ghaebi, Z. Gan, A. George, A. Tur- chanin, G. Cerullo, and G. Soavi, All-optical polariza- tion and amplitude modulation of second-harmonic gen- eration in atomically thin semiconductors, Nat. Photon. 15, 837 (2021)

2021

Show all 46 references

[9]

Dogadov, C

O. Dogadov, C. Trovatello, B. Yao, G. Soavi, and G. Cerullo, Parametric nonlinear optics with layered ma- terials and related heterostructures, Laser Photon. Rev. 16, 2100726 (2022)

2022
[10]

Herrmann, S

P. Herrmann, S. Klimmer, T. Lettau, T. Weickhardt, A. Papavasileiou, K. Mosina, Z. Sofer, I. Paradisanos, D. Kartashov, J. Wilhelm,et al., Nonlinear valley selec- tion rules and all-optical probe of broken time-reversal symmetry in monolayer wse2, Nat. Photon.19, 300 (2025)

2025
[11]

J. A. Wilson and A. D. Yoffe, The transition metal dichalcogenides: Discussion and interpretation of the ob- served optical, electrical and structural properties, Adv. Phys.18, 193 (1969)

1969
[12]

Chhowalla, H

M. Chhowalla, H. S. Shin, G. Eda, L.-J. Li, K. P. Loh, and H. Zhang, The chemistry of two-dimensional layered transition metal dichalcogenide nanosheets, Nat. Chem. 5, 263 (2013)

2013
[13]

Manzeli, D

S. Manzeli, D. Ovchinnikov, D. Pasquier, O. V. Yazyev, and A. Kis, 2D transition metal dichalcogenides, Nat. Rev. Mater.2, 17033 (2017)

2017
[14]

Kumar, S

N. Kumar, S. Najmaei, Q. Cui, F. Ceballos, P. M. Ajayan, J. Lou, and H. Zhao, Second harmonic mi- croscopy of monolayer MoS 2, Phys. Rev. B87, 161403 (2013)

2013
[15]

L. M. Malard, T. V. Alencar, A. P. M. Barboza, K. F. Mak, and A. M. de Paula, Observation of intense second harmonic generation from MoS 2 atomic crystals, Phys. Rev. B87, 201401 (2013)

2013
[16]

Zhang, W

J. Zhang, W. Zhao, P. Yu, G. Yang, and Z. Liu, Second harmonic generation in 2d layered materials, 2D Mater. 7, 042002 (2020)

2020
[17]

L. Du, Z. Huang, J. Zhang, F. Ye, Q. Dai, H. Deng, G. Zhang, and Z. Sun, Nonlinear physics of moir´ e super- lattices, Nat. Mater.23, 1179 (2024)

2024
[18]

Huang, Y

W. Huang, Y. Xiao, F. Xia, X. Chen, and T. Zhai, Second harmonic generation control in 2d layered materials: sta- tus and outlook, Adv. Funct. Mater.34, 2310726 (2024)

2024
[19]

Hsu, Z.-A

W.-T. Hsu, Z.-A. Zhao, L.-J. Li, C.-H. Chen, M.-H. Chiu, P.-S. Chang, Y.-C. Chou, and W.-H. Chang, Sec- ond harmonic generation from artificially stacked transi- tion metal dichalcogenide twisted bilayers, ACS Nano8, 2951 (2014)

2014
[20]

Psilodimitrakopoulos, L

S. Psilodimitrakopoulos, L. Mouchliadis, I. Paradis- anos, G. Kourmoulakis, A. Lemonis, G. Kioseoglou, and E. Stratakis, Twist angle mapping in layered WS 2 by polarization-resolved second harmonic generation, Sci. Rep.9, 14285 (2019)

2019
[21]

Y. Wang, J. Xiao, S. Yang, Y. Wang, and X. Zhang, Second harmonic generation spectroscopy on two- dimensional materials, Opt. Mater. Express9, 1136 (2019)

2019
[22]

L. Zhou, H. Fu, T. Lv, C. Wang, H. Gao, D. Li, L. Deng, and W. Xiong, Nonlinear optical characterization of 2d materials, Nanomater.10, 2263 (2020)

2020
[23]

Paradisanoset al., Optical versus electron diffraction imaging of twist-angle in 2D transition metal dichalco- genide bilayers, npj 2D Mater

I. Paradisanoset al., Optical versus electron diffraction imaging of twist-angle in 2D transition metal dichalco- genide bilayers, npj 2D Mater. Appl.5, 59 (2021)

2021
[24]

Paradisanos, A

I. Paradisanos, A. M. Saiz Raven, T. Amand, C. Robert, P. Renucci, K. Watanabe, T. Taniguchi, I. C. Gerber, X. Marie, and B. Urbaszek, Second harmonic generation control in twisted bilayers of transition metal dichalco- genides, Phys. Rev. B105, 115420 (2022)

2022
[25]

R. W. Boyd,Nonlinear Optics, 3rd ed. (Academic Press, 2008)

2008
[26]

J. Shi, P. Yu, F. Liu, P. He, R. Wang, L. Qin, J. Zhou, X. Li, J. Zhou, X. Sui,et al., 3r mos2 with broken inver- sion symmetry: a promising ultrathin nonlinear optical device, Adv. Mater.29, 1701486 (2017)

2017
[27]

Attaccalite and M

C. Attaccalite and M. Gr¨ uning, Nonlinear optics from an ab initio approach by means of the dynamical Berry phase: Application to second- and third-harmonic gener- ation in semiconductors, Phys. Rev. B88, 235113 (2013)

2013
[28]

Attaccalite, M

C. Attaccalite, M. Palummo, E. Cannuccia, and M. Gr¨ uning, Second-harmonic generation in single-layer monochalcogenides: A response from first-principles real- time simulations, Phys. Rev. Mater.3, 074003 (2019)

2019
[29]

N. A. Pike and R. Pachter, Second-order nonlinear op- tical properties of monolayer transition-metal dichalco- genides by computational analysis, J. Phys. Chem. C 125, 11075 (2021)

2021
[30]

Ruan, Y.-H

J. Ruan, Y.-H. Chan, and S. G. Louie, Exciton enhanced nonlinear optical responses in monolayer h-bn and mos2: insight from first-principles exciton-state coupling for- malism and calculations, Nano Lett.24, 15533 (2024)

2024
[31]

Gr¨ uning and C

M. Gr¨ uning and C. Attaccalite, Second harmonic genera- tion in h-BN and MoS2 monolayers: Role of electron-hole interaction, Phys. Rev. B89, 081102 (2014)

2014
[32]

C. Wang, X. Liu, L. Kang, B.-L. Gu, Y. Yao, and W. Duan, First-principles calculation of nonlinear opti- cal responses by Wannier interpolation, Phys. Rev. B96, 115147 (2017)

2017
[33]

Garcia-Goiricelayaet al., Phys

P. Garcia-Goiricelayaet al., Phys. Rev. B107, 195406 (2023)

2023
[34]

K. F. Mak, C. Lee, J. Hone, J. Shan, and T. F. Heinz, Atomically thin MoS2: A new direct-gap semiconductor, Phys. Rev. Lett.105, 136805 (2010)

2010
[35]

Kresse and J

G. Kresse and J. Hafner, Ab initio molecular dynamics for liquid metals, Phys. Rev. B47, 558 (1993)

1993
[36]

Kresse and J

G. Kresse and J. Hafner, Ab initio molecular-dynamics simulation of the liquid-metal–amorphous-semiconductor transition in germanium, Phys. Rev. B49, 14251 (1994)

1994
[37]

J. P. Perdew, K. Burke, and M. Ernzerhof, Generalized gradient approximation made simple, Phys. Rev. Lett. 77, 3865 (1996)

1996
[38]

J. P. Perdew, K. Burke, and M. Ernzerhof, Generalized gradient approximation made simple [phys. rev. lett. 77, 3865 (1996)], Phys. Rev. Lett.78, 1396 (1997). 14

1996
[39]

P. E. Bl¨ ochl, Projector augmented-wave method, Phys. Rev. B50, 17953 (1994)

1994
[40]

Grimme, J

S. Grimme, J. Antony, S. Ehrlich, and H. Krieg, A con- sistent and accurateab initioparametrization of density functional dispersion correction (DFT-D) for the 94 ele- ments H-Pu, J. Chem. Phys.132, 154104 (2010)

2010
[41]

A. A. Mostofi, J. R. Yates, Y.-S. Lee, I. Souza, D. Van- derbilt, and N. Marzari, Wannier90: A tool for obtaining maximally-localised Wannier functions, Comput. Phys. Commun.178, 685 (2008)

2008
[42]

J. E. Sipe and A. I. Shkrebtii, Second-order optical re- sponse in semiconductors, Phys. Rev. B61, 5337 (2000)

2000
[43]

Korm´ anyos, V

A. Korm´ anyos, V. Z´ olyomi, N. D. Drummond, P. Rakyta, G. Burkard, and V. I. Fal ´ko, Monolayer MoS 2: trigo- nal warping, theγvalley, and spin-orbit coupling effects, Phys. Rev. B88, 045416 (2013)

2013
[44]

Splendiani, L

A. Splendiani, L. Sun, Y. Zhang, T. Li, J. Kim, C.-Y. Chim, G. Galli, and F. Wang, Emerging photolumines- cence in monolayer MoS 2, Nano Lett.10, 1271 (2010)

2010
[45]

Xiao, G.-B

D. Xiao, G.-B. Liu, X. Feng, Wanxiang Xu, and W. Yao, Coupled spin and valley physics in monolayers of MoS 2 and other group-VI dichalcogenides, Phys. Rev. Lett. 108, 196802 (2012)

2012
[46]

Bao, S.-Q

X. Bao, S.-Q. Li, H. Chen, Q. Chen, Y. Wu, L. Zhang, Y. Xie, and Y. Chen, Extremely anisotropic second harmonic generation in in-plane janus transition metal dichalcogenide heterostructures, Phys. Rev. B114, 045418 (2026)

2026

Pith tools

Reviewed August 1, 2026 · model on record in the stance chip above.

[1] [1]

(A2) This constraint reduces the 27 initial components to 3× 6 = 18 independent elements

Intrinsic Permutation Symmetry In the SHG experimental setting, two incoming pho- tons oscillate at the identical frequencyω, making the last two indices ofχ (2) ijk completely interchangeable: χ(2) ijk =χ (2) ikj .(A1) For any fixed indexi∈ {x, y, z}, the remainingjkpair possesses 3×3 = 9 permutations, only 6 of which are linearly independent underj↔kexc...

[2] [2]

Reduction byC 3 Rotation The threefold rotation operationC 3 around the out- of-planez-axis represents a counterclockwise rotation by 120◦ about thez-axis, defined by the standard transfor- mation matrix: R(C3) =   − 1 2 − √ 3 2 0√ 3 2 − 1 2 0 0 0 1   .(A3) Neumann’s principle requires that the tensor must re- main invariant under any valid symmetry o...

[3] [3]

Horizontal Mirror Planeσ h The horizontal mirror planeσ h reflects spatial coordi- nates across the basal plane, mappingx→x,y→y, z→ −z, with matrixR(σ h) = diag(1,1,−1). Under this transformation, eachzindex (n z) contributes a factor of −1: χ(2) ijk σh − − →(−1)nz χ(2) ijk .(A12) Applying Neumann’s principle then forcesχ (2) ijk = 0 for any combination w...

[4] [4]

Here, eachyindex (n y) contributes a factor of−1: χ(2) ijk σv − →(−1)ny χ(2) ijk .(A15) Components with an odd value ofn y vanish

V ertical Mirror Planeσv The vertical mirror planeσ v reflects coordinates across a vertical plane perpendicular to the layer, mappingx→ x,y→ −y,z→z, expressed asR(σ v) = diag(1,−1,1). Here, eachyindex (n y) contributes a factor of−1: χ(2) ijk σv − →(−1)ny χ(2) ijk .(A15) Components with an odd value ofn y vanish. Taking as an example the in-plane compone...

[5] [5]

Inversion SymmetryI A macroscopic spatial inversion operation maps all co- ordinates through the origin,r→ −r, yielding a trans- formation matrixR(I) = diag(−1,−1,−1). For a third- rank polar tensor, this operation introduces an odd trans- formation parity across all indices simultaneously: χ(2) ijk = 0∀i, j, k.(A18) Enforcing Neumann’s principle requires...

[6] [6]

Y. Li, Y. Rao, K. F. Mak, Y. You, S. Wang, C. R. Dean, and T. F. Heinz, Probing symmetry properties of few- layer MoS2 and h-BN by optical second-harmonic gener- ation, Nano Lett.13, 3329 (2013)

2013

[7] [7]

Autere, H

A. Autere, H. Jussila, Y. Dai, Y. Wang, H. Lipsanen, and Z. Sun, Nonlinear optics with 2D layered materials, Adv. Mater.30, 1705963 (2018)

2018

[8] [8]

Klimmer, O

S. Klimmer, O. Ghaebi, Z. Gan, A. George, A. Tur- chanin, G. Cerullo, and G. Soavi, All-optical polariza- tion and amplitude modulation of second-harmonic gen- eration in atomically thin semiconductors, Nat. Photon. 15, 837 (2021)

2021