Piezoelectric Strain Coefficients (d33)

• Piezoelectric strain coefficient d33 is a measure of electric polarization generated per unit mechanical stress along the material’s principal axis, combining intrinsic lattice deformations and extrinsic domain effects.
• It is determined using methods such as the Berlincourt technique, impedance-based resonance, and advanced piezoelectric force microscopy, ensuring accurate characterization from bulk crystals to nanoscale films.
• Optimization of d33 through compositional tuning, domain engineering, and strain management enhances performance in practical applications like sensors, actuators, and energy harvesting devices.

The piezoelectric strain coefficient $d_{33}$ quantifies the induced electric charge per unit mechanical stress, or equivalently, the mechanical strain per unit electric field, measured along a material’s principal (usually poling) axis. It is central to electromechanical energy conversion in piezoelectric materials and serves as a key figure of merit for devices such as actuators, sensors, and energy harvesters. $d_{33}$ encompasses contributions from intrinsic lattice deformation as well as extrinsic domain and phase-change effects, and its value depends critically on crystallography, material composition, microstructure, sample processing, and measurement methodology.

1. Definition and Measurement Principles

The definition of $d_{33}$ corresponds to the response in the so-called longitudinal mode, typically given by

$$d_{33} = \left.\frac{\partial S_3}{\partial E_3}\right|_{E=0} = \frac{Q}{F}$$

where $S_3$ is the strain along the polar (3) axis, $E_3$ the electric field along that axis, $Q$ the charge generated (direct effect), and $F$ the applied force. In converse measurements, $d_{33}$ also describes the induced strain per applied field.

Various experimental approaches are used to quantify $d_{33}$:

• Berlincourt method (quasistatic force application and charge detection) for bulk and single crystals (Guennou et al., 2012).
• Impedance-based resonance methods, employing standards such as EN 50324-1:2002/IEC 60483:1976, relating resonance frequency data and sample geometry to extract $d_{33}$0, dielectric permittivity $d_{33}$1, elastic compliance $d_{33}$2, and electromechanical coupling factor $d_{33}$3 via

$d_{33}$4

(Kumar et al., 2022).

• Piezoelectric force microscopy (PFM) and advanced variants such as Direct Piezoelectric Force Microscopy (DPFM; (Gomez et al., 2016)) or optimized PFM with mitigation of parasitic contributions (Singh et al., 6 May 2025), for nanoscale or thin-film systems.
• For composites or polymer-based systems, local displacement measurements via AFM are common, extracting $d_{33}$5 from the slope of displacement versus bias (Das et al., 2019).

2. Material Chemistry, Structure, and d₃₃ Enhancement

$d_{33}$6 is highly sensitive to the chemical composition, presence of dopants or alloying elements, phase symmetries, and microstructural features:

• Single-crystal and engineered domain states: For Mn-doped NBT-6%BT single crystals, a single domain configuration achieved via [001] poling yields $d_{33}$7 pC/N, while multidomain states formed by poling along [011] or [111] drive $d_{33}$8 down to 225 or 130 pC/N, reflecting the dilution of piezoelectric anisotropy through domain averaging (Guennou et al., 2012).
• Phase coexistence and morphotropic phase boundaries (MPB): The coexistence of multiple ferroelectric phases near MPB, as in PMN-PT, PZT, and Ru-doped BaTiO$d_{33}$9, provides easy polarization rotation pathways and enhances $d_{33}$0 (e.g., up to 530 pC/N in PMN-BT-PT at room temperature (Huang et al., 2013), 269 pC/N in Ru:BaTiO$d_{33}$1 (Brajesh et al., 2021), or 200 pC/N in PMN-0.3PT (Kumar et al., 2022)).
• Alloying strategies: In wurtzite-structure AlN, alloying with transition metals—especially Mo, Sc, Cr—softens the elastic modulus and increases the non-clamped ionic contribution to the piezoelectric stress constant $d_{33}$2, leading to substantial $d_{33}$3 enhancements up to an order of magnitude ($d_{33}$4 up to 58.6 pC/N for optimized Mo$d_{33}$5Al$d_{33}$6N (Zha et al., 2022)).
• Compositional tuning of 2D materials: Layered or monolayer systems such as $d_{33}$7-AsP or g-C$d_{33}$8N$d_{33}$9 can exhibit exceptionally large, strain-enhanced $$d_{33} = \left.\frac{\partial S_3}{\partial E_3}\right|_{E=0} = \frac{Q}{F}$$0/$$d_{33} = \left.\frac{\partial S_3}{\partial E_3}\right|_{E=0} = \frac{Q}{F}$$1 due to highly anisotropic or puckered bonding, and strain engineering can boost $$d_{33} = \left.\frac{\partial S_3}{\partial E_3}\right|_{E=0} = \frac{Q}{F}$$2 by over 300% (Guo et al., 2019, Guo et al., 2020).
• Thin films and substrate constraint: In PMN-PT thin films, mechanical clamping to a substrate severely inhibits field-induced phase transitions, quenching $$d_{33} = \left.\frac{\partial S_3}{\partial E_3}\right|_{E=0} = \frac{Q}{F}$$3 (<100 pm/V), whereas release into a membrane allows the full transition and $$d_{33} = \left.\frac{\partial S_3}{\partial E_3}\right|_{E=0} = \frac{Q}{F}$$4 pm/V—matching bulk values (Brewer et al., 2021).

