Pith. sign in

REVIEW 3 major objections 2 minor 27 references

Effect of the local field and dipole-dipole interaction on the spontaneous ordering of dipole moments in bismuth monolayers with an orthorhombic structure

T0 review · 3 major / 2 minor · reviewed 2026-08-02 · deepseek-v4-flash

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

Pith's one-line read In orthorhombic bismuth monolayers, dipole–dipole interaction—mediated by the only nonzero Born charge component Z_xz—selectively softens the ferroelectric ZO phonon, leaves the antiferroelectric chiral mode untouched, and thereby makes fer

desk verdict Plausible and non-fitted mechanism for the ZO soft mode in bismuth monolayers, worth refereeing; the force-constant stabilization needs a systematic sensitivity analysis before the FE-over-AFE conclusion is bulletproof. read the letter →

arxiv 2607.12542 v2 pith:KUSQHMLP submitted 2026-07-14 cond-mat.mes-hall

D. V. Ponkratova , I. V. Zagorodnev , V. V. Enaldiev This is my paper
classification cond-mat.mes-hall
keywords ferroelectricityantiferroelectricitybismuthmonolayertwo-dimensionalferroelectricdipole-dipoleinteractionlocalfieldsoftopticalphononBorneffectivecharge
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

The paper argues that electrostatics, not just short-range lattice forces, determines how electric dipoles order in a monoatomic bismuth monolayer with an orthorhombic structure. It shows that the local field produced by atomic dipoles (the Lorentz–Lorenz field) admits several ferroelectric and antiferroelectric phases with boundaries in an experimentally plausible range of atomic polarizability. In the phonon picture, the same dipole–dipole interaction acts through the only nonzero Born charge component Z_xz and selectively lowers the frequency of the ZO optical mode that drives ferroelectricity, while the chiral mode that drives antiferroelectric order is unaffected. The conclusion is that ferroelectric order is more favorable than antiferroelectric order because of this dipole-mediated softening. If right, it explains the observed in-plane ferroelectric polarization of these monolayers and identifies the microscopic mechanism behind it.

What carries the argument

The machinery is the interplay between a lattice-dynamics model and a dipole-lattice electrostatics model. The load-bearing object is the Born effective charge tensor of the sublattices, whose only nonzero component Z_xz couples out-of-plane atomic displacements to in-plane polarization; it is this component that turns the ZO optical mode into an oscillating array of in-plane dipoles. The dipole field of that array, computed with Ewald lattice sums in the Lorentz–Lorenz local-field picture, lowers the ZO frequency. The competing chiral C mode at the Brillouin-zone corner involves opposite dipole moments in neighboring cells, so its macroscopic field cancels and its frequency is unchanged. Co

What would settle it

Compute the phonon spectrum of the nonpolar orthorhombic bismuth monolayer using fully relaxed lattice constants and unadjusted force constants; if the ZO mode is not the branch whose frequency is selectively lowered once dipole–dipole interactions are added, or if the C mode also softens, the central claim fails. Experimentally, a Raman or infrared measurement of the ZO phonon frequency under varied dielectric screening could test the predicted dipole-mediated softening.

Watch

Extended reading notes

Core claim

The central claim is that the ferroelectric instability of an orthorhombic bismuth monolayer is driven by dipole–dipole interactions. Working from atomic polarizabilities and the Lorentz–Lorenz local field, the paper finds spontaneous solutions for the sublattice dipole moments in zero external field, giving several ferroelectric and antiferroelectric phases whose boundaries lie at physically relevant polarizabilities (roughly 6.5–7.7 Å^3). Working from the phonon spectrum of the nonpolar lattice, it identifies two soft optical modes: the ZO mode, whose sublattice displacements add up to a uniform out-of-plane deformation with a net in-plane polarization, and a chiral C mode at the Brillouin

Load-bearing premise

The conclusion rests on the phonon model being reliable enough that the relative softness of the ferroelectric and antiferroelectric modes is real; that model uses DFT parameters from a lattice that is not at equilibrium, stabilized only by hand-adjusting two force constants, and the in-plane polarizability is inferred from a macroscopic potential slope.

Editorial extensions

