The core challenge with component-level optimization is computational cost. When you model every single part of a satellite in high fidelity, the number of calculations required to explore the design space becomes astronomical. Traditional genetic programming methods, while powerful for exploring complex solution spaces, struggle to keep up when faced with such granularity. To address this bottleneck, Li and Dong proposed a strategy they call Multi-granularity Genetic Programming and Batch Sequential Design, or MGGPBSD. The name might sound like a mouthful, but the underlying concept is elegant: instead of evaluating every potential design configuration at the highest possible level of detail, the algorithm strategically switches between different levels of modeling fidelity depending on where it is in the search process.

At its heart, MGGPBSD is a memetic algorithm, meaning it combines global search strategies with local refinement techniques. The global search is handled by genetic programming, a technique inspired by biological evolution where candidate solutions are represented as tree structures and iteratively refined through operations like crossover and mutation. In this case, the tree structures represent the hierarchical organization of satellite components. The algorithm explores this component tree to identify promising architectural configurations—decisions about which types of components to use and how they should be interconnected.

But genetic programming alone isn’t enough. Once the algorithm has identified a promising tree structure in a given generation, it needs to optimize the continuous design parameters associated with that configuration—things like power ratings, dimensions, and operating characteristics. This is where the batch sequential design method comes into play. Rather than evaluating parameters one at a time, the batch sequential approach groups evaluations into batches and iteratively refines the parameter estimates based on the results from each batch. This allows for rapid convergence to optimal parameter values within each architectural candidate, dramatically accelerating the inner loop of the optimization.

The researchers also developed a multi-granularity simulation model for satellite components. This means that different components can be modeled at different levels of fidelity—some with simplified analytical representations and others with more detailed numerical simulations. An on-demand switching criterion determines when to use the higher-fidelity model versus the coarser one. During the early stages of optimization, when the algorithm is still exploring broad architectural possibilities, the coarser models are sufficient. As the search converges toward specific candidate solutions, the more detailed models are invoked to ensure accuracy. This adaptive approach significantly reduces the total computational burden without sacrificing the quality of the final design.

To demonstrate the effectiveness of their method, the researchers applied MGGPBSD to a ground observation satellite as a case study. The results were striking: the proposed method achieved a 53 percent reduction in total computational cost compared to a baseline genetic programming approach, all while maintaining optimization accuracy. For satellite engineering teams working under tight schedules and limited computational budgets, a reduction of this magnitude is not merely an incremental improvement—it fundamentally changes what kinds of design problems become tractable.

The implications extend beyond just speed. Component-level optimization produces results that directly guide practical component selection in real engineering applications. Rather than abstract parameter values that must be translated into hardware choices, the algorithm outputs specific component configurations that engineers can implement directly. This bridging of the gap between optimization theory and engineering practice could accelerate satellite development cycles and enable more ambitious mission designs.

The research builds on a growing body of work in multidisciplinary design optimization for spacecraft. Previous efforts have explored collaborative optimization frameworks, metamodel-assisted techniques, and various hybrid approaches that combine different algorithmic strategies. What distinguishes MGGPBSD is its systematic integration of tree-structured genetic programming with batch sequential parameter optimization and multi-granularity modeling—all unified under the memetic algorithm paradigm. This combination allows the method to handle both the discrete architectural decisions and the continuous parameter refinements that characterize real satellite design problems.

Looking forward, the authors suggest that their approach could be extended to even more complex satellite systems and multi-satellite constellations. The key laboratory that supported this research—the Key Laboratory of Spacecraft Design Optimization and Dynamic Simulation Technologies, funded by China’s Ministry of Education—continues to push the boundaries of what computational optimization can achieve in aerospace engineering. As satellites become more sophisticated and mission requirements more demanding, tools like MGGPBSD will likely become indispensable for the next generation of space system designers.

The study, received in January 2026 and published in September 2026, represents a meaningful step toward democratizing access to advanced satellite design optimization. By reducing the computational barriers that have historically limited component-level optimization to well-resourced organizations, methods like MGGPBSD could empower a broader community of engineers and researchers to design the satellites of tomorrow.

The intellectual roots of the new framework reach into several distinct research traditions that have matured over the past two decades. Genetic programming, the evolutionary technique at the core of the method, has been applied to satellite problems before, including work on satellite system topology and parameter optimization and on multi-satellite cooperative task allocation. What those earlier applications demonstrated is that tree-structured representations are naturally suited to design problems where the arrangement of parts matters as much as the numerical values assigned to them. By evolving entire component architectures rather than fixed parameter vectors, the search can discover configurations that a conventional optimizer, locked into a single predefined structure, would never encounter.

The batch sequential design component draws on a separate line of statistical research. Sequential design of experiments has long been used to decide where to sample an expensive simulation next, and batch variants extend this idea by evaluating groups of points simultaneously, which is particularly valuable when simulations can run in parallel or when each evaluation carries fixed overhead. Prior applications of batch sequential designs range from accelerated life testing of polymer composites to minimum-energy design problems in process engineering, and sliced Latin hypercube constructions have provided a mathematical foundation for building such batches coherently. Importing this machinery into the inner loop of an evolutionary search is what allows the new method to refine continuous parameters quickly once an architecture has been selected.

The multi-granularity modeling idea also has a conceptual pedigree. Granular computing, a framework from approximate reasoning research, treats knowledge at multiple levels of detail and selects the granularity appropriate to the decision at hand. In the satellite context, this translates into a modeling hierarchy in which the same physical component can be represented by a fast surrogate when coarse guidance suffices and by a high-fidelity simulation when precision is required. Earlier work on multi-granularity modeling for remote sensing satellite effectiveness evaluation, and on uncertainty quantification for multi-granularity models in satellite optimization, helped establish that such switching can be done without corrupting the optimization results, providing a foundation the present study builds upon.

The memetic framing deserves attention as well. Memetic algorithms pair population-based global search with local improvement procedures, and they have proven effective in domains ranging from large-scale cooperative coevolution to multi-agent dispatching and feature selection. The essential insight is that neither global exploration nor local refinement alone is efficient for problems with mixed discrete and continuous structure. Global search alone wastes evaluations on poorly tuned architectures, while local refinement alone cannot escape the neighborhood of an inferior configuration. The hybrid strategy evaluated here applies that lesson to satellite engineering, using evolutionary operators to navigate the combinatorial space of component trees while delegating parameter tuning to the sequential design procedure within each generation.

The case study choice of a ground observation satellite is also meaningful. Earth observation missions involve tightly coupled subsystems—an imaging payload, its supporting optics, the attitude control system that stabilizes pointing, and the power and thermal systems that sustain operation—so component-level decisions propagate strongly across disciplinary boundaries. This makes such missions a demanding testbed for any optimization framework that claims to handle architectural and parametric decisions simultaneously. The reported 53 percent reduction in computational cost relative to a baseline genetic programming method, achieved without loss of accuracy, suggests that the granularity switching and batch refinement are doing genuine work rather than merely shifting expense between stages.

For practitioners, the practical significance lies in what becomes feasible at given computational budgets. Component-level optimization has historically been reserved for late-stage design or well-funded programs because exhaustive high-fidelity evaluation of every candidate configuration was prohibitively expensive. Methods that adaptively ration fidelity change that calculus, allowing design teams to explore component choices earlier in a program, when the cost of revisiting decisions is lowest. The authors note that data supporting the findings are available from the corresponding author upon reasonable request, and the work was partially supported by the Key Laboratory of Spacecraft Design Optimization and Dynamic Simulation Technologies under China’s Ministry of Education, reflecting the institutional infrastructure behind sustained aerospace optimization research at Beihang University.

Subject of Research: A multi-granularity genetic programming and batch sequential design method for reducing computational cost in component-level satellite optimization

Article Title: Multi-granularity Genetic Programming and Batch Sequential Design Method for Component-Level Satellite Optimization Design

Article References: Li, Z., & Dong, Y. (2026). Multi-granularity Genetic Programming and Batch Sequential Design Method for Component-Level Satellite Optimization Design. International Journal of Aeronautical and Space Sciences. https://doi.org/10.1007/s42405-026-01291-8

Image Credits: AI Generated

DOI: 10.1007/s42405-026-01291-8

Keywords: satellite design, genetic programming, memetic algorithm, batch sequential design, multi-granularity modeling, component-level optimization, aerospace engineering, space systems, optimization, computation cost reduction, Beihang University, multidisciplinary design optimization

The technology at the heart of the study is the Vertical Solar Array Technology, or VSAT, a concept NASA has been advancing for lunar polar missions. Unlike conventional flat panels, a VSAT features an elevated, vertically oriented photovoltaic panel that rotates to always face the Sun. This design is ideally suited to the lunar poles, where the Sun never rises high into the sky. Because the Sun skims just a few degrees above the horizon, a vertical rotating panel can receive sunlight at nearly normal incidence, maximizing the solar energy it captures. A single large VSAT measuring roughly fifty square meters, with a solar conversion efficiency of thirty percent, can produce about fifteen kilowatts of power in full sunlight. That is impressive for one unit, but early lunar activities will quickly outgrow such bespoke power sources, so the real challenge is combining hundreds of panels into a farm capable of generating power in the megawatt range.

