Avalanche statistics suggest Tabby’s star is near a continuous phase transition

Siv Schwink

Tabby's Star in infrared (2MASS survey) and ultraviolet (GALEX).
Tabby's Star in infrared (2MASS survey) and ultraviolet (GALEX).
In its search for extrasolar planets, the Kepler space telescope looks for stars whose light flux periodically dims, signaling the passing of an orbiting planet in front of the star. But the timing and duration of diminished light flux episodes Kepler detected coming from KIC 846852, known as Tabby’s star, are a mystery. These dimming events vary in magnitude and don’t occur at regular intervals, making an orbiting planet an unlikely explanation. The source of these unusual dimming events is the subject of intense speculation. Suggestions from astronomers, astrophysicists, and amateur stargazers have ranged from asteroid belts to alien activity.  

Now a team of scientists at the University of Illinois at Urbana-Champaign—physics graduate student Mohammed Sheikh, working with Professors Karin Dahmen and Richard Weaver—proffer an entirely novel solution to the Tabby’s star puzzle. They suggest the luminosity variations may be intrinsic to the star itself.

Tabby’s star is in most regards a standard F-class star, located in the constellation Cygnus, approximately 1,276 light years from Earth. Its unusual light curve—the graph of its light intensity as a function of time—shows intense dimming events of up to 20 percent, punctuated with smaller irregular dimming events.

The location of Tabby's Star
The location of Tabby's Star
Weaver comments, “There are a few telltale signs of occultation, or dimming by an independent body blocking the view. The most important is periodicity. In Tabby's star, the small and big events are not periodic—they don’t occur at regular intervals—and this is one of the central mysteries of the light curve. “

The Ilinois team applied a statistical analysis to the light curve’s smaller irregular variations. What they found is a mathematical pattern consistent with a well-established avalanche model: the smaller dimming events are the “crackling noise” or small avalanches that are observed during the time intervals between the larger avalanches, equated to the larger dimming events. The small dimming events come in a remarkably broad range of sizes, which are distributed according to a simple scaling law. These results suggest the dimming events may be intrinsic to Tabby’s star and that the star may be near the critical point of an underlying continuous phase transition.

Sheikh performed the calculations for the analysis of the observational data. He explains the mathematical method, which starts by establishing a median dimming threshold across the light curve.

 “The threshold is an artifice we resort to in order to define what an avalanche is in the context of the light curve. In fact, the statistics are pretty robust to where we choose the threshold, so the exact value isn't important. What is important is that we get enough avalanches to do statistics. Once the light curve drops below the threshold, we consider such an event the start of an avalanche. While the light curve remains below the threshold the avalanche continues, and it stops when it increases again to a value above the threshold.”

Avalanches have two main properties, size and duration. The size is the total area enclosed by the light curve (underneath the threshold) and the threshold.

The avalanche size,” continues Sheikh, “is related to the net decrease in energy emitted by the star during the dimming event, when compared to a constant emission rate of the star, or the constant threshold value. The avalanche duration is the length of the event. We also look at the power spectral density, which is related to how much power per unit of frequency is contained in the light curve.

“Basically, we are looking at the statistical distributions of the fluctuations. All of these things have power laws associated with them. This gives us an independent way to interpret the events and check consistency with the model.”

Power laws have the interesting feature that they look the same on different scales. So when you zoom in to small scales and short times you get the same types of statistical distributions as when you zoom out to larger scales and longer times. Power laws reflect self-similarity of the system over a wide range of length and time scales—similar to fractals—that look the same when you zoom in or when you zoom out.

Significantly, the statistics of Tabby star’s smaller dimming events are consistent with the predictions of a scaling theory. In her body of research, Dahmen has established that a small-event scaling pattern punctuated by larger events is typical of systems near a phase transition. She has seen this in the intermittent deformation dynamics of nanocrystals, the event statistics of metallic glasses, rocks, and granular materials, and in earthquakes on much larger scales spanning 12 decades in length. Similar types of avalanches are also seen in neuron firing avalanches in the brain, in magnetic systems, and in many other condensed matter systems.

