Showing posts with label Supercomputer Simulations. Show all posts
Showing posts with label Supercomputer Simulations. Show all posts

Thursday, August 28, 2025

Fresh twist to mystery of Jupiter's core

A high-resolution simulation of a planet colliding with Jupiter, used to study whether this process could be responsible for forming the planet's dilute core. The impact generates striking shock waves and stirs material in Jupiter's interior through turbulent mixing. However, the core material rapidly re-settles, and no dilute core is produced in the simulations. Credit: Jacob Kegerries/Thomas Sandnes

The mystery at Jupiter's heart has taken a fresh twist – as new research suggests a giant impact may not have been responsible for the formation of its core.

It had been thought that a colossal collision with an early planet containing half of Jupiter's core material could have mixed up the central region of the gas giant, enough to explain its interior today.

But a new study published in Monthly Notices of the Royal Astronomical Society suggests its make-up is actually down to how the growing planet absorbed heavy and light materials as it formed and evolved.

Unlike what scientists once expected, the core of the largest planet in our solar system doesn't have a sharp boundary but instead gradually blends into the surrounding layers of mostly hydrogen – a structure known as a dilute core.

How this dilute core formed has been a key question among scientists and astronomers ever since NASA's Juno spacecraft first revealed its existence.

Tn impacting planet collides with Jupiter's core in the simulations, triggering shock waves and turbulent mixing. Credit: Thomas Sandnes/Durham University
Licence type: Attribution (CC BY 4.0)

Using cutting-edge supercomputer simulations of planetary impacts, with a new method to improve the simulation's treatment of mixing between materials, researchers from Durham University, in collaboration with scientists from NASA, SETI, and CENSSS, University of Oslo, tested whether a massive collision could have created Jupiter's dilute core.

The simulations were run on the DiRAC COSMA supercomputer hosted at Durham University using the state-of-the-art SWIFT open-source software.

The study found that a stable dilute core structure was not produced in any of the simulations conducted, even in those involving impacts under extreme conditions.

Instead, the simulations demonstrate that the dense rock and ice core material displaced by an impact would quickly re-settle, leaving a distinct boundary with the outer layers of hydrogen and helium, rather than forming a smooth transition zone between the two regions.

Reflecting on the findings, lead author of the study Dr Thomas Sandnes, of Durham University, said: "It's fascinating to explore how a giant planet like Jupiter would respond to one of the most violent events a growing planet can experience.

"We see in our simulations that this kind of impact literally shakes the planet to its core – just not in the right way to explain the interior of Jupiter that we see today

This image from the simulations shows how the collision of the impactor with Jupiter's core produces striking patterns of fluid instabilities as materials mix. Credit: Jacob Kegerreis/Thomas Sandnes/Durham University
Licence type: Attribution (CC BY 4.0)

The core material rapidly re-settles in the simulations to form a core with a sharp boundary. Credit: Jacob Kegerreis/Thomas Sandnes/Durham University
Licence type: Attribution (CC BY 4.0)

Jupiter isn't the only planet with a dilute core, as scientists have recently found evidence that Saturn has one too.

Dr Luis Teodoro, of the University of Oslo, said: "The fact that Saturn also has a dilute core strengthens the idea that these structures are not the result of rare, extremely high-energy impacts but instead form gradually during the long process of planetary growth and evolution."

The findings of this study could also help inform scientists' understanding and interpretation of the many Jupiter- and Saturn-sized exoplanets that have been observed around distant stars. If dilute cores aren't made by rare and extreme impacts, then perhaps most or all of these planets have comparably complex interiors.

Co-author of the study Dr Jacob Kegerreis said: "Giant impacts are a key part of many planets' histories, but they can't explain everything!

"This project also accelerated another step in our development of new ways to simulate these cataclysmic events in ever greater detail, helping us to continue narrowing down how the amazing diversity of worlds we see in the Solar System and beyond came to be."




Media contacts:

Sam Tonkin
Royal Astronomical Society
Mob: +44 (0)7802 877 700

press@ras.ac.uk

Science contacts:

Dr Thomas Sandnes
Durham University

thomas.d.sandnes@durham.ac.uk



