Astrobites RSS

multi-wavelength image of Messier 51

Editor’s Note: Astrobites is a graduate-student-run organization that digests astrophysical literature for undergraduate students. As part of the partnership between the AAS and astrobites, we occasionally repost astrobites content here at AAS Nova. We hope you enjoy this post from astrobites; the original can be viewed at

Title: X-ray View of Little Red Dots: Do They Host Supermassive Black Holes?
Authors: Tonima Tasnim Ananna (তনিমা তাসনিমঅনন্যা) et al.
First Author’s Institution: Wayne State University
Status: Published in ApJL

Among the fascinating objects observed by JWST are galaxies in the early universe known as little red dots (read more about them here and here!). These galaxies are believed to be rich in dust, which contributes to their distinct red hue. The infrared spectra of little red dots from JWST reveal the presence of strong emission lines. This emission can be a potential indicator of either 1) an active black hole or 2) vigorous star formation. That leaves us with a tantalizing question: which one is it!?

If the little red dots indeed host black holes, research suggests that the masses of these black holes must be between 10 million and 1 billion times the mass of the Sun. This poses a significant challenge — how can these black holes accumulate such large masses in such a short time from their births? This question has left astronomers puzzled, highlighting the need to understand if these little red dots do, in fact, hold such massive black holes.

In astronomy, akin to life, sometimes it helps to take a step back and observe something from a different perspective. Looking at the little red dots at a different wavelength of light can help understand what exactly is happening in these galaxies. Today’s authors set out to harness the power of multiwavelength astronomy and to look at these little red dots at X-ray wavelengths to see if they emit any X-ray signals. Even though both accreting black holes and star formation activity can lead to emission in the X-ray, the difference between the strength of the X-ray emission between the two is greater in the X-ray than in the infrared. So, obtaining the strength of X-ray emission from the little red dots should give us a good idea of whether the X-rays are caused by a black hole or by stellar processes.

The authors select a bunch of little red dots that are hidden behind a big galaxy. Since X-rays can be challenging to detect, a massive galaxy in front can act as a gravitational lens and magnify any X-ray signals from the little red dots. If little red dots indeed host really massive black holes, it is not unreasonable to expect X-ray signals from them.

Lo and Behold! There Are None!

plot of stacked X-ray signals from 21 little red dot galaxies

Figure 1: All 21 little red dots that were stacked to see if they emitted any X-ray signals. There does not appear to be a spike in the X-ray counts anywhere, indicating that there is no strong X-ray emission from the little red dots. [Adapted from Ananna et al. 2024]

The authors do not detect any X-rays from the little red dots. To further improve sensitivity, in case the X-ray signal is too weak to be determined from individual galaxies, the authors also stack the galaxies. X-rays were still undetected with significant confidence (Figure 1).

Although the authors do not detect any X-rays from the little red dots, they can still use the fact that they did not detect any X-rays to obtain conservative limits on the maximum mass of the black hole, if any, that might be present in the little red dots. The authors assume accretion at the Eddington limit, a theoretical maximum rate at which a black hole can grow by accretion of mass, and obtain black hole masses less than 1.5 x 106 M (Figure 2). This is a couple of orders of magnitude lower than the black hole masses in little red dots determined by other works.

plot showing the derived upper limits on the black hole masses of the little red dot galaxies

Figure 2: The derived upper limits to the black hole masses of the little red dots, plotted against their stellar masses. The bottom-most thick dashed green line gives the average value (assuming Eddington limit) of the black hole mass. The black hole masses are not extremely over-massive unless the accretion is a fraction of the Eddington limit. [Ananna et al. 2024]

This mismatch may arise due to the method by which the black hole masses were estimated in these little red dots. Masses were determined using broad emission lines, which have been standardized based on nearby galaxies with active galactic nuclei and thus may not work for galaxies at much higher redshift. It is also possible that the X-rays were not detected because these galaxies are incredibly dusty, but the authors believe that to be very unlikely.

The Mystery of Little Red Dots Continues!

This work suggests that little red dots will not likely host extremely massive black holes. More data coming from JWST, including a detailed study of the spectra of these little red dots, can help us understand them better. This work primarily uses data from the Chandra X-ray Observatory, which has now unfortunately been defunded despite continuing to produce very significant scientific results. (Read here to learn how you can save Chandra!) The mysteries of the universe will ultimately be solved only by looking at it from every possible wavelength, and we need such multiwavelength telescopes to achieve this lofty goal!

Original astrobite edited by Janette Suherli.

About the author, Archana Aravindan:

I am a PhD candidate at the University of California, Riverside, where I study black hole activity in small galaxies. When I am not looking through some incredible telescopes, you can usually find me reading, thinking about policy, or learning a cool language!

spiral galaxy NGC 3432

Editor’s Note: Astrobites is a graduate-student-run organization that digests astrophysical literature for undergraduate students. As part of the partnership between the AAS and astrobites, we occasionally repost astrobites content here at AAS Nova. We hope you enjoy this post from astrobites; the original can be viewed at

Title: CloudFlex: A Flexible Parametric Model for the Small-Scale Structure of the Circumgalactic Medium
Authors: Cameron B. Hummels et al.
First Author’s Institution: California Institute of Technology
Status: Accepted to ApJ

If you or I spent a long time sitting in a bathtub, we’d probably shrivel up and quickly become indistinguishable from a raisin. Galaxies, on the other hand, spend their lives bathed in a reservoir of gas known as the circumgalactic medium. In fact, the presence of the circumgalactic medium is necessary for the long-term growth and survival of the galaxy. More precisely, the circumgalactic medium provides the gas that regulates a galaxy’s ability to form stars — for example, as star formation depletes gas from the interstellar medium, inflows from the circumgalactic medium replenish that supply.

Our understanding of this diffuse gas reservoir has largely been driven by observing its effect on background light (i.e., along so-called quasar sight lines) through absorption features. These measurements have indicated the presence of gas in both hot (T ~ 106 K) and cool (T < 104K) phases. However, the presence of such a stark temperature difference between the various components of this intergalactic gas would produce fluid instabilities that could quickly drive the cold component to fragment into a collection of small “cloudlets.” Self-consistently modeling these structures in large numerical simulations is challenging because they require significant dynamic range — one must be able to resolve clouds on scales smaller than a few light-years while also modeling the large scale motions of gas on galactic (thousands of light-years) scales.

Inspired by this challenge, today’s authors set out to construct a simple model to efficiently describe the distribution and properties of these cool gas cloudlets embedded in broader cloud complexes (larger clouds composed of many cloudlets) within the circumgalactic medium. Leveraging such a model to predict observed line shapes will allow us to begin to constrain the detailed substructure of this galactic bathtub.

Constructing the Model

To construct their model, today’s authors populate these so-called cloud “complexes” with a collection of randomly sampled cloudlets and observe the resulting configuration. In detail, they statistically sample cloudlet properties based on a series of predetermined distributions. For example, the total mass of their cloud is set to be a million times the mass of our Sun. Then, based on their chosen parameters for the minimum mass of an individual cloudlet and the shape of their mass probability distribution (i.e., the function that tells you how likely it is for a cloudlet with a given mass to exist), they randomly generate cloudlets until the total mass of these objects equals the pre-set total mass. They combine this with probability distributions for the distance of the cloudlets from the center of the complex and the turbulent velocity of the cloudlets to produce a mock cloud complex structure, such as the three examples shown in Figure 1.

examples of cloudlet distributions

Figure 1: Three examples of cloudlet distributions generated using the model with the minimum cloudlet mass varied from 10-3 solar masses at left to 105 solar masses at right. Because these are 3D distributions in reality, what is shown here is the cumulative distribution along one direction, with the colors representing the density of the gas along the line of sight at any given point. [Hummels et al. 2024]

