The motivation for the study stems from a growing realization among planetary scientists that many of the smaller moons in the outer Solar System may not be the originals. Several lines of evidence have led researchers to question whether the moons we observe today around Saturn and Uranus, particularly those with radii below 1,000 kilometers, are primordial bodies that formed alongside their parent planets. Instead, a compelling hypothesis holds that most or all of these smaller moons were blown apart in large collisions at some point in their history, with the moons now visible to spacecraft reassembling later from the orbiting debris of their shattered predecessors. For Saturn’s system, some models have proposed disruptive collisions occurring remarkably recently — even within the last 100 million years. That possibility raises an immediate and consequential question: if the candidate ocean moons we see today are, in effect, second-generation bodies stitched together from collisional wreckage, can they still host oceans at all?

To tackle this problem, the research team assembled a sophisticated two-part modeling framework. The first component was a smoothed particle hydrodynamics, or SPH, model, a computational technique widely used in astrophysics to simulate the behavior of fluids and deformable materials under extreme conditions. SPH allowed the team to render catastrophic collisions in fine detail, tracking how a moon’s icy and rocky components scatter, mix, and reaccumulate when struck by an impactor of comparable or substantial size. The second component was a thermal-structural evolution model, which simulated how a moon’s interior evolves over time following such an event — how heat generated by the impact redistributes, how melted material refreezes, and how the boundary between ice and rock migrates as the body cools. By coupling these tools, the scientists could compare three scenarios side by side: a moon before an impact, the same moon after a disruption, and an otherwise identical moon that never experienced a collision at all.

The team’s definition of a disruptive impact was deliberately stringent. As Rhoden explained, the threshold required that the largest surviving fragment be less than half the size of the original target moon. In other words, these were not glancing blows or crater-forming events, but true catastrophic disruptions — collisions energetic enough to break a moon apart entirely, dispersing some fraction of its material into space while the remainder gravitationally reconsolidates into what is essentially a brand-new moon built from the debris. Such events, the researchers note, are far from exotic. Disruptive impacts were common throughout the Solar System’s formative era, and they would fundamentally rework a moon’s interior structure rather than merely scarring its surface.

The intuitive expectation might be that pumping enormous quantities of energy into a frozen moon would be a recipe for ocean creation. After all, a sufficiently violent impact should melt ice, and melted ice is, by definition, liquid water. Rhoden framed the question directly: imagine a small, frozen, geologically quiet moon. If something slams into it with enough force to cause a major collision, would that impart enough energy to generate an ocean? The models delivered a decisive answer — it is actually incredibly difficult. In most simulated scenarios, a small moon that experiences a disruption either loses an existing ocean or fails to form one in the first place. The counterintuitive explanation lies in the physics of heat loss at small scales. Rufu noted that although introducing energy into the system could in principle melt ice and create oceans, that energy dissipates very quickly, which works against ocean formation rather than for it.

Rhoden offered an everyday analogy to illustrate the effect: baking a whole potato versus baking French fries. The smaller pieces heat up and cool down far faster than an intact potato. The same principle applies to a moon blown apart by an impact. Once the body is fragmented into countless pieces of ice and rock, each fragment radiates its heat into space rapidly, and any meltwater generated during the collision freezes again before it can accumulate into a lasting ocean. The very act of breaking the moon apart, which concentrates energy in the moment, undermines the long-term thermal conditions needed to sustain liquid water. The net result, the study found, is that large-scale collisions affect the thickness and longevity of existing oceans but do not contribute meaningfully to creating new ones.

Disruptive collisions do, however, leave a lasting structural imprint on the moons that survive them. In the aftermath of a catastrophic impact, the reassembling body briefly experiences widespread melting, which triggers a process known as ice-rock differentiation. A disrupted moon may begin as a jumbled mixture of rock and ice, but when the ice melts, the denser rocky material sinks toward the center while the lighter water rises. As the water refreezes, it forms a thicker coating of ice around a more substantial core. As Rufu described it, when everything consolidates after a collision, the end product is a moon with a larger core and a thicker ice shell than it had before. This differentiation matters for future exploration, because the thickness of a moon’s ice shell and the structure of its interior strongly influence how scientists interpret gravity measurements, surface features, and potential plume activity observed by spacecraft.

