We present progress incorporating non-LTE (non- Local Thermodynamic Equilibrium) effects into the plasma equation of state (EOS) in the Ares multi-physics code [1] developed at Lawrence Livermore National Laboratory (LLNL). Accounting for non-LTE effects is crucial for obtaining accurate simulations of hohlraums and imploding capsules for inertial confinement fusion applications on the National Ignition Facility. The non-LTE atomic kinetics are tracked via routines originally developed for the LLNL plasma radiative transfer code CRETIN [2]. Past work has utilized the CRETIN routines to generate non-LTE opacities, emissivities, and ionization states in Ares, but with limited self-consistency regarding other EOS quantities, particularly the material specific heat and pressure. We include test problems related to laser-illuminated hohlraums illustrating our progress in achieving this broader self-consistency, including comparisons to KULL [3], another LLNL multi-physics code.

We compare X-ray emission from several general relativistic, multi-frequency, radiation magnetohydrodynamic simulations of thin black hole accretion disks with different accretion rates and spins. The simulations were performed using the M1 closure scheme, resolved with twelve frequency (energy) bins logarithmically spaced from 5 × 10^-3 to 5 × 10^3 keV. We apply a general relativistic Monte Carlo transport code to post-process the simulation data with greater fidelity in frequency resolution and Compton scattering treatment. Despite the relatively few energy bins and Kompaneets approximation to Compton scattering utilized in the M1 method, we find generally good agreement between the methods. Both produce prominent thermal profiles with peaks around 2 - 2.5 keV, where agreement is particularly strong and representative of the soft state. Both also find weaker (lower luminosity) thermally sourced emission extending out to 100 keV due to the hotter innermost regions of the disks. Inverse Compton scattering becomes increasingly effective at hardening spectral outputs with increasing black hole spin, and becomes the dominant mechanism for photons that escape with energies between 10 to several hundred keV. At very high rates of spin the radiation flux in this upscattered component becomes comparable to the thermal flux, a phenomenon typically associated with intermediate states. Beyond 10^4 keV, we observe faint, free-free emission from hot, optically thin coronal regions developing near the horizon, common to both spinning and nonspinning black holes.

We compare X-ray emission from several general relativistic, multifrequency, radiation magnetohydrodynamic simulations of thin black hole accretion disks with different accretion rates and spins. The simulations were performed using the M _1 closure scheme, resolved with 12 frequency (energy) bins logarithmically spaced from 5 × 10 ^−3 to 5 × 10 ^3 ​​​ keV. We apply a general relativistic Monte Carlo transport code to postprocess the simulation data with greater fidelity in frequency resolution and Compton scattering treatment. Despite the relatively few energy bins and Kompaneets approximation to Compton scattering utilized in the M _1 method, we find generally good agreement between the methods. Both produce prominent thermal profiles with peaks around 2–2.5 keV, where agreement is particularly strong and representative of the soft state. Both also find weaker (lower-luminosity) thermally sourced emission extending out to 100 keV due to the hotter innermost regions of the disks. Inverse Compton scattering becomes increasingly effective at hardening spectral outputs with increasing black hole spin, and becomes the dominant mechanism for photons that escape with energies between 10 to several hundred keV. At very high rates of spin, the radiation flux in this upscattered component becomes comparable to the thermal flux, a phenomenon typically associated with intermediate states. Beyond 10 ^4 keV, we observe faint, free–free emission from hot, optically thin coronal regions developing near the horizon, common to both spinning and nonspinning black holes.

We present a set of six general relativistic, multifrequency, radiation magnetohydrodynamic simulations of thin accretion disks with different target mass accretion rates around black holes with spins ranging from nonrotating to rapidly spinning. The simulations use the M _1 closure scheme with 12 independent frequency (or energy) bins ranging logarithmically from 5 × 10 ^−3 keV to 5 × 10 ^3 keV. The multifrequency capability allows us to generate crude spectra and energy-dependent light curves directly from the simulations without a need for special postprocessing. While we generally find roughly thermal spectra with peaks around 1–4 keV, our high-spin cases showed harder-than-expected tails for the soft or thermally dominant state. This leads to radiative efficiencies that are up to five times higher than expected for a Novikov–Thorne disk at the same spin. We attribute these high efficiencies to the high-energy, coronal emission. These coronae mostly occupy the effectively optically thin regions near the inner edges of the disks and also cover or sandwich the inner ∼15 GM / c ^2 of the disks.

