Early theoretical work at Illinois foreshadowed LIGO/Virgo announcement

Siv Schwink
10/16/2017

U of I Professor of Physics and Astronomy Stuart Shapiro
U of I Professor of Physics and Astronomy Stuart Shapiro
The historic October 16 joint announcement by the U.S.-based Laser Interferometer Gravitational-Wave Observatory (LIGO) and the Europe-based Virgo detector of the first detection of gravitational waves produced by colliding neutron stars is doubly noteworthy. It’s also the first cosmic event observed in both gravitational waves and light—some 70 ground- and space-based observatories observed the colliding neutron stars. This is arguably the biggest moment to date in “multi-messenger astronomy.”

In a press release issued by LIGO and Virgo collaborations, National Science Foundation Director France A. Córdova comments, “It is tremendously exciting to experience a rare event that transforms our understanding of the workings of the universe. This discovery realizes a long-standing goal many of us have had, that is, to simultaneously observe rare cosmic events using both traditional as well as gravitational-wave observatories. Only through NSF’s four-decade investment in gravitational-wave observatories, coupled with telescopes that observe from radio to gamma-ray wavelengths, are we able to expand our opportunities to detect new cosmic phenomena and piece together a fresh narrative of the physics of stars in their death throes.”

Well before the development of today’s innovative technologies supporting this simultaneous gravitational-wave and optical observation, early research in numerical relativity at the University of Illinois at Urbana-Champaign helped to lay the theoretical foundation for it. In fact, many features of the discovery had been predicted in the early computational simulations of Professor of Physics and Astronomy Stuart Shapiro and his group.

"This discovery of merging binary neutron stars with counterpart electromagnetic radiation is particularly rewarding to me and to my University of Illinois colleagues. It means that some of the cosmic phenomena that we have been studying and simulating on supercomputers for many years actually occurs in nature, much in the way we predicted," Shapiro comments.

When Shapiro joined the faculty at Illinois in 1996, bringing with him his then-postdoc Thomas Baumgarte, Shapiro was immediately invited to join NCSA as a research scientist by then-Director Larry Smarr.

“Physics Illinois was a hotbed of research on the properties of neutron stars, and NCSA was a leading center for numerical relativity,” remembers Shapiro.

At that time, Shapiro was already an established expert in the physics of compact objects—including black holes, neutron stars, and white dwarfs—and was a co-developer of one of the very first computer codes that could calculate the structure of a spinning neutron star according to Einstein’s theory of general relativity. With Baumgarte and several other collaborators, Shapiro was also one of the builders of the first code in general relativity that modeled binary neutron stars moving in a close circular orbit about each other.

“While binary neutron stars had been observed at that time, all were in widely separated orbits,” Shapiro recalls. “But these orbits were known to be shrinking at the exact rate predicted by general relativity, due to their emission of gravitational waves.”

“Of course, gravitational waves from widely separated stars are too weak to be detected directly on Earth,” he adds. “Only after these orbits would decay, after millions or billions of years, and the stars collide and merge, could the waves be detected by LIGO/Virgo.”

In 1999, using the computer tools they’d devised, Shapiro and Baumgarte, along with another postdoc in the group, Masaru Shibata, constructed models of a single rotating neutron star with mass greatly exceeding the maximum mass of a non-spinning neutron star. They proposed that such a star was a possible short-lived remnant of a binary neutron star merger and coined the name “hypermassive neutron star” for this transient object. Their model suggested that it would be stable over many rotation periods, but predicted it would eventually collapse to form a spinning black hole.

The team predicted this delayed collapse would be attributable to gravitational radiation loss, turbulent viscosity, and/or internal magnetic fields, which would radiate away or transport outwards angular momentum (spin), the hypermassive neutron star’s principal source of support against gravitational collapse.

But solving Einstein's equations of general relativity to simulate the full evolution of a binary neutron star, beginning with a binary in circular orbit and through its late inspiral, merger, and eventual collapse, was proving to be quite a challenge for numerical relativists. The same was true for the inspiral and merger of binary black holes.

So in 2000, Baumgarte and Shapiro, building on earlier work by Shibata and his graduate advisor Takashi Nakamura, both of Kyoto University, designed a new, robust numerical algorithm that enabled simulations of the inspiral and merger of these compact binary stars to be performed in a stable fashion according to general relativity. The scheme they designed is now referred to as the Baumgarte-Shapiro-Shibata-Nakamura (BSSN) numerical relativity scheme and is the most widely implemented approach to solving Einstein's equations on a supercomputer.

Shapiro and his group spent part of the next several years studying exactly how a hypermassive neutron star eventually collapses to a spinning black hole, finding in 2006 that the black hole typically becomes embedded in a low-mass disk of magnetized gas of neutron-star tidal debris that accretes slowly onto the black hole.

“This accretion fuels a central engine for generating outward electromagnetic power,” shares Shapiro. “But a central question remained: could the merger of a binary neutron star, or a binary black hole-neutron star, actually be the engine that generates the powerful short gamma-ray bursts that have been observed by satellites since the 1970s?”

“Such an identification had been hypothesized by theorists, but a firm indication that this might be possible in general relativity was still lacking,” he recounts.

