Pith. sign in

REVIEW 4 major objections 4 minor 53 references

Quadratic piezoelectricity from stacking-engineered interference in multilayer sliding ferroelectrics

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

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

Pith's one-line read Stacking order in multilayer sliding ferroelectrics can cancel the linear piezoelectric response and expose a purely quadratic one, offering a programmable route to nonlinear electromechanics.

desk verdict Genuinely interesting design principle with solid DFT for the two showcase stackings, but the interference mechanism and gamma_e descriptor rest on a non-unique channel partition that needs a robustness test. read the letter →

arxiv 2607.29214 v1 pith:ZSYUDK5S submitted 2026-07-31 cond-mat.mtrl-sci

Jiaxuan Fan , Yuexing Zhao , Xiao-Ping Li , Yurong Yang , Xing-Qiu Chen , Lei Wang This is my paper · ORCID
classification cond-mat.mtrl-sci PACS 77.65.-j77.80.-e
keywords slidingferroelectricspiezoelectricinterferencenonlinearpiezoelectricityquadraticresponsestackingengineeringinterlayerchannelsMoS2NiTe2
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 tries to establish that the out-of-plane piezoelectric response of multilayer sliding ferroelectrics is not a single material property but a sum of stacking-dependent contributions from each interlayer gap. When these contributions have the same sign, they reinforce to give a large linear response; when they have opposite signs, they cancel the linear coefficient while leaving a finite quadratic coefficient, producing a nearly parabolic polarization–strain curve. The authors show this mechanism in BAAC-stacked MoS2 and CBB-stacked NiTe2, and predict that sliding a tetralayer MoS2 sample can switch it between linear, quadratic, and sign-inverted response states. A sympathetic reader would care because it offers a general, reconfigurable design principle for nonlinear piezoelectricity, avoiding the delicate strain or phase-instability tuning previously required.

What carries the argument

The central object is the interface-resolved channel decomposition: the total out-of-plane piezoelectric coefficient e33 is written as the sum of per-gap coefficients e33^(i), each associated with a specific interlayer gap and determined by local stacking. This decomposition, combined with the interference factor γe that measures constructive versus destructive addition, carries the argument. The paper supports it with an energy–polarization landscape in interlayer-distance space (d1, d2), where the linear response is governed by the projection of the strain-relaxation path onto the polarization gradient; near-orthogonality between the two suppresses e33. A generalized Ginzburg–Landau free e

What would settle it

Measure the polarization–strain curve of a ferrielectric CBAB-stacked tetralayer MoS2 sample; if it is not close to parabolic (i.e., if e33 is not strongly suppressed), the predicted destructive interference is not realized. Alternatively, compute e33 and B333 of BAAC-MoS2 using a method that does not rely on per-gap partitioning (e.g., direct Berry-phase strain derivative); if e33 remains large, the interference explanation collapses.

Watch

Extended reading notes

Core claim

The central claim is that each van der Waals gap in a multilayer sliding ferroelectric acts as an independent piezoelectric channel whose sign and magnitude are set by the local stacking environment, so the total out-of-plane linear coefficient is the algebraic sum of channel coefficients. Constructive interference of same-signed channels gives a linear-dominated response, while destructive interference of opposite-signed channels suppresses the linear coefficient without extinguishing the quadratic coefficient. Quantitative first-principles results show BAAC-MoS2 with e33 ≈ +4.68 pC/m and B333 ≈ −1073.83 pC/m, and CBB-NiTe2 with e33 ≈ +1.07 pC/m and B333 ≈ −593.26 pC/m, both approaching the

Load-bearing premise

The macroscopic piezoelectric response factorizes cleanly into independent per-gap channels e33^(i), so that the total linear coefficient is an algebraic sum of these local channel coefficients—a partition of charge redistribution that is a modeling choice rather than a directly measurable quantity.

Editorial extensions

