Blue Large-Amplitude Pulsators (BLAPs)

• BLAPs are hot, compact variable stars characterized by short pulsation periods (from a few to about 75 minutes), high temperatures (around 25,000–35,000 K), and asymmetric sawtooth light curves.
• Their pulsations are driven by the iron-group opacity bump, with period–gravity and period–luminosity relations strongly supporting a dominant fundamental radial mode.
• Multiple evolutionary channels—including binary stripping, helium burning, and merger scenarios—explain the diverse properties and rapid period changes observed in BLAPs.

Searching arXiv for recent and foundational BLAP papers to ground the article in the literature. Searching arXiv for Blue Large-Amplitude Pulsators and closely related formation-channel studies. Querying arXiv for “Blue Large-Amplitude Pulsators BLAP”. Blue Large-Amplitude Pulsators (BLAPs) are a recently identified class of hot, compact variable stars, usually interpreted as radially pulsating objects with effective temperatures $T_{\rm eff}\simeq 20{,}000$–$35{,}000$ K, pulsation periods from a few minutes to about 75 minutes, relatively low surface gravities for their temperature, and large-amplitude, non-sinusoidal light curves (Zhang et al., 26 Sep 2025). First recognized in OGLE data as a new class of short-period blue variables, they occupy a region of the Hertzsprung–Russell diagram between hot subdwarfs and upper main-sequence stars and have become a focal point for work on stellar pulsation, binary stripping, helium burning, and merger remnants (Pietrukowicz et al., 2017). Although their observational coherence is well established, their mode identification and evolutionary origin remain active problems, and multiple formation channels are now supported in the literature (Jadlovský et al., 2024).

1. Discovery and observational definition

The class entered the literature through OGLE as more than a dozen previously unknown short-period variables with periods in the range of 20–40 min and amplitudes of 0.2–0.4 mag in the optical passbands (Pietrukowicz et al., 2017). The original sample quickly proved to be physically distinct from $\delta$ Scuti stars: the light curves were sawtooth-shaped, the colors were very blue, and follow-up spectroscopy showed temperatures near $30{,}000$ K rather than the much cooler atmospheres expected for classical short-period main-sequence pulsators. The prototype objects were also found to be helium enriched and to have surface gravities around the hot-subdwarf regime, but not identical to ordinary sdB or sdOB stars (Pietrukowicz, 2018).

Observationally, BLAPs are defined by a linked set of properties: very short pulsation periods, large optical amplitudes, blue colors, hot atmospheres, and asymmetric phased light curves with a fast rise and slow decline. Their atmospheric parameters typically fall around $T_{\rm eff}\approx 25{,}000$–$34{,}000$ K and $\log g \approx 4.2$–$5.7$, while the original “classical” BLAPs concentrate near 20–40 min and lower gravities, and the later high-gravity subgroup occupies roughly 2–8 min or 3–8 min with $\log g\sim 5.3$–$5.7$ (Pietrukowicz et al., 2024, McWhirter et al., 2022).

The empirical boundaries of the class broadened substantially as searches expanded beyond the initial OGLE sample. A major OGLE re-analysis reported periods from 5.36 minutes to 76.36 minutes, including 14 BLAPs with periods below 10 minutes and five with periods longer than 60 minutes, bringing the total known BLAP sample to almost 200 (Borowicz et al., 23 Oct 2025). This undermined the earlier impression of a clean period gap between the original “classical” BLAPs and the shorter-period, higher-gravity objects. Earlier hints that the gap might not be physically real had already come from transitional objects such as TMTS-BLAP-1 and OW-BLAP-1, which occupy the interval between the two observational groups (Zhang et al., 2023, Lin et al., 2022).

2. Pulsation phenomenology and mode content

BLAP light curves are characteristically asymmetric and non-sinusoidal. The classical morphology is sawtooth-like, with a steep rise to maximum brightness and a slower decline, reminiscent of fundamental-mode Cepheids and RR Lyrae stars, but occurring at far higher temperatures (Pietrukowicz et al., 2017). Subsequent OGLE work showed that the morphology varies systematically with period: shorter-period objects tend to have more rounded, symmetric light curves, whereas many longer-period BLAPs exhibit an additional bump or a dip near maximum light (Pietrukowicz et al., 2024).

The dominant oscillation is generally treated as a single radial mode. High-cadence Liverpool Telescope photometry of OGLE-BLAP-009 and OGLE-BLAP-014 recovered only one astrophysical period in each object, $35{,}000$0 min and $35{,}000$1 min, with no additional short-period variability down to $35{,}000$2 mmag for BLAP-009 and $35{,}000$3 mmag for BLAP-014 (McWhirter et al., 2020). More generally, long OGLE baselines indicate that most BLAPs are single-mode stars, with additional periodicities detected only rarely; OGLE-BLAP-030 is an important exception, with triple-mode behaviour (Pietrukowicz et al., 2024).

The excitation mechanism is usually identified with the $35{,}000$4-mechanism associated with the iron-group opacity bump. Several studies place the relevant driving region near $35{,}000$5 K or $35{,}000$6–5.35, and emphasize that radiative levitation of iron and nickel can be critical for creating sufficient opacity in the envelope (Zhang et al., 2023, Zhang et al., 2024). This links BLAPs to other hot pulsators near the hot-subdwarf branch, but their periods, amplitudes, and global parameters place them in a distinct instability domain.

The preferred mode identification remains slightly more complicated than the observational shorthand suggests. A strong line of work favors the radial fundamental mode on the basis of light-curve shape, radial-velocity behaviour, period relations, and linear stability calculations (Pietrukowicz et al., 2024, Wu et al., 25 Oct 2025). However, low-mass pre-white-dwarf models have also been used to argue that high-order nonradial $35{,}000$7 modes can reproduce the sign and magnitude of observed $35{,}000$8 more naturally than the radial fundamental mode, even while the morphology appears more radial-like (Córsico et al., 2018). The resulting tension has not disappeared, although more recent observational syntheses increasingly treat the dominant mode as the fundamental radial mode in most stars (Pietrukowicz et al., 2024).

3. Spectroscopic properties, empirical relations, and secular variability

Spectroscopy shows that BLAPs form a comparatively homogeneous family in $35{,}000$9 and $\delta$0, but not in surface composition. A MagE-based study of OGLE BLAPs found effective temperatures between about $\delta$1 and $\delta$2 K and surface gravities between about 4.3 and 5.7, while also identifying a split into helium-rich and helium-poor groups (Pietrukowicz et al., 2024). In that sample, 13 of 15 spectroscopic targets were helium enriched with $\delta$3 roughly between $\delta$4 and $\delta$5, whereas two objects had helium abundances below $\delta$6 dex. The He-poor stars cluster preferentially at the hotter, shorter-period end, although the separation is not absolute.

The helium-rich stars are also metal enriched. Co-added spectra reveal strong lines from C, N, O, Mg, Al, and Si, with modeling indicating approximate enrichments of carbon, magnesium, and aluminum by about $\delta$7 solar, nitrogen by about $\delta$8 solar, oxygen by about $\delta$9 solar, and silicon by roughly $30{,}000$0–$30{,}000$1 solar (Pietrukowicz et al., 2024). This composition pattern is important because BLAP pulsation is thought to be driven by the iron-group opacity bump; metal-enriched envelopes are therefore physically consistent with the instability mechanism.

Two empirical relations have become especially prominent. The first is the period–gravity relation,

$30{,}000$2

with

$30{,}000$3

for $30{,}000$4 in minutes (Pietrukowicz et al., 2024). The second is a period–luminosity relation anchored with OGLE-BLAP-009,

$30{,}000$5

Both relations are widely taken as support for radial pulsation, and model grids using MESA-RSP found that the observed period–luminosity, period–gravity, and period–radius relations agree better with the low-mass scenario than with higher-mass alternatives (Jadlovský et al., 2024).

Secular period change is another defining empirical feature. Long-term monitoring shows that typical BLAP period-change rates are of order $30{,}000$6, with both positive and negative signs represented (Pietrukowicz et al., 2017). In one large compilation, 35 objects had measured $30{,}000$7, of which 10 were negative and 25 positive, spanning $30{,}000$8 to $30{,}000$9 (Wu et al., 25 Oct 2025). OGLE-BLAP-030 is an extreme case: one study measured a dominant-period change of about $T_{\rm eff}\approx 25{,}000$0 (Pietrukowicz et al., 2024), while later OGLE work identified three objects with even larger rates of order $T_{\rm eff}\approx 25{,}000$1, including $T_{\rm eff}\approx 25{,}000$2, $T_{\rm eff}\approx 25{,}000$3, and $T_{\rm eff}\approx 25{,}000$4 (Borowicz et al., 23 Oct 2025). These values imply rapid structural evolution in at least a subset of the class.

4. Evolutionary interpretations and formation channels

No single evolutionary channel explains all observed BLAPs. The literature instead supports several viable pathways that converge on similar regions of the $T_{\rm eff}\approx 25{,}000$5–$T_{\rm eff}\approx 25{,}000$6 plane.

Channel	Representative masses	Distinctive result
He-core pre-WD / pre-ELM WD	$T_{\rm eff}\approx 25{,}000$7–$T_{\rm eff}\approx 25{,}000$8	Binary-stripped, low-mass interpretation
Helium-burning binaries	$T_{\rm eff}\approx 25{,}000$9–$34{,}000$0 BLAPs from $34{,}000$1–$34{,}000$2 progenitors	Long residence in BLAP region
HeWD+MS merger	$34{,}000$3–$34{,}000$4	Shell flashes produce both signs of $34{,}000$5
Double low-mass WD merger	total mass $34{,}000$6–$34{,}000$7	Short-lived post-merger BLAPs
Shell-He-burning subdwarf	$34{,}000$8	Large positive $34{,}000$9 in transitional objects

The low-mass pre-white-dwarf channel was the first quantitatively developed explanation. Fully evolutionary pre-ELM WD models place BLAPs at the hot end of binary-stripped He-core tracks, with masses typically near $\log g \approx 4.2$0, and show that both high-order nonradial $\log g \approx 4.2$1 modes and, for the shortest periods, low-order radial modes can fall in the observed period range (Romero et al., 2018). A detailed case study of OGLE-BLAP-009 strengthened that picture: spectroscopy, Gaia parallax, and SED fitting gave $\log g \approx 4.2$2, $\log g \approx 4.2$3, and $\log g \approx 4.2$4, while MESA+GYRE matching favored a first-overtone solution with $\log g \approx 4.2$5 (Bradshaw et al., 2023).

A second major channel invokes helium-burning stars. Binary population synthesis with BPASS finds that helium-burning BLAP progenitors come mainly from intermediate-mass binaries with initial masses $\log g \approx 4.2$6–$\log g \approx 4.2$7, producing BLAPs with final masses $\log g \approx 4.2$8–$\log g \approx 4.2$9 (Zhang et al., 26 Sep 2025). These models occupy the pulsation region for longer than previously studied pre-white-dwarf models and enter the BLAP regime within a narrow age range, $5.7$0–$5.7$1. This channel therefore provides a major Galactic contribution even if it is not dominant in every parameter regime.

Merger channels are also well developed. One study showed that mergers of a helium-core white dwarf with a low-mass main-sequence star evolve through the observed BLAP region between helium-shell ignition and full helium-core burning, with plausible BLAP masses $5.7$2–$5.7$3 (Zhang et al., 2023). Flash-driven expansion and contraction loops in these models naturally reproduce both positive and negative $5.7$4, a long-standing constraint on BLAP evolution. Another study treated mergers of two low-mass white dwarfs and found that BLAPs can form if the total merger mass lies between $5.7$5 and $5.7$6; these remnants become BLAPs within $5.7$7–$5.7$8 yr after coalescence and spend only $5.7$9–$\log g\sim 5.3$0 yr in the BLAP region (Kołaczek-Szymański et al., 2024).

Individual objects can discriminate between channels. TMTS-BLAP-1, with $\log g\sim 5.3$1 min and $\log g\sim 5.3$2, is difficult to reconcile with contracting pre-WD models or ordinary core-helium-burning models, but is consistent with a shell-helium-burning subdwarf in a short-lived shell-helium-ignition phase (Lin et al., 2022). A plausible implication is that the BLAP instability strip is populated by more than one kind of stripped or post-merger star rather than by a single homogeneous evolutionary product.

5. Binary systems, companion channels, and Galactic populations

Binary evidence is central to BLAP formation theory, even though most known BLAPs are not obvious eclipsing or interacting binaries. The first confirmed BLAP in a binary system is HD 133729, a system containing a late B-type main-sequence star and a BLAP secondary (Pigulski et al., 2022). The BLAP pulsates with a period of 32.37 min, and the $\log g\sim 5.3$3 diagram shows a light-travel-time modulation with $\log g\sim 5.3$4 d. Because the B-type primary dominates the optical light, the observed TESS pulsation amplitude is diluted to about 31.6 ppt; after correcting for flux dilution, the intrinsic BLAP amplitude is about $\log g\sim 5.3$5 mag. If the components are coeval and no mass transfer occurred, low-mass progenitor scenarios are excluded for this system.

