New simulation also shows gamma ray bursts are 10 times rarer than previously thought
EVANSTON, Ill. — A Northwestern University-led team of astrophysicists has developed the first-ever full 3D simulation of an entire evolution of a jet formed by a collapsing star, or a “collapsar.”
Because these jets generate gamma ray bursts (GRBs) — the most energetic and luminous events in the universe since the Big Bang — the simulations have shed light on these peculiar, intense bursts of light. Their new findings include an explanation for the longstanding question of why GRBs are mysteriously punctuated by quiet moments — blinking between powerful emissions and an eerily quiet stillness. The new simulation also shows that GRBs are even rarer than previously thought.
The new study will be published on June 29 in Astrophysical Journal Letters. It marks the first full 3D simulation of the entire evolution of a jet — from its birth near the black hole to its emission after escaping from the collapsing star. The new model also is the highest-ever resolution simulation of a large-scale jet.
“These jets are the most powerful events in the universe,” said Northwestern’s Ore Gottlieb, who led the study. “Previous studies have tried to understand how they work, but those studies were limited by computational power and had to include many assumptions. We were able to model the entire evolution of the jet from the very beginning — from its birth by a black hole — without assuming anything about the jet’s structure. We followed the jet from the black hole all the way to the emission site and found processes that have been overlooked in previous studies.”
Gottlieb is a Rothschild Fellow in Northwestern’s Center for Interdisciplinary Exploration and Research in Astrophysics (CIERA). He coauthored the paper with CIERA member Sasha Tchekhovskoy, an assistant professor of physics and astronomy at Northwestern’s Weinberg College of Arts and Sciences.
The most luminous phenomenon in the universe, GRBs emerge when the core of a massive star collapses under its own gravity to form a black hole. As gas falls into the rotating black hole, it energizes — launching a jet into the collapsing star. The jet punches the star until finally escaping from it, accelerating at speeds close to the speed of light. After breaking free from the star, the jet generates a bright GRB.
“The jet generates a GRB when it reaches about 30 times the size of the star — or a million times the size of the black hole,” Gottlieb said. “In other words, if the black hole is the size of a beach ball, the jet needs to expand over the entire size of France before it can produce a GRB.”
Due to the enormity of this scale, previous simulations have been unable to model the full evolution of the jet’s birth and subsequent journey. Using assumptions, all previous studies found that the jet propagates along one axis and never deviates from that axis.
But Gottlieb’s simulation showed something very different. As the star collapses into a black hole, material from that star falls onto the disk of magnetized gas that swirls around the black hole. The falling material causes the disk to tilt, which, in turn, tilts the jet. As the jet struggles to realign with its original trajectory, it wobbles inside the collapsar.
This wobbling provides a new explanation for why GRBs blink. During the quiet moments, the jet doesn’t stop — its emission beams away from Earth, so telescopes simply cannot observe it.
“Emission from GRBs is always irregular,” Gottlieb said. “We see spikes in emission and then a quiescent time that lasts for a few seconds or more. The entire duration of a GRB is about one minute, so these quiescent times are a non-negligible fraction of the total duration. Previous models were not able to explain where these quiescent times were coming from. This wobbling naturally gives an explanation to that phenomenon. We observe the jet when its pointing at us. But when the jet wobbles to point away from us, we cannot see its emission. This is part of Einstein’s theory of relativity.”
These wobbly jets also provide new insights into the rate and nature of GRBs. Although previous studies estimated that about 1% of collapsars produce GRBs, Gottlieb believes that GRBs are actually much rarer.
If the jet were constrained to moving along one axis, then it would only cover a thin slice of the sky — limiting the likelihood of observing it. But the wobbly nature of the jet means that astrophysicists can observe GRBs at different orientations, increasing the likelihood of spotting them. According to Gottlieb’s calculations, GRBs are 10 times more observable than previously thought, which means that astrophysicists are missing 10 times fewer GRBs than previously thought.
“The idea is that we observe GRBs on the sky in a certain rate, and we want to learn about the true rate of GRBs in the universe,” Gottlieb explained. “The observed and true rates are different because we can only see the GRBs that are pointing at us. That means we need to assume something about the angle that these jets cover on the sky, in order to infer the true rate of GRBs. That is, what fraction of GRBs we are missing. Wobbling increases the number of detectable GRBs, so the correction from the observed to true rate is smaller. If we miss fewer GRBs, then there are fewer GRBs overall in the sky.”
If this is true, Gottlieb posits, then most of the jets either fail to be launched at all or never succeed in escaping from the collapsar to produce a GRB. Instead, they remain buried inside.
The new simulations also revealed that some of the magnetic energy in the jets partially converts to thermal energy. This suggests that the jet has a hybrid composition of magnetic and thermal energies, which produce the GRB. In a major step forward in understanding the mechanisms that power GRBs, this is the first time researchers have inferred the jet composition of GRBs at the time of emission.
“Studying jets enables us to ‘see’ what happens deep inside the star as it collapses,” Gottlieb said. “Otherwise, it’s difficult to learn what happens in a collapsed star because light cannot escape from the stellar interior. But we can learn from the jet emission — the history of the jet and the information that it carries from the systems that launch them.”
The major advance of the new simulation partially lies in its computational power. Using the code “H-AMR” on supercomputers at the Oak Ridge Leadership Computing Facility in Oak Ridge, Tennessee, the researchers developed the new simulation, which uses graphical processing units (GPUs) instead of central processing units (CPUs). Extremely efficient at manipulating computer graphics and image processing, GPUs accelerate the creation of images on a display.
The study, “Black hole to photosphere: 3D GRMHD simulations of collapsars reveal wobbling and hybrid composition jets,” was supported by the Fermi Cycle 14 Guest Investigator program (award numbers 80NSSC22K0031, 80NSSC21K1951 and 80NSSC21K1938), the National Science Foundation (award numbers AST-2107839, AST-1815304, AST-1911080, AST-2031997, AST-2107806 and 2020747), NASA (award number 80NSSC18K0565) and the Israel Science Foundation (award number 1657/18).
<iframe width="560" height="315" src="https://www.youtube.com/embed/mbOjHB3wHY4" frameborder="0" allow="accelerometer; autoplay; encrypted-media; gyroscope; picture-in-picture" allowfullscreen></iframe>
Please credit all images to Ore Gottlieb/Northwestern University