{"id":"d3236353-4ae7-4b9d-b670-2cabd20c6f8f","arxiv_id":"2412.14876","paper_version":2,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Edge phase crystals in weak-coupling d-wave superconductors remain the ground state up to roughly half the critical impurity strength, and a uniform edge-current Vorontsov state competes in mesoscopic squares.","lead":"Phase crystals are ripples in the phase of a superconducting wavefunction that drive spontaneous current loops at sample edges. This paper computes how much nonmagnetic impurity disorder these ripples can survive in a d-wave superconductor, finding they persist up to 40 to 50 percent of the impurity strength that destroys superconductivity.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Real-material robustness depends on specular edges; unquantified surface roughness could shrink the 40–50% impurity window.","rationale":"The reader correctly identifies the specular-surface assumption as the main weakness, and I agree that it is the most load-bearing model assumption for the real-material interpretation of the central claim. However, I think the verdict should be moved from a flat ACCEPT to CONDITIONAL: the numerical phase diagram is likely correct for the stated model, but the paper's conclusion that phase crystals 'should be possible to observe... in weak-coupling superconductors' implicitly assumes that surface disorder does not independently destroy the ABS that seed the phase crystal. Since the paper explicitly leaves the critical roughness unquantified, the real-material claim is not yet fully supported. My proposed check is a single, concrete numerical experiment that would settle whether the 40–50% impurity-robustness window survives realistic surface roughness; until then, conditional acceptance is the more precise verdict.","tokens_in":21154,"tokens_out":6664,"duration_ms":62141,"concrete_test":"In SuperConga, implement Nagato et al. generalized boundary conditions with a roughness/specularity parameter R for the triangle-like geometry, and compute T*(Γ = 0) and T*(Γ = 0.4Γ_c) in both the unitary and Born limits at T = 0.05T_c as R varies. If T*(Γ = 0.4Γ_c) collapses toward zero for R values corresponding to modest roughness (e.g. specularity below ~0.5), the headline robustness window is conditional on surface quality.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The headline claim that the phase crystal survives up to ~40–50% of the critical impurity strength is computed for specular, atomically flat [110] surfaces, as stated in Sec. IIF. The zero-energy Andreev bound states that drive the phase crystal are sharp only under ideal specular reflection; surface roughness and finite interface transparency broaden these states, which the authors acknowledge in Sec. VI: 'Our results imply that phase crystals should be stable against such surface roughness but only until some critical roughness where the broadening of surface ABS is too large.' Ref. [41] already showed that surface ruggedness can suppress spontaneously symmetry-broken surface states. Because surface disorder is independent of the bulk impurity scattering treated here, a real sample's phase-crystal window is the intersection of the impurity-scattering window with a surface-quality window. The paper provides no estimate of the critical roughness; if realistic roughness broadens the ABS on the same scale as impurities at Γ ≈ 0.4Γ_c, the 40–50% robustness number does not transfer to real materials. This does not invalidate the self-consistent numerical phase diagram for ideal surfaces, but it makes 'plausible in real materials' conditional rather than established.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper uses quasiclassical Eilenberger theory with t-matrix impurity self-energies, solved self-consistently in the order parameter, impurity self-energies, and vector potential, to study phase crystals at [110] edges of nodal d-wave superconductors. It quantifies impurity-induced broadening of the surface Andreev bound states, computes Tc(Γ) in single-edge and square geometries against the Abrikosov-Gor'kov benchmark, and maps a T–Γ phase diagram for a range of scattering phase shifts. The central findings are that the phase crystal survives up to roughly 40–50% of the critical impurity strength in the Born and unitary limits, and that in mesoscopic squares edge–edge hybridization stabilizes a distinct translationally invariant TRSB (Vorontsov) state.","tokens_in":21344,"tokens_out":11210,"duration_ms":105674,"significance":"The result is significant because it extends phase-crystal predictions to weak-coupling d-wave superconductors with disorder, which is important for experimental feasibility. The paper's numerical scheme is a strength: full self-consistency, benchmarks against analytic ABS broadening and the Abrikosov-Gor'kov Tc curve, and public open-source code and data. The 40–50% window is an emergent output of a free-energy comparison rather than a fitted target. The remaining concerns are about how the ideal-surface result is worded, not about the internal consistency of the calculation.","major_comments":[],"minor_comments":[{"comment":"The sentence beginning 'We find find that the ABS...' contains a duplicated word and should be corrected to 'We find that...'.","section":"Sec. I"},{"comment":"Because the calculations assume specular, atomically flat pair-breaking edges, the abstract and conclusions should explicitly qualify the robustness claim: surface roughness broadens the ABS and the critical roughness is not quantified here. This is not a flaw in the numerical phase diagram, but the real-material inference would be clearer if the caveat appeared in the abstract or conclusions rather than only in the outlook.","section":"Sec. IIF and Sec. VI"},{"comment":"Please state explicitly whether the 40–50% of the critical impurity strength is quoted for the triangular single-edge geometry or the square, and whether the square's low-temperature TRSB state at the largest Γ is the Vorontsov phase rather than the phase crystal; as written, Fig. 4(b) shows that region as Vorontsov, so the headline claim needs this qualification.","section":"Sec. V and Fig. 4"},{"comment":"'edges rotated45◦ degrees' should read 'edges rotated 45°' or 'edges rotated by 45 degrees'.","section":"Fig. 1 caption"}],"recommendation":"minor_revision","confidential_remarks":"The paper is suitable for this journal and the central calculation appears sound. I do not see a novelty or attribution issue; prior phase-crystal work is properly cited. The revision should mainly tighten the wording of the robustness claim and clarify the phase-crystal versus Vorontsov distinction."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Magnus,\n\nThis is the paper that turns the phase crystal prediction into a quantitative disorder-robustness statement. The central number—phase crystal survives up to ~40–50% of the critical impurity strength in both Born and unitary limits—is computed with full self-consistency in the impurity self-energies, the d-wave order parameter, and the vector potential, and it is an emergent output of a standard Eilenberger free energy, not a fit. The benchmarks against the Abrikosov-Gor'kov Tc curve and analytic ABS broadening check out. Data and the SuperConga code are public. That is real evidence, and it holds up.\n\nWhat is genuinely new is the full Γ–T phase diagram with general scattering phase shifts, and the impurity-driven emergence of the Vorontsov phase in a fully confined square geometry. Both follow directly from the numerics and are presented clearly. Prior phase crystal work from this group is cited appropriately; this paper adds the missing weak-coupling disorder map.\n\nThe soft spot is the one the authors admit in Sec. VI: the calculation assumes specular, atomically flat [110] surfaces. The zero-energy ABS that drive the crystal are sharp only under specular reflection, and roughness or finite interface transparency broadens them. The paper says the crystal should be stable up to a critical roughness but does not quantify it. That means the 40–50% window is for ideal surfaces. If realistic roughness broadens the ABS on the same scale as impurity scattering at Γ≈0.4Γc, the robustness number could shrink. The stress-test note puts this fairly; it is a real limitation, not a manufactured one. I would not overturn the paper over it, because the authors explicitly flag it and the qualitative conclusion—disorder alone does not kill the phase crystal—is well supported. But \"plausible in real materials\" is conditional, not established.\n\nThe other approximations, noncrossing t-matrix and circular Fermi surface, are standard and discussed. Nothing in the citation pattern bothers me.\n\nVerdict: send to peer review. It deserves a serious referee, and the referee's main task is to weigh the roughness question against the ideal-surface phase diagram. I would cite this paper if I worked on phase crystals or edge states in d-wave superconductors.","headline":"A self-consistent numerical phase diagram that makes the phase crystal's disorder robustness credible for ideal surfaces; the specular-edge assumption is the main caveat for real materials.","tokens_in":21903,"tokens_out":2849,"would_cite":true,"duration_ms":21696,"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":"Phase crystals survive impurities up to ~50% of critical strength","keywords":["phase crystal","d-wave superconductor","Andreev bound states","impurity scattering","t-matrix approach","time-reversal symmetry breaking","quasiclassical theory","mesoscopic superconductivity"],"falsifier":"Compute the same phase diagram with diffuse or rough surface scattering instead of specular reflection at the pair-breaking edges. If the phase crystal disappears at impurity strengths below roughly 40% of Γc for realistic roughness, the claimed disorder robustness window would close.","tokens_in":20959,"feed_emoji":"🌀","tokens_out":8520,"duration_ms":54730,"temperature":0.7,"pith_summary":"This paper asks whether the phase crystal—a spontaneously modulated superconducting state with periodic currents and magnetic fields—survives the disorder that is always present in real materials. Using a t-matrix treatment of nonmagnetic impurities inside quasiclassical theory, fully self-consistent in impurity self-energies, the d-wave order parameter, and the vector potential, it computes the ground-state phase diagram as a function of impurity strength and temperature. The central result is that the phase crystal at [110] edges of a nodal d-wave superconductor survives up to roughly 40–50% of the critical impurity strength that destroys superconductivity, across the whole range from weak (Born) to strong (unitary) scattering. A second result is that in mesoscopic squares, impurities can shrink the effective system size and convert the phase crystal into a competing state with uniform edge currents. The upshot is that phase crystals should be observable in realistic weakly disordered d-wave samples, not only in strongly correlated models.","feed_headline":"Phase crystals survive impurities up to ~50% of critical strength","feed_subtitle":"Full self-consistent phase diagram puts the periodic-current state within reach of real samples.","key_machinery":"The t-matrix impurity self-energy inside quasiclassical theory, parameterized by a scattering energy Γ and phase shift δ0 (Born limit δ0→0, unitary limit |δ0|→π/2), solved self-consistently with the d-wave order parameter and the vector potential. The phase crystal is driven by zero-energy Andreev bound states at pair-breaking [110] edges: Doppler-shifting these states to finite energies lowers the free energy, while the condensate backflow costs kinetic energy, and the nonlocal superfluid stiffness tensor couples the crystal period to the bound-state decay length. The Vorontsov phase is a competing time-reversal breaking state with uniform edge currents that appears when edge-edge hybridization becomes strong in finite systems.","core_discovery":"The central claim is that the phase crystal—a spontaneously phase-modulated superconducting ground state with periodic currents and magnetic fields—survives nonmagnetic disorder in weak-coupling nodal d-wave superconductors up to roughly 40–50% of the impurity strength that destroys superconductivity. The authors establish the full ground-state phase diagram in temperature–impurity strength space, solving the quasiclassical Eilenberger equation with t-matrix impurity self-energies, the d-wave order parameter, and the vector potential all self-consistent. Across scattering phase shifts from Born to unitary, the phase crystal transition temperature T* decreases roughly linearly with impurity strength but remains nonzero until Γ* ≈ (0.4–0.5)Γc. In mesoscopic squares with multiple pair-breaking edges, impurity-enhanced coherence length shrinks the effective system size and drives hybridization between edge states, suppressing Tc and eventually converting the phase crystal into the Vorontsov phase with translationally invariant edge currents.","pith_inferences":["Beyond the paper's explicit claims, the same self-consistent t-matrix machinery could be run with diffuse surface boundary conditions to produce a roughness–impurity phase diagram; the paper only argues qualitatively that roughness weakens the crystal, leaving the crossover unquantified.","Because impurity concentration changes the number and size of current loops, disorder could serve as a control knob for the topological defect structure of the superflow, consistent with the generalized Poincaré–Hopf constraint invoked in the paper.","The smooth transition between phase crystal and Vorontsov phase at large Γ suggests that coexisting periodic and uniform edge currents may appear in finite samples, which scanning probes could distinguish.","Extending the same phase-diagram computation to other nodal pairings could test whether a roughly 40–50% robustness window is a generic feature of flat-band edge states or specific to d-wave symmetry."],"forward_implications":["Phase crystals are plausible in real, weakly disordered d-wave superconductors, since disorder at 40–50% of the critical level still leaves them as the ground state.","The phase diagram gives experimentalists a concrete temperature–impurity window in which to look for periodic spontaneous currents and fields with scanning probes.","In mesoscopic squares, increasing impurity concentration can switch the time-reversal breaking state from a phase crystal to the uniform edge-current (Vorontsov) state, so sample size and disorder jointly select the broken-symmetry pattern.","Finite-size effects suppress Tc more strongly than the bulk Abrikosov-Gorkov curve in systems with multiple pair-breaking edges, so interpreting Tc measurements in small grains requires accounting for edge-edge hybridization.","The phase crystal is more robust in the unitary scattering limit than in the Born limit, which can be used to identify the dominant scattering regime from the disorder dependence of time-reversal breaking."],"supporting_citations":[{"why":"Establishes that phase crystals form at pair-breaking edges of d-wave superconductors via flat-band Andreev bound states.","marker":"[5]"},{"why":"Derives the phase crystal instability from negative, inhomogeneous superfluid stiffness and couples crystal period to the bound-state decay length.","marker":"[6]"},{"why":"Provides the free-energy argument that shifting zero-energy Andreev bound states to finite energies lowers the free energy, driving spontaneous time-reversal breaking.","marker":"[17]"},{"why":"Showed phase crystals are disorder-robust when strong electron correlations are included, the comparison case the authors extend to weak coupling.","marker":"[35]"},{"why":"Defines the t-matrix impurity self-energy and pair-breaking strength Γ used to map the phase diagram.","marker":"[45]"},{"why":"Gives the analytic impurity broadening of Andreev bound states that the authors reproduce and build on.","marker":"[49]"},{"why":"Proposed the competing translationally invariant edge-current (Vorontsov) phase in a slab, which the mesoscopic square realizes.","marker":"[54]"},{"why":"Supplies the Abrikosov-Gorkov Tc(Γ) curve used to define the critical impurity strength and benchmark the bulk limit.","marker":"[62]"}],"fun_headline_variants":["Phase crystals survive half the impurity threshold","Phase crystals withstand ~50% of critical impurity strength","Impurity-tolerant phase crystals persist to half the limit","Phase crystals remain stable up to 50% critical impurity"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The 40–50% robustness numbers assume atomically flat, specularly reflecting pair-breaking edges; if real surface roughness broadens the zero-energy Andreev bound states as much as impurity scattering does, the phase crystal could vanish well before that window.","fun_headline_variants_meta":{"raw":{"variants":["Phase crystals survive half the impurity threshold","Phase crystals withstand ~50% of critical impurity strength","Impurity-tolerant phase crystals persist to half the limit","Phase crystals remain stable up to 50% critical impurity"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000224,"raw_usage":{"total_tokens":1467,"prompt_tokens":956,"completion_tokens":511,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":572,"completion_tokens_details":{"reasoning_tokens":450}},"tokens_in":572,"tokens_out":511,"duration_ms":4243,"temperature":1.0,"reasoning_tokens":450,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T11:49:59.884805+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the same phase diagram with diffuse or rough surface scattering instead of specular reflection at the pair-breaking edges. If the phase crystal disappears at impurity strengths below roughly 40% of Γc for realistic roughness, the claimed disorder robustness window would close.","supporting_citations":[{"cited_title":"Chakraborty, T","cited_arxiv_id":null,"evidence_quote":"Showed phase crystals are disorder-robust when strong electron correlations are included, the comparison case the authors extend to weak coupling."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the t-matrix impurity self-energy and pair-breaking strength Γ used to map the phase diagram."},{"cited_title":"Poenicke, Y","cited_arxiv_id":null,"evidence_quote":"Gives the analytic impurity broadening of Andreev bound states that the authors reproduce and build on."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Abrikosov-Gorkov Tc(Γ) curve used to define the critical impurity strength and benchmark the bulk limit."}],"review_version":1}