Early theoretical work at Illinois foreshadowed LIGO/Virgo announcement

Siv Schwink

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

  • In the Media

A second solar farm planned in Savoy will put the University of Illinois in the lead among American universities in terms of solar energy, a top campus proponent says.

The campus is moving ahead with a 55-acre solar farm along the north side of Curtis Road, between First and Neil streets in Savoy, about a mile south of the first 21-acre farm on Windsor Road.

Physics Professor Scott Willenbrock, who recently served as a provost's fellow for sustainability, briefed the Academic Senate about the project Monday, saying it will help the campus meet its goal of generating 5 percent of its energy needs from renewable sources. That target was part of the Illinois Climate Action Plan, known as iCap.

  • Research
  • Biological Physics

A previously unappreciated interaction in the genome turns out to have possibly been one of the driving forces in the emergence of advanced life, billions of years ago.

This discovery began with a curiosity for retrotransposons, known as “jumping genes,” which are DNA sequences that copy and paste themselves within the genome, multiplying rapidly. Nearly half of the human genome is made up of retrotransposons, but bacteria hardly have them at all.

Nigel Goldenfeld, Swanlund Endowed Chair of Physics and leader of the Biocomplexity research theme at the IGB, and Thomas Kuhlman, a former physics professor at Illinois who is now at University of California, Riverside, wondered why this is.“We thought a really simple thing to try was to just take one (retrotransposon) out of my genome and put it into the bacteria just to see what would happen,” Kuhlman said. “And it turned out to be really quite interesting.”

  • Research
  • High Energy Physics
  • Particle Physics
The lead ion run is under way. On 8 November at 21:19, the four experiments at the Large Hadron Collider - ALICE, ATLAS, CMS and LHCb - recorded their first collisions of lead nuclei since 2015. For three weeks and a half, the world’s biggest accelerator will collide these nuclei, comprising 208 protons and neutrons, at an energy of 5.02 teraelectronvolts (TeV) for each colliding pair of nucleons (protons and neutrons). This will be the fourth run of this kind since the collider began operation. In 2013 and 2016, lead ions were collided with protons in the LHC.

Anne Sickles is co-convener of the ATLAS Heavy Ion Working Group, which will use these data.
  • Outreach
  • Quantum Information Science
  • Atomic, Molecular, and Optical Physics
  • Quantum Physics
  • Quantum Computing

A two-day summit in Chicago taking place November 8 and 9 has brought together leading experts in quantum information science to advance U.S. efforts in what’s been called the next technological “space race”—and to position Illinois at the forefront of that race. The inaugural Chicago Quantum Summit, hosted by the Chicago Quantum Exchange, includes high-level representation from Microsoft, IBM, Alphabet Inc.’s Google, the National Science Foundation, the U.S. Department of Energy, the U.S. Department of Defense, and the National Institute of Standards and Technology.

The University of Illinois at Urbana-Champaign recently joined the Chicago Quantum Exchange as a core member, making it one of the largest quantum information science (QIS) collaborations in the world. The exchange was formed last year as an alliance between the University of Chicago and the two Illinois-based national laboratories, Argonne and Fermilab.

Representing the U of I at the summit are physics professors Brian DeMarco, Paul Kwiat, and Dale Van Harlingen, who are key players in the planned Illinois Quantum Information Science and Technology Center (IQUIST) on the U of I campus. The U of I news bureau announced last week the university’s $15-million commitment to the new center, which will form a collaboration of physicists, engineers, and computer scientists to develop new algorithms, materials, and devices to advance QIS.