---
title: 'AT 2020zso: Nuclear Tidal Disruption Event'
url: https://www.emergentmind.com/topics/at-2020zso
type: topic
---

# AT 2020zso: Nuclear Tidal Disruption Event

AT 2020zso is a nuclear tidal disruption event (TDE) in the galaxy SDSS J222217.13-071558.9 at \(z = 0.0563\), identified from its UV/optical blackbody evolution, broad He/H/Bowen emission-line phenomenology, and time-dependent spectral development. It is notable for two distinct observational results: optical spectroscopy indicates a newly formed, highly eccentric, nearly edge-on accretion disk, while late-time radio monitoring reveals at least two physically separate non-relativistic outflows. The event therefore connects prompt post-disruption disk formation, strong viewing-angle effects, and multi-episode outflow launching in a single system [2202.08268, 2509.14317].

## 1. Discovery, host environment, and TDE classification

AT 2020zso was discovered in 2020 by wide-field optical surveys as ZTF20acqoiyt, and was also seen by ATLAS and Gaia. Gaia astrometry places the transient at a nuclear offset of \(46 \pm 60\) pc, consistent with the galaxy nucleus. The host is an elliptical galaxy, and late-time narrow-line diagnostics classify the nucleus as a Seyfert AGN. That AGN context is central to the interpretation: the transient occurred in a nucleus already showing narrow-line evidence for activity, yet its broad-line and continuum evolution were argued to be inconsistent with ordinary AGN variability [2202.08268].

The TDE classification rests on the combination of a hot, blue continuum, evolving broad He/H/Bowen emission features, and UV/optical blackbody evolution unlike standard AGN fluctuations. The event is X-ray faint, with a combined Swift/XRT upper limit of \(L_X \lesssim 3.8\times10^{41}\ {\rm erg\ s^{-1}}\) in 0.3–10 keV for a 75 eV blackbody assumption, or \(\lesssim 4.3\times10^{41}\ {\rm erg\ s^{-1}}\) in 3–20 keV for an AGN-like power law. At the same time, it shows Bowen fluorescence lines, implying an ionizing EUV/soft X-ray source hidden from direct view. This observational combination places AT 2020zso naturally within the viewing-angle-dependent picture often invoked for optical/UV TDEs in which the intrinsic ionizing source is reprocessed and partially obscured [2202.08268].

## 2. UV/optical light curve and blackbody evolution

The UV/optical light curves show three phases in the ZTF \(g\) and \(r\) bands: a very steep early rise, a break to a slower rise, and a turnover into decline. The early rise was fit with
\[
L = a + b \times (t-t_0)^{\alpha},
\]
yielding an initial behavior consistent with \(L \propto t^2\), with \(\alpha = 1.9 \pm 0.4\). After a break around phase \(\sim -12\) d, the rise flattens to \(\alpha = 1.55 \pm 0.25\). The optical bands peak later than the UV bands, indicating strong color evolution [2202.08268].

Blackbody fits to the UV+optical spectral energy distributions show that the source cooled by about \(10{,}000\) K during the first part of the light curve. Around peak, the characteristic color temperature is \(\sim 25{,}000 \pm 5{,}000\) K. The blackbody radius initially expands approximately linearly, corresponding to a photospheric expansion speed of \(2900 \pm 300\ {\rm km\ s^{-1}}\), reaches a maximum near \(R_{\rm bb,max} \sim 10^{15}\ {\rm cm}\), plateaus near peak, and later recedes to \(R_{\rm bb} \sim 5\times10^{13}\ {\rm cm}\). For the favored black hole mass range, these scales correspond to \(\sim 3000\)–\(6000\,R_{\rm g}\) at peak and \(\sim 200\)–\(400\,R_{\rm g}\) at late times. The peak bolometric UV/optical luminosity is \(L_{\rm bol,peak} \simeq 7 \pm 1 \times 10^{43}\ {\rm erg\ s^{-1}}\), approximately Eddington for a black hole of \(\sim 5\times10^5\,M_\odot\), or about \(0.4\,L_{\rm Edd}\) for \(10^6\,M_\odot\) [2202.08268].