If this is right

  • If the central claim is correct, the ferroelectric phase observed in orthorhombic bismuth monolayers is stabilized by long-range dipole–dipole forces, not only by short-range force-constant softness.
  • The same electrostatic mechanism predicts a cascade of antiferroelectric and ferroelectric phases as polarizability or lattice parameters vary, so tuning strain or lattice constant could switch between orders.
  • Because the ZO phonon frequency is lowered by the dipolar field, the soft mode should be detectable and sensitive to the dielectric environment.
  • The mechanism ties the direction of ferroelectric polarization (in-plane, along the armchair axis) directly to the nonzero Z_xz Born charge.
  • Other monoatomic monolayers with similar symmetry and charge tensors may show the same selective softening and should be examined for analogous ferroelectric order.

Reading between the lines

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

  • Beyond the paper: because the softening is set by the dipole lattice sums, embedding the monolayer in a higher-permittivity environment should suppress the ZO-mode instability and could switch the ground state toward antiferroelectric order; this is a tunable experimental knob the paper does not discuss.
  • Beyond the paper: the same Z_xz-mediated coupling should make the ZO phonon strongly infrared- or Raman-active only for the ferroelectric-distortion symmetry; comparative spectroscopy above and below the transition could isolate this mechanism from generic lattice anharmonicity.
  • Beyond the paper: one could close the loop between the two approaches by deriving the atomic polarizability tensor from the same force-constant and Born-charge data used in the phonon calculation, giving a parameter-free check of both phase diagrams.
  • Beyond the paper: the claim implies a temperature-dependent force-constant renormalization is needed for consistency with the experimentally stable high-temperature phase; a finite-temperature or anharmonic calculation of the two mode frequencies would test whether the selective softening survives.
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

3 major / 2 minor

Summary. The manuscript studies spontaneous dipole ordering in orthorhombic bismuth monolayers using two approaches: (i) a Lorentz–Lorenz local-field model with atomic polarizabilities from DFT, which yields phase boundaries among paraelectric, antiferroelectric, and ferroelectric orders; and (ii) a nearest-neighbour phonon model with Born effective charges, which identifies a ZO soft mode at the zone center and a chiral (C) soft mode at the M point. The central claim is that the dipole–dipole interaction, mediated by the only nonzero Born charge component Z_xz, selectively reduces the ZO phonon frequency and therefore makes ferroelectric order more favorable than the antiferroelectric C-mode order.

Significance. If the central claim is correct, the paper provides a mechanism-based explanation for the recently observed in-plane ferroelectricity in monatomic bismuth monolayers, connecting local-field electrostatics to phonon instabilities. The paper is careful to report DFT-derived parameters, Ewald lattice sums, and supplementary spectra, and it makes specific, falsifiable predictions about the competition between FE and AFE phases. However, the phonon conclusion rests on a force-constant model that requires ad hoc stabilization, and the displayed dynamical equations contain a potentially load-bearing q→0 approximation for the dipole–dipole interaction. These issues must be resolved before the central claim can be accepted.

major comments (3)
  1. [Section II, Eqs. (10)–(11)] The dipole field from a phonon with wave vector q must be computed with the q-dependent lattice sum Π^{ss'}_{jk}(q), not with the q→0 limit. As written, Eq. (11) uses Π^{ss'}_{jk}(q→0) for every q, so the comparison at the M point (q_M=(π/a,π/b)) is not described correctly. The staggered-dipole cancellation invoked for the C mode is a finite-q effect; it cannot be obtained from a q→0 lattice sum. If the numerical implementation actually used Π(q), the equations should be corrected and the distinction stated. If Π(0) was used, the conclusion that the dipole–dipole interaction leaves the C mode unchanged is not established.
  2. [Section II, Fig. 3 and Table S1] The conclusion that the dipole interaction selectively softens the Γ-point ZO mode relative to the M-point C mode depends sensitively on the force constants. The DFT values in Table I produce an unstable nonpolar lattice even without dipole–dipole interaction (imaginary branches in Fig. 3a). The only stabilized set shown, Table S1, requires flipping the sign of Φ12_zz and reducing Φ12_xy by 5%, precisely the parameters that control the ZO and C modes. The statement that 'the choice of the original or adjusted force constants does not change the main result' is not demonstrated: no comparative soft-mode analysis for the stabilized set is presented in the main text or supplement. Please provide a sensitivity analysis over plausible stabilizing corrections and show the ZO and C mode frequencies with and without dipole–dipole interaction for each set.
  3. [Section S1 and Eq. (4)] The polarizabilities used in the Lorentz–Lorenz condition (4) are extracted in Section S1 using, for in-plane fields, a 'local field' defined by the slope of the macroscopic potential difference in nanoribbons (Fig. S1). This is the macroscopic average field, not the local field at the atomic sites. For a polarizable lattice these differ by the Lorentz field. Using the macroscopic-field slope in Eq. (S1) yields an effective polarizability, and inserting it into the dipole lattice-sum condition (4) is inconsistent with the microscopic polarizability required there. The resulting phase boundaries in Fig. 2(b) and the physically relevant range of α should be re-examined, or the method should be justified explicitly.
