Why Hubble and Webb Just Found 27 Ancient Frozen Worlds Beyond Neptune

Two linked studies published in The Astronomical Journal have established the deepest direct census of the outer solar system ever conducted, using coordinated observations from the James Webb Space Telescope and the Hubble Space Telescope to identify 27 previously unknown trans-Neptunian objects (TNOs). Orbiting in the freezing perimeter beyond Neptune, these newly cataloged bodies represent primordial planetesimals—the foundational blocks of planet building that stalled mid-assembly roughly 4.6 billion years ago. The faintest of these detections reaches a near-infrared magnitude of 29.3, a threshold comparable to resolving the light of a tiny cluster of fireflies on the surface of the Moon from an observation post on Earth.

The survey pushed direct optical and infrared detection down to diameters of five kilometers, penetrating a size regime previously accessible only through statistical inference, crater-counting models on visited outer moons, and stellar occultation surveys. The empirical return from this multi-observatory campaign, led by doctoral researchers Marielle Eduardo at the University of Victoria and Anastasia Morgan at Northern Arizona University, delivers two findings that contradict long-standing assumptions in planetary astrophysics.

First, the survey uncovered an unexpected deficit of sub-ten-kilometer objects. Rather than revealing a dense population of collisional fragments generated by billions of years of mutual impacts, the size distribution follows a remarkably flat power-law slope of 0.29. Second, spectrophotometric analysis across visible and infrared wavelengths shows that these tiny bodies retain the distinct surface colors and chemical signatures of their multi-hundred-kilometer counterparts, defying expectations that impact excavation would have stripped away their weathered outer mantles.

Taken together, these 27 frozen worlds beyond Neptune serve as a rigorous test case for the physical mechanisms governing the early protoplanetary disk. The observations challenge classical bottom-up hierarchical accretion, constrain the nature of high-velocity collisional grinding in the scattered disk, and offer direct physical evidence for rapid, gravity-driven planetesimal formation via aerodynamic pebble concentration.

       SUMMARY OF THE DUAL-OBSERVATORY DISCOVERY PARAMETERS
┌───────────────────────────┬─────────────────────────────────────────────────┐
│ Primary Observatories     │ James Webb Space Telescope (NIRCam)             │
│                           │ Hubble Space Telescope (ACS / WFC3)             │
├───────────────────────────┼─────────────────────────────────────────────────┤
│ Cataloged Objects         │ 27 new Trans-Neptunian Objects (TNOs)           │
│ Hubble Optical Counterparts│ 13 objects with full 0.35–3.2 μm coverage       │
├───────────────────────────┼─────────────────────────────────────────────────┤
│ Survey Footprint          │ 0.05 square degrees (deep pencil-beam profile)  │
├───────────────────────────┼─────────────────────────────────────────────────┤
│ Limiting Sensitivity      │ F150W2 magnitude 28.8 (40% recovery limit)      │
│                           │ Optical r-band equivalent ~ 29.8                │
│                           │ Faintest confirmed detection: mag 29.3          │
├───────────────────────────┼─────────────────────────────────────────────────┤
│ Measured Physical Diameters│ ~5 km to 40 km (assuming ~0.15 geometric albedo)│
├───────────────────────────┼─────────────────────────────────────────────────┤
│ Size-Frequency Slope      │ Power-law cumulative index: 0.29 (+0.08 / -0.07)│
├───────────────────────────┼─────────────────────────────────────────────────┤
│ Core Empirical Findings   │ 1. Steep deficit of sub-10-km bodies            │
│                           │ 2. Color bifurcation preserved at tiny sizes    │
│                           │ 3. Identical size slopes across Hot/Cold belts  │
└───────────────────────────┴─────────────────────────────────────────────────┘

The Deepest Pencil Beam: 27 Planetesimals at the Edge of Detection

Capturing sunlight reflected by a five-kilometer object located beyond thirty astronomical units (AU) presents an extreme photon-budget challenge. Sunlight dims inversely with the square of distance ($1/r^2$) as it travels outward toward the Kuiper Belt, and the minute fraction reflected from an icy surface decays by another factor of $1/r^2$ before reaching detectors near Earth, resulting in an overall flux attenuation scaling with the fourth power of heliocentric distance ($1/r^4$). For objects with diameters ($D$) below ten kilometers and geometric albedos ($p_V$) typically hovering between 0.04 and 0.15, the incoming photon stream is overwhelmed by the zodiacal light background and sensor read noise.

Ground-based optical facilities, including the 8.2-meter Subaru Telescope equipped with Hyper Suprime-Cam, encounter a hard sensitivity limit for moving outer solar system objects at approximately apparent magnitude $r \sim 25.5$ to $26.0$. At that threshold, direct detection is restricted to bodies roughly twenty to thirty kilometers in diameter or larger. To break through this observational wall, the research campaign designed a coordinated pencil-beam survey combining Webb’s 25-square-meter primary mirror and infrared sensitivity with Hubble’s optical resolution.