That is where the geometry turns hostile. The same low Sun angle that makes vertical panels so effective at the poles also creates a serious problem: shadows. Because the Sun sits near the horizon, every VSAT casts an enormous shadow directly behind it, and a dense field of panels can shade one another severely, slashing the total energy output. On a hillside where the terrain itself can also block the Sun, self-shading between panels becomes a critical design constraint. Previous studies, including work by Ross and colleagues, had shown that lunar polar hills offer near-constant solar visibility but suffer abrupt shading from nearby structures, and that a megawatt farm would require hundreds of panels whose interactions had to be modeled carefully.

Earlier illumination studies of the lunar poles, such as those by Noda, Speyerer, Gläser, Mazarico and their colleagues, modeled the local horizon as a slowly changing elevation function around a site. These approaches revealed that a ridgeline on the Connecting Ridge beside Shackleton Crater could receive over eighty percent illumination across an eighteen-year precessional cycle. But they were not designed to handle the abrupt, discrete shadows cast by neighboring solar panels, and they were computationally far too expensive to evaluate the thousands of possible panel arrangements needed to optimize a farm of hundreds of units.

The new study introduces two novel ideas to break that computational bottleneck. The first is a binary sky bitmask technique, borrowed from terrestrial solar energy estimation. On Earth, sky imagers use bitmasks to capture clouds and scattered light, which vary continuously in brightness. The airless Moon is far simpler: sunlight is either blocked or not. The researchers built their terrain model from a cascade of digital elevation models, including a high-resolution five-meter DEM for the Connecting Ridge and coarser regional models from the Lunar Orbiter Laser Altimeter. For each candidate panel location, the narrow corridor of sky within a few degrees of the horizon where the Sun can appear is mapped onto a rectangular bitmask with an angular grid resolution of 0.05 degrees, fine enough that the solar disk spans ten pixels and eighty active bits. Computing the insolation at any moment then reduces to a fast bitwise AND operation between the sun bitmask and the sky bitmask, with bits packed into sixteen-bit integers for speed. The team validated the method by comparing it against a traditional ray-tracing calculation, finding the two power curves for a hundred-panel farm almost perfectly overlapped, and by checking their terrain registration against actual photographs from the Lunar Reconnaissance Orbiter.

With the fast computation in hand, the researchers turned to the design problem. They considered a grid of 6,166 candidate VSAT locations spaced twenty meters apart above the 1,900-meter elevation contour on the Connecting Ridge. Each panel was modeled as five meters wide and ten meters tall, with its bottom edge five meters above the surface, delivering up to fifteen kilowatts in full sun. A greedy optimization algorithm then built the farm panel by panel: at each step, the candidate location that added the most total energy to the farm, after accounting for the shading it cast on all previously placed panels and received from them, was permanently selected. Starting from the single best-lit spot on the ridge, the algorithm generated an ordered sequence of panel positions.

The first hundred panels in the energy-maximizing configuration produced a total of 10,695 megawatt-hours during the year 2030, an average of 1.22 megawatts, comfortably above NASA’s one-megawatt requirement. But the design had a hidden flaw. Twice per lunar day, the Sun’s direction aligned with the long diagonal of the panel layout, and most panels shaded one another simultaneously, causing dramatic dips in output known as brownouts. The energy-optimized arrangement had naturally exploited alignments that concentrated shading at times when other panels were already well lit, but these transient power troughs posed a real problem for a grid that must supply life-support systems without interruption.

The second novel idea addressed the brownouts directly. Instead of maximizing total integrated energy, the researchers optimized a logarithmic value function that rewards marginal power most when the farm is generating little of it. In effect, an extra kilowatt is worth far more during a power trough than during peak production, so the algorithm steers panel placement toward times and places that fill in the gaps. Empirically, the result proved insensitive to the scaling constant chosen for the value function. The reoptimized hundred-panel farm lost only two percent of its total annual output, delivering 10,522 megawatt-hours or an average of roughly 1.2 megawatts, yet the brownout troughs vanished entirely, leaving only the unavoidable blackouts that occur when the entire hill is plunged into darkness by the shadow of Mons Mouton during part of the year.

Those terrain-driven blackouts are the one problem no arrangement of solar panels can solve, and the study is explicit about the implication: survival power during dark periods must come from stored energy or a nuclear source. The researchers also quantified the diminishing returns of farm expansion. While the first panel contributes an average of 13.7 kilowatts, the hundredth adds only 10.8 kilowatts, a 21 percent reduction, and the five-hundredth adds just 7.7 kilowatts, down 44 percent, as self-shading and less favorable sites take their toll. Even so, the Connecting Ridge could sustain an average of 4.8 megawatts with five hundred panels, and the analysis indicates the chosen boundary comfortably accommodated the first five hundred placements, though larger farms would need more space. The economics of panel deployment versus marginal benefit will ultimately determine the optimal size of such an installation.

The study opens several avenues for future work. Real deployments will face slope constraints and other siting restrictions, which can be incorporated by simply removing candidate locations from the grid. More intriguingly, the authors suggest that if individual VSAT panels could move by meters over hours, coordinated motion could further suppress self-shading, a far harder optimization problem that might one day be tamed by deep learning or autonomous distributed computation, or even solved by a simple emergent behavior at each station. For now, the work provides a practical blueprint: by combining a fast bitmask illumination model with value-based optimization, lunar planners can design solar farms that grow gracefully from a handful of panels into a megawatt-class power grid, bringing permanent human habitation at the Moon’s south pole a significant step closer to reality.

The bitmask approach also carries practical implications for how lunar power infrastructure might be planned and expanded over time. Because pre-existing installations can simply be added to the initial sky bitmask at every candidate location, the method naturally accommodates incremental growth: a farm that begins with a handful of bespoke panels can be extended later without redesigning the whole site. This aligns with the staged reality of lunar exploration, in which early missions arrive with modest power needs and only later demand megawatt-scale generation as surface activity intensifies.

The study site itself was chosen for good reason. The Connecting Ridge, on the rim of Shackleton Crater near the south pole, is considered a strong candidate for high-activity landing zones, partly because the hypothesized presence of water ice and other frozen volatiles in permanently shadowed regions nearby makes the area attractive for resource utilization. The Sun’s path there is confined to a narrow swath only a few degrees above the horizon, shifting between solstice extremes and bounded by a margin of about two degrees to ensure no part of the solar disk escapes the modeled corridor from any vantage point in the study area.

There is also a strategic dimension to the work. By estimating how much power a single, well-chosen lunar hill can deliver before self-shading becomes excessive, the analysis offers long-term guidance on the eventual mix of solar and nuclear power that sustained lunar presence will require. If the best hilltops can supply several megawatts but no more, planners can determine at what point stored energy or fission systems must take over, informing architecture decisions years before hardware is deployed.

Subject of Research: Design optimization of a megawatt-scale vertical solar array farm for power generation at the lunar south pole

Article Title: Design optimization of a megawatt scale lunar VSAT solar array farm

Article References: Butler, W. H., & Freudenreich, Y. P. (2026). Design optimization of a megawatt scale lunar VSAT solar array farm. Space and Planetary Resources, 2(1), Article 4. https://doi.org/10.1007/s44461-026-00009-8

Image Credits: AI Generated

DOI: 10.1007/s44461-026-00009-8

Keywords: lunar south pole, VSAT, solar array farm, self-shading, sky bitmask, illumination modeling, Shackleton Crater, Connecting Ridge, lunar power systems, photovoltaics, space exploration, megawatt power

The central challenge the authors tackle is deceptively simple to state: neutron stars are born hot and magnetized, and every observable we measure — their X-ray glow, their spin periods, their sudden outbursts — is shaped by how heat and magnetic field decay and interact over time. Solving this problem requires numerically integrating two coupled families of equations: the heat transfer equation governing the cooling of the stellar interior, and the induction equation describing the evolution of the magnetic field. Both depend sensitively on microphysics such as thermal and electrical conductivities, neutrino emission rates and the behavior of superfluid and superconducting particles at densities found nowhere else in nature.

The cooling story begins dramatically. A newborn neutron star starts at more than ten billion kelvin, transparent only to neutrinos after roughly a minute, and shrinks from about a hundred kilometers to its final radius of ten to fourteen kilometers. Within the first year, its core — containing about ninety-nine percent of the mass — becomes nearly isothermal thanks to enormous thermal conductivity. Long-term cooling then proceeds over hundreds of thousands of years, with temperature gradients confined to the thin crust and envelope. Crucially, the absence of internal convection means that neutron stars cannot operate the kind of self-sustained dynamos that generate magnetic fields in planets and ordinary stars, so their fields must originate elsewhere and be preserved or dissipated over time.

One of the most consequential recent developments highlighted in the review concerns enhanced neutrino cooling. For decades, theorists have debated whether neutron stars cool via ‘minimal’ channels involving modified Urca processes or whether some stars undergo dramatically faster ‘direct Urca’ cooling, possibly triggered by exotic matter such as hyperons or quarks. The Vela pulsar’s anomalously low temperature hinted at such processes, but definitive evidence was elusive. A recent analysis of three young, nearby and extremely cold neutron stars now shows that their properties require enhanced cooling, allowing researchers for the first time to place meaningful constraints on the dense-matter equation of state directly from thermal observations.

