Quantum Paraelectric Materials

• Quantum paraelectric materials are insulating systems near a ferroelectric instability where quantum zero-point fluctuations prevent long-range polar order.
• Experimental studies reveal soft infrared-active phonon modes and Barrett-law dielectric saturation in perovskites like SrTiO₃ and KTaO₃ as key signatures.
• Recent research extends the concept to metallic and strained systems, highlighting tunability via isotope substitution, THz excitation, and pressure for device applications.

Quantum paraelectric materials are insulating crystals, and in some recent extensions also metallic systems near ferroelectric criticality, that lie extremely close to a ferroelectric instability but do not develop long-range ferroelectric order because quantum zero-point fluctuations stabilize the high-symmetry phase. In the canonical formulation, the polar soft mode softens strongly on cooling, the dielectric susceptibility rises, and the material approaches a quantum critical point at $T=0$, yet the transition is suppressed and the low-temperature state remains paraelectric rather than ferroelectric (Anderson et al., 21 Feb 2025). Across the literature summarized here, this regime is identified not only in perovskites such as SrTiO$_3$ and KTaO$_3$, but also in M-type hexaferrites with triangular-lattice dipoles and in proton-bonded iridates where tunneling dipoles remain dynamically disordered (Zhang et al., 2022).

1. Defining characteristics and experimental signatures

Quantum paraelectricity is defined by the coexistence of an incipient ferroelectric instability with the absence of static long-range polar order at low temperature. In SrTiO$_3$, the soft infrared-active phonon modes soften on cooling, but near $T_c \approx 35$ K quantum fluctuations suppress the transition and the material enters a quantum paraelectric ground state (Shen et al., 25 Jul 2025). In KTaO$_3$, the soft ferroelectric mode softens only down to a finite value and then saturates at low temperature, so the paraelectric phase remains the ground state (Ranalli et al., 2022). In CaTiO$_3$, the dielectric constant follows Curie–Weiss behavior at higher temperature and then saturates near 35 K, consistent with a Barrett-law quantum paraelectric response (Bradarić et al., 2019).

A central thermodynamic signature is low-temperature saturation rather than divergence of the dielectric constant. For dipolar systems described in Barrett form, the dielectric response is written as

$$\varepsilon_r = A + \frac{M}{\frac{1}{2}T_1 \coth\!\left(\frac{T_1}{2T}\right)-T_0},$$

with $T_1$ associated with a quantum tunneling scale and $T_0$ reflecting effective dipole-dipole coupling (Zhang et al., 2022). Closely related Barrett-type expressions are used for protonic dipoles in H$_3$0LiIr$_3$1O$_3$2, where the electric susceptibility takes the form

$_3$3

with $_3$4 K and $_3$5 K in mean-field theory (Wang et al., 2018). The phenomenological content is the same: Curie–Weiss-like behavior crosses over to a plateau when zero-point motion suppresses ordering.

The soft-mode viewpoint is equally central. In conventional displacive ferroelectrics, the zone-center transverse optical mode at $_3$6 softens to zero and condenses into a homogeneous macroscopic polarization. In quantum paraelectrics, the mode remains finite or is renormalized upward by quantum and anharmonic effects, so the lattice remains nonpolar on average (Wang et al., 12 Mar 2026). The Lyddane–Sachs–Teller relation,

$_3$7

is used in first-principles studies of KTaO$_3$8 to connect soft-mode behavior directly to the enhanced dielectric response (Ranalli et al., 2022).

Not all signatures are purely static. In H$_3$9LiIr$_3$0O$_3$1 and D$_3$2LiIr$_3$3O$_3$4, dielectric spectroscopy reveals a relaxation step in $_3$5, strong frequency dependence, broad non-Debye dispersion, and a crossover from thermally activated hopping to quantum tunneling on cooling, leading to the conclusion that these materials are simultaneously Kitaev quantum-spin-liquid candidates and quantum paraelectrics (Geirhos et al., 2020). This broadens the operational definition beyond simple dielectric plateaus toward tunneling-dominated dipolar dynamics.

2. Microscopic mechanisms and theoretical descriptions