Further information

The paper ‘No dilute core produced in simulations of giant impacts on to Jupiter’ by T. D. Sandnes, V. R. Eke, J. A. Kegerreis, R. J. Massey and L. F. A. Teodoro, has been published in Monthly Notices of the Royal Astronomical Society. DOI: 10.1093/mnras/staf1105.



Notes for editors

About the Royal Astronomical Society

The Royal Astronomical Society (RAS), founded in 1820, encourages and promotes the study of astronomy, solar-system science, geophysics and closely related branches of science.

The RAS organises scientific meetings, publishes international research and review journals, recognises outstanding achievements by the award of medals and prizes, maintains an extensive library, supports education through grants and outreach activities and represents UK astronomy nationally and internationally. Its more than 4,000 members (Fellows), a third based overseas, include scientific researchers in universities, observatories and laboratories as well as historians of astronomy and others.

The RAS accepts papers for its journals based on the principle of peer review, in which fellow experts on the editorial boards accept the paper as worth considering. The Society issues press releases based on a similar principle, but the organisations and scientists concerned have overall responsibility for their content.

Keep up with the RAS on Instagram, Bluesky, LinkedIn, Facebook and YouTube.

About Durham University

Durham University is a globally outstanding centre of teaching and research based in historic Durham City in the UK.

We are a collegiate university committed to inspiring our people to do outstanding things at Durham and in the world.

We conduct research that improves lives globally and we are ranked as a world top 100 university with an international reputation in research and education (QS World University overwritten to this Rankings 2025).

We are a member of the Russell Group of leading research-intensive UK universities and we are consistently ranked as a top 10 university in national league tables (Times and Sunday Times Good University Guide, Guardian University Guide and The Complete University Guide).

For more information about Durham University visit:
www.durham.ac.uk/about/

Submitted by Sam Tonkin


Thursday, December 01, 2022

Emission lines from the simulated interstellar medium


Gas density and gas temperature of the simulated ISM at high gas surface densities of 100 M⦿/pc2 (left panels) and the corresponding emission line maps for hydrogen (middle), nitrogen and oxygen ions (right panels).© MPA


This plot shows the velocity dispersion (y-axis) as a function of the star formation rate (x-axis) for the warm-ionized medium (WIM, ~ 10.000 K, orange squares) and the cold neutral medium (CNM, ~ 300K, blue squares). The sound speed of the WIM is indicated by the grey band. The simulations follow the observed trends (Leroy et al. 2008, data points with black boxes) very well. The velocity dispersion increases with star formation rate and the warm gas motions become supersonic. © MPA



This map shows the diagnostic line-ratios (BPT diagram) with galaxies from the SDSS (grey) in the background. While the simulation overall (pink star) is in the star formation regime (left), a look at higher resolution (green-yellow dots) reveals significant emission in the AGN regime (right). The two lines (dotted, Kauffman et al. 2003 and solid, Kewley et al. 2001) are predictions separating galaxies with ionization by stellar radiation from radiation from active galactic nuclei (AGN, i.e. accreting supermassive black holes). In this simulation, however, the emission is generated by diffuse ionized gas heated by supernova shocks and not by AGN radiation. If this radiation were excluded the whole region would move in the direction of the arrow. © MPA




All stars in galaxies form in the dense gas of the interstellar medium (ISM). Ionizing radiation from newly born massive stars and supernova explosions lets the gas shine at characteristic wavelengths of certain atoms and ions. The relative strength of such line fluxes is an important observational diagnostic to reveal the internal state and composition of the ISM. However, emission by diffuse ionized gas has different flux ratios making accurate predictions difficult. Scientist at MPA and their European collaborators have used supercomputers to simulate a realistic star forming interstellar medium and to quantify the contribution of the diffuse gas. This finding allows for a more accurate interpretation of observations also at early cosmic times when these extreme conditions are more common than in the local Universe. 

The interstellar medium (ISM) is the backbone of every galaxy that actively forms stars. Gas accreted from the cosmic web or returned from previous generations of stars accumulates here. Dense regions in the ISM cool and collapse into molecular clouds, the nursery of newly forming stars. Most of these stars have masses like our Sun or lower. However, about one in a hundred is more massive than eight solar masses. These massive stars emit strong ionizing radiation and expel material from their surfaces as stellar winds. At the end of their lives, they die as supernovae, dumping enormous amounts of energy into the ISM. These processes are termed as “feedback” and change the conditions of the ISM resulting in a complex interplay between star formation and the turbulent, multiphase structure of the gas. Modern supercomputer simulations can study this in detail.

Scientists at MPA and collaborators at the Universities of Cologne and Heidelberg, the Czech Academy of Sciences, the Institute d’Astrophysique de Paris, and the École Polytechnique Fédérale de Lausanne of the SILCC (Simulating the Life Cycle of molecular Clouds) supercomputing project have used SuperMUC-NG at the Leibnitz Supercomputing Center, one of the world’s fastest supercomputers, to simulate a realistic star forming ISM with all fundamental constituents. This includes gas, stars, radiation, and dust, but also the non-thermal components magnetic fields and cosmic-rays. The complex interaction of the ISM with stellar feedback processes has been simulated with conditions prevailing in a typical spiral galaxy in the local Universe, similar to the solar neighborhood. In addition, models with higher gas density represent more extreme conditions in the local Universe, like in interacting galaxies, or galaxies at earlier cosmic epochs with higher gas fractions.

The simulations are very realistic and explain many details seen in observations. Their multi-phase gas structure, star formation properties, and kinematic properties agree very well with direct ISM observations. In regions with higher star formation rates, the supernova explosions trigger higher turbulent velocities in the cold (< 300 Kelvin) and warm (~ 10000 Kelvin) gas with a similar trend as observed. At high star formation rates, the warm ISM motions become supersonic and magneto-hydrodynamic shocks become ubiquitous.

The ISM is mainly composed of hydrogen and helium but is also enriched with traces of heavier elements such as oxygen or nitrogen. When these elements are ionized and excited by stellar radiation, they start emitting at specific wavelengths when they return to lower energy levels. Some of these transitions, like the [OIII] emission of doubly ionized oxygen, happens in extremely rarefied gas which only exist in space, typically in so-called HII regions around massive stars. It is therefore sometimes referred to as a “forbidden” line. The relative strength of this line emission is an important observational diagnostic to reveal the internal state and composition of the ISM at all cosmic epochs.

For example, the flux ratios of emission lines are used to characterize e.g. the density of the ISM, abundances of heavy elements, or the source of the ionizing radiation. However, also in regions without massive stars and their ionizing radiation, shocks can heat up the gas leading to collisional ionization of the same elements. This diffuse ionized gas (DIG) has very different flux ratios than the gas ionized by stellar radiation in HII regions. The uncertain contribution of the DIG in observations makes accurate predictions about the ISM structure difficult.

The research team has used simulations of a large region of the ISM to construct detailed maps for the expected emission in hydrogen, oxygen and nitrogen lines. The whole region has flux ratio properties similar to the ISM of star forming galaxies in the local Universe. However, a much higher resolution reveals the contribution of the DIG. In the line-ratio diagram (known as the BPT diagram) it populates a region, which is traditionally associated with the presence of an active galactic nucleus. In this case, however, the ionized gas with a temperature of ~ 20.000 Kelvin is created by supersonic turbulence in the warm ionized gas (WIM) triggered by many supernova explosions.

Simulations and data analysis were carried out on supercomputers at MPCDF and SuperMUC-NG at the Leibnitz Supercomputing Center (LRZ).

This work is supported by MPCDF and Gauss Centre for Supercomputing

More Information: The SILCC project



Authors

Thorsten Naab
Scientific Staff
tel.2295

tnaab@mpa-garching.mpg.de

Tim-Eric Rathjen

Original publications:

1.
Rahtjen, T.-E., Naab. T. et al.
SILCC VII - Gas kinematics and multiphase outflows of the simulated
ISM at high gas surface densities
submitted to MNRAS
Source

2.
Rahtjen, T.-E., Naab, T. et al.
Optical emission lines from the simulated interstellar medium in prep



Tuesday, February 16, 2021

Supercomputer Turns Back Cosmic Clock

