Astronomer leverages supercomputers to study black holes, galaxies

Over the last year, two research teams led by Stelios Kazantzidis, a Long-Term Fellow at the Center for Cosmology and Astro-Particle Physics (CCAPP) at The Ohio State University, have used what would average out to nearly 1,000 computing hours each day on the parallel high performance computing systems of the Ohio Supercomputer Center (OSC). To develop their detailed models and resulting simulations, Kazantzidis and his colleagues tapped OSC's flagship system, the Glenn IBM Cluster 1350, which features more than 9,600 Opteron cores and 24 terabytes of memory.

Kazantzidis and University of Zurich student Simone Callegari recently authored a paper, "Growing Massive Black Hole Pairs in Minor Mergers of Disk Galaxies," and submitted it for publication in the Astrophysical Journal. Their study involved a suite of high-resolution, smoothed-particle hydrodynamics simulations of merging disk galaxies with supermassive black holes (SMBHs). These simulations include the effects of star formation and growth of the SMBHs, as well as feedback from both processes.

"Binary SMBHs are very important, because once they form there is always the possibility that the two black holes may subsequently merge," Kazantzidis explained. "Merging SMBHs will produce the strongest signal of gravitational wave emission in the universe. Gravitational waves have not yet been directly detected, although Einstein predicted them in his Theory of General Relativity."

The astronomers found that the mass ratios of SMBH pairs in the centers of merged galaxies do not necessarily relate directly to the ratios they had to their original host galaxies, but are "a consequence of the complex interplay between accretion of matter (stars and gas) onto them and the dynamics of the merger process." As a result, one of the two SMBHs can grow in mass much faster than the other.

Enlarge

In simulations conducted at the Ohio Supercomputer Center, Ohio State astronomer Stelios Kazantzidis found that in minor mergers of disk galaxies, the growth of supermassive black holes is governed by complex factors. This screenshot shows gas density maps that correspond to the moment just before and after the orbit of a secondary supermassive black hole loops into the disk of the primary. Credit: Stelios Kazantzidis, The Ohio State University

Kazantzidis believes simulations of the formation of binary SMBHs have the potential to open a new window into astrophysical and physical phenomena that cannot be studied in other ways and might help to verify general relativity, one of the most fundamental theories of physics.

Kazantzidis and his colleagues also recently developed sophisticated computer models to simulate the formation of dwarf spheroidal galaxies, which are satellites of our own galaxy, the Milky Way. The study concluded that, in a majority of cases, disk-like dwarf galaxies – known in the field as disky dwarfs – experience significant loss of mass as they orbit inside their massive hosts, and their stellar distributions undergo a dramatic morphological, as well as dynamical, transformation: from disks to spheroidal systems.

"These galaxies are very important for astrophysics, because they are the most dark matter-dominated galaxies in the universe," Kazantzidis said. "Understanding their formation can shed light into the very nature of dark matter. Environmental processes like the interactions between dwarf galaxies and their massive hosts we've been investigating should be included as ingredients in future models of dwarf galaxy formation and evolution."

For this project, Kazantzidis, Callegari, Ewa Lokas of the Nicolaus Copernicus Astronomical Center – all of whom utilized the Glenn Cluster – and the rest of the team have submitted to the Astrophysical Journal an article titled, "On the Efficiency of the Tidal Stirring Mechanism for the Origin of Dwarf Spheroidals: Dependence on the Orbital and Structural Parameters of the Progenitor Disky Dwarfs."

Supercomputing centers such as OSC allow astronomers to create extremely sophisticated models that are not feasible on desktop systems. However, even with supercomputers, Kazantzidis and his colleagues find that simulating the multitude of elements involved in these galactic processes remains an enormous challenge.

"Our models can only follow a small subset of, say, the stars in a galaxy," he explained. "For example, a galaxy like our Milky Way contains hundreds of billions of stars, and even the most sophisticated numerical simulations to date can only simulate a tiny fraction of this number. The situation becomes increasingly more difficult in simulations that involve dark matter. This is because the dark matter particle is an elementary particle and, therefore, it is much less massive than a star. A galaxy like the Milky Way contains of the order of 1067 dark matter particles (that is, the number one followed by 67 zeros)."

The goal of Kazantzidis' team is to develop representations of galaxies that are as accurate as possible. Access to the Glenn Cluster increases the number of objects (or simulation particles) that can be depicted in the model, enhancing their ability to perform accurate and meaningful calculations.

"The powerful hardware and software available at OSC are particularly well-suited for cutting-edge astronomy research, such as that being conducted by Dr. Kazantzidis," said Ashok Krishnamurthy, interim co-executive director and director of research at OSC. "The results he and his colleagues have been able to achieve through their research projects are impressive and firmly demonstrate the Center's ability to help accelerate innovation and discovery."

Provided by Ohio Supercomputer Center