REVIEW 3 major objections 4 minor 46 references
Photostriction-tunable Polarization and Structural Dynamics in Interlayer Sliding Ferroelectrics
T0 review · 3 major / 4 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read Photodoping stretches bilayer 3R-MoS2 and strengthens its sliding ferroelectric polarization.
desk verdict A solid but under-specified computational study: the abstract contradicts its own interlayer-spacing data, and the constrained-DFT occupation scheme is never given. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central objects are the rhombohedral (3R) stacking of bilayer MoS2, whose broken inversion symmetry produces a net out-of-plane sliding polarization, and the photostriction mechanism that couples photoexcited carriers to the lattice. The argument is carried by two computational routes: real-time time-dependent density functional theory, which simulates a 2 eV laser pulse and the subsequent coupled electron-ion motion, and orbital-constrained density functional theory, which freezes orbital occupations to model photodoped states and maps their total energy against the in-plane lattice constant. The load-bearing quantities are the energy-versus-lattice-constant curves at fixed carrier densities, the geometric-phase polarization and interlayer distance as functions of biaxial strain, nudged-elastic-band barriers between the two sliding stacking orders, and the real-time interlayer oscillation period of the A1g mode.
What would settle it
Measure the in-plane lattice constant of bilayer 3R-MoS2 with ultrafast electron or X-ray diffraction immediately after a 2 eV pump pulse that creates roughly 0.1 to 0.8 electrons per unit cell: the calculations predict a transient expansion of about 1% to 5%, with the interlayer A1g mode softening from about 405 cm-1 toward 369 cm-1. Seeing no expansion, or an expansion far below 1% at low fluence, would contradict the central claim.
Extended reading notes
Core claim
The central claim is that photodoping in sliding ferroelectric bilayer 3R-MoS2 produces a pronounced photostrictive response, and that the strain it creates feeds back into the ferroelectric and electronic degrees of freedom. Orbital-constrained density functional theory places the energy minimum of the photoexcited system at an in-plane lattice constant of 3.30 Å at 0.8 e/u.c., compared with 3.15 Å in the ground state, a roughly 5% expansion; at the same carrier density the geometric-phase polarization rises monotonically from 0.84 to 1.09 pC/m and the indirect gap falls from 1.30 eV to 0.22 eV. The same expansion softens the interlayer A1g breathing-like mode, lengthening its period from 76.0 fs to 86.6 fs and shifting its frequency from about 405 cm-1 to about 369 cm-1, and it raises the sliding energy barrier only modestly, from 18 to 22 meV/u.c., so polarization switching remains accessible. The mechanism is the strong electromechanical coupling inherent to the non-centrosymmetric rhombohedral stacking: photoexcited carriers redistribute interlayer charge, and the weak van der Waals interlayer bonding makes the out-of-plane response unusually sensitive to in-plane expansion.
Load-bearing premise
The load-bearing premise is that a fixed pattern of orbital occupations in constrained density functional theory faithfully represents the electronic state a real laser pulse creates, and the paper never specifies which orbitals are emptied and filled for each carrier density, so a different occupation pattern could change the predicted expansion and polarization shift.
Editorial extensions
If this is right
- The paper's results imply that a single laser pulse can dynamically shrink the bandgap of bilayer 3R-MoS2 by about 1 eV, moving the optical response from visible to near-infrared.
- The strain-enhanced polarization means an in-plane stretch can be read out as a change in ferroelectric polarization, creating a direct strain-light-memory coupling.
- Because the interlayer A1g mode softens with expansion, photostriction also tunes phonon frequencies and thus vibrational and thermal properties on ultrafast timescales.
- The sliding energy barrier stays low under 5% expansion, so ferroelectric switching remains energetically feasible while the bandgap is being renormalized, making optically programmable ferroelectric memory a plausible target.
- The contrast with 2H stacking, which shows much weaker interlayer relaxation and bandgap response under the same strain, identifies broken inversion symmetry as the ingredient that turns lattice expansion into strong ferroelectric and electronic tuning.
Reading between the lines
- By extension, the same photostriction mechanism should appear in other 3R-stacked transition metal dichalcogenides and in rhombohedral boron nitride, wherever sliding polarization exists; the paper does not compute those cases.
- A concrete experimental test would be ultrafast electron or X-ray diffraction after a 2 eV pump: the calculations predict a transient in-plane expansion of about 1% at low fluence and several percent at high fluence, together with a red-shift of the interlayer Raman mode.
- The paper's abstract states that photodoping increases interlayer spacing, while the main-text numbers show the relaxed interlayer distance decreasing from 3.06 Å to 2.95 Å under 5% expansion; the sign of this sub-effect is one place the manuscript is internally inconsistent.
- The supplemental methods acknowledge a limitation of the coupled electron-ion dynamics: nuclear motion is represented by a single averaged trajectory, which can miss multi-path effects; the authors argue the laser-driven path dominates here.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports first-principles calculations on bilayer 3R-MoS2 as a sliding ferroelectric. Using real-time TDDFT to estimate photodoping carrier densities (nph = 0.1, 0.2, 0.8 e/u.c. for increasing laser intensity) and orbital-constrained DFT (OCDFT) to compute total energies as a function of in-plane lattice constant, the authors claim a pronounced photostrictive response: in-plane expansion up to ~5% at 0.8 e/u.c., a decrease in interlayer spacing, enhanced out-of-plane polarization (0.84 to 1.09 pC/m), a strong bandgap reduction (1.30 to 0.22 eV), and phonon softening of the interlayer A1g mode. The paper further argues that these effects arise from the non-centrosymmetric rhombohedral stacking and proposes applications in optically programmable ferroelectric memories.
Significance. The work addresses an interesting and timely question: whether light-induced carrier injection can dynamically tune both the structural and ferroelectric properties of sliding ferroelectrics. The computational setup is standard for the community (DFT with PBE+optB86b-vdW, 550 eV cutoff, 21x21x1 k-grid; TDDFT with converged real-time propagation) and the inclusion of a comparison with the 2H phase is a strength. If the central claim is substantiated, the predicted 5% photostrictive strain and the associated polarization and bandgap changes would be a notable step toward optical control of slidetronics devices. However, the quantitative predictions depend on an under-specified mapping between the TDDFT carrier density and the OCDFT occupation scheme, and the internal inconsistency between the abstract and Fig. 3(a) currently undermines the phenomenological narrative.
major comments (3)
- [Supplement, Eq. (1)] The orbital-constrained DFT occupation scheme is never specified. Eq. (1) defines rho_OCDFT(r) = sum_i n_i^OCDFT phi_i^*(r) phi_i(r), but the manuscript does not state which orbitals (band index, k-point, spin) are emptied and filled to realize the target carrier densities nph = 0.1, 0.2, and 0.8 e/u.c. The real-time TDDFT simulations determine only the total number of excited electrons, not a unique static occupation pattern. Different occupation choices (e.g., excitations at specific k-points versus a uniform occupation shift, spin-conserving versus spin-flip) can lead to different charge redistributions, forces, and equilibrium lattice constants. Since the central claim of a ~5% lattice expansion at 0.8 e/u.c. is derived from the minimum of the OCDFT energy surface (Fig. 2d), the result is not reproducible without this information. Please specify the occupation pattern used or provide a sensitivity analysis over plausible occupation schemes.
- [Abstract and Fig. 3(a)] The abstract states that electron-hole excitation leads to 'substantial in-plane expansion, increased interlayer spacing, and enhanced ferroelectric polarization,' but Fig. 3(a) and the associated main-text discussion report the opposite: the interlayer spacing decreases from 3.06 Å at zero expansion to 2.95 Å at 5% expansion. This is a direct contradiction in the central physical description. The authors should correct the abstract (or, if the interlayer distance actually increases in the dynamically photoexcited state, clearly explain the difference between the static strained geometries in Fig. 3 and the TDDFT dynamics). As written, the inconsistency makes it difficult to interpret the polarization and phonon results.
- [Main text, paragraph after Fig. 3(a)] There is an internal numerical inconsistency in the reported interlayer distance values. The text first states the relaxed interlayer distance is 3.06 Å and decreases to 2.95 Å at 5% expansion (consistent with the Fig. 3(a) axis label), but later, when comparing 3R and 2H stacking, it states the 3R interlayer distance decreases 'from 3.04 Å to 2.94 Å.' These two sets of values should be reconciled, as the reported polarization and phonon softening depend on the quantitative interlayer relaxation.
minor comments (4)
- [Conclusion] Typo: 'combing' should be 'combining'.
- [Figure 2 caption] Typo: 'The stars denote he energy minima' should be 'the energy minima'.
- [Fig. 4(a) and surrounding text] The text refers to 'the fluctuation reaches 80.6 fs' and later 'the oscillation period under tensile strain,' but it is not clear whether the quoted quantities are oscillation periods or something else. Please label the quantity explicitly.
- [Abstract and main text] The abstract mentions 'increased interlayer spacing' while the main text reports a decrease; beyond the major issue, please ensure all occurrences are consistent after revision.
Circularity Check
No significant circularity: the photostrictive expansion is a first-principles response to a computed carrier density, not a fit or a self-citation.
full rationale
The derivation chain is self-contained: real-time TDDFT yields the photoexcited carrier densities nph = 0.1, 0.2, and 0.8 e/u.c.; orbital-constrained DFT then minimizes the total energy E_KS[rho_OCDFT] over in-plane lattice constants to locate the expanded minima shown in Fig. 2. The lattice expansion is therefore the output of the calculation, not an input fitted to reproduce the expansion. Polarization, interlayer distance, bandgap, and NEB barrier are computed at those expanded lattice constants from the relaxed structures. The self-references are to the TDAP code (Refs. [26], [29]) and to prior work on light-induced polarization reversal (Ref. [21]); these are method/tool citations or contextual background and are not load-bearing for the central numerical result. The main caveat, that the OCDFT occupation scheme is not specified in the paper, is a reproducibility and model-validity concern, not a circular reduction of the prediction to its own input.
Assumptions & free parameters
free parameters (3)
- Photodoping carrier density n_ph =
0.1, 0.2, 0.8 e/u.c.
- Laser field strength E0 =
0.64, 1.25, 2.50 x 10^-4 V/A
- Photon energy =
2 eV
assumptions (5)
- domain assumption Kohn-Sham DFT with the PBE exchange-correlation functional provides a sufficiently accurate description of the electronic ground state and band structure of bilayer 3R-MoS2.
- domain assumption The optB86b-vdW functional correctly captures the interlayer van der Waals interactions and the response to strain.
- ad hoc to paper Orbital-constrained DFT with fixed occupation numbers faithfully represents the non-equilibrium photoexcited state created by a laser pulse.
- domain assumption Ehrenfest dynamics (mean-field ionic trajectories) adequately captures the coupled electron-ion motion after photoexcitation.
- standard math The Berry phase method gives a well-defined polarization for the relaxed slabs.
Cite this review
Pith. "Pith review of Photostriction-tunable Polarization and Structural Dynamics in Interlayer Sliding Ferroelectrics." pith.science (2026). https://pith.science/paper/4ZY3YHN5
@misc{pith2026250524186,
author = {Pith},
title = {Pith review of: Photostriction-tunable Polarization and Structural Dynamics in Interlayer Sliding Ferroelectrics},
year = {2026},
howpublished = {\url{https://pith.science/paper/4ZY3YHN5}},
note = {Machine review of arXiv:2505.24186}
}
read the original abstract
Two-dimensional ferroelectrics with robust polarization offer promising opportunities for non-volatile memory, field-effect transistors, and optoelectronic devices. However, the impact of lattice deformation on polarization and photoinduced structural response remains poorly understood. Here, we employ first-principles calculations to demonstrate photodoping-induced lattice expansion in rhombohedrally stacked bilayer MoS2, revealing a strong coupling between photodoping carrier and lattice structure. We identify a pronounced photostrictive response in sliding ferroelectrics, wherein electron-hole excitation leads to substantial in-plane expansion, increased interlayer spacing, and enhanced ferroelectric polarization. This strain-induced modulation drives significant bandgap renormalization. The photostriction-tunable polarization and structural dynamics arise from the strong electromechanical coupling inherent to the non-centrosymmetric rhombohedral stacking. The findings provide critical insights into the nonthermal lattice expansion governing sliding ferroelectrics at atomic-scale timescales, while simultaneously laying the groundwork for next-generation electronic and memory technologies by leveraging lattice-tunable polarization switching.
Reference graph
Works this paper leans on
-
[1]
R. Bian et al. , Developing fatigue -resistant ferroelectrics using interlayer sliding switching, Science 385, 57 (2024)
work page 2024
-
[2]
Meng et al., Sliding induced multiple polarization states in two-dimensional ferroelectrics, Nat
P. Meng et al., Sliding induced multiple polarization states in two-dimensional ferroelectrics, Nat. Commun. 13, 7696 (2022)
work page 2022
-
[3]
L.-P. Miao, N. Ding, N. Wang, C. Shi, H. -Y. Ye, L. Li, Y. -F. Yao, S. Dong, and Y. Zhang, Direct observation of geometric and sliding ferroelectricity in an amphidynamic crystal, Nat. Mater. 21, 1158 (2022)
work page 2022
-
[4]
B. Qin et al., Interfacial epitaxy of multilayer rhombohedral transition-metal dichalcogenide single crystals, Science 385, 99 (2024)
work page 2024
- [5]
-
[6]
P. Tang and G. E. W. Bauer, Sliding Phase Transition in Ferroelectric van der Waals Bilayers, Phys. Rev. Lett. 130, 176801 (2023)
work page 2023
-
[7]
M. Vizner Stern, S. Salleh Atri, and M. Ben Shalom, Sliding van der Waals polytypes, Nat. Rev. Phys. 7, 50 (2025)
work page 2025
-
[8]
Vizner Stern et al., Interfacial ferroelectricity by van der Waals sliding, Science 372, 1462 (2021)
M. Vizner Stern et al., Interfacial ferroelectricity by van der Waals sliding, Science 372, 1462 (2021)
work page 2021
Show all 46 references
-
[9]
Yasuda et al., Ultrafast high-endurance memory based on sliding ferroelectrics, Science 385, 53 (2024)
K. Yasuda et al., Ultrafast high-endurance memory based on sliding ferroelectrics, Science 385, 53 (2024)
2024
-
[10]
Assavachin and F
S. Assavachin and F. E. Osterloh, Ferroelectric Polarization in BaTiO3 Nanocrystals Controls Photoelectrochemical Water Oxidation and Photocatalytic Hydrogen Evolution, J. Am. Chem. Soc. 145, 18825 (2023)
2023
-
[11]
J. F. Scott, Applications of Modern Ferroelectrics, Science 315, 954 (2007)
2007
-
[12]
Jiang et al., Enabling ultra-low-voltage switching in BaTiO3, Nat
Y. Jiang et al., Enabling ultra-low-voltage switching in BaTiO3, Nat. Mater. 21, 779 (2022)
2022
-
[13]
D. Yang, J. Liang, J. Wu, Y. Xiao, J. I. Dadap, K. Watanabe, T. Taniguchi, and Z. Ye, Non-volatile electrical polarization switching via domain wall release in 3R -MoS2 bilayer, Nat. Commun. 15, 1389 (2024)
2024
-
[14]
T. H. Yang et al. , Ferroelectric transistors based on shear -transformation-mediated rhombohedral-stacked molybdenum disulfide, Nat. Electron. 7, 29 (2024)
2024
-
[15]
Wang et al
L. Wang et al. , Bevel -edge epitaxy of ferroelectric rhombohedral boron nitride single crystal, Nature 629, 74 (2024)
2024
-
[16]
Yasuda, X
K. Yasuda, X. Wang, K. Watanabe, T. Taniguchi, and P. Jarillo -Herrero, Stacking-engineered ferroelectricity in bilayer boron nitride, Science 372, 1458 (2021)
2021
-
[17]
Rogé e, L
L. Rogé e, L. Wang, Y. Zhang, S. Cai, P. Wang, M. Chhowalla, W. Ji, and S. P. Lau, Ferroelectricity in untwisted heterobilayers of transition metal dichalcogenides, Science 376, 973 (2022)
2022
-
[18]
Paillard, B
C. Paillard, B. Xu, B. Dkhil, G. Geneste, and L. Bellaiche, Photostriction in Ferroelectrics from Density Functional Theory, Phys. Rev. Lett. 116, 247401 (2016)
2016
-
[19]
Xiang, Y
Z. Xiang, Y. Chen, Y. Quan, and B. Liao, High -Throughput Search for Photostrictive Materials Based on a Thermodynamic Descriptor, J. Am. Chem. Soc. 146, 33732 (2024)
2024
-
[20]
Tsai et al
H. Tsai et al. , Light -induced lattice expansion leads to high -efficiency perovskite solar cells, Science 360, 67 (2018)
2018
-
[21]
Yang and S
Q. Yang and S. Meng, Light -Induced Complete Reversal of Ferroelectric Polarization in Sliding Ferroelectrics, Phys. Rev. Lett. 133, 136902 (2024)
2024
-
[22]
J. Wang, X. Li, X. Ma, L. Chen, J.-M. Liu, C.-G. Duan, J. Í ñ iguez-Gonzá lez, D. Wu, and Y. Yang, Ultrafast Switching of Sliding Polarization and Dynamical Magnetic Field in van der Waals Bilayers Induced by Light, Phys. Rev. Lett. 133, 126801 (2024)
2024
-
[23]
Gao and L
L. Gao and L. Bellaiche, Large Photoinduced Tuning of Ferroelectricity in Sliding Ferroelectrics, Phys. Rev. Lett. 133, 196801 (2024)
2024
-
[24]
Wang et al., Lattice Expansion Enables Large Surface Carrier Diffusion in WS2 Monolayer, ACS Energy Lett
L. Wang et al., Lattice Expansion Enables Large Surface Carrier Diffusion in WS2 Monolayer, ACS Energy Lett. 10, 1741 (2025)
2025
-
[25]
Runge and E
E. Runge and E. K. U. Gross, Density -Functional Theory for Time -Dependent Systems, Phys. Rev. Lett. 52, 997 (1984)
1984
-
[30]
Ordejó n, E
P. Ordejó n, E. Artacho, and J. M. Soler, Self-consistent order-N density-functional calculations for very large systems, Phys. Rev. B 53, R10441 (1996)
1996
-
[31]
Kresse and J
G. Kresse and J. Furthmü ller, Efficient iterative schemes for ab initio total -energy calculations using a plane-wave basis set, Phys. Rev. B 54, 11169 (1996)
1996
-
[34]
Schonfeld, J
B. Schonfeld, J. J. Huang, and S. C. Moss, Anisotropic mean -square displacements (MSD) in single-crystals of 2H- and 3R-MoS2, Acta Crystallogr., Sect. B 39, 404 (1983)
1983
-
[35]
Henkelman, B
G. Henkelman, B. P. Uberuaga, and H. Jó nsson, A climbing image nudged elastic band method for finding saddle points and minimum energy paths, J. Chem. Phys. 113, 9901 (2000)
2000
-
[36]
Henkelman and H
G. Henkelman and H. Jó nsson, Improved tangent estimate in the nudged elastic band method for finding minimum energy paths and saddle points, J. Chem. Phys. 113, 9978 (2000)
2000
-
[37]
J.-U. Lee, S. Woo, J. Park, H. C. Park, Y.-W. Son, and H. Cheong, Strain-shear coupling in bilayer MoS2, Nat. Commun. 8, 1370 (2017)
2017
-
[38]
R. T. Sam, T. Umakoshi, and P. Verma, Probing stacking configurations in a few layered MoS2 by low frequency Raman spectroscopy, Sci. Rep. 10, 21227 (2020)
2020
-
[39]
J. Park, I. W. Yeu, G. Han, C. S. Hwang, and J. -H. Choi, Ferroelectric switching in bilayer 3R MoS2 via interlayer shear mode driven by nonlinear phononics, Sci. Rep. 9, 14919 (2019)
2019
-
[40]
Photostriction-tunable Polarization and Structural Dynamics in Interlayer Sliding Ferroelectrics
L. Liang, J. Zhang, B. G. Sumpter, Q.-H. Tan, P.-H. Tan, and V. Meunier, Low-Frequency Shear and Layer-Breathing Modes in Raman Scattering of Two -Dimensional Materials, ACS Nano 11, 11777 (2017). Figures and Captions Figure 1. Atomic and electronic structures of bilayer 3R -M...
2017
-
[41]
Runge and E
E. Runge and E. K. U. Gross, Density -Functional Theory for Time -Dependent Systems, Phys. Rev. Lett. 52, 997 (1984). 0 1 2 3 4 5 0.0 0.3 0.6 0.9 1.2Bandgap (eV) Expansion (%) 0 1 2 3 4 5 2.94 2.97 3.00 3.03 3.06Interlayer distance (Å) Expansion (%) 0 1 2 3 4 5 2.94 2.97 3.00 ...
1984
-
[42]
C. Lian, M. Guan, S. Hu, J. Zhang, and S. Meng, Photoexcitation in Solids: First -Principles Quantum Simulations by Real-Time TDDFT, Adv. Theory Simul. 1, 1800055 (2018)
2018
-
[43]
Meng and E
S. Meng and E. Kaxiras, Real -time, local basis -set implementation of time -dependent density functional theory for excited state dynamics simulations, J. Chem. Phys. 129, 054110 (2008)
2008
-
[44]
Sá nchez-Portal, P
D. Sá nchez-Portal, P. Ordejó n, E. Artacho, and J. M. Soler, Density-functional method for very large systems with LCAO basis sets, Int. J. Quantum Chem. 65, 453 (1997)
1997
-
[45]
M. S. José , A. Emilio, D. G. Julian, G. Alberto, J. Javier, O. Pablo, and S.-P. Daniel, The SIESTA method for ab initio order-N materials simulation, J. Phys.: Condens. Matter 14, 2745 (2002)
2002
-
[46]
Kaduk, T
B. Kaduk, T. Kowalczyk, and T. Van Voorhis, Constrained Density Functional Theory, Chem. Rev. 112, 321 (2012)
2012
-
[47]
Kresse and J
G. Kresse and J. Hafner, Ab initio molecular dynamics for liquid metals , Phys. Rev. B 47, 558 (1993)
1993
-
[48]
Kresse and J
G. Kresse and J. Furthmü ller, Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set, Phys. Rev. B 54, 11169 (1996)
1996
-
[49]
J. P. Perdew, J. A. Chevary, S. H. Vosko, K. A. Jackson, M. R. Pederson, D. J. Singh, and C. Fiolhais, Atoms, molecules, solids, and surfaces: Applications of the generalized gradient approximation for exchange and correlation, Phys. Rev. B 46, 6671 (1992)
1992
-
[50]
P. E. Blö chl, Projector augmented-wave method, Phys. Rev. B 50, 17953 (1994)
1994
-
[51]
Klimeš, D
J. Klimeš, D. R. Bowler, and A. Michaelides, Phys. Rev. B 83, 195131 (2011)
2011
-
[52]
R. D. King-Smith and D. Vanderbilt, Theory of polarization of crystalline solids, Phys. Rev. B 47, 1651 (1993)
1993
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.