{"id":"a1c3079c-5765-40b2-959e-3b5b17f7d015","arxiv_id":"2507.05112","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A variational first-principles framework reproduces known small polarons and predicts weakly bound large polarons with symmetry-broken degeneracies in LiF, MgO, Li2O2, and GaN.","lead":"Researchers combined a variational polaron equation with density functional theory inputs to compute self-trapped electron and hole polarons in five polar materials. The method reveals multiple symmetry-broken polaron states and handles very large polarons, matching earlier calculations where comparisons exist.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The GaN hole-polaron result is the load-bearing weakness: Eq. (47) admits a missing LA-phonon correction, and the -7 meV extrapolated binding could change sign if that term is included.","rationale":"The reader's weakest assumption identifies the same soft spot: the Eq. (37) extrapolation for large polarons in MgO/GaN and the uncorrected LA divergence in Eq. (47). I agree with that choice. My stress-test sharpens it: the existence claim for GaN is the most fragile part of the central 'all five materials' statement. The -7 meV value is smaller than the likely uncertainty introduced by the omitted piezoelectric LA correction, and the extrapolation window is narrow. Implementing the LA geff correction, which the authors themselves suggest, would settle the matter. This concern does not invalidate the paper's methodological contribution: the small-polaron results are well validated against independent calculations, and the long-range LO correction is a separately tested ingredient. Therefore the reader's CONDITIONAL verdict remains appropriate: the GaN large-polaron numbers should be treated as provisional until the LA correction or an explicit error bar is provided.","tokens_in":30266,"tokens_out":19654,"duration_ms":217390,"concrete_test":"Recompute the zb-GaN hole polaron including a long-range correction for the piezoelectric LA branch, derived on the same footing as Eq. (25) but using the |q|^{-1/2} acoustic vertex, and redo the Np^{-1/3} extrapolation of Epol. If the corrected extrapolated value becomes positive or shifts by more than about 5 meV, the reported GaN self-trapped hole polaron is not robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The weakest load-bearing point is the claimed self-trapped hole polaron in GaN, whose binding energy is an order of magnitude smaller than the LO phonon scale yet rests on an extrapolation that omits an acknowledged finite-size term. For zb-GaN, Epol = -7 meV is obtained by fitting Eq. (37) to meshes 75^3-85^3, i.e., Np^{-1/3} from 0.0133 to 0.0118, and extrapolating to 0 over roughly eight times the sampled interval. In the same section, Eq. (47) shows that for the piezoelectric LA branch the deformation field B_qnu scales as |q|^{-1/2}, and the text states that a complete treatment would require an extra geff correction for LA phonons. Because that missing LA term has the same Np^{-1/3} scaling as the fitted correction in Eq. (37), it is absorbed into the slope and biases the intercept; with a value as small as -7 meV the intercept could change sign. Thus the central claim that VarPEq yields self-trapped polaron solutions in all five materials is not yet supported for GaN.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript applies the variational polaron equations (VarPEq) framework, implemented in ABINIT, to compute self-trapped electron and hole polarons in LiF, MgO, Li2O2, wurtzite GaN, and zincblende GaN. It introduces an effective long-range electron-phonon correction for finite q-meshes, an energy-filtering procedure for large polarons, and an orthogonalization strategy for enumerating symmetry-broken degenerate solutions. The reported small-polaron binding energies agree closely with earlier calculations for LiF and Li2O2, and more approximately for MgO, while the large MgO electron polaron matches the strong-coupling Fröhlich estimate. The central claim is that the method yields self-trapped polaronic solutions in all five materials, including weakly bound large hole polarons in GaN.","tokens_in":30539,"tokens_out":8378,"duration_ms":94017,"significance":"The paper is potentially significant as a practical variational alternative to supercell-based polaron calculations: it avoids charged supercells, handles degenerate band edges, and reaches large polaron sizes via energy filtering. Its strengths include external benchmarks against Sio et al., Falletta and Pasquarello, Radin and Siegel, and the Fröhlich strong-coupling formula, as well as an open implementation in ABINIT and data deposited on Materials Cloud. The authors also candidly state one central limitation: for piezoelectric GaN, an LA-phonon correction needed for a complete finite-size treatment is missing. Because the GaN binding energies are of order a few meV and rest on a long extrapolation, the claim that self-trapped hole polarons are obtained in both GaN polymorphs is not yet fully supported. With the GaN point addressed, the paper would make a solid contribution.","major_comments":[{"comment":"The claimed large hole polarons in GaN are not yet supported by the evidence presented. The zb-GaN value Epol = -7 meV is obtained by fitting Eq. (37) to data on meshes 75^3-85^3, i.e., Np^{-1/3} from roughly 0.0133 to 0.0118, and extrapolating to Np^{-1/3} = 0 over about eight times the sampled interval. In the same section, Eq. (47) and the surrounding text state that for the piezoelectric LA branch in GaN the deformation field B_{qν} diverges as |q|^{-1/2} and that a complete finite-size treatment requires an additional geff correction for LA phonons, which is not implemented. Because a missing LA-phonon term enters with the same Np^{-1/3} scaling as the fitted slope in Eq. (37), it can bias the intercept; with a reported binding energy of only -7 meV, the sign of Epol is not robust. The same concern applies to the wz-GaN value of -9 meV. The abstract and conclusion therefore overstate the result that self-trapped hole polarons are found in both GaN polymorphs. I ask the authors to either implement the LA correction or present the GaN numbers as tentative, and to provide a quantitative estimate of the extrapolation uncertainty.","section":"Section IV C, Eq. (37), Fig. 9(d), Eq. (47)"},{"comment":"No error bars, residuals, or confidence intervals are reported for any of the extrapolated Epol and εloc values, despite the fact that the large-polaron results in MgO and GaN are tens of meV or less and the extrapolation is made over a long lever arm in Np^{-1/3}. For MgO the agreement with the Fröhlich strong-coupling estimate (-23/-24 meV) provides an external check, but for GaN no such check exists. Please report fit residuals and parameter uncertainties and, for at least one system, test the sensitivity of the intercept to including an O(Np^{-1}) term or to dropping the smallest mesh.","section":"Section IV C, Figs. 8 and 9"}],"minor_comments":[{"comment":"The caption lists panels (a), (b), (d), and (f), but the figure contains four panels and the text refers to Fig. 5(c) and Fig. 5(d); the panel labels in the caption should be corrected to (a)-(d).","section":"Fig. 5 caption"},{"comment":"The exponent in the asymptotic formula is written with the symbol γ, while the subsequent text uses μ for the same exponent (μ = 1/2 for piezoelectric materials); please make the notation consistent.","section":"Eq. (47)"},{"comment":"There are typographical errors: the Section II heading contains 'V ariational' and the text near Eq. (28) contains 'precondioner'; both should be corrected.","section":"Section II, heading and Section II E"},{"comment":"The 30% difference for the MgO hole polaron relative to Ref. 17 is acknowledged and plausibly explained, but I suggest stating more explicitly in the text that this is a known quantitative limitation rather than part of the claimed agreement with previous studies.","section":"Section IV A, Table II"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the scope of the journal, and the open-source implementation plus data availability are clear assets. The GaN issue is significant but appears fixable through additional analysis or carefully qualified wording; I would not reject the paper on this basis. My recommendation of major revision is driven by the load-bearing nature of the GaN large-polaron claim, which rests on an acknowledged incomplete finite-size treatment."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a solid methods paper from the ABINIT team, extending their VarPEq framework with three genuinely useful pieces: an analytic long-range correction geff, a procedure for finding degenerate symmetry-broken polaronic states via orthogonalization against previous solutions and their translated images, and energy filtering to handle large polarons. The small-polaron benchmarks are the strongest part: LiF hole and Li2O2 electron match Sio et al. to about 2%, Li2O2 hole matches Radin and Siegel, and the MgO large electron polaron (-23 meV) matches the Fröhlich strong-coupling estimate. That is real externally benchmarked work, not circular.\n\nThe weak point is the GaN large hole polarons. The -7 meV (zb) and -9 meV (wz) values come from a linear fit in Np^{-1/3} over a very narrow window, with no error bars. More concerning, the authors acknowledge in Eq. (47) that the piezoelectric LA branch makes B diverge as |q|^{-1/2} and that a complete treatment would require an extra geff-type correction for LA phonons. That missing term has the same Np^{-1/3} scaling as the fitted correction, so it gets absorbed into the slope and can bias the intercept. At -7 meV, the intercept could plausibly change sign. The MgO electron polaron is safer because it agrees with the Fröhlich model, but it also lacks error bars. The energy-filtering validation is done only for LiF, though there the extrapolated values agree with the un-filtered calculation. The discussion of the MgO hole discrepancy with Falletta and Pasquarello is honest; the 30% difference is plausibly explained by linear coupling and neglect of Debye-Waller self-energy.\n\nOverall, the method is clearly explained, the implementation is in a public release (ABINIT v10.4), and the data is on Materials Cloud. The core formalism comes from their PRB 2022 paper, but the new pieces are non-trivial and the benchmarks are external. I would send this to peer review. The referee should insist on error bars for the large-polaron energies and either a treatment or an explicit bound on the LA-phonon correction in GaN. If the GaN numbers are presented as qualitative, fine; as quantitative values, they are not yet supported. This paper is for people working on first-principles polaron methods and polaron transport, and it deserves a serious referee.","headline":"Solid extension of VarPEq with strong small-polaron benchmarks; the GaN large-polaron binding energies need error bars and a proper LA-phonon correction before they can be quoted.","tokens_in":686,"tokens_out":843,"would_cite":true,"duration_ms":27408,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["71.38.-k","78.20.Bh"],"model":"deepseek-v4-flash","headline":"Electrons and holes self-trap in all five materials studied, with binding energies between -4.89 eV and -7 meV, and degenerate band edges break symmetry into enumerable partner states such as D4h in LiF and D3d in GaN.","keywords":["polaron","self-trapped polaron","variational polaron equations","electron-phonon coupling","symmetry breaking","density functional theory","long-range correction","energy filtering"],"falsifier":"Recompute the GaN hole polarons with the long-range piezoelectric (LA-mode) correction included on the same footing as the Fröhlich correction, or with larger supercells and a different extrapolation form, and check whether the extrapolated -7 and -9 meV binding energies survive; a shift comparable to the binding energy itself would mean GaN self-trapping is not yet established.","tokens_in":30050,"feed_emoji":"⚛️","tokens_out":13606,"duration_ms":126645,"temperature":0.7,"pith_summary":"A self-trapped polaron is a charge carrier that digs its own potential well by distorting the surrounding lattice, and this paper claims that one variational scheme — the variational polaron equations fed with density-functional band structures, phonons, and electron-phonon couplings — finds such states in every material it tries: hole polarons in LiF and MgO, electron and hole polarons in Li2O2, and hole polarons in both GaN polymorphs. The binding energies span from -4.89 eV for the Li2O2 electron polaron to -7 meV for the zincblende GaN hole polaron, and the method finds that degenerate band edges spontaneously break the crystal symmetry into distinct partner states. A sympathetic reader would care because these binding energies and spatial symmetries decide whether charge moves by band transport or by hopping, and the paper supplies ab initio numbers, a systematic way to enumerate degenerate polarons, and a filtering trick that makes even 100-Bohr-radius polarons affordable. The results match earlier ab initio studies where those exist.","feed_headline":"From -4.89 eV to -7 meV: polarons self-trap in five solids","feed_subtitle":"A variational first-principles method reproduces known polaron states and finds symmetry-broken ones in all five test crystals.","key_machinery":"The load-bearing object is the polaron binding-energy functional, whose coefficients $A_{nk}$ (charge) and $B_{q\\nu}$ (lattice distortion) are optimized by preconditioned conjugate gradients with the single-charge normalization enforced through the Lagrange multiplier $\\varepsilon_{\\mathrm{loc}}$. Three devices carry the argument: the long-range correction $g^{\\mathrm{eff}}_\\nu$, which replaces the divergent $q = \\Gamma$ electron-phonon element by the Fröhlich-type coupling averaged over a sphere with the microzone volume, removing most of the finite-size error; the orthogonalization step, applied to each new trial state against prior solutions and their translated images, which systematically uncovers degenerate symmetry-broken partners; and energy filtering, which keeps only bands within a window $\\Delta\\varepsilon$ of the band edge and makes meshes up to 95×95×95 affordable. The generalized Fröhlich model serves as the benchmark for the large polarons, and the strong-coupling Fröhlich formula reproduces the MgO electron polaron binding energy.","core_discovery":"Minimizing the binding energy functional $E_{\\mathrm{pol}}(A,B) = E_{\\mathrm{el}}(A) + E_{\\mathrm{ph}}(B) + E_{\\mathrm{e-ph}}(A,B)$ by self-consistent gradient descent, with electronic coefficients $A_{nk}$ and lattice-distortion coefficients $B_{q\\nu}$ as independent variables and DFT-supplied bands, phonons, and electron-phonon matrix elements as input, the authors report self-trapped solutions in all five materials, with extrapolated binding energies from -4.89 eV to -7 meV. The central result beyond the numbers is structural: when the relevant band edge is degenerate, the polaron breaks the crystal point group, and the degenerate partner states — D4h in LiF, D3d and D2h in MgO, C2v and D3h in Li2O2, D2h in wurtzite GaN, D3d in zincblende GaN — can be enumerated by orthogonalizing each new solution against all previously found ones and their translations. Two computational additions make this practical: an analytic long-range correction that accounts for the Fröhlich divergence at $q = \\Gamma$, speeding finite-size convergence without changing the extrapolated energies, and an energy-filtering window that restricts the variational space to states near the band edge, reducing cost by up to roughly three orders of magnitude for MgO. The method reproduces earlier theoretical results for LiF, Li2O2, and MgO, including the strong-coupling Fröhlich value for the MgO electron polaron.","pith_inferences":["The reported -7 and -9 meV GaN binding energies rest on a linear $N_p^{-1/3}$ extrapolation over a narrow supercell window while an uncorrected piezoelectric LA-phonon divergence remains in the deformation field; a long-range correction for those LA modes, which the paper sketches but does not implement, could shift the values by an amount comparable to their magnitude.","The Li2O2 electron polaron's 11 eV localization energy exceeds the material's band gap, an artifact of restricting the response to conduction states; extending the formalism to the full conduction-valence manifold might revise the binding energy, not just the localization energy.","Energy-filtering convergence is cross-checked only for LiF; for MgO and GaN the chosen windows are justified by convergence of $\\varepsilon_{\\mathrm{loc}}$ and by the variational guarantee, leaving the possibility of a small systematic bias in the large-polaron numbers.","Because the electron-phonon matrix elements are gauge-dependent and the calculation must disable Brillouin-zone symmetry, extending this approach to high-throughput material scans will likely require gauge-recovery techniques to keep the cost manageable."],"forward_implications":["Self-trapped polarons form in every material tested, from strongly bound small polarons (-4.89 eV in Li2O2) to weakly bound large polarons (-7 to -9 meV in GaN), so the variational framework spans the full small-to-large polaron range.","Degenerate band edges spontaneously break the polaron's symmetry and generate several nearly degenerate states; the orthogonalization procedure can enumerate them, removing a blind spot of self-consistent-field methods.","The long-range correction accelerates finite-size convergence while leaving the extrapolated infinite-size binding energies essentially unchanged.","Energy filtering reproduces the un-filtered LiF electron polaron to within about 1 meV while cutting computational cost by factors of roughly 3 to 900 across the materials, opening very large polarons to ab initio treatment.","Because the weakly bound GaN polarons have binding energies below the most coupled phonon frequencies, their full description will require going beyond self-trapping to dynamical effects."],"supporting_citations":[{"why":"Origin of the reciprocal-space polaron formalism and effective polaron Hamiltonian that the variational method builds on, and source of the LiF hole and Li2O2 electron reference values.","marker":"[11,12]"},{"why":"Defines the variational polaron equations and the gradient-based self-consistent optimization that this paper extends to real materials.","marker":"[14]"},{"why":"Real-space DFT studies of small polarons that supply the MgO hole polaron comparison.","marker":"[17,18]"},{"why":"Generalized Fröhlich model with degenerate bands; predicts the D3d MgO polaron and the trigonal-antiprism zincblende-GaN symmetry breaking.","marker":"[29]"},{"why":"Earlier Li2O2 hole polaron calculation that the PBEsol results closely match.","marker":"[53]"},{"why":"Fourier interpolation of electron-phonon scattering potentials with analytic dipole and quadrupole long-range treatment, used to parametrize the method.","marker":"[49]"},{"why":"Previous report of the D2h MgO polaron solution that this work also finds.","marker":"[59]"},{"why":"Strong-coupling Fröhlich formula that reproduces the MgO large electron polaron binding energy.","marker":"[37,38]"},{"why":"Hosts the open-source implementation of the variational method used for the calculations.","marker":"[35,36]"}],"fun_headline_variants":["Variational method reveals symmetry-broken polarons in five crystals","First-principles variational polarons: five materials, new symmetry states","Polarons self-trap across -4.89 eV to -7 meV in five solids","Variational method predicts self-trapped polarons in five solids","Self-trapped polarons in five solids: variational method works"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the linear finite-size extrapolation $E_{\\mathrm{pol}}(N_p) = E^\\infty_{\\mathrm{pol}} + \\alpha N_p^{-1/3}$ gives accurate infinite-size binding energies for the MgO and GaN large polarons, even though the GaN deformation field still contains an uncorrected long-wavelength piezoelectric divergence and the sampled supercell-size window is narrow.","fun_headline_variants_meta":{"raw":{"variants":["Variational method reveals symmetry-broken polarons in five crystals","First-principles variational polarons: five materials, new symmetry states","Polarons self-trap across -4.89 eV to -7 meV in five solids","Variational method predicts self-trapped polarons in five solids","Self-trapped polarons in five solids: variational method works"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001426,"raw_usage":{"total_tokens":5803,"prompt_tokens":1041,"completion_tokens":4762,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":657,"completion_tokens_details":{"reasoning_tokens":4666}},"tokens_in":657,"tokens_out":4762,"duration_ms":34235,"temperature":1.0,"reasoning_tokens":4666,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T19:32:15.194693+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the GaN hole polarons with the long-range piezoelectric (LA-mode) correction included on the same footing as the Fröhlich correction, or with larger supercells and a different extrapolation form, and check whether the extrapolated -7 and -9 meV binding energies survive; a shift comparable to the binding energy itself would mean GaN self-trapping is not yet established.","supporting_citations":[{"cited_title":"Vasilchenko , author A","cited_arxiv_id":null,"evidence_quote":"Defines the variational polaron equations and the gradient-based self-consistent optimization that this paper extends to real materials."},{"cited_title":"Vasilchenko \\ and\\ author X","cited_arxiv_id":null,"evidence_quote":"Generalized Fröhlich model with degenerate bands; predicts the D3d MgO polaron and the trigonal-antiprism zincblende-GaN symmetry breaking."},{"cited_title":"Kokott , author S","cited_arxiv_id":null,"evidence_quote":"Previous report of the D2h MgO polaron solution that this work also finds."}],"review_version":1}