The survey targeted a concentrated sky footprint of roughly 0.05 square degrees. Webb’s Near-Infrared Camera (NIRCam) gathered deep integrations using wide-band filters, primarily F150W2, which offers exceptional throughput between 1.0 and 2.0 micrometers while minimizing thermal background noise. This enabled the pipeline to achieve a 40 percent detection completeness limit at magnitude 28.8, corresponding to an optical r-band brightness of approximately 29.8.

Simultaneously, the Hubble Space Telescope directed its Advanced Camera for Surveys (ACS) and Wide Field Camera 3 (WFC3) at the identical celestial coordinates. Hubble’s shorter optical bandpasses (spanning 0.35 to 0.8 micrometers) captured the reflected visible sunlight necessary to calculate broad-spectrum color indices. Out of the 27 moving targets isolated by Webb, 13 were bright enough in visible wavelengths to be securely identified in Hubble’s exposures.

APPARENT FLUX ATTENUATION OF DISTANT PLANETESIMALS
Solar Irradiation Flux:     F_inc  ∝ 1 / r^2
Reflected Earthward Flux:   F_obs  ∝ (p_v * D^2) / (r^2 * Δ^2)
When r ≈ Δ (Outer Solar System):  F_obs  ∝ (p_v * D^2) / r^4

Resulting observational penalty: Moving an icy body from 10 AU to 40 AU 
decreases its detected flux by a factor of 256, necessitating extreme integration 
depths and multi-frame stacking to extract sub-10-km planetesimals.

The physical scale of these objects underscores the sensitivity of the dataset. Assuming a median geometric albedo of 15 percent, the faintest targets measure approximately five to ten kilometers across. These bodies represent the smallest planetesimals ever detected through direct reflected sunlight beyond Neptune. They inhabit an elusive middle ground between the massive classical dwarf planets, such as Pluto, Makemake, and Haumea, and the tiny sub-kilometer cometary debris fields whose existence has been recorded primarily through momentary occultations of background stars.

The Shift-and-Stack Architecture: Reconstructing Motion from Noise

The detection of these objects did not occur through traditional single-frame imaging. At magnitude 29, individual exposures captured by Webb or Hubble contain virtually zero statistically significant signal from an object moving across the frame; the individual photons are scattered among background detector read noise and diffuse astrophysical glare. To extract the signal, the research consortium deployed advanced shift-and-stack processing techniques.

                 PRINCIPLE OF SHIFT-AND-STACK PROCESSING
      Frame t_0                  Frame t_1                  Frame t_2
┌──────────────────┐       ┌──────────────────┐       ┌──────────────────┐
│  * (Star)        │       │  *               │       │  *               │
│                  │       │                  │       │                  │
│       · (TNO)    │       │         · (TNO)  │       │               ·  │
│  *               │       │  *               │       │  *               │
└──────────────────┘       └──────────────────┘       └──────────────────┘
                            │                          │
                            ▼ Vector Shift (-v_x, -v_y)▼
                   ┌──────────────────────────────────┐
                   │  *   *   * (Streaked Stars)      │
                   │                                  │
                   │         ● (Coadded TNO Signal)   │
                   │  *   *   *                       │
                   └──────────────────────────────────┘
                   Coadded frame boosts TNO S/N by √N

Shift-and-stack data reduction operates by acquiring a long sequence of consecutive exposures of the same patch of sky. Stationary background sources, including distant galaxies and galactic stars, remain fixed in celestial coordinates, while trans-Neptunian objects drift at rates dictated by their orbital velocity and Earth's (or the space observatories') reflex parallactic motion. For objects located between 30 and 50 AU, this apparent angular velocity typically ranges between 1 and 4 arcseconds per hour along the ecliptic plane.

The computational pipeline executes a systematic search over a multidimensional parameter space:

Velocity Vector Testing: The processing pipeline shifts individual calibrated frames across thousands of discrete velocity vectors $(\dot{\alpha}, \dot{\delta})$ corresponding to kinematically plausible bound solar orbits.
Frame Co-addition: For each trial vector, the pixel arrays are shifted and co-added. When the trial velocity matches the true sky-projected motion of a trans-Neptunian object, the faint photons collected across every frame fall precisely onto the same synthesized pixel coordinates.
Signal-to-Noise Boosting: The signal-to-noise ratio ($S/N$) of the moving target accumulates roughly with the square root of the number of combined exposures ($\sqrt{N}$), elevating a sub-threshold detection into a statistically verifiable peak. Concurrently, stationary background stars and galaxies are smeared into linear streaks, dispersing their localized flux and allowing clean digital subtraction.
Machine Learning Filtering: Searching millions of synthetic orbital trajectories across noisy background arrays generates thousands of spurious candidate detections. To resolve this false-alarm issue, the team integrated a convolutional neural network trained on synthetic point-spread functions (PSFs) embedded into realistic background noise. The classifier evaluated candidate detections based on their spatial profile, flux stability across sub-stacks, and temporal consistency, rejecting false artifacts while preserving genuine moving bodies.