Several black hole mass estimators were considered. MOSFit gives \(\log_{10}(M_{\rm BH}/M_\odot) \approx 5.9\)–6.1, with a representative value \(\sim 6.0 \pm 0.3\), and prefers a low stellar mass with \(\beta \sim 0.63\), suggestive of a partial disruption, though with large systematic uncertainties. TDEMASS yields a mass consistent with \(\sim 10^6\,M_\odot\) and a stellar mass near solar. Independently, the host-galaxy stellar velocity dispersion is \(\sigma = 60 \pm 1\ {\rm km\ s^{-1}}\), corresponding to \(\log M_{\rm BH} = 5.3 \pm 0.4\) using McConnell & Ma and \(\log M_{\rm BH} = 6.2 \pm 0.3\) using Kormendy & Ho. A practical adopted range is therefore \(M_{\rm BH} \sim 5\times10^5 - 10^6\,M_\odot\). From the radius extrapolation, the first observations were taken about 15 days after expansion began, implying a disruption date around MJD 59149 and a rise time of about 35 days to peak; MOSFit instead gives MJD \(59157 \pm 2\) [2202.08268].

## 3. Emission-line evolution and Bowen deblending

The optical spectra show transient broad emission associated with He II \(\lambda4686\), H\(\alpha\), He I \(\lambda5876\), and Bowen N III at \(\lambda\lambda 4100, 4640\). Their evolution is highly structured. At the earliest epochs, \(-14\) and \(-13\) d, the broad features are better described by single broad Gaussians than by disk profiles, with line centroids consistent with rest velocity within uncertainties. In the \(-13\) d X-shooter spectrum, the measured Gaussian FWHM values are \(36{,}000 \pm 7000\ {\rm km\ s^{-1}}\) for H\(\alpha\) and \(31{,}000 \pm 1000\ {\rm km\ s^{-1}}\) for He II, with velocities near systemic: \(800 \pm 400\ {\rm km\ s^{-1}}\) for H\(\alpha\) and \(-170 \pm 135\ {\rm km\ s^{-1}}\) for He II. The He II/H\(\alpha\) equivalent-width ratio is initially very high, dropping from about 11 at \(-14\) d to 7 at \(-13\) d [2202.08268].

By \(-3\) d the line profiles are clearly non-Gaussian and become double-peaked. The complex near 4700 Å is interpreted mainly as He II \(\lambda4686\) with contamination from Bowen N III \(\lambda4640\), while the 6563 Å feature is identified as H\(\alpha\). The red peaks of He II and H\(\alpha\) occur at similar velocities, around \(+8000\ {\rm km\ s^{-1}}\), and the red wings extend to about \(26{,}000\ {\rm km\ s^{-1}}\). By \(+10\) d, H\(\alpha\) develops a triple-peaked structure consisting of a pronounced red peak, a blue peak, and a central component near systemic velocity; the red peak reaches around \(+11{,}500\ {\rm km\ s^{-1}}\). Similar structure appears in He I \(\lambda5876\) and in the feature around 4100 Å, which is identified not as H\(\delta\) but as Bowen N III \(\lambda4100\). The interpretation implies double-peaked Bowen lines, which were argued to be the first such profiles clearly seen in a TDE. Meanwhile, the He II/H\(\alpha\) equivalent-width ratio falls rapidly and then stabilizes near \(\sim 1.5\) after about \(-10\) d [2202.08268].

A key methodological issue is deblending He II from N III. A stable multi-Gaussian decomposition of the 4600–4700 Å complex was not obtained. Instead, the independently visible N III \(\lambda4100\) feature at \(+10\) d was used, assuming that N III 4100 and 4640 have roughly a 1:1 flux ratio and share the same velocity structure. Subtracting the observed N III 4100 profile from the He II complex yields an inferred He II profile that is remarkably similar to H\(\alpha\), supporting the conclusion that He II, H\(\alpha\), and Bowen N III arise from the same disk-like region [2202.08268].

## 4. Relativistic elliptical accretion-disk inference

