The effect at the heart of the study is called limb darkening, and it is exactly what the name suggests: a star appears fainter at its edge than at its center. This happens because stars are not uniformly bright disks. When we look toward the center of a stellar disk, our line of sight penetrates deeper into the star’s atmosphere, reaching hotter and brighter layers. Near the edge, the light we receive originates from shallower, cooler layers, and it must pass through more stellar material on its way out. The result is a gradual fade from a brilliant center to a dimmer rim, a signature that encodes the temperature and density structure of the star’s outer atmosphere.

Lead author Narsireddy Anugu, a staff scientist at Georgia State University’s CHARA Array, emphasized that the survey goes beyond simply sizing up stars. “We are not just measuring how large these stars are,” Anugu said. “We are measuring how their light is distributed across the stellar disk, which directly tests stellar-atmosphere models.” That distinction is crucial. Stellar radii can be measured in many ways, but the distribution of brightness across a stellar surface is a far more demanding test, one that probes the physics of how energy flows through a star’s outermost layers.

The CHARA Array itself is a marvel of engineering that makes such measurements possible. Located on Mount Wilson, the facility combines light from six telescopes positioned at different sites across the observatory grounds. By interfering the light together, CHARA achieves the resolving power of a much larger telescope, allowing astronomers to resolve details on the surfaces of stars that would otherwise remain forever blurred into points of light. Rather than capturing direct images, CHARA measures how the contrast of the interference pattern changes with the spacing between telescopes. These measurements reveal both the size of a star and how brightness changes across its surface, a technique that has made CHARA one of the most productive stellar-imaging facilities in the world.

Gail Schaefer, director of the CHARA Array, highlighted the facility’s unique capabilities in enabling the survey. “This study demonstrates the powerful capabilities of our facility,” Schaefer said. “By combining light from telescopes across the mountaintop, we can image stars with enough detail to see what their surfaces actually look like.” The observations were conducted by observing the same stars simultaneously through two different near-infrared filters, allowing the team to measure how the center-to-edge fading changes with wavelength. This dual-wavelength strategy is what transformed the survey from a collection of stellar portraits into a quantitative test of atmospheric physics.

The sample consisted of 31 bright stars in late stages of their lives, stars whose outer layers have expanded outward far beyond the dimensions they had during the prime of their existence. These evolved subgiant, giant, and supergiant stars provide especially valuable tests of stellar-atmosphere models because the large convective motions churning through their extended atmospheres can produce complex brightness profiles. Unlike the Sun, whose surface is a relatively thin layer governed by well-understood radiative processes, these puffed-up stars have atmospheres where convection dominates, making them natural laboratories for testing whether theoretical models can capture the messy reality of stellar surfaces.

The quantitative result is striking. Across the sample, the limb-darkening strength decreased by about 38 percent going from near-infrared wavelengths at 1.6 microns to 2.2 microns. The atmosphere models tested in the study predicted the same overall trend, weaker limb darkening at longer infrared wavelengths, but with a smaller decrease of only about 17 to 22 percent. In other words, the models capture the broad behavior correctly, but they do not fully reproduce how strongly the wavelength dependence plays out in real stars. The discrepancy means that current models may underestimate how stellar brightness changes across a star’s surface, a systematic error that could propagate into downstream measurements that astronomers rely on every day.

Why does this matter so much? Limb darkening directly affects measurements of stellar diameters, since the apparent edge of a star depends on how its brightness fades. More dramatically, it impacts how astronomers characterize exoplanets that transit across the surface of their parent star. When a planet crosses in front of a star, the amount of light blocked during the transit depends on the brightness across the stellar disk. A planet passing in front of the bright center blocks more light than one skimming the dimmer edge. If models misrepresent that brightness distribution, the inferred sizes and other properties of transiting exoplanets inherit the error. With thousands of confirmed exoplanets and missions like TESS and the upcoming PLATO survey depending on precise transit modeling, an empirical correction to limb darkening is a gift to the exoplanet community.

The survey also delivered a reassuring null result: the team found no evidence for surface features such as large starspots or hidden companion stars anywhere in the sample. That cleanliness matters, because spots and companions can mimic or mask limb-darkening signatures, and their absence confirms that the measured wavelength dependence is a genuine atmospheric effect rather than contamination from surface asymmetries or stellar multiplicity. The result strengthens the case that the discrepancy with models is real and systematic, not an artifact of a few unusual stars.