The microscopic origin of quantum paraelectricity depends on the material class, but the recurring mechanism is competition between a shallow polar instability and quantum motion. In SrTiO$_3$6, the lattice potential becomes strongly anharmonic and can be viewed as a double-well landscape for the soft modes, while quantum fluctuations prevent the ions from freezing into a ferroelectric state (Shen et al., 25 Jul 2025). In M-type BaFe$_3$7O$_3$8, local electric dipoles arise from off-center displacement of Fe$_3$9 ions in the bipyramidal site along the $_3$0-axis; these dipoles remain dynamic because of quantum fluctuations and quantum tunneling, yielding a low-temperature dielectric plateau instead of ferroelectric order (Zhang et al., 2022). In H$_3$1LiIr$_3$2O$_3$3, each proton in an O–H–O bond occupies a double-well potential with minima displaced by about $_3$4 Å from the bond center, creating a local uniaxial dipole whose ordering tendency is overwhelmed by tunneling (Wang et al., 2018).

Several complementary theoretical frameworks are used. Landau–Ginzburg–Devonshire theory is employed to connect dielectric susceptibility, polarization, and nonlinear optical or piezoelectric response in SrTiO$_3$5, with free energy

$_3$6

and dielectric response set by the free-energy curvature (Anderson et al., 21 Feb 2025). In nanoscale electromechanics of SrTiO$_3$7 films, a related LGD-based functional incorporates electrostriction, flexoelectricity, Maxwell stress, deformation-potential coupling, and Vegard strain, with the quadratic coefficient obeying a Barrett-law form

$_3$8

to encode quantum-paraelectric saturation (Morozovska et al., 2011).

A different but widely used approach is to treat the unstable polar mode through quantum lattice dynamics. For BaTiO$_3$9, SrTiO$T_c \approx 35$0, and KTaO$T_c \approx 35$1, a single-particle quantum-mechanical description solves

$T_c \approx 35$2

in a DFT-derived quartic double-well potential

$T_c \approx 35$3

and classifies materials according to whether the zero-point energy lies below the barrier, near the barrier top, or above it (Esswein et al., 2021). In that construction, BaTiO$T_c \approx 35$4 is ferroelectric, SrTiO$T_c \approx 35$5 is the canonical quantum paraelectric, and KTaO$T_c \approx 35$6 lies closer to ordinary paraelectric behavior (Esswein et al., 2021).

For strongly anharmonic perovskites, first-principles work combines the stochastic self-consistent harmonic approximation with machine-learned force fields, and in SrTiO$T_c \approx 35$7 further with random phase approximation electronic structure, to show that the paraelectric phase is stabilized by anharmonic quantum fluctuations even when harmonic phonons predict a ferroelectric instability (Verdi et al., 2022). In KTaO$T_c \approx 35$8, the same SSCHA-based strategy shows that including anharmonic terms is essential to stabilize the spurious imaginary ferroelectric phonon predicted by harmonic DFT, yielding the experimentally observed low-temperature soft-mode plateau (Ranalli et al., 2022).

In protonic systems, the microscopic Hamiltonian is naturally written as a transverse-field Ising model,

$T_c \approx 35$9

where $_3$0 labels the dipole orientation and the transverse field $_3$1 represents proton tunneling (Wang et al., 2018). For H$_3$2LiIr$_3$3O$_3$4, $_3$5 meV greatly exceeds the meV-scale dipolar couplings, yielding a quantum-disordered paraelectric ground state with a dipole excitation gap $_3$6 meV (Wang et al., 2018).

3. Principal material platforms

SrTiO$_3$7 is the canonical quantum paraelectric platform throughout the recent literature. It displays Curie–Weiss-like softening toward an extrapolated ferroelectric temperature scale near 35 K, but long-range ferroelectric order is suppressed by zero-point fluctuations (Li et al., 2018). Its low-frequency dielectric constant exceeds 25,000 in the natural crystal at low temperature, rises to about 42,000 in STO-28 and 82,000 in STO-33 near the tuned quantum critical regime, and under DC bias can exhibit a dynamically tunable linear Pockels coefficient $_3$8 exceeding 500 pm/V at $_3$9 K, with values above 1100 pm/V after oxygen isotope tuning (Anderson et al., 21 Feb 2025). Because of this combination of giant dielectric response, strong soft-mode nonlinearity, and cryogenic tunability, SrTiO$_3$0 underpins work on nonlinear terahertz polaritonics, electro-optics, piezoelectricity, varactors, and proposed microwave three-wave mixers (Shen et al., 25 Jul 2025).

