REVIEW 3 major objections 5 minor 43 references
Barkhausen noise in the columnar hexagonal organic ferroelectric BTA
T0 review · 3 major / 5 minor · reviewed 2026-08-11 · deepseek-v4-flash
Pith's one-line read The paper claims that BTA polarization reversal is a critical avalanche process below about 175 K, with power-law distributed switching events too small to be seen as Barkhausen noise.
desk verdict First BTA Barkhausen kMC study with an honest experimental null result, but the claimed mean-field SOC exponents are compared against a benchmark the paper's own introduction contradicts. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is a kinetic Monte Carlo model of BTA in which each amide group is a flipable dipole on a hexagonal columnar lattice, with electrostatic dipole-dipole interactions computed inside a cutoff sphere and a reaction-field approximation for longer-range contributions. Structural disorder enters as subcolumns of randomly drawn lengths, positional shifts, rotations, and occasional chirality changes between subcolumns; a dipole flips with rate $\nu_0\exp(-\Delta U/k_BT)$ when flipping is uphill and $\nu_0$ otherwise, which lets cascades of flips become avalanches. The argument is carried by comparing the fitted exponents of the event-size and event-energy distributions against the mean-field plasticity predictions, and by the observation that increasing column height, but not lateral box width, enlarges the largest events. The machinery produces the claimed temperature threshold near 175 K, the disorder dependence analogous to a random-field Ising model, and the tiny event sizes that explain the null experimental result.
What would settle it
Run the same simulated event lists through an automated maximum-likelihood power-law fit with threshold as a free parameter: if the fitted exponents drift with threshold or depart from $\tau\approx1.5$ and $\epsilon\approx1.33$ across the 100–175 K range, the self-organized criticality claim is not supported. Experimentally, a device with a much smaller electrode area that reduces the required simultaneous-switch count to a few hundred dipoles should resolve individual BTA avalanches if they exist.
Extended reading notes
Core claim
The authors report that field-driven polarization reversal in the supramolecular ferroelectric BTA is a critical crackling process at low temperature: dipole flips organize into avalanches whose size and energy distributions follow power laws with exponents close to the mean-field plasticity values (about 1.5 for sizes and about 1.33 for energies), and these exponents are stable against disorder strength and modest parameter changes below roughly 175 K. The avalanches propagate essentially one-dimensionally along the supramolecular columns, so lateral coupling between columns is weak. At higher temperatures the power-law exponents increase, indicating creep, and the simulated events are so small that a purpose-built high-sensitivity electrical setup should not detect them; no Barkhausen noise was observed, which the paper takes as consistency with the model rather than as evidence against switching avalanches.
Load-bearing premise
The claim of self-organized criticality rests on the assumption that the fitted power-law exponents reflect the simulated switching physics rather than the manual choices of event threshold and fit range.
Editorial extensions
If this is right
- Below about 175 K, polarization switching in BTA should be scale-invariant: the same avalanche statistics hold across decades of event size, up to the system-size cutoff.
- Raising temperature or lowering sweep frequency should move the system from critical avalanches into creep, visible as larger fitted power-law exponents under the same analysis.
- Barkhausen noise in BTA should be undetectable in conventional and even moderately optimized setups, because single-column avalanches involve hundreds of dipoles rather than the roughly 600,000 simultaneously switching regions needed for a measurable current.
- Switching in the hexagonal columnar phase is quasi-one-dimensional: lateral column coupling is weak, so increasing the lateral box size should not change avalanche sizes, while column height sets the upper cutoff below roughly 300 molecules.
Reading between the lines
- If the low-temperature exponents are truly disorder-independent, one testable extension is to seek the same critical exponents in other columnar supramolecular ferroelectrics with similar hydrogen-bonded stacks, since the mechanism would be morphological rather than chemical.
- The null experimental result could be turned positive using nanoscale electrodes or single-column devices, where the number of switchable dipoles is small enough that individual avalanches would exceed the measurement noise; a power-law distribution of switching currents there would confirm self-organized criticality in real BTA.
- The paper's logic implies that the absence of Barkhausen noise in a ferroelectric is not evidence against jerky switching; it can simply mean the switching entities are too small, so interpreting null results requires a model of event sizes such as the one provided here.
- The claimed temperature crossover near 175 K and its frequency dependence could be probed by measuring coercive-field statistics across the hysteresis loop at different temperatures, since the paper argues the exponents are anti-correlated with coercivity.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper reports kinetic Monte Carlo simulations of polarization reversal in the columnar organic ferroelectric BTA, from which avalanche size and energy distributions are extracted and fitted to power laws. The central claim is that for temperatures below 175 K and sufficiently strong structural disorder the system exhibits self-organized criticality, with exponents matching mean-field plasticity values, while at higher temperatures a creep regime is entered. The authors also report an experimental search for Barkhausen noise in BTA that found none, consistent with their simulated event sizes being well below the measurement noise floor.
Significance. If the self-organized criticality claim were supportable, it would be an interesting advance for organic ferroelectrics, where crackling noise has rarely been studied, and the combination of a detailed kMC model with a null experimental result is a useful attempt at quantitative consistency. The manuscript is careful in describing the simulation protocol and the experimental setup, and it makes explicit, falsifiable statements about the temperature and disorder dependence of the exponents. However, the central claim is currently undermined by an inconsistency between the mean-field benchmark used for comparison and the field-integrated nature of the simulated events, as well as by the lack of robustness checks for the analysis choices that determine the fitted exponents.
major comments (3)
- [Introduction and Figure 6 / SI S5.3] The mean-field benchmark used for the horizontal lines in Fig. 6 and for the values quoted in SI S5.3 is the constant-stress set (τ=1.5, ε=4/3), but the events are collected over a full hysteresis loop (Section 3, Fig. 2). In the Introduction the authors explicitly state that for events measured over a whole (stress-integrated) loop the mean-field values are τ=2 and ε=5/3, citing Ref. [16]. Since the simulated events come from complete field ramps, the correct benchmark under the paper's own framework is the loop-integrated set. The reported agreement with τ≈1.5 and ε≈1.33 therefore does not support the claim of self-organized criticality; it actually indicates a disagreement with the loop-integrated prediction. Please redo the comparison either against τ=2, ε=5/3, or by restricting the event analysis to the coercive-field region and stating clearly which benchmark is used.
- [SI S2 and Section 3 (threshold and fitting procedure)] The extracted exponents, which form the basis of the SOC claim, depend on several manual analysis choices that are not accompanied by robustness checks. The event threshold is set 'just right of the thermal noise peak' (SI S2), the power-law fit range is determined 'by visualizing the event data' (SI S2), and data sets with different thresholds are merged by creating 'weight events' (Section 3). No error bars, bootstrap estimates, or sensitivity analyses are reported for the fitted exponents. It is therefore unclear whether the apparent low-temperature plateau in Fig. 6a is a genuine physical feature or an artifact of these choices. Please provide a systematic variation of the threshold level, the fit range, and the merging procedure, and report the resulting spread in the exponents.
- [SI S5 (disorder dependence)] The claim that the exponents are invariant to disorder is not fully supported by the presented data. Table S5.1 shows that the size exponent τ varies from 1.30 to 1.71 across chirality settings, and Fig. S5.2b shows a clear increase of the exponents for small subcolumn lengths (N ≤ 5) and a breakdown of the power-law form for N = 100. These variations are acknowledged qualitatively, but they are not quantified against the proposed universal regime. If the SOC regime is defined only for an intermediate range of disorder, the boundaries of that range should be stated and the universality claim restricted accordingly.
minor comments (5)
- [Section 4 and Figure 6b] The temperature unit is inconsistent: Fig. 6b and the surrounding text use '160 °C' for the sweeping-frequency study, whereas the text elsewhere uses kelvin (e.g., 'below 175 K' and 'T = 160 K' in SI S4). Please correct the unit throughout.
- [SI S2] The text refers to 'Figure 1d' when describing the merging of data sets; the corresponding figure in the main text is Figure 2d. Please correct the reference.
- [Section 4] There is a typo: 'the size S is defined as he highest number of flipped dipoles per time step' should read 'the highest number'. Please also clarify once whether S is the maximum instantaneous flip rate or the total number of flips, since the summed size SΣ is later defined separately.
- [SI S5.3 caption] The caption says 'the extracted power-law exponents stay at the mean field predictions of τ = 1.5 and ε = 1.33'; the value 1.33 is the decimal approximation of 4/3, but the text elsewhere uses the fraction. Please use one consistent representation.
- [Section 5] The estimate in SI S7 rests on the assumption that the experimental sample can be decomposed into independent ~700-molecule columns and that an event size of 200 dipoles is representative. It would be helpful to state explicitly that this is an order-of-magnitude estimate and to discuss how the conclusion would change if larger lateral avalanches (e.g., spanning several columns) were possible.
Circularity Check
No significant circularity: the SOC exponents are emergent outputs of the kMC model and are tested against external mean-field benchmarks, not fitted inputs.
full rationale
The central claim (self-organized criticality below ~175 K) is an emergent result of the kinetic Monte Carlo simulations. The power-law exponents are extracted from simulated avalanche-size and energy distributions and compared with externally published mean-field plasticity values (refs 12-16), not with any quantity the model was tuned to reproduce. The model parameters are inherited from prior work by the same group (refs 25-27), but those self-citations determine the Hamiltonian and switching rates, not the exponents; the simulations are not fitted to tau=1.5 or epsilon=4/3. The disorder-invariance tests in SI S5 (e.g., Fig. S5.3 and S5.4) explicitly show that the exponents remain near the mean-field values while the subcolumn-length variance is varied, which is the opposite of fitting the claim. The experimental null result is used only as a consistency check (SI S7), where simulated event sizes are compared with the setup noise floor; the absence of Barkhausen noise is not used to define or fit the exponents. The manual threshold and visually chosen fit range described in SI S2 are analysis choices that affect robustness, but they do not make the extracted exponents equal to the inputs by construction. One non-circular concern worth noting: the paper compares full-loop avalanche exponents to the constant-stress mean-field values (tau=1.5, epsilon=4/3) even though its own introduction states that the stress-integrated loop values are tau=2 and epsilon=5/3; this is a benchmark-consistency issue that should be addressed, but it is not a circular reduction of the derivation to its inputs.
Assumptions & free parameters
free parameters (4)
- Event threshold level =
~15 dipoles per time step
- Power-law fit range =
e.g., S = 40 to 250 (Figure S5.3)
- Weight-event correction for merged thresholds =
Not specified
- Subcolumn length N and variance sigma_N =
N=10, sigma_N=1 (typical)
assumptions (4)
- domain assumption The kMC model with fixed dipole positions and only electrostatic interactions (permanent and induced dipoles, reaction field, cutoff r_c=30) captures the essential switching dynamics of BTA.
- domain assumption The disorder in real BTA films is represented by subcolumns with Gaussian length distribution (mean N, variance sigma_N), positional shifts, random rotations, and helicity changes.
- domain assumption The mean-field plasticity / random-field Ising universality class applies to field-driven switching in BTA, so comparison of exponents to mean-field predictions is a valid test of SOC.
- domain assumption The reaction-field cutoff at r_c=30 and the precomputed nearest-neighbor flipping-energy differences remain accurate throughout an avalanche.
Cite this review
Pith. "Pith review of Barkhausen noise in the columnar hexagonal organic ferroelectric BTA." pith.science (2026). https://pith.science/paper/DX5QONYJ
@misc{pith2026241212666,
author = {Pith},
title = {Pith review of: Barkhausen noise in the columnar hexagonal organic ferroelectric BTA},
year = {2026},
howpublished = {\url{https://pith.science/paper/DX5QONYJ}},
note = {Machine review of arXiv:2412.12666}
}
read the original abstract
Upon a polarization reversal within a ferroelectric material, one stable state changes into another which is typically described by a progression of switching events of smaller fractions of the material. These events give rise to crackling or Barkhausen noise and follow a characteristic distribution in their sizes. Barkhausen noise has been studied to better understand the switching processes of ferroelectrics and has been applied for inorganic ferroelectric materials and perovskites. In this work, we present results from kinetic Monte Carlo simulations investigating the switching process of the small organic molecular ferroelectric benzene-1,3,5-tricarboxamides (BTAs). For temperatures below 175 K and sufficiently strong structural disorder, the system exhibits self-organized critical behavior; for higher temperatures, a creep regime is entered. Our extracted power-law exponents are smaller than those typically measured in inorganic crystals and ceramics which indicates that in the more disordered material BTA larger spanning avalanches are possible. The system was experimentally investigated with a high-sensitivity setup. No Barkhausen noise was observed which is consistent with the simulated event sizes, lying several orders beneath the noise threshold of the experimental setup. This finding corroborates the notion that switching in BTA progresses along the 1D columns in the hexagonal liquid crystal lattice, with little coupling between the columns that could give rise to larger lateral avalanches.
Reference graph
Works this paper leans on
-
[16]
Durin, S
G. Durin, S. Zapperi, J. Stat. Mech. 2006, 2006, P01002
2006
-
[1]
V. Ya. Shur, E. L. Rumyantsev, D. V. Pelegov, V. L. Kozhevnikov, E. V. Nikolaeva, E. L. Shishkin, A. P. Chernykh, R. K. Ivanov, Ferroelectrics 2002, 267, 347
2002
-
[2]
I. S. Baturin, M. V. Konev, A. R. Akhmatkhanov, A. I. Lobov, V. Ya. Shur, Ferroelectrics 2008, 374, 136
2008
-
[3]
V. Ya. Shur, A. R. Akhmatkhanov, I. S. Baturin, E. V. Shishkina, Journal of Applied Physics 2012, 111, 014101
2012
-
[5]
Flannigan, C
C. Flannigan, C. D. Tan, J. F. Scott, J. Phys.: Condens. Matter 2019, 32, 055403
2019
-
[7]
Casals, G
B. Casals, G. F. Nataf, E. K. H. Salje, Nat Commun 2021, 12, 345
2021
-
[8]
M. Mai, H. Kliem, Journal of Applied Physics 2013, 114, 224104
2013
-
[10]
J. P. Sethna, K. A. Dahmen, C. R. Myers, Nature 2001, 410, 242
2001
Show all 43 references
-
[11]
H. J. Williams, W. Shockley, Phys. Rev. 1949, 75, 178
1949
-
[12]
K. A. Dahmen, Y. Ben-Zion, J. T. Uhl, Phys. Rev. Lett. 2009, 102, 175501
2009
-
[13]
E. K. H. Salje, K. A. Dahmen, Annu. Rev. Condens. Matter Phys. 2014, 5, 233
2014
-
[14]
K. A. Dahmen, in Avalanches in Functional Materials and Geophysics (Eds.: E. K. H. Salje, A. Saxena, A. Planes), Springer International Publishing, Cham, 2017, pp. 19–30
2017
-
[15]
Zaiser, B
M. Zaiser, B. Marmo, P. Moretti, in Proceedings of International Conference on Statistical Mechanics of Plasticity and Related Instabilities — PoS(SMPRI2005), Sissa Medialab, Indian Institute of Science, Bangalore, India, 2006, p. 053
2006
-
[18]
A. V. Gorbunov, Ferroelectric Switching Phenomena in Organic Dielectrics and Semiconductors, Technische Universiteit Eindhoven, 2018
2018
-
[19]
Urbanaviciute, Multifunctional Supramolecular Organic Ferroelectrics, Linköping University Electronic Press, 2019
I. Urbanaviciute, Multifunctional Supramolecular Organic Ferroelectrics, Linköping University Electronic Press, 2019
2019
-
[20]
Urbanaviciute, X
I. Urbanaviciute, X. Meng, M. Biler, Y. Wei, T. D. Cornelissen, S. Bhattacharjee, M. Linares, M. Kemerink, Mater. Horiz. 2019, 6, 1688
2019
-
[21]
Cornelissen, Switching Kinetics and Charge Transport in Organic Ferroelectrics, Linköping University Electronic Press, Linköping, 2020
T. Cornelissen, Switching Kinetics and Charge Transport in Organic Ferroelectrics, Linköping University Electronic Press, Linköping, 2020
2020
-
[22]
A. V. Gorbunov, T. Putzeys, I. Urbanavičiūtė, R. a. J. Janssen, M. Wübbenhorst, R. P. Sijbesma, M. Kemerink, Phys. Chem. Chem. Phys. 2016, 18, 23663
2016
-
[23]
C. F. C. Fitié, W. S. C. Roelofs, M. Kemerink, R. P. Sijbesma, J. Am. Chem. Soc. 2010, 132, 6892
2010
-
[24]
C. F. C. Fitié, W. S. C. Roelofs, P. C. M. M. Magusin, M. Wübbenhorst, M. Kemerink, R. P. Sijbesma, J. Phys. Chem. B 2012, 116, 3928
2012
-
[25]
Urbanavičiūtė, T
I. Urbanavičiūtė, T. D. Cornelissen, X. Meng, R. P. Sijbesma, M. Kemerink, Nat. Commun. 2018, 9, 4409
2018
-
[27]
T. D. Cornelissen, I. Urbanaviciute, M. Kemerink, Phys. Rev. B 2020, 101, 214301
2020
-
[28]
Kulkarni, E
C. Kulkarni, E. W. Meijer, A. R. A. Palmans, Acc. Chem. Res. 2017, 50, 1928
2017
-
[29]
M. P. Lightfoot, F. S. Mair, R. G. Pritchard, J. E. Warren, Chemical Communications 1999, 0, 1945. 20
1999
-
[30]
P. J. M. Stals, M. M. J. Smulders, R. Martín-Rapún, A. R. A. Palmans, E. W. Meijer, Chemistry – A European Journal 2009, 15, 2071
2009
-
[31]
P. J. M. Stals, J. C. Everts, R. de Bruijn, I. A. W. Filot, M. M. J. Smulders, R. Martín- Rapún, E. A. Pidko, T. F. A. de Greef, A. R. A. Palmans, E. W. Meijer, Chemistry – A European Journal 2010, 16, 810
2010
-
[32]
K. K. Bejagam, G. Fiorin, M. L. Klein, S. Balasubramanian, J. Phys. Chem. B 2014, 118, 5218
2014
-
[33]
B. A. GIL-VILLEGAS, S. C. McGROTHER and GEORGE JACKSON, Molecular Physics 1997, 92, 723
1997
-
[34]
Vopsaroiu, J
M. Vopsaroiu, J. Blackburn, M. G. Cain, P. M. Weaver, Phys. Rev. B 2010, 82, 024109
2010
-
[35]
Vopsaroiu, P
M. Vopsaroiu, P. M. Weaver, M. G. Cain, M. J. Reece, K. B. Chong, IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control 2011, 58, 1867
2011
-
[36]
J. P. Sethna, K. A. Dahmen, O. Perkovic, 2005, DOI 10.48550/arXiv.cond- mat/0406320
2005 doi
- [37]
-
[38]
T. D. Cornelissen, M. Biler, I. Urbanaviciute, P. Norman, M. Linares, M. Kemerink, Phys. Chem. Chem. Phys. 2019, 21, 1375
2019
-
[39]
C. D. Tan, C. Flannigan, J. Gardner, F. D. Morrison, E. K. H. Salje, J. F. Scott, Phys. Rev. Mater. 2019, 3, 034402
2019
-
[40]
E. K. H. Salje, D. Xue, X. Ding, K. A. Dahmen, J. F. Scott, Phys. Rev. Mater. 2019, 3, 014415
2019
-
[41]
Tückmantel, Scanning Probe Studies of Structural and Functional Properties of Ferroelectric Domains and Domain Walls, Springer Nature, 2021
P. Tückmantel, Scanning Probe Studies of Structural and Functional Properties of Ferroelectric Domains and Domain Walls, Springer Nature, 2021
2021
-
[42]
Casals, G
B. Casals, G. F. Nataf, D. Pesquera, E. K. H. Salje, APL Materials 2020, 8, 011105
2020
-
[43]
Tückmantel, I
P. Tückmantel, I. Gaponenko, N. Caballero, J. C. Agar, L. W. Martin, T. Giamarchi, P. Paruch, Phys. Rev. Lett. 2021, 126, 117601
2021
- [44]
-
[45]
Verbakel, S
F. Verbakel, S. C. J. Meskers, R. A. J. Janssen, H. L. Gomes, M. Cölle, M. Büchel, D. M. de Leeuw, Applied Physics Letters 2007, 91, 192103
2007
-
[46]
Cölle, M
M. Cölle, M. Büchel, D. M. de Leeuw, Organic Electronics 2006, 7, 305
2006
-
[47]
Urbanaviciute, X
I. Urbanaviciute, X. Meng, T. D. Cornelissen, A. V. Gorbunov, S. Bhattacharjee, R. P. Sijbesma, M. Kemerink, Advanced Electronic Materials 2017, 3, 1600530
2017
-
[48]
A. V. Gorbunov, X. Meng, I. Urbanaviciute, T. Putzeys, M. Wübbenhorst, R. P. Sijbesma, M. Kemerink, Physical Chemistry Chemical Physics 2017, 19, 3192
2017
Reviewed August 11, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.