{"id":"b65d8116-0e5c-4d8a-af04-09b68d0c15f4","arxiv_id":"2412.12666","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Kinetic Monte Carlo simulations of the organic ferroelectric BTA predict power-law distributed dipole avalanches and self-organized criticality below about 175 K, while measured current noise in BTA films stays below the experimental detection threshold.","lead":"This paper uses computer simulations to study how the organic ferroelectric BTA reverses its polarization, finding jerky avalanche-like switching events whose sizes follow a power-law distribution at low temperatures. The authors also tried to measure these 'Barkhausen' clicks in real BTA samples, found none, and argue this is consistent with their simulated events being too small to detect.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The paper compares full-loop avalanche exponents to constant-stress mean-field values (tau=1.5, epsilon=4/3) instead of the loop-integrated values (tau=2, epsilon=5/3) it cites, so the reported agreement may not support the SOC claim.","rationale":"The reader's conditional verdict is appropriate, but the single most load-bearing concern is not primarily the analysis robustness identified in the reader's weakest_assumption; it is the benchmark mismatch. The reader's rationale does mention the benchmark issue, so there is partial agreement. The paper's introduction correctly distinguishes constant-stress (tau=1.5, epsilon=4/3) from stress-integrated (tau=2, epsilon=5/3) mean-field exponents, yet the results are compared to the constant-stress values despite collecting events over full hysteresis loops. This is an internal inconsistency: if the fitted exponents are 1.5 and 1.33, they contradict the paper's own loop-integrated theory, and the SOC claim as stated is unsupported. The threshold and fit-range choices are a secondary concern because they could bias the exponents, but the benchmark mismatch is sufficient by itself to invalidate the central claim unless the authors clarify that events were actually taken from a narrow coercive-field window or re-benchmark to the loop-integrated values. A concrete reanalysis with maximum-likelihood fitting would settle whether the discrepancy is real or an artifact of the comparison. The experimental null result and the simulation code are not at issue; the concern is purely about the interpretation of the extracted exponents. Therefore the reader's CONDITIONAL verdict remains the right call: the paper should not be accepted as is, but a clarification or reanalysis could resolve the concern.","tokens_in":18237,"tokens_out":7189,"duration_ms":62495,"concrete_test":"Re-analyze the raw event data (or rerun the kMC simulations) to compute the avalanche size and energy distributions from (a) the complete hysteresis loop and (b) a narrow field window around the coercive field, using maximum-likelihood power-law fitting with a fixed lower cutoff. Compare the full-loop exponents to tau=2 and epsilon=5/3, and the coercive-window exponents to tau=1.5 and epsilon=4/3. If the full-loop fit yields tau~2 and epsilon~5/3, the SOC claim is supported; if it yields tau~1.5 and epsilon~4/3 while events are taken from the full loop, the central claim is not supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that BTA exhibits self-organized criticality below ~175 K rests on the extracted power-law exponents matching mean-field predictions. However, the paper's own introduction states that for events measured over a whole hysteresis loop, the stress-integrated mean-field exponents are tau=2 and epsilon=5/3, whereas the constant-stress values are tau=1.5 and epsilon=4/3. The simulations collect events from full field ramps (Figure 2, SI S2), yet the horizontal lines in Figure 6 and the explicit values in Figure S5.3 correspond to the constant-stress values (1.5 and 1.33). If the reported exponents are indeed near 1.5 and 1.33, they do not match the loop-integrated benchmark; they match only the constant-stress theory, which is not the experimental condition of a full hysteresis loop. This internal inconsistency directly undermines the abstract's claim. The additional analysis choices (manual threshold, visually chosen fit range, weight-event merging) compound the issue, but the benchmark mismatch is the primary load-bearing concern because even a perfectly robust fit would be compared against the wrong theory.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":18410,"tokens_out":4656,"duration_ms":43849,"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":[{"comment":"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.","section":"Introduction and Figure 6 / SI S5.3"},{"comment":"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.","section":"SI S2 and Section 3 (threshold and fitting procedure)"},{"comment":"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.","section":"SI S5 (disorder dependence)"}],"minor_comments":[{"comment":"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.","section":"Section 4 and Figure 6b"},{"comment":"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":"SI S2"},{"comment":"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.","section":"Section 4"},{"comment":"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":"SI S5.3 caption"},{"comment":"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.","section":"Section 5"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is well within the scope of the journal and the simulation methodology is described in detail. The main issue is the mismatch between the stress-integrated nature of the simulated events and the constant-stress mean-field exponents used for comparison; this is a fixable but load-bearing error. In addition, the analysis lacks uncertainty quantification, which is essential when the central claim rests on fitted exponents. I would encourage the authors to resubmit after correcting the benchmark and adding robustness checks."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing: this is a genuine first — a kMC study of Barkhausen noise in BTA, with a new experimental null result reported honestly. The SOC claim in the abstract should not be taken at face value. The reader's stress-test note is right about the benchmark: the paper's own introduction gives tau=2 and epsilon=5/3 for whole-loop integrated mean-field events, but Figure 6 and S5.3 compare to 1.5 and 1.33. If the data are full-loop, that's the wrong reference. The only escape would be if all events used in the fits are effectively confined to a narrow field window around the coercive field, so the constant-stress limit applies. The authors mention such a filtering option in SI S2 but don't clearly say that's what the temperature sweep uses. They need to fix that.\n\nWhat's good: the model is described in enough detail to reproduce; there are finite-size checks showing avalanche size saturates with column length and depends weakly on lateral width, which supports quasi-1D propagation; the disorder sweep maps out a regime where exponents are stable, with behaviour reminiscent of RFIM. The experimental null result is handled well: noise floor characterized, filament artifacts identified, and an order-of-magnitude estimate shows the expected currents sit orders below sensitivity. That part is solid and citable on its own. The circularity concern in the reader notes doesn't land: the model parameters are from earlier work, and the exponents are fitted, not imposed.\n\nSoft spots beyond the benchmark: event threshold is set by hand just above thermal noise, fit range by eye, and the merging of runs uses a 'weight-event' correction that duplicates low-end events. There are no error bars or alternative fit tests, and no code or raw event data released. For a power-law claim, those choices are material. They might move the fitted tau by a few tenths, which is exactly the difference between the quoted values and either benchmark. I don't think it's fatal — the distributions do look power-law-like over a decade — but the central claim is more fragile than the writing suggests.\n\nWho this is for: people working on organic ferroelectrics or crackling noise will want to read it, especially for the null result. A competent referee can fix the benchmark and ask for robustness checks; the paper deserves that step. I'd send it to review, with the request that the authors clarify the benchmark, add uncertainty estimates and alternative fits, and publish event-level data.","headline":"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.","tokens_in":19004,"tokens_out":7863,"would_cite":true,"duration_ms":69992,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["77.80.-e","77.80.Fm","05.65.+b"],"model":"deepseek-v4-flash","headline":"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.","keywords":["Barkhausen noise","crackling noise","self-organized criticality","organic ferroelectric","BTA","kinetic Monte Carlo","power-law exponents","ferroelectric switching"],"falsifier":"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.","tokens_in":17995,"feed_emoji":"⚡","tokens_out":6198,"duration_ms":52572,"temperature":0.7,"pith_summary":"The paper claims that polarization reversal in the columnar organic ferroelectric BTA proceeds through dipole avalanches whose sizes and energies are power-law distributed, matching mean-field predictions of self-organized criticality below about 175 K. Above that temperature the extracted exponents grow, which the authors interpret as a thermal creep regime in which many small thermally triggered events replace field-driven cascades. The same kinetic Monte Carlo simulations produce only tiny events, several orders of magnitude below the current noise floor of a sensitive experimental setup, which the authors offer as the reason no Barkhausen noise is seen in real BTA samples. If correct, this identifies the microscopic switching unit of BTA as a subcolumn within a single supramolecular column and explains the silence of the material in Barkhausen-noise measurements.","feed_headline":"BTA ferroelectric switching is a cascade of tiny avalanche events","feed_subtitle":"Kinetic Monte Carlo matches mean-field exponents and explains why no Barkhausen noise is seen.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Supplies the mean-field plasticity model and the predicted exponents 1.5 and 4/3 that the simulated BTA exponents are compared against.","marker":"[12]"},{"why":"Extends the mean-field plasticity picture to switching in soft ferroic materials, giving the framework for interpreting BTA avalanches as crackling.","marker":"[13]"},{"why":"Establishes the subcolumn/hysteron picture of BTA disorder and provides the basis of the kinetic Monte Carlo model the simulations build on.","marker":"[25]"},{"why":"Reaction-field method that approximates long-range electrostatics beyond the interaction cutoff, needed to compute flipping energies.","marker":"[33]"},{"why":"Thermally activated nucleation-limited switching model linking coercive field to temperature and sweep frequency, used to explain the exponent trends.","marker":"[34]"},{"why":"Random-field Ising model classification of disorder regimes, used to interpret why exponents are stable only for intermediate disorder.","marker":"[36]"},{"why":"Successful Barkhausen noise measurement in P(VDF-TrFE) with the same setup, serving as the positive control for the null result in BTA.","marker":"[37]"}],"fun_headline_variants":["Tiny avalanche cascades explain quiet BTA ferroelectric switching","BTA ferroelectric: critical avalanches too small to detect","No Barkhausen crackle in BTA because flips are minuscule","BTA ferroelectric switches by 1D avalanches that stay silent","Critical avalanches in BTA ferroelectric: too quiet for Barkhausen"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Tiny avalanche cascades explain quiet BTA ferroelectric switching","BTA ferroelectric: critical avalanches too small to detect","No Barkhausen crackle in BTA because flips are minuscule","BTA ferroelectric switches by 1D avalanches that stay silent","Critical avalanches in BTA ferroelectric: too quiet for Barkhausen"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000649,"raw_usage":{"total_tokens":2978,"prompt_tokens":947,"completion_tokens":2031,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":563,"completion_tokens_details":{"reasoning_tokens":1935}},"tokens_in":563,"tokens_out":2031,"duration_ms":13331,"temperature":1.0,"reasoning_tokens":1935,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T13:51:34.413590+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reaction-field method that approximates long-range electrostatics beyond the interaction cutoff, needed to compute flipping energies."},{"cited_title":"Barkhausen noise in the organic ferroelectric copolymer P(VDF:TrFE)","cited_arxiv_id":"2412.12671","evidence_quote":"Successful Barkhausen noise measurement in P(VDF-TrFE) with the same setup, serving as the positive control for the null result in BTA."}],"review_version":1}