KTaO$_3$1 is a closely related incipient ferroelectric that retains the cubic perovskite structure to low temperature and is frequently treated as a cleaner quantum paraelectric because it has no structural phase transition from room temperature down to millikelvin temperatures (Cheng et al., 2023). First-principles studies show that the harmonic cubic phase is unstable, but anharmonic quantum renormalization removes the ferroelectric instability and produces a temperature-dependent soft-mode plateau below about 30 K (Ranalli et al., 2022). KTaO$_3$2 also serves as a central system for anomalous thermal transport near the ferroelectric quantum critical point, where optical phonons contribute a thermal conductivity scaling $_3$3 with $_3$4 (Bhalla et al., 2021).

CaTiO$_3$5 appears in the literature as another prominent quantum paraelectric host. Its Barrett-fit parameters from prior work are $_3$6, $_3$7, and $_3$8, and the dielectric constant saturates near 35 K (Bradarić et al., 2019). In mixed CaTi$_3$9Ru$$\varepsilon_r = A + \frac{M}{\frac{1}{2}T_1 \coth\!\left(\frac{T_1}{2T}\right)-T_0},$$0O$$\varepsilon_r = A + \frac{M}{\frac{1}{2}T_1 \coth\!\left(\frac{T_1}{2T}\right)-T_0},$$1, a concentration-independent ferromagnetic-like transition near 35 K is reported and interpreted as evidence for coupling between Ru impurity spins and quantum paraelectric fluctuations in the CaTiO$$\varepsilon_r = A + \frac{M}{\frac{1}{2}T_1 \coth\!\left(\frac{T_1}{2T}\right)-T_0},$$2 host (Bradarić et al., 2019).

BaFe$$\varepsilon_r = A + \frac{M}{\frac{1}{2}T_1 \coth\!\left(\frac{T_1}{2T}\right)-T_0},$$3O$$\varepsilon_r = A + \frac{M}{\frac{1}{2}T_1 \coth\!\left(\frac{T_1}{2T}\right)-T_0},$$4 extends the subject beyond perovskites. It is presented as a new quantum paraelectric with triangular-lattice local-dipole geometry, in which off-center Fe$$\varepsilon_r = A + \frac{M}{\frac{1}{2}T_1 \coth\!\left(\frac{T_1}{2T}\right)-T_0},$$5 displacements on FeO$$\varepsilon_r = A + \frac{M}{\frac{1}{2}T_1 \coth\!\left(\frac{T_1}{2T}\right)-T_0},$$6 bipyramids generate dipoles along the $$\varepsilon_r = A + \frac{M}{\frac{1}{2}T_1 \coth\!\left(\frac{T_1}{2T}\right)-T_0},$$7-axis (Zhang et al., 2022). Because the bipyramids form a perfect two-dimensional triangular lattice, geometric frustration combines with quantum tunneling to motivate a possible quantum-dipole-liquid ground state, although the cited work explicitly treats this as a possibility rather than a demonstrated phase (Zhang et al., 2022).

Hydrogen-bonded iridates provide a further non-oxide variant. In H$$\varepsilon_r = A + \frac{M}{\frac{1}{2}T_1 \coth\!\left(\frac{T_1}{2T}\right)-T_0},$$8LiIr$$\varepsilon_r = A + \frac{M}{\frac{1}{2}T_1 \coth\!\left(\frac{T_1}{2T}\right)-T_0},$$9O$T_1$0, first-principles calculations estimate a dipole moment $T_1$1 for each proton-centered bond dipole (Wang et al., 2018). Dielectric spectroscopy on H$T_1$2LiIr$T_1$3O$T_1$4 and D$T_1$5LiIr$T_1$6O$T_1$7 finds glassy dipolar freezing with relaxation rates in the mHz range at low temperature, so the materials are described as both quantum paraelectrics and Kitaev quantum-spin-liquid candidates (Geirhos et al., 2020).

