Theory of central peak and acoustic anomaly in cubic BaTiO3 close to ferroelectric transition

This paper presents a Ginzburg-Landau theory demonstrating that electrostrictive coupling between polarization and elastic displacement in cubic BaTiO$_3$ generates a central peak in displacement correlations and induces singularities in frequency-dependent elastic moduli, sound velocity, and attenuation near the ferroelectric transition.

The Gist

Here is an explanation of the paper, translated from complex physics jargon into a story you can understand.

The Big Picture: A Crystal That "Sighs" Before It Changes

Imagine a block of Barium Titanate (BaTiO₃). Right now, it's in a "paraelectric" state. Think of this like a crowd of people in a room who are all facing random directions. They are jiggling around, but there is no overall order.

As the room gets cooler (approaching a specific temperature called the transition temperature), something dramatic is about to happen: the crowd is about to suddenly all turn and face the same direction (becoming ferroelectric).

This paper explains two strange things that happen in the room just before everyone turns:

1. The "Central Peak": A mysterious, slow "sigh" or hum that appears in the data.
2. The "Acoustic Anomaly": The sound waves traveling through the room start to behave strangely, getting slower and getting absorbed more than usual.

The author, Akira Onuki, argues that these aren't caused by dirt, impurities, or defects in the crystal. Instead, they are caused by the electricity inside the crystal getting tangled up with the mechanical squeezing of the crystal structure.

---

The Main Characters: The "Electric" and the "Squeeze"

To understand the paper, we need two main characters:

1. The Polarization ($p$): Imagine this as the "mood" of the crowd. It represents the tiny electric arrows inside the atoms. Right now, they are chaotic.
2. The Elastic Displacement ($u$): Imagine this as the physical shape of the room. If the atoms push against each other, the room stretches or squishes.

The Connection (Electrostriction):
In most materials, electricity and shape don't talk much. But in BaTiO₃, they are best friends. This is called Electrostrictive (ES) coupling.

• The Metaphor: Imagine the "mood" (electricity) of the crowd is so strong that it physically forces the walls of the room to stretch or shrink. If the crowd gets excited (fluctuates), the room physically wobbles.

The Mystery: The "Central Peak"

For decades, scientists have looked at these crystals using light or neutrons (like taking a photo with a super-fast camera). They noticed a weird signal right in the middle of the data (at zero frequency). They called it the Central Peak.

• The Old Theory: Scientists used to think this peak was caused by "dust" (impurities) or tiny bubbles (domains) floating around.
• Onuki's Theory: He says, "No, it's not dust. It's the crowd itself."

The Analogy:
Imagine the crowd of people is trying to decide which way to face.

• Fast Fluctuations: Most people are just jiggling their heads quickly. This creates the "Acoustic Peak" (normal sound).
• Slow Fluctuations: But right before they all agree to face North, a few people start to hesitate. They wobble back and forth very slowly, taking a long time to decide.
• The Result: Because the "mood" (electricity) is so tightly linked to the "walls" (the crystal shape), these slow, hesitant wobbles make the walls of the room vibrate slowly too.

In the data, this slow vibration looks like a Central Peak. It's a "sigh" of the crystal as it struggles to decide its new state. The paper proves mathematically that this slow sigh is a natural result of the electricity squeezing the crystal, not a defect.

The Sound Waves: The "Acoustic Anomaly"

When you hit a crystal, sound waves travel through it. Usually, these waves travel at a steady speed. But near the transition temperature, the paper shows that:

1. The Sound Slows Down: The "stiffness" of the crystal changes. Because the electric mood is fluctuating so wildly, it makes the crystal feel "squishier" to the sound waves.
2. The Sound Gets Absorbed: The sound energy gets eaten up. Why? Because the sound wave is trying to push the crystal, but the crystal's electric mood is fighting back, turning that sound energy into heat (friction).

The Analogy:
Imagine running through a hallway.

• Normal Hallway: You run fast and smooth.
• The "Critical" Hallway: The walls are made of rubber bands that are attached to a nervous crowd. As you run, the walls stretch and snap back, but the crowd is hesitating. You lose your energy fighting the rubber bands. You slow down, and you get tired (absorbed) much faster.

The "Secret Sauce": The Time Delay

The paper introduces a concept called the Debye Relaxation Time ($\tau_D$).

• What it is: It's the "thinking time" of the crystal. How long does it take for the electric arrows to react to a squeeze?
• The Magic: The paper shows that the weird "Central Peak" and the "Sound Anomaly" are perfectly synchronized with this thinking time.
  • If you look at the crystal faster than its thinking time, it looks stiff.
  • If you look at it slower, you see the "Central Peak" (the slow sigh).

Why This Matters

For a long time, scientists thought these strange signals were just "noise" or "bad samples." This paper says: "No, this is the crystal talking."

It tells us that the Central Peak is a fundamental feature of how ferroelectric materials work. It's the sound of the material's internal electric forces and its physical shape dancing together right before a major phase change.

Summary in One Sentence

The paper explains that the strange, slow "hum" (Central Peak) and the weird slowing down of sound in Barium Titanate crystals aren't caused by dirt or defects, but are the natural result of the material's electric mood and physical shape getting stuck in a slow, hesitant dance right before they change state.

Technical Summary

Here is a detailed technical summary of the paper "Theory of central peak and acoustic anomaly in cubic BaTiO3 close to ferroelectric transition" by Akira Onuki.

1. Problem Statement

Ferroelectric transitions, particularly in materials like Barium Titanate (BaTiO$_3$), are characterized by critical fluctuations of polarization ($p$) and elastic strains ($e_{\alpha\beta}$). Two major phenomena have long been observed experimentally but lacked a unified theoretical explanation:

• The Central Peak: A narrow, low-frequency peak observed in dynamic scattering experiments (neutron and light scattering) near the transition temperature ($T_c$). Its origin has been debated, with previous attributions to thermal entropy diffusion (Rayleigh scattering), impurities, or small domains.
• Acoustic Anomalies: Significant changes in sound velocity and attenuation (critical sound absorption) as the material approaches the ferroelectric transition.

While piezoelectric systems (like KDP) show strong elastic anomalies due to linear coupling, cubic BaTiO$_3$ in its paraelectric phase is centrosymmetric. Here, the coupling between polarization and strain is electrostrictive (ES) (quadratic in polarization, $p^2 e$). The paper aims to derive a Ginzburg-Landau theory explaining how this ES coupling alone, without impurities or defects, generates the central peak and specific acoustic anomalies.

2. Methodology

The author employs a Ginzburg-Landau (GL) theory combined with linear response theory and mean-field approximation.

• Free Energy Functional: The system is described by a free energy functional $F$ containing:
  • Bulk polarization energy ($f_B$) with Landau expansion up to sixth order (crucial for the weakly first-order transition in BaTiO$_3$).
  • Electric field energy (accounting for depolarization effects).
  • Elastic energy ($f_{el}$).
  • Electrostrictive Coupling ($f_{int}$): A term bilinear in stress and quadratic in polarization ($-\sum \Pi^s_{\alpha\beta} \epsilon_{\alpha\beta}$), where the singular stress $\Pi^s_{\alpha\beta} \propto p_\alpha p_\beta$.
• Dynamics:
  • Polarization: Modeled via a relaxational Langevin equation (Model A dynamics) with a Debye relaxation time $\tau_D$.
  • Displacement: Modeled via the dynamic elasticity equation.
  • Decomposition: The elastic displacement $u$ is decomposed into an acoustic part ($w$) and an electrostrictive part ($m$). The ES part $m$ is slaved to the slow polarization fluctuations, while $w$ represents the fast acoustic modes.
• Correlation Functions: The paper calculates time-correlation functions for polarization and displacement in Fourier-Laplace space. It utilizes the Ginzburg criterion to justify the mean-field approximation, noting that for BaTiO$_3$, the critical region where fluctuations dominate is extremely narrow.

3. Key Contributions and Theoretical Framework

A. Origin of the Central Peak

The paper demonstrates that the central peak arises naturally from the slow relaxation of the electrostrictive stress tensor.