minor comments (2)
  1. [Fig. 3 inset/caption] The notation 'qM = 0' for the ZO mode is confusing; the zone-center point is Γ, while qM is used for the M point (π/a, π/b). Please use standard labels or define the notation.
  2. [Section II and Eq. (S9)] The Born charge signs in the main text, Z(1)_xz = −Z(2)_xz = −Z(3)_xz = Z(4)_xz ≈ 1.5e, differ from Eq. (S9) by an overall sign (S9 has Z(1)=Z(4)=−1.49e and Z(2)=Z(3)=+1.49e). The phonon frequencies are invariant under a global sign, but the text and supplementary material should be made consistent.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the central derivations are self-contained computations from externally computed DFT inputs and Ewald lattice sums.

full rationale

The paper's two derivations each take external DFT-computed inputs—atomic polarizability tensors (Supplementary S1) and Born effective charges plus nearest-neighbor force constants (Supplementary S1/S3)—and combine them with independently derived lattice sums (Ewald summation, Section S2) and a symmetry-parameterized dynamical matrix. The local-field phase boundaries in Fig. 2 are obtained by setting det[α^{-1} − Π] = 0 from Eq. (4), a genuine self-consistency condition; α is a single-atom response property, not a parameter fitted to the FE/AFE phases. The phonon conclusion follows from Eq. (11), in which the dipole term −ZΠZ is evaluated with the computed Born charges and Ewald sums; the ZO-mode softening at q=0 is a calculated consequence, while the C-mode cancellation at the M point is a consequence of the dipole-field geometry. No result is defined in terms of the quantity it purports to predict, and no self-citation supplies a load-bearing premise. The adjusted force constants in Table S1 and the assertion that the main result is unchanged are legitimate robustness/correctness concerns, not circularity: the stabilization corrections are not fitted to the relative ZO-versus-C softening under the dipole interaction, and the original DFT force constants already contain the soft modes. The DFT inputs are externally reproducible from stated first-principles calculations, so a reader could falsify the model without invoking the paper's conclusions. Therefore no circular step is present.

Assumptions & free parameters 4 free parameters · 7 assumptions · 0 invented entities

The central claim rests on DFT-derived inputs (polarizability tensors, Born charges, force constants) plus a classical local-field/phonon model. No new particles or fields are introduced. The main free choices are the scanned polarizability, the ad hoc force-constant corrections, and the environment permittivity.

free parameters (4)
  • Isotropic diagonal polarizability α = varied; phase boundaries at α ≈ 6.5–7.7 ų in Fig. 2
    The determinant condition is scanned over α; the physical relevance relies on this range matching literature/DFT values.
  • Force constant Φ12_zz = sign flipped: +0.409 → −0.409 eV/Ų
    Introduced in Table S1 to stabilize the nonpolar lattice at experimental parameters.
  • Force constant Φ12_xy = 4.54 → 4.313 eV/Ų (5% decrease)
    Adjusted to make the chiral C mode real in the absence of dipole interaction (Table S1).
  • ZO displacement amplitude h0 = 0.94 Å
    Chosen from DFT energy minimum in Section S5 to evaluate ferroelectric-phase stability; not needed for the main soft-mode comparison.