Finally, recent theory extends the language to quantum paraelectric metals, defined as metals near ferroelectric quantum criticality with soft transverse optical phonons. In these systems, a soft TO branch with dispersion

$T_1$8

is taken to generate Rashba-like electron-phonon coupling and anomalous spin transport even for trivial electronic band structures (Kim et al., 2024). This suggests that the concept now spans both insulating and metallic contexts, though the metallic extension remains explicitly tied to proximity to ferroelectric criticality.

4. Tuning, criticality, and phase competition

A recurrent theme is that quantum paraelectrics are exceptionally sensitive to external control parameters. In BaFe$T_1$9O$T_0$0, quantum paraelectricity is tuned by $T_0$1Fe isotope substitution, in-plane compressive strain, and hydrostatic pressure (Zhang et al., 2022). The tunneling scale is modeled as

$T_0$2

so increasing the Fe mass lowers tunneling probability; full $T_0$3Fe substitution is estimated to reduce tunneling probability by about 0.88% (Zhang et al., 2022). Experimentally, 95% $T_0$4Fe replacement and in-plane strain induce a low-temperature peak in the dielectric constant and move the system closer to a critical region, while hydrostatic pressure suppresses the peak and restores plateau-like behavior, pushing the system away from the quantum critical point (Zhang et al., 2022).

In SrTiO$T_0$5, oxygen isotope substitution is an established route toward quantum criticality in the nonlinear optical context. Increasing $T_0$6O content from natural STO to STO-28 and STO-33 raises the dielectric constant peak and enhances both electro-optic and piezoelectric response, with $T_0$7 increasing from $T_0$8 to $T_0$9 to $_3$00 pm/V and $_3$01 from $_3$02 to $_3$03 to $_3$04 pC/N (Anderson et al., 21 Feb 2025). At the same time, the bias field required for maximal response decreases from about 0.080 MV/m in natural STO to 0.036 MV/m and 0.008 MV/m in the isotope-enriched samples (Anderson et al., 21 Feb 2025). A plausible implication is that approaching the ferroelectric quantum critical point amplifies not only dielectric susceptibility but also the coupling of susceptibility into optical and mechanical nonlinearities.

Strain is another especially strong tuning parameter. In epitaxial BaFe$_3$05O$_3$06 films on SrTiO$_3$07, in-plane compressive strain forces Fe ions farther off-center, increases the dipole moment and barrier, and produces a low-temperature dielectric peak around $_3$08 K instead of a simple plateau (Zhang et al., 2022). In the single-particle quantum treatment of perovskites, experimentally accessible strains and volume changes have a much larger effect on the ferroelectric–quantum-paraelectric balance than direct oxygen isotope mass substitution, with a 1% volume increase strongly reshaping the double well in both SrTiO$_3$09 and KTaO$_3$10 (Esswein et al., 2021).

Pressure does not act universally. In BaFe$_3$11O$_3$12, hydrostatic pressure from 0 to 29.2 kbar suppresses the normalized dielectric peak, turns the peak back into a plateau, and at the highest pressure extends plateau behavior to about 15 K (Zhang et al., 2022). This contrasts with the qualitative effect of in-plane strain, showing that the sign of tuning depends on how the perturbation modifies the local dipolar geometry.

Nonequilibrium fields can also traverse or reshape the phase landscape. In SrTiO$_3$13, intense single-cycle THz excitation can dynamically induce a hidden ferroelectric phase, with molecular dynamics simulations giving threshold fields roughly 800 kV/cm along one crystallographic direction and 1100 kV/cm along another, and concluding that a single-cycle THz field can induce a long-lived ferroelectric polarization when the field exceeds about 1 MV/cm (Li et al., 2018). A later study reports a re-entrant sequence in which increasing THz field drives SrTiO$_3$14 from the quantum paraelectric ground state to an intermediate ferroelectric phase and then back to a hidden quantum paraelectric phase above about 500 kV/cm, attributed to coherent population of higher soft-mode eigenstates (Li et al., 2024).

5. Nonequilibrium dynamics, hidden phases, and current controversies