Schematic diagram of the evolution of the Universe from the inflation (left) to the present (right). The “reconstruction method” winds back the evolution from right to left on this illustration to reproduce the primordial density fluctuations from the current galaxy distribution. (Credit: The Institute of Statistical Mathematics) Original size (1.7 MB)

Astronomers have tested a method for reconstructing the state of the early Universe by applying it to 4000 simulated universes using the ATERUI II supercomputer at the National Astronomical Observatory of Japan (NAOJ). They found that together with new observations the method can set better constraints on inflation, one of the most enigmatic events in the history of the Universe. The method can shorten the observation time required to distinguish between various inflation theories.

Just after the Universe came into existence 13.8 billion years ago, it suddenly increased more than a trillion, trillion times in size, in less than a trillionth of a trillionth of a microsecond; but no one knows how or why. This sudden “inflation,” is one of the most important mysteries in modern astronomy. Inflation should have created primordial density fluctuations which would have affected the distribution of where galaxies developed. Thus, mapping the distribution of galaxies can rule out models for inflation which don’t match the observed data.

However, processes other than inflation also impact galaxy distribution, making it difficult to derive information about inflation directly from observations of the large-scale structure of the Universe, the cosmic web comprised of countless galaxies. In particular, the gravitationally driven growth of groups of galaxies can obscure the primordial density fluctuations.

A research team led by Masato Shirasaki, an assistant professor at NAOJ and the Institute of Statistical Mathematics, thought to apply a “reconstruction method” to turn back the clock and remove the gravitational effects from the large-scale structure. They used ATERUI II, the world’s fastest supercomputer dedicated to astronomy simulations, to create 4000 simulated universes and evolve them through gravitationally driven growth. They then applied this method to see how well it reconstructed the starting state of the simulations. The team found that their method can correct for the gravitational effects and improve the constraints on primordial density fluctuations.

“We found that this method is very effective,” says Shirasaki. “Using this method, we can verify of the inflation theories with roughly one tenth the amount of data. This method can shorten the required observing time in upcoming galaxy survey missions such as SuMIRe by NAOJ’s Subaru Telescope.”

These results appeared as Masato Shirasaki et. al. “Constraining Primordial Non-Gaussianity with Post-reconstructed Galaxy Bispectrum in Redshift Space,” in Physical Review D on January 4, 2021.

 Related Link

Source:  NAOJ/Science


Monday, July 20, 2020

Biggest Cosmic Nuclear Bombs -- See the First Supernovae in Cutting-edge Supercomputer Simulations

A 2-D snapshot of a pair-instability supernovae as the explosion waves is about to break through the star's surface. The tiny disturbs represent fluid instability - in a region where different elements interact and mix. Image Credit: ASIAA/Ken Chen

A 3-D profile of a pair-instability supernovae. The blue cube shows the entire simulated space. Orange region is where nickel 56 decays. Image Credit: ASIAA/Ken Chen

A hypernova is a type of supernovae that are 100 times more energetic. Astronomers think these biggest cosmic bombs hold the key to peek into genesis moment of first supernovae birth. However, hypernovae are extremely rare in observation. Therefore, the ASIAA team led by Ke-Jung (Ken) Chen, using the NAOJ’s CfCA supercomputer, has completed high-resolution simulations to tackle this issue. By exploring deeply into the core, they discovered what hypernovae would look like after 300 days of their explosion. The highly innovative work offers an unprecedented conclusion: the effect which gas movement has on the luminosity estimation, has long been overlooked in previous theoretical models. This result boosts our understanding of hypernova formation and may prove to be instrumental in the future hypernova observations.

Right after the big bang, the only elements produced in the universe were hydrogen and helium, all the other natural elements did not come about until after the first stars born and evolved. To understand how the first stars and the elements formed in the first place necessitates the research on supernova. Nearly 50 years of supernova research has simply proven to us: it is not an easy task, many mysteries still open. Thinking that hypernova plays a key role in the breaking-through, the ASIAA Ken Chen team took a deep look into the heart of hypernova by numerical simulations. Because, despite hypernova ejects 100 times more energy than supernovae, observationally, it is in fact extremely rare. And that’s where astronomers start looking for help from good theoretical models and supercomputer simulations.