Observing the Circumgalactic Medium

In this manner the structure of the cloud complexes is determined by a choice of 11 parameters that can be flexibly varied to explore the effects of cloud structure on the resulting observations. The key observational probe that has been historically used to study the circumgalactic medium is the shape of absorption features in the spectra of distant, luminous quasars. Because the shape of these features is a direct result of the interaction of background light with the circumgalactic medium gas and the circumgalactic medium is thought to be a highly disordered, turbulent medium, observed features will in principle be sensitive to where in the medium the light passes and how dense the absorbing objects are. For example, light passing through the center of the cloud in the left panel of Figure 1 will likely interact with more cloudlets than light passing near the edge and thus the observed spectra should reflect this.

simulated absorption lines for 17 cloudlets

Figure 2: Top: A predicted magnesium-II absorption profile (black) generated from the cumulative effect of absorption in the 17 cloudlets intersected along the chosen sight line. The absorption profiles are shifted left and right based on the velocity of the cloudlets. Bottom: A depiction of the chosen sight line along the z axis with the observer at the left. The velocity vectors of the intersected cloudlets are shown with colored arrows. [Hummels et al. 2024]

In Figure 2, today’s authors demonstrate how they can use this model to predict an observed line profile (read: shape) — in this case for a magnesium-II line. This chosen line intersects 17 cloudlets as it passes through the complex (i.e., travels along the z direction in the lower panel). Each of these cloudlets has a different velocity that will cause it to contribute to a slightly different absorption feature in the spectrum because of the Doppler effect (the motion of the cloudlet causes the wavelength of the light it receives to stretch or shrink so that a given absorption line — which should have a fixed wavelength location — appears at a longer or shorter wavelength in the rest frame of the quasar), yielding the left–right offsets of the dips in the upper panel of Figure 2. Each cloudlet’s contribution to the absorption is then modeled by a so-called Voigt profile, which describes how the width of a spectral line will be determined by its internal temperature and density. The cumulative contribution of these profiles then produces the observed spectral feature.

Turning the Knobs

To understand the utility of such a framework, the authors then demonstrate how observed profiles respond to variations in the chosen parameters, such as minimum cloudlet mass. To this end, they generate 10,000 random sight lines through a given cloud complex (i.e., they do the procedure shown in Figure 2 10,000 times) and analyze the properties of the absorption in aggregate (such as by counting the number of intersected absorbers). One key metric to quantify these effects is the observed equivalent width, which basically measures the strength/depth of an absorption feature — a sight line that passes through a high density of absorbers will have a larger equivalent width than one that intersects only a few. For example, from Figure 3, they show that decreasing the minimum cloudlet mass, which should increase the total number of cloudlets (as we saw in the leftmost panel of Figure 1), increases the number of intersected cloudlets, as you might expect. Because these cloudlets are much more common, there are more sight lines that will pass through a lower density column of gas, yielding a broader distribution of column densities. Despite these effects, the overall distribution of observed equivalent widths is roughly independent of minimum cloudlet mass, but the fraction of sight lines that yield a given equivalent width is significantly different between the different minimum masses.

Figure 3: Left three panels: The distribution of the number of intersected cloudlets, column density along the line of sight, and equivalent width for different choices of the minimum cloudlet mass (different colors). Rightmost panel: The fraction of sight lines yielding equivalent widths greater than the value on the horizontal axis. [Adapted from Hummels et al. 2024]

This cloud complex model can be applied in aggregate to “simulate” the distribution of clouds in the halo around the galaxy. That is, several of these complexes will populate the circumgalactic medium around a galaxy, so a given sight line will intersect a number of complexes, each of which is composed of a collection of cloudlets, as we’ve seen. Applying this procedure to such a distribution allows one to predict observed line profiles and thus equivalent width distributions, a key observable used to probe the circumgalactic medium.

Comparing the observed and predicted distributions and the sensitivity of these distributions to underlying parameters will shed light on the detailed structure of this complex medium.

Original astrobite edited by Jessie Thwaites

About the author, Sahil Hegde:

I am an astrophysics PhD student at UCLA working on using semi-analytic models to study the formation of the first stars and galaxies in the universe. I completed my undergraduate at Columbia University, and am originally from the San Francisco Bay Area. Outside of astronomy you’ll find me playing tennis, surfing (read: wiping out), and playing board games/TTRPGs!

NGC 1277 as seen by Hubble

Editor’s Note: Astrobites is a graduate-student-run organization that digests astrophysical literature for undergraduate students. As part of the partnership between the AAS and astrobites, we occasionally repost astrobites content here at AAS Nova. We hope you enjoy this post from astrobites; the original can be viewed at

Title: RUBIES: Evolved Stellar Populations with Extended Formation Histories at z ∼ 7 − 8 in Candidate Massive Galaxies Identified with JWST/NIRSpec
Authors: Bingjie Wang et al.
First Author’s Institution: Pennsylvania State University
Status: Published in ApJL

When JWST was launched, scientists were hopeful that the telescope would be able to answer hundreds of questions about the early universe. Luckily, it has been able to answer many of these questions. And even more luckily, observations of the early universe with JWST have been able to generate even more questions we don’t know the answers to!

One new conundrum scientists are running into is that it looks like some galaxies in the early universe have what we call evolved stellar populations. This means that even though these galaxies are only a couple hundred million years old (super young, right?), they already have ~100-million-year-old stars, meaning that these galaxies had to make a large amount of stars in only a short span of time. We’re not used to seeing this many stars form this fast — in fact, we expect that star formation should be slow this early in the universe. So how do we reconcile this?

The authors of today’s article look in more depth at the spectra of three such galaxies to figure out if these galaxies have really been forming stars for millions of years already. These three galaxies are located at redshifts of z = 6.7–8.4, which places them just 600 million years after the Big Bang.

plot of spectra of three galaxies

Figure 1: Balmer break at rest wavelength 3645 Angstroms (Å) observed in the galaxy spectra. The strength of the break (the relative flux difference before and after the break) is 1.71. [Wang et al. 2024]

Balmer Breaks Create Questions

First, the authors looked to confirm that evolved stellar populations do in fact exist in these galaxies. One foolproof way to do this is to look for a Balmer break in the galaxies’ spectra. The Balmer break comes from the atmospheres of older stars, indicating that the galaxy contains a population of older stars and has had little star formation for at least the past 100 million years. Excitingly, the authors confirm that there is in fact a Balmer break in the spectra of these three galaxies, as shown in Figure 1. Unfortunately, the Balmer break only tells you that there are older stars, not necessarily how old they are or how many of them there are.

Stellar Populations or Active Galactic Nuclei?

The authors then turn to understanding the finer details of the stellar population, specifically its formation history. The authors found evidence of broad lines in their emission-line spectra, signaling that there might be active galactic nucleus activity in at least two of these galaxies. Broad lines come from Doppler broadening, where photons emitted by an active galactic nucleus have a wide distribution of velocities, resulting in a broad emission line. However, the authors aren’t entirely ready to say conclusively that there are active galactic nuclei in the spectra.

In the local, low-redshift universe, we’re pretty good at understanding how these broad lines relate to active galactic nuclei because we have hundreds of spectra to test. But at high redshift, we don’t have a lot of data. We aren’t sure why supermassive black holes act the way they do early in the universe, and we don’t want to assume they act the same as the ones we observe in the local universe. So researchers have to rely on testing theoretical models against high-redshift galaxy observations.

To test how much of the light in these galaxies is coming from active galactic nuclei or evolved stellar populations, the authors test three different models: 1) The light is mostly from stars. 2) It’s from a mixture of stars and active galactic nuclei. 3) It’s mostly from active galactic nuclei. The authors create the model spectra using Prospector and compare the modeled spectra to the observed spectra to see how well they match.

