3D computer simulation sheds new light on colliding stars

Scientists have created a 3D computer model depicting in unprecedented detail the aftermath of a collision between two neutron stars.

The study, published in the journal Monthly Notices of the Royal Astronomical Society, provides a better understanding of how some of the universe's fundamental elements form in cosmic collisions.

"The collision creates heavy elements including gold and lead," said Rodrigo Fernandez from the University of Alberta in Canada.

Fernandez worked with an international research team using supercomputers at the US National Energy Research Scientific Computing Center and data from a collision scientists detected in August of 2017 -- the first such collision ever observed.

"We also saw for the first time a gamma-ray burst from two neutron stars colliding. There's a large amount of science coming out of that discovery," he added, including helping researchers calculate the mass of the neutron stars and even confirm how fast the universe is expanding.

Neutron stars are the smallest and densest stars, packing more mass than Earth's sun into an area the size of a city.

When two of them collide, they merge in a flash of light and debris known as a kilonova, as material explodes outward.

Until now, computer simulations of the collisions haven't been sophisticated enough to account for where all that material ends up.

For example, the new 3D model shows that the accretion disk -- the collection of leftover debris that orbits the combined star -- ejects twice the amount of material and at higher speeds compared with previous 2D models.

"While our results do not fully reconcile all discrepancies, they bring the numbers closer together," Fernandez said, adding that his model provides a better understanding of how heavy elements are created and ejected into space.

By modelling the aftermath of the collision in such detail, Fernandez and the team were also able to account for another way matter is ejected from the collision: on an astrophysical jet, a narrow plume of particles and radiation shot out at nearly the speed of light as the stars collide.

The jet is also thought to be the source of the gamma-ray burst.

"It was expected that we could find jets, but this is the first time we've been able to model this in enough detail to see this effect emerge," said Fernandez.

Modelling the event in 3D was no easy task, he added. Although a neutron star collision happens in just milliseconds, the accretion disk can last for seconds.

Its formation also involves complex physics and many physical processes all happening at once, making it far harder for computers to simulate.

"Among the processes at work, the main culprit is actually the magnetic field acting on the matter," noted Fernandez.

"We know the equations that describe that process, but the only way that we can properly describe them is in 3D. So, not only do you have to run the simulation for a long time, you also have to model it in three dimensions, which is computationally very expensive.

"The simulation's technical aspects are impressive from a scientific standpoint because the interactions are so complex," he said.