{"id":"cad89ed2-2cfc-4af7-8b8e-50748016295e","arxiv_id":"2607.05890","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.5,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Microscopic Usadel calculation of Bc1 for dirty superconducting thin-film strips yields a larger vortex self-energy than Pearl–London and no single width-independent cutoff.","lead":"This paper computes the lower critical field of dirty superconducting thin-film strips from the Usadel equations without a phenomenological vortex-core cutoff. It gives a microscopic Gibbs-energy route to Bc1 for nanostrips used in quantum and sensing devices.","discovery_kind":"new_method","skeptic_critique":{"model":"grok-4.5","headline":"The claimed excess vortex self-energy over Pearl–London hinges on the microscopic Gibbs functional reducing correctly to the Pearl electromagnetic energy; any mismatch in the magnetic or kinetic terms would make the excess (and the single-cutoff claim) an artifact.","rationale":"The Reader correctly flags the dirty-limit 2-D Usadel/thin-film reduction as a modeling assumption that must hold for the core and edge deformation to be trustworthy; that is a real scope limitation. The more load-bearing internal concern, however, is whether the Gibbs functional used to extract the self-energy is thermodynamically consistent with the Pearl–London energy that serves as the benchmark. Without that consistency the numerical finding “larger than Pearl–London and not a single cutoff” cannot be trusted even inside the dirty 2-D regime. Because the abstract already frames the result as a cutoff-free microscopic calculation whose main surprise is the excess self-energy, the free-energy reduction is the single point that most directly supports or undermines the strongest claim. The Reader’s CONDITIONAL verdict already awaits full methods and quantitative benchmarks; the free-energy check is precisely one of those benchmarks, so no change of verdict category is required. If the London-limit test passes, the paper’s claim stands on firmer ground; if it fails, the excess-self-energy statement must be withdrawn or re-derived.","tokens_in":2164,"tokens_out":661,"duration_ms":45230,"concrete_test":"Extract or re-derive the London (rigid-|Δ|, large-κ) limit of the paper’s free-energy functional for both the Meissner state and a Pearl vortex in a wide strip; evaluate the energy difference with the same current density that enters the Pearl–London formula. If the two expressions differ by more than a few percent of the expected core-energy scale, the claimed excess self-energy and the single-cutoff conclusion are unreliable.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper’s central quantitative claim—that the vortex self-energy obtained from the Gibbs difference of self-consistent Usadel solutions is larger than the naive Pearl–London estimate and cannot be absorbed into a London logarithm with one width-independent cutoff—rests on a direct numerical comparison of free energies. The Usadel free-energy functional (Matsubara sum over the anomalous Green’s function plus the electromagnetic contribution) must recover the standard Pearl–London energy when the order parameter is taken rigid outside a core and the thin-film nonlocal kernel is used. If the implementation of the magnetic term at fixed applied field (or the subtraction that isolates the vortex self-energy) omits the applied-field interaction, double-counts the supercurrent kinetic energy already present in the Usadel current, or mishandles the Pearl kernel under the 2-D reduction, the reported excess is not physical. The abstract states that Bc1 follows from the Gibbs-energy difference and that the self-energy exceeds Pearl–London, but does not exhibit the explicit free-energy expression or its London-limit check; that reduction is therefore the least secure link in the argument.","agreement_with_reader":"partial"},"referee_report":{"model":"grok-4.5","summary":"The manuscript formulates a microscopic theory of the lower critical field Bc1 for a dirty superconducting thin-film strip by solving the two-dimensional Usadel equations in the film plane, with the applied field entering the gauge-invariant momentum. Self-consistent Meissner and single-vortex solutions are obtained at fixed applied field, and Bc1 is identified with the field at which their Gibbs free energies cross. The calculation resolves the vortex core and its finite-width deformation without a phenomenological cutoff. The central quantitative claims are that the resulting vortex self-energy exceeds the naive Pearl–London estimate and that this excess cannot, in general, be absorbed into a London logarithm with a single width-independent cutoff. The formulation is stated to apply at any T < Tc and to provide a microscopic basis for Bc1 in superconducting nanostrips.","tokens_in":2375,"tokens_out":1230,"duration_ms":35548,"significance":"If the free-energy comparison and the London-limit reduction are correct, the work supplies a cutoff-free microscopic route to Bc1 in dirty thin-film strips, which is of direct practical interest for vortex-free operation of superconducting nanodevices. Explicit credit is due for (i) obtaining Bc1 from an independent Gibbs comparison of two self-consistent Usadel solutions rather than by fitting a core cutoff to the same data, (ii) working at arbitrary T < Tc within the dirty-limit framework, and (iii) showing that a single width-independent London cutoff is insufficient. Those strengths would make the paper a useful reference for device-oriented estimates of Bc1 once the free-energy implementation is fully documented and checked.","major_comments":[{"comment":"The central claim that the vortex self-energy exceeds the Pearl–London estimate rests on the Gibbs free-energy difference between self-consistent Usadel solutions. The manuscript must exhibit the explicit free-energy (or Gibbs) functional that is evaluated numerically—Matsubara sum over the anomalous Green’s function plus the electromagnetic contribution—and must demonstrate that this functional reduces to the standard Pearl electromagnetic energy when the order parameter is taken rigid outside a core and the thin-film nonlocal kernel is used. Without that reduction check (or an equivalent analytic limit), the reported excess self-energy and the single-cutoff claim could be an artifact of double-counting kinetic energy, mishandling the applied-field interaction at fixed B, or an incorrect Pearl kernel under the 2-D reduction. This is load-bearing for the abstract’s quantitative statement","section":"Theory / free-energy evaluation (Gibbs difference used for Bc1)"},{"comment":"Numerical reliability of the 2-D Usadel solutions is not adequately established for a result that hinges on small free-energy differences. Mesh resolution near the core and edges, convergence with Matsubara cutoff and spatial discretization, and independence of the Meissner–vortex energy crossing on solver tolerances should be documented (tables or appendices). Residual risks on mesh/convergence and free-energy functional details leave the quantitative excess over Pearl–London only moderately secure.","section":"Methods / numerical solution of the 2-D Usadel problem"},{"comment":"The thin-film reduction to a 2-D in-plane Usadel problem (uniform order parameter across thickness, dirty-limit diffusivity, Pearl-type nonlocal kernel) is assumed to capture the core and edge deformation that set the Gibbs difference. The manuscript should state the range of thickness d, mean free path, and strip width W for which this reduction remains controlled, and should indicate how Bc1 would shift if nonlocal electrodynamics beyond the Pearl kernel or thickness-averaged order-parameter variations become important. This bounds the domain of the claimed excess self-energy.","section":"Formulation of the 2-D Usadel problem and thin-film reduction"}],"minor_comments":[{"comment":"Define the precise subtraction that isolates the vortex self-energy from the total Gibbs difference (Meissner reference, applied-field interaction term, and any constant offsets) so that readers can reproduce the comparison to Pearl–London.","section":"Results / vortex self-energy extraction"},{"comment":"State material inputs (Tc, diffusivity or elastic scattering rate, normal-state conductivity / Pearl length, strip geometry W and d) in a single table or paragraph for each figure that reports Bc1 or self-energy.","section":"Figures and parameter lists"},{"comment":"Clarify notation for the gauge-invariant momentum and the boundary conditions at the strip edges (no-current / vacuum) when the applied field is included directly in the Usadel equations.","section":"Theory setup"},{"comment":"A short comparison plot or table of Bc1 versus the Pearl–London form with several conventional cutoffs (e.g., ξ, ξ/2) would make the ‘cannot be represented by a single width-independent cutoff’ claim immediately readable.","section":"Results / comparison to Pearl–London"},{"comment":"Minor typographical and reference-list consistency checks (Usadel / Pearl–London spelling, arXiv or journal citations for the free-energy functional used) would improve presentation.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The reader’s and skeptic’s concern about the free-energy functional’s London-limit reduction is the decisive load-bearing issue; if the authors supply the explicit functional and a clean reduction check, the paper is likely suitable after revision. Scope fits cond-mat.supr-con. No novelty or citation-pattern concerns beyond ordinary self-citation of prior Usadel/strip work."},"author_rebuttal":{"model":"grok-4.5","summary":"We thank the referee for a careful and constructive report. The three major comments correctly identify load-bearing elements of the work: the Gibbs free-energy functional and its London/Pearl reduction, numerical reliability of the 2-D Usadel solutions, and the controlled domain of the thin-film reduction. We address each point below. Where the manuscript was incomplete we will revise; where the existing formulation already answers the concern we explain why. No standing objections remain after these revisions.","responses":[{"response":"We agree that the free-energy functional and its London/Pearl reduction are load-bearing and were not exhibited with sufficient explicitness. In the revised manuscript we will (i) write out the full Gibbs functional that is evaluated numerically: the standard dirty-limit Usadel free-energy density (Matsubara sum involving the anomalous Green’s function F and the self-consistent gap Δ, together with the condensation-energy term) plus the electromagnetic contribution expressed through the thin-film Pearl kernel for the vector potential at fixed applied field B; and (ii) demonstrate analytically and numerically that, when the order parameter is held rigid outside a prescribed core and the same Pearl nonlocal kernel is used, the functional reduces to the standard Pearl electromagnetic energy (plus the usual core contribution that is absorbed into the London cutoff). This reduction check rules out double-counting of kinetic energy and mishandling of the applied-field interaction. With that documentation in place, the reported excess self-energy relative to the naive Pearl–London estimate and the failure of a single width-independent cutoff remain genuine microscopic results rather than numerical artifacts.","revision_made":"yes","referee_comment":"[Theory / free-energy evaluation (Gibbs difference used for Bc1)] The central claim that the vortex self-energy exceeds the Pearl–London estimate rests on the Gibbs free-energy difference between self-consistent Usadel solutions. The manuscript must exhibit the explicit free-energy (or Gibbs) functional that is evaluated numerically—Matsubara sum over the anomalous Green’s function plus the electromagnetic contribution—and must demonstrate that this functional reduces to the standard Pearl electromagnetic energy when the order parameter is taken rigid outside a core and the thin-film nonlocal kernel is used. Without that reduction check (or an equivalent analytic limit), the reported excess self-energy and the single-cutoff claim could be an artifact of double-counting kinetic energy, mishandling the applied-field interaction at fixed B, or an incorrect Pearl kernel under the 2-D reducti"},{"response":"We accept this criticism. The present manuscript does not document mesh resolution, Matsubara cutoff, spatial discretization, or solver-tolerance independence at the level required for a result that rests on small Gibbs differences. In revision we will add an appendix (or supplementary tables) that reports: (a) local mesh refinement near the vortex core and the strip edges, with measured residual of the Usadel equations; (b) convergence of both the Meissner and single-vortex free energies (and of their crossing field Bc1) with Matsubara cutoff Nω and with successive spatial refinements; and (c) independence of the Meissner–vortex energy crossing on the nonlinear solver tolerances within the range used for production runs. These checks will place the quantitative excess over Pearl–London on a firmer numerical footing.","revision_made":"yes","referee_comment":"[Methods / numerical solution of the 2-D Usadel problem] Numerical reliability of the 2-D Usadel solutions is not adequately established for a result that hinges on small free-energy differences. Mesh resolution near the core and edges, convergence with Matsubara cutoff and spatial discretization, and independence of the Meissner–vortex energy crossing on solver tolerances should be documented (tables or appendices). Residual risks on mesh/convergence and free-energy functional details leave the quantitative excess over Pearl–London only moderately secure."},{"response":"We agree that the domain of validity of the 2-D reduction should be stated explicitly. The reduction assumes (i) dirty-limit diffusivity (ℓ ≪ ξ0), (ii) film thickness d small enough that the order parameter and the Green’s functions are essentially uniform across the thickness (d ≲ ξ(T) and d ≪ λ), and (iii) the standard Pearl nonlocal kernel for the in-plane electrodynamics (valid when d ≪ λ and the strip is treated as infinitesimally thin for the magnetic field). Under these conditions the core and edge deformations that determine the Gibbs difference are captured by the 2-D Usadel problem. In the revised text we will list the corresponding ranges of d, ℓ, and W (relative to ξ and λ) and will note that if nonlocal electrodynamics beyond the Pearl kernel (e.g., finite-thickness corrections to the kernel) or thickness-averaged order-parameter variations become important, Bc1 would receive additional corrections of relative order d/λ or d/ξ; those corrections lie outside the present framework and would systematically shift the absolute scale of Bc1 while leaving the qualitative conclusion—that a single width-independent London cutoff is insufficient—intact within the stated domain.","revision_made":"yes","referee_comment":"[Formulation of the 2-D Usadel problem and thin-film reduction] The thin-film reduction to a 2-D in-plane Usadel problem (uniform order parameter across thickness, dirty-limit diffusivity, Pearl-type nonlocal kernel) is assumed to capture the core and edge deformation that set the Gibbs difference. The manuscript should state the range of thickness d, mean free path, and strip width W for which this reduction remains controlled, and should indicate how Bc1 would shift if nonlocal electrodynamics beyond the Pearl kernel or thickness-averaged order-parameter variations become important. This bounds the domain of the claimed excess self-energy."}],"tokens_in":1932,"tokens_out":1227,"duration_ms":15669,"standing_objections":[]},"desk_editor":{"model":"grok-4.5","letter":"The one thing to know: this paper replaces the Pearl–London Bc1 estimate (which needs a phenomenological core cutoff) with a self-consistent 2D Usadel calculation for a dirty thin-film strip. Vortex and Meissner solutions at fixed applied field give Bc1 from their Gibbs difference, so the core and edge deformation are resolved without a cutoff. That is the actual new result.\n\nWhat it does well is straightforward. The route is standard and defensible: Usadel in the film plane, applied field in the gauge-invariant momentum, free-energy comparison of two self-consistent states. Circularity is low—Bc1 is not fitted via a core radius. Material inputs are external (Tc, diffusivity, geometry). It claims the vortex self-energy comes out larger than the naive Pearl–London value and cannot in general be written as a London log with one width-independent cutoff. If that holds, it is useful for nanostrip devices (SNSPDs, KIDs) where the thermodynamic scale for vortex entry actually matters. The formulation is meant to work at any T < Tc.\n\nThe soft spot is real but proportionate. The excess self-energy claim lives or dies on whether the Usadel free-energy functional (Matsubara sum plus electromagnetic piece) correctly recovers the Pearl electromagnetic energy when the order parameter is rigid outside a core. If the magnetic term at fixed field, the kinetic-energy counting, or the 2D Pearl kernel is mishandled, the reported excess is an artifact. The abstract states the Gibbs difference and the excess without exhibiting the functional or a London-limit check; that reduction is the least secure link. The dirty-limit / thin-film reduction itself is the other assumption: if thickness or mean-free-path corrections matter, the numbers move. Neither issue is a load-bearing collapse of the method; both are checkable in the methods and numerics.\n\nThis is for people who design or model superconducting nanostrips and need a microscopic Bc1 rather than a cutoff-tuned London log. It deserves a serious referee. I would send it to peer review: the calculation is a legitimate extension with clear device relevance, and the free-energy reduction is something a referee can demand and verify. Not a desk reject.","headline":"Cutoff-free microscopic Bc1 for dirty thin-film strips via 2D Usadel + Gibbs comparison; method is sound, free-energy reduction to Pearl–London is the check that matters.","tokens_in":3088,"tokens_out":567,"would_cite":false,"duration_ms":17304,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["74.25.Op","74.78.-w","74.20.Fg","74.25.Ha"],"model":"grok-4.5","headline":"A microscopic theory of Bc1 for dirty superconducting thin-film strips is obtained by solving the 2D Usadel equations without a phenomenological core cutoff.","keywords":["lower critical field","Bc1","Usadel equations","superconducting thin-film strip","vortex self-energy","Pearl–London theory","dirty limit","Gibbs free energy"],"falsifier":"Measure Bc1 of a dirty superconducting nanostrip of known width and thickness at fixed temperature and compare the measured value (and its width dependence) against the Usadel Gibbs-energy prediction; systematic deviation that cannot be absorbed by material parameters would falsify the reduction.","tokens_in":2966,"feed_emoji":"🧲","tokens_out":889,"duration_ms":15944,"temperature":0.7,"pith_summary":"The lower critical field Bc1 of a narrow superconducting thin-film strip is the thermodynamic field at which a vortex first becomes favorable, setting the scale for vortex-free operation in a perpendicular field. Standard Pearl–London estimates leave the vortex core unresolved and therefore require an ad-hoc cutoff whose value is uncertain. This paper replaces that estimate with a microscopic calculation for a dirty strip: the two-dimensional Usadel equations are solved self-consistently in the film plane with the applied field included directly in the gauge-invariant momentum. Vortex and Meissner solutions are obtained at fixed field, and Bc1 is read from the crossing of their Gibbs free energies. Because the core and its deformation by the strip edges are resolved, the resulting vortex self-energy is larger than the naive Pearl–London value and cannot be written as a London logarithm with a single width-independent cutoff. The framework works at any temperature below Tc and supplies a parameter-light basis for predicting Bc1 in nanostrips and related thin-film devices.","feed_headline":"Microscopic Bc1 for dirty thin-film strips, no core cutoff","feed_subtitle":"2D Usadel solutions give a larger vortex self-energy than Pearl–London estimates allow","key_machinery":"Self-consistent numerical solutions of the 2D Usadel equations in the film plane, with the applied field entering the gauge-invariant momentum; the Gibbs free-energy difference between the vortex and Meissner solutions at fixed field supplies Bc1 and resolves the core and edge deformation.","core_discovery":"By solving the two-dimensional Usadel equations self-consistently for a dirty superconducting strip, the lower critical field Bc1 is obtained from the Gibbs-energy difference between vortex and Meissner states without introducing a phenomenological core cutoff; the vortex self-energy exceeds the naive Pearl–London estimate and cannot in general be represented by a London logarithm with a single width-independent cutoff.","pith_inferences":["The same 2D Usadel Gibbs construction could be extended to multi-vortex or edge-barrier configurations to map the full low-field phase diagram of a strip.","If the strip is only moderately dirty, a quasiclassical Eilenberger treatment would be needed to check how much the self-energy shifts away from the Usadel result.","Quantitative comparison with existing Bc1 data on Nb or NbN nanostrips would immediately test whether the excess self-energy is observed.","The method supplies a microscopic route to the effective Pearl length and core size that enter circuit models of superconducting nanowire detectors and resonators."],"forward_implications":["Bc1 of dirty nanostrips can be predicted from material parameters and geometry without an adjustable core cutoff.","The vortex self-energy is systematically larger than Pearl–London estimates that use a conventional cutoff of order the coherence length.","A single width-independent London cutoff cannot reproduce the microscopic self-energy across different strip widths.","The same Usadel framework yields Bc1 at any temperature below Tc, including near Tc where the core size diverges.","Design of vortex-free superconducting thin-film devices can use the computed Gibbs crossing rather than phenomenological estimates."],"fun_headline_variants":["2D Usadel yields larger vortex self-energy without core cutoff","Microscopic Bc1 for dirty strips exceeds Pearl-London estimate","Usadel equations set strip Bc1 from Gibbs energy free of cutoff","Vortex self-energy in dirty strips exceeds single-cutoff London form","Self-consistent Usadel resolves strip vortex core without cutoff"],"cache_read_input_tokens":128,"weakest_assumption_plain":"The dirty-limit Usadel description reduced to two dimensions (thin-film limit, order parameter uniform across thickness, and nonlocal electrodynamics captured only by the Pearl kernel) is assumed to capture the core and edge physics that set the energy difference.","fun_headline_variants_meta":{"raw":{"variants":["2D Usadel yields larger vortex self-energy without core cutoff","Microscopic Bc1 for dirty strips exceeds Pearl-London estimate","Usadel equations set strip Bc1 from Gibbs energy free of cutoff","Vortex self-energy in dirty strips exceeds single-cutoff London form","Self-consistent Usadel resolves strip vortex core without cutoff"]},"model":"grok-4.5","cost_usd":0.016166,"raw_usage":{"total_tokens":3251,"prompt_tokens":739,"num_sources_used":0,"completion_tokens":91,"cost_in_usd_ticks":161660000,"prompt_tokens_details":{"text_tokens":739,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":2421,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":739,"tokens_out":91,"duration_ms":24641,"temperature":1.0,"reasoning_tokens":2421,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-08T19:32:20.232883+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"Measure Bc1 of a dirty superconducting nanostrip of known width and thickness at fixed temperature and compare the measured value (and its width dependence) against the Usadel Gibbs-energy prediction; systematic deviation that cannot be absorbed by material parameters would falsify the reduction.","supporting_citations":[],"review_version":1}