The line profiles at \(-3\) and \(+10\) d were modeled with the relativistic elliptical disk prescription of Eracleous et al. (1995), applied to H\(\alpha\) and the deblended He II line. The earliest \(-13\) d spectrum was not fit with the disk model because the lines were interpreted as arising primarily in an outflowing envelope. The model uses seven main physical parameters: the emissivity power-law index \(q\), local line broadening \(\sigma\), disk inclination \(i\), eccentricity \(e\), orientation angle \(\phi_0\), inner pericentre radius \(r_1\), and outer pericentre radius \(r_2\). Broad priors were adopted: \(q=2\)–3, \(\sigma = 500\)–4500 km s\(^{-1}\), \(i = 0^\circ\)–\(90^\circ\), \(e=0\)–1, \(\phi_0 = 0^\circ\)–\(360^\circ\), \(r_1 = 100\)–550 \(R_g\), and \(r_2 = 750\)–4750 \(R_g\). A weak extra Gaussian component was also included to represent a possible outflow/wind or central component [2202.08268].

The best-fit solutions are mutually consistent in implying a highly inclined, highly eccentric, relatively compact accretion flow. For H\(\alpha\) at \(-3\) d, the fit gives \(q = 2.41 \pm 0.04\), \(\sigma = 1018 \pm 273\ {\rm km\ s^{-1}}\), \(i = 85 \pm 6^\circ\), \(e = 0.96 \pm 0.02\), \(\phi_0 = 256 \pm 2^\circ\), \(r_1 = 302 \pm 13\ R_g\), and \(r_2 = 4095 \pm 315\ R_g\). For He II at \(-3\) d, the fit gives \(q = 2.12 \pm 0.05\), \(\sigma = 1570 \pm 205\ {\rm km\ s^{-1}}\), \(i = 89 \pm 1^\circ\), \(e = 0.97 \pm 0.01\), \(\phi_0 = 219 \pm 2^\circ\), \(r_1 = 121 \pm 12\ R_g\), and \(r_2 = 2250 \pm 375\ R_g\). For H\(\alpha\) at \(+10\) d, the fit gives \(q = 2.28 \pm 0.04\), \(\sigma = 1270 \pm 145\ {\rm km\ s^{-1}}\), \(i = 78 \pm 5^\circ\), \(e = 0.96 \pm 0.02\), \(\phi_0 = 242 \pm 2^\circ\), \(r_1 = 270 \pm 11\ R_g\), and \(r_2 = 1453 \pm 74\ R_g\). For He II at \(+10\) d, the fit gives \(q = 2.92 \pm 0.06\), \(\sigma = 1330 \pm 100\ {\rm km\ s^{-1}}\), \(i = 86 \pm 3^\circ\), \(e = 0.97 \pm 0.01\), \(\phi_0 = 208 \pm 1^\circ\), \(r_1 = 530 \pm 10\ R_g\), and \(r_2 = 2100 \pm 90\ R_g\) [2202.08268].

Averaging across lines and epochs yields \(i \simeq 85 \pm 5^\circ\), \(e \simeq 0.97 \pm 0.01\), \(\phi_0 \sim 210^\circ\)–\(260^\circ\) with an average near \(240^\circ\), an inner radius of several hundred \(R_g\), an outer radius of several thousand \(R_g\), local broadening of \(\sim 1000\)–\(1500\ {\rm km\ s^{-1}}\), and emissivity slope \(q \sim 2\)–3. The line asymmetry is essential to the eccentricity inference: in a standard relativistic circular disk, the blue peak is usually at least as strong as the red peak because of Doppler boosting, whereas in AT 2020zso the red peak is stronger. That strongly favors a non-axisymmetric eccentric disk over a circular one [2202.08268].

## 5. Envelope, viewing angle, and spin constraints

The preferred physical picture is a two-component structure consisting of an optically thick outflowing/reprocessing envelope at early times and an accretion disk that becomes spectroscopically visible near peak. Before peak, the photosphere expands and the continuum plus lines are dominated by a broad, quasi-spherical reprocessing region. The very broad, nearly Gaussian He II seen at \(-14\) to \(-13\) d is attributed to this envelope rather than to the disk. The rapid drop in the He II/H\(\alpha\) ratio is interpreted as an optical-depth effect: when the envelope is compact and dense, H\(\alpha\) is more self-absorbed than He II, making the spectrum appear He-rich; as the envelope expands and cools, H\(\alpha\) becomes less suppressed [2202.08268].