The same cooling physics reaches into particle physics beyond the Standard Model. Because hypothetical light particles such as axions would drain energy from a neutron star’s interior, cooling measurements can bound their properties. Studies cited in the review use the thermal behavior of young neutron stars to constrain the QCD axion mass, and more recent work incorporating envelope structure has tightened those limits further. Meanwhile, the famous case of the neutron star in the Cassiopeia A supernova remnant — once thought to show a rapid temperature decline signaling the onset of neutron superfluidity — illustrates the difficulty of the enterprise. The inferred cooling rate has shrunk from about four percent per decade to between roughly 1.6 and 2.2 percent as data and detector calibrations improved, and recent critical reanalyses question whether the decline is robust at all.

On the magnetic side, the review lays out the three dominant interior processes. Ohmic dissipation, the slow diffusion of currents through the highly conductive crustal lattice, acts on timescales of hundreds of thousands to millions of years. The Hall drift, a nonlinear advection of magnetic field by the electron fluid, becomes dominant in magnetars and can transfer magnetic energy from large scales to small scales, accelerating dissipation. Ambipolar diffusion, the coupled drift of charged particles relative to the neutron superfluid, may dominate in the cores of strongly magnetized stars during their first hundred thousand years, though its exact role in superfluid and superconducting conditions remains actively contested. Newer ingredients — including the chiral magnetic effect, in which a tiny imbalance between left- and right-handed electrons can amplify fields, and crustal failure mechanisms ranging from brittle starquakes to plastic flow — are given detailed treatment, with the authors cautioning that some widely used failure criteria may overestimate how often the crust actually cracks.

Numerical methodology occupies a substantial portion of the review, and for good reason: the equations are stiff, nonlinear and span many orders of magnitude in the relevant coefficients. Older simulations relied on spectral methods elegant in accuracy but fragile near the discontinuities that the Hall term inevitably produces. Recent codes have shifted toward finite-difference and finite-volume schemes borrowed from high-resolution shock-capturing techniques in computational fluid dynamics. The new three-dimensional code MATINS, developed by members of the review team, goes further by adopting a cubed-sphere grid that sidesteps the numerical pathologies of the polar axis, while incorporating realistic microphysics throughout the crust. Benchmark tests against analytical solutions confirm accuracy at the sub-percent level.

What happens when such codes are initialized with realistic, turbulent magnetic fields inherited from the proto-neutron star dynamo phase is striking. Fully three-dimensional simulations starting from complex fields dominated by small and intermediate scales show that the tangled structure persists for hundreds of thousands of years, with the Hall term continuously feeding energy from large scales while Ohmic dissipation erodes the small scales. These turbulent initial conditions naturally reproduce the observed properties of central compact objects and the so-called low-field magnetars — stars that erupt in magnetar-like bursts despite possessing relatively weak large-scale dipoles. However, the simulations still struggle to generate the ultra-strong, dominant dipole that classical magnetars require, leaving the origin of magnetar-strength dipolar fields an open question. One provocative possibility explored in recent work is that the chiral magnetic effect could grow the dipolar component to magnetar strengths within fifty to a hundred years of birth, driven by an extraordinarily small but persistent chiral imbalance sustained by the star’s magnetic helicity.

The final piece of the puzzle is rotation. Spin period and its derivative are the most precisely measured neutron star observables, and they encode the electromagnetic torque exerted by the magnetosphere. The review emphasizes that the familiar textbook formula for inferring magnetic fields from timing data, based on a vacuum dipole, systematically overestimates the true field strength. Plasma-filled force-free magnetospheres exert torques even on aligned rotators, and general relativistic effects near the stellar surface further amplify the spin-down luminosity. A properly relativistic formula yields dipolar field estimates roughly a factor of four lower than the classical expression, potentially reshaping how the entire neutron star population is classified — including how many objects truly qualify as magnetars.

As instruments such as the Square Kilometre Array Observatory prepare to detect many thousands of new pulsars and X-ray observatories continue to monitor magnetar outbursts, the theoretical framework assembled in this review provides the essential bridge between microphysics and observation. The authors’ outlook is clear: the future lies in fully three-dimensional, self-consistent magneto-thermal simulations coupled to dynamical magnetospheres, with machine-learning solvers such as physics-informed neural networks emerging as a computationally efficient complement to traditional schemes. What emerges from the review as a whole is a picture of neutron stars not as static lighthouses but as dynamic, evolving systems in which magnetism, heat and rotation are inseparably entangled — a picture that is only now becoming possible to compute.

Beyond the theoretical machinery, the review situates itself within a rapidly growing observational census. Nearly four thousand rotation-powered radio pulsars are now catalogued, a number set to multiply as the Square Kilometre Array Observatory comes online, while high-energy instruments have identified several hundred gamma-ray pulsars and roughly a hundred X-ray sources. Only a few dozen neutron stars show the soft, thermal surface emission in X-rays that directly probes interior cooling, making each detection disproportionately valuable for constraining the underlying microphysics.

Magnetars, the extreme end of this population, remain scarce — a few dozen confirmed candidates — but their energetic output is extraordinary. Their persistent X-ray luminosities of 10^33 to 10^35 erg per second routinely exceed their rotational energy loss rates by orders of magnitude, implying that magnetic energy, not spin-down, powers their emission. The most violent events, giant flares, release up to 10^46 erg in under a second; only three have ever been recorded. This magnetic reservoir interpretation, first advanced by Thompson and Duncan in the mid-1990s, now underpins most theoretical modeling of magnetar activity.

Perhaps the most significant conceptual shift the review documents is the blurring of historical subclass boundaries. Low-field magnetars — nominally ordinary pulsars that occasionally erupt in bursts — and some central compact objects demonstrate that a strong surface dipole is not the decisive ingredient for magnetar-like behavior. Instead, attention has turned to the hidden internal field architecture: how magnetic energy is partitioned between toroidal and poloidal components and across spatial scales, where electrical currents reside, and how helicity is transported outward into the magnetosphere. These questions set the agenda for the three-dimensional simulations and future code development efforts the review outlines.

Subject of Research: Magneto-thermal and rotational evolution models of isolated neutron stars

Article Title: Magnetic, thermal and rotational evolution of isolated neutron stars

Article References: Pons, J. A., Dehman, C., & Viganò, D. (2026). Magnetic, thermal and rotational evolution of isolated neutron stars. Living Reviews in Computational Astrophysics, 12(1), Article 5. https://doi.org/10.1007/s41115-026-00028-4

Image Credits: AI Generated

DOI: 10.1007/s41115-026-00028-4

Keywords: neutron stars, magnetars, pulsars, magnetic field decay, neutron star cooling, neutrino emission, Hall drift, Ohmic dissipation, ambipolar diffusion, magneto-thermal simulations, dense matter equation of state, numerical astrophysics

The surface of Mercury tells a remarkable story of planetary contraction. Like a grape transforming into a raisin, the planet has wrinkled as its interior slowly cooled. When a planet’s interior loses heat, its rocky outer shell is forced to accommodate the shrinking volume beneath it. The crust responds by crumpling, forming dramatic geological features such as ridges, cliffs, and steep scarps that scar the planet’s surface. These features, known to scientists as shortening structures, are direct evidence of the immense compressive forces that have acted on Mercury’s crust as the planet gradually lost its internal heat.

For decades, planetary scientists have relied on mapping these shortening structures to estimate exactly how much Mercury has shrunk. The logic is elegant: the more the crust has crumpled and folded, the more the planet’s radius must have decreased. By measuring the total accumulated strain expressed in these tectonic features, researchers can work backward to calculate the planet’s loss of radius. Previous estimates based on this method, informed by data from NASA’s MESSENGER mission, suggested that Mercury’s radius had decreased by approximately 8.3 kilometers since the planet’s early history.

However, the new research reveals a fundamental flaw in that approach. According to the study, published in the journal Geophysical Research Letters, the visible record of Mercury’s tectonic activity is far from complete. The planet’s rough and pockmarked surface, bombarded by impacts over billions of years, has obscured many of the very structures scientists use to measure its contraction. In essence, the planet’s own violent history has hidden the evidence of its shrinkage.

“Mercury’s surface preserves a record of how the planet has cooled and contracted, but we found that this record is incomplete,” explains Gaku Nishiyama, lead author of the study and a researcher at Hokkaido University. “We compared a global map of Mercury’s surface roughness with maps of shortening structures and contraction. Once we account for the effect of rough terrain, Mercury appears to have shrunk considerably more than what the visible tectonic record alone suggested.”

The research team’s methodology centered on a systematic comparison between surface roughness and the distribution of tectonic features across the entire planet. By overlaying global maps of surface texture with maps documenting the locations of ridges and scarps, the researchers identified a striking pattern: shortening structures were far more likely to be identified in smooth terrain than in rough regions. This correlation was too strong to be coincidental. Instead, it pointed to a systematic bias in the geological record, where the roughest terrain was concealing tectonic structures from view.

A particularly compelling piece of evidence emerged when the researchers examined areas surrounding relatively young impact craters. When a large asteroid or comet strikes a planetary surface, the explosion excavates material from beneath the crust and throws it outward in a debris field known as ejecta. These rough, jumbled deposits can blanket the surrounding landscape, burying older geological features beneath layers of impact debris. Around the Rachmaninoff crater, a prominent impact basin on Mercury, the researchers observed that shortening structures were noticeably less common in areas covered by rough ejecta. Additionally, the structures became progressively less visible closer to the crater itself, suggesting they were partially buried under impact-related material. This spatial pattern provided a natural experiment, demonstrating exactly how impact debris can erase tectonic evidence from the scientific record.

By mathematically correcting for this obscuration effect, the researchers recalculated Mercury’s total radial contraction. The revised estimate increased from the previously accepted value of 8.3 kilometers to approximately 11.6 kilometers, an increase of roughly 10 to 30 percent compared to earlier figures. Importantly, the team acknowledges that even this new number may still represent an underestimate. Because the relationship between surface roughness and hidden tectonic structures cannot be perfectly quantified, the true magnitude of Mercury’s shrinkage could be even greater than the revised calculation suggests.

The implications of this finding extend far beyond Mercury itself. The same principle—that rough surfaces can obscure tectonic records and lead scientists to underestimate planetary contraction—could apply to other rocky worlds in our solar system. The Moon, whose surface is even rougher than Mercury’s due to its long history of impacts, is a prime candidate for similar reassessment. If the lunar tectonic record has been similarly compromised by impact ejecta and rough terrain, scientists may need to revise estimates of how much the Moon has contracted as its interior cooled. The finding could also influence our understanding of other airless, cratered bodies such as Mars, whose ancient surface preserves billions of years of geological history buried under layers of impact debris.

The timing of this research is particularly significant given the upcoming BepiColombo mission, a joint endeavor between the European Space Agency and the Japan Aerospace Exploration Agency (JAXA). Currently en route to Mercury, the spacecraft carries a sophisticated instrument called the BepiColombo Laser Altimeter, or BELA, which will measure the planet’s surface topography with unprecedented precision. BELA is expected to characterize surface roughness at much finer scales than previous instruments, potentially revealing subtle relationships between younger geological events and tectonic structures that have been difficult to detect. The new findings from Hokkaido University provide a conceptual framework that will help scientists interpret the data BepiColombo sends back, ensuring that the mission’s measurements can be used to build a more complete picture of Mercury’s thermal evolution.

The study also touches on deeper questions about planetary cooling processes. All terrestrial planets—Mercury, Venus, Earth, and Mars—have been losing heat since their formation roughly 4.5 billion years ago. The rate and magnitude of that cooling, expressed through planetary contraction, provides crucial insights into a planet’s internal structure, its thermal history, and even its potential for past geological activity such as volcanism. Mercury, lacking a substantial atmosphere and plate tectonics, preserves its tectonic record with exceptional fidelity, making it an ideal laboratory for studying how rocky planets evolve. But as this study demonstrates, even the best-preserved records contain gaps, and accounting for those gaps can fundamentally change our understanding of a planet’s history.

For the scientific community, the research serves as a cautionary tale about the limitations of surface observation. When scientists estimate planetary contraction from visible tectonic features, they are essentially reading an incomplete book. Pages have been torn out by impacts, buried by ejecta, and hidden by rough terrain. The Hokkaido University team’s approach—systematically comparing terrain characteristics with tectonic feature distribution—offers a template for identifying and correcting such biases on other worlds. As missions to Mercury, the Moon, and Mars continue to return increasingly detailed data, methods like this will be essential for transforming raw observations into accurate reconstructions of planetary history.

Mercury, it turns out, has been hiding its true transformation in plain sight. The little planet that has shrunk so dramatically over billions of years has done so even more extensively than anyone realized. With future observations from BepiColombo poised to refine these measurements further, scientists are now one step closer to understanding the full story of how Mercury—and potentially other rocky worlds—cooled, contracted, and evolved into the planets we see today.

The research team, led by Xiangtong Wang and Wei Li of Sichuan University together with colleagues from the China Academy of Space Technology and other institutions, focuses on the two dynamics that most severely degrade routing in operational low-Earth-orbit systems. The first is the persistent churn of the satellite-ground topology: as spacecraft orbit and rotate with the planet, the map of who can talk to whom is rewritten continuously. The second is the intermittent failure of inter-satellite links, the point-to-point connections that carry traffic between neighboring satellites. Each of these effects is troublesome on its own; together, the authors argue, they simultaneously flood the network’s control plane with global update traffic and carve the connectivity graph into sparse, fragmented islands, degrading both the speed of route computation and the network’s ability to survive faults.

Conventional routing protocols were simply not built for this environment. Interior gateway protocols such as OSPF assume a relatively stable topology and rely on network-wide flooding of link-state updates whenever anything changes. In a mega-constellation, where link failures may occur somewhere in the network almost constantly, that flooding translates into an enormous and perpetual control overhead, consuming the very satellite bandwidth and processing capacity needed for user traffic. Geographic routing approaches, which forward packets based on physical position, avoid some of this overhead but can struggle with the polar regions, seam discontinuities and the sparse connectivity left behind by link failures. The BlockFLEX design borrows strategically from both traditions while insulating them from the worst of the churn.

The core idea is a robust virtual overlay built from anonymous blocks. Instead of treating every satellite as an individually addressable, constantly moving node, BlockFLEX groups satellites into clusters whose membership changes are hidden behind a stable abstraction. To the routing layer above, each block behaves as a coherent unit, masking the underlying motion of individual spacecraft and the comings and goings of their links. This is the same general principle that has long been used in mobile ad hoc networks, where clustering reduces the frequency with which routes must be recomputed, but the authors refine it for the specific geometry and failure patterns of LEO constellations, ensuring that the overlay remains stable even as its members change.

On top of this overlay, BlockFLEX deploys a two-tier hybrid routing strategy that splits responsibilities by scale. Between blocks, the architecture relies on convergence-free geographic forwarding: packets are passed from block to block based on spatial position, a method that requires no global route computation and therefore never stalls waiting for routing tables to converge after a change. Within each block, by contrast, a conventional convergence-isolated routing protocol manages fine-grained delivery. Because each block contains only a small fraction of the constellation’s satellites, the control traffic needed to maintain routes inside it stays confined to that block, localizing the propagation of updates and preventing a single link failure from triggering a constellation-wide signaling storm.

This division of labor is where the architecture derives most of its resilience. When an inter-satellite link fails, the blast radius of the disruption is contained. Packets in transit between blocks are rerouted by the geographic layer without any protocol negotiation, while the affected block recomputes its internal paths quietly and locally. The authors also incorporate complementary mechanisms aimed at scalability and efficiency, ensuring that the approach continues to perform as constellations grow from hundreds toward many thousands of satellites. The design draws on established tools such as graph partitioning for block construction and satellite positioning knowledge, including GPS-based ephemeris data, to keep the geographic layer accurate as orbits evolve.

The evaluation is notable for being grounded in models of currently operational LEO mega-constellations rather than purely synthetic topologies, and the team has released its source code and data publicly on GitHub to allow independent reproduction. In stress tests that randomly disabled up to thirty percent of all inter-satellite links, a failure rate far beyond what most operational scenarios would expect, BlockFLEX substantially outperformed state-of-the-art routing schemes in both resilience and efficiency. In practical terms, that means a higher fraction of packets still reach their destinations when links fail, and the network spends far less of its capacity on control-plane bookkeeping instead of user traffic.

The stakes of this engineering problem are rising quickly. Starlink already operates thousands of satellites with laser inter-satellite links, OneWeb and Amazon’s Kuiper are deploying their own fleets, and Chinese constellations such as Guowang and Qianfan are moving toward large-scale service. These systems aim to deliver broadband to aircraft, ships, rural communities and military users, and their commercial viability depends on squeezing maximum throughput out of expensive orbital infrastructure. Researchers have previously shown that adversarial or accidental failures can degrade LEO network performance dramatically, and studies of attack surfaces such as ICARUS have highlighted how exposed these networks can be. A routing layer that degrades gracefully under heavy link loss is therefore not merely an academic nicety but a competitive necessity.

BlockFLEX also fits into a broader architectural conversation about how the space network of the future should be organized. Some proposals advocate software-defined satellite networking, in which a centralized controller computes all routes; others favor entirely distributed or purely geographic schemes; standards bodies have begun formalizing routing architectures for satellite networks in documents such as recent RFCs on satellite routing. The hybrid, hierarchical philosophy behind BlockFLEX suggests a middle path: keep the global layer dumb and stable, and push the dynamic, failure-sensitive work to small, self-contained neighborhoods where it can be handled quickly and cheaply. That philosophy echoes design lessons from the terrestrial Internet itself, where hierarchy and locality have historically been the keys to scaling.

Open questions remain. The anonymous block abstraction must be maintained somehow, and the cost of regrouping satellites as orbits evolve, particularly across the constellation’s seams and polar regions, will matter at full scale. Interoperability with existing ground segment operations and with the proprietary scheduling algorithms used by commercial operators is another practical hurdle. Nonetheless, the demonstration that a hybrid block-based design can withstand the loss of nearly a third of its inter-satellite links while outperforming existing schemes marks a meaningful step toward routing infrastructure that is genuinely worthy of the mega-constellation era. As humanity wires the sky, the networks that survive will be the ones designed, like BlockFLEX, to expect failure as the norm rather than the exception.

Beyond the headline results, the study sits within a well-established lineage of research on grouping nodes in highly dynamic networks. Clustering was first explored in the late 1990s for reconfigurable wireless and mobile ad hoc networks, where grouping mobile nodes reduced how often routes needed repair. Applying that principle to orbiting spacecraft is attractive because the predictable geometry of satellite orbits makes membership changes computable in advance, unlike the random motion of terrestrial mobile devices.

The geographic forwarding tier likewise draws on decades of prior work. Early proposals for LEO constellations, including datagram routing algorithms from the INFOCOM literature of 2000 and geographic packet routing demonstrated in experimental satellite systems, showed that physical position is a remarkably stable routing signal in orbit, since satellite motion is deterministic and ephemeris data can be predicted far into the future. BlockFLEX inherits this determinism: the coarse decision about which block should receive a packet next can be made from positional knowledge alone, requiring no signaling exchange across a shaken topology.

Another contextual thread is the failure profile of the links themselves. Inter-satellite links in current constellations rely on free-space optical terminals with narrow beams, which must continuously acquire, track and point at one another as relative geometry changes. Acquisition and pointing errors, thermal effects and blockage all contribute to intermittent outages, making a thirty percent random-failure stress test a plausible extreme rather than pure fiction. The authors’ decision to evaluate on models of operational constellations, alongside publicly released code, positions the work for scrutiny as deployment scales and real-world link telemetry accumulates.

Subject of Research: Efficient and resilient routing for low-Earth-orbit mega-constellation satellite networks

Article Title: Towards an Efficient, Resilient and Scalable Routing for Mega-Constellation Networks

Article References: Wang, X., Li, W., Yang, M., Han, S., Jiang, M., & Wang, Y. (2026). Towards an Efficient, Resilient and Scalable Routing for Mega-Constellation Networks. International Journal of Aeronautical and Space Sciences. https://doi.org/10.1007/s42405-026-01289-2

Image Credits: AI Generated

DOI: 10.1007/s42405-026-01289-2

Keywords: mega-constellation networks, LEO satellites, inter-satellite link failures, hybrid routing, routing resilience, virtual overlay, geographic forwarding, network scalability, BlockFLEX, satellite communication, control-plane overhead, Sichuan University

The raw material in question is lunar regolith, the powdery, shattered blanket of rock that covers the Moon’s surface. Forged by billions of years of micrometeoroid impacts and space weathering, regolith is a complex mixture of silicate minerals and metal oxides, including silicon dioxide, aluminum oxide, calcium oxide, iron oxide, titanium dioxide, and magnesium oxide. Locked inside these oxides are exactly the elements a lunar base needs most: oxygen for breathing and rocket propellant, and metals for construction, fabrication, and power infrastructure. But the chemistry is unforgiving. These elements are bound tightly in mineral matrices, and liberating them demands either extreme heat or aggressive electrochemical processing. The Moon’s own environment compounds the challenge, offering high vacuum, reduced gravity, abrasive dust, and a fourteen-day night that strains any power system.

Where you dig matters enormously. The dark volcanic plains of the lunar maria are enriched in iron, titanium, and magnesium, making them attractive for processes that target those oxides. The brighter highlands, by contrast, are dominated by anorthosite, a rock built largely from the aluminum- and calcium-rich feldspar anorthite. Because aluminum is prized for structural applications, electrical conductors, and recyclability, the research team focused its comparison on anorthite as a representative highlands feedstock, giving all three extraction methods an identical starting material so that differences in performance could be attributed to the processes themselves rather than to geochemistry.

The first contender is carbothermal reduction, or CTR, a high-temperature chemical process that uses carbon, typically as graphite or methane, to strip oxygen from metal oxides. In theory, certain oxides can be reduced near 850 degrees Celsius, but practical operation generally demands temperatures above 1500 degrees Celsius, where partial melting of the feedstock improves reactant mobility and conversion. CTR splits into two regimes: a solid-gas pathway below roughly 1200 degrees Celsius with slower kinetics and incomplete conversion, and a molten-phase pathway above about 1400 degrees Celsius that delivers more extensive reduction. The process is especially attractive for silicon and iron production from iron-rich mare regolith, and it pairs naturally with concentrated solar thermal heating, reducing the demand for electrical power. Its Achilles heel is carbon. The Moon has essentially none, so large-scale CTR would require either constant resupply from Earth or in-situ carbon production, both of which add logistical complexity. Carbon can also react with silicon products to form silicon carbide, threatening both product purity and reactor longevity. Prior studies of carbothermal reduction on lunar regolith report oxygen yields of only 2.4 to 7.9 percent of the feedstock mass, corresponding to just 5.4 to 17.5 percent of the total available oxygen.

The second contender is molten salt electrolysis, or MSE, an electrochemical approach in which metal oxides are reduced inside a molten salt electrolyte, most commonly calcium chloride. Operating at a comparatively moderate 800 to 1000 degrees Celsius, MSE follows the logic of the FFC Cambridge process, in which oxygen ions are pulled directly out of solid oxides. Metal cations migrate to the cathode and deposit as metal, while oxygen gas evolves at the anode. The study highlights a lunar-adapted configuration called LISAP-MSE, for Lunar In-Situ Aluminum Production via Molten Salt Electrolysis, which first leaches anorthite with hydrochloric acid, thermally decomposes the resulting aluminum chloride hydrate into alumina, and then electrolytically reduces that alumina to metallic aluminum and oxygen. The calcium chloride electrolyte may even be derivable from processed regolith itself. MSE’s great strength is selectivity: because only electrochemically active oxides are reduced, beneficiation steps such as acid leaching, flotation, or magnetic separation can be used upstream to boost feedstock quality and deliver high-purity metals. The trade-offs are system complexity, slower oxygen production rates, and the need to manage electrode degradation and electrolyte stability over long durations in a dusty environment.

The third contender, molten regolith electrolysis, or MRE, is the most direct of all. Rather than beneficiating feedstock or adding an electrolyte, MRE simply melts raw regolith above roughly 1300 degrees Celsius, where the silicate melt becomes ionically conductive, and applies an electric potential across it. Metal cations such as iron, titanium, silicon, aluminum, and calcium are reduced at the cathode into mixed alloys, while oxygen streams off at the anode. This brute-force simplicity makes MRE appealing for early-stage missions where mass efficiency matters and oxygen is the priority, and it can in principle process the complete range of oxides present in lunar soil. But the same extreme temperatures impose severe engineering burdens: refractory containment materials must resist corrosion by molten silicate, electrodes must survive repeated thermal cycling and oxidation, and continuous power delivery is essential to keep the melt conductive and the reaction progressing.

To compare the three fairly, the researchers built a rigorous thermodynamic framework using temperature-adjusted enthalpy, entropy, and Gibbs free energy, drawing on NIST thermochemical data supplemented by specialized datasets for anorthite, steam, and graphitic carbon, whose high-temperature heat capacities required alternative formulations. Each process was modeled as complete conversion of anorthite to its products, establishing an upper-bound theoretical limit for aluminum and oxygen extraction. The results were normalized per kilogram of aluminum and oxygen produced, allowing a direct, first-principles comparison of inherent energy requirements at each process’s realistic operating temperature.

The verdict was unambiguous. All three processes are endothermic and require external energy input, but LISAP-MSE posted the lowest enthalpy and Gibbs free energy values of the three, making it the most energetically efficient route to simultaneous aluminum and oxygen production under the study’s assumptions. CTR landed in the middle, with higher thermal demands compounded by its consumable carbon feedstock and logistical overhead. MRE fared worst on energy, displaying the highest enthalpy and Gibbs free energy values, a reflection of the enormous cost of maintaining a fully molten silicate phase and driving electrochemistry within it. Notably, none of the reactions were thermodynamically spontaneous at the evaluated temperatures and one bar pressure, confirming that all of these processes will depend on sustained external heating and electrical power. The authors caution, however, that the ranking is framework-dependent: if bulk oxygen were the dominant goal, or if iron- and titanium-rich mare regolith were the feedstock, CTR and MRE could climb back up the rankings, and hydrogen reduction, excluded from this comparison, might also become attractive.

Perhaps the study’s most consequential insight is that the three technologies are not competitors at all, but natural partners. CTR and MRE share similar high-temperature operating conditions and produce high-throughput oxygen along with mixed-metal alloys of limited purity, while MSE operates cooler and refines selectively. The researchers sketch integrated architectures in which calcium aluminate byproducts from CTR and MRE are fed into an MSE module as beneficiated feedstock, yielding a high-purity calcium-aluminum alloy that then acts as a metallothermic reducing agent to liberate additional silicon and iron from leftover oxides, with the reformed calcium aluminates looping back to regenerate the alloy and release more oxygen. A second synergy uses the MSE-produced alloy as a pre-treatment for MRE, preferentially displacing iron, magnesium, and titanium from the melt to leave a silicon-enriched mixture that electrolyzes more cleanly. On the infrastructure side, the authors map out a phased concept of operations: deploy modular reactors with imported spares, pilot CTR or MRE for early oxygen, integrate autonomous excavation and hauling at rates of at least 2.54 kilograms of regolith per hour for a modest 10,000-kilogram-per-year oxygen plant, then scale up MSE modules as aluminum demand for construction and shielding grows. Power infrastructure must reach tens to hundreds of kilowatts, with nuclear fission surface power flagged as the most practical near-term backbone, supplemented by solar arrays during the lunar day. Autonomy will be essential given communication delays, and dust mitigation, through sealed pneumatic transfer, dust-tolerant seals, and electrodynamic dust shields, will decide whether reactors survive the abrasive lunar surface. The team’s recommendations for future research center on three fronts: materials durability in molten environments, continuous autonomous feedstock handling, and system-level integration of thermal, electrical, and mass flows. The message is clear: no single process can carry a lunar economy alone, and the first self-sustaining Moon base will be built on the clever coupling of all of them.

Subject of Research: Comparative thermodynamic analysis of carbothermal reduction, molten salt electrolysis, and molten regolith electrolysis for extracting aluminum and oxygen from lunar regolith

Article Title: A review and analysis of extraction methods for aluminum and oxygen from lunar regolith

Article References: Ortega, J. N., Rupp, B., & Han, F. D. (2026). A review and analysis of extraction methods for aluminum and oxygen from lunar regolith. Space and Planetary Resources, 2(1), Article 5. https://doi.org/10.1007/s44461-026-00010-1

Image Credits: AI Generated

DOI: 10.1007/s44461-026-00010-1

Keywords: lunar regolith, in-situ resource utilization, aluminum extraction, oxygen production, carbothermal reduction, molten salt electrolysis, molten regolith electrolysis, anorthite, lunar ISRU, space resources, thermodynamics, Moon base

Magnetic twist in the solar corona is no exotic hypothetical. Observations from instruments aboard TRACE and SOHO have revealed rotating sunspots, braided loop structures and helical signatures in extreme-ultraviolet images, all pointing to the presence of azimuthal magnetic field components superimposed on the dominant longitudinal field threading coronal loops. Such twisted configurations are expected to arise when magnetic flux tubes emerge from beneath the photosphere, and they are intimately linked to some of the Sun’s most dramatic behavior, including kink instabilities that may trigger flares. Yet when it comes to modeling wave propagation in these structures, most theoretical studies have relied on simplifications—either treating the twist as confined to a thin surface layer or ignoring it entirely. The Shandong team set out to relax that simplification, constructing a straight magnetic cylinder in which both an axis-aligned longitudinal field and a distributed azimuthal field coexist throughout the equilibrium.

The technical heart of the study lies in the derivation of the dispersion relation for the fast sausage modes, obtained by solving the eigenvalue problem that emerges from the one-dimensional resistive magnetohydrodynamic equations. Resistive MHD, which incorporates a finite electrical resistivity into the otherwise idealized description of conducting plasma, is essential here for a subtle reason: in a continuous equilibrium distribution, the natural resonant surfaces—locations where the wave’s phase speed matches a local Alfvén or continuum speed—cannot be treated as sharp discontinuities. Instead, the resistive formulation allows the singular behavior at these resonant positions to be handled numerically, following a long tradition in magnetohydrodynamic wave theory that stretches back to classic work on continuous spectra of cylindrical MHD equilibria in the 1970s. The authors employed the numerical code PDE2D to solve the resulting boundary value problem, carefully choosing the outer computational boundary to be far enough away that it did not contaminate the eigenfrequencies—a methodological point they address in detail in an appendix, demonstrating that the numerical results become insensitive to the outer boundary once it extends beyond roughly ten loop radii.

One of the most consequential findings of the study is a confirmation of a previous result with far-reaching observational meaning: the principal fast sausage mode exhibits no cutoff. In wave physics, a cutoff is a critical wavelength or wavenumber below which disturbances cannot propagate as trapped oscillations within the guiding structure, instead leaking energy into the surrounding environment. For ordinary fast sausage modes in straight, untwisted cylinders, such cutoffs exist for all but the fundamental mode, meaning that only perturbations with sufficiently short wavelengths remain trapped within the loop. The new calculation shows that when the magnetic field is twisted with a continuous distribution, the principal sausage mode retains its ability to propagate at arbitrary wavelengths. In practical terms, this means the fundamental mode can carry information along a twisted coronal loop regardless of how long the wavelength of the disturbance is, a property that greatly expands the range of oscillations that seismologists might hope to detect in solar observations.

The paper goes further, systematically mapping how the oscillation frequencies and damping rates of fast sausage modes depend on two key parameters: the strength of the magnetic twist and the longitudinal wavenumber, which characterizes how rapidly the wave varies along the length of the loop. The analysis reveals that twist modifies the dispersive properties of the modes, shifting the relationship between frequency and wavelength in ways that depend on the detailed equilibrium structure. These dependencies matter because they form the theoretical foundation of coronal seismology—the technique of inferring otherwise unmeasurable plasma and magnetic parameters of the corona from observed oscillation properties, in much the same way that geologists infer the Earth’s interior structure from seismic waves. If oscillation periods, damping times and phase speeds can be measured precisely, the theoretical dispersion relations allow researchers to work backward to quantities such as magnetic field strength, density contrast and, in this case, the degree of twist itself.

Resonant absorption, the mechanism by which wave energy is transferred from the fast sausage oscillation into the continuous Alfvén spectrum at resonant surfaces, was a particular focus of the investigation. This damping mechanism has been extensively studied in the context of kink oscillations of coronal loops, where it is widely regarded as a leading candidate for explaining the rapid decay of transverse loop oscillations observed after flares and coronal mass ejections. For the fast sausage modes in twisted loops considered here, however, the outcome is strikingly different. For the flux tube parameters explored in the study, the damping induced by resonant absorption turns out to be far too weak to produce observable signatures in actual solar observations. The damping rates fall below the threshold of detectability for current instrumentation, suggesting that if observed sausage oscillations are found to decay rapidly, some other damping mechanism—or some other structural property of the loop—must be responsible.

This result carries a cautionary message for observers. Quasi-periodic pulsations with periods ranging from fractions of a second to several minutes are routinely detected in solar flare emission across radio, extreme-ultraviolet and hard X-ray wavelengths, and fast sausage modes are frequently invoked as one of the physical interpretations, alongside mechanisms involving MHD oscillations of the flare loops, periodic reconnection, or wave-particle interactions. Spatially resolved microwave observations of flare loops have revealed oscillation patterns consistent with standing sausage modes, and these interpretations have been used to diagnose flare loop parameters such as magnetic field strengths and density scales. The new study does not invalidate such diagnostics—indeed, the confirmed absence of a cutoff for the principal mode strengthens the theoretical basis for interpreting long-wavelength, global oscillations in terms of the fundamental sausage mode. But it does suggest that damping-based diagnostics, which would exploit the temporal decay of sausage signals to infer transverse density structuring, may be futile for twisted loops of the type modeled, since the resonant damping is simply too feeble to matter.

The modeling framework itself deserves attention. By treating the coronal loop as a structured straight cylinder with a continuous equilibrium distribution—meaning the plasma density, longitudinal field and azimuthal field all vary smoothly with radius rather than jumping abruptly at the loop boundary—the authors sidestep the mathematical pathologies associated with sharp discontinuities while retaining the essential physics of magnetic twist. The azimuthal field component, which is what distinguishes a twisted loop from a simple potential field, enters the wave equations through the equilibrium magnetic field curvature and the associated current distribution, coupling the compressive sausage motions to the torsional degrees of freedom. Solving the full eigenvalue problem numerically, rather than relying on asymptotic approximations valid only in limiting regimes, allows the study to capture the behavior of the modes across the full range of longitudinal wavenumbers, from the long-wavelength fundamental regime through the cutoff wavenumbers of the higher overtones.

The broader significance of the work lies in its contribution to a half-century-long effort to understand MHD waves in structured solar plasmas. Since the first suggestions in 1970 that MHD pulsations might be observable in the corona, and the landmark theoretical treatment of wave propagation in magnetic cylinders in 1983, sausage modes have been studied in configurations of increasing realism: with transverse density structuring, with surface currents, with asymmetric environments and, as here, with magnetic twist. Each increment in realism sharpens the tools available for interpreting the growing archive of high-resolution solar observations from missions such as SDO, IRIS and the Daniel K. Inouye Solar Telescope. The Shandong team’s results—confirming the cutoff-free nature of the principal mode while demonstrating the observational insignificance of resonant damping in their parameter regime—provide both reassurance and a boundary marker, telling the seismology community which inferences can be trusted and which lines of diagnostic reasoning are unlikely to bear fruit.

For now, the twisted loops of the solar corona keep their secrets reluctantly. But studies of this kind, by mapping precisely how waves encode the properties of their host structures, are steadily turning the Sun’s oscillating atmosphere into a readable text—one pulsation at a time.

To appreciate the significance of the work, it helps to understand what Lorentzian distance measures. In general relativity, spacetime is modeled as a four-dimensional manifold equipped with a Lorentzian metric, a geometric object that defines causal relationships between events. The Lorentzian distance between two causally related points is the length of the longest timelike curve—the worldline of a permissible observer—connecting them. Unlike ordinary Riemannian distance, which is smooth and well-behaved, Lorentzian distance is only continuous in general, and it vanishes whenever the two points cannot be causally connected. This inherent roughness makes Lorentzian geometry a delicate subject, and many of its deepest theorems rely on careful control of how these distance functions behave.

The classical Lorentzian splitting theorems, which trace back to work in the late 1980s by researchers including Richard Newman, J.-H. Eschenburg, and Galloway himself, are a case in point. These theorems assert, roughly, that if a spacetime satisfies the strong energy condition—the statement that gravity, as encoded in the Einstein equations, tends to focus matter—and contains a complete timelike line, meaning an inextendible geodesic that maximizes the time separation between every pair of points along it, then the spacetime splits as a product. In physical terms, the universe must decompose into a static spatial part crossed with ordinary time, much like the flat spacetimes of special relativity. Splitting theorems are rigidity statements: they show that under physically reasonable hypotheses, the universe cannot have an exotic global shape.

