Next Article in Journal / Special Issue
Progress on the Haystack Observatory Postprocessing System
Previous Article in Journal
Introduction to Electroweak Baryogenesis
Previous Article in Special Issue
Measuring Photon Rings with the ngEHT
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling Reconstructed Images of Jets Launched by SANE Super-Eddington Accretion Flows around SMBHs with the ngEHT

1
Black Hole Initiative, Harvard University, 20 Garden Street, Cambridge, MA 02138, USA
2
Center for Astrophysics, Harvard & Smithsonian, 60 Garden Street, Cambridge, MA 02138, USA
3
Department of Physics & Astronomy, The University of Texas at San Antonio, One UTSA Circle, San Antonio, TX 78249, USA
*
Author to whom correspondence should be addressed.
Galaxies 2022, 10(6), 117; https://doi.org/10.3390/galaxies10060117
Submission received: 14 November 2022 / Revised: 6 December 2022 / Accepted: 7 December 2022 / Published: 13 December 2022

Abstract

:
Tidal disruption events (TDEs) around supermassive black holes (SMBHs) are a potential laboratory to study super-Eddington accretion disks and sometimes result in powerful jets or outflows which may shine in the radio and sub-millimeter bands. In this work, we modeled the thermal synchrotron emission of jets by general relativistic radiation magneto-hydrodynamics (GRRMHD) simulations of a BH accretion disk/jet system which assumed the TDE resulted in a magnetized accretion disk around a BH accreting at ∼12–25 times the Eddington accretion rate. Through synthetic observations with the Next Generation Event Horizon Telescope (ngEHT) and an image reconstruction analysis, we demonstrate that TDE jets may provide compelling targets within the context of the models explored in this work. In particular, we found that jets launched by a SANE super-Eddington disk around a spin a * = 0.9 reach the ngEHT detection threshold at large distances (up to 100 Mpc in this work). A two-temperature plasma in the jet or weaker jets, such as a spin a * = 0 model, requires a much closer distance, as we demonstrate detection at 10 Mpc for limiting cases of a * = 0 , R = 1 or a * = 0.9 , R = 20 . We also demonstrate that TDE jets may appear as superluminal sources if the BH is rapidly rotating and the jet is viewed nearly face on.

1. Introduction

