Papers
Topics
Authors
Recent
Search
2000 character limit reached

CSS161010: Luminous Fast Blue Optical Transient

Updated 9 July 2026
  • CSS161010 is a fast blue optical transient located in a dwarf galaxy, characterized by a rapid optical rise, high peak luminosity, and persistent blue colors.
  • Multi-wavelength observations reveal an evolving shock front with mildly relativistic ejecta, detailed via synchrotron self-absorption and equipartition modeling.
  • Competing models—including engine-driven explosions, tidal disruption by an intermediate-mass black hole, and CEJSN jet interactions—offer diverse interpretations of its origin.

CSS161010 is a nearby luminous fast blue optical transient (LFBOT), or fast blue optical transient (FBOT), at redshift z0.034z\approx0.034 in a dwarf galaxy at a luminosity distance of about 150150\,Mpc. It is distinguished by an extremely rapid optical rise, a high peak bolometric luminosity, persistently blue colors, luminous radio and X-ray emission, and an unusual spectroscopic evolution from a blue featureless continuum to very broad, entirely blueshifted hydrogen emission. The event has therefore become a reference case for several competing compact-object-powered interpretations, including an engine-driven aspherical stellar explosion, a tidal disruption event (TDE) by an intermediate-mass black hole (IMBH), and a Common-Envelope Jets Supernova (CEJSN)-like jet interaction scenario (Gutiérrez et al., 2024, Coppejans et al., 2020, Akashi et al., 2020).

1. Discovery, localization, and host-galaxy context

CSS161010 was discovered by CRTS on 2016 Oct 10 and independently by ASAS-SN, with the outburst start defined as JD 2457669.92±2.002457669.92 \pm 2.00. Its reported coordinates are consistent with a source in the dwarf galaxy WISEA J045834.37–081804.4, and its host redshift was measured from narrow H and [O III] emission to be z=0.0340±0.0006z = 0.0340 \pm 0.0006, corresponding to DL149.4D_L \approx 149.4\,Mpc (Gutiérrez et al., 2024).

The host environment is itself an important part of the classification problem. Prospector fitting of the GTC/OSIRIS spectrum plus BgVriz spectral energy distribution yielded log(M/M)=8.1350.079+0.087\log(M_*/M_\odot)=8.135^{+0.087}_{-0.079}, a recent star-formation rate of about 0.0150.008+0.010Myr10.015^{+0.010}_{-0.008}\,M_\odot\,{\rm yr}^{-1}, stellar metallicity Z0.0260.010+0.007ZZ_* \approx 0.026^{+0.007}_{-0.010}\,Z_\odot, and gas-phase metallicity 0.23Z\approx 0.23\,Z_\odot (Gutiérrez et al., 2024). A separate fit based on deep Keck/LRIS spectroscopy and multi-band photometry gave M107MM_*\approx10^7\,M_\odot, 150150\,0, and 150150\,1, placing the host in the extreme dwarf-starburst regime (Coppejans et al., 2020). These differing host-mass estimates reflect modeling choices rather than a settled inconsistency.

The host has also been used to motivate an IMBH interpretation. Extrapolating 150150\,2–150150\,3 relations suggests 150150\,4, i.e. an intermediate-mass black hole (Gutiérrez et al., 2024). Another analysis noted that CSS161010 lies within 150150\,5 of its host nucleus and found no prior variability or AGN lines, while host narrow-line diagnostics place the galaxy in the star-forming regime rather than the AGN regime (Coppejans et al., 2020).

2. Optical and spectroscopic phenomenology

The optical evolution is unusually fast even within the LFBOT class. The V-band light curve reaches an absolute peak of 150150\,6mag in 150150\,7days, and the bolometric light curve peaks at 150150\,8 at about the same epoch, with total radiated energy 150150\,9erg (Gutiérrez et al., 2024). After maximum, the optical bands follow a single power law,

2457669.92±2.002457669.92 \pm 2.000

and the intrinsic colors remain very blue for about 2457669.92±2.002457669.92 \pm 2.001days, with 2457669.92±2.002457669.92 \pm 2.002mag and 2457669.92±2.002457669.92 \pm 2.003mag (Gutiérrez et al., 2024).