A crucial ingredient in all known proofs of these splitting theorems is the local equi-Lipschitz continuity of families of Lorentzian distance functions—often constructed as Busemann functions—associated with a complete timelike line. In essence, equi-Lipschitzness means that an entire family of functions shares a common bound on how fast they can change: there is a single Lipschitz constant that controls them all in a neighborhood of the line. This uniform control allows mathematicians to extract smoothly converging subsequences, take limits of Busemann functions, and ultimately produce the preferred time coordinate whose existence forces the splitting. Without it, the limiting machinery at the heart of the proofs would collapse.

Galloway, McCann, and Ohanyan asked a natural next question: does the same uniform regularity hold when the distance functions are taken not to a single point, but to the level sets of a Cauchy temporal function? A Cauchy temporal function is a smooth time function that increases along every future-directed timelike curve and whose level sets are Cauchy surfaces—spatial slices that every inextendible timelike curve crosses exactly once. Such functions are the gold standard for imposing a clean global notion of time on a spacetime, and they exist in all globally hyperbolic spacetimes, the class of spacetimes most physicists consider reasonable models of our universe. One might therefore expect the distances to these well-behaved foliations to inherit the same regularity as distances to points.

The authors show that this expectation is false. In general, families of Lorentzian distances to and from the level sets of a Cauchy temporal function fail to be locally equi-Lipschitz. The failure is not an artifact of pathological metrics or exotic causal structures; it is an intrinsic feature of how Lorentzian geometry treats spatial slices. Whereas a complete timelike line provides a rigid scaffold—its maximizing property propagates uniform control through the neighborhood—the foliation by Cauchy surfaces offers no such mechanism. The distance to a surface can concentrate its variation in ways that no single Lipschitz constant can tame, even locally. This negative result matters because several proposed approaches to splitting conjectures for Cauchy surfaces have implicitly assumed or hoped for exactly this kind of regularity, and the new theorem rules out the naive strategy of extending the classical point-based arguments wholesale to foliations.

The stakes become clear when the authors connect their findings to Bartnik’s splitting conjecture. In 1988, Robert Bartnik conjectured that a cosmological spacetime—one that admits a compact Cauchy surface—satisfying the strong energy condition should admit a splitting of a related kind, with deep ties to the existence of constant mean curvature surfaces. Decades of partial progress have been recorded; the authors note in their notes that specific classes of Cauchy surfaces have been treated in earlier works, yielding splitting results subject to additional conditions. But a general proof remains elusive. Galloway, McCann, and Ohanyan formulate new conjectures based on the existence of Cauchy temporal functions in cosmological spacetimes and in timelike geodesically complete spacetimes, conjecturing that the Lorentz distances to the level sets of such functions are equi-Lipschitz precisely in the circumstances that matter. Strikingly, they prove that these conjectures are equivalent to Bartnik’s splitting conjecture.

This equivalence is the conceptual heart of the paper. It transforms an analytic question about the regularity of distance functions into a statement about the global causal structure of the cosmos, and vice versa. On the one hand, if Bartnik’s conjecture holds, then the equi-Lipschitz property must follow in the conjectured settings, providing a new handle on spacetime rigidity. On the other hand, anyone seeking to prove Bartnik’s conjecture can now target the equi-Lipschitz property directly, using tools from analysis, partial differential equations, and geometric measure theory that were developed for related regularity problems. The authors’ negative theorem serves as a warning sign along the way: the property does not come for free, so any successful proof must identify precisely which additional structure restores uniform control.

The work also sits within a broader modern program to extend Lorentzian geometry beyond smooth manifolds. Recent research, including the authors’ own collaborations with Mathias Braun, Nicola Gigli, Clemens Sämann, and others, has developed splitting theorems and causal calculus on nonsmooth and metric-measure spacetimes, inspired in part by optimal transport and by the synthetic treatment of curvature. In such settings, where the metric may be only continuous or the spacetime may carry a measure-theoretic weight, questions of Lipschitz regularity become even more delicate. The new negative result calibrates expectations across this program, indicating that techniques tied to distance functions along lines—which do enjoy local equi-Lipschitzness—cannot be transplanted blindly to settings built around temporal functions and their foliations.

For physicists, the paper is a reminder that the geometry underlying cosmological models is governed by rigid mathematical constraints that remain only partially mapped. Splitting theorems constrain the possible global shapes of universes obeying the energy conditions that dominate classical cosmology, and Bartnik’s conjecture sits at the frontier of this understanding. By demonstrating exactly where the standard toolkit fails—and by proving that a specific regularity conjecture is equivalent to that frontier problem—Galloway, McCann, and Ohanyan have redrawn the map of what is known and what must be proven. The full paper, including the precise counterexamples and the proofs of the equivalences, is available in General Relativity and Gravitation.

The starting point is a thought experiment that has haunted quantum foundations since Eugene Wigner first posed it in the 1960s. Wigner’s friend sits inside a sealed laboratory and measures a quantum system, obtaining a definite outcome. Wigner, standing outside, treats the entire laboratory — friend and system included — as a single quantum system in an entangled superposition. The friend says the outcome is definite; Wigner’s quantum formalism says no outcome yet exists. The recent no-go theorems, developed by researchers including Časlav Brukner, Daniela Frauchiger and Renato Renner, and a team led by Kok-Wei Bong, sharpen this paradox into rigorous impossibility results. They target what philosophers call the Absoluteness of Observed Events, or AOE: the claim that any measurement outcome has a unique, definite value that does not depend on which observer measured it. The theorems show that AOE cannot be reconciled with quantum theory plus other reasonable-sounding assumptions about locality and free choice.

Pienaar’s contribution begins with a careful dissection of what AOE actually asserts. It bundles together two distinct claims: that an outcome is unique, and that it is absolute — the same for everyone. One can reject the package either way. Many-worlds interpretations reject uniqueness: every outcome occurs, each in its own branch, and there is arguably no single classical space-time threading them together, only something like a “block multiverse.” But a second family of interpretations, which Pienaar calls perspectival, retains uniqueness while denying absoluteness. Each outcome really happens, with a single definite value — but only relative to the observer who performs the measurement. Prominent members of this camp include QBism, the quantum-Bayesian approach developed by Christopher Fuchs and colleagues, and Carlo Rovelli’s Relational Quantum Mechanics.

For these perspectival views, the question Pienaar tackles is brutally simple to state: can observers who reject absoluteness still agree that they live inside one shared space-time? The concern was voiced most crisply by Eric Cavalcanti in a 2021 paper on QBism’s response to Wigner’s friend. If we reject AOE, Cavalcanti argued, then the classical notion of an event must be challenged too: outcomes that are definite for the friend but not for Wigner cannot be located in “Wigner’s space-time” at all. They occur, in effect, inside a “Wigner bubble” — a region of reality that exists for one observer but not for another. Fuchs himself has long embraced a related sentiment, describing space-time in QBism as an “abstract diagram” each agent uses to organize their own expectations, and writing in 2011 that the lesson of Wigner’s friend is that the world is more truly a “pluriverse” than a universe.

Yet Cavalcanti’s argument, Pienaar noticed, contains a hidden weakness. It works only if the space-time point where the friend’s outcome occurs genuinely cannot be matched to any point in Wigner’s manifold. In a seemingly devastating counterexample — which Pienaar names “Wigner’s diamond” — that matching looks easy. Inside a sealed, quantum-controlled chamber, the friend opens a box containing a diamond whose nitrogen-vacancy center encodes a qubit. She measures it by shining a green laser, at a wavelength of roughly 637 nanometers, onto the NV-center and watching for fluorescence: light on means outcome “1,” no light means “0.” Wigner watches through a window. A tiny screen blocks his view of the diamond, so no information about the fluorescence escapes, and Wigner’s entangled superposition survives intact. But the chamber is fitted with rigid measuring rods and embedded clocks — a physical coordinate grid — so Wigner can see exactly where the diamond sits and exactly when the laser fires. Both observers agree on the where and the when. The only disagreement is whether an outcome with a definite value actually happened there. For the friend it did; for Wigner the light field is entangled with the friend and there is simply no event to speak of. Does this not defeat the bubble conclusion?

The answer, Pienaar argues, hinges on deep and contested questions in the philosophy of space-time — specifically on Einstein’s famous “hole argument” from general relativity. General relativity’s equations are invariant under smooth transformations of the coordinate system called diffeomorphisms, which means that radically different-looking distributions of matter and metric fields over the space-time manifold can be observationally identical. Physics alone, it seems, cannot tell you which mathematical manifold point “really is” the point where a given material event happens. Einstein’s own solution, appealing to Leibniz’s identity of indiscernibles, was a relational view: space-time points just are the material events and field configurations located there, with no further “container” behind them. On this relational view, Wigner’s diamond is no counterexample at all. The friend’s measurement outcome is not merely located at a space-time point; it helps define that point. Since the outcome does not exist as a material event for Wigner, the point it defines cannot exist for Wigner either — and Cavalcanti’s bubble returns.

Two rival metaphysics, however, open an escape hatch. Space-time substantivalism holds that manifold points exist independently of any matter or fields, so the same bare point can be shared by Wigner and his friend even if the events occurring there differ. More intriguingly, Pienaar discusses a recent “reference frame” view, developed in work on quantum reference frames and the so-called quantum hole argument, according to which space-time points are identified with the material events of a specially chosen subset of matter — the reference frame. Once the rods-and-clocks grid inside the chamber is chosen as the reference frame, Wigner and his friend can operationally agree that the laser strike happens at one and the same juncture of the grid. This yields the almost self-contradictory-sounding sentence that there is a unique space-time point shared by both observers where the outcome occurred for the friend but not for Wigner. On either non-relational view, Cavalcanti’s original argument fails to force the bubble conclusion.