Looking ahead, the team plans to extend the comparison to broader wavelength coverage, stretching from visible light into the near-infrared, which will further constrain how well models reproduce the full spectral behavior of stellar atmospheres. The researchers also intend to image smaller main-sequence stars, which are far more difficult to resolve but are especially important for understanding transiting exoplanets around Sun-like stars, the primary targets in the search for habitable worlds. For now, the CHARA Array’s survey of 31 fading disks stands as a reminder that even the most familiar objects in the sky still hold surprises, and that pushing observational precision to its limits remains one of the most reliable ways to find where our theories fall short.

Subject of Research: Near-infrared limb darkening measurements of evolved stars using optical interferometry

Article Title: Georgia State Telescope array reveals how nearby stars fade from center to edge

Article References: Georgia State Telescope array reveals how nearby stars fade from center to edge. (n.d.). Original publication

Image Credits: AI Generated

DOI: Not provided

Keywords: CHARA Array, limb darkening, stellar atmospheres, optical interferometry, evolved stars, near-infrared astronomy, exoplanet transits, stellar diameters, Mount Wilson, supergiant stars, stellar convection, The Astronomical Journal

The scale of the problem is difficult to overstate. AGB stars shine with luminosities of 100 to 1000 times that of the Sun and have effective temperatures below 3000 kelvin, giving them radii of several hundred solar radii. Red supergiants are even more extreme: hotter, between 3450 and 4100 kelvin, with radii that can exceed 1000 solar radii. Despite their vast size, their atmospheres are astonishingly thin and dynamic. Convection, the process that transports heat from the interior to the visible surface by bulk motion of hot gas, is not a local phenomenon in these stars. In the Sun, convection manifests as millions of small granules covering the photosphere. In giant envelopes, the pressure scale height is so large that individual convective plumes grow to sizes comparable to the stellar radius itself. Observations with optical interferometers have now directly imaged these giant cells on stars such as Betelgeuse, Antares and the AGB star π1 Gruis, revealing a handful of enormous bright and dark patches that evolve on weekly to yearly timescales.

Understanding such convection demands a fully three-dimensional, time-dependent treatment, because the process is non-local, non-linear, and spans wildly disparate length scales. Traditional one-dimensional stellar evolution models rely on mixing length theory, a prescription invented by Böhm-Vitense in 1958 that compresses all of convection into a single tunable parameter. For the Sun and solar-like stars, calibrated mixing lengths work reasonably well. But in the coolest and most luminous stars, where the convective plumes span a large fraction of the star, turbulence becomes supersonic, and radiative opacity varies by orders of magnitude, the one-dimensional approximation breaks down dramatically. The review emphasizes that only global simulations, often called star-in-a-box models, which encompass the majority of the convective envelope plus the near circumstellar environment, can capture the physics at play.

Two numerical codes anchor this effort. The first, CO5BOLD, developed over decades and first applied to global giant-star simulations by Freytag and collaborators in 2002, solves the coupled equations of compressible hydrodynamics and non-local radiative energy transport on a three-dimensional Cartesian grid embedded in a fixed, spherically symmetric gravitational field. It uses a finite-volume approach with directional splitting and an approximate Riemann solver of Roe type, along with a short-characteristics method for radiative transfer built on Rosseland mean opacity tables that merge high-temperature OPAL data with low-temperature PHOENIX data. The code can be extended to magnetohydrodynamics and can handle rotation through a corotating frame. The second code, Athena++, developed by Stone and colleagues, takes a different route: it employs adaptive mesh refinement, task-based dynamic execution, and a Godunov-type shock-capturing scheme, solving the time-dependent, frequency-integrated radiative transfer equation directly over discrete ordinates, with radiation fully coupled to the gas energy and momentum equations. For red supergiants, Athena++ simulations use spherical polar coordinates covering 70.6 percent of the face-on hemisphere, with an accurate gravitational potential that includes the mass of the simulated envelope itself, a detail that matters because these loosely bound atmospheres are anything but negligible in mass.

What do these simulations actually show? The picture that emerges is spectacularly turbulent. The surfaces of simulated AGB and red supergiant stars are dominated by a few giant convective cells, with wide upflows surrounded by narrow, turbulent downdrafts, evolving on timescales from weeks to years. Non-stationary convection in the deep interior generates sound waves. These waves travel outward, slow down and compress as the temperature and sound speed drop, and their amplitudes grow with the steep decline in density until they steepen into shocks. Just above the visible surface, the resulting shock network is extraordinary: turbulent pressure can exceed twenty times the gas pressure, shock waves propagate with average Mach numbers of 3 to 4, and local peaks reach Mach numbers of 8 to 10. Individual shocks above the photosphere can span 200 to 250 solar radii, a width comparable to the star’s own radius, at temperatures around 2500 kelvin. These shocks intermittently levitate gas to distances where dust can condense, providing the crucial first step in the wind-acceleration mechanism that drives the massive mass-loss rates of 10 to the minus 6 to 10 to the minus 4 solar masses per year observed in AGB stars.