Tidal disruptions of stars by supermassive black holes (SMBHs), or tidal disruption events (TDEs), have recently become a regularly observed transient phenomenon. Stars which enter the tidal radius
R t 7 × 10 12 M BH 10 6 M 1 / 3 M * M 1 / 3 R * R [ cm ]
of the central SMBH in their host galaxy will be disrupted [1,2], either partially or fully, depending on the orbit and equation of state of the star [3,4]. The bound stream of gas returns towards the BH delivering mass at the fallback rate ( M ˙ fb ). Apsidal precession of the returning stream leads to self-intersection with material that has yet to pass through pericenter and leads to dissipation and disk formation. Dissipation due to self-intersection may also be a source of early emission in a TDE [5]. After the initial rise to peak, the fallback rate follows a power law behavior which can be approximated as
M ˙ fb = M ˙ fb , peak t t fb 5 / 3 .
Here, M ˙ fb , peak is the peak mass fallback rate:
M ˙ fb , peak M ˙ Edd 133 M BH 10 6 M 3 / 2 M * M 2 R * R 3 / 2
making the “frozen in” approximation as in Stone et al. [6], and
t fb = 3.5 × 10 6 s M BH 10 6 M 1 / 2 M * M 1 R * R 3 / 2 ,
is the fallback time, which is the orbital time of the most bound part of the stream. Of note is the fact that the mass fallback rate can greatly exceed the Eddington mass accretion rate M ˙ Edd . Since the mass accretion rate is expected to be similar in magnitude to the mass fallback rate, it is possible that TDEs result in a super-Eddington accretion disk, which are geometrically thick, radiatively inefficient accretion disks [7,8]. The exact power law behavior varies with the properties and orbit of the star [9].
TDEs are typically seen as optical/X-ray transients [10,11], but several TDEs have resulted in outflows or jets which shine in the radio bands [12]. In the most common case in which no relativistic jet is launched (commonly referred to as “non-jetted” TDEs), the X-ray and optical/UV luminosity follows a roughly t 5 / 3 decline, similar to the fallback rate. If the TDE leads to prompt disk formation, the X-rays are thought to arise from an accretion disk, and the optical/UV emission arises from a large scale reprocessing layer [13,14]. TDEs have also been observed to launch relativistic X-ray jets in a few cases. These jetted TDEs have been argued to arise due to a magnetically arrested disk (MAD) [15] forming around the BH during the TDE [13,16,17], which leads to jet production via the Blandford–Znajek (BZ) mechanism [18] extracting spin energy from the BH. Alternatively, these powerful relativistic jets may be produced thanks to radiative acceleration of gas through a narrow funnel region [19,20].
It is important to note that the BZ mechanism is the strong field limit of the magnetic Penrose process (MPP), but the MPP can operate with weaker magnetic field strengths [21,22]. Note as well that the magnetic Penrose process, as its name suggests, is itself the magnetic flavor of the Penrose process (PP) [23]. While the PP relies on the existence of negative energy orbits in the vicinity of the BH, the MPP relies on negative energy orbits and quadrupole electric fields produced by twisting of the magnetic field lines threading the BH horizon. A key distinction, and why the BZ mechanism is often cited while interpreting astrophysical jets, is that the energy extraction from the MPP can exceed 100 percent, whereas the maximum energy extraction from a PP is ∼20% [21,22].
A handful of non-jetted TDEs have been observed to produce radiowaves peaking at tens of GHz with L radio 10 37 39 erg s 1 . This emission is thought to arise from an outflow launched by the TDE with velocity v 0.1 c shocking on the gas surrounding the BH. Meanwhile, jetted TDEs produce bright radio emission peaking at L radio 10 40 42 erg s 1 . The appearance of radio emission is often delayed by several weeks from the initial appearance of the optical/UV/X-ray emission in non-jetted TDEs, which hints at some connection to the disk formation process to the occurrence of outflows.
When M ˙ fb rises to peak, it has previously been assumed that by this stage a circularized accretion disk has formed [13,17]. The first direct demonstration of circularization near the peak fallback rate was recently demonstrated in a numerical simulation by Steinberg and Stone [5]. Multiple authors have argued in favor of a picture in which an inner accretion flow is surrounded by a quasi-spherical reprocessing layer, since this naturally explains the sometimes delayed appearance of X-ray emission in optical/UV-discovered TDEs [13,14,24]. This picture naturally arises if the accretion flow is actually super-Eddington [13], which has motivated multiple studies of GRRMHD simulations magnetized of super-Eddington disks. Motivated by this fact and the demonstration that a “standard and normal evolution” (SANE) [15] super-Eddington accretion disk also leads to viewing angle effects which may explain the behavior in non-jetted TDEs [17], we conducted a study of the outflows launched by SANE models in Curd et al. [14] and studied their radio-submm emission. We focus on these SANE models in this work as well. The BZ mechanism also acts in the SANE model [14,25]; however, the disk and BH have a much lower magnetic flux and thus a substantially lower jet power. Although BZ driven outflows in SANE disks can be relativistic, they are weakly comptonized and do not appear as powerful, beamed X-ray sources like the presently known jetted TDEs [14].
The Next Generation Event Horizon Telescope (ngEHT) will provide more baseline coverage and faster response times than the previous mission. An estimate for ngEHT is to double the antenna sites [26] of its 20 μ as predecessor EHT, and thus the number of possible baseline pairs and triads of sites available for imaging jet/accretion flow/black hole systems will scale combinatorially (the number of baselines grows with the number of antennae as (N(N − 1)/2). This could allow for interesting sources, such as jets from nearby tidal disruption events, to be imaged directly. In our previous work [14], we provided the first demonstration that SANE super-Eddington accretion flows can produce radio emission which is bright enough at 230 GHz to be detected and resolved. Here we take things a step further and produce reconstructed images assuming such jets happen in the nearby universe.
The detection rate of TDEs in the optical/UV/X-ray is expected to grow rapidly once the Large Synoptic Survey Telescope comes online [27,28]. Assuming rapid follow-up of TDEs in radio-submm bands finds detectable emission, this could provide a large number of targets for the ngEHT. As we demonstrated in Curd et al. [14], some models produced detectable emission, even at ∼180 Mpc. At this distance, a conservative estimate of the volume integrated TDE rate suggests more than 200 TDEs per year assuming volumetric TDE rates based on Stone and Metzger [29]. Even at < 40 Mpc, we estimate that several TDEs should occur per year (see Figure 2 in Curd et al. [14]), which suggests some nearby TDEs may become targets of opportunity during the ngEHT mission.
We stress that jets such as those in our first work on the subject of jets from SANE models of TDE accretion disks [14] do not resemble any previously detected radio TDEs. This may suggest most, or even all, TDEs do not form accretion disks which resemble SANE models to begin with. However, the number of TDEs that have appeared in the radio-submm in the first place is extremely small as of this writing; fewer than twenty radio TDEs have been reported. Furthermore, magnetic fields were certainly present in the forming disk, albeit dynamically subdominant to hydrodynamic effects early in the disk formation [30,31]. Nevertheless, it is possible that after the disk circularizes, which Steinberg and Stone [5] suggests may take tens of days, the magnetic field builds up in a dynamo effect similar to Sadowski et al. [30]. In this case, one would almost certainly expect the magnetic field to become dynamically important, in which case a magnetized outflow may be launched as in Curd et al. [14]. In addition, if the currently known jetted TDEs are indeed MAD disks around rapidly rotating BHs, it stands to reason that a sub-class of less strongly magnetized disks which could be described by a SANE model should exist. TDEs continue to surprise observers in terms of the range of behavior, so such a jet formation channel may yet be discovered.

2. Numerical Methods

2.1. GRRHMD Simulations

Throughout this work, we often use gravitational units to define length and time. In particular, we use the gravitational radius
r g = G M BH c 2
and the gravitational time
t g = G M BH c 3 ,
where M BH is the mass of the black hole (BH). We also adopt the following definition for the Eddington mass accretion rate:
M ˙ Edd = L Edd η NT c 2 ,
where L Edd = 1.25 × 10 38 ( M BH / M ) erg s 1 is the Eddington luminosity and η NT is the radiative efficiency of a Novikov–Thorne thin disk around a BH with spin parameter a * [32].
We conducted an imaging analysis of GRRMHD simulations presented in Curd et al. [14]. In particular, we analyzed the most massive BH models m7a0.0-HR and m7a0.9-HR, which are M BH = 10 7 M BHs of spin a * = 0 and 0.9 BHs. We specify the simulation diagnostics relevant for this work in Table 1. The simulations were conducted in 2D ( r , ϑ ) coordinates on a N r × N ϑ = 640 × 256 grid with added resolution near the poles to adequately resolve both the disk and jet. The radial grid cells were logarithmically spaced with a maximum domain radius of R max = 10 5 r g to capture the large-scale features of the jet.
On horizon scales, gas is flowing across the BH horizon in an accretion disk due to angular momentum transport driven by the magneto-rotational instability. The disk is optically thick and turbulent, and gas inside of the disk is advected with the gas across the BH horizon. However, an optically thin funnel above and below the disks exists. Here, radiation can escape freely and pushes on gas, accelerating a significant outflow. In addition, the funnel is magnetized and sometimes exhibits magnetization parameter σ = b 2 / ρ c 2 > 1 , where b is the magnetic field strength and ρ is the mass density of the plasma. In the jet, magnetic energy is partially converted into kinetic energy, as it contributes to accelerating gas into an outflow.
In both simulations, radiative and Poynting acceleration drove fast outflows. The jet reached relativistic speeds with Lorentz factor γ > 5 for model m7a0.9-HR, which was likely due to the BZ effect extracting spin energy from the BH, which may produce roughly ∼1 percent of the jet efficiency, even though the magnetic flux threading the black hole was well below the MAD limit [14,25]. The primary sites of dissipation were the jet head and internal shocks inside of the jet. Internal shocks were due to fast- and slow-moving gas interacting downstream of the jet head, in addition to recollimation shocks. As we show in Figure 1, this resulted in a hot, magnetized jet which reached large scales ( r > 30 , 000 r g ) by the end of the simulation. The a * = 0.9 model had a significantly more magnetized jet and also produced a more powerful jet by a significant fraction (see Table 1).
The simulations were conducted in 2D to overcome the substantial numerical requirements of resolving the accretion disk and jet in a global GRRMHD simulation. While SANE super-Eddington disks can be described by an axisymmetric flow quite well [33], the jet may undergo 3D instabilities which cannot be captured in our numerical setup. In particular, the jet may become kink-unstable as it propagates [34]. This will certainly change the morphological features, and additional dissipation of magnetic energy along the jet may occur. If the dissipation is extreme, especially at large scales such as those we consider here, the synchrotron emission which we describe in the next section may differ from our analysis of 2D jets.

2.2. 230 GHz Emission

We post-processed the KORAL simulation data with the general relativistic ray tracing (GRRT) code ipole [35,36,37], which includes synchrotron and Bremsstrahlung emission and absorption. Since the KORAL simulation was conducted in 2D, we first copied the 2D data onto a full 3D ( r , ϑ , ϕ ) grid with 32 cells in azimuth which spanned ϕ = 0 2 π before ray tracing. We did this by simply copying the original 2D data at each ϕ cell such that the new 3D data maintained the assumed axisymmetry in ϕ . The electron distribution function was assumed to be thermal. Ohmura et al. [38,39] demonstrated that large-scale active galactic nuclei (AGN) jets can produce a two-temperature plasma. Motivated by their findings and the possibility that a two-temperature plasma will be produced due to shocks at the jet head and within the jet itself, we tested a simple two-temperature jet model by scaling the electron temperature relative to the ion temperature via the plasma temperature ratio:
R = T i T e ,
where T i is the temperature of the ions and T e is the temperature of the electrons. Note that T i was obtained directly from the KORAL simulation by setting T i = T gas .
The peak of the radio-submm spectra in m7a0.0-HR is lower than that of m7a0.9-HR, so increasing R has a much more significant impact on the 230 GHz emission and can make the jet undetectable even at 10 Mpc for values of R > 2 [14]. It is possible that a non-thermal electron distribution will have greater high-energy emission even as R increases, but we saved an exploration of non-thermal electron models for a future analysis.
Each model was imaged at 230 GHz. For both models, we imaged the simulation at times t = 38 , 000 t g and t = 78 , 000 t g for a difference in observing times of ∼23 days. We chose a distance D = 10 Mpc, R = 1 , observing angles relative to the jet axis (z in Figure 1) of θ = 10 , 45 , and 90 , respectively. Note that we use θ for the observer angle, and ϑ is the polar angle in the KORAL grid coordinates. For model m7a0.9-HR, we also tested limiting cases D = 100 Mpc, R = 1 , and D = 10 Mpc, R = 20 , imaged at θ = 90 . The total 230 GHz flux of each ray-traced model is tabulated in Table 2.
We show a full library of each of the ipole images convolved with a Gaussian beam with a 20 μ as full width at half maximum (FWHM) in Figure A1 and Figure A3.

2.3. Synthetic ngEHT Observations and Image Reconstruction

In order to test to what extent the jet features in our models can be observed, we simulated observations with a potential ngEHT array, consisting of the 2022 EHT stations and 11 additional stations, selected from Raymond et al. [40] and similar to the ngEHT reference array used in the ngEHT Analysis Challenges [41]. The new dishes were assumed to have a diameter of 10 m and a receiver temperature of 50 K, with the array operating at a bandwidth of 8 GHz. For each image, we simulated a 24-hour observation with a 50% duty cycle with this array, using the ngehtsim 1 library, which makes use of eht-imaging [42,43]; see also [44]. The atmospheric opacity was set to reflect a good day in April, using the top 1 σ quantile from the MERRA-2 data interpolated and integrated for each site on a 3 h cadence for a 10-year period [45,46]. Thermal noise was added to the complex visibilities, visibility phases were randomized, and no systematic visibility amplitude errors were added to the data.
We subsequently used the regularized maximum likelihood framework in eht-imaging to produce image reconstructions, with maximum entropy and (squared) total variation regularizers, fitting to visibility amplitudes and closure phases; see [42,43,47]. After establishing a set of well-performing imaging parameters on the m7a0.9-HR, R = 1 model at t g = 78 , 000 and a distance of 10 Mpc (Figure 2), we applied the same script to all other simulated datasets.

3. Results

In this section, we comment on the detectability of our models and then compare the ipole images with the reconstructed images. We comment on features which may be of interest in terms of the broader study of astrophysical jets.

3.1. Reconstructed Images

In the ray-traced image (i.e., see the left panel in Figure 2), the jet head produces bright emission as it shocks on the circumnuclear medium (CNM). In addition, various shocks occur within the jet due to both slow- and fast-moving components colliding radially and due to recollimation shocks. This leads to dissipation within the jet and bright “bubbles” of emission at 230 GHz. As we show in the right panel of Figure 2, the jet head and the structures in the jet are faithfully reproduced in the reconstruction for favorable viewing angles ( θ = 45 and 90 ) as long as the source is nearby ( D = 10 Mpc). Jets viewed near θ = 10 are dominated by emission from the jet head, and distinguishing internal jet features would be unlikely. This can be seen by comparing the base images with the full library of reconstructed images for models m7a0.0-HR and m7a0.9-HR (Figure A1Figure A4).
For distant sources ( D = 100 Mpc), distinguishing internal features is impossible, and only the jet head can be fully distinguished in the reconstruction (right panel in Figure 3). This would still allow for the jet motion to be tracked, but detailed information is lost.

3.2. Tracking Jet Motion

In this subsection, we demonstrate that the original ray-traced images and the reconstructed images allow for the jet motion to be tracked and yield similar results for the time evolution of the jet. The jet features are approximately Lorentzian, so we fit Lorentzian profiles to the image to find the position of the top and bottom jet in both the base images and the reconstructed images. We detail the peak finding algorithm in Appendix B. Since we cannot properly center the jet in the reconstructed images (there is no bright, central radiation from the near BH), we only measured the distance between the two jet peaks y 1 and y 2 , respectively. We define the apparent jet length as:
l jet = | y 2 y 1 | .
Note that we have not differentiated the “top” or “bottom” jet here, as we are only concerned with the total distance between the jet heads. We obtained errors on the jet length from the error estimates of the jet head locations using standard error propagation analysis:
δ l jet = δ y 1 2 + δ y 2 2 .
We computed the relative difference between the jet lengths for the ray-traced ( l jet ) and reconstructed ( l jet , rec ) images in order to quantify the extent to which measurements of the jet length and velocity agree. We found that, generally, the ray-traced and reconstructed images yield similar jet lengths within 3 standard deviations (Figure 4). In general, there are much larger errors on the fit for the Lorentzian profile’s center at steep angles (i.e., see the relative difference for θ = 10 ). The agreement is also effected by how bright the source is, as illustrated by the shift to the right for model m7a0.9-HR when D = 100 Mpc. We compiled estimates of the jet lengths for each image in Table 3.
Since we lack centering information in the reconstructed images, we chose to estimate the jet velocity perpendicular to the line of sight by assuming both the top and bottom jet have the same speed. Then, the velocity of the jet in the source’s frame is:
v = 1 2 | l jet ( t 2 ) l jet ( t 1 ) | t 2 t 1 .
Note that we used the same expression to derive the velocity in the reconstructed images ( v rec ) but replaced l jet with l jet , rec in Equation (11). Very long baseline interferometry (VLBI) observations can provide the apparent motion of the jet, which may be superluminal due to relativistic effects. To account for this, we assumed the jet has a velocity of v rec and then estimated the apparent velocity ( v app , rec ) via the time in the observer’s frame ( t 1 , t 2 ):
v app , rec = 1 2 | l jet , rec ( t 2 ) l jet , rec ( t 1 ) | t 2 t 1 = v rec 1 v rec c cos θ sin θ 1 .
We used the relationship
t 2 t 1 = ( t 2 t 1 ) 1 v rec c cos θ sin θ
in the last expression. Note that the division by sin θ was to account for the fact the v rec measures the velocity parallel to the line of sight with no time delay effects while we require an estimate of the velocity along the jet axis (which we can obtain since the geometry is fully known). We only present the apparent velocity for the reconstructed images ( v app , rec ), since this represents an estimate of what VLBI observations would truly be.
We used the data at t 1 = 38 , 000 t g and t 2 = 78 , 000 t g for each model to estimate the jet velocity. We tabulated the estimated velocity and apparent velocity for each model in Table 4. We found excellent agreement between the velocities derived from the ray-traced and reconstructed images. It is interesting to note the apparently faster jet for R = 20 . We suspect the increased speed is due to the emitting material being dominated by material near the jet axis rather than some of the slower moving material around the jet head, which does not produce much emission at 230 GHz as the electron temperature is reduced. The more powerful jet model m7a0.9-HR demonstrates that such jets may appear as superluminal sources, as we found a maximum v app , rec 1.526 c at θ = 10 .

4. Discussion

4.1. Extracting Jet Physics from VLBI Images

A key feature of the jets in our models is the bright “bubbles” (or knots) of 230 GHz emission, which appear to correlate with recollimation shocks. Such structures have been seen in VLBI images of various AGN jets [48,49,50]. Previous simulations of jets in various astrophysical contexts have demonstrated that recollimation occurs when there is a pressure mismatch between the jet and the surrounding medium, which could be a static atmosphere or a slower-moving jet sheath [51,52,53,54]. The number of recollimation shocks along the jet axis is dependent on the properties of the jet and medium. It is therefore possible that direct VLBI of TDE jets will allow in-depth modeling of jet launching and could also aid in constraining the properties of the surrounding medium. For instance, one work successfully applied simulations of MHD jets to constrain properties of BL Lacartae, which is a blazar jet with recollimation features [55].
We suggest that a similar approach may be applied in TDE jets. With a suitable exploration of the parameter space, it is conceivable that an analysis similar to that of [55] could be applied to TDE jets in cases where VLBI is possible. A broader exploration of TDE jets through various simulation methodologies is strongly suggested. We plan to explore the effects of the ambient medium, magnetic field strength, and disk accretion rate on the jet properties in the case of a SANE, super-Eddington disk in a future work.

4.2. Proposed Observational Methodology

Our synthetic ngEHT observations demonstrate that a 24-hour observation may be sufficient to study both the structure and/or motion of newly born TDE jets. This is much shorter than the fallback time, which is on the order of month(s), and the duration of radio emissions in several TDEs, which can sometimes be visible for years [12]. Our suggested observational strategy is conducting rapid followup of newly discovered optical/X-ray TDEs when they are near the peak of their emission in order to study both the early- and late-time properties of their jets (if present). Observations with a single telescope at 230 GHz could be conducted to search for TDEs emitting in the radio. If emission at 230 GHz is detected, we suggest that the ngEHT conduct VLBI follow-up of targets within no more than a month. Our imaging simulations were done assuming a full ngEHT array consisting of the 2022 EHT stations plus 11 additional sites, but depending on the target, not all sites may need to be available in order to obtain a high-fidelity image reconstruction.
Unlike many other EHT/ngEHT targets, TDEs will appear randomly across the sky, and the ngEHT will need to be capable of follow-up observations in the order of a week to weeks. Our current modeling of jets and outflows from super-Eddington disks is too sparse to let us make predictions regarding how long the jets will be visible at 230 GHz. However, if radio TDE observations are any indicator, emissions may persist for many months [12].
TDEs provide an excellent laboratory for studying jet/accretion/black hole systems across a wide range of accretion states over a relatively short period of time (∼1-a few years). In several cases, TDEs have shown state transitions after several hundred days in the X-ray range, which are likely associated with the evolution of the disk as the mass accretion rate declines. Stone and Metzger [29] argue, for instance, that the transition from a thick, super-Eddington disk to a thin disk can explain the jet shut-off in jetted TDEs, such as Swift J1644+57 and Swift J2058+05 [56,57], but recent simulations [58,59] demonstrated that MAD is possible even for thinner accretion disks. As such, long-term VLBI monitoring is strongly suggested. as this would allow for (1) the radio-submm emission of the outflows to be characterized and compared to the behavior of the accretion flow, and (2) the direct study of how the jet evolves morphologically as the disk state changes.
Another attractive potential target which we have yet to attribute a self-contained study to is jetted TDEs. These TDEs are extremely rare, and current observations suggest only about 1% of all TDEs will produce powerful relativistic jets. These jets will produce extremely bright radiation in the X-ray and the radio-submm regions. However, most have been distant due to the lower probability of their occurrence. Should a jetted TDE occur nearby enough for VLBI to resolve the jet, we strongly suggest such jets be treated as targets for the ngEHT.
Lastly, a recent TDE AT2018hyz showed a late outflow (∼3 years after the initial outburst) and brightened in the radio over several hundred days [60]. Unlike many other radio TDEs, a relatively bright 240 GHz emission was detected. If AT2018hyz is in fact a jet instead of a spherical outflow, Cendes et al. [60] estimates that the velocity could reach ≲0.6c. AT2018hyz is a relatively nearby TDE at ∼204 Mpc, but the flux density at the time of detection (∼0.2 mJy at 240 GHz) makes it too dim for ngEHT follow-up. However, placing AT2018hyz at ∼50 Mpc would shift the flux density to ∼5 mJy, which is the minimum estimated flux density required for an ngEHT VLBI detection. Future TDEs will likely be monitored across the radio-submm, so nearby targets of opportunity such as late radio TDEs like AT2018hyz should be considered, should they show significant radio emission.
While the EHT operated using dedicated time, the ngEHT will include more arrays and some telescopes may be sequestered for targeting of opportune VLBI observations (Doeleman et al., in prep). Development and proposals for the ngEHT are still underway, but we expect that fast (within 1-several weeks) response coverage of sources of interest in order to conduct VLBI of TDEs (and other transient sources) should be possible (Shepard Doeleman, private communication). The Atacama Large Millimeter Array (or ALMA), the most sensitive antenna, will maintain sustained coverage [40]. As such, we anticipate that at least unresolved radio data can be captured, which would prove equally invaluable to the study of radio TDEs.

5. Conclusions

In this work, we have demonstrated through a synthetic imaging analysis that TDE jets resembling the GRRMHD models presented in Curd et al. [14] are compelling ngEHT targets. We also confirmed that the detection limits considered in Curd et al. [14] are roughly applicable, as m7a0.0-HR did not produce detectable emission at a distance of 100 Mpc in our imaging analysis.
Various shock features in the jet are visible for the 10 Mpc images we consider, and studying the jet morphology in these cases could aid in characterizing the environment of the BH. Most TDEs that occur during the ngEHT mission will be farther away, but the apparent motion, which may be superluminal, can be extracted in such cases.
We suggest that the ngEHT be utilized for radio follow-up of TDEs. Our models study the birth of a TDE jet in the first 48 days after the disk forms under the assumption that the disk is SANE, super-Eddington, and threaded by a dynamically important magnetic field. However, TDEs which occur nearby and have jet properties similar to those of jetted TDEs, such as Swift J1644+57 or AT2018hyz, may provide interesting targets.

Author Contributions

Conceptualization, R.E. and R.A.; Formal analysis, B.C. and F.R.; Investigation, B.C.; Writing—original draft, B.C., F.R. and R.A.; Writing—review & editing, B.C. All authors have read and agreed to the published version of the manuscript.

Funding

Brandon Curd was supported by NSF grant AST-1816420, and made use of computational support from NSF via XSEDE/ACCESS resources (grant TG-AST080026N). Razieh Emami acknowledges the support by the Institute for Theory and Computation at the Center for Astrophysics and grants 21-atp21-0077, NSF AST-1816420, and HST-GO-16173.001-A for very generous support. Freek Roelofs was supported by NSF grants AST-1935980 and AST-2034306. This work was supported by the Black Hole Initiative at Harvard University, made possible through the support of grants from the Gordon and Betty Moore Foundation and the John Templeton Foundation. The opinions expressed in this publication are those of the authors and do not necessarily reflect the views of the Moore or Templeton Foundations.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

We graciously acknowledge Koushik Chatterjee, Lani Oramas, Joaquin Duran, and Hayley West for fruitful discussions that were helpful in the preparation of this work.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Full Image Library

In Figure A1Figure A4, we show the full library of images analyzed in this work. All of the base images were ray traced at ν = 230 GHz and then convolved with a Gaussian beam with a FWHM of 20 μ as . Similarly, we blurred the reconstructed image using the same beam for comparison. Note that for Figure A2 and Figure A4, we also shifted the reconstructed image to be approximately centered for comparison with the base image. In general, bright features are represented quite well in the reconstruction; however, some noise is introduced in dimmer sources (i.e., see the reconstruction of m7a0.9-HR at D = 100 Mpc).
Figure A1. Here we show the full library of base ipole images for model m7a0.0-HR. The time of each column is indicated at the top; the distance D and plasma temperature ratio R for each row are indicated on the right. The angle of the observer relative to the jet axis is indicated above each set of two rows. Each image spans 550 × 770 μ as 2 and is blurred by convolving the base image with a Gaussian beam with a FWHM of 20 μ as. The color scale is logarithmic, spanning three orders of magnitude, and each image uses the same maximum for the intensity scale.
Figure A1. Here we show the full library of base ipole images for model m7a0.0-HR. The time of each column is indicated at the top; the distance D and plasma temperature ratio R for each row are indicated on the right. The angle of the observer relative to the jet axis is indicated above each set of two rows. Each image spans 550 × 770 μ as 2 and is blurred by convolving the base image with a Gaussian beam with a FWHM of 20 μ as. The color scale is logarithmic, spanning three orders of magnitude, and each image uses the same maximum for the intensity scale.
Galaxies 10 00117 g0a1
Figure A2. The same as Figure A1 but showing the reconstructed images. Note that the intensity scale for each panel is the same as the corresponding panel in Figure A1 for comparison.
Figure A2. The same as Figure A1 but showing the reconstructed images. Note that the intensity scale for each panel is the same as the corresponding panel in Figure A1 for comparison.
Galaxies 10 00117 g0a2
Figure A3. Here we show the full library of base ipole images for model m7a0.9-HR. Each image in the top and middle rows spans 1000 × 1400 μ as 2 , and the bottom row spans 170 × 238 μ as 2 . Each image was blurred by convolving the base image with a Gaussian beam with a FWHM of 20 μ as (indicated by the white circle in the bottom right panel). The time of each column is indicated at the top; the distance D and plasma temperature ratio R for each row are indicated on the right. The angle of the observer relative to the jet axis is indicated above each set of two rows. The color scale is logarithmic, spanning three orders of magnitude in each image. We used the same color scale for the D = 10 Mpc images, but reduced the maximum by an order of magnitude in the D = 100 Mpc images to better show the image features.
Figure A3. Here we show the full library of base ipole images for model m7a0.9-HR. Each image in the top and middle rows spans 1000 × 1400 μ as 2 , and the bottom row spans 170 × 238 μ as 2 . Each image was blurred by convolving the base image with a Gaussian beam with a FWHM of 20 μ as (indicated by the white circle in the bottom right panel). The time of each column is indicated at the top; the distance D and plasma temperature ratio R for each row are indicated on the right. The angle of the observer relative to the jet axis is indicated above each set of two rows. The color scale is logarithmic, spanning three orders of magnitude in each image. We used the same color scale for the D = 10 Mpc images, but reduced the maximum by an order of magnitude in the D = 100 Mpc images to better show the image features.
Galaxies 10 00117 g0a3
Figure A4. The same as Figure A3 but showing the reconstructed images. Note that the intensity scale for each panel is the same as the corresponding panel in Figure A3 for comparison.
Figure A4. The same as Figure A3 but showing the reconstructed images. Note that the intensity scale for each panel is the same as the corresponding panel in Figure A3 for comparison.
Galaxies 10 00117 g0a4

Appendix B. Fitting Procedure for Jet Head Position

Here we describe the algorithm implemented to estimate the positions of the top and bottom jet heads in each image. In order to track the jet’s motion, we first smoothed the data (either the base image or the reconstructed image) with a Gaussian beam which assumed an angular resolution for the VLBI observations of Δ θ = 20 μ as . We then binned the data along the symmetry axis of the jet (y) by summing the flux along each row (x).
Since the images appear to be roughly Lorentzian, we first attempted to fit a double Lorentzian function of the form:
f d ( y ) = a 1 w 1 2 ( y y cen , 2 ) 2 + w 1 2 + a 2 w 2 2 ( y y cen , 2 ) 2 + w 2 2 ,
where a is the amplitude, w is the width, and y cen is the center defining the curve. We took y cen as a measure of the jet head location. If this fitting procedure did not produce a good fit, we found the peaks by performing a two-step fitting procedure in which we fit a single Lorentzian:
f 1 ( y ) = a 1 w 1 2 ( y y cen , 1 ) 2 + w 1 2 ,
and then subtracted the fit f 1 ( y ) from the data and then fit the second peak with:
f 2 ( y ) = a 2 w 2 2 ( y y cen , 2 ) 2 + w 2 2 .
We implemented the Python package SciPy [61] to optimize the curve(s) and estimate the jet head positions and errors. We show an example of the data and the Lorentzian fit for both the base image and the reconstruction in Figure A5. The reconstruction tends to be a bit broader, but the fitting procedure works equally well for all images and reconstructed images.
Figure A5. We demonstrate the fit performance using a ray-traced image of model m7a0.9-HR at t = 78 , 000 t g , θ = 90 , D = 10 Mpc, and R = 20 . We show the x-binned data and a double Lorentzian (Equation (A1)) fit for (a) the blurred base image and (b) the blurred reconstruction.
Figure A5. We demonstrate the fit performance using a ray-traced image of model m7a0.9-HR at t = 78 , 000 t g , θ = 90 , D = 10 Mpc, and R = 20 . We show the x-binned data and a double Lorentzian (Equation (A1)) fit for (a) the blurred base image and (b) the blurred reconstruction.
Galaxies 10 00117 g0a5

Note

1
https://github.com/Smithsonian/ngehtsim (accessed on 1 June 2022).

References

  1. Hills, J.G. Possible power source of Seyfert galaxies and QSOs. Nature 1975, 254, 295–298. [Google Scholar] [CrossRef] [Green Version]
  2. Rees, M.J. Tidal disruption of stars by black holes of 106–108 solar masses in nearby galaxies. Nature 1988, 333, 523–528. [Google Scholar] [CrossRef]
  3. Guillochon, J.; Ramirez-Ruiz, E. Hydrodynamical Simulations to Determine the Feeding Rate of Black Holes by the Tidal Disruption of Stars: The Importance of the Impact Parameter and Stellar Structure. Astrophys. J. 2013, 767, 25. [Google Scholar] [CrossRef] [Green Version]
  4. Mainetti, D.; Lupi, A.; Campana, S.; Colpi, M.; Coughlin, E.R.; Guillochon, J.; Ramirez-Ruiz, E. The fine line between total and partial tidal disruption events. Astron. Astrophys. 2017, 600, A124. [Google Scholar] [CrossRef]
  5. Steinberg, E.; Stone, N.C. The Origins of Peak Light in Tidal Disruption Events. arXiv 2022, arXiv:2206.10641. [Google Scholar] [CrossRef]
  6. Stone, N.; Sari, R.; Loeb, A. Consequences of strong compression in tidal disruption events. MNRAS 2013, 435, 1809–1824. [Google Scholar] [CrossRef] [Green Version]
  7. Abramowicz, M.A.; Calvani, M.; Nobili, L. Thick accretion disks with super-Eddington luminosities. Astrophys. J. 1980, 242, 772–788. [Google Scholar] [CrossRef]
  8. Abramowicz, M.A.; Czerny, B.; Lasota, J.P.; Szuszkiewicz, E. Slim Accretion Disks. Astrophys. J. 1988, 332, 646. [Google Scholar] [CrossRef]
  9. Golightly, E.C.A.; Nixon, C.J.; Coughlin, E.R. On the Diversity of Fallback Rates from Tidal Disruption Events with Accurate Stellar Structure. Astrophys. J. 2019, 882, L26. [Google Scholar] [CrossRef]
  10. Komossa, S. Tidal disruption of stars by supermassive black holes: Status of observations. J. High Energy Astrophys. 2015, 7, 148–157. [Google Scholar] [CrossRef]
  11. Gezari, S. Tidal Disruption Events. ARA&A 2021, 59, 21–58. [Google Scholar] [CrossRef]
  12. Alexander, K.D.; van Velzen, S.; Horesh, A.; Zauderer, B.A. Radio Properties of Tidal Disruption Events. Space Sci. Rev. 2020, 216, 81. [Google Scholar] [CrossRef]
  13. Dai, L.; McKinney, J.C.; Roth, N.; Ramirez-Ruiz, E.; Miller, M.C. A Unified Model for Tidal Disruption Events. Astrophys. J. 2018, 859, L20. [Google Scholar] [CrossRef] [Green Version]
  14. Curd, B.; Emami, R.; Anantua, R.; Palumbo, D.; Doeleman, S.; Narayan, R. Jets from SANE Super-Eddington Accretion Disks: Morphology, Spectra, and Their Potential as Targets for ngEHT. arXiv 2022, arXiv:2206.06358. [Google Scholar]
  15. Gammie, C.F.; McKinney, J.C.; Tóth, G. HARM: A Numerical Scheme for General Relativistic Magnetohydrodynamics. Astrophys. J. 2003, 589, 444–457. [Google Scholar] [CrossRef]
  16. Tchekhovskoy, A.; Metzger, B.D.; Giannios, D.; Kelley, L.Z. Swift J1644+57 gone MAD: The case for dynamically important magnetic flux threading the black hole in a jetted tidal disruption event. MNRAS 2014, 437, 2744–2760. [Google Scholar] [CrossRef] [Green Version]
  17. Curd, B.; Narayan, R. GRRMHD simulations of tidal disruption event accretion discs around supermassive black holes: Jet formation, spectra, and detectability. MNRAS 2019, 483, 565–592. [Google Scholar] [CrossRef] [Green Version]
  18. Blandford, R.D.; Znajek, R.L. Electromagnetic extraction of energy from Kerr black holes. MNRAS 1977, 179, 433–456. [Google Scholar] [CrossRef]
  19. Abramowicz, M.A.; Ellis, G.F.R.; Lanza, A. Relativistic Effects in Superluminal Jets and Neutron Star Winds. Astrophys. J. 1990, 361, 470. [Google Scholar] [CrossRef]
  20. Coughlin, E.R.; Begelman, M.C. Structured, relativistic jets driven by radiation. MNRAS 2020, 499, 3158–3177. [Google Scholar] [CrossRef]
  21. Dadhich, N.; Tursunov, A.; Ahmedov, B.; Stuchlík, Z. The distinguishing signature of magnetic Penrose process. MNRAS 2018, 478, L89–L94. [Google Scholar] [CrossRef]
  22. Stuchlík, Z.; Kološ, M.; Kovář, J.; Slaný, P.; Tursunov, A. Influence of Cosmic Repulsion and Magnetic Fields on Accretion Disks Rotating around Kerr Black Holes. Universe 2020, 6, 26. [Google Scholar] [CrossRef] [Green Version]
  23. Penrose, R. Gravitational Collapse: The Role of General Relativity. Nuovo Cim. Riv. Ser. 1969, 1, 252. [Google Scholar]
  24. Thomsen, L.L.; Kwan, T.M.; Dai, L.; Wu, S.C.; Roth, N.; Ramirez-Ruiz, E. Dynamical Unification of Tidal Disruption Events. Astrophys. J. 2022, 937, L28. [Google Scholar] [CrossRef]
  25. Penna, R.F.; Narayan, R.; Sądowski, A. General relativistic magnetohydrodynamic simulations of Blandford-Znajek jets and the membrane paradigm. MNRAS 2013, 436, 3741–3758. [Google Scholar] [CrossRef] [Green Version]
  26. Doeleman, S.; Blackburn, L.; Dexter, J.; Gomez, J.L.; Johnson, M.D.; Palumbo, D.C.; Weintroub, J.; Farah, J.R.; Fish, V.; Loinard, L.; et al. Studying Black Holes on Horizon Scales with VLBI Ground Arrays. Bull. Am. Astron. Soc. 2019, 51, 256. [Google Scholar]
  27. Ivezić, Ž.; Kahn, S.M.; Tyson, J.A.; Abel, B.; Acosta, E.; Allsman, R.; Alonso, D.; AlSayyad, Y.; Anderson, S.F.; Andrew, J.; et al. LSST: From Science Drivers to Reference Design and Anticipated Data Products. Astrophys. J. 2019, 873, 111. [Google Scholar] [CrossRef]
  28. Bricman, K.; Gomboc, A. The Prospects of Observing Tidal Disruption Events with the Large Synoptic Survey Telescope. Astrophys. J. 2020, 890, 73. [Google Scholar] [CrossRef]
  29. Stone, N.C.; Metzger, B.D. Rates of stellar tidal disruption as probes of the supermassive black hole mass function. MNRAS 2016, 455, 859–883. [Google Scholar] [CrossRef] [Green Version]
  30. Sadowski, A.; Tejeda, E.; Gafton, E.; Rosswog, S.; Abarca, D. Magnetohydrodynamical simulations of a deep tidal disruption in general relativity. Mon. Not. Roy. Astron. Soc. 2016, 458, 4250–4268. [Google Scholar] [CrossRef] [Green Version]
  31. Curd, B. Global simulations of tidal disruption event disc formation via stream injection in GRRMHD. Mon. Not. Roy. Astron. Soc. 2021, 507, 3207–3227. [Google Scholar] [CrossRef]
  32. Novikov, I.D.; Thorne, K.S. Astrophysics of black holes. In Black Holes (Les Astres Occlus); Gordon & Breach: New York, NY, USA, 1973; pp. 343–450. [Google Scholar]
  33. Sądowski, A.; Narayan, R.; Tchekhovskoy, A.; Abarca, D.; Zhu, Y.; McKinney, J.C. Global simulations of axisymmetric radiative black hole accretion discs in general relativity with a mean-field magnetic dynamo. MNRAS 2015, 447, 49–71. [Google Scholar] [CrossRef] [Green Version]
  34. Barniol Duran, R.; Tchekhovskoy, A.; Giannios, D. Simulations of AGN jets: Magnetic kink instability versus conical shocks. MNRAS 2017, 469, 4957–4978. [Google Scholar] [CrossRef]
  35. Mościbrodzka, M.; Gammie, C.F. IPOLE - semi-analytic scheme for relativistic polarized radiative transport. MNRAS 2018, 475, 43–54. [Google Scholar] [CrossRef] [Green Version]
  36. Yarza, R.; Wong, G.N.; Ryan, B.R.; Gammie, C.F. Bremsstrahlung in GRMHD Models of Accreting Black Holes. Astrophys. J. 2020, 898, 50. [Google Scholar] [CrossRef]
  37. Wong, G.N.; Prather, B.S.; Dhruv, V.; Ryan, B.R.; Mościbrodzka, M.; Chan, C.k.; Joshi, A.V.; Yarza, R.; Ricarte, A.; Shiokawa, H.; et al. PATOKA: Simulating Electromagnetic Observables of Black Hole Accretion. Astrophys. J. 2022, 259, 64. [Google Scholar] [CrossRef]
  38. Ohmura, T.; Machida, M.; Nakamura, K.; Kudoh, Y.; Asahina, Y.; Matsumoto, R. Two-Temperature Magnetohydrodynamics Simulations of Propagation of Semi-Relativistic Jets. Galaxies 2019, 7, 14. [Google Scholar] [CrossRef] [Green Version]
  39. Ohmura, T.; Machida, M.; Nakamura, K.; Kudoh, Y.; Matsumoto, R. Two-temperature magnetohydrodynamic simulations for sub-relativistic active galactic nucleus jets: Dependence on the fraction of the electron heating. MNRAS 2020, 493, 5761–5772. [Google Scholar] [CrossRef] [Green Version]
  40. Raymond, A.W.; Palumbo, D.; Paine, S.N.; Blackburn, L.; Córdova Rosado, R.; Doeleman, S.S.; Farah, J.R.; Johnson, M.D.; Roelofs, F.; Tilanus, R.P.J.; et al. Evaluation of New Submillimeter VLBI Sites for the Event Horizon Telescope. Astrophys. J. 2021, 253, 5. [Google Scholar] [CrossRef]
  41. Roelofs, F.; et al. [Black Hole Initiative at Harvard University], 2022; in preparation.
  42. Chael, A.A.; Johnson, M.D.; Narayan, R.; Doeleman, S.S.; Wardle, J.F.C.; Bouman, K.L. High-resolution Linear Polarimetric Imaging for the Event Horizon Telescope. Astrophys. J. 2016, 829, 11. [Google Scholar] [CrossRef]
  43. Chael, A.A.; Johnson, M.D.; Bouman, K.L.; Blackburn, L.L.; Akiyama, K.; Narayan, R. Interferometric Imaging Directly with Closure Phases and Closure Amplitudes. Astrophys. J. 2018, 857, 23. [Google Scholar] [CrossRef]
  44. Doeleman, S.; et al. [Harvard-Smithsonian Center for Astrophysics], 2022; in preparation.
  45. Gelaro, R.; McCarty, W.; Suárez, M.J.; Todling, R.; Molod, A.; Takacs, L.; Randles, C.A.; Darmenov, A.; Bosilovich, M.G.; Reichle, R.; et al. The Modern-Era Retrospective Analysis for Research and Applications, Version 2 (MERRA-2). J. Clim. 2017, 30, 5419–5454. [Google Scholar] [CrossRef]
  46. Paine, S. The Am Atmospheric Model. 2019. Available online: https://zenodo.cern.ch/record/3406496#.Y5fRyH1ByUk (accessed on 1 June 2022).
  47. Event Horizon Telescope Collaboration; Akiyama, K.; Alberdi, A.; Alef, W.; Asada, K.; Azulay, R.; Baczko, A.K.; Ball, D.; Baloković, M.; Barrett, J.; et al. First M87 Event Horizon Telescope Results. IV. Imaging the Central Supermassive Black Hole. Astrophys. J. 2019, 875, L4. [Google Scholar] [CrossRef]
  48. Jorstad, S.G.; Marscher, A.P.; Lister, M.L.; Stirling, A.M.; Cawthorne, T.V.; Gear, W.K.; Gómez, J.L.; Stevens, J.A.; Smith, P.S.; Forster, J.R.; et al. Polarimetric Observations of 15 Active Galactic Nuclei at High Frequencies: Jet Kinematics from Bimonthly Monitoring with the Very Long Baseline Array. Astrophys. J. 2005, 130, 1418–1465. [Google Scholar] [CrossRef]
  49. Lister, M.L.; Aller, M.F.; Aller, H.D.; Homan, D.C.; Kellermann, K.I.; Kovalev, Y.Y.; Pushkarev, A.B.; Richards, J.L.; Ros, E.; Savolainen, T. MOJAVE. X. Parsec-scale Jet Orientation Variations and Superluminal Motion in Active Galactic Nuclei. Astrophys. J. 2013, 146, 120. [Google Scholar] [CrossRef] [Green Version]
  50. Cohen, M.H.; Meier, D.L.; Arshakian, T.G.; Homan, D.C.; Hovatta, T.; Kovalev, Y.Y.; Lister, M.L.; Pushkarev, A.B.; Richards, J.L.; Savolainen, T. Studies of the Jet in Bl Lacertae. I. Recollimation Shock and Moving Emission Features. Astrophys. J. 2014, 787, 151. [Google Scholar] [CrossRef]
  51. Kohler, S.; Begelman, M.C.; Beckwith, K. Recollimation boundary layers in relativistic jets. MNRAS 2012, 422, 2282–2290. [Google Scholar] [CrossRef] [Green Version]
  52. Lazzati, D.; Morsony, B.J.; Blackwell, C.H.; Begelman, M.C. Unifying the Zoo of Jet-driven Stellar Explosions. Astrophys. J. 2012, 750, 68. [Google Scholar] [CrossRef]
  53. Mizuno, Y.; Gómez, J.L.; Nishikawa, K.I.; Meli, A.; Hardee, P.E.; Rezzolla, L. Recollimation Shocks in Magnetized Relativistic Jets. Astrophys. J. 2015, 809, 38. [Google Scholar] [CrossRef] [Green Version]
  54. Hervet, O.; Meliani, Z.; Zech, A.; Boisson, C.; Cayatte, V.; Sauty, C.; Sol, H. Shocks in relativistic transverse stratified jets. A new paradigm for radio-loud AGN. Astron. Astrophys. 2017, 606, A103. [Google Scholar] [CrossRef] [Green Version]
  55. Gómez, J.L.; Lobanov, A.P.; Bruni, G.; Kovalev, Y.Y.; Marscher, A.P.; Jorstad, S.G.; Mizuno, Y.; Bach, U.; Sokolovsky, K.V.; Anderson, J.M.; et al. Probing the Innermost Regions of AGN Jets and Their Magnetic Fields with RadioAstron. I. Imaging BL Lacertae at 21 Microarcsecond Resolution. Astrophys. J. 2016, 817, 96. [Google Scholar] [CrossRef] [Green Version]
  56. Zauderer, B.A.; Berger, E.; Margutti, R.; Pooley, G.G.; Sari, R.; Soderberg, A.M.; Brunthaler, A.; Bietenholz, M.F. Radio Monitoring of the Tidal Disruption Event Swift J164449.3+573451. II. The Relativistic Jet Shuts Off and a Transition to Forward Shock X-Ray/Radio Emission. Astrophys. J. 2013, 767, 152. [Google Scholar] [CrossRef]
  57. Pasham, D.R.; Cenko, S.B.; Levan, A.J.; Bower, G.C.; Horesh, A.; Brown, G.C.; Dolan, S.; Wiersema, K.; Filippenko, A.V.; Fruchter, A.S.; et al. A Multiwavelength Study of the Relativistic Tidal Disruption Candidate Swift J2058.4+0516 at Late Times. Astrophys. J. 2015, 805, 68. [Google Scholar] [CrossRef]
  58. Curd, B.; Narayan, R. GRRMHD Simulations of MAD Accretion Disks Declining from Super-Eddington to Sub-Eddington Accretion Rates. arXiv 2022, arXiv:2209.12081. [Google Scholar]
  59. Liska, M.T.P.; Musoke, G.; Tchekhovskoy, A.; Porth, O.; Beloborodov, A.M. Formation of Magnetically Truncated Accretion Disks in 3D Radiation-transport Two-temperature GRMHD Simulations. Astrophys. J. 2022, 935, L1. [Google Scholar] [CrossRef]
  60. Cendes, Y.; Berger, E.; Alexander, K.D.; Gomez, S.; Hajela, A.; Chornock, R.; Laskar, T.; Margutti, R.; Metzger, B.; Bietenholz, M.F.; et al. A Mildly Relativistic Outflow Launched Two Years after Disruption in Tidal Disruption Event AT2018hyz. Astrophys. J. 2022, 938, 28. [Google Scholar] [CrossRef]
  61. Virtanen, P.; Gommers, R.; Oliphant, T.E.; Haberland, M.; Reddy, T.; Cournapeau, D.; Burovski, E.; Peterson, P.; Weckesser, W.; Bright, J.; et al. SciPy 1.0: Fundamental Algorithms for Scientific Computing in Python. Nat. Methods 2020, 17, 261–272. [Google Scholar] [CrossRef]
Figure 1. Here we show snapshots of the GRRMHD KORAL simulations that we post-processed with ipole. All data are shown for t = 78 , 000 t g in both m7a0.0-HR (top) and m7a0.9-HR (bottom). The colors indicate the gas density ρ (left), gas temperature T gas (middle), and magnetic field strength | B | (right); the yellow contours indicate the σ = 1 boundary, in which we set ρ = 0 in the ray tracing step to prevent emission.
Figure 1. Here we show snapshots of the GRRMHD KORAL simulations that we post-processed with ipole. All data are shown for t = 78 , 000 t g in both m7a0.0-HR (top) and m7a0.9-HR (bottom). The colors indicate the gas density ρ (left), gas temperature T gas (middle), and magnetic field strength | B | (right); the yellow contours indicate the σ = 1 boundary, in which we set ρ = 0 in the ray tracing step to prevent emission.
Galaxies 10 00117 g001
Figure 2. Model m7a0.9-HR at t = 78 , 000 t g imaged at D = 10 Mpc with θ = 90 and R = 1 . We show the base ipole image with no blurring (left), the base ipole image blurred via convolution with a 20 μ as FWHM Gaussian beam (middle), and the reconstructed image blurred using the same Gaussian beam (right).
Figure 2. Model m7a0.9-HR at t = 78 , 000 t g imaged at D = 10 Mpc with θ = 90 and R = 1 . We show the base ipole image with no blurring (left), the base ipole image blurred via convolution with a 20 μ as FWHM Gaussian beam (middle), and the reconstructed image blurred using the same Gaussian beam (right).
Galaxies 10 00117 g002
Figure 3. The same as Figure 2 but for model m7a0.9-HR at t = 78 , 000 t g imaged at D = 100 Mpc with θ = 90 and R = 1 .
Figure 3. The same as Figure 2 but for model m7a0.9-HR at t = 78 , 000 t g imaged at D = 100 Mpc with θ = 90 and R = 1 .
Galaxies 10 00117 g003
Figure 4. Here we show the relative difference (open circles) and the error within 3 standard deviations (horizontal bars) between the jet lengths obtained from the ray-traced images ( l jet ) and the reconstructed images ( l jet , rec ) . Models m7a0.0-HR (left panel) and m7a0.9-HR (right panel) are shown. For each choice of viewing angle θ , distance to the source D, and R , we show the data at t = 38 , 000 t g (blue) and t = 78 , 000 r g . We indicate the viewing angle next to each pairing of error bars. For each model, except where explicitly indicated to be otherwise, we show D = 10 Mpc and R = 1 . In general, there is a very small relative difference when the jet is viewed edge on. As the viewing angle approaches face on or the jet is placed at a larger distance, the relative difference increases along with the error. In all but one image, even when the relative difference shifts away from zero, the zero relative difference line is within 3 standard deviations, which suggests excellent agreement between jet lengths derived from the raw ipole and reconstructed images.
Figure 4. Here we show the relative difference (open circles) and the error within 3 standard deviations (horizontal bars) between the jet lengths obtained from the ray-traced images ( l jet ) and the reconstructed images ( l jet , rec ) . Models m7a0.0-HR (left panel) and m7a0.9-HR (right panel) are shown. For each choice of viewing angle θ , distance to the source D, and R , we show the data at t = 38 , 000 t g (blue) and t = 78 , 000 r g . We indicate the viewing angle next to each pairing of error bars. For each model, except where explicitly indicated to be otherwise, we show D = 10 Mpc and R = 1 . In general, there is a very small relative difference when the jet is viewed edge on. As the viewing angle approaches face on or the jet is placed at a larger distance, the relative difference increases along with the error. In all but one image, even when the relative difference shifts away from zero, the zero relative difference line is within 3 standard deviations, which suggests excellent agreement between jet lengths derived from the raw ipole and reconstructed images.
Galaxies 10 00117 g004
Table 1. We tabulate the mass accretion rate M ˙ , jet efficiency η jet L jet / M ˙ c 2 (where L jet is the jet power as defined in [14]), and total simulation duration t sim for each KORAL simulation. Note that M ˙ and η jet are time averaged over the final 50,000 t g of each simulation.
Table 1. We tabulate the mass accretion rate M ˙ , jet efficiency η jet L jet / M ˙ c 2 (where L jet is the jet power as defined in [14]), and total simulation duration t sim for each KORAL simulation. Note that M ˙ and η jet are time averaged over the final 50,000 t g of each simulation.
Model M ˙ η jet t sim
( M ˙ Edd ) ( t g )
m7a0.0-HR120.24%83,000
m7a0.9-HR251.15%81,200
Table 2. Here we tabulate the 230 GHz flux density for each model given a specific time, viewing angle θ , distance D, and temperature ratio R .
Table 2. Here we tabulate the 230 GHz flux density for each model given a specific time, viewing angle θ , distance D, and temperature ratio R .
ModelTimeDistance R F 230 GHz
( t g ) (Mpc) (Jy)
θ = 10 θ = 45 θ = 90
m7a0.0-HR38,0001010.2190.2140.074
78,0001010.0140.0130.006
m7a0.9-HR38,0001012.0014.4526.036
78,00010111.96826.78035.092
38,0001020--0.190
78,0001020--0.485
38,0001001--0.060
78,0001001--0.351
Table 3. Here we tabulate the estimated jet length for each model at each time for various choices of the distance D, viewing angle θ , and plasma temperature ratio R . We compare the jet length as computed from the base ipole image ( l jet ) and the reconstructed image ( l jet , rec ).
Table 3. Here we tabulate the estimated jet length for each model at each time for various choices of the distance D, viewing angle θ , and plasma temperature ratio R . We compare the jet length as computed from the base ipole image ( l jet ) and the reconstructed image ( l jet , rec ).
ModelTimeDistance θ R l jet l jet , rec
( t g ) (Mpc) ( r g ) ( r g )
m7a0.0-HR38,00010 10 1 5091 272 + 272 4927 317 + 317
78,00010 10 1 9042 356 + 356 9104 407 + 407
38,00010 45 1 22 , 006 123 + 123 21 , 774 731 + 731
78,00010 45 1 39 , 832 318 + 318 39 , 678 501 + 501
38,00010 90 1 28 , 658 + 148 28 , 463 173 + 173
78,00010 90 1 57 , 068 401 + 401 56 , 844 237 + 237
m7a0.9-HR38,00010 10 1 8457 446 + 446 4191 726 + 726
78,00010 10 1 17 , 585 233 + 233 16 , 837 726 + 726
38,00010 45 1 33 , 128 336 + 336 31 , 941 683 + 683
78,00010 45 1 73 , 535 890 + 890 71 , 855 1058 + 1058
38,00010 90 1 48 , 468 306 + 306 46 , 972 1028 + 1028
78,00010 90 1 111 , 503 1085 + 1085 110 , 327 2657 + 2657
38,00010 90 20 49 , 457 150 + 150 48 , 930 358 + 358
78,00010 90 20 119 , 897 173 + 173 119 , 256 415 + 415
38,000100 90 1 46 , 715 400 + 400 40 , 244 2714 + 2714
78,000100 90 1 109 , 091 959 + 959 98 , 822 2075 + 2075
Table 4. Here we tabulate the estimated jet velocity for each model for different choices of the distance D, viewing angle θ , and plasma temperature ratio R . The velocities shown were calculated using the base ipole images (v), the reconstructed images ( v rec ), and the reonstructed images, while accounting for the possibility of superluminal motion ( v app , rec ).
Table 4. Here we tabulate the estimated jet velocity for each model for different choices of the distance D, viewing angle θ , and plasma temperature ratio R . The velocities shown were calculated using the base ipole images (v), the reconstructed images ( v rec ), and the reonstructed images, while accounting for the possibility of superluminal motion ( v app , rec ).
ModelDistance θ R v v rec v app , rec
(Mpc) ( c ) ( c ) ( c )
m7a0.0-HR10 10 1 0 . 049 0.006 + 0.006 0 . 052 0.006 + 0.006 0 . 074 0.009 + 0.009
10 45 1 0 . 223 0.005 + 0.005 0 . 224 0.011 + 0.011 0 . 288 0.014 + 0.014
10 90 1 0 . 355 0.005 + 0.005 0 . 355 0.004 + 0.004 0 . 355 0.004 + 0.004
m7a0.9-HR10 10 1 0 . 114 0.006 + 0.006 0 . 158 0.011 + 0.011 1 . 526 0.102 + 0.102
10 45 1 0 . 505 0.012 + 0.012 0 . 499 0.016 + 0.016 0 . 996 0.031 + 0.031
10 90 1 0 . 788 0.014 + 0.014 0 . 792 0.036 + 0.036 0 . 792 0.036 + 0.036
10 90 20 0 . 88 0.002 + 0.002 0 . 879 0.007 + 0.007 0 . 879 0.007 + 0.007
100 90 1 0 . 780 0.013 0.013 0 . 732 0.043 + 0.043 0 . 732 0.043 + 0.043
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Curd, B.; Emami, R.; Roelofs, F.; Anantua, R. Modeling Reconstructed Images of Jets Launched by SANE Super-Eddington Accretion Flows around SMBHs with the ngEHT. Galaxies 2022, 10, 117. https://doi.org/10.3390/galaxies10060117

AMA Style

Curd B, Emami R, Roelofs F, Anantua R. Modeling Reconstructed Images of Jets Launched by SANE Super-Eddington Accretion Flows around SMBHs with the ngEHT. Galaxies. 2022; 10(6):117. https://doi.org/10.3390/galaxies10060117

Chicago/Turabian Style

Curd, Brandon, Razieh Emami, Freek Roelofs, and Richard Anantua. 2022. "Modeling Reconstructed Images of Jets Launched by SANE Super-Eddington Accretion Flows around SMBHs with the ngEHT" Galaxies 10, no. 6: 117. https://doi.org/10.3390/galaxies10060117

APA Style

Curd, B., Emami, R., Roelofs, F., & Anantua, R. (2022). Modeling Reconstructed Images of Jets Launched by SANE Super-Eddington Accretion Flows around SMBHs with the ngEHT. Galaxies, 10(6), 117. https://doi.org/10.3390/galaxies10060117

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop