Pith. sign in

REVIEW 3 major objections 5 minor 54 references

Dipolar Nematic State in Relaxor Ferroelectrics

T0 review · 3 major / 5 minor · reviewed 2026-08-05 · deepseek-v4-flash

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

Pith's one-line read Molecular dynamics across Pb-, Bi-, and Ba-based relaxors shows a universal dipolar nematic state—long-range orientational order of local polarizations without cluster alignment—whose destruction is tracked by a single skewness-based order

desk verdict A genuinely ambitious computational paper that reproduces a remarkable range of relaxor phenomenology and proposes a new universal nematic picture, but the load-bearing claim of equilibrium orientational order rests on 1 ns trajectories that are not validated against slower director fluctuations. read the letter →

arxiv 2509.01464 v1 pith:E377EAIR submitted 2025-09-01 cond-mat.mtrl-sci physics.comp-ph

Yuan-Jinsheng Liu , Tyler C. Sterling , Shi Liu This is my paper
classification cond-mat.mtrl-sciphysics.comp-ph
keywords relaxorferroelectricsdipolarnematicstatemoleculardynamicsmachine-learnedinteratomicpotentialorientationalorderpolarizationautocorrelationmastercurvelead-freerelaxors
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 claims that relaxor ferroelectrics are governed by a universal 'dipolar nematic' state: local polarizations keep long-range orientational coherence even though neighboring dipoles are not aligned into clusters. Using large-scale molecular dynamics driven by a first-principles machine-learned interatomic potential, the authors reproduce experimental signatures of three chemically distinct relaxors—Pb-based PIN-PMN-PT, Bi-based BNT, and Ba-based BZT—then show that the orientational order is visible in the statistics of how each unit-cell polarization persists in time rather than in conventional cluster analysis. They define an order parameter from the skewness of the distribution of local polarization autocorrelation values; plotted against T/Tm, all three systems fall on one master curve that decays to zero at the dielectric peak. If true, this replaces the polar-nanoregion and polar-nanodomain pictures with a single statistical mechanism that also explains the diffuse phase transition, frequency-dependent dielectric dispersion, and reversible giant piezoelectricity.

What carries the argument

The central tool is a sliding-window analysis of each unit cell's local polarization trajectory. For each cell, the trajectory is divided into 200 ps windows; at the end of each window one records the autocorrelation value a_k between the polarization at the window start and at the window end. The distribution of these values gives a mean mu_j (directional memory) and standard deviation sigma_j (temporal variability), which sorts cells into Type-III persistent, Type-II diffuse, and Type-I stochastic dynamics. The ensemble polarization sphere—all instantaneous local polarization vectors projected onto a unit sphere—then reveals whether the preferred directions have global orientational order.

What would settle it

Run the same analysis on 10 ns trajectories or multiple independent seeds at 300 K for PIN-PMN-PT: if the ensemble polarization sphere becomes isotropic or the preferred axis drifts continuously across the sphere over time, the claimed nematic order is a finite-time artifact. Alternatively, a macroscopic probe sensitive to orientational order—such as optical birefringence or anisotropy of diffuse neutron scattering—that stays nonzero above Tm while the mu-distribution skewness order parameter decays would contradict the claimed one-to-one link to the dielectric peak.

Watch

Extended reading notes

Core claim

On the paper's own terms, the discovery is that relaxor behavior does not come from nanoscale clusters of aligned dipoles. Instead, in every system simulated, each unit cell's polarization fluctuates around one or more preferred directions, and those preferred directions are collectively aligned along a global axis (or, in BNT, eight symmetry-related axes) even when the macroscopic polarization is zero and the spatial polarization correlation decays to zero. This is the signature of a dipolar nematic: orientational order without local alignment. The paper further claims that the thermal destruction of this order, quantified by the skewness S_mu of the distribution of per-cell autocorrelation

Load-bearing premise

The load-bearing assumption is that each 1 ns simulation trajectory is long enough for time averages to equilibrate, so the anisotropic polarization sphere reflects a true thermodynamic nematic state rather than a slowly wandering global director in an otherwise disordered system.

Editorial extensions

If this is right

  • The relaxor diffuse transition is two-stage: macroscopic polarization disappears first, while orientational order survives until Tm; the dielectric peak marks the loss of nematic order, not the disappearance of polar domains.
  • Frequency-dependent dielectric dispersion follows from the coexistence of three dynamically distinct cell populations with different relaxation times; high frequencies freeze out the persistent Type-III cells.
  • Giant reversible piezoelectricity arises because the nematic ground state is orientationally soft: an electric field biases the population balance among cell types, and removing the field restores the original mu-distribution.
  • The normalized skewness (1 - S_mu) can serve as a practical order parameter for comparing relaxor behavior across chemistries, with the dielectric maximum predicted where it crosses zero.
  • Anisotropic diffuse neutron scattering features—butterfly shapes, ellipsoidal intensity, and 1/q^2 tails—are signatures of the Pb/O sublattice orientational correlations, not of compact polar clusters.

Reading between the lines

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

  • Not claimed by the paper, but the same autocorrelation-skewness order parameter could be applied to other disordered polar or magnetic materials—dipolar glasses, supercooled molecular liquids, spin ices—where hidden orientational order may coexist with local disorder.
  • If the master curve is truly universal, it makes a concrete prediction: for any new relaxor composition, knowing Tm alone would locate its normalized temperature on the same curve, so one could infer the microscopic order parameter evolution directly from dielectric measurements.
  • The structural-memory picture suggests a design rule for high-performance relaxors—maximize the population of persistent Type-III cells with narrow directional spread rather than trying to grow large aligned domains—an inversion of conventional domain-engineering strategy.
  • A field-cycling experiment with variable field strength and rate, measuring how quickly the mu-distribution recovers after field removal, could distinguish the nematic restoration mechanism from ordinary domain back-switching.
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 / 5 minor

Summary. The manuscript reports large-scale NPT molecular dynamics simulations of three relaxor families (PIN-PMN-PT, BNT, and BZT) using a modified attention-free UniPero machine-learning potential trained on DFT data. It reproduces several experimental signatures: diffuse dielectric maxima, anisotropic diffuse neutron scattering, the waterfall effect in the dynamic structure factor, and effective d33 values above 1200 pC/N. From the trajectories, the authors perform a sliding-window analysis of local polarization autocorrelations and derive per-cell persistence and variance measures (μ, σ). They interpret the resulting anisotropic ensemble polarization spheres as evidence of a universal dipolar nematic state: long-range orientational order of local polarizations without local collinear alignment. They further propose that the skewness Sμ of the μ-distribution, normalized and expressed as 1−Sμ, collapses the thermal evolution of all three systems onto a single master curve when plotted against T/Tm, and they connect this nematic order to reversible giant piezoelectricity and to the diffuse phase transition.

Significance. If established, the proposed dipolar nematic state would be a significant conceptual shift in relaxor physics, potentially unifying lead-based and lead-free systems and explaining the coexistence of local disorder with global orientational coherence. The paper's strengths are substantial: the MLIP is stated to be publicly available; the simulations reproduce multiple independent experimental observables; and the statistical analysis of local autocorrelation distributions is inventive and clearly presented. I find no circularity in deriving the nematic order from MD trajectories trained on DFT data. However, the central claim of equilibrium long-range orientational order rests on 1 ns trajectories without ergodicity checks, and the universal master curve relies on normalization choices that are not fully specified. These points must be addressed before the universality claim is convincing.

major comments (3)
  1. [Methods C / Fig. 2f / Fig. S7] The load-bearing evidence for long-range orientational order is the anisotropic ensemble polarization sphere and the nematic order parameter S, both computed from a single 1 ns NPT trajectory per temperature (Methods C). The quasi-static pair correlations for τ>2 ps (Fig. S6) constrain local pair reorientation, not collective director rotation. If the director performs a slow random walk on the unit sphere with correlation time ≫1 ns, the trajectory samples only a small patch of orientations, producing an artificially anisotropic sphere and nonzero S even if the equilibrium ensemble is isotropic. No independent runs, block averages, error bars, or director time series are reported. Please provide (i) S computed from a traceless nematic tensor as a function of time for trajectories several times longer, or from multiple independent 1 ns runs; (ii) the director trajectory (e.g., the princi
  2. [Discussion / Fig. 5d] The universal master curve is constructed from 1−Sμ with Sμ 'normalized with respect to its high-temperature value,' and each system is rescaled by its own Tm. As written, the normalized quantity is zero at the high-temperature reference by construction, so the collapse in Fig. 5d does not by itself demonstrate universality unless the raw skewness values are shown to have the same functional form before normalization. Please define Sμ explicitly (formula), report raw skewness for every system and temperature, state the reference temperature used for each normalization, and show the sensitivity of the collapse to the choice of reference. The stretched-exponential fit, β exp(−(T/Tm)^n), introduces additional parameters β and n; these should be reported with uncertainties. With per-system Tm, a high-temperature reference, β, and n, the fit has enough degrees of freedom that a visually good
  3. [Fig. 1f / Fig. 2f] The manuscript should reconcile the claim of long-range orientational order with the fact that C_p(r) decays to zero at 300 K. A nematic state with ±n symmetry can indeed have vanishing vector correlation at long distances while retaining quadrupolar order; the text hints at this but does not state it. More importantly, the nematic order parameter S used in Fig. S7 is never defined in the main text, and its connection to C_p(r) is not given. Please state the definition of S (e.g., largest eigenvalue of the traceless second-moment tensor of the normalized polarization orientation distribution) and show that S is robust to trajectory length and system size. A finite-size check (for instance 16^3 vs 32^3 supercells) would materially strengthen the claim that the observed anisotropy is true long-range order rather than a finite-size or finite-time fluctuation.
minor comments (5)
  1. [Sec. II A] There is a typo: 'the the butterfly-shaped' should read 'the butterfly-shaped'. Also, Eq. (1) would benefit from a cleaner typesetting of the B-site term so that the reader can immediately see the summation convention.
  2. [Methods C] The sliding-window construction is not fully specified: please state the stride between successive windows and whether the windows overlap. This determines the effective number of statistically independent samples in the {a_k} distribution and is relevant to the reported μ and σ values.
  3. [Fig. 5d / text] The caption labels the ordinate ΔSμ = 1 − Sμ, while the text states that Sμ is normalized with respect to its high-temperature value. Please align the notation and state explicitly what is plotted: raw 1−Sμ, or 1−Sμ/Sμ(ref)?
  4. [Data Availability / Methods A] Methods A states that the training dataset, hyperparameters, and model are publicly available, but the Data Availability section says data are available from the corresponding author upon reasonable request. Please make the two statements consistent.
  5. [Sec. II C / Fig. 2c] The thresholds separating Type-I, Type-II, and Type-III cells are described only qualitatively ('low μ', 'high μ', etc.). Please report the exact classification criteria, either in the main text or in the Supplementary Methods, so that the dynamic-type fractions can be reproduced.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the claimed nematic state and master curve are obtained from MD trajectories generated by a DFT-trained potential and validated against external experimental observables; normalization and empirical fitting do not constitute definitional circularity.

full rationale

The central structural claim (dipolar nematic state) is a direct analysis of MD trajectories produced by the UniPero machine-learned potential (Ref. [17]) trained on DFT data. Although Ref. [17] shares authors with this paper, it functions as an independent computational model: the present work validates it against experimental diffuse scattering, phonon dispersion, dielectric anomalies, and piezoelectric response, and the nematic interpretation is not an input to the potential. The universal order parameter is defined as a normalized skewness of the local-polarization autocorrelation distribution (Methods C, Fig. 5a–d). Normalizing by each system's high-temperature skewness and by its own T_m is a standard data-collapse procedure; it does not force the observed shape of the curve or the shared stretched-exponential decay, which is explicitly described as empirical. The concern about 1 ns trajectories and possible slow director reorientation is a statistical equilibration issue, not a circularity of the derivation chain. No equation in the paper reduces by construction to a fitted input, and no load-bearing uniqueness claim rests solely on author self-citation.

Assumptions & free parameters 5 free parameters · 8 assumptions · 1 invented entities

The central claim rests on (i) the accuracy and transferability of a neural-network potential fitted to DFT data, (ii) the ergodicity of 1 ns MD trajectories, (iii) the supercell size, (iv) the definition of local polarization via averaged Born charges, and (v) the validity of skewness as an order parameter. The universal curve adds the stretched-exponential form and per-system Tm rescaling as fitted elements.

free parameters (5)
  • Modified UniPero neural network weights = Unknown (network weights trained on 19,773 DFT configurations)
    All MD trajectories depend on this ML potential; it is fitted to first-principles data, not derived analytically. The claim 'first-principles-based' rests on this surrogate.
  • Sliding-window duration Delta T = 200 ps
    Chosen by hand in Methods C; it directly shapes the mu and sigma distributions and therefore the skewness order parameter.
  • Stretched exponential parameters beta and n = Not reported
    The universal master curve is described as an 'empirical stretched exponential function' beta exp(-(T/Tm)^n); beta and n are fitted to the collapsed data.
  • Per-system Tm values = 340 K (PIN-PMN-PT), 640 K (BNT), 210 K (BZT)
    Tm is extracted from each simulated dielectric peak and used to rescale temperature, which contributes to producing the apparent universal collapse.
  • High-temperature normalization reference for S_mu = Not specified
    S_mu is normalized with respect to its high-temperature value, a choice that trivially forces the order parameter to start near unity and decay to zero.
assumptions (8)
  • domain assumption The modified UniPero neural network potential accurately represents the potential energy surface of the three relaxors and extrapolates to compositions not explicitly in the training set.
    Invoked in Methods A; all MD results depend on this transferability.
  • domain assumption 1 ns NPT MD trajectories reach equilibrium and time averages equal ensemble averages.
    Invoked in Methods C and Fig. 2f; the polarization sphere anisotropy is interpreted as equilibrium order.
  • domain assumption A 32x32x32 supercell (163,840 atoms) is large enough to avoid finite-size artifacts in long-range correlations.
    Methods A; the supercell size sets the maximum correlation length probed.
  • domain assumption Local polarization from Eq. (1) using averaged Born effective charges captures the essential polar degrees of freedom.
    Methods B; all subsequent correlation and skewness analyses rely on this definition.
  • ad hoc to paper Skewness of the mu-distribution is a valid order parameter for the dipolar nematic state.
    Introduced in Discussion; no independent derivation shows skewness is the correct thermodynamic order parameter.
  • ad hoc to paper The stretched exponential form is the correct universal functional form for the master curve.
    Discussion and Fig. 5d; it is an empirical fit, not derived from the model.
  • domain assumption The experimental BNT dielectric data provided by S. Deng are reliable.
    Used in Fig. 3a inset as validation for the BNT simulations.
  • domain assumption The classical approximation and Gaussian energy resolution used in neutron scattering calculations are adequate for reproducing the experimental diffuse scattering.
    Methods D; affects the comparison in Figs. 1b-c.
invented entities (1)
  • Dipolar nematic state (including the multipolar nematic variant for BNT)
    purpose: Proposed new phase of relaxor ferroelectrics characterized by long-range orientational order without local dipolar alignment; used to explain diffuse transitions, frequency dispersion, and reversible giant piezoelectricity.
    The state is inferred from MD trajectories and time-averaged polarization spheres. No direct experimental signature unique to this state is proposed that would distinguish it from polar nanoregion or slush models.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Dipolar Nematic State in Relaxor Ferroelectrics." pith.science (2026). https://pith.science/paper/E377EAIR

@misc{pith2026250901464,
  author       = {Pith},
  title        = {Pith review of: Dipolar Nematic State in Relaxor Ferroelectrics},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/E377EAIR}},
  note         = {Machine review of arXiv:2509.01464}
}
read the original abstract