If this is right

  • Stacking order alone can suppress the linear piezoelectric coefficient without requiring critical strain or material-specific instability, making nonlinear piezoelectricity designable rather than accidental.
  • The interference factor γe serves as a compact descriptor: systems with γe < 0.2 exhibit vertex strains of only a few percent, enabling screening for near-parabolic piezoelectric candidates.
  • In a single tetralayer MoS2 sample, lateral sliding is predicted to connect linear-dominated, quadratic-dominated, and sign-inverted response states, providing a non-volatile, reconfigurable electromechanical element.
  • The quadratic-dominated state offers a route to frequency doubling: a sinusoidal vertical ac field would produce an output at twice the frequency, extending beyond ordinary ferroelectric polarization reversal.
  • The mechanism naturally generalizes to other layered systems with multiple comparable interfacial polar units, such as moiré superlattices and oxide heterostructures, though the present demonstration is in transition-metal dichalcogenides.

Reading between the lines

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

  • The same interference logic could apply to other stacking-dependent order parameters beyond polarization—for instance, nonlinear optical susceptibilities or spin–orbit-coupled responses—where per-gap sign cancellation might tune higher-order coefficients while suppressing the leading term.
  • Because the quadratic-dominated state is reached at zero strain, the vertex strain itself becomes a stack-tunable engineering parameter: one could design a material whose harmonic output peaks at zero bias, potentially useful for energy harvesting or mechanical logic where amplitude is not the distinguishing feature.
  • An experimental test could use artificially assembled heterostructures with one chemically modified gap (e.g., an intercalant): if the per-gap channel decomposition is faithful, the modified gap's contribution should add linearly to the rest, allowing direct measurement of individual channel coefficients.
  • If sliding-programmable quadratic response is confirmed, it suggests a mechanical 'frequency-doubling switch' that could be toggled by a lateral field, which would be a new functionality for van der Waals devices beyond simple polarization reversal.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

4 major / 4 minor

Summary. The manuscript proposes a stacking-engineered 'piezoelectric interference' mechanism in multilayer sliding ferroelectrics. It decomposes the out-of-plane linear piezoelectric coefficient into per-interlayer-gap channels e33^(i), argues that same-signed channels give constructive interference and opposite-signed channels give destructive interference, and uses this to design stackings with suppressed linear response but finite quadratic response. Direct DFT calculations on representative MoS2 and NiTe2 multilayer stackings show near-parabolic P_z(ε) curves for BAAC-MoS2 and CBB-NiTe2, with vertex strains reduced by factors of ~25 and ~60 relative to CBA-MoS2. The paper also presents a Ginzburg–Landau model, a 12-material materials map based on an interference factor γ_e, and predicted sliding pathways in tetralayer MoS2 that connect linear-dominated, quadratic-dominated, and sign-inverted response states.

Significance. If the interference mechanism is robust, the work offers a genuinely useful design principle for nonlinear piezoelectricity: instead of delicate critical tuning, stacking order would provide a reversible, layer-by-layer control of the balance between linear and quadratic electromechanical response. The direct DFT data are internally consistent, and the use of four vdW treatments with a consistent CI/DI classification is a clear strength; the classical-dipole polarization method is also benchmarked against Berry-phase results for semiconducting references. The sliding-pathway predictions for tetralayer MoS2 are experimentally testable. However, the central explanatory claim rests on a per-gap channel decomposition that is not uniquely defined and is not yet shown to be invariant under reasonable alternative partitions. The geometric orthogonality argument in Fig. 2 provides a more partition-independent route, but it is demonstrated for only one system. Thus the significance is conditional: the numerical findings are valuable, but the 'interference' mechanism and the γ_e descriptor need additional validation before they can be accepted as a general design principle.

