Stellar Remnants

Nuclear Pasta

A teaspoon of neutron-star matter weighs 6 billion tons — here's what holds it together

Nuclear pasta is the exotic state of matter ~1 km below a neutron star's surface, where nucleons rearrange into shapes named after Italian pasta — gnocchi (spheres), spaghetti (rods), lasagna (sheets), bucatini (tubes), and Swiss cheese (sphere holes) — at densities 1014–1015 g/cm3. Discovered theoretically in 1983 by Ravenhall, Pethick & Wilson; the strongest known material in the universe (~1010× stronger than steel) according to 2018 simulations by Caplan et al. A 5-mm cube weighs ~50 million tons. The crust transitions from a normal nuclear lattice (gnocchi) at 0.5 ρnuclear to uniform neutron fluid at ρnuclear ≈ 2.6 × 1014 g/cm3, with pasta phases between. Beneath, a superfluid core spins frictionlessly at up to 716 Hz (PSR J1748-2446ad).

  • Density1014–1015 g/cm3
  • Phasesgnocchi, spaghetti, lasagna, bucatini, Swiss cheese
  • Strength~1010× steel
  • Teaspoon mass~6 billion tons
  • Crust thickness~1 km
  • Superfluid spin record716 Hz (PSR J1748-2446ad)

Interactive visualization

Press play, or step through manually. The visualization is yours to drive — try it before reading on.

Open visualization fullscreen ↗

Watch the 60-second explainer

A condensed visual walkthrough — narrated, captioned, under a minute.

Why nuclear pasta matters

Nuclear pasta is not a curiosity — it sits at the center of three open problems in physics. Each one couples laboratory measurements to astrophysical observations through this strange phase of matter.

  • The nuclear equation of state. The relationship between pressure and density at supranuclear densities is unknown. We can squeeze atomic nuclei in colliders to ~3 ρ0 for 10-22 seconds; we cannot reach the 5–10 ρ0 of a neutron-star core. Pasta sits between the well-understood lattice regime and the uniform-fluid regime, and its phase boundaries depend sensitively on the symmetry energy slope L (somewhere between 40 and 80 MeV). Pin down where pasta begins and ends, and you constrain L; constrain L, and you predict neutron-star radii to within ~0.5 km.
  • Gravitational wave signatures. When two neutron stars merge — as observed in GW170817 on August 17 2017 — the inspiral waveform's late phase carries an imprint of tidal deformability Λ. Soft equations of state (smaller Λ) give more compact stars; stiff ones (larger Λ) give larger radii. LIGO-Virgo measured Λ1.4 < 800, ruling out the stiffest pre-merger models. Pasta-phase contributions to bulk and shear viscosity also damp post-merger oscillations at 2–4 kHz, a band Cosmic Explorer and Einstein Telescope will reach.
  • Pulsar glitches. Crustquakes alone cannot account for the spin-up signature of large glitches like Vela's. Two-component models invoke crustal superfluid that decouples from the lattice and stores angular momentum, then releases it when vortices unpin. The pasta layer expands the moment-of-inertia reservoir by ~50% in some equation-of-state families, reconciling the inferred superfluid fraction with theoretical predictions.
  • Cooling curves. Pasta phases conduct heat anisotropically. Lasagna sheets channel phonons within the plane and block them across; spaghetti rods do the opposite. X-ray observations of cooling young neutron stars (Cassiopeia A, ~340 years old) constrain thermal conductivity in this layer to roughly 1022–1023 erg/(cm s K), informing transport-coefficient calculations.
  • Magnetic field decay. Hall drift and Ohmic dissipation in the inner crust set magnetar lifetimes. Pasta's elongated structures bend resistive timescales by factors of 10 to 100, explaining why magnetars retain ~1014–1015 G fields for 104–105 years before becoming radio-quiet.
  • r-process nucleosynthesis. Crustal matter ejected during binary mergers — a few percent of a solar mass at velocities ~0.2c — synthesizes half the periodic table beyond iron. The mass distribution and lanthanide fraction depend on whether the pre-ejecta sat in pasta or uniform fluid; kilonova spectra (AT2017gfo) confirmed strontium and constrained the opacity, anchoring r-process models.

The five pasta phases in order of depth