Model Comparisons

The results of the modeling are shown in Figure 2. The most important thing the authors are testing here is what these different levels of contribution from the active galactic nucleus mean for the formation history of the stellar population. For instance, the model with the maximal stellar contribution means that the stellar population had to form earlier, because there need to be enough stars to create all the light seen. The minimal stellar model can form stars a bit later because the active galactic nucleus contributes more to the light seen, so you don’t need as many stars producing a lot of light. However, there is no way to verify which model is more accurate at this stage, so the authors turn to theorizing how each model could fit into the evolutionary stages of later galaxies.

comparison of modeled spectra to observed spectra

Figure 2: Left: Each model compared to the spectrum for one of the observed galaxies. Right: Predicted star formation rate as a function of time since Big Bang for each model. As seen, models with higher active galactic nucleus contribution predict later star formation. [Wang et al. 2024]

Progenitors of Massive Quiescent Galaxies?

To find out how their models fit into the evolution of galaxies, the authors then compare the different suggested star formation histories of the medium and maximal stellar contribution models with the projected histories of massive quiescent galaxies that exist at z ~ 3-5. These massive quiescent galaxies lack an explanation for how they got so massive yet have such slow star formation now. As can be seen in Figure 3, the star formation rate traced across the age of the universe for the massive galaxies (black, blue, green, and yellow distributions) matches much better in the maximal stellar contribution model. The authors posit that these galaxies could create a lot of stellar mass in a short time, then have their star formation slowed greatly to become these massive quiescent galaxies. With this explanation, the authors note that these galaxies would be a relatively small fraction of all the galaxies in the early universe.

comparison of star formation histories of massive quiescent galaxies and the RUBIES sample galaxies

Figure 3: Comparison of medium (purple) and maximal (red) stellar models with the projected formation histories of a sample of massive quiescent galaxies at z ~ 3-5. [Wang et al. 2024]

Of course, there is still the complication that there’s convincing evidence that there are active galactic nuclei in this sample of galaxies. Also, having this many stars form this early is still an uncomfortable thought for many astronomers. The jury is still out on these galaxies and how much star formation they had early on, but future studies looking into the far infrared can hopefully distinguish between light from stars and active galactic nuclei and lend more insight into what makes up these galaxies. And hopefully with JWST we can find some more answers while exploring even more questions.

Original astrobite edited by Amaya Sinha.

About the author, Caroline von Raesfeld:

I’m a second-year PhD student at Northwestern University. My research explores how we can better understand high-redshift galaxy spectra using observations and modeling. In my free time, I love to read, write, and learn about history.

Illustration of gravitational waves emanating from a pair of black holes approaching a merger

Editor’s Note: Astrobites is a graduate-student-run organization that digests astrophysical literature for undergraduate students. As part of the partnership between the AAS and astrobites, we occasionally repost astrobites content here at AAS Nova. We hope you enjoy this post from astrobites; the original can be viewed at

Title: The Metallicity Dependence and Evolutionary Times of Merging Binary Black Holes: Combined Constraints from Individual Gravitational-Wave Detections and the Stochastic Background
Authors: Kevin Turbang et al.
First Author’s Institution: Free University of Brussels and University of Antwerp
Status: Published in ApJ

The Gravitational Wave Ocean

Before the final few seconds of the gravitational wave signal from a merger washes up on LIGOVirgoKAGRA detectors on Earth, the black holes or neutron stars that caused the powerful waves were actually orbiting around each other for much, much longer — possibly for hundreds of millions of years. The gravitational waves emitted prior to the pair’s final hurrah are too quiet to be picked up via ground-based interferometry with current detectors, and they instead wash together in a background hum like waves out at sea.

The length of time between two black holes forming and their eventual merging — otherwise known as delay time — is a useful signature for telling us where the black holes came from and the lives of the stars from which they were born. Fortunately, we don’t have to sit around for a billion years watching the cosmos play at spinning black holes around each other (though it would be fun to help if you were an interstellar cat). The distribution of time delays can be calculated using the number of merging black holes per given time at a certain redshift (called the merger rate), something that we can measure with current detections of merging black holes with ground-based detectors; as long as we know the rate of star formation and how this varies with the metallicity of the forming stars.

Instead of waiting for one binary system to merge somewhere in the cosmos, we can understand the time-delay distribution by detecting many black holes that have formed and merged at different points in time. The different information that multiple detections give us allows us to paint a picture of this time-delay distribution without waiting the lifetime of a black hole binary. There’s one crux, however: with the limited sensitivity of current ground-based gravitational wave detectors, we only really resolve merging black holes out to a redshift of around 1, about 8 billion years ago, when the universe had already started slowing its pace of forming stars. The black holes merging closer than this distance make loud enough gravitational waves that they stand out from the noise, making them resolvable, while beyond this distance gravitational waves that are quieter and harder to detect are lost in detector noise, and hence are unresolved. This means we can only understand the merger rate at small redshifts with resolved detections, while minimising the information we can get about the time-delay distribution.

Heading for Choppy Waters

Fortunately, today’s authors have a trick up their sleeve. By combining the information from the current detections loud enough to individually identify and the information from the background hum of merging black holes too weak to detect, the so-called stochastic background, they can help bypass this problem. The stochastic background is made up of unresolved signals from many more inspiralling black holes — a random sea of gravitational waves — and provides vastly more information about the black hole population at higher redshifts. The authors show that by using information about the stochastic background, we can extract information about the shape of the time-delay distribution (i.e., whether most black holes take a long time to merge or whether most merge relatively quickly after forming), the minimum time a black hole binary takes to merge, and the maximum metallicity environment that can create merging black holes. However, this background has so far been too quiet to detect. Today’s authors show that predicting these parameters using only resolved signals doesn’t tell us much about the time-delay distribution, but things change significantly when including different models of the stochastic background.

Even though we haven’t detected this stochastic background, we can still make predictions about the time-delay distribution because we know that the stochastic signal needs to be quieter than we can currently hear in the detectors. Therefore, there should be fewer merging binaries making up the background, making it not loud enough to meet the detection threshold. This means we have an estimate of the energy per unit frequency of the stochastic background, without having even detected it.

The authors’ model assumes the time-delay probability distribution is a power law with slope κ. A steeply positively sloped time-delay distribution would mean that there is very little probability of any objects merging quickly after they have formed and it’s very likely that black holes take a long time to merge. A steeply negative slope would mean that it’s very likely that black holes merge quickly after forming. Using the resolved detections and this stochastic background estimate, the authors find κ is strongly constrained to be negative. Whether the slope would be positive or negative is determined by the low-redshift black hole population — a positive slope would mean that since every greater time delay is more probable than the smaller time delay before, every day there would be more and more black holes merging. Fortunately, we don’t need to explain an ever-increasing number of merging black holes, as measurements at low redshifts show that there are fewer merging black holes closer to us as far as our resolved detections reach. Meanwhile, they can’t place any constraints on the minimum delay time or the maximum sensitivity. This is a consequence of the estimated stochastic background being much weaker than what we can detect with the current observations.

Our Next Great Voyages

With the next big set of improvements to the LIGO detectors currently being implemented comes two possibilities: either this stochastic background is detected, or it remains too quiet to be picked up. The authors show that either of these cases will be more informative about the time-delay distribution than what the current observations can say. By simulating a stochastic background for each of these future cases and adding it to the current resolved detections, they can say more about the astrophysics of the time delay and metallicities of the black hole mergers in the universe.

For the case of a future detection of the stochastic background, the slope of the time-delay distribution is likely to be highly negative, while the metallicity parameter and the minimum time delay are both constrained to smaller values, as seen in Figure 1. This means that a universe with a loud stochastic background has shorter delay times with lots of black holes merging in the early universe and at low metallicities.