The fate of an existing ocean, meanwhile, depends heavily on the moon’s size. The simulations showed that when moons possessed subsurface oceans prior to impact, larger bodies — those with radii of 1,000 kilometers or more — might retain their oceans through a disruption and reassembly. Smaller moons, by contrast, tend to lose their oceans entirely. This size-dependent outcome carries significant implications for how scientists interpret the moon systems of Saturn and Uranus. If a small moon is suspected of harboring an ocean today, and if it also shows signs of having been disrupted and reassembled, the new results suggest the two histories may be difficult to reconcile — unless the disruption occurred long enough ago, or under conditions the models have not yet fully explored.

That caveat points to the remaining uncertainties in the study. Outer Solar System moons span a wide range of physical and orbital characteristics, and the parameters governing disruptive collisions — the frequency, scale, and timing of such events — are not well constrained by current observations. The authors emphasize that more work is needed to determine whether any set of conditions could promote ocean formation in smaller reassembled moons. The modeling framework itself, combining SPH collision simulations with long-term thermal evolution, offers a template for such follow-up studies, and future refinements could probe a broader range of impact scenarios, moon compositions, and orbital environments.

For now, the findings offer a measure of reassurance for scientists excited about some of the most famous ocean world candidates. Rhoden noted that it seems unlikely that moons like Saturn’s Enceladus or Dione — both thought to harbor oceans today — would have been disrupted and reassembled within the past 100 million years. Enceladus, in particular, has captivated the planetary science community since NASA’s Cassini spacecraft observed plumes of water vapor and icy particles erupting from its south polar region, providing direct evidence of a subsurface sea in contact with a rocky core. If those oceans are to be taken as long-lived and potentially habitable environments, the new simulations suggest their host moons must have avoided catastrophic disruption over geologically recent timescales — a conclusion that, if confirmed, strengthens the case that these small, icy worlds are genuine ocean worlds rather than transient products of violent collisions.

Subject of Research: The effects of disruptive impacts on the formation and longevity of subsurface oceans in icy moons of the outer Solar System

Article Title: SwRI simulations indicate large impacts may limit ocean formation in icy moons

Article References: SwRI simulations indicate large impacts may limit ocean formation in icy moons. (n.d.). Original publication

Image Credits: AI Generated

DOI: Not provided

Keywords: icy moons, subsurface oceans, disruptive impacts, Saturn, Uranus, Enceladus, smoothed particle hydrodynamics, ice-rock differentiation, Nature Astronomy, Southwest Research Institute, ocean worlds, planetary science

The study, published in The Astrophysical Journal Letters, is the first to incorporate temperature-dependent material strength into high-resolution simulations of the moon-forming giant impact. Previous simulations, including the foundational work published in 2001 by Robin Canup and Erik Asphaug, treated the colliding bodies as fluids. That assumption seemed reasonable at the time: the impact was believed to be energetic enough to melt and vaporize large portions of both Earth and Theia, and surely molten or vaporized rock cannot hold its shape. But the new research, led by Adeene Denton, a former postdoctoral researcher at the University of Arizona’s Lunar and Planetary Laboratory who is now at SwRI, demonstrates that this simplification was hiding something fundamental about how the collision actually unfolded.

“We discovered that the preexisting geology of the Mars-sized proto-moon matters,” Denton said. “When you simulate the Earth and the moon as colliding bodies with geologic properties, it changes how the moon forms out of that impact, that’s something we considered unnecessary before.” The insight grew out of Denton’s earlier work on the Pluto-Charon system, where she found that material strength played a crucial role in how the dwarf planet and its large moon came to be. She and her colleagues wondered whether the same physics might apply to Earth’s own moon, even though the moon-forming impact involved far larger bodies and far greater energies.

To test the idea, the team employed an advanced version of smoothed particle hydrodynamics, or SPH, a computational technique that represents planetary bodies as millions of interacting particles and tracks how they deform, heat, and flow during a collision. The version used in this study, developed at the University of Arizona and the University of Bern in Switzerland, includes a strength model that gives the simulated bodies the kind of resistance to deformation expected of real geologic materials, such as the rock and metal that would have composed Theia and the proto-Earth, or solid ice. In practical terms, the modeled planets were no longer infinitely deformable blobs; they had internal structure, rigidity, and a memory of their thermal state, and those properties shaped how the impact played out.

The results were striking. Hotter bodies are weaker than colder ones, and the team found that the outcome of the moon-forming impact is highly sensitive to the temperatures of the colliding bodies at the moment of contact. In some scenarios, the impact destroyed Theia and produced a massive protolunar disk of debris around Earth, which would then need to cool and coalesce into the moon over time, much as in the classic models. But when the researchers used the same parameters as the original impact modeling, including the equal temperature structures inside both bodies, something remarkable happened: within around five hours, an intact moon emerged directly from the collision. “Depending on how hot the Earth and moon are prior to the collision, the impact can destroy Theia and produce this massive disk of debris that eventually forms the moon,” Denton said. “But when we used the same parameters as original impact modeling, down to the equal temperature structures inside both bodies, within around five hours, an intact moon emerged.”

While intact-moon outcomes have appeared in prior simulations under different assumptions, this work is the first to show that material strength and temperature are central controls on which formation pathway the impact takes. The distinction is far from academic. A moon that forms intact from the collision would preserve material in a very different state than one assembled from a disk of vaporized and molten debris, with potentially different histories of heating, differentiation, and loss of volatile elements. That means the physical and chemical properties of the moon we see today, including perhaps its volatile content, could act as a fingerprint of the thermal conditions prevailing at the moment of the giant impact.

“These surprising and exciting new results imply a potential connection between the physical properties of the moon today, including perhaps its volatile content, and the thermal state of the Earth and Theia at the time of the giant impact,” said Canup, who was not involved in the new study. “This in turn might help scientists better constrain when the moon-forming event occurred.” Because protoplanets generally begin their lives hot and cool as they age, the thermal states of Earth and Theia encode information about how long each body existed before the collision. If the moon’s characteristics can be tied to those thermal states, the moon itself becomes a kind of clock, offering a new avenue for dating one of the most consequential events in the solar system’s history.

The findings do not resolve every mystery surrounding the moon. One long-standing conundrum remains the striking similarity between the materials that make up Earth and the moon, which is difficult to explain if the moon formed largely from the remains of a distinct impactor. “Because Earth and Mars formed in the same neighborhood of the solar system, they are like siblings,” Denton said. “The moon and Earth are more like fraternal twins.” One possible explanation is that Theia and the proto-Earth formed from a common reservoir of material in the same region of the solar system, while Mars, which is compositionally distinct from both Earth and the moon, formed farther out. The new strength-aware simulations may ultimately help test such scenarios by tracking more precisely which material ends up where.

“We now know that the geophysical state of Earth and Theia play a fundamental role in shaping the outcome of the collision,” said study co-author Namya Baijal, a doctoral student in Asphaug’s group at the Lunar and Planetary Laboratory. “This gives us a new way to explore the conditions of the impact and what they might reveal about the moon’s origin.” Asphaug, a co-author of both the original 2001 work and the current study, put the shift in perspective simply: because the collision was considered violent enough to melt and vaporize large portions of Earth and Theia, earlier papers assumed it was acceptable to approximate them as fluids. “Based on our new results, however, we think that it is time to reconsider that,” he said.

The broader lesson of the study extends beyond the Earth-moon system. Material strength, long dismissed as irrelevant for planetary-scale collisions, turns out to be essential for understanding impacts among smaller bodies such as asteroids, dwarf planets, and moons, and its importance for Earth’s moon suggests that thermal and structural properties deserve a place in every giant-impact model going forward. As computational techniques continue to improve, researchers can now ask not just whether a collision can build a moon, but what kind of moon it builds, and what that moon can tell us about the conditions under which it was born. For a companion world that has hung in our sky for four and a half billion years, the moon, it seems, still has plenty of secrets about its own origin story left to reveal.