major comments (4)
  1. [Results — 'Stacking-controlled piezoelectric interference'; Eq. (5), Eq. (9), Supplementary Note 3] The central classification into constructive versus destructive interference, and the descriptor γ_e, is built on the per-gap decomposition e33^(i). This decomposition is a modeling choice and is not unique: the total e33 is a single ground-state response, and assigning it to individual vdW gaps requires a partition convention. The isolated-gap variation used in Supplementary Note 3 assumes that strain-induced charge redistribution localizes to one gap and ignores collective relaxation across gaps. If an equally reasonable partition (e.g., Bader/Wannier decomposition of the uniform-strain response, or simultaneous independent gap variations) changes the sign of any e33^(i), then 'interference' is a property of the bookkeeping rather than a physical mechanism. The authors should demonstrate sign invariance of the e33^(i) under at least one alternative partition, or reframe the claims with
  2. [Results — 'Geometric origin of piezoelectric interference'; Fig. 2d,e] The geometric orthogonality argument (relaxation path nearly parallel to polarization contours in the (d1,d2) plane) is the most convincing, partition-independent evidence for destructive interference, but it is presented only for CBB-NiTe2. For BAAC-MoS2 — the other crucial near-parabolic system — no such landscape is shown, and the 12-material map in Fig. 3c relies on γ_e, not on direct orthogonality. To make the design principle general, the authors should either show the gradient–relaxation-path orthogonality for BAAC-MoS2 and CBAB-MoS2 as well, or provide a quantitative orthogonality criterion applied to all surveyed systems.
  3. [Results — 'Ginzburg–Landau framework and materials design space'; Eq. (6), Supplementary Note 6] The Ginzburg–Landau model is presented as an explanatory framework, but its parameters (A_i, β_i, g_i, λ, C) are not fitted to the DFT data, and Eq. (6) is not used to compute e^(1), e^(2) from microscopic inputs. Since Eq. (5) is just the definition of the decomposition, the model can only restate the cancellation algebra; it cannot independently validate the constructive/destructive classification. The authors should either fit the GL parameters to the DFT-computed channels and show predictive agreement, or explicitly label the model as a phenomenological illustration rather than a derivation.
  4. [Methods and Results — reported coefficients] The reported values of e33 and B333 (e.g., +24.05, +4.68, +1.07 pC/m and −218.15, −1073.83, −593.26 pC/m) are given without uncertainties, fit windows, or convergence tests. Because the central claim of 'suppressed linear response' depends on comparing these numbers to zero and to one another, the reader cannot assess whether e33 ≈ +4.68 pC/m is significantly different from zero or from the CBA value. The authors should report fitting strain ranges, standard errors from the quadratic fits, and convergence checks with respect to k-point density and strain steps.
minor comments (4)
  1. [Methods — First-principles calculations] The text contains typographical spacing in 'V ASP' and 'DA TA A V AILABILITY'; these should be corrected to 'VASP' and 'DATA AVAILABILITY'.
  2. [Throughout] The stacking notation (CBA, BAAC, CBB, CBAC, CBAB, ACAB, BCAB) is used without a general definition in the main text. A short convention statement or a table mapping each label to the layer stacking sequence would improve readability.
  3. [Methods — Polarization and piezoelectricity] For semimetallic NiTe2, the classical dipole method is used because Berry phase is ill-defined; the benchmark against Berry phase is said to be in Supplementary Note 1. A one-sentence summary of the validation for the semimetallic case in the main text would help the reader trust the NiTe2 numbers, which are central to one of the two main parabolic examples.
  4. [Fig. 3c] The materials map should explicitly label axes and define the color/symbol coding in the caption; in the provided text it is not clear what the axes are or which symbol corresponds to which vdW treatment/material. If this is a rendering issue, the figure must be self-contained.

Circularity Check

0 steps flagged · score 1.0 of 10

No load-bearing circularity; the central DFT predictions are self-contained and the interference framing is interpretive rather than derived from itself.

full rationale

The paper's quantitative claims—e33, B333, the factor-of-25 crossover reduction, and the four-state sliding response—are obtained by direct DFT fits to Pz(epsilon) and are not generated by the Ginzburg-Landau decomposition. The interface-resolved e33^(i) are defined as separately calculated gap-response channels, and Eq. (5) is a model-level projection (P=P1+P2), not a definition that forces the DFT total; the non-uniqueness of the partition is a modeling caveat, not a circular step. The interference factor gamma_e in Eq. (9) is a normalized expression of the e33^(i) signs, so the gamma_e-vs-epsilon_v map is partly a restatement of the smallness of |e33|, but this descriptor is used for organization and design rationale rather than as the source of the first-principles predictions. The polarization method is benchmarked against Berry-phase results, and the classification is checked across four vdW treatments; no load-bearing self-citation or imported uniqueness theorem appears. The central derivation is therefore self-contained, with a mild interpretive tautology in the gamma_e descriptor that does not rise to circularity.

Assumptions & free parameters 2 free parameters · 4 assumptions · 2 invented entities