Blackbody fits show roughly constant 2457669.92±2.002457669.92 \pm 2.004K from 2457669.92±2.002457669.92 \pm 2.005 to 2457669.92±2.002457669.92 \pm 2.006days. Over the same interval, the inferred photospheric radius rises from about 2457669.92±2.002457669.92 \pm 2.007cm at 2457669.92±2.002457669.92 \pm 2.008days to about 2457669.92±2.002457669.92 \pm 2.009cm at peak, then declines to about z=0.0340±0.0006z = 0.0340 \pm 0.00060cm by z=0.0340±0.0006z = 0.0340 \pm 0.00061days (Gutiérrez et al., 2024). Another observational summary likewise emphasized blue spectra with an effective temperature z=0.0340±0.0006z = 0.0340 \pm 0.00062K and a decay timescale z=0.0340±0.0006z = 0.0340 \pm 0.00063days (Gao et al., 25 Aug 2025).

The spectroscopic sequence is the most unusual aspect of the event. At z=0.0340±0.0006z = 0.0340 \pm 0.00064days the spectrum is a blue, featureless continuum with z=0.0340±0.0006z = 0.0340 \pm 0.00065K. At z=0.0340±0.0006z = 0.0340 \pm 0.00066–z=0.0340±0.0006z = 0.0340 \pm 0.00067days, broad He II z=0.0340±0.0006z = 0.0340 \pm 0.00068 and He II z=0.0340±0.0006z = 0.0340 \pm 0.00069 emerge at DL149.4D_L \approx 149.4\,0. From DL149.4D_L \approx 149.4\,1days onward, the spectra become dominated by very broad, entirely blueshifted HDL149.4D_L \approx 149.4\,2 and HDL149.4D_L \approx 149.4\,3 emission, with no rest-wavelength or redshifted component ever seen (Gutiérrez et al., 2024).

The hydrogen-line kinematics are extreme. The emission peak remains blueshifted by at least DL149.4D_L \approx 149.4\,4 at all epochs; the maximum ejecta velocity traced by the blue wing declines from about DL149.4D_L \approx 149.4\,5 (DL149.4D_L \approx 149.4\,6) at DL149.4D_L \approx 149.4\,7days to about DL149.4D_L \approx 149.4\,8 (DL149.4D_L \approx 149.4\,9) at log(M/M)=8.1350.079+0.087\log(M_*/M_\odot)=8.135^{+0.087}_{-0.079}0days; and the FWHM drops from about log(M/M)=8.1350.079+0.087\log(M_*/M_\odot)=8.135^{+0.087}_{-0.079}1 to about log(M/M)=8.1350.079+0.087\log(M_*/M_\odot)=8.135^{+0.087}_{-0.079}2 (Gutiérrez et al., 2024). The Hlog(M/M)=8.1350.079+0.087\log(M_*/M_\odot)=8.135^{+0.087}_{-0.079}3 luminosity rises to about log(M/M)=8.1350.079+0.087\log(M_*/M_\odot)=8.135^{+0.087}_{-0.079}4 at about log(M/M)=8.1350.079+0.087\log(M_*/M_\odot)=8.135^{+0.087}_{-0.079}5days, while log(M/M)=8.1350.079+0.087\log(M_*/M_\odot)=8.135^{+0.087}_{-0.079}6 declines from about log(M/M)=8.1350.079+0.087\log(M_*/M_\odot)=8.135^{+0.087}_{-0.079}7 to about log(M/M)=8.1350.079+0.087\log(M_*/M_\odot)=8.135^{+0.087}_{-0.079}8 between log(M/M)=8.1350.079+0.087\log(M_*/M_\odot)=8.135^{+0.087}_{-0.079}9 and 0.0150.008+0.010Myr10.015^{+0.010}_{-0.008}\,M_\odot\,{\rm yr}^{-1}0days (Gutiérrez et al., 2024). The persistent blueshifted emission and the lack of any emission at the rest wavelength were described as unique features not seen in any transient before CSS161010 (Gutiérrez et al., 2024).

