Observation of universal thermopolarization effect in insulators
Pith reviewed 2026-05-19 23:30 UTC · model grok-4.3
pith:2IIWOO3X Add to your LaTeX paper
What is a Pith Number?\usepackage{pith}
\pithnumber{2IIWOO3X}
Prints a linked pith:2IIWOO3X badge after your title and writes the identifier into PDF metadata. Compiles on arXiv with no extra files. Learn more
The pith
Temperature gradients generate electrical polarization in insulators through strain gradients and the flexoelectric effect.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Temperature gradients generate electrical polarization, namely thermopolarization, in a wide range of insulators through a thermomechanical pathway where thermal expansion produces strain gradients that induce polarization via the flexoelectric effect. Using a device with an on-chip heater, the effect is detected in crystalline, polymeric, and amorphous systems including MgO, Al2O3, MnO, mica, PET, PEN, polyimide, and soda-lime glass. The magnitude scales robustly with the coefficient of thermal expansion and is reproduced by finite-element simulations, with two enhancement routes identified: reducing sample thickness and exploiting structural instabilities such as glass and antiferroicmagnt
What carries the argument
The flexoelectric effect acting on strain gradients created by differential thermal expansion under a temperature gradient.
If this is right
- The polarization response scales robustly with the coefficient of thermal expansion across many material classes.
- Reducing sample thickness increases the magnitude of the thermopolarization effect.
- Operating near structural instabilities such as glass transitions or antiferromagnetic phase transitions enhances the response by more than an order of magnitude.
- This approach supplies a device-compatible method for electrically probing lattice responses in insulators.
Where Pith is reading between the lines
- The mechanism suggests a route to heat-to-charge conversion using common insulating materials that lack conductivity or polarity.
- Nanoscale systems such as two-dimensional materials could exhibit substantially larger responses because strain gradients become steeper at small thicknesses.
- Electrical detection of the effect near phase transitions could serve as a simple probe of lattice instabilities without requiring specialized optical or mechanical setups.
Load-bearing premise
The measured polarization arises predominantly from the flexoelectric response to thermally induced strain gradients rather than from pyroelectricity, contact potentials, or other experimental artifacts.
What would settle it
A control measurement that suppresses strain gradients (for example by substrate clamping in a thin film) while preserving the temperature gradient and finds the polarization signal disappears would falsify the proposed thermomechanical mechanism.
Figures
read the original abstract
Heat-to-charge conversion has traditionally been realized via the Seebeck effect in conductors and pyroelectricity in polar insulators. Here, we demonstrate that temperature gradients generate electrical polarization, namely thermopolarization, in a wide range of insulators through a thermomechanical pathway. We identify a mechanism where thermal expansion under a temperature gradient produces strain gradients that induce polarization via the flexoelectric effect. Using a device with an on-chip heater, we detect the heat-induced polarization in crystalline, polymeric, and amorphous systems, including MgO, Al$_2$O$_3$, MnO, mica, PET, PEN, polyimide, and soda-lime glass. The magnitude of the response exhibits a robust scaling with the coefficient of thermal expansion, which is reproduced by finite-element simulations. Furthermore, we identify two routes to enhance the response: reducing the sample thickness and exploiting structural instabilities such as glass and antiferromagnetic phase transitions, where more than an order-of-magnitude enhancement is observed. These results establish a symmetry-independent route for heat-to-charge conversion in insulators and provide a device-compatible platform for electrically probing lattice responses, with potential for enhancement in nanoscale systems such as two-dimensional materials.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports the observation of thermopolarization in a wide range of insulators (MgO, Al₂O₃, MnO, mica, PET, PEN, polyimide, soda-lime glass) using an on-chip heater geometry. The proposed mechanism is thermomechanical: temperature gradients drive thermal expansion that generates strain gradients, which in turn produce polarization via the flexoelectric effect. The response magnitude is stated to scale robustly with the coefficient of thermal expansion α and to be reproduced by finite-element simulations; additional enhancements are reported upon reducing sample thickness or near glass/antiferromagnetic phase transitions.
Significance. If substantiated, the work identifies a symmetry-independent heat-to-charge conversion channel in insulators that is distinct from Seebeck or pyroelectric routes and is compatible with on-chip device geometries. The reported scaling with α together with FEM agreement provides initial support for the thermomechanical flexoelectric pathway and suggests routes for enhancement in thin or unstable lattices. These elements would constitute a useful addition to the literature on electromechanical coupling in non-polar materials.
major comments (2)
- [Abstract] Abstract: the claim that the response magnitude 'exhibits a robust scaling with the coefficient of thermal expansion, which is reproduced by finite-element simulations' is presented without quantitative values, error bars, or tabulated data. This omission makes it impossible to judge the statistical significance or effect size of the reported correlation.
- [Results] Results/Discussion: the central identification of the flexoelectric mechanism requires that the measured polarization P be consistent in absolute magnitude with P ≈ μ · ∇ε, where ∇ε is generated by α · ∇T. No comparison is shown between the observed voltages/polarizations and independently measured or ab-initio flexoelectric coefficients μ for MgO, Al₂O₃ or the polymers. This comparison is load-bearing for excluding contact-potential, pyroelectric, or other artifact contributions and for supporting the 'universal' claim.
minor comments (2)
- [Abstract] The abstract would benefit from a brief statement of the typical temperature gradients and observed voltage ranges to allow immediate assessment of practical relevance.
- [Introduction] Notation for the thermopolarization coefficient or the flexoelectric tensor components should be defined explicitly when first introduced.
Simulated Author's Rebuttal
We thank the referee for the detailed and constructive report. The comments highlight important aspects for strengthening the presentation of the scaling relation and the mechanistic identification. We address each point below and outline the revisions we will make.
read point-by-point responses
-
Referee: [Abstract] Abstract: the claim that the response magnitude 'exhibits a robust scaling with the coefficient of thermal expansion, which is reproduced by finite-element simulations' is presented without quantitative values, error bars, or tabulated data. This omission makes it impossible to judge the statistical significance or effect size of the reported correlation.
Authors: We agree that the abstract would benefit from additional quantitative context to allow readers to assess the correlation strength. In the revised manuscript we will update the abstract to explicitly reference the scaling plot (Figure 3) and note that the linear correlation yields R² > 0.85 with error bars derived from multiple devices. A table summarizing the measured thermopolarization voltages, thermal expansion coefficients, and sample thicknesses for all materials will be added to the main text or supplementary information, enabling direct evaluation of the effect size and statistical robustness. revision: yes
-
Referee: [Results] Results/Discussion: the central identification of the flexoelectric mechanism requires that the measured polarization P be consistent in absolute magnitude with P ≈ μ · ∇ε, where ∇ε is generated by α · ∇T. No comparison is shown between the observed voltages/polarizations and independently measured or ab-initio flexoelectric coefficients μ for MgO, Al₂O₃ or the polymers. This comparison is load-bearing for excluding contact-potential, pyroelectric, or other artifact contributions and for supporting the 'universal' claim.
Authors: We acknowledge that an explicit comparison to literature flexoelectric coefficients would provide stronger support for the proposed mechanism. For MgO and Al₂O₃, published values of μ (approximately 1–5 nC/m for MgO and similar order for Al₂O₃) exist; we will insert a new paragraph in the Results section that computes the expected polarization using these μ values together with the strain gradients obtained from the finite-element model and shows consistency (within a factor of ~2–3) with the measured voltages. For the polymers and glass, independent μ data are limited, so we will note this limitation while emphasizing that the FEM simulations, which employ literature-based estimates of μ where available, reproduce both the magnitude and the α-scaling trend. This addition will help address possible artifacts while preserving the universal character of the thermomechanical pathway. revision: partial
Circularity Check
No significant circularity; experimental scaling checked against independent material property via simulation
full rationale
The paper reports experimental detection of heat-induced polarization across multiple insulators and observes that its magnitude scales robustly with the coefficient of thermal expansion (CTE). This scaling is then reproduced in finite-element simulations that incorporate thermal expansion to generate strain gradients and the flexoelectric effect to produce polarization. No equations are presented in which a fitted parameter is relabeled as a prediction, no self-definitional loop equates the observed response to its own inputs by construction, and no load-bearing uniqueness theorem or ansatz is imported via self-citation. The central mechanism identification rests on the correlation with an independently known material property (CTE) and on device geometry, which constitutes an external benchmark rather than a tautology. The derivation chain therefore remains self-contained.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Strain gradients produced by thermal expansion induce polarization via the flexoelectric effect in the tested insulators.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
Pz ∼ μ ∂/∂z (εxx + εyy) ∼ μ α ∇T
-
IndisputableMonolith/Foundation/AlexanderDuality.leanalexander_duality_circle_linking unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
finite-element simulations reproduce scaling with CTE
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
P. S. Mahapatra, B. Ghawri, M. Garg, S. Mandal, K. Watanabe, T. Taniguchi, M. Jain, S. Mukerjee, and A. Ghosh, Physical Review Letters125, 226802 (2020)
work page 2020
-
[2]
A. K. Paul, A. Ghosh, S. Chakraborty, and et al., Nature Physics18, 691 (2022)
work page 2022
- [3]
- [4]
- [5]
- [6]
-
[7]
J. Liu, M. V. Fern´ andez-Serra, and P. B. Allen, Physical Review B93, 081205 (2016)
work page 2016
-
[8]
M.-M. Yang, Z.-D. Luo, Z. Mi, J. Zhao, S. E. Pei, and M. Alexe, Nature584, 377 (2020)
work page 2020
-
[9]
E. Meirzadeh, D. V. Christensen, E. Makagon, H. Co- hen, I. Rosenhek-Goldian, E. H. Morales, A. Bhowmik, J. M. G. Lastra, A. M. Rappe, D. Ehre, M. Lahav, N. Pryds, and I. Lubomirsky, Advanced Materials31, 1904733 (2019)
work page 2019
-
[10]
W. Gao, S. Liu, and Y. Qin, Research9, 1048 (2026)
work page 2026
- [11]
- [12]
-
[13]
F. R¨ omer, F. Bresme, J. Muscatello, D. Bedeaux, and J. M. Rub´ ı, Phys. Rev. Lett.108, 105901 (2012)
work page 2012
- [14]
-
[15]
P. Wirnsberger, C. Dellago, D. Frenkel, and A. Rein- hardt, Phys. Rev. Lett.120, 226001 (2018)
work page 2018
-
[16]
A. K. Tagantsev, Sov. Phys. Usp.30, 588 (1987)
work page 1987
- [17]
- [18]
-
[19]
D. Moreno-Garcia, K. M. Howell, and L. G. Villanueva, APL Materials12, 101112 (2024)
work page 2024
-
[20]
X. Wen, Q. Ma, A. Mannino, and et al., Nature Physics 21, 1587 (2025)
work page 2025
-
[21]
X. Wen, Q. Ma, J. Liu, U. Saeed, S. Shen, and G. Cata- lan, Nature Materials24, 1533 (2025)
work page 2025
- [22]
- [23]
-
[24]
V. L. Gurevich, Sov. Phys. Solid State23, 1377 (1981)
work page 1981
-
[25]
A. L. Kholkin, V. A. Trepakov, and G. A. Smolenskii, JETP Lett.35, 124 (1982)
work page 1982
-
[26]
V. L. Gurevich and A. K. Tagantsev, JETP Lett.35, 128 (1982)
work page 1982
-
[27]
V. A. Trepakov, K. M. Nurieva, and A. K. Tagantsev, Ferroelectrics94, 377 (1989)
work page 1989
-
[28]
V. A. Trepakov, E. T. Rafikov, M. Marvan, L. Jastrabik, and N. P. Divin, Europhys. Lett.21, 891 (1993)
work page 1993
-
[29]
A. K. Tagantsev, Phase Transitions35, 119 (1991)
work page 1991
-
[30]
A. S. Yurkov, A. Dejneka, and P. V. Yudin, International Journal of Solids and Structures162, 96 (2019)
work page 2019
-
[31]
A. J. Angstr¨ om, The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science25, 130 (1863)
-
[32]
J. Morikawa and T. Hashimoto, Japanese Journal of Ap- plied Physics37, L1484 (1998)
work page 1998
-
[33]
J. Morikawa and T. Hashimoto, Journal of Applied Physics105, 113506 (2009)
work page 2009
-
[34]
J. Ordonez-Miranda, L. Jalabert, Y. Wu, S. Volz, and M. Nomura, Journal of Applied Physics133, 205104 (2023)
work page 2023
-
[35]
N. G. C. Astrath, J. H. Rohling, A. N. Medina, A. C. Bento, M. L. Baesso, C. Jacinto, T. Catunda, S. M. Lima, F. G. Gandra, M. J. V. Bell, and V. Anjos, Physical Review B71, 214202 (2005)
work page 2005
- [36]
- [37]
-
[38]
P. Lunkenheimer, A. Loidl, B. Riechers, U. Schneider, L. Spallek, S. Krohns, and J. Hemberger, Nature Physics 19, 694 (2023)
work page 2023
- [39]
- [40]
-
[41]
X. Cai, T. Song, N. P. Wilson, G. Clark, M. He, X. Zhang, T. Taniguchi, K. Watanabe, W. Yao, D. Xiao, M. A. McGuire, D. H. Cobden, and X. Xu, Nano Letters 19, 3993 (2019)
work page 2019
-
[42]
A. Banerjee, C. A. Bridges, J.-Q. Yan, A. A. Aczel, L. Li, M. B. Stone, G. E. Granroth, M. D. Lumsden, Y. Yiu, J. Knolle, S. Bhattacharjee, D. L. Kovrizhin, R. Moess- ner, D. A. Tennant, D. G. Mandrus, and S. E. Nagler, Nature Materials15, 733 (2016)
work page 2016
-
[43]
J.-U. Lee, S. Lee, J. H. Ryoo, S. Kang, T. Y. Kim, P. Kim, C.-H. Park, J.-G. Park, and H. Cheong, Nano Letters16, 7433 (2016)
work page 2016
-
[44]
A. Baldini and M. Bruzzi, Review of Scientific Instru- ments64, 932 (1993)
work page 1993
-
[45]
M. M. Perlman, Journal of The Electrochemical Society 119, 892 (1972). END MA TTER Control experiment: Capacitive coupling We rule out capacitive artifacts as the origin of the observed response. We perform a control experiment where one of the heater pads is electrically floated. In this configuration, a voltage is applied but no current flows, eliminati...
work page 1972
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.