The paper's design principle rests on an additive channel decomposition that is definitionally exact once the partition is chosen, plus standard DFT/GL assumptions. The GL parameters are not quantified, the NiTe2 polarization is unbenchmarked, and the gamma_e descriptor is built from the same fitted coefficients it is used to explain.

free parameters (2)
  • GL free-energy parameters A_i, beta_i, g_i, lambda, C
    Eq (3) parameters are not fitted or reported; the CI/DI classification is obtained by choosing signs of coupling/driving terms rather than from an independent parameter-free derivation.
  • Quadratic-fit strain range for P_z(epsilon)
    e33 and B333 are obtained by fitting DFT data to Eq (1); the fit window and uncertainties are not given, so vertex strains depend on an unstated choice.
assumptions (4)
  • ad hoc to paper Total out-of-plane polarization and piezoelectric response decompose exactly into interlayer-gap channels (P=sum P_i; e_tot=sum e33^(i))
    Central to the interface-resolved analysis (Figs 1d-f, Eq 5); the partition in Supplementary Note 3 is a modeling choice not uniquely defined and not experimentally isolated.
  • domain assumption Classical dipole moment computed from slab charge density gives quantitatively reliable P_z for semimetallic NiTe2
    Methods: Berry phase is ill-defined for semimetals; the benchmark is only against semiconducting bilayers, not NiTe2.
  • domain assumption Ginzburg-Landau expansion to fourth order in P and linear in epsilon is sufficient
    Eq (3) truncates at O(epsilon P_i^2,...); no justification is given that omitted terms are small in the strain range considered.
  • domain assumption DFT with PBE + vdW corrections correctly captures stacking energetics and polarization for MoS2 and NiTe2
    Standard assumption; four vdW treatments are used but there is no experimental benchmark for the predicted electromechanical coefficients.
invented entities (2)
  • Interface-resolved piezoelectric channel e33^(i)
    purpose: Decomposes total piezoelectric response into per-gap contributions to diagnose constructive vs destructive interference
    Not directly measurable; defined via a chosen partition (Supplementary Note 3), with no experimental probe that isolates a single gap's response.
  • Interference factor gamma_e
    purpose: Descriptor to classify stacking states as constructive vs destructive and to screen for near-parabolic piezoelectricity
    A derived scalar from fitted e33^(i) values; useful as a bookkeeping metric but not an independent physical quantity.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Quadratic piezoelectricity from stacking-engineered interference in multilayer sliding ferroelectrics." pith.science (2026). https://pith.science/paper/ZSYUDK5S

@misc{pith2026260729214,
  author       = {Pith},
  title        = {Pith review of: Quadratic piezoelectricity from stacking-engineered interference in multilayer sliding ferroelectrics},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/ZSYUDK5S}},
  note         = {Machine review of arXiv:2607.29214}
}
abstract

Designing nonlinear piezoelectricity requires suppressing the linear piezoelectric coefficient without extinguishing higher-order electromechanical response, yet a general and reconfigurable route remains lacking. Here we introduce stacking-engineered piezoelectric interference as such a mechanism in multilayer sliding ferroelectrics. Combining first-principles calculations with a generalized Ginzburg--Landau framework, we show that each interlayer gap acts as a local piezoelectric channel whose sign and magnitude are determined by stacking. Constructive interference between same-signed channels produces a linear-dominated response, whereas destructive interference between oppositely signed channels suppresses the linear coefficient while preserving a finite quadratic response. Representative MoS$_2$ and NiTe$_2$ multilayers approach the parabolic limit, with BAAC-stacked MoS$_2$ reducing the linear-to-quadratic crossover strain by a factor of 25 relative to CBA-stacked MoS$_2$. Experimentally accessible tetralayer MoS$_2$ sliding pathways further connect linear-dominated, quadratic-dominated and sign-inverted states. Here, we identify stacking-engineered interference as a design principle for programmable nonlinear electromechanics in layered materials.

Figures

Figures reproduced from arXiv: 2607.29214 by the authors.