assumptions (7)
  • domain assumption Each atom responds as a point dipole with a polarizability tensor; higher multipoles are neglected.
    Used in Eq. (3) to relate induced dipole to local field.
  • standard math The local field of the periodic dipole array is given by the Ewald-summed lattice sum Π (Eq. 2).
    Standard electrostatics for periodic dipoles.
  • domain assumption The Born effective charge tensor Z couples displacements to macroscopic polarization and to the dipole field (Eqs. 9-10).
    Standard phonon/polarization model.
  • domain assumption Only nearest-neighbor force constants matter; the dynamical matrix is parametrized by Φ11, Φ12, Φ14.
    Eqs. (12)-(18); validated against Phonopy in Fig. S2, but long-range force constants neglected.
  • domain assumption DFT with LDA/GGA-PBE and no spin-orbit coupling gives accurate polarizabilities, Born charges, and force constants for Bi monolayers.
    Section S1/S3; SOC is known to be significant in bismuth.
  • domain assumption The static permittivity of the environment is ε=1.
    Used for Fig. 3(b); screening by substrate or encapsulation would weaken dipole interaction.
  • domain assumption The experimental lattice parameters define the nonpolar phase even though they are not equilibrium positions in DFT, with stability restored by temperature or ad hoc force-constant corrections.
    Section II and S4; underlies the reality of the soft modes.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Effect of the local field and dipole-dipole interaction on the spontaneous ordering of dipole moments in bismuth monolayers with an orthorhombic structure." pith.science (2026). https://pith.science/paper/KUSQHMLP

@misc{pith2026260712542,
  author       = {Pith},
  title        = {Pith review of: Effect of the local field and dipole-dipole interaction on the spontaneous ordering of dipole moments in bismuth monolayers with an orthorhombic structure},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/KUSQHMLP}},
  note         = {Machine review of arXiv:2607.12542}
}
read the original abstract

Using two complementary approaches, the instability of bismuth monolayers with an orthorhombic structure with respect to the emergence of spontaneous electric dipole moments at the crystal lattice sites has been studied. The first approach, based on determining the dipole moments of lattice sites through the polarizability of bismuth atoms and taking into account the local Lorentz-Lorenz field, suggests the possibility of the existence of several antiferroelectric and ferroelectric phases in orthorhombic bismuth monolayers. Using the second approach, the vibrational spectrum of a nonpolar symmetric lattice has been analyzed, and two types of soft optical modes corresponding to the ferroelectric and antiferroelectric instabilities of the monolayers have been revealed. The relationship of the Born effective charges to the polarization direction in the ferroelectric phase, as well as their role in the reduction of the frequency of the polar optical phonon, which characterizes the lattice deformation in the ferroelectric phase due to the dipole-dipole interaction, has been shown.

Figures

Figures reproduced from arXiv: 2607.12542 by the authors.