Around peak, the envelope reaches \(R_{\rm env} \sim 10^{15}\ {\rm cm}\), corresponding to \(\sim 3000\)–\(6000\,R_g\) for \(M_{\rm BH} \sim 5\times10^5 - 10^6\,M_\odot\), and then becomes optically thin. At that point, disk emission dominates the broad-line profiles. The late-time blackbody radius of \(\sim 5\times10^{13}\ {\rm cm}\) overlaps the spectroscopically inferred disk scale, suggesting that the UV continuum may become directly disk-dominated at late epochs. The inferred eccentricity is also close to the characteristic value for returning tidal debris,
\[
e \sim 0.98 \left[\frac{(M_{\rm BH}/M_\star)}{10^6}\right]^{1/3},
\]
which supports the interpretation that the debris has not yet circularized strongly. This suggests that an accretion disk can form within about a month after disruption while remaining highly eccentric [2202.08268].

AT 2020zso also bears directly on TDE unification arguments. It is Bowen-strong and X-ray dim, yet line-profile modeling gives \(i \sim 85^\circ\). The event was therefore presented as the first quantitative confirmation of a high-inclination geometry in a Bowen TDE, consistent with models in which viewing angle largely determines whether the EUV/soft X-ray source is directly observed. Alternative explanations were examined and disfavored. A binary SMBH interpretation was rejected because binary-driven eccentric AGN disk scenarios evolve on thousands-of-years timescales, whereas AT 2020zso shows major line-profile changes on weeks timescales; the narrow-line region also appears more consistent with rotation-dominated gas than with a disturbed dual-AGN system. Turn-on or changing-look AGN, bipolar outflows, and spiral-arm interpretations were also considered less satisfactory than a fresh TDE disk [2202.08268].

A further inference concerns black hole spin. The modeled disk parameters, particularly inclination and orientation, do not change significantly over a baseline of about 15 days, so the empirical constraint \(t_{\rm align} > 15\ {\rm days}\) was adopted. The only explicit timescale formula quoted in this context is
\[
t_{\rm dyn} = 2 M_6 \xi^{3/2}\ {\rm days},
\]
which gives roughly 5–15 days for AT 2020zso. Using published alignment and precession calculations, the absence of measurable precession-driven inclination change was taken to rule out high spins, with a quoted upper limit \(a < 0.8\) for \(M_{\rm BH}\sim10^6\,M_\odot\); in the abstract this is phrased as excluding high spin values for disk viscosity \(\alpha \gtrsim 0.1\). This is an indirect, theory-informed constraint rather than a direct measurement of precession frequency [2202.08268].

## 6. Long-lived radio emission and multiple outflows

Multi-year radio monitoring later showed that AT 2020zso also produced unusually rich long-lived radio emission. The campaign extended for more than three years, primarily with the VLA and supplemented by GMRT. VLA coverage spans \(L/S/C/X/Ku\) bands, roughly 1.3–17.4 GHz, while GMRT adds 0.65 and 1.26 GHz. The first VLA detection occurred at \(\delta t=40\) d with \(F_\nu(15\,{\rm GHz})\approx 22\,\mu{\rm Jy}\). The radio luminosity peaks at \(\sim 2\times 10^{38}\) erg s\(^{-1}\), well below the \(\gtrsim 10^{40}\) erg s\(^{-1}\) characteristic of on-axis relativistic TDE jets. High-frequency emission rises and fades earlier than low-frequency emission, as expected for an expanding self-absorbed synchrotron source whose peak frequency moves downward with time [2509.14317].