Figure 1
Figure 1. presents a comparative first-principles analysis of three representative multilayer systems. We first con￾sider CBA-stacked trilayer MoS2 (CBA-MoS2), whose out￾of-plane polarization varies almost linearly with vertical strain (Fig. 1a), giving a large positive linear coefficient of e (tot) 33 ≈ [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2 [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. b, λ > 0 favors antiparallel alignment of the two polar subunits. In an approximately antipolar ferrielectric￾like configuration, P1,0P2,0 < 0. When the effective strain￾driving terms are opposite-signed, g eff 1 g eff 2 < 0, and their stiffness-weighted contributions nearly cancel, the local linear responses become opposite-signed and strongly compensate: e (1)e (2) < 0, |e (tot)| ≪ |e (1)| + |e (2)|. (8) This dest… view at source ↗
Figures from the paper (1 more)
Figure 4
Figure 4. Figure 4: a shows the minimum-energy sliding pathway con￾necting four stacking states, denoted S1–S4, corresponding to CBAC, CBAB, ACAB, and BCAB stackings. These states are close in energy and separated by moderate calculated slid￾ing barriers, comparable to energy scales repor…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

53 extracted references

  1. [1]

    Wei, X.-K., Domingo, N., Sun, Y ., Balke, N., Dunin- Borkowski, R. E. & Mayer, J. Progress on emerging ferroelec- tric materials for energy harvesting, storage and conversion. Adv. Energy Mater.12, 2201199 (2022)

    2022

  2. [2]

    & Lee, P

    Li, T. & Lee, P. S. Piezoelectric energy harvesting technol- ogy: from materials, structures, to applications.Small Struct. 3, 2100128 (2022)

    2022

  3. [3]

    & Koc ¸, M

    Sezer, N. & Koc ¸, M. A comprehensive review on the state- of-the-art of piezoelectric energy harvesting.Nano Energy80, 105567 (2021)

    2021

  4. [4]

    & Lin, L

    Liu, H., Zhong, J., Lee, C., Lee, S.-W. & Lin, L. A compre- hensive review on piezoelectric energy harvesting technology: materials, mechanisms, and applications.Appl. Phys. Rev.5, 041306 (2018)

    2018

  5. [5]

    V ., Karthikeyan, P

    Sabarianand, D. V ., Karthikeyan, P. & Muthuramalingam, T. A review on control strategies for compensation of hysteresis and creep on piezoelectric actuators based micro systems.Mech. Syst. Signal Process.140, 106634 (2020)

    2020

  6. [6]

    Lines, M. E. & Glass, A. M.Principles and Applications of Fer- roelectrics and Related Materials(Oxford Univ. Press, Oxford, 2001)

    2001

  7. [7]

    Ferroelectric, dielectric and piezoelectric prop- erties of ferroelectric thin films and ceramics.Rep

    Damjanovic, D. Ferroelectric, dielectric and piezoelectric prop- erties of ferroelectric thin films and ceramics.Rep. Prog. Phys. 61, 1267–1324 (1998)

    1998

  8. [8]

    & Vanderbilt, D

    Bernardini, F., Fiorentini, V . & Vanderbilt, D. Spontaneous po- larization and piezoelectric constants of III-V nitrides.Phys. Rev. B56, R10024–R10027 (1997)

    1997

Show all 53 references
  1. [9]

    &´I˜niguez, J

    Dutta, S., Buragohain, P., Glinsek, S., Richter, C., Aramberri, H., Lu, H., Schroeder, U., Defay, E., Gruverman, A. &´I˜niguez, J. Piezoelectricity in hafnia.Nat. Commun.12, 7301 (2021)

    2021

  2. [10]

    & Dong, S

    Wang, Z. & Dong, S. Large in-plane negative piezoelectricity and giant nonlinear optical susceptibility in elementary ferro- electric monolayers.Phys. Rev. B108, 235423 (2023)

    2023

  3. [11]

    & Liu, J

    Hong, Y ., Deng, J., Kong, Q., Yin, Y ., Ding, X., Sun, J. & Liu, J. Z. Origin of different piezoelectric responses in elemental Sb and Bi monolayers.Phys. Rev. B109, 035204 (2024)

    2024

  4. [12]

    & Zunger, A

    Bester, G., Wu, X., Vanderbilt, D. & Zunger, A. Importance of second-order piezoelectric effects in zinc-blende semiconduc- tors.Phys. Rev. Lett.96, 187602 (2006)

    2006

  5. [13]

    & Bester, G

    Prodhomme, P.-Y ., Beya-Wakata, A. & Bester, G. Nonlinear piezoelectricity in wurtzite semiconductors.Phys. Rev. B88, 121304 (2013)

    2013

  6. [14]

    A., Schulz, S

    Caro, M. A., Schulz, S. & O’Reilly, E. P. Origin of nonlin- ear piezoelectricity in III-V semiconductors: internal strain and bond ionicity from hybrid-functional density functional theory. Phys. Rev. B91, 075203 (2015)

    2015

  7. [15]

    N., Landahl, E

    Grigoriev, A., Sichel, R., Lee, H. N., Landahl, E. C., Adams, B., Dufresne, E. M. & Evans, P. G. Nonlinear piezoelectricity in epitaxial ferroelectrics at high electric fields.Phys. Rev. Lett. 100, 027604 (2008)

    2008

  8. [16]

    & Yang, Y

    Cheng, H., Jiao, P., Wang, J., Qing, M., Deng, Y ., Liu, J.-M., Bellaiche, L., Wu, D. & Yang, Y . Tunable and parabolic piezo- electricity in hafnia under epitaxial strain.Nat. Commun.15, 394 (2024)

    2024

  9. [17]

    & Resta, R

    de Gironcoli, S., Baroni, S. & Resta, R. Piezoelectric properties of III-V semiconductors from first-principles linear-response theory.Phys. Rev. Lett.62, 2853–2856 (1989)

    1989

  10. [18]

    First-principles responses of solids to atomic dis- placements and homogeneous electric fields: implementation of a conjugate-gradient algorithm.Phys

    Gonze, X. First-principles responses of solids to atomic dis- placements and homogeneous electric fields: implementation of a conjugate-gradient algorithm.Phys. Rev. B55, 10337– 10354 (1997)

    1997

  11. [19]

    B., Kjær, K

    Katsouras, I., Asadi, K., Li, M., van Driel, T. B., Kjær, K. S., Zhao, D., Lenz, T., Gu, Y ., Blom, P. W. M., Damjanovic, D., Nielsen, M. M. & de Leeuw, D. M. The negative piezoelec- tric effect of the ferroelectric polymer poly(vinylidene fluoride). Nat. Mater.15, 78–84 (2016)

    2016

  12. [20]

    & Cohen, R

    Liu, S. & Cohen, R. E. Origin of negative longitudinal piezo- electric effect.Phys. Rev. Lett.119, 207601 (2017)

    2017

  13. [21]

    & Zhu, H

    Jing, H., Guan, C. & Zhu, H. Understanding and tuning nega- tive longitudinal piezoelectricity in hafnia.npj Comput. Mater. 10, 160 (2024)

    2024

  14. [22]

    Li, L. & Wu, M. Binary compound bilayer and multilayer with vertical polarizations: two-dimensional ferroelectrics, multifer- roics, and nanogenerators.ACS Nano11, 6382–6388 (2017)

    2017

  15. [23]

    Yang, Q., Wu, M. & Li, J. Origin of two-dimensional vertical ferroelectricity in WTe2 bilayer and multilayer.J. Phys. Chem. Lett.9, 7160–7164 (2018)

    2018

  16. [24]

    & Jarillo-Herrero, P

    Yasuda, K., Wang, X., Watanabe, K., Taniguchi, T. & Jarillo-Herrero, P. Stacking-engineered ferroelectricity in bi- layer boron nitride.Science372, 1458–1462 (2021)

    2021

  17. [25]

    & Ben Shalom, M

    Vizner Stern, M., Waschitz, Y ., Cao, W., Nevo, I., Watanabe, K., Taniguchi, T., Sela, E., Urbakh, M., Hod, O. & Ben Shalom, M. Interfacial ferroelectricity by van der Waals sliding.Science 372, 1462–1466 (2021)

    2021

  18. [26]

    Weston, A., Castanon, E. G., Enaldiev, V ., Ferreira, F., Bhat- tacharjee, S., Xu, S., Corte-Le ´on, H., Wu, Z., Clark, N., Sum- merfield, A., Hashimoto, T., Gao, Y ., Wang, W., Hamer, M., Read, H., Fumagalli, L., Kretinin, A. V ., Haigh, S. J., Kaza- kova, O., Geim, A. K., Fa...

    2022

  19. [27]

    Kim, J., Rabe, K. M. & Vanderbilt, D. Negative piezoelectric response of van der Waals layered bismuth tellurohalides.Phys. Rev. B100, 104115 (2019)

    2019

  20. [28]

    A., Liu, F., Chang, L., Ichinose, D., Funakubo, H., Hu, W., Wu, T., Liu, Z., Dong, S

    You, L., Zhang, Y ., Zhou, S., Chaturvedi, A., Morris, S. A., Liu, F., Chang, L., Ichinose, D., Funakubo, H., Hu, W., Wu, T., Liu, Z., Dong, S. & Wang, J. Origin of giant negative piezo- electricity in a layered van der Waals ferroelectric.Sci. Adv.5, eaav3780 (2019)

    2019

  21. [29]

    & Rappe, A

    Qi, Y . & Rappe, A. M. Widespread negative longitudinal piezo- electric responses in ferroelectric crystals with layered struc- tures.Phys. Rev. Lett.126, 217601 (2021)

    2021

  22. [30]

    & Liu, F

    Meng, P., Wu, Y ., Bian, R., Pan, E., Dong, B., Zhao, X., Chen, J., Wu, L., Sun, Y ., Fu, Q., Liu, Q., Shi, D., Zhang, Q., Zhang, Y .-W., Liu, Z. & Liu, F. Sliding induced multiple polariza- tion states in two-dimensional ferroelectrics.Nat. Commun.13, 7696 (2022)

    2022

  23. [31]

    Ouyang, T., Cha, S., Sun, Y ., Taniguchi, T., Watanabe, K., Ga- bor, N. M. & Lui, C. H. Electrically switching ferroelectric or- der in 3R-MoS2 layers.Nano Lett.25, 1459–1465 (2025)

    2025

  24. [32]

    C., Erturk, A., Mann, B

    Stanton, S. C., Erturk, A., Mann, B. P. & Inman, D. J. Nonlin- ear piezoelectricity in electroelastic energy harvesters: model- ing and experimental identification.J. Appl. Phys.108, 074903 (2010)

    2010

  25. [33]

    & Takagi, K

    Yamamoto, K., Kokubo, A., Sakai, K. & Takagi, K. Nonlinear piezoelectricity in PZT ceramics for generating ultrasonic phase conjugate waves.Ultrasonics38, 830–833 (2000)

    2000

  26. [34]

    Yang, L., Ding, S., Gao, J. & Wu, M. Atypical sliding and moir´e ferroelectricity in pure multilayer graphene.Phys. Rev. Lett.131, 096801 (2023)

    2023

  27. [35]

    Zhang, L., Ding, J., Xiang, H., Liu, N., Zhou, W., Wu, L., Xin, N., Watanabe, K., Taniguchi, T. & Xu, S. Electronic ferroelec- tricity in monolayer graphene moir ´e superlattices.Nat. Com- mun.15, 10905 (2024)

    2024

  28. [36]

    A., Ashoori, R., Gedik, N., Fu, L., Xu, S.-Y

    Zheng, Z., Ma, Q., Bi, Z., de la Barrera, S., Liu, M.-H., Mao, N., Zhang, Y ., Kiper, N., Watanabe, K., Taniguchi, T., Kong, J., Tisdale, W. A., Ashoori, R., Gedik, N., Fu, L., Xu, S.-Y . & Jarillo-Herrero, P. Unconventional ferroelectricity in moir ´e heterostructures.Nature5...

    2020

  29. [37]

    B., G ¨ohler, C., Liu, S., Gruver- man, A

    Hu, Y ., Lu, H., Masood, S. B., G ¨ohler, C., Liu, S., Gruver- man, A. & Alexe, M. A 2D hybrid perovskite ferroelectric with switchable polarization and photoelectric robustness down to monolayer.Nat. Commun.16, 3028 (2025)

    2025

  30. [38]

    M., Santiso, J., Cairney, J

    Liu, Y ., Niu, R., Uriach, R., Pesquera, D., Caicedo Roque, J. M., Santiso, J., Cairney, J. M., Liao, X., Arbiol, J. & Catalan, G. Coexistence of ferroelectric and ferrielectric phases in ultrathin antiferroelectric PbZrO3 thin films.Microstructures4, 2024045 (2024)

    2024

  31. [39]

    S., Manj ´on-Sanz, A., Liu, J., Ullah, S., Kong, J., Babori, C., Br ´ezard-Oudot, A., Daniel, L., Misture, S

    Pramanik, T., Htet, C. S., Manj ´on-Sanz, A., Liu, J., Ullah, S., Kong, J., Babori, C., Br ´ezard-Oudot, A., Daniel, L., Misture, S. T. & Pramanick, A. An atomistic structural description of the ferrielectric polar phase involving non-coplanar cation dis- placements.Scr. Mater...

    2025

  32. [40]

    & Joubert, D

    Kresse, G. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method.Phys. Rev. B59, 1758–1775 (1999)

    1999

  33. [41]

    & Furthm ¨uller, J

    Kresse, G. & Furthm ¨uller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B54, 11169–11186 (1996)

    1996

  34. [42]

    & Furthm ¨uller, J

    Kresse, G. & Furthm ¨uller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set.Comput. Mater. Sci.6, 15–50 (1996)

    1996

  35. [43]

    P., Burke, K

    Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradi- ent approximation made simple.Phys. Rev. Lett.77, 3865–3868 (1996)

    1996

  36. [44]

    Semiempirical GGA-type density functional con- structed with a long-range dispersion correction.J

    Grimme, S. Semiempirical GGA-type density functional con- structed with a long-range dispersion correction.J. Comput. Chem.27, 1787–1799 (2006)

    2006

  37. [45]

    & Krieg, H

    Grimme, S., Antony, J., Ehrlich, S. & Krieg, H. A consistent and accurate ab initio parametrization of density functional dis- persion correction (DFT-D) for the 94 elements H-Pu.J. Chem. Phys.132, 154104 (2010)

    2010

  38. [46]

    & Goerigk, L

    Grimme, S., Ehrlich, S. & Goerigk, L. Effect of the damping function in dispersion corrected density functional theory.J. Comput. Chem.32, 1456–1465 (2011)

    2011

  39. [47]

    Klime ˇs, J., Bowler, D. R. & Michaelides, A. Van der Waals density functionals applied to solids.Phys. Rev. B83, 195131 (2011)

    2011

  40. [48]

    & Perdew, J

    Sun, J., Ruzsinszky, A. & Perdew, J. P. Strongly constrained and appropriately normed semilocal density functional.Phys. Rev. Lett.115, 036402 (2015)

    2015

  41. [49]

    & de Gironcoli, S

    Sabatini, R., Gorni, T. & de Gironcoli, S. Nonlocal van der Waals density functional made simple and efficient.Phys. Rev. B87, 041108 (2013)

    2013

  42. [50]

    Henkelman, G., Uberuaga, B. P. & J ´onsson, H. A climbing im- age nudged elastic band method for finding saddle points and minimum energy paths.J. Chem. Phys.113, 9901–9904 (2000)

    2000

  43. [51]

    & Zhu, W

    Ding, W., Zhu, J., Wang, Z., Gao, Y ., Xiao, D., Gu, Y ., Zhang, Z. & Zhu, W. Prediction of intrinsic two-dimensional ferro- electrics in In 2Se3 and other III 2-VI3 van der Waals materials. Nat. Commun.8, 14956 (2017)

    2017

  44. [52]

    & Wang, Z

    Sheng, H., Fang, Z. & Wang, Z. Ferroelectric metals in 1T/1T ′ - phase transition metal dichalcogenide MTe 2 bilayers (M=Pt, Pd, and Ni).Phys. Rev. B108, 104109 (2023)

    2023

  45. [53]

    Grassland Talents

    Shi, J., Huan, Y ., Xiao, M., Hong, M., Zhao, X., Gao, Y ., Cui, F., Yang, P., Pennycook, S. J., Zhao, J. & Zhang, Y . Two- dimensional metallic NiTe 2 with ultrahigh environmental sta- bility, conductivity, and electrocatalytic activity.ACS Nano14, 9011–9020 (2020). ACKNOWLED...

    2020

Pith tools

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