Walking down a 1-km column from outer crust to neutron-fluid interface, you cross five distinct geometries. Each forms because, at densities where Coulomb repulsion between protons no longer dominates, the system minimizes surface energy by changing topology.

  • Gnocchi (~0.5 ρ0, depth ~700 m). Spherical nuclei in a body-centered cubic lattice, immersed in a sea of dripped neutrons. Lattice spacing ~30 fm. Each "gnocchi" contains roughly 200–1000 nucleons.
  • Spaghetti (~0.6 ρ0, depth ~800 m). Cylindrical rods of nuclear matter arranged on a hexagonal lattice. The rods can be tens of fm long. Coulomb repulsion is reduced at the cost of higher surface area than spheres of equal volume — net energy savings come from packing.
  • Lasagna (~0.7 ρ0, depth ~900 m). Parallel slabs of nuclear matter separated by neutron fluid. Slab thickness ~5 fm, separation ~10 fm. Caplan et al. (2018) showed that breaking these sheets requires shear stresses of ~1030 erg/cm3 — the source of the "strongest material" claim.
  • Bucatini / anti-spaghetti (~0.8 ρ0, depth ~1000 m). Cylindrical tubes of neutron fluid running through a continuous nuclear-matter background. Topologically inverse to spaghetti.
  • Swiss cheese / anti-gnocchi (~0.9 ρ0, depth ~1050 m). Spherical bubbles of low-density neutron fluid embedded in nuclear matter. The last identifiable structure before the smooth uniform phase.

The transition density, 0.5 to 1.0 ρ0, brackets a thin region but corresponds to a finite physical column because pressure rises slowly with depth in this regime. Total pasta column thickness is several hundred meters in most equation-of-state models.

Why pasta is the strongest material in the universe