3. Radio, X-ray, and mildly relativistic ejecta

CSS161010 was monitored in X-rays with Chandra at rest-frame times 0.0150.008+0.010Myr10.015^{+0.010}_{-0.008}\,M_\odot\,{\rm yr}^{-1}1, 0.0150.008+0.010Myr10.015^{+0.010}_{-0.008}\,M_\odot\,{\rm yr}^{-1}2, and 0.0150.008+0.010Myr10.015^{+0.010}_{-0.008}\,M_\odot\,{\rm yr}^{-1}3days. At 0.0150.008+0.010Myr10.015^{+0.010}_{-0.008}\,M_\odot\,{\rm yr}^{-1}4days the inferred unabsorbed 0.0150.008+0.010Myr10.015^{+0.010}_{-0.008}\,M_\odot\,{\rm yr}^{-1}5–0.0150.008+0.010Myr10.015^{+0.010}_{-0.008}\,M_\odot\,{\rm yr}^{-1}6keV luminosity was 0.0150.008+0.010Myr10.015^{+0.010}_{-0.008}\,M_\odot\,{\rm yr}^{-1}7; at 0.0150.008+0.010Myr10.015^{+0.010}_{-0.008}\,M_\odot\,{\rm yr}^{-1}8days it brightened to 0.0150.008+0.010Myr10.015^{+0.010}_{-0.008}\,M_\odot\,{\rm yr}^{-1}9; and by Z0.0260.010+0.007ZZ_* \approx 0.026^{+0.007}_{-0.010}\,Z_\odot0days it had faded below a Z0.0260.010+0.007ZZ_* \approx 0.026^{+0.007}_{-0.010}\,Z_\odot1 upper limit of Z0.0260.010+0.007ZZ_* \approx 0.026^{+0.007}_{-0.010}\,Z_\odot2 (Coppejans et al., 2020). The sparse X-ray light curve suggests a relatively slow decay or plateau to about Z0.0260.010+0.007ZZ_* \approx 0.026^{+0.007}_{-0.010}\,Z_\odot3days followed by a rapid fade (Coppejans et al., 2020).

Radio observations with the VLA and GMRT span Z0.0260.010+0.007ZZ_* \approx 0.026^{+0.007}_{-0.010}\,Z_\odot4–Z0.0260.010+0.007ZZ_* \approx 0.026^{+0.007}_{-0.010}\,Z_\odot5days at Z0.0260.010+0.007ZZ_* \approx 0.026^{+0.007}_{-0.010}\,Z_\odot6–Z0.0260.010+0.007ZZ_* \approx 0.026^{+0.007}_{-0.010}\,Z_\odot7GHz. The centimeter light curve peaks at about Z0.0260.010+0.007ZZ_* \approx 0.026^{+0.007}_{-0.010}\,Z_\odot8days with a peak spectral luminosity Z0.0260.010+0.007ZZ_* \approx 0.026^{+0.007}_{-0.010}\,Z_\odot9 at 0.23Z\approx 0.23\,Z_\odot0GHz (Coppejans et al., 2020). At 0.23Z\approx 0.23\,Z_\odot1days the spectrum is synchrotron self-absorbed (SSA), with optically thick and thin slopes

0.23Z\approx 0.23\,Z_\odot2

implying an electron index 0.23Z\approx 0.23\,Z_\odot3 (Coppejans et al., 2020).

The radio evolution is steep. After peak, the spectral turnover evolves as 0.23Z\approx 0.23\,Z_\odot4 and 0.23Z\approx 0.23\,Z_\odot5, corresponding to a fixed-frequency decline such as 0.23Z\approx 0.23\,Z_\odot6 (Coppejans et al., 2020). This was interpreted as a strongly decelerating blast wave in a steep density gradient (Coppejans et al., 2020).