The nonequilibrium physics of quantum paraelectrics is marked by a persistent controversy over whether ultrafast THz driving produces genuine ferroelectric order, hidden paraelectric states, or only enhanced local dipolar correlations. In SrTiO$_3$15, early THz-pump experiments interpreted the growth of non-oscillatory TFISH and the appearance of new low-frequency phonon peaks as evidence for a dynamically induced hidden ferroelectric phase (Li et al., 2018). Subsequent work complicated that picture by finding that intense THz excitation can produce a re-entrant hidden quantum paraelectric phase above $_3$16 kV/cm, characterized by activated antiferrodistortive phonon modes and explained as a coherent superposition of ground, first-order, and second-order soft-mode eigenstates (Li et al., 2024).

Another later study argues that the mechanism of THz-induced symmetry breaking in SrTiO$_3$17 is spatially inhomogeneous and strongly defect dependent. Under single-cycle fields up to $_3$18 MV/cm, short-lived coherent antiferrodistortive modes suppress dipole correlations within $_3$19 ps, while vacancy-rich regions show heavily damped soft and AFD modes together with a defect-induced low-frequency mode at $_3$20–$_3$21 THz that prevents long-range ferroelectric coherence (Cheng et al., 1 Dec 2025). In vacancy-sparse regions, the same work interprets a non-monotonic temperature dependence peaking at $_3$22 K and a softening-then-hardening of collective modes as evidence for transient ferroelectric order, whereas vacancy-rich regions show only monotonic, defect-dominated behavior (Cheng et al., 1 Dec 2025). This directly addresses earlier conflicting interpretations by making defects an active regulator rather than a secondary complication.

KTaO$_3$23 provides an instructive counterexample. THz-pump/SHG experiments observe a long-lived SHG relaxation lasting up to 20 ps at 10 K, but detailed analysis of the coherent soft-mode oscillation finds hardening with fluence well described by a single-well potential,

$_3$24

leading to the conclusion that intense THz pulses up to 500 kV/cm do not drive a global ferroelectric phase (Cheng et al., 2023). Instead, the long-lived SHG background is attributed to moderate dipolar correlation between defect-induced local polar structures or polar nanoregions (Cheng et al., 2023). This has become a central caution in the field: long-lived SHG alone is not sufficient evidence for transient ferroelectric order.

A related equilibrium controversy concerns the microscopic nature of the low-temperature state in SrTiO$_3$25. The standard picture takes it to be an incipient ferroelectric whose $_3$26 TO mode is suppressed by quantum fluctuations. A recent finite-momentum x-ray study instead reports that uniaxial tensile strain at 20 K does not produce a strong zone-center ferroelectric response but stabilizes a hidden polar-acoustic phase with nanoscale polarization modulation, signaled by a vertical diffuse streak extending to $_3$27 rlu and a dramatically renormalized finite-$_3$28 transverse acoustic branch with an additional peak near 0.33 THz (Wang et al., 12 Mar 2026). The authors argue that this hidden phase can mimic ferroelectric-like thermodynamic signatures while differing fundamentally in its collective excitations (Wang et al., 12 Mar 2026). This suggests that some phenomena traditionally attributed to suppressed homogeneous ferroelectricity may instead involve competition with a modulated polar state.

Even proposals for cavity engineering have yielded unexpected results. A multimode continuum treatment of a quantum paraelectric slab in a Fabry–Perot cavity finds that, once the full transverse-mode continuum is included, the cavity suppresses ferroelectric correlations near the metallic boundaries rather than enhancing them, producing a surface-confined blue shift of the soft transverse phonon frequency (Curtis et al., 2023). The effect vanishes at high temperature, indicating a purely quantum mechanical origin (Curtis et al., 2023). This stands in explicit tension with earlier single-mode expectations and underscores that boundary conditions and mode counting are not merely technical details.

6. Functional responses, transport, and device implications

Quantum paraelectricity has become a materials design principle for large cryogenic nonlinearities. In SrTiO$_3$29, the nonlinear electro-optic susceptibility is connected to the product $_3$30 through relations such as

$_3$31