In 2018 Caplan, Schneider and Horowitz ran molecular-dynamics simulations on Blue Waters at the National Center for Supercomputing Applications, evolving 5 million nucleons through shear deformation. They measured the breaking strain εbreak at which the lasagna and spaghetti structures rupture, finding values up to 0.3 — meaning pasta deforms by 30% before fracture. Combined with its enormous shear modulus (~1030 erg/cm3, ten orders larger than steel's ~1011 erg/cm3), this gives a yield strength of roughly 1029 erg/cm3 ≈ 1010× steel.

The mechanism is geometric. Lasagna sheets must reconnect across a stiff matrix to accommodate strain; the topology blocks slip planes that ordinary metals exploit. Defects propagate through narrow channels and pin under their own self-energy. Comparable structures in liquid crystals and block copolymers show the same enhancement — pasta is just a relativistically dense version.

That strength has observational consequences: a maximum mountain on a neutron-star surface is bounded by crustal yield. Pasta lets the crust support deformations of order 10-6 in radius before cracking — corresponding to a continuous-wave gravitational signal LIGO has searched for in known pulsars without detection. The non-detections give upper limits on ellipticity below 10-7 for some sources, consistent with — but not yet probing — pasta-strength predictions.

The superfluid core beneath

Below the pasta layer, density exceeds ρ0 and matter becomes a uniform fluid of overwhelmingly neutrons (~95%) with a small admixture of protons (~5%) and electrons. Within minutes of the supernova, neutrinos cool the star below ~109 K, and pairing transitions occur:

  • Neutron 1S0 superfluid in the inner crust. Critical temperature Tc ≈ 1010 K at densities below ~0.7 ρ0. The pairing gap peaks at ~3 MeV.
  • Neutron 3P2 superfluid in the outer core. Tc ≈ 109 K at densities of 1–3 ρ0. The triplet pairing is more uncertain — gaps span 0.05 to 1 MeV in different many-body calculations.
  • Proton 1S0 superconductor. Type-II at most densities; magnetic field threads the core in flux tubes of ~1015 Mx each.

The superfluid carries angular momentum in quantized vortices, each holding a circulation quantum of h/(2mn) ≈ 2 × 10-3 cm2/s. A neutron star spinning at 100 Hz has ~1017 vortices per cm2, on average a few microns apart. As the pulsar spins down through magnetic torque, vortices must migrate outward — but they pin to lattice nuclei or to flux tubes. Pinning and unpinning drive the glitches we observe.

Record-holding pulsars

  • PSR J1748-2446ad. Discovered 2005 in Terzan 5 by Hessels et al. Spin frequency 716 Hz (period 1.396 ms). Equatorial speed ~24% c. Recycled by accretion in a binary; current companion is a 0.14 solar-mass star.
  • PSR J0740+6620. Mass 2.08 +/- 0.07 solar masses (Cromartie et al. 2019, refined by Fonseca et al. 2021). Radius 12.39 (+1.30/-0.98) km via NICER X-ray timing. Currently the most stringent upper bound: the equation of state must support at least 2.08 solar masses without collapse.
  • SGR 1806-20. Magnetar with surface field ~2 × 1015 G. Released 2 × 1046 erg in 0.2 seconds on December 27 2004 — the brightest extrasolar gamma-ray event ever recorded, detected through a partial reflection off Earth's ionosphere.
  • RX J1856.5-3754. Closest known isolated neutron star at 123 pc. Surface temperature ~700,000 K, no detected pulsations, used as a thermal-emission benchmark.
  • The Crab pulsar (PSR B0531+21). Born in 1054 CE. Spin frequency 30.225 Hz dropping by 3.86 × 10-10 Hz/s — the energy lost equals 5 × 1038 erg/s and powers the entire Crab Nebula.
  • The Vela pulsar (PSR B0833-45). Spin frequency 11.19 Hz. Glitches every ~3 years with relative spin-up ΔΩ/Ω ~ 10-6 — the canonical "large glitch" used to test crustal superfluid models.

What might lie deeper

The innermost core of a 2-solar-mass neutron star may not be ordinary nucleons at all. Above ~5 ρ0, several exotic phases compete:

  • Hyperons. Strange-quark-bearing baryons (Λ, Σ, Ξ) become energetically favorable. The "hyperon puzzle": their appearance softens the equation of state and tends to violate the 2-solar-mass observed minimum. Resolution requires repulsive hyperon-nucleon interactions still being mapped at JLab and J-PARC.
  • Quark matter. If chiral symmetry restores, deconfined u-d-s quark matter may form in the deepest core — possibly as a color superconductor in the color-flavor-locked (CFL) phase. Hybrid stars with quark cores and nucleon shells are observationally indistinguishable from pure-nucleon stars in current data.
  • Pion or kaon condensates. Bose-condensed mesons could appear at intermediate densities, lowering pressure and making compact stars more compact still.

None of these is confirmed. Pasta, by contrast, is well-established in simulation — its existence depends on physics (Coulomb plus surface tension) we already understand. The deeper exotic phases may or may not exist; the pasta layer almost certainly does.

Common misconceptions

  • "Neutron stars are pure neutrons." The crust contains protons, electrons, and neutron-rich nuclei (zirconium-122, iron-78, krypton-118 are predicted at successive depths). The core is ~95% neutrons but retains ~5% protons and electrons in beta equilibrium. The deepest interior may even contain quarks or hyperons.
  • "They're solid." Only the outer layer (km-thick crust including pasta) is a rigid solid. Below that is superfluid, which flows without viscosity, threaded by magnetic flux tubes from a type-II proton superconductor. The "surface" you'd land on is a kilometer-thin shell over a sea.
  • "All the same density." Density varies by six orders of magnitude from atmosphere (~106 g/cm3) to inner core (~4 × 1015 g/cm3). The label "1014 g/cm3" applies only to the saturation point at the crust-core boundary.
  • "Smaller = less massive." Mass and radius are inverted in neutron stars: stiffer equations of state give larger but less compact stars; softer ones give smaller and denser ones. All known neutron stars cluster between 1.2 and 2.1 solar masses but radii span 10–14 km — most of the variation is structure, not mass.
  • "Pasta is metaphorical." The shapes are literal geometric configurations confirmed in molecular dynamics. The names are descriptive, not analogical.
  • "You'd be crushed instantly on the surface." Surface gravity is ~2 × 1014 cm/s2, ~2 × 1011 g. You'd be flattened to atomic scale in microseconds — but not by the density, by the gravitational field. The surface itself is a 1-cm atmosphere of carbon or hydrogen at ordinary atomic density.
  • "Nuclear pasta and spaghettification are the same." Different phenomena. Spaghettification is the tidal stretching of objects falling into black holes. Nuclear pasta is a static phase of matter inside neutron stars.
  • "Neutron stars don't have magnetic fields." They have the strongest magnetic fields known. Ordinary radio pulsars: 1011–1013 G. Magnetars: 1014–1015 G. For comparison, Earth's field is ~0.5 G, a refrigerator magnet ~100 G, the strongest sustained lab field ~45 tesla = 4.5 × 105 G.
  • "They emit light because they're hot." Young neutron stars do glow thermally (~106 K surface), but most observed pulsars shine via beamed magnetospheric emission powered by rotational kinetic energy, not heat. Magnetars instead tap their magnetic field's reservoir.

Frequently asked questions

What is nuclear pasta and why is it called that?

Nuclear pasta is a phase of matter found about 1 km below a neutron star's surface where nucleons (protons and neutrons) rearrange into non-spherical shapes due to competition between nuclear attraction and Coulomb repulsion. As density rises from 1013 to 1014 g/cm3, individual spherical nuclei (gnocchi) merge into rods (spaghetti), then sheets (lasagna), then inverse tubes (bucatini), then bubble-riddled fluid (Swiss cheese) before becoming uniform nuclear matter. Ravenhall, Pethick and Wilson named the phases after Italian pasta in their 1983 paper because the geometric shapes match — the analogy stuck because it captures the morphology precisely.

How dense is a neutron star really?

Densities span six orders of magnitude. The outer atmosphere is ~106 g/cm3, similar to white dwarf material. The outer crust (a few hundred meters thick) sits at 109 to 1011 g/cm3, lattices of neutron-rich iron and zirconium nuclei. The inner crust including pasta runs 1011 to 1014 g/cm3. At the crust-core boundary, density reaches saturation density ρ0 = 2.6 × 1014 g/cm3 (also written as 0.16 nucleons per fm3). The core averages 5 to 10 ρ0; the deepest center may reach 10 to 15 ρ0 ~ 4 × 1015 g/cm3. A teaspoon (5 mL) of core matter weighs ~6 billion tons — comparable to the mass of Mount Everest compressed into a sugar cube.

Why doesn't it collapse into a black hole?

Two pressures hold gravity at bay. First, neutron degeneracy pressure: the Pauli exclusion principle forbids two neutrons from occupying the same quantum state, producing a quantum-mechanical pressure independent of temperature. Second, the strong nuclear repulsion at short range: nucleons resist compression below ~0.5 fm separation. Above the Tolman-Oppenheimer-Volkoff (TOV) limit — roughly 2.2 to 2.5 solar masses depending on the equation of state — even these pressures fail and the star collapses to a black hole. The current observed maximum is PSR J0740+6620 at 2.08 +/- 0.07 solar masses, with radius 12.39 km measured by NICER in 2021.

What is a superfluid core?

Below ~109 K — quickly reached after birth — neutrons in the core pair up via the strong attractive force at low momentum (the same mechanism behind Cooper pairing in superconductors). Paired neutrons form a Bose condensate that flows without viscosity. Crust neutrons enter a superfluid 1S0 state; core neutrons form 3P2 pairs at higher density. Protons in the core form a type-II superconductor, threading the magnetic field into quantized flux tubes. Vortices in the neutron superfluid carry the angular momentum: each vortex pins to a flux tube or crustal lattice site, releasing it suddenly during glitches.

How fast can a neutron star spin?

The fastest known is PSR J1748-2446ad in the globular cluster Terzan 5, rotating at 716 Hz (716 times per second, period 1.396 ms). Its equator moves at ~24% the speed of light. The theoretical maximum (mass shedding limit) is ~1500 Hz for a typical 1.4 solar mass neutron star with a stiff equation of state — at higher rates centrifugal force exceeds gravity at the equator and material flies off. Many fast pulsars are 'recycled' millisecond pulsars: spun up by accreting matter from a binary companion, transferring orbital angular momentum to the star over ~108 years. Birth spin rates for isolated pulsars are typically 10 to 100 Hz.

What happens during a starquake?

As a young pulsar spins down, the equatorial bulge shrinks. The rigid solid crust resists deformation until stress exceeds its yield strength — then it cracks, releasing energy and angular momentum suddenly. Observers see a 'glitch': the rotation rate jumps upward by 10-9 to 10-6, then relaxes over days to months as crustal superfluid couples to the new spin. The Vela pulsar glitches every ~3 years with Δω/ω ~ 10-6. Magnetar starquakes are far more violent: SGR 1806-20 released ~1046 erg in 0.2 seconds on December 27 2004 — equal to the Sun's total output over 250,000 years — by cracking its crust under magnetic stress of ~1015 G.