The heart of the paper is a formal proof that this escape route is illusory in general. Pienaar sets up a precise framework: each observer A has a local space-time manifold and a “localization map” that assigns each possible measurement outcome value to the space-time region where it occurs for that observer. The “classical background” condition says all observers’ manifolds can be embedded as submanifolds of one single shared manifold, with every outcome occurring in exactly one region of it. He then defines an “absolute localization map,” a global assignment of outcomes to regions in the shared manifold that is consistent with every observer’s individual localization map, and proves that such a map is a necessary condition for the shared space-time to exist. The key input is a deliberately weak assumption: that it is possible to design a quantum measurement whose different outcome values occur in disjoint space-time regions — as in a photon-polarization measurement where a beamsplitter routes the photon to one of two spatially separated detectors. Rejecting even this assumption, Pienaar notes, would imply that the locality loophole in Bell tests could not be closed in principle, since Alice’s and Bob’s outcomes could never be confined to space-like separated regions — a radically unpalatable position.

With that in place, the main theorem follows. Suppose a set of assumptions in some extended Wigner’s friend scenario forces you to reject AOE. Now redesign the experiment so that every observer’s distinct outcomes land in distinct space-time regions — always possible by the weak assumption. If a shared classical background space-time existed, the corresponding absolute localization map would, in every run, single out a consistent global assignment of one region to each outcome. But the consistency requirement means the tuple of regions uniquely determines a complete set of outcome values — a “cross-section” of all possible results. Assumptions that prohibit AOE prohibit exactly such simultaneously definite value assignments. Contradiction: no absolute localization map can exist, and the classical background condition must fail. In other words, any no-go theorem that topples AOE can be engineered into a version that also shatters the single space-time block. Rejecting absoluteness of observed events commits you, whether you like it or not, to a fragmented space-time in which some points that exist for one observer cannot belong to another’s manifold.

The implications ripple outward. Pienaar shows that perspectival interpretations cannot be “local” in any traditional sense, because notions like space-like separation and light cones presuppose the shared background that rejection of AOE undermines. QBism and RQM have both claimed some form of locality; critics have disputed those claims precisely because locality becomes ambiguous without absoluteness. Pienaar suggests the way forward is to define locality observer by observer — “relative locality,” a phrase that connects his result to an active research program in quantum gravity pioneered by Giovanni Amelino-Camelia, Laurent Freidel, Jerzy Kowalski-Glikman and Lee Smolin. Perhaps, he speculates, the fragmentation demanded by quantum foundations is not a bug but a clue to the quantum structure of space-time itself.

The result also draws a sharp line between two kinds of perspectivism. Special relativity is perspectival but benign: there exists a perspective-free description — Minkowski space-time — from which any observer’s view can be derived. Radical perspectivism, as embodied by QBism and RQM, denies that such a view-from-nowhere is even possible in principle. Philosophers like Dean Rickles had hoped that participatory-realist views could be accommodated within a block universe, much as free will can be lived from within a deterministic block. Pienaar’s theorem closes that door. Whether the resulting picture — observers each inhabiting their own space-time, overlapping only where their events are mutually defined — is a reductio ad absurdum of perspectival quantum theory, or the first sketch of a genuinely new geometry for quantum gravity, is now an open and newly urgent question.

Gamma-ray bursts are the most violent explosions in the universe since the Big Bang, unleashing in seconds the energy that the Sun radiates over billions of years. Yet despite decades of observation, fundamental questions about how these blasts launch their ultra-relativistic jets remain open. Polarization of the prompt gamma-ray emission is among the most powerful diagnostics available, because different theoretical pictures, from ordered magnetic fields to dissipative photospheres to electromagnetically dominated outflows, predict distinct polarization signatures. Measuring that polarization, however, requires accurate knowledge of where the burst is on the sky and what its spectrum looks like, and it is precisely this information that the BSD is designed to provide for POLAR-2’s dedicated High-energy Polarimetry Detector.

The HPD itself follows in the footsteps of the original POLAR instrument, which between 2016 and 2017 performed the first accurate, systematic polarization measurements of gamma-ray burst prompt emission using an array of plastic scintillator bars read out by multi-anode photomultiplier tubes. In POLAR-2, the polarimeter has been substantially upgraded, and a related wide-field photoelectric polarimeter, the LPD, is also being developed within the broader mission framework. But a Compton-scattering polarimeter needs an independent spectrometer companion to constrain the burst geometry and spectral shape, because systematic uncertainties in these quantities can otherwise bias the inferred polarization degree and angle. The BSD fills this role, and the new paper demonstrates quantitatively that it can do so to the accuracy required.

At the heart of the BSD lies a classical but powerful imaging concept: the coded-aperture mask. Positioned above the detector plane, the mask casts a shadow pattern onto the detector that shifts predictably with the direction of the incoming gamma rays, allowing the source position to be reconstructed from the observed shadowgram. This technique, famously employed by instruments such as INTEGRAL’s IBIS and the ECLAIRs camera on SVOM, gives the BSD a remarkably wide field of view. The half-coded field of view measures approximately 132 degrees by 125 degrees, meaning the instrument can monitor a huge swath of the sky simultaneously, an essential property for an experiment that must catch transient events that appear without warning and fade within seconds.

The detector itself is built from pixelated GAGG crystals, cerium-doped gadolinium aluminium gallium garnet, a dense, high-light-yield scintillator that has become a workhorse material for next-generation hard X-ray and gamma-ray instruments. The GAGG arrays were supplied by the 26th Research Institute of China Electronics Technology Group Corporation, formerly known as the Sichuan Institute of Piezoelectric and Acousto-optic Technologies. GAGG’s combination of high stopping power, good energy resolution and non-hygroscopic behaviour makes it well suited to the demanding environment of low Earth orbit, where the instrument must contend with a continuous barrage of trapped charged particles, cosmic rays and atmospheric albedo radiation. The custom front-end electronics, designed and developed by Merlin Kole and Nicolas De Angelis with funding from the Swiss National Science Foundation’s Ambizione program, provide the readout chain for the pixelated crystals.

Operating across an energy range of 10 to 1000 kiloelectronvolts, the BSD covers the full band in which gamma-ray burst prompt emission delivers the bulk of its photons, from the soft X-ray regime through the hard gamma-ray band. This broad coverage allows the instrument to fit spectral models, such as the familiar Band function parameterization used since the era of the Compton Gamma Ray Observatory’s BATSE experiment, and to extract the spectral parameters, peak energy and flux that feed directly into the polarization analysis of the HPD. The simulations carried out by the team, performed with the Geant4 toolkit developed at CERN, model the full detector response, including the in-orbit particle background environment that was characterized using established models of cosmic-ray-induced atmospheric neutron and cosmic-ray fluxes.

The headline result of the performance study concerns localization accuracy. Simulations indicate that the BSD can localize faint gamma-ray bursts similar to GRB 170817A, the famous short burst accompanying the binary neutron star merger detected by LIGO-Virgo in August 2017, to an accuracy of roughly 1.5 degrees. That benchmark matters enormously: GRB 170817A, despite being exceptionally close and bright in gravitational waves, was a faint gamma-ray emitter, and pinning down such events quickly is critical for enabling rapid follow-up by ground- and space-based telescopes. For POLAR-2’s polarimetry goals, a localization of about 1.5 degrees satisfies the core requirements, because the systematic uncertainty in the polarization measurement scales with how well the burst direction and spectrum are known. Precise input from the BSD therefore translates directly into tighter constraints on the polarization degree measured by the plastic scintillator array of the HPD.

Beyond its role as a supporting spectrometer, the BSD carries a scientific payload of its own. According to the simulations presented in the paper, the instrument possesses moderate intrinsic capability for gamma-ray burst polarimetry, particularly at energies of several hundred kiloelectronvolts, where Compton scattering kinematics within the GAGG crystals can encode polarization information. This means that even on its own, the BSD can contribute independent polarization measurements for the brightest events, providing a cross-check on HPD results and potentially extending the mission’s polarization science to a wider class of transients. The authors also note the instrument’s relevance to multi-messenger astronomy in the era of gravitational-wave detectors, where rapid, accurate localization of short gamma-ray bursts is a prized commodity.

The path from design to flight hardware has involved a demanding calibration campaign. The team acknowledges the European Synchrotron Radiation Facility in Grenoble for providing beam time and facilities, with particular thanks to the ID15A beamline team for their assistance during a one-week experimental campaign. Synchrotron beams offer finely tuned, monochromatic X-ray beams that allow instrument response to be mapped with high precision across the energy range, and preliminary ground-based calibration tests have now been combined with the Monte Carlo predictions to yield an overall evaluation of the instrument’s expected scientific performance. The data underpinning the study are available from the corresponding authors upon reasonable request.

With the design validated and performance benchmarks met in simulation and early testing, the Broad-band Spectrometer Detector moves POLAR-2 closer to launch-readiness on the China Space Station. If the instrument performs in orbit as the simulations predict, the mission will be positioned to deliver the next generation of high-precision polarization measurements for gamma-ray bursts, transforming a decades-old observational challenge into a precision test of the physics of relativistic jets, magnetic fields and explosive stellar death.