{"id":"a965cfc1-b0bc-4bfc-9529-0722c56f126f","arxiv_id":"2501.14563","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"The net edge current and orbital magnetic moment in mesoscopic chiral d-wave superconductors can be engineered and strongly enhanced by choosing pentagon, hexagon, or disk shapes, with predicted magnetic fields up to 0.01-0.5 mT.","lead":"The authors show that the shape of a small chiral d-wave superconductor can strongly change the size and even direction of its spontaneous edge currents, with pentagons, hexagons, and disks giving the largest magnetic signals. This suggests that patterning mesoscopic samples is a practical way to experimentally detect chiral d-wave superconductivity, which has so far lacked clear fingerprints.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Surface roughness could flip or wash out the shape-dependent edge-current sign; the paper's central design claim needs a roughness/anisotropy sensitivity test.","rationale":"The reader's weakest assumption matches the load-bearing risk I identify. The numerical Eilenberger results are self-consistent, with stated convergence parameters and public data, so the internal calculations are not in question. What is in question is whether the shape-dependent enhancement is a property of a real mesoscopic sample, since the mechanism relies on ballistic trajectories reflecting specularly off multiple edges. Published studies cited by the authors (Refs. [107,108,116,117]) show that surface roughness is not a small perturbation for chiral d-wave edge currents: it can suppress or reverse them even in the semi-infinite limit. Since the finite-size edge-edge interference is an additional, more delicate contribution on top of these known roughness effects, the central experimental predictions inherit this fragility. The paper's own disclaimer in Sec. IIC and the future-work paragraph in Sec. V make the limitation explicit, which is honest, but it means the headline claims should not be read as robust predictions without a sensitivity analysis. My recommended test, adding diffuse scattering or corrugation to the same framework, directly probes whether the sign and magnitude of the currents survive. I also note the minor triangle-geometry classification error in Sec. IIIA (the triangle is called 'concave'), but this does not affect the central argument. The reader's CONDITIONAL verdict is appropriate; no verdict change is needed.","tokens_in":30020,"tokens_out":11859,"duration_ms":129200,"concrete_test":"Using the same SuperConga Eilenberger solver, implement a mixed specular/diffuse boundary condition (fraction p of diffuse reflections) or a corrugated boundary of amplitude ~0.1 xi_0, and recompute the net current I (Eq. 17), orbital magnetic moment m_z (Eq. 8), and flux Phi_ind (Eq. 9) for pentagon, hexagon, and disk at R = 5-10 xi_0, T = 0.1 T_c. The concern is settled if the relative sign of I between disk and pentagon, and the >50% enhancement over the semi-infinite value, survive for p = 0.1 diffuse scattering (or corrugation amplitude 0.1 xi_0). If the sign flips or the enhancement drops below ~50%, the design claim must be substantially qualified. Compare with the known semi-infinite rough-surface current inversion of Refs. [107,108] to calibrate the roughness model.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that sample shape controls the sign and magnitude of the net chiral edge current through edge-edge interference of ballistic quasiparticle trajectories (Sec. IIIA, Fig. 3). This mechanism presumes specular superconductor-vacuum surfaces and a circular Fermi surface (Sec. IIC). The authors explicitly defer disorder, surface roughness, and normal-state anisotropy to future work (Sec. IIC; Conclusions). That deferral is load-bearing because published work shows that surface roughness alone can suppress or invert chiral d-wave edge currents even in semi-infinite geometries without any edge-edge interference: Refs. [107,108,116,117] include spontaneous-current suppression and current inversion induced by rough surfaces and flat-band Andreev bound states. If diffuse scattering randomizes the phase of reflected quasiparticles, the interference that produces the shape-dependent sign (pentagon/hexagon negative vs disk positive for nu=-2) will be weakened, shifted, or reversed. The headline experimental figures (0.01-0.5 mT fields, mu_B/2 per Cooper pair, pentagon-vs-disk sign reversal) are therefore not yet supported for realistic etched or grown interfaces. A quantitative sensitivity test is needed to know whether the design principle survives at all; absent that test, the central experimental prediction is conditional.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper uses self-consistent quasiclassical Eilenberger theory, implemented in the open-source package SuperConga, to study chiral d-wave superconductors in finite two-dimensional geometries: regular polygons and circular disks. The central results are that, for system sizes of tens of coherence lengths, the net edge current, orbital magnetic moment, and induced flux develop a strong and systematic shape dependence: pentagons and hexagons develop large currents opposite to the bulk chirality, circular disks develop large currents along the chirality, while squares and triangles have much smaller currents. The authors attribute this behavior to edge-edge interference between ballistic quasiparticle trajectories and support it with spatial current-density profiles and net-current, OMM, and flux data. They also report that geometries which break the degeneracy between the two d-wave components, such as squares, exhibit an additional second-order transition from a chiral to a nodal state, with a heat-capacity jump of order 10% of the bulk transition jump. The estimated magnetic fields are 0.01–0.5 mT and the magnetic moment approaches μB/2 per Cooper pair, and the paper concludes that mesoscopic patterning is a viable route to detect chiral d-wave superconductivity.","tokens_in":30163,"tokens_out":7340,"duration_ms":72074,"significance":"If the predictions survive contact with realistic materials, this is a significant contribution to the longstanding problem of detecting chiral d-wave superconductivity, because semi-infinite chiral d-wave edges are predicted in many models to have almost vanishing net currents. The study is carefully executed: the numerical solutions are fully self-consistent with explicit convergence criteria (tolerance 10^-7, 20 points per coherence length, 256 Fermi-surface points), the code and data are public, and the known vanishing-current limit is reproduced. The predictions are falsifiable and quantitatively specific, including a shape-dependent sign reversal and a calorimetric signature. The main caveat, which the authors acknowledge, is that the load-bearing model assumptions of a perfectly cylindrical Fermi surface and specular superconductor-vacuum surfaces have not yet been tested for robustness against surface roughness and Fermi-surface anisotropy; these are known to affect chiral d-wave edge currents, and the manuscript's design claim is therefore currently conditional on the idealized model.","major_comments":[{"comment":"The model assumes a perfectly cylindrical Fermi surface and specular superconductor-vacuum surfaces (Sec. IIC), and the Conclusions explicitly defer surface roughness (Refs. [108,116,117]) and normal-state anisotropy (Ref. [203]) to future work. This is a load-bearing simplification because the shape-dependent sign reversal is attributed in Sec. IIIA to ballistic edge-edge interference along straight quasiparticle trajectories; diffuse surface scattering that randomizes the reflected quasiparticle phase will directly affect that interference. Published results for semi-infinite chiral d-wave edges show that surface roughness can suppress or even invert the spontaneous edge current, so the predicted pentagon/hexagon-versus-disk sign reversal and the 0.01–0.5 mT field range are not yet supported for realistic etched or grown interfaces. I request a quantitative sensitivity test (e.g., varying the surface specularity or adding a disorder self-energy in the same self-consistent scheme) for at least the pentagon, hexagon, and disk geometries, or a clearly stated bound on the domain of validity.","section":"Sec. IIC; Conclusions"},{"comment":"The central mechanism of destructive versus constructive edge-edge interference is inferred from the geometric trend in the current-density profiles rather than demonstrated directly. Because Eq. (6) yields only the total current density, the decomposition into bulk-edge and edge-edge contributions shown schematically in Fig. 3 is not verified in the numerics. A direct check—for example, extracting the edge-edge part of the surface propagator, or varying the healing length r0 at fixed shape—would convert this explanation from a qualitative association into a testable design principle. This matters because the sign change between C5/C6 and C∞ is the headline prediction; if the interference interpretation is incorrect, the suggested extrapolation to other shapes is not justified.","section":"Sec. IIIA; Figs. 3, 9, 10"}],"minor_comments":[{"comment":"The statement that 'the triangular system is the only concave geometry' is incorrect; a regular triangle is convex. The intended distinguishing feature is likely the acute internal angle or the sharpness of the corners, and the wording should be corrected.","section":"Sec. IIIA"},{"comment":"The sentence 'Squares, on the other hand, have the the lowest currents' contains a duplicated article 'the the' and should be edited.","section":"Conclusions"},{"comment":"The abstract promises a 'full phase diagram as a function of temperature and system size for different geometries', but the main text presents the phase diagram only for square and disk geometries; other shapes are relegated to Appendix C. Either move part of the Appendix C discussion into the main text or soften the wording in the abstract.","section":"Abstract; Sec. IV"},{"comment":"The estimate of 0.01–0.5 mT assumes a coherence length of 2–10 nm and R ~ 5–10 ξ0, but no explicit layer thickness or stray-field calculation is given. A brief statement of the assumed experimental geometry would make the claimed measurability more concrete.","section":"Sec. IIIB"},{"comment":"The sentence 'The only assumption we make is to assume the same pairing strength in both d-wave channels' is too broad, since the model also assumes a cylindrical Fermi surface and specular boundaries. The sentence should be rephrased to refer specifically to assumptions about the pairing channels.","section":"Sec. IIB"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is technically sound and the self-consistent numerics are credible; the public code and data are a clear strength. My main concern is that the experimental-viability claims outrun the idealized model, and the roughness/anisotropy sensitivity test should be part of the revision rather than deferred. The incorrect 'concave triangle' wording is a minor distraction that should be fixed. This is a suitable paper for cond-mat.supr-con once the robustness question is addressed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Holmvall and Black-Schaffer have done something genuinely useful here. This is the natural follow-up to their 2023 PRB on size-dependent enhancement, and it delivers: a systematic shape dependence from triangles to disks, a sign reversal between low-symmetry polygons and disks, and a new heat-capacity jump marking a chiral-to-nodal transition in squares. The numerics are as solid as it gets for Eilenberger: fully self-consistent, stated tolerances (10^-7), 20 points per coherence length, 256 Fermi-surface points, reproducible with the public SuperConga code and Zenodo data. They also reproduce the known vanishing net current in the semi-infinite limit, which is a good consistency check. Credit where due: the shape-dependence itself, the pentagon/hexagon versus disk sign reversal, and the calorimetric signature are new relative to their prior work.\n\nSoft spots, in proportion. The triangle geometry classification is simply wrong: a regular triangle is convex, not concave. That is minor and easily fixed, but it should be corrected. The larger issue is the clean-limit idealization. The entire mechanism is edge-edge interference along ballistic specular trajectories on a circular Fermi surface. The authors explicitly defer roughness, disorder, and anisotropy to future work, and they cite the relevant literature showing that roughness alone can suppress or invert chiral d-wave currents even without finite-size effects. That does not kill the paper — it is a legitimate clean-limit calculation, and the limitations are stated honestly — but it does mean the headline 0.01–0.5 mT fields and the sign design rule are conditional on surfaces that are specular on the scale of the coherence length. Real etched or grown samples may behave differently. I would ask for a sensitivity test or at least a clearer statement in the abstract that the estimates are clean-limit and material-dependent.\n\nAlso minor: the 0.01–0.5 mT estimate depends on coherence length and is not a universal number; it should be labeled as such in the abstract.\n\nAll that said, this is a careful paper that deserves serious refereeing. The reader's condition is fair: fix the triangle error, add or reference a roughness/anisotropy sensitivity analysis, and qualify the field estimate. The central argument holds up within its stated assumptions. I would recommend sending it to review rather than desk rejecting, with the expectation that the authors address those points.","headline":"A numerically solid, honest mapping of how sample shape controls chiral edge currents in mesoscopic d-wave superconductors, with the main caveat that the headline fields and sign reversal rest on clean, specular surfaces and a circular Fermi surface.","tokens_in":30771,"tokens_out":2311,"would_cite":true,"duration_ms":21069,"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":"Mesoscopic patterning of chiral d-wave superconductors into pentagons, hexagons, or disks can amplify otherwise vanishing edge currents into magnetic fields of 0.01–0.5 mT, with a magnetic moment up to $\\mu_B/2$ per Cooper pair.","keywords":["chiral d-wave superconductivity","mesoscopic patterning","edge currents","orbital magnetic moment","edge-edge interference","heat capacity jump","scanning SQUID magnetometry","quasiclassical theory"],"falsifier":"A scanning SQUID or cantilever magnetometry measurement on patterned mesoscopic films of a candidate chiral d-wave superconductor, with pentagon-, hexagon-, and disk-shaped islands of size tens of $\\xi_0$, looking for spontaneous magnetic fields of order 0.01–0.5 mT whose sign depends on shape as predicted; equivalently, a numerical simulation that adds surface roughness or an anisotropic Fermi surface to the same model and checks whether the shape-dependent sign and the $1/R$ divergence survive.","tokens_in":29742,"feed_emoji":"🧲","tokens_out":6862,"duration_ms":60042,"temperature":0.7,"pith_summary":"Chiral d-wave superconductivity is hard to verify experimentally because the spontaneous chiral edge currents that should signal it are predicted to nearly vanish in bulk samples. This paper argues that the shape of a mesoscopic sample can be used as a design knob: in finite-sized superconductors, edge-edge interference between ballistic quasiparticle trajectories strongly enhances or suppresses the net edge current depending on geometry. Pentagons and hexagons (low rotation symmetry) produce large currents opposite to the bulk chirality, while circular disks produce comparably large currents along the chirality, with estimated magnetic fields of 0.01–0.5 mT and an orbital magnetic moment approaching $\\mu_B/2$ per Cooper pair. The authors also find an additional heat capacity jump in square-shaped samples, marking a chiral-to-nodal transition, which could serve as an indirect calorimetric signature.","feed_headline":"Shaped superconductors can amplify edge currents to 0.5 mT","feed_subtitle":"Pentagons, hexagons, and disks produce the strongest magnetic signal, offering a practical route to detect chiral d-wave superconductivity.","key_machinery":"The mechanism is edge-edge interference encoded in the quasiclassical surface propagator $\\hat{g}(\\mathbf{p}_F,\\mathbf{R};\\varepsilon)$: in a finite sample, straight ballistic trajectories along the Fermi velocity connect different portions of the edge that lie within a correlation length $r_0 \\sim 2$–$10\\xi_0$ of each other. The relative orientation of connected edges decides the sign of the interference—edges meeting at more perpendicular angles (low rotation symmetry) give destructive interference against the naive edge-mode current, while more parallel orientations (approaching the circular disk) give constructive interference. The calculations solve the Eilenberger equation with full self-consistency for both the superconducting order parameter and the vector potential, using a numerically stable Riccati formulation, and obtain the current, magnetic moment, free energy, and heat capacity from the resulting propagators.","core_discovery":"The central claim is that mesoscopic finite-size effects, far from being a nuisance, can be engineered to greatly enhance the spontaneous chiral edge current and its magnetic fingerprints in chiral $d$-wave superconductors. Working within self-consistent quasiclassical theory with a cylindrical Fermi surface and specular edges, the authors find that the net current and the associated orbital magnetic moment and induced flux grow as the system size shrinks, with the net current scaling roughly as $I \\sim 1/R$ until it is cut off by strong suppression of the superconducting order parameter below roughly $D \\sim 10$–$20\\xi_0$. The sign and magnitude depend systematically on the sample's rotation symmetry: for negative bulk chirality, pentagons ($C_5$) and hexagons ($C_6$) develop a large net current antiparallel to the chirality, the circular disk develops a comparably large current parallel to the chirality, and squares and triangles host the smallest currents. They estimate fields of 0.01–0.5 mT and magnetic moments up to $\\mu_B/2$ per Cooper pair, within reach of scanning-probe and cantilever magnetometry. In squares, the geometric suppression of the $d_{xy}$ component creates an intermediate nodal state and a second-order transition to the chiral state at $T^*(D)$, with a heat capacity jump as large as 10% of the bulk normal-superconducting jump.","pith_inferences":["The same edge-edge interference mechanism should operate in other chiral paired states, including chiral $p$-wave superconductors and superfluids, so shape-dependence could be a general route to amplifying edge currents wherever the semi-infinite net current is suppressed.","A practical diagnostic: fabricating several shapes on the same film and comparing the sign and magnitude of their magnetic signals would isolate the geometric interference effect from material-specific parameters, since the ratios between shapes are the geometric prediction.","If surface roughness or Fermi-surface anisotropy reverses the predicted sign for one shape, the shape-dependence could become a sensitive probe of edge quality; conversely, observing the predicted $C_n$-dependent sign could certify a clean ballistic edge.","The predicted $1/R$ divergence suggests that at intermediate sizes the current is dominated by edge-edge trajectories rather than bulk condensate backflow, so a local measurement of the current density near corners versus flat edges could separate the two contributions."],"forward_implications":["Patterning a chiral d-wave film into pentagons, hexagons, or disks of radius $R \\approx 5$–$50\\xi_0$ should produce spontaneous magnetic fields of 0.01–0.5 mT, detectable with current scanning-probe and cantilever magnetometry.","The net edge current and orbital magnetic moment grow roughly as $1/R$ as the sample shrinks, so smaller mesoscopic samples give stronger signals down to the size where order-parameter suppression cuts off the enhancement near $D \\sim 10$–$20\\xi_0$.","The sign of the net current is selected by rotation symmetry: low-symmetry polygons (pentagons, hexagons) yield current opposite to the bulk chirality, while disks yield current along the chirality, providing a shape-controlled test of the mechanism.","In square-shaped samples, an intermediate nodal $d$-wave state appears below a size- and temperature-dependent $T^*(D)$, accompanied by a second-order heat capacity jump up to about 10% of the bulk jump, offering a calorimetric signature.","In the semi-infinite limit all shapes recover the previously established vanishing net current, so the predicted enhancement is inherently a finite-size effect that disappears as $D \\to \\infty$, with the disk converging slowest because of finite edge curvature."],"supporting_citations":[{"why":"The prior mesoscopic finite-size result that this work extends to shape dependence.","marker":"[90]"},{"why":"Establishes surface states, edge currents, and edge-edge contributions to the propagator that drive the shape dependence.","marker":"[88]"},{"why":"Shows vanishing net edge currents in semi-infinite non-p-wave chiral superconductors, the baseline this work aims to overcome.","marker":"[102]"},{"why":"Extends the vanishing-current result and clarifies its non-topological origin, motivating finite-size engineering.","marker":"[103]"},{"why":"Explains the cancellation between chiral edge modes and condensate backflow that the shape effect counteracts.","marker":"[113]"},{"why":"Supplies the self-consistent numerical framework used for all calculated currents, order parameters, and phase diagrams.","marker":"[97]"},{"why":"Provides the comparison scale for net currents in higher-chirality superconductors.","marker":"[110]"},{"why":"Introduces edge-edge interactions in confined chiral superfluids that underpin the interference picture.","marker":"[139]"}],"fun_headline_variants":["Shaped samples amplify chiral d-wave edge currents","Pentagon and hexagon shapes maximize chiral edge currents","Mesoscopic patterning enhances chiral d-wave edge currents","Sample shape controls chiral superconductor edge currents","Low-symmetry shapes sharpen chiral superconductor signals"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is the idealized sample model of Sec. IIC: a perfectly cylindrical Fermi surface, specular superconductor-vacuum surfaces, and no disorder, roughness, or band anisotropy; if real surfaces are rough or the Fermi surface is anisotropic, the predicted edge-edge interference—and with it the sign and magnitude of the currents—could weaken, shift, or reverse.","fun_headline_variants_meta":{"raw":{"variants":["Shaped samples amplify chiral d-wave edge currents","Pentagon and hexagon shapes maximize chiral edge currents","Mesoscopic patterning enhances chiral d-wave edge currents","Sample shape controls chiral superconductor edge currents","Low-symmetry shapes sharpen chiral superconductor signals"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001002,"raw_usage":{"total_tokens":4357,"prompt_tokens":1183,"completion_tokens":3174,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":799,"completion_tokens_details":{"reasoning_tokens":3100}},"tokens_in":799,"tokens_out":3174,"duration_ms":21893,"temperature":1.0,"reasoning_tokens":3100,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T15:01:16.880210+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A scanning SQUID or cantilever magnetometry measurement on patterned mesoscopic films of a candidate chiral d-wave superconductor, with pentagon-, hexagon-, and disk-shaped islands of size tens of $\\xi_0$, looking for spontaneous magnetic fields of order 0.01–0.5 mT whose sign depends on shape as predicted; equivalently, a numerical simulation that adds surface roughness or an anisotropic Fermi surface to the same model and checks whether the shape-dependent sign and the $1/R$ divergence survive.","supporting_citations":[{"cited_title":"Wu and J","cited_arxiv_id":null,"evidence_quote":"Introduces edge-edge interactions in confined chiral superfluids that underpin the interference picture."}],"review_version":1}