• Because the ES stress is proportional to $p^2$, and $p$ relaxes slowly near $T_c$ (with time scale $\tau_D$), the stress tensor also relaxes slowly.
• This slow relaxation creates a "quasi-static" component in the displacement correlation function.
• The Fourier-Laplace transform of the displacement time-correlation function yields a peak at $\omega = 0$ (the central peak) with a width proportional to the relaxation rate ($1/\tau_D$).
• Crucial Finding: This mechanism does not require impurities, vacancies, or pre-existing domains; it is an intrinsic property of the ES coupling in the paraelectric phase.

B. Frequency-Dependent Elastic Moduli

The author derives the complex, frequency-dependent elastic moduli $C^*_{ij}(\omega)$:
$$C^*_{ij}(\omega) = C_{ij} - M_{ij} I_e A^{-1/2} G_c(i\Omega)$$
Where:

• $A \propto (T/T_0 - 1)$ is the reduced temperature coefficient.
• $G_c(i\Omega)$ is a scaling function dependent on the scaled frequency $\Omega = \omega \tau_D / 2$.
• $M_{ij}$ are coefficients derived from electrostrictive constants ($g_{ij}$) and elastic constants ($c_{ij}$).

The theory predicts that:

• $C^*_{11}(\omega)$ and $C^*_{12}(\omega)$ exhibit strong critical anomalies (softening and frequency dependence).
• $C^*_{44}(\omega)$ remains nearly non-singular, consistent with experimental observations in cubic BaTiO$_3$.

C. Scaling and Universality

The dynamics are governed by a universal scaling function $K_c(s)$ (related to the Weber function) and its Laplace transform $F_c(i\Omega)$.

• The central peak height scales as $A^{-3/2}$.
• The peak width scales as $A$.
• The sound attenuation scales as $A^{-3/2}$ at low frequencies.

4. Key Results

1. Quantitative Agreement with Experiments:

   • Using experimental data from Kashida et al. (ultrasonic velocity) and Ko et al. (Brillouin scattering), the author determines the gradient coefficient $C \approx (1.5 \, \text{\AA})^2$ and the microscopic time $\tau_0 \approx 1.1 \times 10^{-2}$ ps.
   • The calculated Debye relaxation time $\tau_D \approx 0.62 / (T/T_0 - 1)$ ps matches the broadening of the central peak observed in light scattering experiments.
   • The theory successfully reproduces the magnitude of the elastic softening ($\Delta C_{11}/C_{11}$) and the specific frequency dependence of sound absorption.

2. Dynamic Structure Factor:

   • The calculated dynamic structure factor $I_D(q, \omega)$ shows two distinct features:
     • A Brillouin peak (acoustic mode) at finite frequency.
     • A Central peak at $\omega=0$.
   • The ratio of the central peak height to the acoustic peak height is derived. While the central peak grows as $T \to T_c$, the theory predicts it remains lower than the acoustic peak for the specific wave vectors used in neutron scattering experiments, consistent with observations.

3. Anisotropy:

   • The theory predicts an anisotropy factor $D(\hat{q})$ for the central peak intensity, which is maximum along the [001] direction and minimum along [111], matching the symmetry of the cubic crystal.

5. Significance and Conclusion

• Resolution of a Long-standing Mystery: The paper provides a definitive theoretical origin for the central peak in cubic ferroelectrics, attributing it to the intrinsic electrostrictive coupling between polarization and strain, rather than extrinsic defects.
• Unified Description: It offers a unified framework describing both the static elastic anomalies (softening) and dynamic anomalies (central peak, sound absorption) using a single set of parameters derived from the GL free energy.
• Predictive Power: The derived scaling functions ($K_c, F_c, G_c$) allow for precise predictions of frequency-dependent behavior near the transition, which can be tested against high-resolution scattering data.
• Implications for Thin Films: The theory is extended to include boundary conditions (fixed charge vs. fixed potential), offering insights into the dielectric response of ferroelectric thin films where surface effects are significant.

In summary, Akira Onuki's work establishes that the "central peak" in BaTiO$_3$ is a fundamental consequence of the slow relaxation of electrostrictive stress in the paraelectric phase, providing a robust Ginzburg-Landau explanation for critical acoustic and scattering anomalies in cubic ferroelectrics.