Flux-tube bouquets in the superfluid interiors of neutron stars ✛
Context
This work sits among studies of superfluid and superconducting phases in neutron star cores. Nearby works address vortices in neutron superfluidity, supersolidity at liquid interfaces, invisible quantum vortices and their collective behavior. The closest arXiv neighbors study magnetohydrodynamics of superfluid and superconducting cores, mixing of charged and neutral Bose condensates, superconducting phases in two-component models, and minimal models for pulsar glitches. Unlike most neighboring works, here for the first time both neutron vortices and proton flux tubes are present simultaneously in one computational cell, allowing extraction of pinning energies.
Methods
A two-dimensional Ginzburg–Landau model for two coupled condensates, neutron and proton, was used, without phase dragging. The free energy was minimized numerically in two ways: with quasiperiodic boundary conditions fixing the total numbers of vortices and flux tubes, and with fixed phases to compute interaction energy at fixed separation. Calculations were performed on grids of about 5 points per unit length using the Adam optimizer in PyTorch. Model parameters were fixed at reference values from Wood and Graber: κ=1.44, R=0.4, ε=0.1, while couplings α and g were varied within |α|≤0.1, |g|≤0.1.
Results
With attractive local couplings (α<0 or g<0), overlap of a vortex and a flux tube is energetically favorable: the pinning energy increases with distance, and the minimum occurs at d=0. The binding energy per unit length for typical parameters is about 0.13 MeV/fm. In configurations with one vortex and nine flux tubes, bouquets form: 6 to 9 flux tubes gather around the vortex, and in a dense bouquet (α=−0.1) individual flux tubes are indistinguishable in condensate amplitude or magnetic field, yet their nine phase singularities remain topologically independent. For flux tubes without a vortex, a type-1.5 regime is confirmed: the interaction force changes sign at a finite distance, and clusters self-assemble with a characteristic spacing of about 7 (≈220 fm), consistent with the prediction from a two-particle calculation.
Implications
This work complements neighboring studies by explicitly identifying a new type of bound state, the bouquet, and showing that the type-1.5 regime persists even without phase dragging. Unlike earlier works that considered vortices and flux tubes separately or with dragging, here pinning energies are obtained for the first time from the joint relaxation of both defects in a single cell. This provides a quantitative basis for models of dissipative interaction between the superfluid and charged components.
Future development
A natural extension is to include phase dragging, a three-dimensional geometry with curved and misaligned defects, and microscopic calibration of parameters α and g for a realistic equation of state. One could move from two-dimensional cross-sections to full three-dimensional simulations with quasiperiodic boundary conditions, as in the works of Wood et al. This would allow testing the stability of bouquets and clusters in volume and linking them to observable phenomena, such as glitches or thermal evolution.
Impact
The results may affect models of dissipative interaction in neutron star cores, particularly the calculation of mutual friction and the interpretation of roto-chemical heating. They also provide benchmarks for laboratory experiments with multicomponent superconductors and superfluid gases.
Next steps
The next step is to include phase dragging and perform three-dimensional calculations with realistic parameters. It is also necessary to compute effective transport coefficients (e.g., electron scattering off a bouquet) and check the impact on Urca reaction rates.
Key open problems
In neutron star astrophysics, the connection between superfluid vortices and the magnetic field, and the mechanism of glitches, remain open questions. This work directly addresses the problem of vortex pinning on flux tubes and its role in dissipation, as well as the question of the type of superconductivity in the core (type-I, type-II, or intermediate type-1.5).
The knowledge machine's view
This section is addressed to the author of the paper and their colleagues. Neighbouring works were found by meaning-based search across our archive, not by word overlap. We do not review: we show the field around the work.
What the work does
The authors built a two-dimensional Ginzburg–Landau model with local couplings for neutron and proton condensates in the superfluid core of a neutron star. By numerically minimizing the free energy, they studied the interaction between neutron vortices and proton flux tubes, discovering bound states called 'bouquets,' where a vortex holds several flux tubes. They also confirmed a type-1.5 regime for flux tubes without a vortex, where clusters self-assemble at a finite distance.
Where its strength is
The strength of the work lies in simultaneously considering both types of defects in a single computational cell and extracting quantitative pinning energies, which provides a foundation for models of dissipative coupling between the superfluid and charged components.
Where it could go next
- It might be interesting to see how bouquets of vortices and flux tubes affect large-scale magnetohydrodynamics of the core, for example, turbulent amplification of magnetic fields in young neutron stars.based on: Cosmic Kneading: Turbulent Amplification of Magnetic Fields in Newborn Neutron Stars arXiv:2607.21990v1
- Closely related is work on the role of mass, rather than rotation, in generating strong magnetic fields; perhaps one could investigate how bouquets alter magnetic field evolution depending on stellar mass.based on: Neutron Star's Sleeping Volcano: How Mass, Not Spin, Gives Birth to Magnetic Monsters arXiv:2607.17556v1
- Another direction is the influence of bouquets on pulsar glitches: nearby minimal models where vortex pinning plays a key role exist, and one could try to incorporate flux tube effects there.based on: Neutrinos: Born in Black Hole Coronas arXiv:2601.01533v2
Worth trying
- Test how pinning energies change when accounting for entrainment by adding the corresponding term to the free energy.
- Calculate the pinning force of a bouquet on a three-dimensional flux tube lattice to assess collective effects.
- Use the obtained binding energies to estimate the dissipation rate during vortex motion in a glitch model.
Why this may matter
Such structures can affect angular momentum transport and magnetic flux in neutron stars, linking microscopic superconductivity physics to observable phenomena such as glitches and magnetic field evolution.
Nearby works in our archive
- Cosmic Kneading: Turbulent Amplification of Magnetic Fields in Newborn Neutron Starssimilarity 0.54
- A Star with a Double Heartsimilarity 0.54
- Magnetar Traffic Controllers: How Spin Ignites Super-Bright Lightsimilarity 0.54
- Stellar Ashsimilarity 0.53
- A Neutron Star as a Musical Instrumentsimilarity 0.53
- Neutrinos: Born in Black Hole Coronassimilarity 0.52
- Neutron Star's Sleeping Volcano: How Mass, Not Spin, Gives Birth to Magnetic Monsterssimilarity 0.52
The archive is sparse around this work: the closest review is only 0.54, against a usual 0.64. It looks like few have drilled nearby.
A dense bouquet of nine flux tubes around a single vortex creates a magnetic flux several times larger than that of a single flux tube, although topologically it consists of nine separate flux quanta. The spacing between flux tubes in a cluster is about 220 fm, which is thousands of times smaller than the size of an atom.
Key numbers
- pinning energy per unit length: up to ~0.13 MeV/fm
- characteristic spacing between flux tubes in a cluster: ~220 fm
- number of flux tubes in a bouquet: 6–9
- computational cell size: 50×50 (in units of coherence length)