There are currently two theoretical models of how hypernovae formed, and Chen's team chose to build their simulations on the pair-instability supernovae model -- the one that is highly anticipated and relatively more robust (the other is called the core-collapse supernovae model, or, “the black hole model”). The difference between them is that one (the later) leaves a black hole, and the other (the former) doesn't even leave a black hole when it completely blows itself up. Usually when massive stars explode, they leave something behind – either a dense core called a neutron star or a black hole. But for the massive stars - the first stars in the universe - there was only hydrogen and helium, no traces of other elements yet. These very massive first stars can begin making pairs of electron-positrons in the end of their evolution, causing a runaway effect where the pressure drops in the star’s core, triggering a collapse, leading to an enormous explosion that completely disrupts the star, leaving nothing behind, not even a black hole.

"A star must be 140-260 times the mass of the Sun to die in such a manner" Chen said. Astronomers call stars that explode in this way the “pair-instability supernovae” (where the “pair” means the electron – positron pair.)

Such an explosion produces a large amount of radioactive isotope Ni 56, which according to Chen, is “the most important element in a supernova, because its decay energy accounted for most of the visible light of a supernova, and without it, many supernovae would have been too dark to observe".

The international team led by Ken Chen has used the NAOJ's CfCA supercomputer to run their high resolution hydrodynamical simulation for hypernova. Describing the code and the running “extremely challenging”, Chen explains, “larger the simulation scale, to keep the resolution high, the entire calculation will become very difficult and demand much more computational power, not to mention that the physics involved is also complicated.” To combat these, Chen said, their best advantage is their “well-craft code and a robust program structure.”

While previous simulations run for pair-instability supernovae model have only done 30 days after the explosion, Chen's team has run the simulation up to 300 days -- which allow them to study the entire decay process of Ni 56 (which has a half-life of 70 days, so the simulation had to be long enough). They are the first team who has done this. With extensive experience in simulating large scale supernovae, the team probed the relationship between the gas movement and energy radiation inside the supernova. What they found is that during the initial decay of Nickel 56, the heated gas expanded and formed thin-shell structures.

Chen said, "the temperature inside the gas shell is extremely high, from calculation we understand that there should be ~ 30% energy used in gas movement, then the remaining ~ 70% energy can likely become the supernova luminosity. Earlier models have ignored the gas dynamic effects, so the supernova luminosity results were all overestimated.”

Therefore, in the field of pair instability supernovae study, these results will certainly contribute to the further understanding of its radiation mechanism and observational characteristics.Several studies showed that the mass of first stars in the universe would be 100 to 300 solar masses, somewhat hint that the chances for first supernovae to be pair-instability supernovae could really be high. On the other hand, the first stars may be detectable by the James Webb Space Telescope (JWST) - the successor to the Hubble Space Telescope - making the observation and theoretical work of pair-instability supernovae an important subject in the near future.




Terminology explained:

Hypernova: A hypernova is a type of stellar explosion which ejects material with an unusually high kinetic energy, an order of magnitude higher than most supernovae.

Massive Star: A massive star is a star that is larger than eight solar masses during its regular main sequence lifetime.

Pair-instability supernovae model: Massive stars, between about 130 and 250 solar masses, are thought to lead to a pair-instability supernovae (PISN). In these stars, electron-positron pairs are created in the core. This leads the star to become dynamically unstable and leads to the collapse then explosion of these stars.

Core-collapse supernovae model: Core-collapse supernovae are dramatic explosions of giant stars at the end of their thermonuclear evolution giving birth to neutron stars and black holes.



More Information:


Paper:
Gas Dynamics of the Nickel-56 Decay Heating in Pair-Instability Supernovae published on July 14th. by The Astrophysical Journal

Authors are: KE-JUNG CHEN,S. E. WOOSLEY, AND DANIEL J. WHALEN



Media Contact:

Assistant Research Fellow Ken Chen
Email:
kjchen@asiaa.sinica.edu.tw
Tel: (02)2366-5457

Article written by: Ken Chen & Lauren Huang

Webpage Editor: Lauren Huang