Shapiro and his group spent the next several years simulating the merger of binary neutron stars and binary black hole-neutron stars with magnetized neutron stars and spinning black holes, becoming one of the very few groups worldwide capable of simulating such an event according to the laws of general relativistic magnetohydrodynamics.

The key breakthrough came in 2015 when Shapiro and his postdocs Vasilis Paschalidis and Milton Ruiz were able to evolve a magnetized binary black hole-neutron star from its late inspiral, through tidal disruption of the neutron star and merger, to the formation of a magnetized gaseous disk orbiting the remnant black hole.

They found that in addition to the burst of gravitational waves emitted during this event, a jet of plasma was eventually launched above the poles of the black hole, and this jet also shot a collimated beam of electromagnetic radiation outward into space.

“Such a jet, confined by the walls of a tightly wound, twisted, helical magnetic field—picture field lines twisted as in a barber-pole—above the black hole poles, was long believed to be a crucial ingredient in generating a short gamma-ray burst. The jet produced in our simulations had the expected duration and typical luminosity of observed short gamma-ray bursts. This finding was a real first and generated a lot of excitement in the field,” Shapiro comments.

Next the group simulated the merger of magnetized binary neutron stars and, in 2016, reported that those systems that formed a transient hypermassive neutron star following merger and then underwent delayed collapse to a black hole also produced a jet.

“We determined that this model represented a good candidate for the short gamma-ray burst central engine and a counterpart to the gravitational wave burst,” Shapiro states. “This too was all very exciting.”

Most recently, in 2017, the team performed simulations for more massive neutron star binaries that had a total mass above a critical threshold that forced them to undergo prompt collapse following merger.

“These systems do not have time to form a hypermassive neutron star and, as a result, do not launch a jet. Thus, they are apparently not good short gamma-ray burst candidates,” Shapiro continues. “Determining whether or not a binary neutron star can become a short gamma-ray burst may thus depend on its total mass, as well as the nuclear equation of state that provides pressure support in a neutron star and determines the critical threshold mass below which delayed collapse can occur.”

According to Shapiro, details of the discovery announced today will dictate many further questions and future simulations. Shapiro and his U of I group hope to continue to play a significant role in this very exciting and ongoing scientific investigation.

 

Stuart Shapiro's research has been funded by NASA and by the National Science Foundation.

 

Recent News

  • Research

An international team of researchers led by Paul Scherrer Institute postdoctoral researcher Niels Schröter now provide an important benchmark for how "strong" topological phenonena can be in a real material. Writing in Science, the team reports experiments in which they observed that, in the topological semimetal palladium gallium (PdGa), one of the most common classifiers of topological phenomena, the Chern number, can reach the maximum value that is allowed in any metallic crystal. That this is possible in a real material has never been shown before. Moreover, the team has established ways to control the sign of the Chern number, which might bring new opportunities for exploring, and exploiting, topological phenomena. Illinois Physics Professor Barry Bradlyn contributed to the theoretical work elucidating the team's experiments.

At the European Organization for Nuclear Research (CERN), over 200 physicists across dozens of institutions are collaborating on a project called COMPASS. This experiment (short for Common Muon and Proton Apparatus for Structure and Spectroscopy) uses CERN’s Super Proton Synchrotron to tear apart protons with a particle beam, allowing researchers to see the subatomic quarks and gluons that make up these building blocks of the universe. But particle beams aren’t the only futuretech in play – the experiments are also enabled by a heavy dose of supercomputing power.

New findings from physicists at the University of Illinois, in collaboration with researchers at The University of Tokyo and others, clarify the physics of coupling topological materials with simple, conventional superconductors.

Through a novel method they devised to fabricate bulk insulating topological insulator (TI) films on superconductor (SC) substrates, the researchers were able to more precisely test the proximity effect, or coupling when two materials contact one another, between TIs and SCs. They found that when the TI film is bulk insulating, no superconductivity is observed at the top surface, but if it is a metal, as in prior work, strong, long-range superconducting order is seen. The experimental efforts were led by physics Professor Tai-Chang Chiang and Joseph Andrew Hlevyack, postdoctoral researcher in Professor Chiang’s group, in collaboration with Professor James N. Eckstein’s group including Yang Bai, Professor Kozo Okazaki’s Lab at The U. of Tokyo, and five other institutes internationally. The findings are published in Physical Review Letters, which has been highlighted as a PRL Editors’ Suggestion.

  • Accolades

Illinois Physics Assistant Professor Barry Bradlyn has been selected for a 2020 National Science Foundation CAREER (Faculty Early Career Development) Award. This award is conferred annually in support of junior faculty who excel in the role of teacher-scholars by integrating outstanding research programs with excellent educational programs. Receipt of this award also reflects great promise for a lifetime of leadership within the recipients’ respective fields.

Bradlyn is a theoretical condensed matter physicist whose work studying the novel quantum properties inherent in topological insulators and topological semimetals has already shed new light on these extraordinary systems. Among his contributions, he developed a real-space formulation of topological band theory, allowing for the prediction of many new topological insulators and semimetals.