Without this digital synthetic tracking infrastructure, identifying 27 faint moving targets would have required a single aperture far exceeding the physical dimensions of Webb. The success of this algorithmic integration demonstrates that computational synthesis is now just as critical to pushing the observational edge as physical mirror surface area.

The Missing Fragment Paradox: Challenging the Collisional Cascade

The size distribution of small planetary bodies functions as a physical ledger of their dynamical history. When a population of planetesimals is left undisturbed, its size-frequency distribution (SFD) mirrors the primordial outcome of the accretion process. Conversely, if that population undergoes extensive collisional evolution, the exchange of kinetic energy during mutual hypervelocity impacts systematically fractures larger parents into smaller progeny, establishing an equilibrium size distribution known as a collisional cascade.

                     SIZE DISTRIBUTION POWER-LAW SLOPES
Log Cumulative Number [N(>D)]
  │                                
  │       \   Classical Dohnanyi Cascade Slope (q ≈ 3.5, α ≈ 0.50)
  │        \  [Predicted: Dense sea of tiny shattered fragments]
  │         \
  │          \                  
  │           \           
  │            \          Observed Webb-Hubble Survey Slope (α = 0.29)
  │             \───────- [Observed: Pronounced deficit of small bodies]
  │                      \
  │                       \
  │                        \  Primordial Break Scale (~50-100 km)
  │                         \
  └──────────────────────────┴──────────────────────────── Log Diameter (D)
   1 km                     10 km                       100 km

In 1969, astrophysicist John S. Dohnanyi formulated the foundational mathematical framework for collisional cascades, demonstrating that under idealized conditions—where material strength remains scale-invariant and catastrophic fragmentation governs impacts—the differential size distribution $dN(D) \propto D^{-q} dD$ inevitably relaxes to a power-law index of $q \approx 3.5$. Expressed in terms of cumulative logarithmic sky-plane magnitude distributions:

$$\Sigma(m) \propto 10^{\alpha m}$$

where the magnitude slope relates directly to the diameter index via:

$$\alpha = \frac{q - 1}{5}$$

For a standard Dohnanyi fragmentation cascade where $q = 3.5$, the corresponding magnitude slope is $\alpha = 0.50$. Extended numerical simulations of the outer solar system, incorporating velocity dispersions generated during giant planet migrations, consistently suggested that bodies below twenty kilometers should exhibit steep slopes ranging from $\alpha \sim 0.40$ to $\alpha \sim 0.60$. Models predicted that a deep pencil-beam survey piercing the 28th-magnitude floor should encounter a dense population of shattered debris.

The empirical data from Marielle Eduardo and colleagues delivered the opposite result. The measured power-law slope for the newly detected sample settled at:

$$\alpha = 0.29^{+0.08}_{-0.07}$$

This slope translates to a differential size index of roughly $q \approx 2.45$, far below the baseline threshold for an active collisional cascade. The survey revealed a pronounced absence of sub-ten-kilometer planetesimals.

┌─────────────────────────────────┬───────────────────┬───────────────────────┐
│ Dynamic Reservoir               │ Characteristic SFD│ Physical Mechanism    │
│                                 │ Slope Index (α)   │ Driving Distribution  │
├─────────────────────────────────┼───────────────────┼───────────────────────┤
│ Dohnanyi Equilibrium Model      │ α ≈ 0.50 (q=3.50) │ Catastrophic fragment-│
│ (Theoretical Baseline)          │                   │ ation equilibrium     │
├─────────────────────────────────┼───────────────────┼───────────────────────┤
│ Large TNOs (D > 100 km)         │ α ≈ 0.70 to 1.10  │ Primordial run-away   │
│ (Observed via ground surveys)   │ (q=4.50 to 6.50)  │ gravitational capture │
├─────────────────────────────────┼───────────────────┼───────────────────────┤
│ Sub-10-km Deep Webb/Hubble Survey│ α = 0.29 (+0.08/-0.07)│ Streaming instability │
│ (Eduardo et al. 2026 Sample)    │ (q ≈ 2.45)        │ with low collision rate│
└─────────────────────────────────┴───────────────────┴───────────────────────┘

This deficit yields vital dynamical insights into the outer edges of the solar system. For the size distribution to remain this shallow, the volume of space occupied by these objects must have possessed a substantially lower spatial mass density during its developmental epochs than classical planet-formation models assumed. Alternatively, the relative impact velocities between these bodies must have been insufficient to drive full catastrophic disruption.

If relative velocities remain low—as is characteristic of dynamically undisturbed orbits—impacts frequently result in bouncing, mild localized cratering, or gentle sticking, rather than fragmentation. The 27 objects identified by the space observatories demonstrate that the outer fringe of the planetary system avoided extreme collisional grinding, preserving an early architecture of accretion that escaped the violent restructuring seen across the asteroid belt.

The Retention of Primordial Skin: Spectrophotometric Uniformity Across Scale

While the size-frequency distribution measures the physical dimensions of a population, color spectrophotometry reveals its surface composition and historical exposure to space weathering. Trans-Neptunian objects generally fall into distinct dynamical classes, which trace different evolutionary paths:

Dynamically "Cold" Classical TNOs: These bodies inhabit orbits with low inclinations ($i < 4^\circ$) and low eccentricities ($e < 0.1$), remaining close to the original, near-circular plane of the protoplanetary disk. They are considered pristine relics that formed in situ beyond 40 AU, largely protected from major orbital disruptions. Spectroscopically, cold classicals are characterized by distinct red to ultra-red surface colors, caused by the prolonged irradiation of volatile hydrocarbons.
Dynamically "Hot" TNOs: These objects occupy orbits with elevated inclinations ($i > 15^\circ$) and high eccentricities, often crossing the orbital path of Neptune. They originated closer to the Sun—within the primordial giant-planet region between 15 and 30 AU—before being scattered outward during the rapid orbital divergence and migration of Jupiter, Saturn, Uranus, and Neptune. Their surfaces display more neutral, grayish-blue to moderately reddish colors, reflecting distinct formation conditions and different volatile ice concentrations.

                 SPECTRAL REFLECTANCE DIVERGENCE (0.35 - 3.2 μm)
Relative Reflectance (Normalized at 1.0 μm)
  │                                    
3 │                                 / Cold Classical Relics (Ultra-Red Tholins)
  │                                /  [Preserved identical across 10 km to 500 km]
2 │                        . - - -'   
  │                   . - '           
1 │  - - - - - - - - '─────────────── Hot Scattered Population (Neutral / Grey)
  │  [Hubble Optical] │ [Webb Infrared]
0 └───────────────────┴──────────────┴──────────────────────── Wavelength (μm)
 0.35                0.8            2.0                      3.2

Before this multi-observatory survey, planetary scientists widely hypothesized that size would fundamentally alter these spectral relationships. Because small bodies have tiny gravitational potentials, any impact event should easily breach their outer, radiation-weathered layers, casting fragments into space and uncovering unweathered interior ice.

Over billions of years of exposure to solar wind ions, solar ultraviolet radiation, and galactic cosmic rays (GCRs), surface methane ($CH_4$), methanol ($CH_3OH$), and nitrogen ($N_2$) undergo radiolytic decomposition, building complex, carbon-rich macromolecular polymers known as tholins. These organic crusts give unperturbed surfaces their distinct red coloration.

Impacts were expected to shatter this thin skin. A five-kilometer body, according to the standard model, should be an excavated chunk of pure subsurface water ice ($H_2O$) mixed with basic volatiles, giving it a neutral or blue-sloped reflectance spectrum that washed out the classic "hot" versus "cold" color distinctions.

The spectrophotometric analysis led by Anastasia Morgan overturned this assumption. By matching Hubble’s optical photometry with Webb’s infrared filter responses, the research team assembled spectral energy distributions across the 0.35 to 3.2 micrometer range for 13 of the detected bodies.

EXPECTED VS. OBSERVED SPECTRAL BEHAVIOR FOR SUB-10-KM TNOS

EXPECTED (Impact Stripping Hypothesis):
  Small Body -> Severe Collisional Churn -> Stripped Radiation Mantle 
  -> Exposed Pristine Subsurface Ice -> Bluer, Homogenized Spectra Across All Classes

OBSERVED (Morgan et al. 2026 Empirical Data):
  Cold Small Bodies -> Retain Ultra-Red Surface Slopes (Identical to Large Siblings)
  Hot Small Bodies  -> Retain Neutral/Moderately Red Slopes (Class Integrity Intact)
  -> Planetesimals Remain Primordial Down to Single-Digit Kilometer Scales

The small dynamically cold objects occupied the exact same narrow, ultra-red reflectance sequence that characterizes large cold classicals. Meanwhile, the dynamically hot small bodies displayed the wider, more neutral color distribution typical of large scattered disk objects.

This preservation of surface characteristics down to single-digit-kilometer scales reveals that these frozen worlds beyond Neptune are not random chunks chipped off larger parents during subsequent epochs. Instead, they appear to be original planetesimals that have preserved their primordial surfaces largely intact since the birth of the planetary system.

The Streaming Instability Triumphs: Planetesimal Formation Across Variable Disk Regimes

The empirical evidence—a flat size distribution combined with preserved surface chemistry—provides fresh support for modern aerodynamic accretion models over traditional bottom-up assembly.

For decades, the standard explanation for planet formation relied on hierarchical accretion, a theoretical framework where microscopic dust grains slowly collide and adhere through van der Waals forces, gradually growing into pebbles, boulders, kilometer-scale planetesimals, and eventually planetary embryos. This bottom-up pathway, however, faces significant physical barriers:

The Radial Drift Barrier: As dust grains grow into centimeter- to meter-sized boulders, they encounter head-on aerodynamic drag from the surrounding gaseous disk. Because the gas feels an outward radial pressure gradient, it orbits at a slightly sub-Keplerian velocity. Solid boulders, moving at unattenuated Keplerian speeds, lose orbital angular momentum to the slower gas, causing them to spiral into the central star in just centuries—far faster than pairwise collisions can build them into stable, larger bodies.
The Fragmentation and Bouncing Barriers: Laboratory experiments show that when silicates and icy grains reach millimeter- to centimeter-scale dimensions, impact velocities between them often exceed their cohesive strength. Rather than sticking together, colliding particles shatter or bounce off one another, bringing growth to an abrupt halt.

                  STREAMING INSTABILITY ACCRETION PATHWAY
1. Dust Settles to Midplane      2. Aerodynamic Clumping        3. Gravitational Collapse
┌─────────────────────────┐    ┌─────────────────────────┐    ┌─────────────────────────┐
│ Gas Disk: Sub-Keplerian │    │ Gas Drag Feedback Loops │    │ Local Pebble Cloud      │
│ ═══════════════════════ │    │ · Pebbles decelerate gas│    │ Density Exceeds Roche   │
│ · · · · · · · · · · · · │ -> │ · Pressure traps form   │ -> │ Limit: Sudden Collapse  │
│ Pebble Layer Settles    │    │ · Massive overdensities │    │ Into 10-100 km Bodies   │
│ ═══════════════════════ │    │   assemble rapidly      │    │ [No meter-scale drift]  │
└─────────────────────────┘    └─────────────────────────┘    └─────────────────────────┘

The streaming instability model, first formulated by Andrew Youdin and Jeremy Goodman in 2005 and refined through magnetohydrodynamic simulations by Anders Johansen and collaborators, bypasses these growth barriers entirely.

In this framework, particles do not slowly grow one step at a time through pairwise collisions. Instead, when aerodynamic pebbles settle into the midplane of a protoplanetary disk, they begin to push back on the surrounding gas. This collective feedback triggers an aerodynamic instability that concentrates the pebbles into dense, filamentary clouds.

When the local density of these pebble clouds surpasses the Roche density threshold, the cloud’s self-gravity overwhelms the local tidal shear of the central star. The entire filament then collapses directly under its own weight, transforming a loose cloud of pebbles into a gravitationally bound planetesimal in a single step.

The size-distribution data gathered by Eduardo and colleagues aligns closely with the predictions of the streaming instability. High-resolution numerical simulations of gravitational collapse driven by this instability consistently produce a distinctive size distribution:

At larger diameters (typically $D \gtrsim 50$ to $100$ kilometers), the mass distribution follows a steep power-law slope reflecting the broad structure of the collapsing clouds.
Below this transition scale, the size distribution flattens dramatically, producing a shallow power-law slope ($q \lesssim 2.5$) that yields a noticeable deficit of smaller bodies.

THE DISK-INVARIANCE PRINCIPLE IN PLANETESIMAL ACCRETION

Dynamically Cold Classicals          Dynamically Hot Objects
(Formed at 40 - 50 AU)               (Formed at 15 - 30 AU)
  · Low gas density                    · High gas density
  · Extended orbital periods           · Faster dynamical timescales
  · Low ambient temperature (~30 K)    · Higher ambient temperature (~50 K)
                   \                  /
                    ▼                ▼
          Observed Size-Distribution Power-Law Slope:
                       α ≈ 0.29 (Both Populations)
                                  │
                                  ▼
      Conclusion: Planetesimal formation via aerodynamic cloud 
      collapse produces an invariant initial mass function that 
      is largely independent of local protoplanetary disk conditions.

The survey uncovered an even broader insight: both the dynamically cold classical objects and the dynamically hot scattered bodies share the same shallow size distribution slope of 0.29.

This shared slope provides compelling evidence for the universal behavior of planetesimal accretion. The cold classicals assembled in situ in a cold, diffuse, outer environment beyond 40 AU, whereas the hot scattered population grew within the dense, warmer, and more dynamically turbulent neighborhood of the giant planets (15 to 30 AU).

Despite these stark differences in local temperature, gas surface density, and orbital velocity, both regions produced an identical size distribution. This shared trait indicates that streaming instability collapse operates via self-regulating mechanics, producing a characteristic size distribution that is largely independent of local disk conditions.

Dual-Observatory Synthesis: The Optical-Infrared Diagnostic Framework

The characterization of these small planetesimals illustrates a wider principle: high-precision astrophysics at the outer limits of detection increasingly relies on the synchronized pairing of complementary space-based facilities. Neither observatory could have unraveled this mystery on its own.

                     DUAL-OBSERVATORY SYNERGY MATRIX
┌────────────────────────────────┬───────────────────────────────┬──────────────────────────────┐
│ Operational Parameter          │ James Webb Space Telescope    │ Hubble Space Telescope       │
├────────────────────────────────┼───────────────────────────────┼──────────────────────────────┤
│ Primary Optical Regime         │ Near-to-Mid Infrared          │ Ultraviolet to Visible       │
│                                │ (0.6 μm to 28.0 μm)           │ (0.2 μm to 1.7 μm)           │
├────────────────────────────────┼───────────────────────────────┼──────────────────────────────┤
│ Collecting Area                │ 25.4 m² (Segmented Beryllium) │ 4.5 m² (Monolithic Glass)    │
├────────────────────────────────┼───────────────────────────────┼──────────────────────────────┤
│ Discovery Role                 │ Primary Detection Engine:     │ Diagnostic Characterization: │
│                                │ High infrared sensitivity     │ Visible optical bandpasses   │
│                                │ pushed limiting detection     │ measured volatile-dependent  │
│                                │ down to magnitude 29.3        │ spectral slopes              │
├────────────────────────────────┼───────────────────────────────┼──────────────────────────────┤
│ Primary Diagnostic Yield       │ 1. Orbital velocity vectors   │ 1. Tholin absorption metrics │
│                                │ 2. Limiting detection flux    │ 2. Color-class verification  │
│                                │ 3. Deep size-frequency index  │ 3. Optical-to-infrared SED   │
└────────────────────────────────┴───────────────────────────────┴──────────────────────────────┘

Operated alone, the James Webb Space Telescope functions as an extraordinary detection platform. Its large, gold-coated beryllium primary mirror and cryogenically cooled NIRCam instrument easily deliver the sensitivity needed to capture the faint infrared reflections from tiny, distant objects. Webb alone, however, has an observational blind spot: its optical throughput drops sharply at wavelengths below 0.6 micrometers.

Because the key spectral markers of solar system space weathering—specifically the steep absorption slopes produced by complex tholins—are concentrated across ultraviolet and blue-visible wavelengths (0.35 to 0.60 micrometers), relying solely on Webb's infrared data makes it difficult to separate an ancient, weathered, hydrocarbon-rich surface from a younger, impact-excavated, iron-bearing mineral mixture.

Hubble provides the missing piece of the puzzle. While its 2.4-meter primary mirror lacks the collecting area to discover the faintest targets on its own, its ACS and WFC3 instruments offer unmatched optical resolving power. Hubble's blue- and visible-light measurements (F435W, F606W, and F814W) set the baseline for an object’s overall spectral slope.

When researchers anchored Webb’s infrared photometry against Hubble’s visible-light measurements, they created continuous Spectral Energy Distributions (SEDs) covering 0.35 to 3.2 micrometers. This unified optical-to-infrared window enabled the team to clearly verify that these small bodies maintain the chemical signatures of their larger family members, a diagnostic confirmation that neither facility could have achieved on its own.

                 UNIFIED SPECTRAL COVERAGE (0.35 - 3.2 μm)
0.35 μm         0.6 μm           1.0 μm          2.0 μm           3.2 μm
  ├─── Hubble ACS / WFC3 ──────────┤
  │    (Optical Surface Slopes)    │
  │    [Tholin Absorption Edge]    │
  │                                ├───────── Webb NIRCam ──────────┤
  │                                │         (Infrared Thermal /    │
  │                                │          Volatile Ice Bands)   │
  └────────────────────────────────┴────────────────────────────────┘

This coordinated methodology provides an operational blueprint for future outer-system surveys. As astronomical investigations push toward ever-fainter targets, physical mirror size must be paired with wide, cross-wavelength coverage to turn bare positional tracks into rich chemical profiles.

Survey Strategy Realignment: From Rubin’s Wide Net to Roman’s Depth

The discovery of these 27 distant, icy bodies provides crucial context for upcoming wide-field astronomical facilities, primarily the Vera C. Rubin Observatory in Chile and the Nancy Grace Roman Space Telescope.

       NEXT-GENERATION SOLAR SYSTEM SURVEY ARCHITECTURE
