{"id":"cfae149a-5f45-468d-b448-112b89266388","arxiv_id":"2501.09248","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"The energy functions of three repeating FRBs show a universal break around 10^38 erg, which the authors interpret as evidence for starquake triggering in magnetars.","lead":"Three hyperactive repeating FRBs all show a break in their burst energy distribution near 10^38 erg, with different power-law slopes below and above that break. The authors argue this break matches the starquake model, where weak bursts grow on the magnetar surface and strong bursts are limited by crust thickness.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"A single power-law plate-size distribution (Eq. 7) predicts αE2 = 1.5 αE1, but the fitted FRB slopes give αE2/αE1 ≈ 2.5–4.5; the starquake model therefore requires an unstated additional break in N(lp), weakening the claimed dimensionality explanation.","rationale":"The paper's observational core, a common break near 10^38 erg in three large FAST samples with a similar break in SGR 1806-20, is a plausible result, and the comparison with Arecibo data adds credibility. My concern targets the theoretical conversion from the fitted slopes to the starquake prediction. The model as written assumes one plate-size power-law index β but then fits two indices β1 and β2; this is not a harmless notation change. A single β fixes the ratio of the two cumulative slopes to 3/2, and the measured ratio is incompatible with that. The same-β consistency check is exactly what makes the earthquake analogy work, so applying it to the FRB data shows the analogy fails quantitatively unless an additional, unexplained break is inserted into the plate-size distribution. This does not overturn the empirical finding, but it strips the central 'starquakes triggered FRBs' conclusion of one of its main quantitative supports. The reader's correctness-risk label and conditional verdict remain appropriate; I do not see a need to move to reject. My concrete test is a model-comparison check that could settle whether the extra break is required.","tokens_in":18203,"tokens_out":8486,"duration_ms":99888,"concrete_test":"Fit the three FRB energy datasets with an unbinned likelihood using (A) the free broken power law of Eq. (2) and (B) a single-β broken power law that enforces αE2 = 1.5 αE1, which is the prediction of Eq. (10) with one β. Use the same completeness thresholds and data as in the paper, and compare models with ΔAIC or a likelihood-ratio test. If model B is rejected at high significance, the data require a break in the plate-size distribution in addition to the dimensional crossover, and the claimed starquake explanation loses quantitative support.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The theoretical interpretation is internally inconsistent at the level of §4.2, Eqs. (7)–(10). If the fracture-plate sizes follow one power law N(lp) ∝ lp^{-β}, the cumulative energy slopes below and above the break must satisfy αE2 = 1.5 αE1, because αE1 = (β−1)/3 and αE2 = (β−1)/2. For the three FRBs in Table 1 (αE1 = 0.56, 0.35, 0.36), this predicts αE2 = 0.84, 0.53, 0.54, respectively, whereas the fitted values are 1.40, 1.51, 1.61, a discrepancy of many sigma in each case. The paper avoids this by labeling the indices β1 and β2 in Eq. (9) and fitting them separately, which effectively assumes that the plate-size distribution itself changes slope at the same break energy. That extra break in N(lp) is not derived from the starquake model and is not stated as an assumption. The claim that 'difference in dimensionality causes a break' is therefore incomplete: for the observed FRB slopes, the dimensional crossover alone cannot produce the required slope change. The earthquake comparison actually highlights this: earthquakes have αE2/αE1 ≈ 1.47, consistent with a single β, whereas the FRB ratios are ≈2.5–4.5. Thus the universal break may be real, but the starquake interpretation is not uniquely supported and requires an ad hoc second break.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper analyzes the cumulative energy distributions of three hyperactive repeating FRBs (FRB 20121102A, FRB 20201124A, FRB 20220912A) using FAST data, fitting each with a broken power law via MCMC. The authors report a universal break around 10^38 erg in all three sources, with indices αE1 ≈ 0.35–0.56 below the break and αE2 ≈ 1.40–1.61 above it. They compare these with a global earthquake catalog and with SGR 1806-20 X-ray bursts, finding breaks at similar energies. The paper interprets the break as evidence for starquakes in magnetars, arguing that weak bursts grow in three dimensions while strong bursts are confined to two dimensions by the crust thickness, producing a dimensionality-driven change in the energy-function slope.","tokens_in":18544,"tokens_out":6356,"duration_ms":55185,"significance":"The observational identification of a common break at ~10^38 erg in three large, uniformly processed FRB samples is a potentially important contribution to the study of repeating FRB energetics and magnetar physics. The paper makes a commendable effort to homogenize the energy calculations and completeness thresholds across the three sources, and the earthquake and SGR comparisons provide useful physical context. However, the theoretical interpretation is not well supported by the data as presented: the model requires an unstated break in the plate-size distribution, and the break energy is fitted rather than predicted. If the observational result holds, it would motivate further theoretical work, but the current manuscript overstates the support for the starquake trigger.","major_comments":[{"comment":"The model assumes a single power-law plate-size distribution N(lp) ∝ lp^{-β} in Eq. (7), but then fits different exponents β1 and β2 for weak and strong FRBs in Eqs. (9) and (10). If β is the same across the break, the cumulative slopes must satisfy αE2 = 1.5 αE1. For the fitted αE1 values in Table 1 (0.56, 0.35, 0.36), this predicts αE2 = 0.84, 0.53, and 0.54, respectively, far below the fitted values of 1.40, 1.51, and 1.61. The paper does not state or motivate a change in β at the break; without such a change, the dimensional crossover alone cannot produce the observed slope difference. This internal inconsistency makes the proposed explanation ad hoc, and the earthquake comparison actually highlights the problem, because earthquakes have an index ratio of ≈1.47, consistent with a single β, whereas the FRB ratios are ≈2.5–4.5.","section":"§4.2, Eqs. (7)–(10)"},{"comment":"The likelihood function in Eq. (3) treats binned cumulative counts as independent with σcum,i = √Ni. Cumulative bins are strongly correlated, so this underestimates the uncertainties and can bias the parameter estimates; the 1σ errors quoted in Table 1 are likely too small. Moreover, the Kolmogorov-Smirnov test is applied to the same cumulative data used for fitting, so the reported p-values of 0.9999 do not provide independent evidence of goodness of fit. The authors should use a likelihood based on independent energy bins (e.g., Poisson counts in differential bins) or otherwise account for the correlation, and should report a goodness-of-fit test that does not reuse the fitted cumulative curve.","section":"§3.1, Eq. (3)"},{"comment":"The break energy Eb is not predicted by the starquake model; it is fitted from the FRB data and then inserted into Eq. (11) to derive a fracture length L ≈ 2.7×10^4 cm. This is a consistency check rather than a prediction, since the derived length is not independently verified. The abstract's claim that the break 'can be well understood' and that the result 'strongly supports' the starquake trigger therefore overstates the evidential weight. A more convincing case would predict Eb from crustal properties or demonstrate that the derived L and the assumed crust thickness Rc are consistent with independent constraints without tuning η and B.","section":"§4.2, Eq. (11)"},{"comment":"The reduced χ² values for the three FRB fits are 2.06, 2.45, and 1.79, all noticeably above unity, indicating that the broken power law does not fully describe the binned cumulative distributions. The paper does not discuss this discrepancy and instead emphasizes the K-S p-values. The authors should address whether the residuals are systematic (e.g., a smooth curvature) and should report the number of bins and degrees of freedom for each fit, so that the quality of the broken-power-law description can be properly assessed.","section":"Table 1"}],"minor_comments":[{"comment":"The completeness thresholds for FRB 20201124A (2.0×10^36 erg) and FRB 20220912A (1.0×10^36 erg) are stated without derivation; please provide the calculation or cite the original papers for these specific values.","section":"§2"},{"comment":"The exponential cutoff model is written as N(>E) ∝ E^{-αE} e^{-γ}, but γ must multiply E or a cutoff energy for dimensional consistency; as written, the expression is not well defined. Please correct the formula.","section":"§3.2, Eq. (5)"},{"comment":"The notation switches from β in Eq. (7) to β1 and β2 in Eq. (9) without explanation; please clarify that the model allows different power-law exponents for the plate-size distribution below and above the break.","section":"§4.2, Eq. (7)"},{"comment":"There is a typo: 'triggerde' should be 'triggered'.","section":"§6, point (1)"},{"comment":"The break point for FRB 20121102A is quoted as 1.3×10^38 erg in the caption, but Table 1 gives 1.05×10^38 erg; please correct the inconsistency.","section":"Figure 1 caption"},{"comment":"The reduced χ² values should be accompanied by the number of bins and the number of degrees of freedom for each fit.","section":"Table 1"}],"recommendation":"major_revision","confidential_remarks":"The observational dataset is strong and the universal break is an interesting result, but the theoretical interpretation is currently overreaching. The internal inconsistency in the plate-size distribution (β1 ≠ β2 without justification) is a serious flaw that requires a substantial revision of Section 4.2 or a significant softening of the starquake claims. The statistical treatment of the cumulative distribution also needs revision. I recommend major revision; if the authors can either provide a physical mechanism for a break in N(lp) or reframe the paper as an observational finding with a tentative interpretation, it could become publishable."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First: the observational claim is the part to take seriously. A uniform FAST analysis of three bright repeaters shows a break in the cumulative energy function at ~10^38 erg, with indices below and above the break consistent across sources. That is new, and it strengthens the earlier Arecibo hints. The completeness treatment is careful, and the comparison with the SGR 1806-20 X-ray burst break is useful.\n\nThe soft spots are all in the theory section. The stress-test note lands: if the plate-size distribution is a single power law (Eq. 7), the cumulative slopes must satisfy αE2 = 1.5 αE1. The fitted values give ratios 2.5–4.5, so the claimed dimensional crossover alone cannot produce the observed slope change. The paper gets around this by silently switching from β to β1 and β2 in Eq. (9), effectively assuming a second break in the plate-size distribution. That assumption is not stated, not derived, and not tested. The 'explanation' is therefore a consistency check with two new free parameters, not a prediction. The break energy is also fitted, then used to back out a fracture length; nothing predictive is going on there.\n\nA few smaller statistical points: the binned cumulative likelihood in Eq. (3) treats the bins as independent, which is wrong for a cumulative distribution; reduced χ² values of 2.0–2.5 are not great, and the K-S p-values of 0.9999 are suspicious, as if the uncertainties are overestimated. These are fixable in revision, though.\n\nSo: the paper deserves a serious referee. The observational result is likely to stand and will be cited. The theoretical section needs major restructuring before the starquake claim can be taken at face value. I'd send it to review with a clear request to address the β1/β2 issue.","headline":"The three-source universal break near 10^38 erg is a real observational result worth following up; the starquake interpretation has a load-bearing logical gap.","tokens_in":19231,"tokens_out":3139,"would_cite":true,"duration_ms":29992,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Three repeating FRBs share a universal energy break near 10^38 erg, evidence for starquake triggers.","keywords":["fast radio bursts","repeating FRBs","energy function","broken power law","starquake","magnetar crust","Gutenberg-Richter law","FRB 20121102A"],"falsifier":"Measure the energy function of another hyperactive repeater with at least a thousand bursts: if its break is not near $10^{38}$ erg, or if its slopes below and above the break differ from the reported values beyond the quoted uncertainties, the universal-break claim is falsified. A second decisive test would be to check whether the break energy shifts systematically with the assumed magnetic field or crust thickness, as the dimensional argument would predict.","tokens_in":17954,"feed_emoji":"⚡","tokens_out":7830,"duration_ms":63731,"temperature":0.7,"pith_summary":"This paper analyzes the energy distributions of the three most active repeating fast radio bursts — FRB 20121102A, FRB 20201124A, and FRB 20220912A — each with more than a thousand bursts detected by FAST. In all three, the cumulative energy function is not a single power law but a broken power law with a break near $10^{38}$ erg, and the slopes below and above the break are consistent across the three sources. The authors argue that this break mirrors the frequency-magnitude relation of earthquakes, and they show that a starquake model can reproduce it: weak bursts grow in both surface area and crust depth, so their energy scales as size cubed, while strong bursts are confined by crustal thickness and scale as size squared. If correct, this turns the apparent break into evidence that repeating FRBs are triggered by starquakes in magnetar crusts and gives a new way to probe neutron-star crust properties.","feed_headline":"10^38 erg marks a universal break in three repeating FRBs","feed_subtitle":"The break matches starquake predictions and earthquake statistics, a sign that crust cracking triggers these bursts.","key_machinery":"The argument rests on a broken power-law fit to cumulative burst energy, $N(\\ge E) \\propto E^{-\\alpha_{E1}}$ below and $E^{-\\alpha_{E2}}$ above a break energy $E_b$. The physical mechanism is a dimensional-scaling argument adopted from earthquake seismology: fracture plates of size $l_p$ are assumed to follow the power-law distribution $N(l_p) \\propto l_p^{-\\beta}$, and burst energy is taken to scale as $E \\propto l_p^3$ for weak bursts, which can grow in both area and depth, and $E \\propto l_p^2$ for strong bursts, which are confined by the crust thickness $R_c$. This converts the plate-size power law into two different energy-function slopes, $(1-\\beta_1)/3$ and $(1-\\beta_2)/2$, matching the fitted $\\alpha_{E1}$ and $\\alpha_{E2}$. The break occurs where the starquake penetration depth reaches the crust thickness, and the fitted break energy is used to estimate a fracture length, giving a consistency check against expected crust properties.","core_discovery":"The central discovery is that the cumulative energy functions $N(\\ge E)$ of FRB 20121102A, FRB 20201124A, and FRB 20220912A are each well described by a broken power law with break energies of $1.05 \\times 10^{38}$ erg, $1.13 \\times 10^{38}$ erg, and $1.09 \\times 10^{38}$ erg, respectively. The faint-end slopes $\\alpha_{E1}$ take values between $0.35$ and $0.56$, while the bright-end slopes $\\alpha_{E2}$ are close to $1.5$, and the paper finds these agree within uncertainties after systematic effects are considered. The same broken-power-law shape appears in a global earthquake catalog, with a break at magnitude $7.6$, and in the X-ray bursts of SGR 1806-20, with a break near $1.48 \\times 10^{38}$ erg. The paper concludes that the break is intrinsic to the sources, that the dimensional difference between shallow and crust-confined starquakes explains it, and that this supports the starquake trigger hypothesis for FRBs.","pith_inferences":["A testable prediction follows: any newly discovered hyperactive repeater with a comparable burst count should show a break within a factor of a few of $10^{38}$ erg, or the claimed universality would be weakened.","The break energy is fitted rather than predicted, and converting it into a fracture length depends on the assumed beaming factor, conversion efficiency, magnetic field, and breaking strain, so the inferred scale can shift by orders of magnitude under different assumptions.","Because the same dimensional argument predicts breaks in both FRB and magnetar X-ray-burst energy functions, simultaneous radio and X-ray monitoring of a single active repeater could test whether the two breaks trace the same crust physics."],"forward_implications":["Extrapolations of the energy function for bright repeating FRBs must use the steep above-break slope; single-power-law extrapolations would overpredict the rate of the most energetic bursts.","The break near $10^{38}$ erg is robust against telescope sensitivity and completeness effects, since the independent Arecibo sample of FRB 20121102A shows a break at a similar energy.","If repeating FRBs are starquake-triggered, their energy functions should resemble those of magnetar X-ray bursts, and the observed SGR 1806-20 break at $1.48\\times10^{38}$ erg is consistent with that expectation.","The starquake interpretation connects FRB statistics to neutron-star crust physics, including the maximum breaking strain and the crustal thickness."],"supporting_citations":[{"why":"FAST sample of 1652 bursts from FRB 20121102A, the primary dataset for its energy function.","marker":"Li et al. 2021"},{"why":"FAST sample of 1863 bursts from FRB 20201124A, with completeness estimates and energy calibration.","marker":"Xu et al. 2022"},{"why":"FAST sample of 1076 bursts from FRB 20220912A, used for its energy function.","marker":"Zhang et al. 2023"},{"why":"Independent Arecibo analysis showing a similar break around $10^{38}$ erg for FRB 20121102A.","marker":"Jahns et al. 2023"},{"why":"Worldwide earthquake catalog complete for magnitude $\\ge 7.0$, the comparison dataset for the Gutenberg-Richter break.","marker":"Pacheco et al. 1992"},{"why":"Predicted the earthquake energy-function break from a change in rupture dimensionality.","marker":"Rundle 1989"},{"why":"Provides the starquake energy scaling used to relate fracture size to released energy.","marker":"Wang et al. 2018"},{"why":"Provides the FRB-to-X-ray energy efficiency $\\eta \\sim 10^{-4}$ used to infer fracture length.","marker":"Mereghetti et al. 2020"}],"fun_headline_variants":["Universal energy break links FRBs to starquakes","Three repeating FRBs share starquake-like energy cutoff","FRB energy breaks at 10^38 erg, mimicking earthquakes","Starquakes implicate universal FRB energy break","Repeating FRBs show earthquake-like energy limit"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The explanation assumes that fracture-plate sizes follow a pure power-law distribution and that the weak-to-strong transition is exactly the change from $l^3$ to $l^2$ energy scaling at the fitted break energy, while the break value itself is not derived from first principles.","fun_headline_variants_meta":{"raw":{"variants":["Universal energy break links FRBs to starquakes","Three repeating FRBs share starquake-like energy cutoff","FRB energy breaks at 10^38 erg, mimicking earthquakes","Starquakes implicate universal FRB energy break","Repeating FRBs show earthquake-like energy limit"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000772,"raw_usage":{"total_tokens":3468,"prompt_tokens":1047,"completion_tokens":2421,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":663,"completion_tokens_details":{"reasoning_tokens":2357}},"tokens_in":663,"tokens_out":2421,"duration_ms":16648,"temperature":1.0,"reasoning_tokens":2357,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T20:06:45.578233+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the energy function of another hyperactive repeater with at least a thousand bursts: if its break is not near $10^{38}$ erg, or if its slopes below and above the break differ from the reported values beyond the quoted uncertainties, the universal-break claim is falsified. A second decisive test would be to check whether the break energy shifts systematically with the assumed magnetic field or crust thickness, as the dimensional argument would predict.","supporting_citations":[{"cited_title":"F., Scholz , C","cited_arxiv_id":null,"evidence_quote":"Worldwide earthquake catalog complete for magnitude $\\ge 7.0$, the comparison dataset for the Gutenberg-Richter break."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Predicted the earthquake energy-function break from a change in rupture dimensionality."}],"review_version":1}