Remarkably, the simulations show that radial pulsations are not imposed on these models but emerge spontaneously from the convective motions themselves, a result first demonstrated for AGB stars by Freytag and colleagues in 2017 and confirmed in later work by Ahmad and coauthors in 2023. The extracted pulsation periods agree with those observed in long-period variable stars. This interplay between convection and pulsation leaves an observable fingerprint in stellar spectra. Spectral lines form over a range of atmospheric depths, and each line carries the imprint of the local velocity field through its shift, width and asymmetry. The tomographic method, originally conceived to study shock waves in pulsating Cepheids and Mira variables, exploits this by sorting spectral lines according to their formation depth and cross-correlating them with observed spectra. In Mira stars, the technique beautifully reproduces the classic Schwarzschild scenario, in which a shock front marching upward through the atmosphere produces a sequence from red-shifted absorption to line doubling and finally blue-shifted emission.

Red supergiants, it turns out, play by different rules. Applied to years of high-resolution spectra of the prototypical red supergiants μ Cephei and Betelgeuse, the tomographic method revealed no such clean shock scenario. Instead, the atmospheric layers exhibit a phase shift between velocity and temperature variations, tracing a hysteresis loop in the temperature-velocity plane. Crucially, synthetic spectra computed from three-dimensional CO5BOLD simulations show the same hysteresis behavior with timescales matching the observed photometric variability. The interpretation is elegant: the loop is the signature of convective turnover, in which hot rising material reaches the upper atmosphere, cools, and falls back, modulated by acoustic waves generated by the turbulent convective flow below. The method even provided insight into Betelgeuse’s famous Great Dimming of early 2020, when the star faded to magnitude 1.6. Tomography identified two successive shocks, in February 2018 and January 2019, that expanded part of the atmosphere and, through increased molecular opacity in the cooler gas, obscured the star, an explanation consistent with, and complementary to, the dust-veil picture from direct imaging.

The Athena++ simulations have yielded equally consequential results on the theoretical side. Goldberg and colleagues in 2022 performed the first three-dimensional calibration of mixing length theory in the red supergiant regime, recovering surprisingly high convective efficiencies with a mixing length parameter of roughly 3 to 4, and showing that radiation pressure provides about one-third of the support against gravity in the envelope. The simulations also revealed an inverse correlation between density and radial velocity in the outer layers, uncharacteristic of ordinary convection, and an atmosphere far more extended than one-dimensional models predict. In these luminous layers, where the opacity locally exceeds the Eddington opacity, the usual correlation between dense sinking material and convection flips entirely: cold, opaque, dense regions move outward. This has direct implications for how red supergiants die. When such a star explodes as a Type II-P supernova, its shock must traverse this wildly inhomogeneous envelope, whose surface undulates over radial ranges of 10 to 20 percent of the stellar radius. The shock therefore breaks out at different radii and different times around the star, smearing the ultraviolet flash of shock breakout over roughly 3 to 6 hours instead of the sub-hour duration predicted by spherical models, and reducing peak luminosities by factors of 2 to 10. This means astronomers cannot use breakout rise times to measure progenitor radii via light-travel-time effects, but they can instead use breakout durations to probe the scale of convective inhomogeneities themselves.

The convergence of three-dimensional radiation-hydrodynamics modeling, interferometric imaging, ALMA observations of molecular envelopes, spectropolarimetric monitoring of surface magnetism, and tomographic spectroscopy is transforming our understanding of the late lives of stars. Convection, once a fudge factor in one-dimensional codes, is now a quantitatively constrained, observable, three-dimensional process that shapes winds, dust formation, pulsation, magnetism and the very first light of supernovae. As future ultraviolet missions such as ULTRASAT catch hundreds of shock breakouts at high cadence, and as ever more powerful simulations push toward even greater realism, the giant convective cells churning on the surfaces of dying stars will continue to yield their secrets, one shock wave, one spectral line, and one hysteresis loop at a time.

Subject of Research: Three-dimensional radiation-hydrodynamics simulations of convection in the atmospheres of asymptotic giant branch and red supergiant stars

Article Title: Signatures of convection in the atmospheres of cool evolved stars

Article References: Signatures of convection in the atmospheres of cool evolved stars. (n.d.). https://doi.org/10.1007/s41115-024-00020-w

Image Credits: AI Generated

DOI: 10.1007/s41115-024-00020-w

Keywords: stellar convection, red supergiants, AGB stars, radiation hydrodynamics, CO5BOLD, Athena++, tomography, Betelgeuse, shock breakout, supernovae, mixing length theory, stellar mass loss