corner plot of the three astrophysical parameters of the time-delay distribution

Figure 1: A corner plot of the three astrophysical parameters of the time-delay distribution (from left to right: the time delay power-law slope κ, log of the maximum metallicity in solar metallicities, and log of the minimum time delay in billions of years) measured with current resolved gravitational wave observations, and a stochastic background detected with future LIGO–Virgo–KAGRA detectors. The three panels at the top of each column show the 1D probability distributions of the three parameters, where the dotted line is the current constraints on these parameters. The three panels in the bottom left represent the joint probability distributions between pairs of parameters, showing the correlations between the parameters. [Turbang et al. 2024]

On the other hand, a non-detection at higher sensitivities means that the time-delay slope is more likely to be weakly negative, with larger values of maximum metallicity and minimum time delay. In this case, there are fewer probable values of the time-delay power-law slope (we can understand this as a narrower probability distribution of κ in Figure 2, meaning a smaller uncertainty on this value), while there is a larger uncertainty on the maximum metallicity and minimum time-delay parameters. Essentially, this result is because a shallower slope of the time delay (a more even distribution of delay times across cosmic time) is the only way to get a stochastic background that is too quiet to be detected with future sensitivities. A quiet background means that there are fewer farther away black holes merging to contribute to the background, and so there can’t have been many more black holes merging at the start of the universe as there are now.

corner plot of the three astrophysical parameters of the time-delay distribution in the case of a non-detected stochastic background

Figure 2: The same as Figure 1 but considering a non-detected stochastic background at future detector sensitivity. In this case the value of κ is still negative but closer to zero, meaning a less steep distribution of black hole time delays. There is a wider range of possible maximum metallicities and minimum time delays. [Turbang et al. 2024]

So how do we tell how long it takes for two black holes in the universe to merge? No matter if we see it or not, the stochastic background will give us more information when the next sets of gravitational wave observations are gathered. With the currently ongoing fourth observing run of the LIGO–Virgo–KAGRA detectors, hopefully it won’t be long before we can sail the stochastic seas, so get to your seats for the great black hole race of the universe!

Original astrobite edited by Evan Lewis.

About the author, Storm Colloms:

Storm is a postgraduate researcher at the University of Glasgow, Scotland. They work on understanding populations of binary black holes and neutron stars from the gravitational wave signals emitted when they merge, and what that tells us about the lives and deaths of massive stars. Outwith astrophysics they spend their time taking digital and film photos and making fun doodles of their research.

emission shell surrounding the X-ray binary CI Camelopardalis

Editor’s Note: Astrobites is a graduate-student-run organization that digests astrophysical literature for undergraduate students. As part of the partnership between the AAS and astrobites, we occasionally repost astrobites content here at AAS Nova. We hope you enjoy this post from astrobites; the original can be viewed at

Title: Discovery of an [O III] Emission Shell Around the X-ray Binary CI Cam
Authors: Robert Fesen et al.
First Author’s Institution: Dartmouth College
Status: Published in RNAAS

From Be Star to X-ray Binary

Among the spectral classes of normal stars that we know and love, there is one particularly eccentric type called the Be stars. Unfortunately, no, the “e” there does not stand for eccentric. Much like the B-type stars, Be stars are incredibly hot and massive objects. However, what sets them apart is the notable presence of excess emission features in their spectral profiles, in contrast with the absorption lines typically seen in the spectra of regular B stars. Thus, the “e” stands for emission.

Our star today, CI Camelopardalis (CI Cam), was initially cataloged as a Be star, but when it was later found that the spectrum was particularly rich in forbidden emission lines, it was reclassified as a B[e]star. The brackets designate the forbidden lines that come from an improbable transition between metastable energy levels in atoms or ions, which can only be observed in low-density astrophysical environments. As an example of this notation, the doubly ionized oxygen forbidden line is written as [O III].

Artist’s illustration of an X-ray binary system

Figure 1: Artist’s illustration of an X-ray binary system that consists of a normal star and a compact object. The gravity of the compact object pulls matter off the normal star creating an accretion disk, shown in blue. [European Space Agency, NASA and Felix Mirabel (the French Atomic Energy Commission & the Institute for Astronomy and Space Physics/Conicet of Argentina)]

CI Cam underwent a dramatic X-ray outburst in 1998, which not only revealed the presence of an unseen compact companion, but also firmly established the system as an X-ray binary. An X-ray binary is a fascinating duo of a normal star and a compact object, which could be a white dwarf, a neutron star, or a black hole (see Figure 1). The two objects orbit closely around a common center of mass, emitting intense X-ray radiation. This high-energy emission is produced by the accretion process in which material from the star gets pulled onto the compact object by its immense gravitational field. As this material spirals inward, it heats up to temperatures in the millions of degrees, causing it to radiate very brightly in the X-ray part of the electromagnetic spectrum. Interestingly, the exact nature of CI Cam’s compact companion remains unclear to this day.

The Discovery

In a serendipitous finding, today’s authors uncovered a previously unknown shell of emitting gas surrounding the CI Cam system. This shell, spanning about 8 arcminutes by 12 arcminutes on the sky, is brightest in the forbidden [O III] emission line at 5007 Angstroms (within the optical or visible light regime). The discovery unfolded from the optical data of the large supernova remnant G150.3+4.5 (G150) that was observed as part of a limited wide-field imaging survey of optical emission associated with galactic remnants. This survey was led by Robert Fesen from Dartmouth College (the first author of today’s article). The emission shell was spotted adjacent to the remnant’s filaments along its northwestern edge (see Figure 2).

Zoomed out and detailed views of the newly discovered emission shell

Figure 2: Composite image showing the emission shell around CI Cam and the western portion of the galactic supernova remnant G150.3+4.5. Blue represents the [O III] emission and red represents the Hα emission. [Adapted from Fesen et al. 2024]

The presence of [O III] emission indicates the existence of extremely hot ionizing sources that can provide the high-energy photons needed to doubly ionize oxygen atoms. Typically, strong [O III] emissions are seen to have filamentary or arc-like structures tracing the interfaces between ionized and neutral gas regions. These emissions can also be enhanced by the presence of shocks, in addition to photoionization by hot sources. An emission shell featuring strong [O III] emission around CI Cam with faint emissions (see Figure 3) implies the presence of high-velocity shocks of about 90 to 150 km/s, potentially driven by CI Cam’s outflows into the surrounding interstellar medium.

Two grayscale views of the shell in O III and H alpha filters

Figure 3: Monochromatic images of the emission shell around CI Cam as seen in [O III] (top) and in Hα (bottom). [Adapted from Fesen et al. 2024]

By adopting CI Cam’s newly refined Gaia distance of 4.1 kiloparsecs (about 13,000 light-years), the physical size of the shell, assuming it is indeed the result of CI Cam’s outflows, spans approximately 9.5 x 14.3 parsecs (31 x 47 light-years). To put this into perspective, the longer side of the nebula is roughly equivalent to a hundred billion times the distance between New York City and Los Angeles via the I-40. This dimension is comparable to the dimensions of nebulae observed around Wolf–Rayet stars, which are known to generate such large structures through their powerful stellar winds and outflows.

Due to the similarity in [O III] brightness between the shell and the supernova remnant G150’s filaments, their close proximity raised the question of a potential physical association between the two objects. However, the authors found no direct evidence in the morphology linking the shell to the remnant on the imaging data. They also took into account the available distance estimates for G150 derived from radio and γ-ray observations, which suggested a much closer distance compared to CI Cam. The authors argued that if G150 were actually at CI Cam’s distance, it would have to be a much larger and more ancient supernova remnant, which seems highly unlikely given its bright optical filaments that suggest a relatively young age. The authors then concluded that the apparent proximity of the shell to G150 is likely just a chance alignment with Cl Cam.

While X-ray binaries are primarily studied in the X-ray, optical observations offer a complementary perspective, helping us to understand the details of the binary’s interactions with its surrounding medium. That the size of the emission shell around CI Cam is comparable to those created by powerful Wolf–Rayet stars provides evidence for the strong stellar winds and the transient outburst that CI Cam has unleashed over an extended period. This, in turn, promises to shed light on the mass-loss mechanisms and accretion processes in this class of objects. This discovery not only highlights the significance of wide-field imaging surveys in unveiling previously unknown structures around peculiar objects, but also opens up a new observational window for advancing our understanding of how such systems enrich and shape their cosmic neighborhoods.

Original astrobite edited by Ivey Davis.

About the author, Janette Suherli:

Janette is a PhD student at University of Manitoba in Winnipeg, Canada. Her research focuses on the utilization of integral field spectroscopy for the studies of supernova remnants and their compact objects in the optical. She is also the current chair of Graduate Student Committee for the Canadian Astronomical Society (CASCA). She grew up in Indonesia where it is summer all year round! Before pursuing her PhD in astrophysics, Janette worked as a data analyst for a big Indonesian tech company, combating credit card fraud.

photograph of the Sun with plasma streaming from its surface

Editor’s Note: Astrobites is a graduate-student-run organization that digests astrophysical literature for undergraduate students. As part of the partnership between the AAS and astrobites, we occasionally repost astrobites content here at AAS Nova. We hope you enjoy this post from astrobites; the original can be viewed at

Title: The Death of Vulcan: NEID Reveals That the Planet Candidate Orbiting HD 26965 Is Stellar Activity
Authors: Abigail Burrows et al.
First Author’s Institution: Dartmouth College and NASA’s Jet Propulsion Laboratory
Status: Published in AJ

It’s possible you may have heard of HD 26965, otherwise known as 40 Eridani A, the stellar host for the fictional planet Vulcan (homeworld of Spock, of Star Trek fame). Back in 2018, two teams of astronomers announced the likely presence of a super-Earth — Vulcan, perhaps — orbiting HD 26965, based on radial-velocity measurements from a variety of instruments. These astronomers measured the planet’s mass to be equivalent to more than 8 Earths and the orbital period to be roughly 42 days. But before Trekkies could break out the Romulan ale, researchers noted that more work was needed to fully separate the radial-velocity signals from those of the exoplaneteer’s eternal enemy: stellar activity, which can mimic planetary signals.

The Fate of Vulcan

Since then, study after study has highlighted stellar activity as the likely primary source of Vulcan’s radial-velocity signal, raising further uncertainty about the existence of the planet. What might bring us closer to a definite answer would be to combine detailed radial-velocity analysis techniques with what we already know about the effects of solar activity on radial-velocity measurements. This knowledge can be applied to high-resolution spectra of HD 26965 taken by NEID at the Kitt Peak National Observatory. NEID represents the latest and greatest in ground-based spectrometers, ensuring clearer, more frequent data than previous studies had access to.

Important to determining the status of HD 26965 b (the standard name for our Vulcan) is the concept of phase lag. Based on what we know from studying solar activity, the observational parameters used to study stellar activity can become offset in time from the effects of stellar activity on radial-velocity measurements. This would mean that stellar activity might not have been fully corrected for in early studies of this planet. If a decaying starspot or plage were present on HD 26965, correcting for it would lead to a reduction in the radial-velocity signal.

Combining Radial Velocity Analysis Techniques

To account for this phase lag, the authors first compute what are essentially smoothed radial-velocity signals by taking the sum of every radial velocity corresponding to every spectral line, weighted by their corresponding errors. These are corrected for NEID systematic quirks (such as by converting measurements into the stellar frame from the observatory frame NEID usually works in). This produces a “template” radial velocity that the authors compare each observed spectral line to, noting differences between observed radial velocity and “template” radial velocity. This comparison allows them to view how phase lag may affect the amplitude and location of the radial-velocity signal.

The Hits Just Keep On Coming

The observed NEID radial velocities have a similar 42-day period to the planet model proposed by an earlier study, but Figure 1 shows that they are out of phase by 30–40%! Straight off the bat, it’s not looking so good for Vulcan with strike number one. In addition, comparing these radial velocities to activity indicator lines Ca Ⅱ, H, and K reveals that they share a similar period but do not seem to be correlated. Strike number two!

observed radial velocities compared to the model including a planet in the system

Figure 1: The NEID-observed radial velocities of HD 26965 (marked in purple) compared to the planet model proposed by an earlier study (Ma et al. 2018, marked in yellow). The NEID radial velocities are out of phase with the planet model by 30–40%. [Burrows et al. 2024]

Going back to phase lag, the authors calculate phase offsets for all activity indicators using a Gaussian process, finding a consistent phase lag between 4.65 and 6.67 days, more than 10% of the star’s rotation period. Shifting data with this phase lag leads to significant reduction of the radial-velocity signal strength when modeled with a 42-day period, and overall over a variety of periods, seen in the periodograms in Figure 2. Strike number three!

radial-velocity measurements and periodograms per stellar activity

Figure 2: Left: A variety of corrected radial-velocity measurements per activity metric, marked with purple circles, accompanied by original radial velocities, marked with gray squares. Right: Periodograms for corrected radial velocities, marked in solid purple, and for original radial velocities, marked in gray dashes. The red diamonds denote periods with the highest remaining power. At the supposed 42-day period, the probability of that being the actual period (called the power) drops significantly for most activity indicators. [Burrows et al. 2024]

Stop, Stop, He’s Already Dead!

These are just a small sampling of demonstrations in the article — ultimately, they all seem to point towards stellar activity probably being the source of the radial-velocity signal from HD 26965, unfortunately for Vulcan.

Fortunately for us, this article aims to act as a Swiss army knife of sorts. It accomplished several things — first, the authors showed that phase lags may be important in determining the relationship between radial velocity and whatever may have affected it. Second, a bundling of analyses makes for a stronger case compared to a single analysis technique on its own, giving more solid evidence for the stellar activity hypothesis.

Last but not least, the methodologies used in this article can be used on radial velocities from other stars, including our own Sun! In addition to testing this multi-pronged analysis on other types of stars, the authors hope to use it on our Sun during especially turbulent activity periods. By refining our understanding of how our Sun and other stars fluctuate through time, we better our exoplanet detection techniques and our chances of finding habitable worlds.

Original astrobite edited by Dee Dunne.

About the author, Diana Solano-Oropeza:

I’m a first-year astronomy PhD student at Cornell University, where I study exoplanets, stars, and habitability using Gaia data. I earned my BS in physics at Drexel University before entering the Bridge to the PhD in STEM program at Columbia University. There, I researched TESS-detected exoplanets for two years. My hobbies include practicing Muay Thai, writing fiction, and playing video games. You can check out my website at

Illustration of the exoplanet WASP-12 b.

Editor’s Note: Astrobites is a graduate-student-run organization that digests astrophysical literature for undergraduate students. As part of the partnership between the AAS and astrobites, we occasionally repost astrobites content here at AAS Nova. We hope you enjoy this post from astrobites; the original can be viewed at

Title: Doomed Worlds I: No New Evidence for Orbital Decay in a Long-Term Survey of 43 Ultra-Hot Jupiters
Authors: Elisabeth R. Adams et al.
First Author’s Institution: Planetary Science Institute
Status: Accepted to PSJ

Ultra-hot Jupiters: Fleeting Beauties?

