{"id":"e5d5d422-6d1e-4a9e-a09c-714f7815f96b","arxiv_id":"1908.06215","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":8.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Microscopic calculations show the o vortex is the most stable axisymmetric vortex in 3P2 superfluids under strong magnetic fields and hosts two zero-energy Majorana fermions in its core.","lead":"This paper computes quantized vortex structures in 3P2 superfluids with a microscopic theory (Eilenberger and Bogoliubov-de Gennes equations), going beyond the usual Ginzburg-Landau approach. It finds that in a strong magnetic field the axisymmetric o vortex is most stable and carries two zero-energy Majorana fermions in its core, which could matter for neutron star cooling and glitches.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Stability claim is conditional on the axisymmetric ansatz of Eq. (18); the paper's own open question about nonaxisymmetric vortices is the load-bearing gap for the physical conclusion.","rationale":"The reader's weakest_assumption is exactly right, and my independent reading agrees. The paper is careful to qualify the main stability statement as 'axisymmetric,' and the Majorana part rests on P3 symmetry plus a direct BdG diagonalization, both of which are solid within the stated sector. The conditional verdict is therefore appropriate; my concern does not change that verdict. I considered whether a stronger internal objection exists, such as the non-self-consistent magnetization feedback into the order parameter, but the authors explicitly flag this (end of Sec. III.A) and it does not bear on the topological zero-mode existence, only on quantitative profiles and possibly subtle free-energy shifts. The nonaxisymmetric gap is the single load-bearing limitation because it controls whether the o vortex, rather than some uncomputed vortex, is the physical ground state in a rotating neutron-star core. The paper's speculation in Sec. IV that the symmetric traceless tensor suppresses double-core vortices is reasonable but not a computation; Ref. 50's half-quantized vortex pairs are a concrete nonaxisymmetric alternative. Thus the reader's CONDITIONAL verdict should stand unchanged.","tokens_in":23693,"tokens_out":6652,"duration_ms":76537,"concrete_test":"Perform a two-dimensional Eilenberger + gap-equation calculation on a full (ρ,θ) grid without imposing Eq. (18), at T=0.4Tc and VZ/Tc=0.9, 1.4, 1.5, using as initial conditions (i) the D4-BN-o2 vortex with small quadrupolar e±2iθ perturbations, (ii) a 3He-B-type double-core ansatz adapted to the symmetric traceless 3P2 tensor, and (iii) a pair of half-quantum vortices from Ref. 50. Compare the converged Luttinger-Ward free energies (Eq. (27)) to the axisymmetric values. If any nonaxisymmetric solution converges to lower Jsn, or if the o vortex is linearly unstable to the e±2iθ deformation, the statement that the o vortex is the most stable state under strong fields fails beyond the axisymmetric sector.","verdict_should_be":"UNCHANGED","load_bearing_attack":"All self-consistent free-energy comparisons (Figs. 3 and 6, Sec. III) are confined to order parameters satisfying the axisymmetry condition (J3−κ)A=0, Eq. (18). Within that subspace the o vortex at strong field is the most stable solution, and the two k3=0 zero modes are confirmed numerically (Fig. 7(f)) and by the P3 chiral winding argument, Eqs. (41)-(42). The load-bearing gap is not internal to this calculation but its reach: the physical conclusion that o vortices govern strong-field vortex matter requires that no nonaxisymmetric vortex has lower free energy. The paper does not test this, and explicitly lists nonaxisymmetric vortices as future work (Sec. IV). The analogy is suggestive but not a proof: in 3He-B the nonaxisymmetric double-core vortex remains stable even when axisymmetric v vortices are destabilized by strong fields (Refs. 61, 62, 82), and Ref. 50 predicts half-quantized vortex pairs in the D4-BN phase. The paper's argument that the symmetric traceless tensor of 3P2 may suppress these states is qualitative. Because the zero modes are protected only in the o vortex's P3-symmetric sector, a lower-energy nonaxisymmetric configuration would remove the claimed physical relevance even though the axisymmetric theorem stands.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a microscopic study of singly quantized axisymmetric vortices in a $^{3}P_{2}$ superfluid, modeled by a zero-range spin-triplet $p$-wave interaction with Zeeman coupling. The authors solve the Eilenberger equation and the gap equation self-consistently for several boundary conditions corresponding to UN, D2-BN, and D4-BN uniform phases, classify the resulting vortices as $o$ or $v$ according to whether the magnetic $\\pi$-rotation symmetry $P_{3}$ is preserved, and compare their free energies at $T=0.4T_{c}$ with and without a magnetic field. They then solve the Bogoliubov-de Gennes (BdG) equation on the self-consistent order-parameter profiles to obtain quasiparticle spectra, spin-polarized bound-state branches, and local spin densities. The main results are that within the axisymmetric ansatz the $o$ vortex becomes the most stable configuration in a strong magnetic field, that it hosts two spin-degenerate zero-energy Majorana bound states at $k_{3}=0$ protected by a $P_{3}$-symmetry winding number, and that the BdG magnetization profiles differ substantially from the GL-order-parameter-based estimates.","tokens_in":23924,"tokens_out":11224,"duration_ms":120267,"significance":"The result, if it holds, is a meaningful advance: it appears to be the first microscopic (quasiclassical plus BdG) description of vortex cores in $^{3}P_{2}$ superfluids, and it identifies a concrete strong-field regime in which topologically protected Majorana zero modes can exist in the cores of neutron-star vortices. The paper earns credit for doing genuine self-consistent Eilenberger calculations, for checking the zero modes both numerically and through the $P_{3}$ winding-number argument without parameter fitting, and for comparing BdG magnetizations against earlier GL estimates. The main caveat is the explicit restriction to axisymmetric configurations in Eq. (18); the paper itself lists nonaxisymmetric vortices as important future work. Therefore the stability and zero-mode statements are rigorously established only within that subspace, and the physical relevance to vortex matter in neutron stars remains conditional on the absence of a lower-energy nonaxisymmetric vortex such as the double-core vortex known in $^{3}$He-B.","major_comments":[],"minor_comments":[{"comment":"The BdG spin density in Eq. (39) is computed with an energy cutoff of $15T_{c}$, but no convergence check with respect to this cutoff is reported; please add a sentence confirming convergence or a short convergence test, since the continuum contribution is important for the total magnetization profile.","section":"Sec. III.A, Eq. (39)"},{"comment":"The citation number 82 in Sec. IV corresponds to an unnumbered footnote placed at the end of the reference list; this should be formatted as a proper numbered reference, or the in-text citation should be changed to the footnote marker, to avoid confusion.","section":"Sec. IV and reference list"},{"comment":"The word 'nonaxisymetric' in the first paragraph of Sec. IV should be corrected to 'nonaxisymmetric'.","section":"Sec. IV"},{"comment":"The abstract and conclusion are properly careful about the axisymmetric restriction, but the neutron-star framing in the introduction suggests a global stability claim; adding one sentence in the introduction stating explicitly that global stability against nonaxisymmetric perturbations is not addressed would make the scope fully transparent.","section":"Abstract and Sec. IV"}],"recommendation":"minor_revision","confidential_remarks":"This is a solid, carefully scoped paper. The nonaxisymmetric-vortex concern is real but is acknowledged by the authors and does not invalidate the axisymmetric theorem; I would not block publication on it. The main requested changes are cosmetic and one convergence statement. Please also check the formatting issue with reference 82."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Let me cut to it: this is the first microscopic Eilenberger-plus-BdG calculation of vortices in a 3P2 superfluid, and the central result is the o vortex at strong field carrying two zero-energy Majorana bound states. The Majorana claim is not numerology: the P3-symmetry winding number predicts two zero modes, and the BdG spectra show them at k3=0 with no fitted parameters. The self-consistent Eilenberger solutions also reveal a v-vortex analog that GL theory missed, and the free-energy comparison is a genuine variational one within the axisymmetric ansatz. The v vortex wins at zero field, while the o vortex becomes most stable as the field grows, with the D4-BN-o2 branch taking over near VZ ~ 0.9Tc. That is a clean, internally consistent story. The magnetization comparison between BdG and GL is also a useful warning that GL core profiles are unreliable.\n\nThe soft spots are both acknowledged in the paper, and neither kills the core calculation. First, everything lives in the axisymmetric sector of Eq. (18). The abstract and Sec. IV are careful to say the stability statement is for axisymmetric vortices. The stress-test worry is about the reach, not the internal logic: if a nonaxisymmetric double-core vortex (as in 3He-B) or a pair of half-quantized vortices wins energetically, the o vortex is no longer the physically realized state and the P3-protected zero modes lose their immediate relevance. The paper's argument that the symmetric traceless tensor may suppress the double-core deformation is qualitative; that is the main gap between the calculation and a neutron-star conclusion. Second, the core-magnetization profile is computed without feeding the magnetization back into the order parameter. The authors say exactly this after Eq. (44); it makes the BdG spin-density plots indicative rather than final, though the topological zero modes do not depend on that feedback.\n\nOn the literature: the paper leans on earlier work from the same group plus the 3He-B vortex classification; the citations are appropriate, not padding. The numerics use fixed parameters (T=0.4Tc, omega_c=10Tc, kF xi0=4, R0/xi0=80) and a BdG energy cutoff; that is normal for this kind of calculation, and the consistency check against quasiclassical LDOS in Appendix D is a nice touch.\n\nWho gets value: anyone working on 3P2 pairing, neutron-star vortex microphysics, or Majorana zero modes in multicomponent superfluids. I agree with the conditional verdict. Send it to a serious referee rather than desk-rejecting; the referee should ask for the nonaxisymmetric caveat to be stated more prominently and for the magnetization-feedback limitation to be kept in view.","headline":"First microscopic Eilenberger-plus-BdG calculation for 3P2 vortices; the axisymmetric o-vortex/Majorana result is solid, but the neutron-star stability conclusion still waits on nonaxisymmetric vortices.","tokens_in":24520,"tokens_out":3227,"would_cite":true,"duration_ms":33988,"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":"The paper claims that in a strong magnetic field the o vortex is the most stable axisymmetric vortex in a $^3P_2$ superfluid and that it binds two zero-energy Majorana fermions in its core.","keywords":["3P2 superfluid","neutron star","vortex","Majorana zero mode","Bogoliubov-de Gennes equation","Eilenberger equation","o vortex","spin-triplet pairing"],"falsifier":"Solve the order-parameter equations without imposing the axisymmetric ansatz of Eq. (18), allowing elliptic or double-core deformations, and compare free energies at fields $V_Z \\approx 0.9$–$1.5\\,T_c$; if any such nonaxisymmetric solution has lower free energy than the D4-BN-o2 vortex, the most-stable claim fails. A direct numerical check that the two $\\ell = 0$ modes at $k_3 = 0$ remain exactly degenerate under arbitrary small P3-breaking perturbations would also test the protection claimed for the zero modes.","tokens_in":23454,"feed_emoji":"🌀","tokens_out":5412,"duration_ms":52460,"temperature":0.7,"pith_summary":"This paper argues that quantized vortices in the spin-triplet $^3P_2$ superfluid thought to exist inside neutron-star cores should be described by a fully microscopic theory, not just by Ginzburg–Landau order-parameter profiles. Solving the Eilenberger equation self-consistently and then the Bogoliubov–de Gennes equation, the authors identify a family of axisymmetric vortex solutions whose stability depends on magnetic-field strength. Their central result is that in a strong magnetic field the so-called o vortex is the most stable axisymmetric configuration, and it binds two zero-energy Majorana fermions in its core, protected by a discrete magnetic rotation symmetry called P3. The paper also finds that the magnetization profile around the core computed microscopically differs sharply from earlier order-parameter-only estimates. If correct, this identifies the vortex microphysics that governs neutron-star cores in strong fields and gives a concrete astronomical setting where topologically protected Majorana zero modes could occur.","feed_headline":"Strong-field vortex binds two Majorana fermions","feed_subtitle":"In 3P2 superfluids the o vortex beats the v vortex and carries protected zero modes.","key_machinery":"The work combines the quasiclassical Eilenberger equation with the self-consistent gap equation to determine the axisymmetric order parameter $A(R) = \\sum_{M=-2}^{2} \\gamma_M(\\rho) e^{i(\\kappa-M)\\theta} \\Gamma_M$, and the Bogoliubov–de Gennes equation to obtain quasiparticle eigenenergies and core magnetization. The o vortex is defined as the axisymmetric configuration that preserves the P3 magnetic $\\pi$-rotation symmetry, meaning the components $\\gamma_0$ and $\\gamma_{\\pm 2}$ are real and $\\gamma_{\\pm 1} = 0$; this symmetry, combined with particle-hole conjugation, defines a chiral operator $\\check{\\Gamma}$ and a one-dimensional winding number $w_{\\rm 1d} = 2$, which guarantees two zero-energy Majorana modes at $k_3 = 0$. Stability is compared through the Luttinger–Ward free-energy functional, and the BdG spin density is contrasted with the Ginzburg–Landau order-parameter formula.","core_discovery":"Within axisymmetric vortex configurations in $^3P_2$ superfluids, the paper demonstrates that the o vortex—the configuration preserving all three discrete symmetries P1, P2, and P3—is energetically the most stable vortex in the presence of a strong magnetic field, specifically in the D4-BN phase. The o vortex contains two spin-degenerate zero-energy Majorana bound states at $k_3 = 0$, protected by the P3 symmetry, which yields a one-dimensional winding number $w_{\\rm 1d} = 2$. In contrast, the v vortex spontaneously breaks P3, so its would-be zero modes mix and split away from zero energy. The paper further shows that self-consistent microscopic order parameters include induced components with $M = \\pm 2$ that are absent from the earlier Ginzburg–Landau ansatz, and that the local spin density computed from Bogoliubov–de Gennes quasiparticles is finite at the o-vortex core where the order-parameter-only picture gives zero. These results constitute the first microscopic calculation of a single vortex in a multicomponent superfluid with a finite Zeeman field.","pith_inferences":["The paper leaves nonaxisymmetric vortices open; if the symmetric-traceless tensor structure of the $^3P_2$ order parameter indeed suppresses the double-core vortex that destabilizes the o vortex in $^3$He-B, then the o vortex with its two Majorana modes could be the global ground-state vortex in magnetar interiors, not merely the best axisymmetric one.","Two Majorana zero modes per vortex could support non-Abelian statistics beyond the single-mode case; the paper raises this question but does not answer it, so a concrete next step is to compute braiding properties of the o-vortex pair.","A testable extension is to feed the self-consistent vortex profiles and BdG spectra into a vortex-dynamics calculation and ask whether the spectral-flow force they produce changes predicted glitch relaxation timescales compared with Ginzburg–Landau-based estimates.","The same microscopic machinery could be applied to the cyclic and ferromagnetic $^3P_2$ phases, where 1/3-quantized non-Abelian vortices have been predicted, to check whether those cores also host zero-energy fermions."],"forward_implications":["In strong magnetic fields of the size relevant to magnetars, the stable vortex state is the o vortex, so each vortex line there should carry a topologically protected pair of Majorana zero modes.","The v vortex, which has no protected zero modes, is stable only at weaker fields, so the presence or absence of Majorana fermions depends on magnetic-field strength through the vortex configuration.","Core magnetization computed from fermionic quasiparticles differs drastically from Ginzburg–Landau estimates, so neutron-star cooling and vortex dynamics calculations that use order-parameter-only profiles miss a substantial contribution.","At sufficiently strong fields the $M = \\pm 1$ components of the v vortex vanish, the vortex continuously recovers the P3 symmetry, and the two zero modes reappear, making the strongly magnetized regime the natural place to look for these modes.","If these vortices exist in neutron-star cores, the fermion bound states contribute to the spectral-flow force and therefore to vortex unpinning and glitch dynamics, not only to the static structure."],"supporting_citations":[{"why":"Supplies the quasiclassical Eilenberger transport equation used to determine the order parameters self-consistently.","marker":"[71]"},{"why":"Provides the Luttinger–Ward free-energy functional used to compare vortex stability.","marker":"[73]"},{"why":"Gives the discrete-symmetry classification of o and v vortices in $^3$He-B that the paper adapts to $^3P_2$ superfluids.","marker":"[57]"},{"why":"Establishes that the P3 magnetic $\\pi$ rotation protects two spin-degenerate Majorana zero modes in the o-vortex core.","marker":"[65]"},{"why":"Reviews and generalizes the P3-protection argument for Majorana fermions in vortex cores.","marker":"[66]"},{"why":"Shows that the uniform $^3P_2$ superfluid is a topological class-DIII superfluid and supplies the microscopic starting point for the present calculations.","marker":"[52]"},{"why":"Provides the Ginzburg–Landau phase diagram and vortex boundary conditions that the microscopic solutions are compared against.","marker":"[30]"},{"why":"Identifies fermion bound states in vortex cores as relevant to neutron-star cooling and vortex dynamics.","marker":"[67]"}],"fun_headline_variants":["o Vortex in 3P2 Superfluid Binds Two Majoranas","Microscopic Theory Reveals Stable o Vortex with Majorana Modes","First Microscopic Vortex Study in 3P2 Superfluids","o Vortex Most Stable in 3P2 Superfluid, Hosts Majorana Pair","Two Majorana Fermions Found in 3P2 Superfluid Vortex"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculations only consider vortices that keep circular symmetry around the vortex line, with boundary conditions fixed by the uniform UN or D4-BN phases; if a nonaxisymmetric vortex has lower free energy, the o vortex would not be the global most-stable state and the P3 symmetry protecting the two Majorana zero modes could be broken.","fun_headline_variants_meta":{"raw":{"variants":["o Vortex in 3P2 Superfluid Binds Two Majoranas","Microscopic Theory Reveals Stable o Vortex with Majorana Modes","First Microscopic Vortex Study in 3P2 Superfluids","o Vortex Most Stable in 3P2 Superfluid, Hosts Majorana Pair","Two Majorana Fermions Found in 3P2 Superfluid Vortex"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000614,"raw_usage":{"total_tokens":2851,"prompt_tokens":940,"completion_tokens":1911,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":556,"completion_tokens_details":{"reasoning_tokens":1808}},"tokens_in":556,"tokens_out":1911,"duration_ms":11880,"temperature":1.0,"reasoning_tokens":1808,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:52:29.282601+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Solve the order-parameter equations without imposing the axisymmetric ansatz of Eq. (18), allowing elliptic or double-core deformations, and compare free energies at fields $V_Z \\approx 0.9$–$1.5\\,T_c$; if any such nonaxisymmetric solution has lower free energy than the D4-BN-o2 vortex, the most-stable claim fails. A direct numerical check that the two $\\ell = 0$ modes at $k_3 = 0$ remain exactly degenerate under arbitrary small P3-breaking perturbations would also test the protection claimed for the zero modes.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the Luttinger–Ward free-energy functional used to compare vortex stability."},{"cited_title":"Tsutsumi , author T","cited_arxiv_id":null,"evidence_quote":"Establishes that the P3 magnetic $\\pi$ rotation protects two spin-degenerate Majorana zero modes in the o-vortex core."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Identifies fermion bound states in vortex cores as relevant to neutron-star cooling and vortex dynamics."}],"review_version":1}