{"id":"fba92eea-27bc-494a-9799-6848c2ea13be","arxiv_id":"2412.18754","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"In a simulated s±d wave nematic superconductor, vortices at a twin boundary split into half-quantum fractional vortices identified as core-up and core-down merons, forming skyrmion pairs.","lead":"This paper uses computer simulations of a two-component Ginzburg-Landau model to study vortices in a model of the iron-based superconductor FeSe. It predicts that vortices trapped on a nematic twin boundary split into fractional half-quantum vortices with meron and skyrmion character, a signature that could be seen in vortex-flow experiments.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Half-quantum-vortex identification is inferred from order-parameter weights, not from a computed flux integral; the central observable consequence is therefore unverified.","rationale":"The reader's weakest assumption was that the chosen GL parameter regime is representative of FeSe and that the barrierless twin-boundary model is adequate. That is a legitimate concern, but the more load-bearing issue is internal to the paper's argument: the half-quantum-vortex classification is asserted from the near-equal amplitudes of the s and d components, not from the quantity that physically defines a half-quantum vortex, namely the magnetic flux per fractional vortex. The authors do present direct TDGL solutions showing separated vortex cores and a reduced B(r) peak, so the core-splitting mechanism is credible; credit is also due for checking robustness against ad = 0.8 and Tcd = 0.8 and for a different twin-boundary orientation. However, the paper's own caveat about parameter fine-tuning and the hedged phrase 'can be treated as' signal that the half-flux statement is not a verified output of the simulation. Since the proposed experimental signatures (vortex flow splitting/recombination, flux fraction measurements) are all tied to the magnitude of the fractional flux, this missing flux integral is the single most decisive check. The concern is addressable and does not warrant rejection; the original conditional verdict remains appropriate, so I recommend no change to the reader's verdict.","tokens_in":11296,"tokens_out":3854,"duration_ms":35542,"concrete_test":"From the relaxed H = 0.4 TDGL configuration used for Fig. 2(e), compute the magnetic flux Phi = integral(B_z dA) over a contour around an isolated green (Delta_s-core) and red (Delta_d-core) fractional vortex on the twin boundary, using the same numerical grid, and compare with Phi_0/2. Then repeat with ad = 0.8 and Tcd = 0.8 to see whether Phi deviates from Phi_0/2 as the s/d weights become unequal. If Phi/Phi_0 > 0.6 or < 0.4 in the symmetric case, the 'half-quantum vortex' headline needs revision; if it stays near 0.5 across asymmetric parameters, the claim is robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that twin-boundary vortices are half-quantum vortices rests on the assertion in Section 3: 'Since both order parameters give almost equal weight to the contribution in the present choice of parameters, the fractional vortices can be treated as half-quantum vortices.' The paper shows separated cores and a lower B(r) maximum, but never integrates the local magnetic flux around a green/red fractional vortex. In two-component GL theory, the flux carried by a vortex in one component depends on the superfluid-density ratio, the gradient coefficients Kd and K_tilde, and the coupling terms c_epsilon and gamma2, not on amplitude equality alone. The chosen parameter set (ad = bd = Kd = Tcd = 1, gamma1 = 1.2, gamma2 = 1.0, K_tilde = 1/sqrt(2)) is deliberately symmetric, and the authors themselves note 'fine-tuning of parameters is necessary for quantitative comparison with experimental observation.' Without a flux integral, the observable prediction (half-flux vortices detectable by flow or scanning probes) is not established. The topologically assigned Q = -1/2 per meron is likewise an inference from the n-field texture, not a numerically integrated skyrmion number.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper studies vortex states near a nematic twin boundary in a two-component Ginzburg-Landau model for superconducting FeSe. Using time-dependent GL simulations with a sign-changing Josephson coupling c_epsilon across the boundary, the authors find that the orientation of the nematic vortex core rotates by 90 degrees between s+d and s-d domains and that a vortex placed on the twin boundary splits into two fractional vortices with separated s-wave and d-wave cores. They characterize these objects as core-up and core-down merons, assign a skyrmion number Q = -1 per fractional-vortex pair, and propose that vortex-flow experiments can observe the trapping, splitting, and recombination dynamics.","tokens_in":11495,"tokens_out":4065,"duration_ms":40365,"significance":"If confirmed, the prediction of alternating half-quantum vortices at FeSe twin boundaries would be a concrete and falsifiable manifestation of fractional vortices in a multi-component superconductor, with direct relevance to ongoing STM and scanning-probe experiments. The manuscript's strength is that the vortex splitting and recombination are emergent outputs of a direct numerical solution of the stated TDGL equations, not constructed by inputting the desired vortex texture, and the model inherits the established s+d GL framework. Its weakness is that the half-quantum and meron/skyrmion identifications are inferred from order-parameter weights and local field suppression rather than computed topological or flux integrals, so the quantitative new claims are not yet established.","major_comments":[{"comment":"The central claim that twin-boundary vortices are half-quantum is inferred from the lower internal magnetic field and from the near-equal weights of the two order-parameter components, but no direct flux integral is reported. In the two-component GL functional of Eq. (3), the flux carried by a vortex in one component depends on Kd, K̃, gamma2, and c_epsilon through the coupled Maxwell/GL equations, not on amplitude equality alone. I request a direct evaluation of Φ = ∮ A·dl or ∫ B dS around one red and one green fractional vortex, with the integration contour enclosing only that core, and a comparison with Φ0/2. Without this, the 'half-quantum vortex' label and the proposed experimental signatures are not demonstrated.","section":"§3, Fig. 3(b)"},{"comment":"The assignment Q = -1/2 to each meron is stated on the basis of the sign of the winding of the n-field, but the skyrmion number is not evaluated numerically. The integration region in Eq. (11) is not specified, and a vector-field plot alone does not establish that Q is exactly -1/2 rather than merely having that sign. Please compute Q over the region shown in Fig. 3(e), state how the integration domain is chosen around each isolated vortex, and report the numerical value.","section":"§3, Eq. (11)"},{"comment":"The proposed half-quantum regime relies on nearly equal s-wave and d-wave weights and on a barrierless twin boundary. The manuscript states that 'the barrier potential at the twin-boundary is not considered' and later notes that 'fine-tuning of parameters is necessary for quantitative comparison with experimental observation.' As written, the paper does not show that the chosen parameter set (ad=bd=Kd=Tcd=1, |c_epsilon|=1.0, gamma1=1.2, gamma2=1.0, K̃=1/√2, kappa=5, T=0.1) is representative of FeSe, nor that the vortex splitting survives a moderate s/d imbalance or a finite boundary barrier. I ask for either a robustness scan over the most sensitive coefficients or a clearly scoped statement that the result is a proof of principle for the chosen model rather than a quantitative FeSe prediction.","section":"Parameters around Eq. (3) and §2"}],"minor_comments":[{"comment":"The text defines η† = (Δs*, Δd*), so η†η should equal |Δs|^2 + |Δd|^2, but the manuscript writes |Δs(r)|^2 + |Δd(r)| (missing the square on the second term) and Fig. 3(d) writes n_z = (|Δs|^2 - |Δd|)/η†η (also missing a square). Please correct these expressions.","section":"Eq. (10) and Fig. 3(d)"},{"comment":"The prefactor 1/12 in the TDGL relaxation term of Eq. (4) is not explained; please state its origin or indicate that it is a chosen dimensionless relaxation coefficient.","section":"Numerical method, Eq. (4)"},{"comment":"The manuscript reports only the system size 60×60; for reproducibility it would be helpful to state the grid spacing, time step, and convergence criterion used in the TDGL evolution.","section":"Numerical setup"},{"comment":"The circular arrows indicating phase winding in Fig. 2(d) would be easier to interpret if the color correspondence between red/green fractional vortices and the two winding directions were made explicit in the caption or in the main text.","section":"Fig. 2(d)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the scope of a condensed-matter theory journal and the numerical approach appears sound, but the central quantitative claims—half-quantum flux and exact meron/skyrmion charges—require direct computation before publication. I would not recommend rejection because the missing checks are straightforward extensions of the existing simulation: evaluating the local flux integral around each fractional vortex and numerically integrating the skyrmion density over the region shown in Fig. 3(e)."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"I think the reader's take is about right. The new thing here is concrete: in a two-component s±d GL model with a nematic twin boundary, a vortex right on the boundary splits into two cores, one in each component, and the two fractional vortices form an alternating core-up/core-down meron array with a skyrmion number of -1 per pair. The 90-degree reorientation of the vortex core across the boundary is also clearly demonstrated. The TDGL simulations look internally consistent, and the flow snapshots (trapping and splitting, then recombination on escape) give the experimental signature a concrete form.\n\nCredit where it's due: the paper is clearly written, the figures are informative, and the authors are honest about the model's limitations. They mention the neglected barrier potential, the need for fine-tuning, and they explicitly list flux-quantization measurement as future work. The reference list covers the relevant fractional-vortex and meron literature.\n\nThe soft spot is exactly what the stress-test note identifies. 'Half-quantum' is asserted from the near-equal weights of the two components and a lower B(r) maximum, never from a flux integral through the vortex. That matters: in two-component GL the flux carried depends on gradient terms and couplings, not just amplitude ratios. The skyrmion number is similarly read off the n-field texture rather than computed. These are the load-bearing quantities for the title's claim, and they are inferences. Also, the parameter set is highly symmetric and the twin boundary is idealized; the authors concede quantitative comparison would need fine-tuning. None of this is fatal to the basic phenomenon, but 'realized' in the title overstates what is a numerical prediction.\n\nI'd send it to a referee. The topic is of real interest to the FeSe vortex community and to anyone working on fractional vortices. A referee should ask for a computed flux per vortex, an explicit skyrmion-number integral, and a bit more parameter variation. Those are straightforward requests, and the paper would be stronger with them.","headline":"Solid GL simulation of vortex splitting at a FeSe twin boundary, but the half-quantum label is inferred, not computed.","tokens_in":12111,"tokens_out":2398,"would_cite":true,"duration_ms":34446,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper claims that an ordinary vortex parked on a nematic twin boundary in an s-plus-d superconductor like FeSe splits into two fractional vortices, which behave as half-quantum vortices and form a meron pair with skyrmion number Q =…","keywords":["fractional vortex","half-quantum vortex","nematic superconductor","twin boundary","FeSe","meron","skyrmion","Ginzburg-Landau theory"],"falsifier":"A local probe that resolves sub-flux-quantum magnetic fields, such as scanning SQUID or high-resolution vortex imaging, could settle the claim: the paper predicts that a vortex on a FeSe twin boundary should show a magnetic-field peak roughly half the height of a bulk vortex, two spatially separated fractional cores, and a flow process in which a bulk vortex splits on entering the boundary and recombines on leaving; observing a single unsplit core of full flux at the boundary would refute the claim.","tokens_in":10987,"feed_emoji":"🌀","tokens_out":8937,"duration_ms":80045,"temperature":0.7,"pith_summary":"This paper asks what happens to an ordinary magnetic vortex when it sits exactly on a nematic twin boundary in a superconductor such as FeSe, where the superconducting pairing changes from s+d to s-d across the boundary. Using two-component Ginzburg-Landau theory, it finds that the vortex can no longer keep its s-wave and d-wave components locked together: the two cores separate, producing a pair of fractional vortices that can be treated as half-quantum vortices when the components have nearly equal weight. It then shows that this separated pair is topologically nontrivial, forming a core-up and core-down meron pair with skyrmion number Q = -1, and that the splitting and recombination should be visible in the time evolution of vortex flow along the boundary. If correct, this gives FeSe twin boundaries as a candidate real-material platform for fractional vortices and skyrmion-like textures that have long been sought in multi-component superconductors.","feed_headline":"Twin-boundary vortex splits into two half-quantum pieces","feed_subtitle":"Simulations of FeSe show each trapped vortex becomes two half-quantum vortices, with a flow signature to test.","key_machinery":"The machinery is a two-component Ginzburg-Landau free energy for s-wave and d-wave order parameters, with a Josephson coupling term whose coefficient changes sign across the twin boundary and a cross-gradient term that produces the nematic anisotropy. The sign change forces the real part of the d-wave order parameter to flip while the s-wave component stays roughly constant, creating a region of s-plus-id character with locally broken time-reversal symmetry at the boundary. The object that carries the topological argument is the unit vector n = (n_x, n_y, n_z) built from the two order parameters; the rotation of its in-plane part and the sign change of n_z across the separated cores define the core-up and core-down merons, and integrating n dot (curl n) over a pair gives Q = -1. Time-dependent Ginzburg-Landau equations are then evolved to produce vortex-flow snapshots in which a conventional vortex splits on entry into the boundary and recombines on exit.","core_discovery":"The paper's central claim is that the nematic twin boundary acts as a line where a conventional vortex fractionalizes. In the s+d and s-d domains on either side, the relative phase between the s-wave and d-wave order parameters is locked near 0 or pi, and an ordinary vortex has a single core in both components. On the boundary itself the relative phase rotates to plus or minus pi/2, so a vortex cannot keep both components zero at the same point; instead the s-wave core and the d-wave core separate along the boundary. With the parameter choice used here, the two components contribute almost equally, so each separated core carries roughly half of a flux quantum and can be regarded as a half-quantum vortex. The separated pair is not just a flux split: the map of the unit vector n built from the two order parameters shows a core-up meron and a core-down meron whose combined skyrmion number is Q = -1, and the sequence of vortices along the boundary forms a skyrmion lattice.","pith_inferences":["Beyond the paper: if the s-wave and d-wave weights are tuned to be unequal, the boundary vortices should interpolate continuously between half-quantum and nearly conventional behavior, so the twin boundary could act as a knob for the flux fraction per vortex; the paper's brief runs with asymmetric parameters hint at this but do not quantify it.","Beyond the paper: the same sign-flip mechanism should produce fractional vortices at any domain wall in a two-component superconductor where the Josephson coupling changes sign, so the prediction is not limited to FeSe and could be tested in other nematic or multi-component systems.","Beyond the paper: because each meron pair carries Q = -1, a transport current or field sweep along the boundary may drive collective skyrmion-like motion of the vortex array, analogous to meron lattices in magnets; the paper only demonstrates trapping and escape of single vortices."],"forward_implications":["A conventional vortex crossing the twin boundary splits into two fractional vortices, one centered in the s-wave component and one in the d-wave component, and this splitting is visible in the order-parameter profiles.","Each fractional vortex carries about half of a flux quantum in the studied parameter range, so a vortex on the boundary has a measurably lower local magnetic field than a bulk vortex.","The split pair is a meron-anti-meron texture with skyrmion number Q = -1, and an alternating array of such pairs along the boundary constitutes a skyrmion lattice.","The nematic vortex core orientation rotates by 90 degrees across the twin boundary.","Vortex flow parallel to the boundary shows a repeatable process: a bulk vortex is trapped, splits into two fractional vortices, and recombines when it escapes, giving a time-dependent signature for experiments."],"supporting_citations":[{"why":"Supplies the twin-boundary Ginzburg-Landau model with sign-changing c_epsilon and the s-plus-id state near FeSe twin boundaries that the paper adopts as its starting point.","marker":"[12]"},{"why":"Provides the two-component Ginzburg-Landau treatment of the elliptical vortex and oblique vortex lattice in FeSe, including the nematic cross-gradient term used here.","marker":"[15]"},{"why":"Establishes the concept of vortices carrying fractional flux in multi-component superconductors, the class of objects this paper seeks in FeSe.","marker":"[31]"},{"why":"Supplies the time-dependent Ginzburg-Landau framework for vortex states in multi-component superconductors used to compute the vortex structures.","marker":"[32]"},{"why":"Introduces the unit-vector n and skyrmion-number description of nematic skyrmions in odd-parity superconductors that the paper applies to the fractional vortex pair.","marker":"[36]"},{"why":"Analyzes half-quantum vortices in a nematic superconductor, giving the half-quantum interpretation used for the separated s/d cores.","marker":"[37]"},{"why":"Reports a vortex carrying a temperature-dependent fraction of the flux quantum in an Fe-based superconductor, the experimental precedent for fractional vortices.","marker":"[40]"},{"why":"Supplies the meron and skyrmion topological-texture framework used to classify the core-up and core-down vortices as merons.","marker":"[48]"}],"fun_headline_variants":["Vortex fractionalizes at nematic twin boundary","Half-quantum vortices emerge on twin boundary","FeSe twin boundary splits vortices into merons","Skyrmion lattice from vortex splitting in FeSe","Twin boundary turns vortices into half-quantum pairs"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that real FeSe twin boundaries behave like the idealized, barrier-free boundary in the simulation, with the s-wave and d-wave components contributing almost equally; if the actual parameters sit in a different regime, the vortex splitting into half-quantum pairs could disappear.","fun_headline_variants_meta":{"raw":{"variants":["Vortex fractionalizes at nematic twin boundary","Half-quantum vortices emerge on twin boundary","FeSe twin boundary splits vortices into merons","Skyrmion lattice from vortex splitting in FeSe","Twin boundary turns vortices into half-quantum pairs"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000538,"raw_usage":{"total_tokens":2538,"prompt_tokens":858,"completion_tokens":1680,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":474,"completion_tokens_details":{"reasoning_tokens":1617}},"tokens_in":474,"tokens_out":1680,"duration_ms":10544,"temperature":1.0,"reasoning_tokens":1617,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T04:30:25.729893+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A local probe that resolves sub-flux-quantum magnetic fields, such as scanning SQUID or high-resolution vortex imaging, could settle the claim: the paper predicts that a vortex on a FeSe twin boundary should show a magnetic-field peak roughly half the height of a bulk vortex, two spatially separated fractional cores, and a flow process in which a bulk vortex splits on entering the boundary and recombines on leaving; observing a single unsplit core of full flux at the boundary would refute the claim.","supporting_citations":[{"cited_title":"Song, Y .-L","cited_arxiv_id":null,"evidence_quote":"Supplies the twin-boundary Ginzburg-Landau model with sign-changing c_epsilon and the s-plus-id state near FeSe twin boundaries that the paper adopts as its starting point."},{"cited_title":"Hashimoto, Y","cited_arxiv_id":null,"evidence_quote":"Provides the two-component Ginzburg-Landau treatment of the elliptical vortex and oblique vortex lattice in FeSe, including the nematic cross-gradient term used here."},{"cited_title":"Sigrist and D","cited_arxiv_id":null,"evidence_quote":"Establishes the concept of vortices carrying fractional flux in multi-component superconductors, the class of objects this paper seeks in FeSe."},{"cited_title":"Babaev, V ortices with Fractional Flux in Two-Gap Super- conductors and in Extended Faddeev Model, Phys","cited_arxiv_id":null,"evidence_quote":"Supplies the time-dependent Ginzburg-Landau framework for vortex states in multi-component superconductors used to compute the vortex structures."},{"cited_title":"Tanaka, H","cited_arxiv_id":null,"evidence_quote":"Introduces the unit-vector n and skyrmion-number description of nematic skyrmions in odd-parity superconductors that the paper applies to the fractional vortex pair."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Analyzes half-quantum vortices in a nematic superconductor, giving the half-quantum interpretation used for the separated s/d cores."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reports a vortex carrying a temperature-dependent fraction of the flux quantum in an Fe-based superconductor, the experimental precedent for fractional vortices."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the meron and skyrmion topological-texture framework used to classify the core-up and core-down vortices as merons."}],"review_version":1}