┌───────────────────────────┬───────────────────────────┬───────────────────────────┐
│ Metric / Facility         │ Vera C. Rubin (LSST)      │ Nancy Grace Roman Space   │
│                           │                           │ Telescope                 │
├───────────────────────────┼───────────────────────────┼───────────────────────────┤
│ Observation Strategy      │ Ground-based wide synoptic│ Space-based wide-field    │
│                           │ (Southern sky coverage)   │ infrared surveys          │
├───────────────────────────┼───────────────────────────┼───────────────────────────┤
│ Limiting Magnitude        │ r ~ 24.5 (single visit)   │ F184 ~ 26.7 (wide tiers)  │
│                           │ r ~ 27.5 (10-yr co-add)   │ ~ 29.5 (targeted deeps)   │
├───────────────────────────┼───────────────────────────┼───────────────────────────┤
│ Expected TNO Yield        │ ~ 10,000 to 40,000 bodies │ Thousands of faint TNOs;  │
│                           │ (primarily D > 30-50 km)  │ bridges wide/deep realms  │
├───────────────────────────┼───────────────────────────┼───────────────────────────┤
│ Primary Role in           │ Establishes the regional  │ Probes the small-body     │
│ Small-Body Science        │ population architecture   │ size distribution down to │
│                           │ across the sky            │ single-digit kilometers   │
└───────────────────────────┴───────────────────────────┴───────────────────────────┘

The Vera C. Rubin Observatory’s Legacy Survey of Space and Time (LSST) is set to reshape the study of the outer solar system by scanning the southern sky every few nights with an 8.4-meter telescope and a 3.2-gigapixel camera. Over its planned ten-year run, LSST is projected to catalog between 10,000 and 40,000 new trans-Neptunian objects, expanding the known architectural framework of the Kuiper Belt.

Yet the Hubble-Webb survey highlights an important physical limit for ground-based observatories. Even under ideal atmospheric conditions, LSST’s single-exposure limit of $r \sim 24.5$ restricts its real-time moving object pipeline to bodies larger than roughly thirty to fifty kilometers across.

While co-adding years of imaging will push its detection limit toward $r \sim 27.5$, ground-based data cannot reliably pierce the sub-ten-kilometer domain where the shallow size distribution was just identified. LSST will excel at cataloging the macroscopic components of the outer solar system, but the delicate, low-mass planetesimals will remain hidden from ground-based view.

                 COMPLEMENTARY SURVEY REGIMES
Spatial Sky Coverage
  │  ┌────────────────────────────────────────────────────────┐
  │  │ Rubin Observatory (LSST)                               │
  │  │ Broad Synoptic Net (Tens of thousands of degrees)      │
  │  │ Focus: Macro-architecture (D > 30-50 km)               │
  │  └────────────────────────────────────────────────────────┘
  │                 ┌─────────────────────────────────────────┐
  │                 │ Nancy Grace Roman Space Telescope       │
  │                 │ Wide-Field Infrared Imaging             │
  │                 │ Bridges macro surveys and deep pencils  │
  │                 └─────────────────────────────────────────┘
  │                                     ┌─────────────────────┐
  │                                     │ Webb / Hubble Pencils│
  │                                     │ Micro-scale deep-cut│
  │                                     │ (D down to ~5 km)   │
  │                                     └─────────────────────┘
  └───────────────────────────────────────────────────────────── Limiting Depth
  Mag 24                                                      Mag 30

That crucial physical regime will be probed by NASA’s Nancy Grace Roman Space Telescope. Equipped with a primary mirror matching Hubble’s 2.4-meter diameter and a Wide Field Instrument that delivers a field of view one hundred times larger, Roman merges wide-angle scanning with space-based infrared sensitivity.

Roman’s high-latitude surveys and dedicated microlensing fields near the galactic bulge will image through deep space filters down to near-infrared magnitudes of 27 to 29. The telescope will discover thousands of faint trans-Neptunian bodies, bridging the gap between Rubin’s broad, sky-wide mapping and Webb’s ultra-deep, pencil-beam observations.

Data from Roman will test whether the flat power-law slope of 0.29 holds true across the entire sweep of the Kuiper Belt, or if the initial 27 bodies detected by Webb and Hubble reflect localized structures in the outer solar system.

These astronomical findings also provide vital context for direct spacecraft exploration. In January 2019, NASA’s New Horizons probe conducted an in-situ flyby of 486958 Arrokoth, a 36-kilometer contact binary preserved within the cold classical Kuiper Belt. Arrokoth surprised planetary scientists with its pristine, uncratered lobes and its uniform, ultra-red surface tholin blanket.

At the time, researchers debated whether Arrokoth represented a unique, gentle merger or a typical relic of the outer solar system. The discovery that tiny frozen worlds beyond Neptune systematically share these identical surface colors and show an absence of impact-shattered fragments confirms that Arrokoth’s smooth appearance reflects the prevailing conditions of early planetesimal accretion, establishing its history as standard across the pristine outer belt.