A dedicated binary-evolution study of HD 133729 instead found that a pre-white-dwarf Roche-lobe-overflow channel can reproduce the observed configuration. Using MESA, the best-fitting case had $\log g\sim 5.3$6, $\log g\sim 5.3$7, $\log g\sim 5.3$8, $\log g\sim 5.3$9, and $5.7$0 (Zhang et al., 2024). The model predicts helium and nitrogen enhancement on the surface of the B-type companion and implies that the system will later enter a common-envelope phase and merge. HD 133729 therefore functions as a benchmark for BLAPs with main-sequence companions, but it does not force the entire class into a single channel.

Population synthesis predicts substantial Galactic numbers for more than one binary pathway. For low-mass pre-white-dwarf BLAPs, BPASS-based synthesis found 11,931 systems in the Milky Way under a constant star formation rate of $5.7$1, with 78.5 per cent forming after stable RLOF and 21.5 per cent after common-envelope evolution (Byrne et al., 2021). About 75 per cent of these systems have main-sequence companions and about 25 per cent evolved or compact companions. For helium-burning BLAPs, BPASS predicts approximately 14,351 Galactic systems at solar metallicity, including 12,799 with main-sequence companions and 1,551 with evolved or compact-object companions (Zhang et al., 26 Sep 2025).

The helium-burning study emphasizes a sharp binary dichotomy. BLAPs with main-sequence companions form mainly through Roche lobe overflow and tend to have orbital periods around $5.7$2 days, with companion masses peaking near $5.7$3 (Zhang et al., 26 Sep 2025). BLAPs with evolved or compact companions form mainly through common-envelope evolution and have shorter orbital periods peaking around $5.7$4 days; their companions are often white dwarfs with masses below $5.7$5, while NS/BH companions are rare. This RLOF–CEE period dichotomy is one of the clearest observational tests proposed for the helium-burning channel.

6. Surveys, Galactic distribution, and unresolved issues

BLAP detectability is limited by both rarity and selection effects. Known BLAPs are usually faint and concentrated in the Galactic plane, where interstellar extinction can mask their intrinsically blue colors (McWhirter et al., 2022). The Gaia+ZTF search strategy therefore relies on dereddening Gaia photometry with a 3D dust map and then searching for short-period variability. In one such search, 162,377,536 Gaia sources were processed, 608,219 intrinsically blue subluminous stars were retained after color–magnitude cuts, 44,847 variable candidates survived ZTF filtering, and 22 candidate BLAPs were identified after periodicity analysis and visual inspection (McWhirter et al., 2022).

OGLE remains the dominant source of BLAPs because of its long baselines and large source counts. A later OGLE survey searched more than 400 million $5.7$6-band light curves in the inner Galactic bulge and recovered all previously known inner-bulge BLAPs while adding 88 new objects (Borowicz et al., 23 Oct 2025). The OmegaWhite survey likewise demonstrated that BLAP-like objects can be recovered efficiently in high-cadence Galactic-plane imaging, classifying four stars as BLAPs and showing that spectroscopic follow-up at $5.7$7 is needed to separate them from hot subdwarfs and other short-period blue variables (Ramsay et al., 2022).

The Galactic census is still incomplete. The helium-burning BPASS study predicts that most BLAPs are very faint, with apparent magnitudes $5.7$8, because they lie in or near the Galactic plane and suffer strong extinction (Zhang et al., 26 Sep 2025). Under those assumptions, only future wide-field surveys such as WFST and VRO LSST are expected to detect roughly $5.7$9 and $35{,}000$00 BLAPs, respectively. This large difference between observed and predicted populations is not taken as evidence against BLAP formation models; rather, it is treated as a consequence of survey depth, cadence, extinction, and crowding.

Environmental constraints add another unresolved dimension. A dedicated OGLE search through 77.9 million stellar detections in the Magellanic System found no BLAPs (Pietrukowicz, 2018). The paper connected this non-detection to the idea that BLAP instability in hydrogen-deficient atmospheres may prefer high metallicity, which is consistent with their concentration in the Galactic disk and bulge and their absence in the halo and globular clusters. This remains a suggestion rather than a closed result, because completeness, crowding, and distance effects are also relevant.

A further complication is that some BLAPs may be magnetic oblique pulsators rather than strictly single-period radial pulsators in the simplest sense. OGLE-BLAP-001 and ZGP-BLAP-08 both show three frequencies that are exactly equally spaced in frequency, a pattern interpreted as an oblique pulsator signature with inferred rotation periods of $35{,}000$01 d and $35{,}000$02 d, respectively (Pigulski et al., 2024). The same work argues that these stars are likely magnetic and that merger scenarios could provide a natural origin for such fields. This does not overturn the broader view that BLAPs are predominantly single-mode radial pulsators, but it shows that the class probably contains specialized subpopulations as well as multiple formation pathways.