Equipartition-based SSA modeling inferred an apparent transverse speed of about 0.23Z\approx 0.23\,Z_\odot7 at 0.23Z\approx 0.23\,Z_\odot8days, about 0.23Z\approx 0.23\,Z_\odot9 at M107MM_*\approx10^7\,M_\odot0days, and a deceleration to about M107MM_*\approx10^7\,M_\odot1 by M107MM_*\approx10^7\,M_\odot2days, implying an initial outflow with M107MM_*\approx10^7\,M_\odot3 (Coppejans et al., 2020). The lower-limit internal energy at M107MM_*\approx10^7\,M_\odot4days, M107MM_*\approx10^7\,M_\odot5erg, implies M107MM_*\approx10^7\,M_\odot6erg in the fast ejecta, while the mass coupled to the mildly relativistic component is M107MM_*\approx10^7\,M_\odot7–M107MM_*\approx10^7\,M_\odot8 (Coppejans et al., 2020). This combination was summarized as a mildly relativistic, high-energy, baryon-loaded outflow or “dirty fireball,” consistent with the absence of detected prompt M107MM_*\approx10^7\,M_\odot9-rays (Coppejans et al., 2020).

The X-ray and radio components were not found to be coextensive in a simple synchrotron extrapolation. Once the synchrotron cooling break is included, the observed X-ray flux lies orders of magnitude above the radio extrapolation, requiring a second emission component (Coppejans et al., 2020).

4. Shock interpretation, revised synchrotron inference, and high-energy particle acceleration

A later theoretical treatment revised the standard equipartition framework by incorporating relativistic corrections and a thermal Maxwell–Jüttner electron population in addition to a non-thermal power law. Applied to CSS161010 at 150150\,00days, with 150150\,01GHz and 150150\,02, the event falls in the thermal-SSA regime rather than the classical pure power-law SSA regime (Margalit et al., 2024). In that interpretation, the shock velocity is 150150\,03 or 150150\,04, the circumsource density parameter is 150150\,05, the physical density at the shock is about 150150\,06, and the post-shock energy is about 150150\,07erg (Margalit et al., 2024).

The same work emphasized that the conventional non-relativistic equipartition formula would overestimate the velocity, density, and energy for a source like CSS161010. For the same observed 150150\,08, the standard treatment gives 150150\,09, 150150\,10, and 150150\,11erg, i.e. about a 150150\,12 overestimate in velocity, a factor of about 150150\,13 in density, and an order-of-magnitude overestimate in energy (Margalit et al., 2024). This suggests that some of the physical tension in early modeling depends on the adopted synchrotron microphysics.

A separate particle-acceleration study combined PIC simulations with Monte Carlo transport and magnetic-field amplification to model non-thermal emission in mildly relativistic FBOT shocks. For CSS161010, the best-fit model at 150150\,14days used 150150\,15cm, 150150\,16G, 150150\,17, proton kinetic energy 150150\,18erg, and electron kinetic energy 150150\,19erg, reproducing the observed self-absorbed radio spectra at 150150\,20, 150150\,21, and 150150\,22days within the measurement errors (Romansky et al., 2024). In the same framework, a synchrotron model gives a 150150\,23–150150\,24keV flux 150150\,25, in agreement with Chandra, and predicts proton energies up to 150150\,26GeV, with a possible upper limit of 150150\,27PeV (Romansky et al., 2024).

5. Competing origin models

The physical origin of CSS161010 remains unsettled. One line of interpretation treats it as an engine-driven aspherical stellar explosion from a hydrogen-rich progenitor. On this view, the inferred kinetic energy and mildly relativistic velocity exceed the expectations of a spherically symmetric core-collapse shock, the large baryon load is inconsistent with classical GRB jets, and the progenitor likely retained a substantial hydrogen envelope (Coppejans et al., 2020). Possible central engines include accretion onto a newly formed black hole or spin-down of a nascent magnetar embedded within a hydrogen-rich envelope (Coppejans et al., 2020).

A second line of interpretation argues that the combined optical and spectroscopic properties favor a TDE by an IMBH. In this reading, stellar-explosion scenarios fail because a radioactive Arnett model requires 150150\,28 and 150150\,29, which is unphysical, while Type Ibn-like circumstellar interaction and Wolf–Rayet/compact-object merger scenarios do not reproduce the pure hydrogen spectra, weak narrow lines, radius evolution, or entirely blueshifted profiles (Gutiérrez et al., 2024). The TDE model is motivated by the host mass scale, the fallback-timescale estimate 150150\,30days with corresponding 150150\,31, the requirement of super-Eddington accretion given 150150\,32, and the possibility that the hydrogen lines arise in an optically thick outflow with 150150\,33 where self-absorption or geometry removes the red wing (Gutiérrez et al., 2024).

