MPA Simulates Violent Birth of Neutron Stars

A team of researchers at the Max Planck Institute for Astrophysics conducted the most expensive and most elaborate supercomputer simulations so far to study the formation of neutron stars at the center of collapsing stars with unprecedented accuracy. These worldwide first three-dimensional models with a detailed treatment of all important physical effects confirm that extremely violent, hugely asymmetric sloshing and spiral motions occur when the stellar matter falls towards the center. The results of the simulations thus lend support to basic perceptions of the dynamical processes that are involved when a star explodes as supernova.

Stars with more than eight to ten times the mass of our Sun end their lives in a gigantic explosion, in which the stellar gas is expelled into the surrounding space with enormous power. Such supernovae belong to the most energetic and brightest phenomena in the universe and can outshine a whole galaxy for weeks. They are the cosmic origin of chemical elements like carbon, oxygen, silicon, and iron, of which the Earth and our bodies are made of, and which are bred in massive stars over millions of years or freshly fused in the stellar explosion.

Supernovae are also the birth places of neutron stars, those extraordinarily exotic, compact stellar remnants, in which about 1.5 times the mass of our Sun is compressed to a sphere with the diameter of Munich. This happens within fractions of a second when the stellar core implodes due to the strong gravity of its own mass. The catastrophic collapse is stopped only when the density of atomic nuclei - gargantuan 300 million tons in a sugar cube - is exceeded.


What, however, causes the disruption of the star? How can the implosion of the stellar core be reversed to an explosion? The exact processes are still a matter of intense research. According to the most widely favored scenario, neutrinos, mysterious elementary particles, play a crucial role. These neutrinos are produced and radiated in tremendous numbers at the extreme temperatures and densities in the collapsing stellar core and nascent neutron star. Like the thermal radiation of a heater they heat the gas surrounding the hot neutron star and thus could "ignite" the explosion. In this scenario the neutrinos pump energy into the stellar gas and build up pressure until a shock wave is accelerated to disrupt the star in a supernova. But does this theoretical idea really work? Is it the explanation of the still enigmatic mechanism driving the explosion?


Unfortunately (or luckily!) the processes in the center of exploding stars cannot be reproduced in the laboratory and many solar masses of intransparent stellar gas obscure our view into the deep interior of supernovae. Research is therefore strongly dependent on most sophisticated and challenging computer simulations, in which the complex mathematical equations are solved that describe the motion of the stellar gas and the physical processes that occur at the extreme conditions in the collapsing stellar core. For this task the most powerful existing supercomputers are used, but still it has been possible to conduct such calculations only with radical and crude simplifications until recently. If, for example, the crucial effects of neutrinos were included in some detailed treatment, the computer simulations could only be performed in two dimensions, which means that the star in the models was assumed to have an artificial rotational symmetry around an axis.


Thanks to support from the Rechenzentrum Garching (RZG) in developing a particularly efficient and fast computer program, access to most powerful supercomputers, and a computer time award of nearly 150 million processor hours, which is the greatest contingent so far granted by the "Partnership for Advanced Computing in Europe (PRACE)" initiative of the European Union, the team of researchers at the Max Planck Institute for Astrophysics (MPA) in Garching could now for the first time simulate the processes in collapsing stars in three dimensions and with a sophisticated description of all relevant physics.


"For this purpose we used nearly 16,000 processor cores in parallel mode, but still a single model run took about 4.5 months of continuous computing", says PhD student Florian Hanke, who performed the simulations. Only two computing centers in Europe were able to provide sufficiently powerful machines for such long periods of time, namely CURIE at Très Grand Centre de calcul (TGCC) du CEA near Paris and SuperMUC at the Leibniz-Rechenzentrum (LRZ) in Munich/Garching.


Many Terabytes of simulation data (1 Terabyte are thousand billion bytes) had to be analysed and visualized before the researchers could grasp the essence of their model runs. What they saw caused excitement as well as astonishment. The stellar gas did not only exhibit the violent bubbling and seething with the characteristic rising mushroom-like plumes driven by neutrino heating in close similarity to what can be observed in boiling water. (This process is called convection.) The scientists also found powerful, large sloshing motions, which temporarily switch over to rapid, strong rotational motions. Such a behavior had been known before and had been named "Standing Accretion Shock Instability", or SASI. This term expresses the fact that the initial sphericity of the supernova shock wave is spontaneously broken, because the shock develops large-amplitude, pulsating asymmetries by the oscillatory growth of initially small, random seed perturbations. So far, however, this had been found only in simplified and incomplete model simulations.


"My colleague Thierry Foglizzo at the Service d' Astrophysique des CEA-Saclay near Paris has obtained a detailed understanding of the growth conditions of this instability", explains Hans-Thomas Janka, the head of the research team. "He has constructed an experiment, in which a hydraulic jump in a circular water flow exhibits pulsational asymmetries in close analogy to the shock front in the collapsing matter of the supernova core." This phenomenon was named "SWASI" ("Shallow Water Analogue of Shock Instability") and allows one to demonstrate dynamical processes in the deep interior of a dying star by a relatively simple and inexpensive experimental setup of table size, of course without accounting for the important effects of neutrino heating. For this reason many astrophysicists had been sceptical that this instability indeed occurs in collapsing stars.


The Garching team could now demonstrate for the first time unambiguously that the SASI also plays an important role in the so far most realistic computer models. "It does not only govern the mass motions in the supernova core but it also imposes characteristic signatures on the neutrino and gravitational-wave emission, which will be measurable for a future Galactic supernova. Moreover, it may lead to strong asymmetries of the stellar explosion, in course of which the newly formed neutron star will receive a large kick and spin", describes team member Bernhard Müller the most significant consequences of such dynamical processes in the supernova core.


The researchers now plan to explore in more detail the measurable effects connected to the SASI and to sharpen their predictions of associated signals. Moreover, they plan to perform more and longer simulations to understand how the instability acts together with neutrino heating and enhances the efficiency of the latter. The goal is to ultimately clarify whether this conspiracy is the long-searched mechanism that triggers the supernova explosion and thus leaves behind the neutron star as compact remnant.

Publications:

Hanke F., Müller B., Wongwathanarat A., Marek A., Janka H.-Th.,"SASI Activity in Three-Dimensional Neutrino-Hydrodynamics Simulations of Supernova Cores", Astrophysical Journal 770, 66 (2013); http://arxiv.org/abs/1303.6269 

Foglizzo T., Masset F., Guilet J., Durand G., "Shallow Water Analogue of the Standing Accretion Shock Instability: Experimental Demonstration and a Two-Dimensional Model", Physical Review Letters 108, 051103 (2012); http://arxiv.org/abs/1112.3448 

Acknowledgments:

Elena Erastova and Markus Rampp (Rechenzentrum Garching) are acknowledged for the visualization of the results of the three-dimensional simulations, and Thierry Foglizzo for kindly providing the images and movies. This research project was supported by the Deutsche Forschungsgemeinschaft through Transregional Collaborative Research Center SFB/TR7 "Gravitational Wave Astronomy" and the Cluster of Excellence EXC153 "Origin and Structure of the Universe". The simulations were made possible by supercomputing resources (Tier-0) provide by PRACE on CURIE TN (GENCI@CEA, France) and SuperMUC (GCS@LRZ, Garching). The postprocessing of the simulation data was conducted on the IBM iDataPlex System hydra of the Rechenzentrum Garching.