Relaxor ferroelectrics exhibit exceptional dielectric and electromechanical properties, yet their microscopic origins remain elusive due to the interplay of hierarchical polar structures and chemical complexity. While models based on polar nanoregions or nanodomains offer valuable phenomenological insights, they often lack the first-principles predictive capability necessary for quantitatively describing functional properties such as piezoelectric coefficients. Here, we use large-scale molecular dynamics simulations, enabled by a universal first-principles-based machine-learning interatomic potential, to investigate atomic-scale polar dynamics in canonical Pb-, Bi-, and Ba-based relaxors. Across all systems, we uncover a universal dipolar nematic state, characterized by long-range orientational order of local polarizations without local alignment, challenging conventional polar cluster-based paradigms. We introduce a universal order parameter, derived from the skewness of the distributions of the local polarization autocorrelation functions, that captures the thermal evolution of both lead-based and lead-free systems within a single master curve. This nematic order, and its robust structural memory under electric field cycling, underpins key relaxor phenomena, including diffuse phase transition, frequency-dependent dielectric dispersion, and reversible giant piezoelectricity. Our findings establish a unified microscopic framework for relaxors and present a broadly applicable statistical approach to understanding complex disordered materials.

Figures

Figures reproduced from arXiv: 2509.01464 by the authors.

Figure 1
Figure 1. MD simulations of the lead-based relaxor 0.24PIN–0.42PMN–0.34PT. a, Simu￾lated static dielectric permittivity (ϵr) and total polarization as a function of temperature. The relaxor exhibits a broad dielectric peak at Tm ≈ 340 K and a gradual polarization decay, characteristic of a diffuse phase transition. b, Simulated elastic neutron diffuse scattering patterns in the (H, K, 0) plane at 300 K. The patterns reproduce… view at source ↗
Figure 2
Figure 2. Dipolar nematic state in the Pb-based relaxor 0.24PIN–0.42PMN–0.34PT. a, Schematic of the sliding-window analysis used to construct the distribution {Ak(∆T)} from the local polarization autocorrelation function, based on the polarization trajectory of unit cell j. This distri￾bution yields the mean (µj , dashed line) and standard deviation (σj , shaded region), which together characterize the polarization dynamics o… view at source ↗
Figure 3
Figure 3. Multipolar nematic state in the Bi-based relaxor Bi0.5Na0.5TiO3. a, Simulated static dielectric permittivity and total polarization of Bi0.5Na0.5TiO3 as functions of temperature. The simulation reproduces key experimental anomalies (inset): a dielectric hump at Td ≈ 320 K and a broad peak at Tm ≈ 640 K. b, Scatter plot of the {µj , σj} distributions at 300 K, showing a dominance of unit cells with Type-III dynamics,… view at source ↗
Figures from the paper (2 more)
Figure 4
Figure 4. Figure 4: Dipolar nematic state in the Ba-based relaxor BaZr0.3Ti0.7O3. a, Simulated dielectric permittivity and total polarization of BaZr0.3Ti0.7O3 as functions of temperature, showing a broad dielectric peak at Tm ≈ 210 K. b, Scatter plot of the {µj , σj} distributions at 140…
Figure 5
Figure 5. Figure 5: Universal nematic ordering and its role in functional properties. a–c, Evolution of the µ-distribution upon heating for PIN-PMN-PT, BNT, and BZT, respectively. All three systems exhibit a universal trend where the peak shifts to lower µ values and broadens, correspondi…

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

54 extracted references · 54 canonical work pages

  1. [1]

    multipolar nematic state

    crystallographic axis (Fig. 2f). This coexistence of local polar disorder and sustained ori- entational coherence is indicative of a dipolar nematic state, analogous to the nematic phase in liquid crystals, where molecules lack positional order but align collectively along a preferred axis [29]. Our calculations of the nematic order parameter S in 0.24–0....

    work page

  2. [2]

    Park and T

    S.-E. Park and T. R. Shrout, Ultrahigh strain and piezoelectric behavior in relaxor based ferro- electric single crystals, J. Appl. Phys. 82, 1804 (1997)

    work page 1997

  3. [3]

    Pandya, J

    S. Pandya, J. Wilbur, J. Kim, R. Gao, A. Dasgupta, C. Dames, and L. W. Martin, Pyroelectric energy conversion with large energy and power density in relaxor ferroelectric thin films, Nat. Mater. 17, 432 (2018)

    work page 2018

  4. [4]

    Kutnjak, J

    Z. Kutnjak, J. Petzelt, and R. Blinc, The giant electromechanical response in ferroelectric relaxors as a critical phenomenon, Nature 441, 956 (2006)

    work page 2006

  5. [5]

    F. Li, B. Wang, X. Gao, D. Damjanovic, L.-Q. Chen, and S. Zhang, Ferroelectric materials toward next-generation electromechanical technologies, Science 389, eadn4926 (2025)

    work page 2025

  6. [6]

    F. Li, M. J. Cabral, B. Xu, Z. Cheng, E. C. Dickey, J. M. LeBeau, J. Wang, J. Luo, S. Taylor, W. Hackenberger, L. Bellaiche, Z. Xu, L.-Q. Chen, T. R. Shrout, and S. Zhang, Giant piezoelec- tricity of Sm-doped Pb(Mg 1/3Nb2/3)O3-PbTiO3 single crystals, Science 364, 264 (2019)

    work page 2019

  7. [7]

    Rojac, Piezoelectric response of disordered lead-based relaxor ferroelectrics, Commun

    T. Rojac, Piezoelectric response of disordered lead-based relaxor ferroelectrics, Commun. Mater. 4, 12 (2023)

    work page 2023

  8. [8]

    Burns and B

    G. Burns and B. Scott, Index of refraction in ‘dirty’ displacive ferroelectrics, Solid State Commun. 13, 423 (1973)

    work page 1973

Show all 54 references
  1. [9]

    Burns and F

    G. Burns and F. Dacol, Glassy polarization behavior in ferroelectric compounds Pb(Mg1/3Nb2/3)O3 and Pb(Zn1/3Nb2/3)O3, Solid State Commun. 48, 853 (1983)

    1983

  2. [10]

    Takenaka, I

    H. Takenaka, I. Grinberg, S. Liu, and A. M. Rappe, Slush-like polar structures in single-crystal relaxors, Nature 546, 391 (2017)

    2017

  3. [11]

    J. Kim, A. Kumar, Y. Qi, H. Takenaka, P. J. Ryan, D. Meyers, J.-W. Kim, A. Fernandez, Z. Tian, A. M. Rappe, J. M. LeBeau, and L. W. Martin, Coupled polarization and nanodomain evolution underpins large electromechanical responses in relaxors, Nat. Phys. 18, 1502 (2022)

    2022

  4. [12]

    Westphal, W

    V. Westphal, W. Kleemann, and M. D. Glinchuk, Diffuse phase transitions and random-field- induced domain states of the “relaxor” ferroelectric PbMg 1/3Nb2/3O3, Phys. Rev. Lett. 68, 847 (1992)

    1992

  5. [13]

    Phelan, C

    D. Phelan, C. Stock, J. A. Rodriguez-Rivera, S. Chi, J. Le˜ ao, X. Long, Y. Xie, A. A. Bokov, Z.-G. Ye, P. Ganesh, and P. M. Gehring, Role of random electric fields in relaxors, Proc. Natl. Acad. Sci. USA 111, 1754 (2014). 18

    2014

  6. [14]

    Fu and R

    H. Fu and R. E. Cohen, Polarization rotation mechanism for ultrahigh electromechanical response in single-crystal piezoelectrics, Nature 403, 281 (2000)

    2000

  7. [15]

    Noheda, Structure and high-piezoelectricity in lead oxide solid solutions, Curr

    B. Noheda, Structure and high-piezoelectricity in lead oxide solid solutions, Curr. Opin. Solid State Mater. Sci. 6, 27 (2002)

    2002

  8. [16]

    Y. M. Jin, Y. U. Wang, A. G. Khachaturyan, J. F. Li, and D. Viehland, Adaptive ferroelectric states in systems with low domain wall energy: Tetragonal microdomains, J. Appl. Phys. 94, 3629 (2003)

    2003

  9. [17]

    Y. M. Jin, Y. U. Wang, A. G. Khachaturyan, J. F. Li, and D. Viehland, Conformal miniaturization of domains with low domain-wall energy: Monoclinic ferroelectric states near the morphotropic phase boundaries, Phys. Rev. Lett. 91, 197601 (2003)

    2003

  10. [18]

    J. Wu, J. Yang, Y.-J. Liu, D. Zhang, Y. Yang, Y. Zhang, L. Zhang, and S. Liu, Universal interatomic potential for perovskite oxides, Phys. Rev. B 108, L180104 (2023)

    2023

  11. [19]

    F. Li, S. Zhang, T. Yang, Z. Xu, N. Zhang, G. Liu, J. Wang, J. Wang, Z. Cheng, Z.-G. Ye, J. Luo, T. R. Shrout, and L.-Q. Chen, The origin of ultrahigh piezoelectricity in relaxor-ferroelectric solid solution crystals, Nat. Commun. 7, 13807 (2016)

    2016

  12. [20]

    Eremenko, V

    M. Eremenko, V. Krayzman, A. Bosak, H. Y. Playford, K. W. Chapman, J. C. Woicik, B. Ravel, and I. Levin, Local atomic order and hierarchical polar nanoregions in a classical relaxor ferroelectric, Nat. Commun. 10, 2728 (2019)

    2019

  13. [21]

    Kumar, J

    A. Kumar, J. N. Baker, P. C. Bowes, M. J. Cabral, S. Zhang, E. C. Dickey, D. L. Irving, and J. M. LeBeau, Atomic-resolution electron microscopy of nanoscale local structure in lead-based relaxor ferroelectrics, Nat. Mater. 20, 62 (2020)

    2020

  14. [22]

    J. Kim, Y. Qi, A. Kumar, Y.-L. Tang, M. Xu, H. Takenaka, M. Zhu, Z. Tian, R. Ramesh, J. M. LeBeau, A. M. Rappe, and L. W. Martin, Size-driven phase evolution in ultrathin relaxor films, Nat. Nanotechnol. 20, 478 (2025)

    2025

  15. [23]

    Zhang, J

    S. Zhang, J. Luo, W. Hackenberger, and T. R. Shrout, Characterization of Pb(In1/2Nb1/2)O3–Pb(Mg1/3Nb2/3)O3–PbTiO3 ferroelectric crystal with enhanced phase transi- tion temperatures, J. Appl. Phys. 104, 064106 (2008)

    2008

  16. [24]

    X. Liu, S. Zhang, J. Luo, T. R. Shrout, and W. Cao, Complete set of material constants of Pb(In1/2Nb1/2)O3–Pb(Mg1/3Nb2/3)O3–PbTiO3 single crystal with morphotropic phase boundary composition, J. Appl. Phys. 106, 074112 (2009)

    2009

  17. [25]

    Q. Li, S. Danilkin, G. Deng, Z. Li, R. L. Withers, Z. Xu, and Y. Liu, Soft phonon modes and diffuse 19 scattering in Pb(In 1/2Nb1/2)O3-Pb(Mg1/3Nb2/3)O3-PbTiO3 relaxor ferroelectrics, J. Materiomics 4, 345 (2018)

    2018

  18. [26]

    Bosak, D

    A. Bosak, D. Chernyshov, S. Vakhrushev, and M. Krisch, Diffuse scattering in relaxor ferroelectrics: true three-dimensional mapping, experimental artefacts and modelling, Acta Crystallogr. A68, 117 (2011)

    2011

  19. [27]

    M. J. Krogstad, P. M. Gehring, S. Rosenkranz, R. Osborn, F. Ye, Y. Liu, J. P. C. Ruff, W. Chen, J. M. Wozniak, H. Luo, O. Chmaissem, Z.-G. Ye, and D. Phelan, The relation of local order to material properties in relaxor ferroelectrics, Nat. Mater. 17, 718 (2018)

    2018

  20. [28]

    Chetverikov, A

    Y. Chetverikov, A. Naberezhnov, S. Vakhrushev, B. Dorner, and A. Ivanov, Inelastic and critical neutron scattering in the ergodic phase of the relaxor ferroelectric Pb(Mg 1/3Nb2/3)O3, Appl. Phys. A 74, s989 (2002)

    2002

  21. [29]

    G. Xu, G. Shirane, J. R. D. Copley, and P. M. Gehring, Neutron elastic diffuse scattering study of Pb(Mg1/3Nb2/3)O3, Phys. Rev. B 69, 064112 (2004)

    2004

  22. [30]

    Khoo, Liquid crystals (John Wiley & Sons, 2022)

    I.-C. Khoo, Liquid crystals (John Wiley & Sons, 2022)

    2022

  23. [31]

    V. V. Shvartsman and D. C. Lupascu, Lead-free relaxor ferroelectrics, J. Am. Ceram. Soc. 95, 1 (2011)

    2011

  24. [32]

    X. Zhou, G. Xue, H. Luo, C. R. Bowen, and D. Zhang, Phase structure and properties of sodium bismuth titanate lead-free piezoelectric ceramics, Prog. Mater. Sci. 122, 100836 (2021)

    2021

  25. [33]

    Sakata and Y

    K. Sakata and Y. Masuda, Ferroelectric and antiferroelectric properties of (Na0.5Bi0.5)TiO3-SrTiO3 solid solution ceramics, Ferroelectrics 7, 347 (1974)

    1974

  26. [34]

    Hiruma, H

    Y. Hiruma, H. Nagata, and T. Takenaka, Thermal depoling process and piezoelectric properties of bismuth sodium titanate ceramics, J. Appl. Phys. 105, 084112 (2009)

    2009

  27. [35]

    B. N. Rao, R. Datta, S. S. Chandrashekaran, D. K. Mishra, V. Sathe, A. Senyshyn, and R. Ranjan, Local structural disorder and its influence on the average global structure and polar properties in Na0.5Bi0.5TiO3, Phys. Rev. B 88, 224103 (2013)

    2013

  28. [36]

    I. P. Pronin, P. P. Syrnikov, V. A. Isupov, V. M. Egorov, and N. V. Zaitseva, Peculiarities of phase transitions in sodium-bismuth titanate, Ferroelectrics 25, 395 (1980)

    1980

  29. [37]

    Takenaka, K.-i

    T. Takenaka, K.-i. M. Kei-ichi Maruyama, and K. S. Koichiro Sakata, (Bi 1/2Na1/2)TiO3-BaTiO3 system for lead-free piezoelectric ceramics, Jpn. J. Appl. Phys. 30, 2236 (1991)

    1991

  30. [38]

    C.-S. Tu, I. G. Siny, and V. H. Schmidt, Sequence of dielectric anomalies and high-temperature relaxation behavior in Na 1/2Bi1/2TiO3, Phys. Rev. B 49, 11550 (1994). 20

    1994

  31. [39]

    Suchanicz, K

    J. Suchanicz, K. Roleder, A. Kania, and J. Ha´ naderek, Electrostrictive strain and pyroeffect in the region of phase coexistence in Na 0.5Bi0.5TiO3, Ferroelectrics 77, 107 (1988)

    1988

  32. [40]

    Suchanicz and J

    J. Suchanicz and J. Kwapulinski, X-ray diffraction study of the phase transitions in Na0.5Bi0.5TiO3, Ferroelectrics 165, 249 (1995)

    1995

  33. [41]

    B. N. Rao, A. N. Fitch, and R. Ranjan, Ferroelectric-ferroelectric phase coexistence in Na1/2Bi1/2TiO3, Phys. Rev. B 87, 060102 (2013)

    2013

  34. [42]

    Zhang, H

    D. Zhang, H. Bi, F.-Z. Dai, W. Jiang, X. Liu, L. Zhang, and H. Wang, Pretraining of attention- based deep learning potential model for molecular simulation, npj Comput. Mater. 10, 94 (2024)

    2024

  35. [43]

    Zhang, H

    Y. Zhang, H. Wang, W. Chen, J. Zeng, L. Zhang, H. Wang, and W. E, DP-GEN: A concur- rent learning platform for the generation of reliable deep learning based potential energy models, Comput. Phys. Commun. 253, 107206 (2020)

    2020

  36. [44]

    J. Wu, J. Yang, L. Ma, L. Zhang, and S. Liu, Modular development of deep potential for complex solid solutions, Phys. Rev. B 107, 144102 (2023)

    2023

  37. [45]

    Plimpton, Fast parallel algorithms for short-range molecular synamics, J

    S. Plimpton, Fast parallel algorithms for short-range molecular synamics, J. Comput. Phys. 117, 1 (1995)

    1995

  38. [46]

    Liu and R

    S. Liu and R. E. Cohen, Origin of stationary domain wall enhanced ferroelectric susceptibility, Phys. Rev. B 95, 094102 (2017)

    2017

  39. [47]

    Boresch, P

    S. Boresch, P. H¨ ochtl, and O. Steinhauser, Studying the dielectric properties of a protein solution by computer simulation, J. Phys. Chem. B 104, 8743 (2000)

    2000

  40. [48]

    Ponomareva, L

    I. Ponomareva, L. Bellaiche, and R. Resta, Dielectric anomalies in ferroelectric nanostructures, Phys. Rev. Lett. 99, 227601 (2007)

    2007

  41. [49]

    Ponomareva, L

    I. Ponomareva, L. Bellaiche, and R. Resta, Relation between dielectric responses and polarization fluctuations in ferroelectric nanostructures, Phys. Rev. B 76, 235403 (2007)

    2007

  42. [50]

    Ponomareva, L

    I. Ponomareva, L. Bellaiche, T. Ostapchuk, J. Hlinka, and J. Petzelt, Terahertz dielectric response of cubic BaTiO 3, Phys. Rev. B 77, 012102 (2008)

    2008

  43. [51]

    G. L. Squires, Introduction to the Theory of Thermal Neutron Scattering, 3rd ed. (Cambridge University Press, 2012)

    2012

  44. [52]

    N. J. Weadock, T. C. Sterling, J. A. Vigil, A. Gold-Parker, I. C. Smith, B. Ahammed, M. J. Krogstad, F. Ye, D. Voneshen, P. M. Gehring, et al., The nature of dynamic local order in CH3NH3PbI3 and CH3NH3PbBr3, Joule 7, 1051 (2023)

    2023

  45. [53]

    V. F. Sears, Neutron scattering lengths and cross sections, Neutron news 3, 26 (1992). 21

    1992

  46. [54]

    T. C. Sterling, pynamic-structure-factor (2023). 22 Figure 1. MD simulations of the lead-based relaxor 0.24PIN–0.42PMN–0.34PT. a , Simu- lated static dielectric permittivity ( ϵr) and total polarization as a function of temperature. The relaxor exhibits a broad dielectric peak...

    2023

Pith tools

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