Subject of Research: Giant impact simulations of moon formation incorporating material strength and temperature

Article Title: Soft and stretchy or rock solid? New modeling study sheds light on how the Earth got its moon

Article References: Soft and stretchy or rock solid? New modeling study sheds light on how the Earth got its moon. (n.d.). Original publication

Image Credits: AI Generated

DOI: Not provided

Keywords: moon formation, giant impact, Theia, smoothed particle hydrodynamics, material strength, protolunar disk, proto-Earth, planetary science, lunar origin, computational modeling, Southwest Research Institute, University of Arizona

The stakes could hardly be higher. Gravitational wave detections of merging black holes and neutron stars have made it abundantly clear that some mechanism must shrink stellar orbits by orders of magnitude before the final, feeble gravitational-wave-driven inspiral can begin. For most systems, that mechanism is the common envelope. The concept dates back to 1976, when Bohdan Paczyński proposed it to explain V 471 Tauri, a compact white dwarf and K dwarf binary whose tight orbit could only be understood if the two stars had once shared an envelope. Since then, the common envelope has become the linchpin of compact binary formation theory, invoked to explain cataclysmic variables, binary pulsars, short gamma-ray burst progenitors, and the close white dwarf pairs that will one day merge as Type Ia supernovae. Practically all massive stars are born in multiple systems, and up to seventy percent of them experience binary interaction during their lives, making this phase not an exotic curiosity but a routine chapter in stellar history.

The challenge for theorists is that the common envelope phase is fast, asymmetric, and wildly multi-scale. The review organizes the interaction into three stages: a pre-envelope phase in which unstable mass transfer begins, a dynamical inspiral in which the companion plunges into the giant’s envelope, and a post-inspiral phase in which the system either ejects the envelope or settles toward merger. The dynamical inspiral, the heart of the problem, unfolds on timescales of roughly fifty days for the envelopes of red giant and asymptotic giant branch stars. Yet the cores of those same stars, with sound speeds exceeding a thousand kilometers per second, evolve on dynamical timescales of just twenty seconds. This five-order-of-magnitude temporal gap, combined with spatial scales spanning up to eight orders of magnitude when companions range from planets to black holes, means that fully resolved three-dimensional simulations remain out of reach for the foreseeable future.

For decades, the field relied on parametric shortcuts. The energy formalism, introduced by van den Heuvel and Webbink, balances the binding energy of the envelope against the orbital energy lost during inspiral, mediated by an efficiency parameter called alpha. Population synthesis models that predict gravitational wave event rates and Type Ia supernova progenitor channels depend critically on this parameter, yet its value remains deeply uncertain. Observational calibrations based on post-common-envelope white dwarf binaries have yielded conflicting results, with some studies favoring a constant alpha and others a value that varies with system parameters. A parallel angular momentum formalism, the so-called gamma prescription, suffers from its own lack of predictive power. Compounding the problem, the very definition of envelope binding energy is ambiguous: whether recombination energy released as ionized hydrogen and helium recombine can be tapped to eject the envelope is one of the most contested questions in the field, and the answer may push the effective efficiency above unity.

One-dimensional mechanical models attempted something more physical, integrating an equation of motion for the inspiralling companion under gravity and an assumed drag force while feeding the released energy into a stellar evolution code. The drag itself traces back to classical work on gravitational focusing: Hoyle and Lyttleton’s treatment of accretion onto a moving point mass, refined by Bondi and Hoyle into the famous accretion column picture, and complemented by Bondi’s spherical limit. Dynamical friction, first analyzed by Chandrasekhar for collisionless systems and extended to gaseous media by Dokuchaev and later Ostriker, produces the backward pull that shrinks the orbit. These analytic expressions, involving the Mach number of the companion’s motion and a Coulomb logarithm whose integration limits remain arbitrary, capture the functional dependencies but not the full nonlinear reality of a turbulent, shock-laden envelope.