“We know from other systems near nonequilibrium phase transitions that a system can have small events that show power law scaling and large events that have different dynamics,” Dahmen explains. “Examples of such transitions are magnetic systems that are slowly driven with a magnetic field, or the slow deformation of somewhat brittle materials where there is often first little crackling that gets louder and louder until there is a big snap when the material breaks.”

“The small events in our star analysis would be like the little crackles while the large events would be the analogue of the big snap,” she continues. “Our mean field model is actually able to account for both, small events and large ones. It has a built in ‘weakening’ mechanism that explains why there should be two types of avalanches.”

If the dimming events are associated with an oncoming phase transition, what would the star be transitioning toward and in what time frame? Weaver explains, “As more data is being analyzed we hope it will be possible to identify exactly what type of transition this is. We don't have a deep enough understanding to get a definitive answer, and more observations are required. We can only speculate on what such a transition would be.”

“It’s important to note—lack of periodicity alone is not enough to rule out occultation. That's part of the reason why theories such as comets or planetary debris are so popular. We can't definitely rule out these things with our findings, but we can say that the power laws we have obtained are more consistent with intrinsic variation.”

Dahmen adds, “Our work gives a framework on how to analyze the data and perhaps even to classify stars into how close or how far away the stars are from such a transition. These statistical analysis tools have been tested and successfully applied to avalanche noise in magnetic systems and plastic deformation. We are transporting these tools to astrophysics to learn more about the dynamics of stars and eventually to compare different stars.

“As a next step we believe the same type of analysis should be applied to other stars to see how universal these fluctuation statistics are among the stars that are already known. In other words we would use the statistics of the noise in the light curves in these stars to learn something about the dynamical processes that are going on inside the star.”

This research is published in the December 19 issue of Physical Review Letters.

Recent News

  • Research

In a paper published today in Physical Review Letters, a team of researchers led by Illinois Physics Professor Nicolás Yunes present two conclusions. First, using NICER’s observation datamand universal relations among various properties of neutron stars, the authors infer the moment of inertia; the tidal Love number; the quadrupole moment; and the surface eccentricity of neutron star PSR J0030+045. Next, they use their inferences to propose and implement a novel test of GR.

  • Outreach

Illinois Physics Professor Smitha Vishveshwara has been elected to the Executive Committee of the American Physical Society (APS) Division of Condensed Matter Physics (DCMP). By special election, Vishveshwara fills the seat of Richard L. Greene of University of Maryland, who stepped down from the position of chair-elect.

The DCMP chair line typically represents a four-year commitment of service, and each newly elected member generally starts as vice chair for the first year, then serves as chair-elect, chair, and past chair. Vishveshwara will serve as chair-elect starting immediately, and will take over as chair in March 2022.

A new textbook from Cambridge University Press entitled Numerical Relativity: Starting from Scratch, coauthored by Bowdoin College Physics Professor Thomas W. Baumgarte and Illinois Physics and Astronomy Professor Stuart L. Shapiro, explicates this esoteric subfield of physics for today’s students and scientists. The textbook makes heavy use of analogies from Newtonian gravity, scalar fields, and electromagnetic fields. In this way, it introduces key concepts of numerical relativity in a context familiar to readers without prior expertise in general relativity. Readers can explore the concepts presented by working through textbook exercises, and can see them first-hand by experimenting with the accompanying Python sample codes.

  • Accolades

Three University of Illinois Urbana-Champaign professors have been elected to the National Academy of Sciences, one of the highest professional honors a scientist can receive.

Physics professor Nadya Mason and chemistry professors Ralph Nuzzo and Wilfred van der Donk are among 120 newly elected U.S. members – 59 of whom are women, the most elected in a single year – and 30 international members in recognition of their distinguished and continuing achievements in original research.