Using flow models based on axisymmetric general relativistic magnetohydrodynamics (GRMHD) simulations, we construct radiative models for Sgr A*. Spectral energy distributions that include the effects of thermal synchrotron emission and absorption, and Compton scattering, are calculated using a Monte Carlo technique. Images are calculated using a ray-tracing scheme. All models are scaled so that the 230 GHz flux density is 3.4 Jy. The key model parameters are the dimensionless black hole spin a * , the inclination i, and the ion-to-electron temperature ratio T i /T e . We find that: (1) models with T i /T e = 1 are inconsistent with the observed submillimeter spectral slope; (2) the X-ray flux is a strongly increasing function of a * ; (3) the X-ray flux is a strongly increasing function of i; (4) 230 GHz image size is a complicated function of i, a * , and T i /T e , but the T i /T e = 10 models are generally large and at most marginally consistent with the 230 GHz VLBI data; (5) for models with T i /T e = 10 and i = 85 deg the event horizon is cloaked behind a synchrotron photosphere at 230 GHz and will not be seen by VLBI, but these models overproduce NIR and X-ray flux; (6) in all models whose SEDs are consistent with observations the event horizon is uncloaked at 230 GHz; (7) the models that are most consistent with the observations have a * ∼ 0.9. We finish with a discussion of the limitations of our model and prospects for future improvements.

Multi-dimensional instabilities have become an important ingredient in core-collapse supernova (CCSN) theory. Therefore, it is necessary to understand the driving mechanism of the dominant instability. We compare our parameterized three-dimensional CCSN simulations with other buoyancy-driven simulations and propose scaling relations for neutrino-driven convection. Through these comparisons, we infer that buoyancy-driven convection dominates post-shock turbulence in our simulations. In support of this inference, we present four major results. First, the convective fluxes and kinetic energies in the neutrino-heated region are consistent with expectations of buoyancy-driven convection. Second, the convective flux is positive where buoyancy actively drives convection, and the radial and tangential components of the kinetic energy are in rough equipartition (i.e. K r ∼ K θ + K φ ). Both results are natural consequences of buoyancy-driven convection, and are commonly observed in simulations of convection. Third, buoyant driving is balanced by turbulent dissipation. Fourth, the convective luminosity and turbulent dissipation scale with the driving neutrino power. In all, these four results suggest that in neutrino-driven explosions, the multi-dimensional motions are consistent with neutrino-driven convection.

The multidimensional character of the hydrodynamics in core-collapse supernova (CCSN) cores is a key facilitator of explosions. A major goal over the last decade has been to elucidate how multidimensionality produces more favorable conditions for explosions. Unfortunately, much of this work has necessarily been performed assuming axisymmetry and it remains unclear whether or not this compromises those results. In this work, we present analyses of simplified two-and three-dimensional CCSN models with the goal of comparing the multidimensional hydrodynamics in setups that differ only in dimension. Not surprisingly, we find many differences between 2D and 3D models. While some differences are subtle and perhaps not crucial to understanding the explosion mechanism, others are quite dramatic and make interpreting 2D CCSN models problematic. In particular, we find that imposing axisymmetry artificially produces excess power at the largest spatial scales, power that has been deemed critical in the success of previous explosion models and has been attributed solely to the standing accretion shock instability. Nevertheless, our 3D models, which have an order of magnitude less power on large scales compared to 2D models, explode earlier. Since we see explosions earlier in 3D than in 2D, the vigorous sloshing associated with the large scale power in 2D models is either not critical in any dimension or the explosion mechanism operates differently in 2D and 3D. Possibly related to the earlier explosions in 3D, we find that about 25% of the accreted material spends more time in the gain region in 3D than in 2D, being exposed to more integrated heating and reaching higher peak entropies, an effect we associate with the differing characters of turbulence in 2D and 3D. Finally, we discuss a simple model for the runaway growth of buoyant bubbles that is able to quantitatively account for the growth of the shock radius and predicts a critical luminosity relation.

We have conducted nineteen state-of-the-art 3D core-collapse supernova simulations spanning a broad range of progenitor masses. This is the largest collection of sophisticated 3D supernova simulations ever performed. We have found that while the majority of these models explode, not all do, and that even models in the middle of the available progenitor mass range may be less explodable. This does not mean that those models for which we did not witness explosion would not explode in Nature, but that they are less prone to explosion than others. One consequence is that the "compactness" measure is not a metric for explodability. We find that lower-mass massive star progenitors likely experience lower-energy explosions, while the higher-mass massive stars likely experience higher-energy explosions. Moreover, most 3D explosions have a dominant dipole morphology, have a pinched, wasp-waist structure, and experience simultaneous accretion and explosion. We reproduce the general range of residual neutron-star masses inferred for the galactic neutron-star population. The most massive progenitor models, however, in particular visà vis explosion energy, need to be continued for longer physical times to asymptote to their final states. We find that while the majority of the inner ejecta have Y e = 0.5, there is a substantial proton-rich tail. This result has important implications for the nucleosynthetic yields as a function of progenitor. Finally, we find that the non-exploding models eventually evolve into compact inner configurations that experience a quasi-periodic spiral SASI mode. We otherwise see little evidence of the SASI in the exploding models.

We investigate the structure of the stalled supernova shock in both 2D and 3D and explore the differences in the effects of neutrino heating and the standing accretion shock instability (SASI). We find that early on the amplitude of the dipolar mode of the shock is factors of ∼2−3 smaller in 3D than in 2D. However, later in both 3D and 2D the monopole and dipole modes start to grow until explosion. Whereas in 2D the (l, m) = (1, 0) mode changes sign quasi-periodically, producing the "up-and-down" motion always seen in modern 2D simulations, in 3D this almost never happens. Rather, in 3D when the dipolar mode starts to grow, it grows in magnitude and wanders stochastically in direction until settling before explosion to a particular patch of solid angle. Furthermore, in 2D we find that the amplitude of the dipolar shock deformation separates into two classes. For the first, identified with the SASI and for a wide range of "low" neutrino luminosities, this amplitude remains small and roughly constant. For the other, identified with higher luminosities and neutrino-driven convection, the dipolar amplitude grows sharply. Importantly, it is only for this higher luminosity class that we see neutrino-driven explosions within ∼1 second of bounce. Moreover, for the "low" luminosity runs (including zero), the power spectra of these dipolar oscillations peak in the 30-50 Hz range associated with advection timescales, while for the high-luminosity runs the power spectra at lower frequencies are significantly more prominent. We associate this enhanced power at lower frequencies with slower convective effects and the secular growth of the dipolar shock amplitude. Though our study involves a simplified, parametrized approach, on the basis of it we hypothesize that neutrino-driven buoyant convection should almost always dominate the SASI when the supernova explosion is neutrino-driven.