Ultra-hot Jupiters are gas giants orbiting close to their stars, with orbital periods less than roughly three days. Because these planets are large and close to their stars, they produce large signals, making them promising targets for detection and characterization. But, you know what they say: all good things must come to an end. These planets are expected to experience large tidal effects from their stars, resulting in a loss of angular momentum, orbital decay, and, eventually, the star engulfing the planet.

Several lines of evidence support the picture that ultra-hot Jupiters are subject to orbital decay over long timescales. For instance, stars hosting hot Jupiters tend to be younger than the average exoplanet host star, and ultra-hot Jupiters are rarer around older host stars. A recent research article even reports a direct detection of a planetary engulfment event from the sudden, short-lived increase in brightness of a faint star. While this evidence paints a compelling picture, it is difficult to estimate how quickly we expect ultra-hot Jupiters to experience orbital decay given theoretical uncertainty in stellar tidal effects.

Keeping Time: Working Hard or Hardly Working?

Because we expect orbital decay to occur and we know of thousands of transiting exoplanets, some of which have been observed for decades, several teams have searched for orbital decay and found two promising detections: WASP-12 b and Kepler-1658 b. Searching for orbital decay relies on the detection of transit-timing variations. This is when a planet passes in front of its star along our line of sight earlier or later than expected. There are many sources of transit-timing variations in addition to orbital decay, including precession, perturbations from companion planets or stars, or acceleration of the host star toward Earth.

Let’s say we observe the transit of a planet at time t = 0 and know its period, P. We expect to observe transits at time P, 2P, 3P, etc. In the case of orbital decay, the period of the planet is getting shorter as time goes on, meaning we need to factor in an additional quadratic term encoding the rate at which the period shrinks. Then, to detect a statistically significant signal of orbital decay, we need to show that the quadratic model fits the data better than the constant-period linear model. The authors of today’s article attempt to do exactly this but with an impressive level of care and attention to detail.

This science depends upon accurate and precise measurements of transit times for each planet in the authors’ sample, most of which have been observed by several teams with various instruments and methodologies over years or decades. Moreover, each transit time must be reported in one unified timing system (click here for more info on one of the most common timing systems). Not every transit observation properly identifies its timing system or accurately converts between timing systems, meaning any historical inaccuracies complicate such studies.

Statistical Methods: Comparing Models

The authors of today’s article compile transit times for 43 ultra-hot Jupiters and take new transit data for six of those planets to extend the temporal baseline of observations. To assess whether the linear (constant period) or quadratic (changing period) model fits the data better, the authors use the Bayesian information criterion (BIC), a model selection criterion that awards a good fit but penalizes additional parameters to avoid overfitting. The authors calculate the difference in the BIC (ΔBIC) between the linear and quadratic models, with a larger ΔBIC suggesting the quadratic model is preferred.

The authors additionally perform a variety of steps to ensure the quality of the data. They perform omit-one tests, where individual transit times are removed from the analysis and flagged if they alter the ΔBIC result by more than 25%. This step is essential since one transit time recorded inaccurately or in the wrong system could result in a spurious detection of orbital decay. The authors additionally perform a “rescaling test,” where the error bars are scaled up to account for unrealistically small error bars in reported transit times.


As shown in Figure 1, four planets out of the sample of 43 had a ΔBIC above the detection threshold, including WASP-12 b, which had been found previously to show orbital decay. The authors measure WASP-12 b’s period to be shrinking by 30 milliseconds per year, matching previously reported values. The planets WASP-121 b and WASP-46 b show tentative period increases, but these results are highly dependent on one or a few data points, warranting further observations. The planet TrES-1 b has prior tentative claims of its period decreasing, and the authors find a tentative period decrease of 18 milliseconds per year. However, this rate of period shrinkage suggests stellar tidal effects that would differ greatly from theoretical predictions, perhaps suggesting a cause of period decrease other than orbital decay.

plot of delta BIC for each of the 43 planets studied

Figure 1: The value of each planet’s ΔBIC shown relative to a threshold of ΔBIC = 30 (top), zoomed in results (middle), and rescaled results (bottom) with scaled up uncertainties, indicating only WASP-12 b definitively shows signs of orbital decay. [Adams et al. 2024]

Only one planet, WASP-12 b, was found to have a clear period decrease after rescaling error bars, as shown in the bottom panel of Figure 1. The authors predict that if the orbits of the other planets in the sample were decaying as rapidly as the orbit of WASP-12 b, they could have found significant detections of period decrease in roughly half the sample. There is thus no evidence that orbital decay is common among ultra-hot Jupiters, which is possibly confounding considering the other lines of evidence that suggest ultra-hot Jupiters are subject to decay. Though patience is required, as time goes on, it will be possible to search for orbital decay around more planets at higher precision, helping us ascertain the ultimate fate of close-in planets.

Original astrobite edited by Ivey Davis.

About the author, Kylee Carden:

I am a first-year PhD student at The Ohio State University, where I am an observer of planets outside the solar system. I’m involved with the transiting exoplanet survey of the upcoming Roman Space Telescope and working with high-resolution spectroscopic observations of exoplanet atmospheres. I am a huge fan of my cat Piccadilly, cycling, and visiting underappreciated tourist sites.

Images of six "little red dot" galaxies from JWST

Editor’s Note: Astrobites is a graduate-student-run organization that digests astrophysical literature for undergraduate students. As part of the partnership between the AAS and astrobites, we occasionally repost astrobites content here at AAS Nova. We hope you enjoy this post from astrobites; the original can be viewed at

Title: A Census of Photometrically Selected Little Red Dots at 4 < z < 9 in JWST Blank Fields
Authors: Vasily Kokorev et al.
First Author’s Institution: Kapteyn Astronomical Institute
Status: Published in ApJ

Imagine you peek into a kindergarten class and, to your shock, you see that all the children are well over 6 feet tall. That is precisely how astronomers felt when data from JWST showed galaxies with massive black holes just a few hundred million years after the birth of the universe. Some of these black holes have been measured to be a million times more massive than the mass of our Sun, and astronomers are puzzled as to how they could have gained so much mass in such a short time.

The earliest galaxies likely to host these black holes show up as little red dots in the images JWST took of the early universe (as seen in the long exposure images; see the bottom image in Figure 1). They are believed to be compact (with a small radius) galaxies with a Type I active galactic nucleus and obscured (covered by dust), which accounts for their red color and why they are easily observed in the infrared. The spectrum of the galaxy has a “V” shape, a blue continuum from the unobscured part of the active galactic nucleus in the galaxy, and a red continuum from the obscured part (see Figure 1). These little red dots have appeared in several images taken by JWST, hinting that plenty of massive black holes are lurking in the early universe.

spectrum and images of a little red dot galaxy

Figure 1: The characteristic spectrum of little red dots (on the left), with the compact source contributing to the dust-free blue color in the continuum (top right) and the dust-reddened part (bottom right). [Adapted from Kokorev et al. 2024]

Where You Look Matters!

It is vital to systematically look at these little red dot galaxies to understand how many massive black holes were in the early universe. Two factors can introduce biases in the counting of these galaxies. One is the phenomenon of cosmic variance: is JWST just preferentially looking in a direction with many little red dots, or should we expect the same number even when it looks at different parts of the universe? The other is how crowded it is in the direction in which you are observing: if you have a lot of stars in the direction you are looking, they could be misclassified as little red dots or vice versa. If you happen to have plenty of massive galaxy clusters in the direction you are looking, they may create an illusion that more of these little red dots exist than their true numbers (a phenomenon known as gravitational lensing).