Principles of Early Accretion and Outer Solar System Architecture

Using this dual-observatory discovery as an analytical lens reveals three core physical principles that help refine our understanding of how planetary systems form:

┌─────────────────────────────────┬────────────────────────────────────────────────────┐
│ Principle                       │ Core Physical Meaning                              │
├─────────────────────────────────┼────────────────────────────────────────────────────┤
│ 1. Disk Invariance Principle    │ Aerodynamic pebble concentration builds planetes-  │
│                                 │ imals with a characteristic mass function that is  │
│                                 │ largely independent of local disk temperature, gas │
│                                 │ density, or radial distance from the central star. │
├─────────────────────────────────┼────────────────────────────────────────────────────┤
│ 2. Surface Memory Principle     │ In low-density outer disks, small icy bodies avoid │
│                                 │ disruptive collisions, preserving their radiation- │
│                                 │ altered surfaces over billions of years.           │
├─────────────────────────────────┼────────────────────────────────────────────────────┤
│ 3. Synthetic Depth Principle    │ Pushing solar system observational limits relies on│
│                                 │ multi-frame motion reconstruction and cross-       │
│                                 │ wavelength analysis rather than aperture alone.    │
└─────────────────────────────────┴────────────────────────────────────────────────────┘

1. The Disk Invariance Principle

Planetesimal formation via aerodynamic pebble concentration operates as a self-regulating physical process that is remarkably robust against changing environmental conditions.

Whether collapsing in the warm, dense inner protoplanetary disk where the dynamically hot population originated, or in the cold, thin, distant regions where the classical belt formed, the gravitational collapse of pebble clouds produces an identical initial size-distribution slope ($q \approx 2.45$, corresponding to $\alpha = 0.29$). The basic physics of streaming instability collapse depends primarily on the aerodynamic coupling between solids and gas—defined by the dimensionless Stokes number—rather than on the total ambient gas density of the disk.

2. The Surface Memory Principle

In dynamically undisturbed regions of the outer solar system, small icy bodies avoid severe, surface-stripping impact cascades.

Because the initial assembly of the outer disk left behind a low spatial mass density of small fragments, impact rates across the trans-Neptunian belt remained too low to erode the radiation-processed surfaces of small planetesimals. As a result, bodies measuring just five to ten kilometers across preserve their primordial, radiation-altered crusts, maintaining their original chemical signatures over billions of years without being reduced to pulverized, impact-excavated ice.

3. The Synthetic Depth Principle

Pushing the boundary of solar system exploration requires combining deep multi-frame orbital reconstruction with broad, cross-wavelength spectroscopy.

Aperture size alone cannot overcome the steep $1/r^4$ flux loss of distant planetesimals. Instead, modern astronomical surveys must coordinate optical and infrared instruments to extract faint moving targets from noise, turning isolated, moving points of light into detailed records of our solar system's earliest days.

The Unresolved Frontiers of Outer Solar System Research

The characterization of these 27 frozen worlds beyond Neptune establishes a new empirical foundation, but it also raises new questions that will guide future planetary science:

The True Scale of the Collisional Break: If the power-law slope remains flat at $\alpha = 0.29$ down to diameters of five kilometers, at what threshold does catastrophic disruption finally become evident? Future deep integrations using the Habitable Worlds Observatory or extremely sensitive stellar occultation arrays will be needed to determine whether the size distribution plunges toward zero below one kilometer, or if a secondary cascade of shattered debris appears at sub-kilometer scales.
The Full Volatile Inventory: While broadband optical and infrared photometry confirms that these small bodies retain distinct surface colors, existing signal-to-noise ratios cannot resolve individual spectroscopic absorption bands. Measuring precise surface volatile ratios—such as the balance of methane, water ice, carbon dioxide, and complex tholins—will require ultra-deep spectroscopy from Webb's NIRSpec instrument or future thirty-meter-class ground telescopes equipped with advanced adaptive optics.
The Far Boundary of the Planetesimal Disk: The 27 detected objects reside primarily in the core Kuiper Belt and inner scattered disk, between 30 and 50 AU. Yet ongoing trajectory and telemetry analysis from the New Horizons spacecraft, now cruising past 66 AU, indicates persistent impacts from fine interplanetary dust particles. This raises the intriguing possibility that the planetesimal belt extends much farther into deep space than currently mapped. Determining whether these distant dust grains are released by undiscovered populations of five-kilometer icy bodies remains a vital question for outer solar system astronomy.

The discovery of these 27 distant worlds confirms that the outer solar system retains a clear, uncorrupted memory of its primordial state. By looking past Neptune with our most sensitive space observatories, planetary scientists are no longer relying on theoretical models to guess how planetesimals first gathered in the early dark; they are reading the surviving record left behind in the frozen architecture of the solar system itself.