Material System	Peak $$d_{33} = \left.\frac{\partial S_3}{\partial E_3}\right|_{E=0} = \frac{Q}{F}$$5 (pC/N or pm/V)	Notes
Mn-doped NBT-6%BT [001], single	570	Single-domain, [001] poled (Guennou et al., 2012)
PMN-BT-PT, $$d_{33} = \left.\frac{\partial S_3}{\partial E_3}\right|_{E=0} = \frac{Q}{F}$$6	530	Large, temperature-stable (Huang et al., 2013)
Sr-doped PMS-PZT	281	Near MPB, $$d_{33} = \left.\frac{\partial S_3}{\partial E_3}\right|_{E=0} = \frac{Q}{F}$$7, $$d_{33} = \left.\frac{\partial S_3}{\partial E_3}\right|_{E=0} = \frac{Q}{F}$$8 (Brajesh et al., 2013)
Ru-doped BaTiO$$d_{33} = \left.\frac{\partial S_3}{\partial E_3}\right|_{E=0} = \frac{Q}{F}$$9, 2% Ru	269	Phase coexistence at RT (Brajesh et al., 2021)
PMN-0.3PT ceramic	200 ± 21	$S_3$0, $S_3$1 F/m (Kumar et al., 2022)
Polycrystalline PZT (MPB)	~220–230	With strong extrinsic contributions
Oriented PZT (MPB)	46.4 ± 4.4	Intrinsic, c-axis oriented (Fu et al., 2024)
Mo$S_3$2Al$S_3$3N (metastable config.)	58.6	12.3× pure AlN (Zha et al., 2022)

3. Intrinsic and Extrinsic Contributions

$S_3$4 is typically a sum of lattice-level (intrinsic) and microstructural (extrinsic) effects:

• Intrinsic component: Due to atomic displacement within the unit cell under electric field or stress. First-principles-derived $S_3$5 for aligned tetragonal PZT near MPB is 50–55 pC/N (Fu et al., 2024).
• Extrinsic contributions: From irreversible (non-180$S_3$6) domain wall motion, field-induced phase transformations (especially in MPB compositions), and intergranular interactions. In PZT ceramics, $S_3$7 is amplified %%%%46$E_3$47%%%% over the intrinsic value due to these extrinsic mechanisms.
• Direct measurement distinguishing intrinsic/extrinsic: DPFM can quantitatively map $E_3$0 at the nanoscale, providing values of $E_3$1 pC/N for PPLN and $E_3$2 pC/N for BiFeO$E_3$3 consistent with literature and substantially isolating intrinsic effects (Gomez et al., 2016).
• Mechanisms at MPB: In systems like PbTiO$E_3$4-BiScO$E_3$5, the dominant mechanism for high $E_3$6 is a field-induced inter-ferroelectric transformation rather than enhanced domain wall mobility (V. et al., 2015); i.e., large $E_3$7 does not necessarily correlate with higher lattice strain or domain switching but with phase conversion.

4. Dependence on Microstructure, Orientation, and Strain

Crystal orientation, microstructure, and engineered strain play decisive roles:

• Domain engineering: Achieving a single domain state (e.g., via [001] poling in NBT-6%BT) maximizes intrinsic $E_3$8. Multidomain configurations dilute the response proportionally to tensor rotation projections (Guennou et al., 2012).
• Strain and domain orientation in thin films: In PZT films, compressive misfit strain induces 30$E_3$9 domain alignment, enhancing out-of-plane $Q$0 (to values $Q$1593 pm/V); La$Q$2-doping reduces misfit and switches domain alignment to the $Q$3-axis, nearly doubling the in-plane $Q$4 (A et al., 2023).
• Clamping vs. membrane in PMN-PT thin films: Only free-standing configurations permit full polarization rotation (R$Q$5T), yielding giant $Q$6 ($Q$71000 pm/V), while substrate clamping restricts transformation to (R$Q$8M) with $Q$9 pm/V (Brewer et al., 2021).
• Strain engineering in 2D materials: Biaxial or uniaxial strain in materials such as g-C$F$0N$F$1 and $F$2-AsP dramatically increases $F$3 or $F$4 by reducing the elastic modulus denominator and enhancing ionic polarization response; similar mechanisms may be leveraged for out-of-plane $F$5 (Guo et al., 2020, Guo et al., 2019).

5. Measurement Methodology, Calibration, Artifacts, and Standardization

Reliable extraction of $F$6 from experiments requires careful attention to parasitic and extrinsic effects:

• Mitigation of electrostatic and nonlocal contributions: In PFM, the combination of optimal optical lever sensitivity (by adjusting laser spot position), application of KPFM-determined DC offsets, and edge scanning minimizes both local and nonlocal electrostatic contributions. The effective $F$7 is then confidently extracted from the linear response to AC bias (Singh et al., 6 May 2025).
• Resonance vs. direct methods: The impedance resonance approach yields effective $F$8 values that include both intrinsic and extrinsic contributions, while DPFM or optimized PFM (with artifact subtraction) can isolate more intrinsic values (Brajesh et al., 2013, Gomez et al., 2016).
• Standardization and comparison: Equations such as $F$9 and figure-of-merit metrics $d_{33}$0 enable benchmarking for device applications and facilitate cross-comparison (Kumar et al., 2022).

6. Composition, Anisotropy, and Design Implications

High-throughput and ab initio studies of perovskite and related structures underscore complex dependencies:

• Tetragonality and suppression/saturation: Excessive tetragonality ($d_{33}$1) leads to polarization saturation, with $d_{33}$2 becoming suppressed despite high $d_{33}$3 ratios; such anisotropy metrics may not correlate with practical performance (Jia et al., 2023).
• Design trade-offs: Targeting high $d_{33}$4, substantial polarization, and soft elastic response along the field axis is favored, but large anisotropy ratios can reflect vanishing $d_{33}$5 rather than exceptional $d_{33}$6.
• Distribution/sublayer engineering: In AlN alloys, distributing Mo (or other TM) over multiple Al sublayers and avoiding dimer formation along $d_{33}$7 maximizes $d_{33}$8 (Zha et al., 2022).

7. Applications and Energy Conversion Performance

$d_{33}$9 is a principal parameter in electromechanical conversion efficiency, underpinning:

• Actuators, sensors, and ultrasonic devices, requiring both high $d_{33}$0 and temperature-stable response (e.g., $d_{33}$1 pC/N, $d_{33}$2 pC/N with near-constant $d_{33}$3 over 185–360 K in PMN-BT-PT (Huang et al., 2013)).
• Energy harvesters, where $d_{33}$4 directly relates to generated open-circuit voltage and short-circuit current under dynamic mechanical loads; e.g., PMN-0.3PT delivers $d_{33}$5 pC/N with maximum current density 95 nA/cm$d_{33}$6 and $d_{33}$7 V/cm under 3 N load (Kumar et al., 2022).
• Polymer–nanoparticle composites: Through nanoparticle inclusion (PMNPT into SU8), $d_{33}$8 pm/V is achieved, with compatibility for standard microfabrication and efficient RF acoustic transduction up to 1.5 GHz (Das et al., 2019).

In summary, $d_{33}$9 serves as a fundamental descriptor of a piezoelectric material’s direct electromechanical conversion capability, with its magnitude and functional behavior determined by intrinsic lattice effects, extrinsic phase and domain activity, microstructure, and precise measurement methodologies. Careful engineering of all these facets enables optimization for a broad spectrum of device applications ranging from high-precision actuators to energy harvesters and integrated nanoscale electronics.