To minimize the errors caused by these effects in determining the number of little red dots in the early universe, the authors of today’s research article specifically look at large areas (640 arcminutes2) of the sky by combining JWST data from various programs. This would minimize the effects of cosmic variance as you can measure the numbers over a bigger area of the sky. They also look specifically at data in blank fields (defined here as areas on the sky without galaxy clusters), which helps them determine the accurate number of objects per unit volume. All these galaxies are photometrically selected (i.e., chosen only from looking at images rather than spectra), meaning there is limited spectroscopic data to help confirm what kind of objects they are. Galaxies are determined to be little red dots based on their red colors and how compact they are in the images. Using the obtained fluxes, the authors then construct spectral energy distributions to determine the redshifts (z) of the sample. Limiting the sample to z > 4 (for the early universe), the authors end up with 260 little red dot galaxies

Do Not Judge a Galaxy by Its Size (in Your JWST Image!)

On calculating the total luminosity of the little red dots and comparing it to their redshifts, the authors find a large number of bright little red dots at redshift of z = 5 (around a billion years after the beginning of the universe). The number of little red dots is almost 100 times more than the number of ultraviolet-selected quasars, which are active galactic nuclei identified using another method. They also find that computer models are unable to reproduce the high fraction of the bright galaxies they uncover at high redshifts (left side of Figure 2). The authors derive the mass of the black holes at redshifts of z = 4.5–6.5 to be around 106–108 solar masses, indicating that these black holes were already massive a few hundred million years after the Big Bang. They find deviations from the predicted number densities of massive black holes at these redshifts from galaxy simulations. This is likely because the galaxies that host more gigantic black holes are very dusty, and thus, their spectra do not have any blue continuum. They may then be missed from selections of little red dots as one of the factors it depends on is the characteristic “V”-shaped spectrum, which would need a contribution from the blue continuum (right side of Figure 2).

plots showing the number density of little red dots as a function of luminosity and black hole mass

Figure 2: Left: The number density of the little red dots as a function of luminosity at 6.5 < z < 8.5 with the predicted values from simulations indicated by the blue solid line. Right: The number density as a function of black hole mass at 4.5 < z < 6.5. The observed number density of more massive black holes is lower than the values predicted by simulations. [Adapted from Kokorev et al. 2024]

While spectroscopy is a more reliable method to identify massive black holes, many galaxies that host black holes can still be picked out using near-infrared colors and photometry, which is a much less expensive technique. The challenge lies in ensuring that the photometrically selected sources are reliable, and the authors of today’s article made great use of this technique. Follow-up spectroscopic studies of these photometrically selected samples can help us further understand the exact nature of the black holes. Such studies are already underway, and we can look forward to finally making sense of how these black holes became so massive in such a short time after the formation of the universe!

Original astrobite edited by Delaney Dunne.

About the author, Archana Aravindan:

I am a PhD candidate at the University of California, Riverside, where I study black hole activity in small galaxies. When I am not looking through some incredible telescopes, you can usually find me reading, thinking about policy, or learning a cool language!

two images of Dracula's Chivito

Editor’s Note: Astrobites is a graduate-student-run organization that digests astrophysical literature for undergraduate students. As part of the partnership between the AAS and astrobites, we occasionally repost astrobites content here at AAS Nova. We hope you enjoy this post from astrobites; the original can be viewed at

Title: Dracula’s Chivito: Discovery of a Large Edge-On Protoplanetary Disk with Pan-STARRS
Authors: Ciprian T. Berghea et al.
First Author’s Institution: US Naval Observatory
Status: Published in ApJL

Where Planets Are Born

Studying protoplanetary disks helps us understand how planets, including those in our solar system, are born. These disks are vast and flared structures, consisting of dust and gas orbiting a young star. Protoplanetary disks contain the remnants of the stellar birth process, in which a collapsing molecular cloud gives rise to a central star surrounded by a swirling disk of material. Protoplanetary disks are vital to observe as they are the birth sites of planets. The tiny dust particles come together, sticking to each other and forming larger bodies. This process, influenced by gravity, gas, and radiation, leads to the birth of planets in developing planetary systems.

Meet the Vampire Sub

This research article features a large edge-on protoplanetary disk that was stumbled upon when going through images from the Pan-STARRS research project as a part of  a study of active galactic nucleus candidates. This disk is one of the largest known disks in the sky and is oriented edge on, completely obscuring its central star. Associated with a source of infrared light in the same region of the sky, IRAS 23077+6707, the disk spans approximately 11 in apparent size, with a very faint structure in the disk’s northern part extending out to about 17. It is possibly the largest protoplanetary disk (by angular extent) discovered to date. The structure of this disk is reminiscent of the popular Gomez’s Hamburger, which is not associated with any star-forming region, just like the subject of this article. The similarity to a sandwich, along with the fang-like structures in the northern part of the disk as seen in the images in Figure 1, earned the IRAS 23077+6707 protoplanetary disk the name “Dracula’s Chivito” (chivito is a type of sandwich and the national dish of Uruguay, where one of the co-authors is from).

Left: Pan-STARRS image of Dracule’s Chivito with a bright blue-purple light in the center, as a structure resembling a sandwich. It has faint filaments extending out like “fangs”. Right: A model of the disk on the left, which looks smoother and the fangs are shown to be a part of larger disk envelope, which is a hazy blue light emanating from the disk in the center which has a pink hue

Figure 1: Pan-STARRS image of Dracula’s Chivito (left) and the associated radiative transfer model (right). In the northern part of the disk, very faint filaments are seen extending 17″ out from the edges of the disk. These filaments have been nicknamed “fangs.” These are interpreted as the edges of the disk envelope, as reproduced in the model image. These fangs are likely obscured by a cloud in the southern part of the disk. [Adapted from Berghea et al. 2024]

Decoding DraChi

Analysis: Images of Dracula’s Chivito (henceforth referred to as DraChi) were obtained in the grizy filters of Pan-STARRS (Figure 1). These images and data from the Galaxy Evolution Explorer (GALEX) (ultraviolet), the Two Micron All-Sky Survey (2MASS) (infrared), the Infrared Astronomical Satellite (IRAS) (infrared), and AKARI (infrared) were used for photometric analysis, i.e., flux or brightness measurements. These data were used to construct the spectral energy distribution of the disk (Figure 2), together with a radiative transfer model generated using the code HOCHUNK3D. Radiative transfer models help us understand the disk geometry and how light is scattered by the dust grains in the disk. This light, which we see along our line of sight, is plotted as a spectral energy distribution. Spectral energy distributions give information about how much energy an object gives off at different wavelengths.

A plot of the spectral energy distribution of the disk as obtained from brightness measurements from different instruments and compared to the spectrum from the model. The different colored dots representing data are mostly fit well by the model except in near and mid infrared wavelengths (from about 10 to 50 microns)

Figure 2: The spectral energy distribution of the disk, using photometric data from the image, the model, and other sources as mentioned in the article. Brightness is plotted as a function of wavelength. The data (colored dots) match the model (the solid red line shows the model without extinction, and the dashed red line shows the model with extinction) in most places except in near- and mid-infrared wavelengths, which could be due to discrepancies between instruments or possible variability in the disk’s luminosity. [Berghea et al. 2024]

Distance: It is hard to estimate DraChi’s distance because it is not associated with a known star-forming region. Accurate distances to local molecular clouds are vital to locating protoplanetary disks and comprehending planet formation processes. Therefore, using the Gaia DR3 data for nearby stars, the extinction of the disk was estimated to find the distance to the nearest interstellar clouds, and DraChi is hence estimated to be 978 light-years away.

Basic Properties: The spectral energy distribution suggests that the host star of the disk is a pre-main-sequence star of type A with a temperature of about 6500–8500K. The images of the disk and the resultant radiative transfer model constrain the disk inclination to be between 80° and 84°. The scale height of the disk is about 25–50 au at a radius of 500 au (astronomical units), where 1 au is the distance from Earth to the Sun. Using the distance and the angular extent of the disk in the sky, the disk’s radius is estimated to be 1,650 au. The radiative transfer model, based on the scattering of the light, quantifies the mass of the disk to be about 0.2 times that of the Sun.

The Fangs: The authors noticed two “fang-like” features in the northern part of the disk, and these features were also reproduced in the model of the disk. The “fangs” closely resemble the “edge” of the shadow created by the disk in the bright surrounding envelope. They could be filaments due to a possible outflow from the central part of the disk, which is characteristic of a young disk at the end of the Class I phase (~0.5 million years old). It is possible that the fangs are present in the south, but this region is likely obscured in the images from Pan-STARRS and could be perhaps seen in infrared (longer-wavelength) imaging.

Is DraChi the Only One?

The short answer is no. DraChi is certainly different from most other protoplanetary disks, given its size and large distance from any known star-forming regions. But the existence of Gomez’s Hamburger proves there are more such disks awaiting discovery. A disk as young as DraChi is vital to understanding planet formation in its earlier stages, and its large size makes for interesting future observations using more sensitive instruments.

Original astrobite edited by Kylee Carden and Jessie Thwaites.

About the author, Maria Vincent:

Maria is a PhD candidate in astronomy at the Institute for Astronomy, University of Hawai’i at Manoa. Her research focuses on adaptive optics and high-contrast imaging science and instrumentation with ground-based telescopes. Driven by a fascination with planet formation and the intricate processes shaping our solar system, she uses the Subaru Coronagraphic Extreme Adaptive Optics suite to observe and study morphological features of protoplanetary disks in near-infrared wavelengths, aiming to understand disk structure and processes governing planet formation. On the instrumentation side, she is working on designing and constructing an optical testbed to test and characterize a new deformable mirror as part of the upcoming High-order Advanced Keck Adaptive Optics upgrade. Outside of work, she enjoys blogging, mystery, historical and science fiction literature and cinemedia, photography, hiking, and travel.

Illustration of a binary system undergoing mass transfer

Editor’s Note: Astrobites is a graduate-student-run organization that digests astrophysical literature for undergraduate students. As part of the partnership between the AAS and astrobites, we occasionally repost astrobites content here at AAS Nova. We hope you enjoy this post from astrobites; the original can be viewed at

Title: Expansion of Accreting Main-Sequence Stars During Rapid Mass Transfer
Authors: Mike Y. M. Lau et al.
First Author’s Institution: Heidelberg Institute for Theoretical Studies
Status: Published in ApJL

In our universe, most stars live in binaries. These stars live together, orbiting around each other, evolving and influencing each other’s lives. For many stars, this can be a tumultuous relationship with periods of mass transfer, altering the stars’ lives from the typical single-star evolution from a main-sequence star to a red giant to a stellar remnant. Mass transfer allows stars to grow and shrink in mass through their evolution, changing their identities and speeding up or slowing down their evolution. Mass transfer can also cause stars to go supernova.

diagram of the stages of Roche lobe filling

Figure 1: An Illustration of Roche lobe filling of binary stars in three different scenarios, A, B, and C. The left column shows the gravitational potential wells for each star, while the right column outlines the stars and their Roche lobes (the teardrop shapes around them). In scenario A, the matter of both stars is within each potential well and the stars are detached. In scenario B, the left star has expanded and filled its Roche lobe, causing matter to accrete onto the right star. In scenario C, both stars have filled their Roche lobes and orbit in a contact binary. [Storm Colloms, adapted from Pringle and Wade, Interacting Binary Stars (1985, Cambridge University Press)]

Mass transfer occurs when one of the stars in a binary system expands in radius during its evolution. If this expansion pushes the outside layers of the star beyond the gravitational influence of the star, the outer layers will be pulled from the expanded donor star onto the other star in the binary. The region of gravitational influence of each star is called its Roche lobe, illustrated in Figure 1.

As common as mass transfer is in the lives of stars, there are a lot of unknowns about mass transfer processes. Astronomers make simulations of stars transferring mass to try to understand these unknowns better. Simulating a fully accurate model of these processes requires many computational resources to model the hydrodynamic processes in all dimensions, and so models always make some simplifying assumptions.

One big simplifying assumption is how much of the donor star’s mass can be packed onto the accreting star. Previous models assumed that the accreting star does not change with time when accepting mass from the donor star. However, today’s authors are actually able to investigate how the properties of the accreting star change through the mass-transfer period using a stellar evolution code called Modules for Experiments in Stellar Astrophysics (MESA).

The authors show that how the accretor’s mass and radius respond to mass transfer depends on two timescales: 1) the thermal timescale of the accretor — this is the timescale on which the accretor can radiate away all of its thermal energy, and 2) the timescale of the accretion — this is the rate of mass transfer from the donor to the accretor. The thermal timescale of the accreting star is generally much longer than the mass-transfer timescale, meaning that during mass transfer the accreting star is fed mass much more quickly than it can “swallow” it, or have it thermally cool and settle onto the star. However, the accreting star can store this accreted matter in an envelope around it, stuffing it in its cheeks much like a hamster, as long as this does not spill over the edge of the Roche lobe — hence why the authors coin the name “hamstars” for these stars. This “cheek stuffing” makes the accreting star grow in radius and cool much like a typical main-sequence star becoming a red giant, as shown in Figure 2. As the star accretes more and more and its luminosity increases, its thermal timescale decreases, and it can radiate away the energy of accreting mass and “swallow” the matter faster. The star then shrinks from its expanded radius, having grown several times in mass.

HR diagrams showing the evolution of the accreting star during mass transfe

Figure 2: Hertzsprung–Russell diagrams showing the evolution of the accreting star during mass transfer. The three panels show the temperature versus luminosity for accretors of different initial masses, and each track in the panels is a different rate of mass transfer in the simulation. Most of the stars cool at some point in evolution, expanding in radius, before the accreted matter can cool onto the star and it shrinks and heats again, resuming its usual evolution having grown several times in mass. [Lau et al. 2024]

How much the radius expands during mass transfer changes the fraction of mass lost by the donor that is added to the accretor. This is called the accretion efficiency, and this quantity is often assumed to be an arbitrary constant. With the authors’ prescription for the change in radius of the accretor, they describe a new way of calculating the accretion efficiency that considers the accretor’s expansion. This depends on the initial mass of the accreting star, and it’s found to be higher than typical assumptions for low-mass (less than 15 solar masses) accreting stars, and lower for high-mass (more than 20 solar masses) stars.

The accretion efficiency determines how much of the mass accreted actually makes it onto the accreting star, versus how much might be lost during mass transfer. While low-mass stars can pack a lot of accreting matter in their hamster cheeks, high-mass stars have a lower limit for how much matter they can chew. These higher-mass accreting stars make up observed populations of X-ray binaries and the progenitors of pairs of black holes and/or neutron stars. Compared to current models, the authors’ prescription for accretion efficiency is in better agreement with the observed population of X-ray binaries and predicts fewer merging black hole binaries that could be observed with gravitational waves, something we could investigate further with more observations.

While today’s article considers the accretor’s response to mass transfer (pinning down the similarities between stars and rodents), it does not account for the accreting star’s spin or mass loss. If a star is spinning fast, it could limit how much the accreting star can expand. Some mass can also be swept away by stellar winds, meaning that not all of the mass lost by one star is accreted onto the other. Future work is needed to incorporate these additional factors before we can look up into the cosmos for a great hamstar hunt!

Original astrobite edited by Lucie Rowland.

About the author, Storm Colloms:

Storm is a postgraduate researcher at the University of Glasgow, Scotland. They work on understanding populations of binary black holes and neutron stars from the gravitational wave signals emitted when they merge, and what that tells us about the lives and deaths of massive stars. Outwith astrophysics they spend their time taking digital and film photos, and making fun doodles of their research.

1 2 3 39