Figure 1
Figure 1. FIG. 1. Lattice structure of an orthorhombic bismuth mono [PITH_FULL_IMAGE:figures/full_fig_p001_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. (a, b) Parametric dependences of [PITH_FULL_IMAGE:figures/full_fig_p002_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3. (a) Phonon spectrum in the orthorhombic bismuth [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗

Review history (2 revisions) →

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

27 extracted references · 1 linked inside Pith

  1. [1]

    Valasek, Phys

    J. Valasek, Phys. Rev.17, 475 (1921)

    1921

  2. [2]

    Jaffe, R

    B. Jaffe, R. S. Roth, and S. Marzullo, J. Appl. Phys.25, 809 (1954)

    1954

  3. [3]

    Acosta, N

    M. Acosta, N. Novak, V. Rojas, S. Patel, R. Vaish, J. Ko- ruza, G. A. Rossetti, and J. Rödel, Appl. Phys. Rev.4, 041305 (2017)

    2017

  4. [4]

    N. Wang, Z. Shen, W. Luo, H.-K. Li, Z.-J. Xu, C. Shi, H.-Y. Ye, S. Dong, and L.-P. Miao, Nat. Commun.15, 10160 (2024)

    2024

  5. [5]

    J. Gou, H. Bai, X. Zhang, Y. L. Huang, S. Duan, A. Ar- iando, S. A. Yang, L. Chen, Y. Lu, and A. T. S. Wee, Nature (London)617, 67 (2023)

    2023

  6. [6]

    S. A. Mikhailov, Phys. Rev. B88, 195410 (2013)

    2013

  7. [7]

    armchair

    or Group IV monochalcogenides (such as lead sul- fide) [8] is associated with a lattice deformation dictated by the softness of the corresponding out-of-plane opti- cal mode of lattice vibrations (ZO) [9] (see Fig. 3), the physical reason for its softening remains unclear. At the same time, it is well known that the frequency of po- lar phonons correspond...

    arXiv 2026

  8. [8]

    C. Xiao, F. Wang, S. A. Yang, Y. Lu, Y. Feng, and S. Zhang, Adv. Funct. Mater.28, 1707383 (2018)

    2018

Show all 27 references
  1. [9]

    Sutter, H

    P. Sutter, H. P. Komsa, H. Lu, A. Gruverman, and E. Sutter, Nano Today37, 101082 (2021)

    2021

  2. [10]

    Aktürk, O

    E. Aktürk, O. Ü. Aktürk, and S. Ciraci, Phys. Rev. B 94, 014115 (2016)

    2016

  3. [11]

    V. N. Murzin, R. E. Pasynkov, and S. P. Solov’ev, Sov. Phys. Usp.10, 453 (1968)

    1968

  4. [12]

    P. P. Ewald, Ann. Phys.369, 253 (1921)

    1921

  5. [13]

    Kittel, Phys

    C. Kittel, Phys. Rev.82, 729 (1951)

    1951

  6. [14]

    J. Gou, L. Kong, X. He, Y. L. Huang, J. Sun, S. Meng, K. Wu, L. Chen, and A. T. S. Wee, Sci. Adv.6, eaba2773 (2020)

    2020

  7. [15]

    Maroulis,Atoms, Molecules and Clusters in Elec- tric Fields: Theoretical Approaches to the Calculation of Electric Polarizability, Vol

    G. Maroulis,Atoms, Molecules and Clusters in Elec- tric Fields: Theoretical Approaches to the Calculation of Electric Polarizability, Vol. 1 (World Scientific/Imperial College Press, London, 2006)

    2006

  8. [16]

    Effect of the local field and dipole–dipole interaction on the spontaneous ordering of dipole moments in bismuth monolayers with an orthorhombic structure

    P. Schwerdtfeger and J. K. Nagle, Mol. Phys.117, 1200 (2019). 6 Supplementary Material to the article “Effect of the local field and dipole–dipole interaction on the spontaneous ordering of dipole moments in bismuth monolayers with an orthorhombic structure” S1. ESTIMA TE OF B...

    2019

  9. [17]

    Kresse, J

    G. Kresse, J. Hafner, Phys. Rev. B47, 558 (1993)

    1993

  10. [18]

    Kresse and J

    G. Kresse and J. Furthmüller, Phys. Rev. B54, 11169 (1996)

    1996

  11. [19]

    D. M. Ceperley and B. J. Alder, Phys. Rev. Lett. 45, 566 (1980)

    1980

  12. [20]

    J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett.77, 3865 (1996)

    1996

  13. [21]

    P. E. Blöchl, Phys. Rev. B50, 17953 (1994)

    1994

  14. [22]

    Neugebauer and M

    J. Neugebauer and M. Scheffler, Phys. Rev. B46, 16067 (1992)

    1992

  15. [23]

    J. Gou, L. Kong, X. He, Y. L. Huang, J. Sun, S. Meng, K. Wu, L. Chen, A. T. S. Wee, Sci. Adv.6, eaba2773 (2020)

    2020

  16. [24]

    A. Togo, I. Tanaka, Scr. Mater.108, 1 (2015)

    2015

  17. [25]

    A. Togo, J. Phys. Soc. Jpn.92, 012001 (2023)

    2023

  18. [26]

    H. J. Monkhorst, J. D. Pack, Phys. Rev. B13, 5188 (1976)

    1976

  19. [27]

    Aktürk, O

    E. Aktürk, O. Ü. Aktürk, S. Ciraci, Phys. Rev. B 94, 014115 (2016). 9 FIG. S2. Phonon spectra in orthorhombic bismuth monalayers obtained with dynamical matrix from the main text (left) (without account of dipole correction) and using the Phonopy package [8, 9] with 3x3x1 supe...

    2016

Pith tools

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