The radio light curves require two flares. The first rises soon after optical peak, reaches maximum around \(\delta t\sim 1\) yr, and then fades. A second flare begins around \(\sim 800\) d and becomes dominant after \(\sim 900\) d. Broken-power-law fits give rise and decay indices \(a_1=2.0\pm0.2\) and \(b_1=2.5\pm0.2\) for the first flare, and \(a_2=5.6\pm1.8\) and \(b_2=2.0\pm0.5\) for the second. The decisive evidence for two physical components is spectral: at \(\delta t=882\) d and especially \(\delta t=926\) d, the radio spectral energy distribution shows two distinct peaks. At \(\delta t=882\) d, the first component has \(F_{\nu_{a,1}} = 0.95\pm0.16~{\rm mJy}\) and \(\nu_{a,1}=0.89\pm0.21~{\rm GHz}\), while the second has \(F_{\nu_{a,2}} = 0.27\pm0.05~{\rm mJy}\) and \(\nu_{a,2}=4.95\pm0.73~{\rm GHz}\). At \(\delta t=926\) d, the first is constrained by \(F_{\nu_{a,1}}>0.67~{\rm mJy}\) and \(\nu_{a,1}<0.81~{\rm GHz}\), while the second has \(F_{\nu_{a,2}} = 0.25\pm0.04~{\rm mJy}\) and \(\nu_{a,2}=3.95\pm0.42~{\rm GHz}\). The \(\delta t=926\) d spectrum was stated to be incompatible with any single-region synchrotron-break ordering of the Granot & Sari type, so the emission was interpreted as the sum of two independent synchrotron sources rather than a refreshed shock or simple energy injection into one blast wave [2509.14317].

The radio analysis adopts a non-relativistic synchrotron-emitting outflow with
\[
N(\gamma)\propto \gamma^{-p}, \qquad \gamma>\gamma_m,
\]
and spectral ordering \(\nu_m < \nu_a < \nu_c\). MCMC fitting yields \(p_1 = 2.92\pm0.08\) and \(p_2=3.22\pm0.07\). Physical parameters are then inferred using the equipartition/minimum-energy formalism of Barniol Duran et al. for a non-relativistic spherical blast wave, with an emitting shell of thickness \(0.1R\), filling factors \(f_A=1\) and \(f_V\simeq0.36\), and microphysical assumptions \(\epsilon_e=\epsilon_B=0.1\) and \(\gamma_m=2\) [2509.14317].

For the first flare, the inferred launch date is \(t_{0,R_1}=74^{+12}_{-12}\ {\rm d}\) relative to optical discovery, about 56 d after the bolometric peak. The outflow remains non-relativistic, with \(v\simeq 0.11c\) and epoch-by-epoch \(\beta\) values around 0.10–0.12. The radius evolves from \(R \sim 10^{16.04}\ {\rm cm}\) at \(\delta t=101\) d to \(R \sim 10^{17.42}\ {\rm cm}\) at \(\delta t=882\) d, with \(R>10^{17.39}\) cm at \(\delta t=926\) d. The internal energy rises and then plateaus around \(E \sim 2.2\times10^{49}\ {\rm erg}\), more specifically \(\log_{10}E\simeq 49.38\)–49.53 at late epochs. The magnetic field declines from \(B\sim 10^{0.27}\ {\rm G}\) to \(B\sim 10^{-0.81}\ {\rm G}\), and the ambient density from \(n_{\rm ext}\sim 10^{3.13}\ {\rm cm^{-3}}\) to \(n_{\rm ext}\sim 10^{0.98}\ {\rm cm^{-3}}\). The kinetic mass estimate is \(M_k\sim 2\times 10^{-3}M_\odot\) [2509.14317].

For the second flare, the inferred launch date is \(t_{0,R_2}=836^{+10}_{-13}\ {\rm d}\) relative to optical discovery. This second outflow is likewise non-relativistic and approximately freely coasting, with characteristic velocity \(v\approx0.16c\) and epoch-by-epoch \(\beta\approx0.13\)–0.23. Its radius grows from \(R\sim 10^{16.44}\ {\rm cm}\) at \(\delta t=882\) d to \(R\sim 10^{17.33}\ {\rm cm}\) at \(\delta t=1275\) d, with \(R>10^{17.41}\) cm at \(\delta t=1382\) d. The energy rises to \(E\approx 3.8\times10^{49}\ {\rm erg}\), with \(\log_{10}E\) increasing from 48.34 at \(\delta t=882\) d to 49.65 at \(\delta t=1275\) d and \(>49.75\) at \(\delta t=1382\) d. The magnetic field falls from \(B\sim 10^{0.09}\ {\rm G}\) to \(B\sim 10^{-0.61}\ {\rm G}\), and the ambient density from \(n_{\rm ext}\sim 10^{2.73}\ {\rm cm^{-3}}\) to \(n_{\rm ext}\sim 10^{1.33}\ {\rm cm^{-3}}\). Its kinetic mass is also of order \(M_k\sim 2\times10^{-3}\,M_\odot\) [2509.14317].