We report on a new capability added to our general relativistic radiation-magnetohydrodynamics code, Cosmos++: an implicit Monte Carlo (IMC) treatment for radiation transport. The method is based on a Fleck-type implicit discretization of the radiation-hydrodynamics equations, but generalized for both Newtonian and relativistic regimes. A multiple reference frame approach is used to geodesically transport photon packets (and solve the hydrodynamics equations) in the coordinate frame, while radiation-matter interactions are handled either in the fluid or electron frames then communicated via Lorentz boosts and orthonormal tetrad bases attached to the fluid. We describe a method for constructing estimators of radiation moments using path-weighting that generalizes to arbitrary coordinate systems in flat or curved spacetime. Absorption, emission, scattering, and relativistic Comptonization are among the matter interactions considered in this report. We discuss our formulations and numerical methods, and validate our models against a suite of radiation and coupled radiation-hydrodynamics test problems in both flat and curved spacetimes.

A star crossing the tidal radius of a supermassive black hole will be spectacularly ripped apart with an accompanying burst of radiation. A few tens of such tidal disruption events have now been identified in optical wavelengths, but the exact origin of the strong optical emission remains inconclusive. Here we report polarimetric observations of three tidal disruption events. The continuum polarization appears independent of wavelength, while emission lines are partially depolarized. These signatures are consistent with photons being scattered and polarized in an envelope of free electrons. An almost axisymmetric photosphere viewed from different angles is in broad agreement with the data, but there is also evidence for deviations from axial symmetry before the peak of the flare and significant time evolution at early times, compatible with the rapid formation of an accretion disk. By combining a super-Eddington accretion model with a radiative transfer code, we simulate the polarization degree as a function of disk mass and viewing angle and we show that the predicted levels are compatible with the observations for extended reprocessing envelopes of ~1,000 gravitational radii. Spectropolarimetry therefore constitutes a new observational test for tidal disruption event models, and opens an important new line of exploration in the study of tidal disruption events. Spectropolarimetric observations of three tidal disruption events reveal that they are optically polarized at the 1–2% level by a cloud of electrons surrounding the black hole (in good agreement with theory).

We present the first results of a comprehensive supernova (SN) radiative-transfer (RT) code-comparison initiative (StaNdaRT), where the emission from the same set of standardized test models is simulated by currently-used RT codes. A total of ten codes have been run on a set of four benchmark ejecta models of Type Ia supernovae. We consider two sub-Chandrasekhar-mass ($M_\mathrm{tot} = 1.0$ M$_\odot$) toy models with analytic density and composition profiles and two Chandrasekhar-mass delayed-detonation models that are outcomes of hydrodynamical simulations. We adopt spherical symmetry for all four models. The results of the different codes, including the light curves, spectra, and the evolution of several physical properties as a function of radius and time, are provided in electronic form in a standard format via a public repository. We also include the detailed test model profiles and several python scripts for accessing and presenting the input and output files. We also provide the code used to generate the toy models studied here. In this paper, we describe in detail the test models, radiative-transfer codes and output formats and provide access to the repository. We present example results of several key diagnostic features.

We study the initiation of thermonuclear detonations in tidally disrupted white dwarf stars by intermediate-mass (10 3 M ⊙ ) black holes. The length scales required to resolve the initiation mechanism are not easily reached in 3D, so instead we have devised 2D proxy models, which, together with a logarithmic gridding strategy, can adequately capture detonation wave fronts as material undergoes simultaneous compression and stretching from tidal forces. We consider 0.15 and 0.6 M ⊙ white dwarf stars parameterized by tidal strengths in the range β = 4–23. High spatial resolution elucidates the manner and conditions leading to thermonuclear detonation, linking the initiation sequence to stellar composition and tidal strength. All of our models suffer sustained detonations triggered by a combination of adiabatic compression, mild thermonuclear preconditioning, and collisional heating, in degrees depending primarily on tidal strength. We find that many diagnostics, such as temperature, total released energy, and iron-group products, are fairly well converged (better than 10%) at resolutions below 10 km along the scale height of the orbital plane. The exceptions are intermediate-mass transients like calcium, which remain uncertain up to factors of 2, even at 1 km resolution.

About a hundred tidal disruption events (TDEs) have been observed and they exhibit a wide range of emission properties both at peak and over their lifetimes. Some TDEs peak predominantly at X-ray energies while others radiate chiefly at UV and optical wavelengths. While the peak luminosities across TDEs show distinct properties, the evolutionary behavior can also vary between TDEs with similar peak emission properties. At late time, some optical TDEs rebrighten in X-rays, while others maintain strong UV/optical emission components. In this Letter, we conduct three-dimensional general relativistic radiation magnetohydrodynamics simulations of TDE accretion disks at varying accretion rates ranging from a few to a few tens of the Eddington accretion rate. We make use of Monte Carlo radiative transfer simulations to calculate the reprocessed spectra at various inclinations and at different evolutionary stages. We confirm the unified model proposed by Dai et al. (2018), which predicts that the observed emission largely depends on the viewing angle of the observer with respect to the disk orientation (X-ray strong when viewed face-on and UV/optically strong when viewed disk-on). What is more, we find that disks with higher accretion rates have elevated wind and disk densities, which increases the reprocessing of the high-energy radiation and thus generally augments the optical-to-X-ray flux ratio along a particular viewing angle. This implies that at later times, as the accretion level declines, we expect the funnel to effectively open up and allow more X-rays to leak out along intermediate viewing angles. Such dynamical model for TDEs can provide a natural explanation for the diversity in the emission properties observed in TDEs at peak and along their temporal evolution.

Detections of the tidal disruption flares (TDFs) of stars by supermassive black holes (SMBHs) are rapidly accumulating as optical surveys improve. These detections may provide constraints on SMBH demographics, stellar dynamics, and stellar evolution in galaxies. To maximize this scientific impact, we require a better understanding of how astrophysical parameters interact with survey selection effects in setting the properties of detected flares. We develop a framework for modeling the distributions of optical TDF detections in surveys across attributes of the host galaxies and the flares themselves. This model folds in effects of the stellar disruption rate in each galaxy, the flare luminosity and temperature distributions, the effects of obscuration and reddening by dust in the host galaxy, and survey selection criteria. We directly apply this model to the sample of TDFs detected by the Zwicky Transient Facility and find that the overall flare detection rate is in line with simple theoretical expectation. The model can also reproduce the distribution of total stellar mass and redshift of the host galaxies, but fails to match all details of the detected flares, such as their luminosity and temperature distributions. We also find that dust obscuration likely plays an important role in suppressing the TDF detection rate in star-forming galaxies. While we do not find that the unusual preference of TDFs to have hosts in post-starburst galaxies in the "green valley" can be entirely explained by selection effects, our model can help to quantify the true rate enhancement in those galaxies.

While tidal disruption events (TDEs) have long been heralded as laboratories for the study of quiescent black holes, the small number of known TDEs and uncertainties in their emission mechanism have hindered progress towards this promise. Here present 17 new TDEs that have been detected recently by the Zwicky Transient Facility along with Swift UV and X-ray follow-up observations. Our homogeneous analysis of the optical/UV light curves, including 22 previously known TDEs from the literature, reveals a clean separation of light curve properties with spectroscopic class. The TDEs with Bowen fluorescence features in their optical spectra have smaller blackbody radii, as well as longer rise times and higher disruption rates compared to the rest of the sample. The Bowen fluorescence mechanism requires a high density which can be reached at smaller radii, which in turn yields longer diffusion timescales. Thus, the difference in rise times suggests the pre-peak TDE light curves are governed not by the fallback timescale, but instead by the diffusion of photons through the tidal debris. The small subset of TDEs that show only helium emission lines in their spectra have the longest rise times, the highest luminosities and the lowest rates. We also report, for the first time, the detection of soft X-ray flares from a TDE on day timescales. Based on the fact the flares peak at a luminosity similar to the optical/UV blackbody luminosity, we attribute them to brief glimpses through a reprocessing layer that otherwise obscures the inner accretion flow.

We study the properties of the galaxies hosting the first 19 tidal disruption events (TDEs) detected with the Zwicky Transient Facility (ZTF) within the context of a carefully constructed, representative host galaxy sample. We find that the ZTF sample of TDE hosts is dominated by compact "green valley" galaxies. After we restrict the comparison sample to galaxies with a similar concentration, as measured by the Sersic index, we find this green valley overrepresentation is even larger. That is, concentrated red sequence galaxies are not producing TDEs at elevated levels. We present host galaxy spectra that show that E + A galaxies are overrepresented in the ZTF sample by a factor of 22, which is lower than previous TDE host galaxy studies have found. We find that this overrepresentation can be fully accounted for when taking into account the masses, colors, and Sersic indices of the ZTF TDE hosts. The combination of both green colors and high Sersic index of the typical TDE host galaxy could be explained if the TDE rate is temporarily enhanced following a merger that leads to a higher central concentration of stars.

Recent claimed detections of tidal disruption events (TDEs) in multi-wavelength data have opened potential new windows into the evolution and properties of otherwise dormant supermassive black holes (SMBHs) in the centres of galaxies. At present, there are several dozen TDE candidates, which share some properties and differ in others. The range in properties is broad enough to overlap other transient types, such as active galactic nuclei (AGN) and supernovae (SNe), which can make TDE classification ambiguous. A further complication is that "TDE signatures" have not been uniformly observed to similar sensitivities or even targeted across all candidates. This chapter both reviews those events that are unusual relative to other TDEs, including the possibility of TDEs in pre-existing AGN, and summarises those characteristics thought to best distinguish TDEs from continuously accreting AGN, strongly flaring AGN, SNe, and Gamma-Ray Bursts (GRBs), as well as other potential impostors like stellar collisions, "micro-TDEs," and circumbinary accretion flows. We conclude that multiple observables should be used to classify any one event as a TDE. We also consider the TDE candidate population as a whole, which, for certain host galaxy or SMBH characteristics, is distinguishable statistically from non-TDEs, suggesting that at least some TDE candidates do in fact arise from SMBH-disrupted stars.

We describe how the various outcomes of stellar tidal disruption give rise to observable radiation. We separately consider the cases where gas circularizes rapidly into an accretion disc, as well as the case when shocked debris streams provide the observable emission without having fully circularized. For the rapid circularization case, we describe how outflows, absorption by reprocessing layers, and Comptonization can cause the observed radiation to depart from that of a bare disc, possibly giving rise to the observed optical/UV emission along with soft X-rays from the disc. If, instead, most of the debris follows highly eccentric orbits for a significant time, many properties of the observed optical/UV emission can be explained by the scale of those eccentric orbits and the shocks embedded in the debris flow near orbital apocenter. In this picture, soft X-ray emission at early times results from the smaller amount of debris mass deflected into a compact accretion disc by weak shocks near the stellar pericenter. A general proposal for the near-constancy of the ultraviolet/optical color temperatures is provided, by linking it to incomplete thermalization of radiation in the atmosphere of the emitting region. We also briefly discuss the radio signals from the interaction of unbound debris and jets with the black hole environment.

We describe how the various outcomes of stellar tidal disruption give rise to observable radiation. We separately consider the cases where gas circularizes rapidly into an accretion disc, as well as the case when shocked debris streams provide the observable emission without having fully circularized. For the rapid circularization case, we describe how outflows, absorption by reprocessing layers, and Comptonization can cause the observed radiation to depart from that of a bare disc, possibly giving rise to the observed optical/UV emission along with soft X-rays from the disc. If, instead, most of the debris follows highly eccentric orbits for a significant time, many properties of the observed optical/UV emission can be explained by the scale of those eccentric orbits and the shocks embedded in the debris flow near orbital apocenter. In this picture, soft X-ray emission at early times results from the smaller amount of debris mass deflected into a compact accretion disc by weak shocks near the stellar pericenter. A general proposal for the near-constancy of the ultraviolet/optical color temperatures is provided, by linking it to incomplete thermalization of radiation in the atmosphere of the emitting region. We also briefly discuss the radio signals from the interaction of unbound debris and jets with the black hole environment.

The field of time-domain astrophysics has entered the era of Multi-messenger Astronomy (MMA). One key science goal for the next decade (and beyond) will be to characterize gravitational wave (GW) and neutrino sources using the next generation of Extremely Large Telescopes (ELTs). These studies will have a broad impact across astrophysics, informing our knowledge of the production and enrichment history of the heaviest chemical elements, constrain the dense matter equation of state, provide independent constraints on cosmology, increase our understanding of particle acceleration in shocks and jets, and study the lives of black holes in the universe. Future GW detectors will greatly improve their sensitivity during the coming decade, as will near-infrared telescopes capable of independently finding kilonovae from neutron star mergers. However, the electromagnetic counterparts to high-frequency (LIGO/Virgo band) GW sources will be distant and faint and thus demand ELT capabilities for characterization. ELTs will be important and necessary contributors to an advanced and complete multi-messenger network.

We present light curves and spectra of the tidal disruption event (TDE) ASASSN-18pg/AT 2018dyb spanning a period of one year. The event shows a plethora of strong emission lines, including the Balmer series, He ii, He i, and metal lines of O iii λ3760 and N iii λλ4100, 4640 (blended with He ii). The latter lines are consistent with originating from the Bowen fluorescence mechanism. By analyzing literature spectra of past events, we conclude that these lines are common in TDEs. The spectral diversity of optical TDEs is thus larger than previously thought and includes N-rich events besides H- and He-rich events. We study how the spectral lines evolve with time, by means of their width, relative strength, and velocity offsets. The velocity width of the lines starts at ∼13,000 km s−1 and decreases with time. The ratio of He ii to N iii increases with time. The same is true for ASASSN-14li, which has a very similar spectrum to AT 2018dyb but its lines are narrower by a factor of >2. We estimate a black hole mass of MBH = 3.3 − 2.0 + 5.0 × 10 6 M⊙ by using the M–σ relation. This is consistent with the black hole mass derived using the MOSFiT transient fitting code. The detection of strong Bowen lines in the optical spectrum is an indirect proof for extreme ultraviolet and (reprocessed) X-ray radiation and favors an accretion origin for the TDE optical luminosity. A model where photons escape after multiple scatterings through a super-Eddington thick disk and its optically thick wind, viewed at an angle close to the disk plane, is consistent with the observations.

We report the discovery of nonstellar hydrogen Balmer and metastable helium absorption lines accompanying a transient, high-velocity (0.05c) broad absorption line (BAL) system in the optical spectra of the tidal disruption event (TDE) AT2018zr (z = 0.071). In the Hubble Space Telescope UV spectra, absorption of high-and low-ionization lines is also present at this velocity, making AT2018zr resemble a low-ionization BALQSO. We conclude that these transient absorption features are more likely to arise in fast outflows produced by the TDE than absorbed by the unbound debris. In accordance with the outflow picture, we are able to reproduce the flat-topped Ha emission in a spherically expanding medium without invoking the typical prescription of an elliptical disk. We also report the appearance of narrow (similar to 1000 kms(-1)) N III lambda 4640, He II lambda 4686, H alpha, and H beta emission in the late-time optical spectra of AT2018zr, which may be a result of UV continuum hardening at late times, as observed by Swift. Including AT2018zr, we find a high association rate (three out of four) of BALs in the UV spectra of TDEs. This suggests that outflows may be ubiquitous among TDEs and less sensitive to viewing angle effects compared to QSO outflows.