Three-dimensional hydrodynamic simulations, begun tentatively in the late 1980s and maturing through waves of increasingly capable codes, have transformed the picture. They revealed just how violently non-spherical the interaction is: the envelope deforms into a toroidal shape, spiral shocks issue from the core binary, and shear instabilities shred the outflow. Crucially, they showed that ejecting the envelope using orbital energy alone is extraordinarily difficult. The inspiral tends to stall while a significant fraction of the envelope remains bound, a phenomenon attributed to several conspiring effects: the expansion and dilution of gas around the cores weakens the drag, the gas can be dragged into co-rotation with the binary until the velocity contrast vanishes, and the growing gravitational attraction between the approaching cores demands ever larger forces to sustain the inspiral. Simulations that include recombination energy through realistic equations of state eject substantially more material, though whether that energy thermalizes or radiates away remains contested, making radiation transport an urgent addition to the modeling toolkit.

The numerical techniques themselves form a fascinating landscape. Smoothed particle hydrodynamics, a Lagrangian particle method, excels at conserving angular momentum and avoiding advection errors, making it naturally suited to two orbiting stars embedded in vacuum, but it struggles to resolve shocks and low-density flows without enormous particle counts. Eulerian grid-based finite-volume schemes capture shocks and instabilities with superior accuracy through Riemann solvers, but suffer advection errors on fixed grids and require adaptive mesh refinement to concentrate resolution where it matters. Moving-mesh codes such as arepo, which evolved from cosmological simulation technology, thread the needle by advecting an unstructured Voronoi mesh with the flow, combining Lagrangian flexibility with Godunov-type accuracy. The first three-dimensional magnetohydrodynamic simulations of common envelope evolution, performed with arepo, showed that magnetic fields are amplified enormously by the magnetorotational instability in the accretion flow around the companion, though they remain dynamically subdominant during the inspiral itself.

Equally delicate is the art of setting up these simulations. The core of a giant star cannot be resolved without triggering fatal timestep restrictions, so modelers replace it with a gravitating point particle and reconstruct the envelope in hydrostatic equilibrium using a modified Lane-Emden equation. Gravitational softening must be applied to prevent the point masses from generating unphysical singular forces, yet the choice of softening length measurably alters the inspiral rate and envelope unbinding. Mapping one-dimensional stellar evolution models onto three-dimensional grids introduces noise and spurious velocities that must be damped away over several dynamical timescales before the companion is released. Even the pseudo-vacuum used to fill empty grid cells can contaminate predictions of observable light curves and colors. Energy and angular momentum conservation, monitored to the ten-percent level in early work, remains a critical diagnostic, because the envelope is so loosely bound that numerical errors can masquerade as physical ejection.

What emerges from the review is a field in rapid ascent. The past decade has seen a proliferation of global simulations, wind-tunnel experiments probing drag forces at high resolution, and growing recognition that the pre-envelope and post-inspiral phases, which unfold on thermal timescales of thousands of years, must eventually be coupled to the dynamical calculations. Observational anchors abound: post-common-envelope binaries provide statistical constraints, planetary nebulae bear the morphological fingerprints of ejected envelopes, and luminous red novae such as V1309 Sco, the merger that was caught in pre-outburst survey data, offer direct glimpses of the interaction in action. The ultimate prize is a predictive connection between the parameters of a binary entering the common envelope phase and the properties of the remnant that emerges, whether a tight double white dwarf, an X-ray binary, or a future gravitational wave source. Röpke and De Marco conclude that with rapidly advancing computational power, refined numerical techniques, and a clearer grasp of the relevant physics, one of the last fundamental unsolved problems of stellar astrophysics may finally yield.

Subject of Research: Three-dimensional hydrodynamic simulations of common-envelope evolution in binary stellar systems

Article Title: Simulations of common-envelope evolution in binary stellar systems: physical models and numerical techniques

Article References: Simulations of common-envelope evolution in binary stellar systems: physical models and numerical techniques. (n.d.). https://doi.org/10.1007/s41115-023-00017-x

Image Credits: AI Generated

DOI: 10.1007/s41115-023-00017-x

Keywords: common envelope evolution, binary stars, stellar astrophysics, hydrodynamic simulations, gravitational waves, compact objects, smoothed particle hydrodynamics, adaptive mesh refinement, recombination energy, stellar mergers, computational astrophysics, binary evolution