The external environment inferred from the radio data is consistent with a Bondi-like circumnuclear medium, \(n_{\rm ext}(r)\propto r^{-k}\), with \(k \simeq 1.53\) for the first outflow and \(k \simeq 1.51\) for the second, close to the \(r^{-3/2}\) expectation for spherical Bondi accretion. Using \(M_{\rm BH}\approx10^6\,M_\odot\), \(T=10^7\) K, \(\gamma=5/3\), and \(\mu=0.6\), the Bondi radius is estimated as \(R_B \approx 1.2\times10^{17}\ {\rm cm}\). No density break near \(R_B\) is required to explain the second flare, although magnetic-field inhomogeneities were suggested by an unusually shallow optically thick slope \(\alpha\simeq1.05\) at \(\delta t\approx365\) d, compared with the canonical self-absorbed synchrotron value \(\alpha=5/2\) [2509.14317].

The preferred physical interpretation links the two radio outflows to two accretion states. The first is favored to be accretion-driven, launched while the source was still accreting at a relatively high Eddington fraction; the optical luminosity implies \(L/L_{\rm Edd}\sim 0.5\), consistent with a slim disk capable of driving winds. The second is associated with a later transition to an advection-dominated accretion flow, after an intermediate regime in which strong outflow production is not expected below \(\dot M \lesssim 0.3 \dot M_{\rm Edd}\), and then may resume when \(\dot M \lesssim 0.03 \dot M_{\rm Edd}\). Off-axis relativistic jet and delayed-accretion explanations were considered less compelling, and “energy injection into Outflow 1” was explicitly rejected because two simultaneous spectral peaks isolate two emitting regions [2509.14317].

## 7. Uncertainties, alternatives, and significance

Several limitations qualify the interpretation. On the optical side, the MOSFit light-curve fit is formally poor because the event shows rapid temperature evolution and a broken rise that the model does not capture well, so disruption parameters and stellar-mass estimates carry substantial systematics. The He II profile is blended with Bowen N III, and the deblending assumes a 1:1 flux ratio between N III 4100 and 4640 and identical velocity structure, which likely introduces additional systematics; for this reason the H\(\alpha\) fits were treated as primary. The disk model is idealized, not a full radiative-transfer or hydrodynamical calculation, and some parameters—especially the inner and outer radii—remain degenerate across lines and epochs. The spin constraint is indirect, relying on comparison to published precession and alignment theory rather than on a direct measurement of precession [2202.08268].

The radio interpretation likewise depends on modeling assumptions. The inferred \(R\), \(E\), \(B\), and \(n_{\rm ext}\) values rely on equipartition and quasi-spherical shell geometry, and the launch dates depend on the assumption of free expansion; if the outflows decelerate, the inferred \(t_{0,R}\) values would shift. The first outflow’s launch time disfavors unbound debris and pure collisionally induced outflow scenarios, but does not absolutely exclude some contribution from such processes. The similar density profiles inferred for the two outflows are somewhat surprising if both were quasi-spherical, because the second would be expected to encounter gas modified by the first; this suggests, plausibly, that the two ejections may have different geometries, such as one more polar and one more toroidal. A repeated partial-TDE interpretation was also briefly considered, since optical modeling suggested a partial disruption and the separation between the two radio outflows is comparable to recurrence times in repeating partial TDE candidates, but the lack of an observed second optical flare—during a sun-constrained interval—left that possibility unresolved [2509.14317].

Within these caveats, AT 2020zso is significant for several reasons. Optical spectroscopy provides unusually direct evidence for a compact, edge-on, highly eccentric TDE accretion flow, closer to the expected geometry of newly returned debris than the near-circular disks inferred in some earlier events. The event also supports the unification picture in which Bowen-strong, X-ray-faint TDEs are preferentially seen at high inclination. Finally, the radio data show that a thermal TDE can launch at least two distinct non-relativistic outflows over a multi-year interval, with the separation of components demonstrated spectrally rather than inferred solely from late-time rebrightening. AT 2020zso therefore serves as an important observational link between TDE disk formation, incomplete circularization, anisotropic reprocessing, and episodic accretion-driven outflows in an AGN host [2202.08268, 2509.14317].

Source: https://www.emergentmind.com/topics/at-2020zso