which explain why proximity to a ferroelectric quantum critical point yields giant bias-tunable Pockels and piezoelectric coefficients at cryogenic temperature (Anderson et al., 21 Feb 2025). The same underlying softness makes SrTiO$_3$32 suitable for rf varactors functional at 6 mK, where the capacitance of a quantum-paraelectric parallel-plate device is tuned from about 40 pF to 15 pF, with a maximum around 52 pF near $_3$33 V and an inferred dielectric constant $_3$34 (Apostolidis et al., 2020). Those varactors enable perfect impedance matching and resonator tuning between 167 MHz and 182 MHz, supporting charge sensitivities of $_3$35 and capacitance sensitivities of 0.04 aF/$_3$36 in quantum-dot readout (Apostolidis et al., 2020).

At microwave frequencies, the same large field-tunable permittivity motivates a proposed quantum paraelectric nonlinear dielectric amplifier, or PANDA, based on STO or KTO (Rosenthal et al., 18 Oct 2025). For a nanofabricated parallel-plate capacitor, the effective driven Hamiltonian contains a three-wave mixing term with strength $_3$37 and a quartic Kerr term $_3$38,

$_3$39

The proposal estimates $_3$40 MHz for STO and $_3$41 MHz for KTO, while $_3$42 remains of order $_3$43 Hz, giving $_3$44 (Rosenthal et al., 18 Oct 2025). This suggests that quantum paraelectric dielectrics can function not only as tunable capacitors but as strongly nonlinear yet weak-Kerr mixing elements.

In the terahertz regime, quantum paraelectrics support hybrid light–matter excitations with unusually strong nonlinearity. In SrTiO$_3$45, 2D time-resolved THz Kerr imaging directly tracks bulk phonon-polaritons and shows that the harmonic Drude–Lorentz dielectric response,

$_3$46

describes high-temperature propagation, whereas below about 30 K the response becomes nonperturbative and at 5 K the low-temperature quantum paraelectric phase supports a thresholded, soliton-like transport channel (Shen et al., 25 Jul 2025). The threshold is identified around $_3$47 externally, corresponding to an internal threshold of order 10 kV/cm, and interpreted באמצעות self-induced transparency of an effective two-level soft mode (Shen et al., 25 Jul 2025). This establishes quantum paraelectrics as a platform for nonlinear THz polaritonics and all-optical THz signal control.

Transport phenomena likewise reflect critical lattice fluctuations. In SrTiO$_3$48 and KTaO$_3$49, Kubo-formalism calculations show that a nearly soft transverse optical phonon contributes directly to thermal conductivity and, through TA–TO scattering, produces a low-temperature power law $_3$50 with $_3$51, rather than the usual $_3$52 behavior of acoustic-phonon-dominated insulators (Bhalla et al., 2021). In quantum paraelectric metals, ferroelectric fluctuations generate a Rashba-like phonon-mediated spin-orbit coupling and a spin conductivity quadratic in the wave vector of an inhomogeneous electric field,

$_3$53

with explicit quadrupolar symmetry in components such as $_3$54 (Kim et al., 2024). This extends the influence of ferroelectric quantum criticality from dielectric phenomena into spin transport.

At the nanoscale, quantum paraelectric films exhibit strong electromechanical responses despite being non-piezoelectric in bulk symmetry. In SrTiO$_3$55 thin films with mobile carriers and vacancies, theory predicts that the surface displacement contains contributions from concentration-strain, flexoelectricity, electrostriction, and Maxwell stress,

$_3$56

with a pronounced crossover from linear to quadratic to sub-linear $_3$57 voltage dependence due to dielectric nonlinearity (Morozovska et al., 2011). The response is strongly size dependent through the ratio of film thickness or tip radius to the Debye screening radius (Morozovska et al., 2011).

Quantum paraelectricity also couples to other collective sectors. In CaTi$_3$58Ru$_3$59O$_3$60, impurity moments are proposed to couple to quantum paraelectric fluctuations through dynamical multiferroicity, schematically $_3$61, yielding a concentration-independent ferromagnetic-like transition near 35 K (Bradarić et al., 2019). In H$_3$62LiIr$_3$63O$_3$64, proton dipole fluctuations renormalize magnetic exchange interactions and are argued to help stabilize Kitaev quantum spin liquid behavior (Wang et al., 2018). A plausible implication is that quantum paraelectricity should be viewed less as an isolated dielectric anomaly than as a low-energy sector capable of mediating optical, mechanical, magnetic, and transport phenomena across multiple classes of quantum materials.