A third model invokes jet interaction after a stellar-envelope phase. Three-dimensional hydrodynamical simulations of “late” opposite jets launched into core-collapse supernova ejecta produce hot, low-density bipolar bubbles that cross the photosphere earlier in the polar directions, enhancing radiation leakage and creating viewing-angle-dependent light-curve structure (Akashi et al., 2020). When scaled from jet launch at 150150\,34days to 150150\,35–150150\,36days, as proposed for a neutron star accreting in a CEJSN, the polar bubble reaches the photosphere at 150150\,37–150150\,38days; the high-energy, long-duration run predicts 150150\,39–150150\,40days, 150150\,41, a rapid post-breakout drop for equatorial viewing with 150150\,42days, and a transient blue color 150150\,43–150150\,44K (Akashi et al., 2020). That work states that these simulations naturally explain the timing, luminosity, color, and rapid post-peak decline of CSS161010 within the CEJSN framework (Akashi et al., 2020).

The coexistence of these interpretations is itself a central feature of the source. CSS161010 satisfies defining FBOT phenomenology, but its hydrogen-rich, entirely blueshifted line profiles and its luminous radio/X-ray behavior make it atypical even within that class.

6. Later tests, population context, and multimessenger constraints

CSS161010 has also been used as a target for post-explosion fast radio burst (FRB) searches motivated by newborn-magnetar FBOT models. FAST observed the source twice in late 2022 for a total on-source time of 150150\,45hr at 150150\,46MHz, with a single-pulse search over 150150\,47–150150\,48, boxcar widths up to about 150150\,49ms, and threshold 150150\,50; no FRB-like pulses were found (Gao et al., 25 Aug 2025). The best-case sensitivity reached 150150\,51mJy, and synthetic injections confirmed near-150150\,52 completeness above 150150\,53mJy for DMs up to 150150\,54 (Gao et al., 25 Aug 2025).

Assuming a log-normal luminosity function analogous to repeating FRB 121102, the nondetection implies a burst rate 150150\,55 (Gao et al., 25 Aug 2025). The same study estimated a free–free escape time 150150\,56yr for fiducial FBOT ejecta parameters, shorter than the 150150\,57yr interval between explosion and observation, so the lack of detection was not attributed to bulk ejecta opacity (Gao et al., 25 Aug 2025). The proposed alternatives were weaker bursts, a lower repetition rate than FRB 121102, geometric beaming, propagation effects in the immediate magnetar environment, or an FRB-active phase outside the observed 150150\,58–150150\,59yr window (Gao et al., 25 Aug 2025).

Within the transient population, CSS161010 belongs to the luminous, fast-rising end of the FBOT class alongside AT 2018cow and ZTF18abvkwla. A rate estimate for such luminous FBOTs gave 150150\,60, i.e. 150150\,61 of the local core-collapse supernova rate (Coppejans et al., 2020). One analysis concluded that CSS161010 establishes the existence of a new class of rare H-rich transients that can launch mildly relativistic outflows (Coppejans et al., 2020). Another concluded that its host, light-curve timescale, super-Eddington luminosity, and unique hydrogen-line phenomenology strongly support a partial TDE by an IMBH (Gutiérrez et al., 2024). The CEJSN late-jet interpretation adds explicit predictions of viewing-angle-dependent light-curve breaks and blue spectral excess at breakout (Akashi et al., 2020). As a result, CSS161010 remains a benchmark object for distinguishing among compact-object engines, anisotropic outflows, and shock-powered versus accretion-powered FBOT channels.

Topic to Video (Beta)

No one has generated a video about this topic yet.

Whiteboard

No one has generated a whiteboard explanation for this topic yet.

Follow Topic

Get notified by email when new papers are published related to CSS161010.