We model a compact black hole-accretion disk system in the collapsar scenario with full transport, frequency dependent, general relativistic radiation magnetohydrodynamics. We examine whether or not winds from a collapsar disk can undergo rapid neutron capture (r-process) nucleosynthesis and significantly contribute to solar r-process abundances. We find the inclusion of accurate transport has significant effects on outflows, raising the electron fraction above $Y_e \sim 0.3$ and preventing third peak r-process material from being synthesized. We compare our model to semi-analytic expectations and argue that accurate neutrino transport and realistic initial and boundary conditions are required to capture the dynamics and nucleosynthetic outcome of a collapsar.

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 \`a 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 present a conservative numerical method for radiation magnetohydrodynamics with frequency-dependent full transport in stationary spacetimes. This method is stable and accurate for both large and small optical depths and radiation pressures. The radiation stress-energy tensor is evolved in flux-conservative form, and closed with a swarm of samples that each transport a multigroup representation of the invariant specific intensity along a null geodesic. In each zone, the enclosed samples are used to efficiently construct a Delaunay triangulation of the unit sphere in the comoving frame, which in turn is used to calculate the Eddington tensor, average source terms, and adaptively refine the sample swarm. Radiation four-fources are evaluated in the moment sector in a semi-implicit fashion. The radiative transfer equation is solved in invariant form deterministically for each sample. Since each sample carries a discrete representation of the full spectrum, the cost of evaluating the transport operator is independent of the number of frequency groups, representing a significant reduction of algorithmic complexity for transport in frequency dependent problems. The major approximation we make in this work is performing scattering in an angle-averaged way, with Compton scattering further approximated by the Kompaneets equation. Local adaptivity in samples also makes this scheme more amenable to nonuniform meshes than a traditional Monte Carlo method. We describe the method and present results on a suite of test problems. We find that MOCMC converges at least as $\sim N^{-1}$, rather than the canonical Monte Carlo $N^{-1/2}$, where $N$ is the number of samples per zone. On several problems we demonstrate substantial improvement over Eddington and M1 closures and gray opacities.

The 2017 detection of the inspiral and merger of two neutron stars in gravitational waves and gamma rays was accompanied by a quickly-reddening transient. Such a transient was predicted to occur following a rapid neutron capture (r-process) nucleosynthesis event, which synthesizes neutron-rich, radioactive nuclei and can take place in both dynamical ejecta and in the wind driven off the accretion torus formed after a neutron star merger. We present the first three-dimensional general relativistic, full transport neutrino radiation magnetohydrodynamics (GRRMHD) simulations of the black hole-accretion disk-wind system produced by the GW170817 merger. We show that the small but non-negligible optical depths lead to neutrino transport globally coupling the disk electron fraction, which we capture by solving the transport equation with a Monte Carlo method. The resulting absorption drives up the electron fraction in a structured, continuous outflow, with electron fraction as high as $Y_e\sim 0.4$ in the extreme polar region. We show via nuclear reaction network and radiative transfer calculations that nucleosynthesis in the disk wind will produce a blue kilonova.

Our understanding of the engines and progenitors of gamma-ray bursts has expanded through the ages as a broader set of diagnostics has allowed us to test our understanding of these objects. Here we review the history of the growth in our understanding, focusing on 3 leading engines and 9 potential progenitors. The first gravitational wave detection of a short burst is the latest in a series of breakthrough observations shaping this understanding and we study the importance of multi-diagnostic, multi-messenger observations on these engines and their progenitors. Our understanding based on a detailed study of nearby bursts can be applied to make predictions for the trends expected as we begin to observe high redshift bursts and we discuss these trends.

There are now $\sim$20 multi-dimensional core-collapse supernova (CCSN) simulations that explode. However, these simulations have explosion energies that are a few times $10^{50}$ erg, not $10^{51}$ erg. In this manuscript, we compare the inferred explosion energies of these simulations and observations of 38 SN~IIP. Assuming a log-normal distribution, the mean explosion energy for the observations is $\mu_{\rm obs} = -0.13\pm 0.05$ ($\log_{10}(E/10^{51}\, {\rm erg})$) and the width is $\sigma_{\rm obs} = 0.21^{+0.05}_{-0.04}$. Only three CCSN codes have sufficient simulations to compare with observations: CHIMERA, CoCoNuT-FMT, and FORNAX. Currently, FORNAX has the largest sample of simulations. The two-dimensional FORNAX simulations show a correlation between explosion energy and progenitor mass, ranging from linear to quadratic, $E_{\rm sim} \propto M^{1-2}$; this correlation is consistent with inferences from observations. In addition, we infer the ratio of the observed-to-simulated explosion energies, $\Delta=\log_{10}(E_{\rm obs}/E_{\rm sim})$. For the CHIMERA set, $\Delta=0.33\pm0.06$; for CoCoNuT-FMT, $\Delta=0.62\pm0.05$; for FORNAX2D, $\Delta=0.73\pm0.05$, and for FORNAX3D, $\Delta=0.95\pm0.06$. On average, the simulations are less energetic than inferred energies from observations ($\Delta \approx 0.7$), but we also note that the variation among the simulations (max($\Delta$)-min($\Delta$) $\approx 0.6$) is as large as this average offset. This suggests that further improvements to the simulations could resolve the discrepancy. Furthermore, both the simulations and the observations are heavily biased. In this preliminary comparison, we model these biases, but to more reliably compare the explosion energies, we recommend strategies to un-bias both the simulations and observations.

The 2017 detection of the in-spiral and merger of two neutron stars was a landmark discovery in astrophysics. We now know that such mergers are central engines of short gamma ray bursts and sites of r-process nucleosynthesis, where the heaviest elements in our universe are formed. In the coming years, we expect many more such mergers. Modeling such systems presents a significant computational challenge along with the observational one. To meet this challenge, we present $\nu\texttt{bhlight}$, a scheme for solving general relativistic magnetohydrodynamics with energy-dependent neutrino transport in full (3+1)-dimensions, facilitated by Monte Carlo methods. We present a suite of tests demonstrating the accuracy, efficacy, and necessity of our scheme. We demonstrate the potential of our scheme by running a sample calculation in a domain of interest---the dynamics and composition of the accretion disk formed by a binary neutron star merger.

This paper describes the design and implementation of our new multi-group, multi-dimensional radiation hydrodynamics (RHD) code Fornax and provides a suite of code tests to validate its application in a wide range of physical regimes. Instead of focusing exclusively on tests of neutrino radiation hydrodynamics relevant to the core-collapse supernova problem for which Fornax is primarily intended, we present here classical and rigorous demonstrations of code performance relevant to a broad range of multi-dimensional hydrodynamic and multi-group radiation hydrodynamic problems. Our code solves the comoving-frame radiation moment equations using the M1 closure, utilizes conservative high-order reconstruction, employs semi-explicit matter and radiation transport via a high-order time stepping scheme, and is suitable for application to a wide range of astrophysical problems. To this end, we first describe the philosophy, algorithms, and methodologies of Fornax and then perform numerous stringent code tests, that collectively and vigorously exercise the code, demonstrate the excellent numerical fidelity with which it captures the many physical effects of radiation hydrodynamics, and show excellent strong scaling well above 100k MPI tasks.

Most one-dimensional core-collapse simulations fail to explode, yet multi-dimensional simulations often explode. A dominant multi-dimensional effect aiding explosion is neutrino-driven convection. We incorporate a convection model in approximate one-dimensional core-collapse supernova (CCSN) simulations. This is the 1D+ method. This convection model lowers the neutrino luminosity required for explosion by 30%, similar to the reduction observed in multi-dimensional simulations. The model is based upon the global turbulence model of Mabanta & Murphy (2018) and models the mean-field turbulent flow of neutrino-driven convection. In this preliminary investigation, we use simple neutrino heating and cooling algorithms to compare the critical condition in the 1D+ simulations with the critical condition observed in two-dimensional simulations. Qualitatively, the critical conditions in the 1D+ and the two-dimensional simulations are similar. The assumptions in the convection model affect the radial profiles of density, entropy, and temperature, and comparisons with the profiles of three dimensional simulations will help to calibrate these assumptions. These 1D+ simulations are consistent with the profiles and explosion conditions of equivalent two-dimensional CCSN simulations but are ~100 times faster, and the 1D+ prescription has the potential to be ~100,000 faster than three-dimensional CCSN simulations. The 1D+ technique will be ideally suited to test the explodability of thousands of progenitor models.

We present axisymmetric two-temperature general relativistic radiation magnetohydrodynamic (GRRMHD) simulations of the inner region of the accretion flow onto the supermassive black hole M87. We address uncertainties from previous modeling efforts through inclusion of models for (1) self-consistent dissipative and Coulomb electron heating (2) radiation transport (3) frequency-dependent synchrotron emission, self-absorption, and Compton scattering. We adopt a distance $D=16.7$ Mpc, an observer angle $\theta = 20^{\circ}$, and consider black hole masses $M/M_{\odot} = (3.3\times10^{9}, 6.2\times10^{9})$ and spins $a_{\star} = (0.5, 0.9375)$ in a four-simulation suite. For each $(M, a_{\star})$, we identify the accretion rate that recovers the 230 GHz flux from VLBI measurements. We report on disk thermodynamics at these accretion rates ($\dot{M}/\dot{M}_{\mathrm{Edd}} \sim 10^{-5}$). The disk remains geometrically thick; cooling does not lead to a thin disk component. While electron heating is dominated by Coulomb rather than dissipation for $r \gtrsim 10 GM/c^2$, the accretion disk remains two-temperature. Radiative cooling of electrons is not negligible, especially for $r \lesssim 10 GM/c^2$. The Compton $y$ parameter is of order unity. We then compare derived and observed or inferred spectra, mm images, and jet powers. Simulations with $M/M_{\odot} = 3.3\times10^{9}$ are in conflict with observations. These simulations produce mm images that are too small, while the low-spin simulation also overproduces X-rays. For $M/M_{\odot} = 6.2\times10^{9}$, both simulations agree with constraints on radio/IR/X-ray fluxes and mm image sizes. Simulation jet power is a factor $10^2-10^3$ below inferred values, a possible consequence of the modest net magnetic flux in our models.

In this paper, we present the results of our three-dimensional, multi-group, multi-neutrino-species radiation/hydrodynamic simulation using the state-of-the-art code F{\sc{ornax}} of the terminal dynamics of the core of a non-rotating 16-M$_{\odot}$ stellar progenitor. The calculation incorporates redistribution by inelastic scattering, a correction for the effect of many-body interactions on the neutrino-nucleon scattering rates, approximate general relativity (including the effects of gravitational redshifts), velocity-dependent frequency advection, and an implementation of initial perturbations in the progenitor core. The model explodes within $\sim$100 milliseconds of bounce (near when the silicon-oxygen interface is accreted through the temporarily-stalled shock) and by the end of the simulation (here, $\sim$677 milliseconds after bounce) is accumulating explosion energy at a rate of $\sim$2.5$\times$10$^{50}$ ergs s$^{-1}$. The supernova explosion resembles an asymmetrical multi-plume structure, with one hemisphere predominating. The gravitational mass of the residual proto-neutron star at $\sim$677 milliseconds is $\sim$1.42 M$_{\odot}$. Even at the end of the simulation, explosion in most of the solid angle is accompanied by some accretion in an annular fraction at the wasp-like waist of the debris field. The ejecta electron fraction (Y$_e$) is distributed from $\sim$0.48 to $\sim$0.56, with most of the ejecta mass proton-rich. This may have implications for supernova nucleosynthesis, and could have a bearing on the p- and $\nu$p-processes and on the site of the first peak of the r-process. The ejecta spatial distributions of both Y$_e$ and mass density are predominantly in wide-angle plumes and large-scale structures, but are nevertheless quite patchy.

We present results of 2D axisymmetric core-collapse supernova simulations, employing the FORNAX code, of nine progenitor models spanning 12 to 25 M$_{\odot}$ and evolved over a 20,000-km grid. We find that four of the nine models explode with inelastic scattering off electrons and neutrons as well as the many-body correction to neutrino-nucleon scattering opacities. We show that these four models feature sharp Si-O interfaces in their density profiles, and that the corresponding dip in density reduces the accretion rate around the stalled shock and prompts explosion. The non-exploding models lack such a steep feature, suggesting that Si-O interface is one key to explosion. Furthermore, we show that all of the non-exploding models can be nudged to explosion with modest changes to macrophysical inputs, including moderate rotation and perturbations to infall velocities, as well as to microphysical inputs, including changes to neutrino-nucleon interaction rates, suggesting that all the models are perhaps close to criticality. Exploding models have energies of few $\times$10$^{50}$ ergs at the end of our simulation, and are rising, suggesting the need to continue these simulations over larger grids and for longer times to reproduce the energies seen in Nature. We find that the morphology of the explosion contributes to the explosion energy, with more isotropic ejecta producing larger explosion energies. However, we do not find evidence for the Lepton-number Emission Self-Sustained Asymmetry. Finally, we look at PNS properties and explore the role of dimension in our simulations. We find that convection in the proto-neutron star (PNS) produces larger PNS radii as well as greater "$\nu_\mu$" luminosities in 2D compared to 1D.

We explore with self-consistent 2D F{\sc{ornax}} simulations the dependence of the outcome of collapse on many-body corrections to neutrino-nucleon cross sections, the nucleon-nucleon bremsstrahlung rate, electron capture on heavy nuclei, pre-collapse seed perturbations, and inelastic neutrino-electron and neutrino-nucleon scattering. Importantly, proximity to criticality amplifies the role of even small changes in the neutrino-matter couplings, and such changes can together add to produce outsized effects. When close to the critical condition the cumulative result of a few small effects (including seeds) that individually have only modest consequence can convert an anemic into a robust explosion, or even a dud into a blast. Such sensitivity is not seen in one dimension and may explain the apparent heterogeneity in the outcomes of detailed simulations performed internationally. A natural conclusion is that the different groups collectively are closer to a realistic understanding of the mechanism of core-collapse supernovae than might have seemed apparent.

We present new 1D (spherical) and 2D (axisymmetric) simulations of electron-capture (EC) and low-mass iron-core-collapse supernovae (SN). We consider six progenitor models: the ECSN progenitor from Nomoto (1984, 1987); two ECSN-like low-mass low-metallicity iron core progenitors from Heger (private communication); and the 9-, 10-, and 11-$M_\odot$ (zero-age main sequence) progenitors from Sukhbold et al. (2016). We confirm that the ECSN and ESCN-like progenitors explode easily even in 1D with explosion energies of up to a 0.15 Bethes ($1 {\rm B} \equiv 10^{51}\ {\rm erg}$), and are a viable mechanism for the production of very low-mass neutron stars. However, the 9-, 10-, and 11-$M_\odot$ progenitors do not explode in 1D and are not even necessarily easier to explode than higher-mass progenitor stars in 2D. We study the effect of perturbations and of changes to the microphysics and we find that relatively small changes can result in qualitatively different outcomes, even in 1D, for models sufficiently close to the explosion threshold. Finally, we revisit the impact of convection below the protoneutron star (PNS) surface. We analyze, 1D and 2D evolutions of PNSs subject to the same boundary conditions. We find that the impact of PNS convection has been underestimated in previous studies and could result in an increase of the neutrino luminosity by up to factors of two.

The mechanism driving core-collapse supernovae is sensitive to the interplay between matter and neutrino radiation. However, neutrino radiation transport is very difficult to simulate, and several radiation transport methods of varying levels of approximation are available. We carefully compare for the first time in multiple spatial dimensions the discrete ordinates (DO) code of Nagakura, Yamada, and Sumiyoshi and the Monte Carlo (MC) code Sedonu, under the assumptions of a static fluid background, flat spacetime, elastic scattering, and full special relativity. We find remarkably good agreement in all spectral, angular, and fluid interaction quantities, lending confidence to both methods. The DO method excels in determining the heating and cooling rates in the optically thick region. The MC method predicts sharper angular features due to the effectively infinite angular resolution, but struggles to drive down noise in quantities where subtractive cancellation is prevalent, such as the net gain in the protoneutron star and off-diagonal components of the Eddington tensor. We also find that errors in the angular moments of the distribution functions induced by neglecting velocity dependence are sub-dominant to those from limited momentum-space resolution. We briefly compare directly computed second angular moments to those predicted by popular algebraic two-moment closures, and find that the errors from the approximate closures are comparable to the difference between the DO and MC methods. Included in this work is an improved Sedonu code, which now implements a fully special relativistic, time-independent version of the grid-agnostic Monte Carlo random walk approximation.

We present axisymmetric numerical simulations of radiatively inefficient accretion flows onto black holes combining general relativity, magnetohydrodynamics, self-consistent electron thermodynamics, and frequency-dependent radiation transport. We investigate a range of accretion rates up to $10^{-5} \dot{M}_{\mathrm{Edd}}$ onto a $10^8 M_{\odot}$ black hole with spin $a_{\star} = 0.5$. We report on averaged flow thermodynamics as a function of accretion rate. We present the spectra of outgoing radiation and find that it varies strongly with accretion rate, from synchrotron-dominated in the radio at low $\dot{M}$ to inverse Compton-dominated at our highest $\dot{M}$. In contrast to canonical analytic models, we find that by $\dot{M} \approx 10^{-5} \dot{M}_{\mathrm{Edd}}$, the flow approaches $\sim 1\%$ radiative efficiency, with much of the radiation due to inverse Compton scattering off Coulomb-heated electrons far from the black hole. These results have broad implications for modeling of accreting black holes across a large fraction of the accretion rates realized in observed systems.

We derive an integral condition for core-collapse supernova (CCSN) explosions and use it to construct a new diagnostic of explodability. The fundamental challenge in CCSN theory is to explain how a stalled accretion shock revives to explode a star. In this manuscript, we assume that the shock revival is initiated by the delayed-neutrino mechanism and derive an integral condition for spherically symmetric shock expansion, $v_s > 0$. One of the most useful one-dimensional explosion conditions is the neutrino luminosity and mass-accretion rate ($L_{\nu}-\dot{\mathcal{M}}$) critical curve. Below this curve, steady-state stalled solutions exist, but above this curve, there are no stalled solutions. Burrows & Goshy suggested that the solutions above this curve are dynamic and explosive. In this manuscript, we take one step closer to proving this supposition; we show that all steady solutions above this curve have $v_s > 0$. Assuming that these steady $v_s > 0$ solutions correspond to explosion, we present a new dimensionless integral condition for explosion, $\Psi > 0$. $\Psi$ roughly describes the balance between pressure and gravity, and we show that this parameter is equivalent to the $\tau$ condition used to infer the $L_{\nu}-\dot{\mathcal{M}}$ critical curve. The illuminating difference is that there is a direct relationship between $\Psi$ and $v_s$. Below the critical curve, $\Psi$ may be negative, positive, and zero, which corresponds to receding, expanding, and stalled-shock solutions. At the critical curve, the minimum $\Psi$ solution is zero; above the critical curve, $\Psi_{\rm min} > 0$, and all steady solutions have $v_s > 0$. Using one-dimensional simulations, we confirm our primary assumptions and verify that $\Psi_{\rm min} > 0$ is a reliable and accurate explosion diagnostic.

We perform the first self-consistent, time-dependent, multi-group calculations in two dimensions (2D) to address the consequences of using the ray-by-ray+ transport simplification in core-collapse supernova simulations. Such a dimensional reduction is employed by many researchers to facilitate their resource-intensive calculations. Our new code (F{\sc{ornax}}) implements multi-D transport, and can, by zeroing out transverse flux terms, emulate the ray-by-ray+ scheme. Using the same microphysics, initial models, resolution, and code, we compare the results of simulating 12-, 15-, 20-, and 25-M$_{\odot}$ progenitor models using these two transport methods. Our findings call into question the wisdom of the pervasive use of the ray-by-ray+ approach. Employing it leads to maximum post-bounce/pre-explosion shock radii that are almost universally larger by tens of kilometers than those derived using the more accurate scheme, typically leaving the post-bounce matter less bound and artificially more "explodable." In fact, for our 25-M$_{\odot}$ progenitor, the ray-by-ray+ model explodes, while the corresponding multi-D transport model does not. Therefore, in two dimensions the combination of ray-by-ray+ with the axial sloshing hydrodynamics that is a feature of 2D supernova dynamics can result in quantitatively, and perhaps qualitatively, incorrect results.

We present a model for the temporal variability of long gamma-ray bursts during the prompt phase (the highly variable first 100 seconds or so), in the context of a magnetically arrested disk (MAD) around a black hole. In this state, sufficient magnetic flux is held on to the black hole such that it stalls the accretion near the inner region of the disk. The system transitions in and out of the MAD state, which we relate to the variable luminosity of the GRB during the prompt phase, with a characteristic timescale defined by the free fall time in the region over which the accretion is arrested. We present simple analytic estimates of the relevant energetics and timescales, and compare them to gamma-ray burst observations. In particular, we show how this model can reproduce the characteristic one second time scale that emerges from various analyses of the prompt emission light curve. We also discuss how our model can accommodate the potentially physically important correlation between a burst quiescent time and the duration of its subsequent pulse (Ramirez-Ruiz & Merloni 2001).

We calculate the distance-dependent performance of a few representative terrestrial neutrino detectors in detecting and measuring the properties of the $\nu_e$ breakout burst light curve in a Galactic core-collapse supernova. The breakout burst is a signature phenomenon of core collapse and offers a probe into the stellar core through collapse and bounce. We examine cases of no neutrino oscillations and oscillations due to normal and inverted neutrino-mass hierarchies. For the normal hierarchy, other neutrino flavors emitted by the supernova overwhelm the $\nu_e$ signal, making a detection of the breakout burst difficult. For the inverted hierarchy (IH), some detectors at some distances should be able to see the $\nu_e$ breakout burst peak and measure its properties. For the IH, the maximum luminosity of the breakout burst can be measured at 10 kpc to accuracies of $\sim$30% for Hyper-Kamiokande (Hyper-K) and $\sim$60% for the Deep Underground Neutrino Experiment (DUNE). Super-Kamiokande (Super-K) and Jiangmen Underground Neutrino Observatory (JUNO) lack the mass needed to make an accurate measurement. For the IH, the time of the maximum luminosity of the breakout burst can be measured in Hyper-K to an accuracy of $\sim$3 ms at 7 kpc, in DUNE $\sim$2 ms at 4 kpc, and JUNO and Super-K can measure the time of maximum luminosity to an accuracy of $\sim$2 ms at 1 kpc. Detector backgrounds in IceCube render a measurement of the $\nu_e$ breakout burst unlikely. For the inverted hierarchy, a measurement of the maximum luminosity of the breakout burst could be used to differentiate between nuclear equations of state.

We present bhlight, a numerical scheme for solving the equations of general relativistic radiation magnetohydrodynamics (GRRMHD) using a direct Monte Carlo solution of the frequency-dependent radiative transport equation. bhlight is designed to evolve black hole accretion flows at intermediate accretion rate, in the regime between the classical radiatively efficient disk and the radiatively inefficient accretion flow (RIAF), in which global radiative effects play a sub-dominant but non-negligible role in disk dynamics. We describe the governing equations, numerical method, idiosyncrasies of our implementation, and a suite of test and convergence results. We also describe example applications to radiative Bondi accretion and to a slowly accreting Kerr black hole in axisymmetry.

We present new two-dimensional (2D) axisymmetric neutrino radiation/hydrodynamic models of core-collapse supernova (CCSN) cores. We use the CASTRO code, which incorporates truly multi-dimensional, multi-group, flux-limited diffusion (MGFLD) neutrino transport, including all relevant $\mathcal{O}(v/c)$ terms. Our main motivation for carrying out this study is to compare with recent 2D models produced by other groups who have obtained explosions for some progenitor stars and with recent 2D VULCAN results that did not incorporate $\mathcal{O}(v/c)$ terms. We follow the evolution of 12, 15, 20, and 25 solar-mass progenitors to approximately 600 milliseconds after bounce and do not obtain an explosion in any of these models. Though the reason for the qualitative disagreement among the groups engaged in CCSN modeling remains unclear, we speculate that the simplifying ``ray-by-ray' approach employed by all other groups may be compromising their results. We show that ``ray-by-ray' calculations greatly exaggerate the angular and temporal variations of the neutrino fluxes, which we argue are better captured by our multi-dimensional MGFLD approach. On the other hand, our 2D models also make approximations, making it difficult to draw definitive conclusions concerning the root of the differences between groups. We discuss some of the diagnostics often employed in the analyses of CCSN simulations and highlight the intimate relationship between the various explosion conditions that have been proposed. Finally, we explore the ingredients that may be missing in current calculations that may be important in reproducing the properties of the average CCSNe, should the delayed neutrino-heating mechanism be the correct mechanism of explosion.

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_{\theta} + K_{\phi}). 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.

We present the galaxy two-point angular correlation function for galaxies selected from the seventh data release of the Sloan Digital Sky Survey. The galaxy sample was selected with $r$-band apparent magnitudes between 17 and 21; and we measure the correlation function for the full sample as well as for the four magnitude ranges: 17-18, 18-19, 19-20, and 20-21. We update the flag criteria to select a clean galaxy catalog and detail specific tests that we perform to characterize systematic effects, including the effects of seeing, Galactic extinction, and the overall survey uniformity. Notably, we find that optimally we can use observed regions with seeing $< 1\farcs5$, and $r$-band extinction < 0.13 magnitudes, smaller than previously published results. Furthermore, we confirm that the uniformity of the SDSS photometry is minimally affected by the stripe geometry. We find that, overall, the two-point angular correlation function can be described by a power law, $\omega(\theta) = A_\omega \theta^{(1-\gamma)}$ with $\gamma \simeq 1.72$, over the range $0\fdg005$--$10\degr$. We also find similar relationships for the four magnitude subsamples, but the amplitude within the same angular interval for the four subsamples is found to decrease with fainter magnitudes, in agreement with previous results. We find that the systematic signals are well below the galaxy angular correlation function for angles less than approximately $5\degr$, which limits the modeling of galaxy angular correlations on larger scales. Finally, we present our custom, highly parallelized two-point correlation code that we used in this analysis.

We present a formulation for multigroup radiation hydrodynamics that is correct to order $O(v/c)$ using the comoving-frame approach and the flux-limited diffusion approximation. We describe a numerical algorithm for solving the system, implemented in the compressible astrophysics code, CASTRO. CASTRO uses an Eulerian grid with block-structured adaptive mesh refinement based on a nested hierarchy of logically-rectangular variable-sized grids with simultaneous refinement in both space and time. In our multigroup radiation solver, the system is split into three parts, one part that couples the radiation and fluid in a hyperbolic subsystem, another part that advects the radiation in frequency space, and a parabolic part that evolves radiation diffusion and source-sink terms. The hyperbolic subsystem and the frequency space advection are solved explicitly with high-order Godunov schemes, whereas the parabolic part is solved implicitly with a first-order backward Euler method. Our multigroup radiation solver works for both neutrino and photon radiation.

The multidimensional character of the hydrodynamics in core-collapse supernova (CCSN) cores is a key facilitator of 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 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 to 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. 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, 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. On the basis of our study, we hypothesize that neutrino-driven buoyant convection should almost always dominate the SASI when the supernova explosion is neutrino-driven.

Monte Carlo approaches to radiation transport have several attractive properties such as simplicity of implementation, high accuracy, and good parallel scaling. Moreover, Monte Carlo methods can handle complicated geometries and are relatively easy to extend to multiple spatial dimensions, which makes them potentially interesting in modeling complex multi-dimensional astrophysical phenomena such as core-collapse supernovae. The aim of this paper is to explore Monte Carlo methods for modeling neutrino transport in core-collapse supernovae. We generalize the Implicit Monte Carlo photon transport scheme of Fleck & Cummings and gray discrete-diffusion scheme of Densmore et al. to energy-, time-, and velocity-dependent neutrino transport. Using our 1D spherically-symmetric implementation, we show that, similar to the photon transport case, the implicit scheme enables significantly larger timesteps compared with explicit time discretization, without sacrificing accuracy, while the discrete-diffusion method leads to significant speed-ups at high optical depth. Our results suggest that a combination of spectral, velocity-dependent, Implicit Monte Carlo and discrete-diffusion Monte Carlo methods represents a robust approach for use in neutrino transport calculations in core-collapse supernovae. Our velocity-dependent scheme can easily be adapted to photon transport.

In accretion-based models for Sgr A* the X-ray, infrared, and millimeter emission arise in a hot, geometrically thick accretion flow close to the black hole. The spectrum and size of the source depend on the black hole mass accretion rate $\dot{M}$. Since Gillessen et al. have recently discovered a cloud moving toward Sgr A* that will arrive in summer 2013, $\dot{M}$ may increase from its present value $\dot{M}_0$. We therefore reconsider the "best-bet" accretion model of Moscibrodzka et al., which is based on a general relativistic MHD flow model and fully relativistic radiative transfer, for a range of $\dot{M}$. We find that for modest increases in $\dot{M}$ the characteristic ring of emission due to the photon orbit becomes brighter, more extended, and easier to detect by the planned Event Horizon Telescope submm VLBI experiment. If $\dot{M} \gtrsim 8 \dot{M}_0$ this "silhouette of the black hole will be hidden beneath the synchrotron photosphere at 230 GHz, and for $\dot{M} \gtrsim 16 \dot{M}_0$ the silhouette is hidden at 345 GHz. We also find that for $\dot{M} > 2 \dot{M}_0$ the near-horizon accretion flow becomes a persistent X-ray and mid-infrared source, and in the near-infrared Sgr A* will acquire a persistent component that is brighter than currently observed flares.

We report transient quasi-periodic oscillations (QPOs) on minute timescales in relativistic, radiative models of the galactic center source Sgr A*. The QPOs result from nonaxisymmetric $m=1$ structure in the accretion flow excited by MHD turbulence. Near-infrared (NIR) and X-ray power spectra show significant peaks at frequencies comparable to the orbital frequency at the innermost stable circular orbit (ISCO) $f_o$. The excess power is associated with inward propagating magnetic filaments inside the ISCO. The amplitudes of the QPOs are sensitive to the electron distribution function. We argue that transient QPOs appear at a range of frequencies in the neighborhood of $f_o$ and that the power spectra, averaged over long times, likely show a broad bump near $f_o$ rather than distinct, narrow QPO features.

Global, general relativistic magnetohydrodynamic (GRMHD) simulations of nonradiative, magnetized disks are widely used to model accreting black holes. We have performed a convergence study of GRMHD models computed with HARM3D. The models span a factor of 4 in linear resolution, from 96x96x64 to 384x384x256. We consider three diagnostics of convergence: (1) dimensionless shell-averaged quantities such as plasma \beta; (2) the azimuthal correlation length of fluid variables; and (3) synthetic spectra of the source including synchrotron emission, absorption, and Compton scattering. Shell-averaged temperature is, except for the lowest resolution run, nearly independent of resolution; shell-averaged plasma \beta\ decreases steadily with resolution but shows signs of convergence. The azimuthal correlation lengths of density, internal energy, and temperature decrease steadily with resolution but show signs of convergence. In contrast, the azimuthal correlation length of magnetic field decreases nearly linearly with grid size. We argue by analogy with local models, however, that convergence should be achieved with another factor of 2 in resolution. Synthetic spectra are, except for the lowest resolution run, nearly independent of resolution. The convergence behavior is consistent with that of higher physical resolution local model (shearing box) calculations and with the recent nonrelativistic global convergence studies of Hawley et al. (2011).

Electron-positron pairs may be produced near accreting black holes by a variety of physical processes, and the resulting pair plasma may be accelerated and collimated into a relativistic jet. Here we use a self-consistent dynamical and radiative model to investigate pair production by \gamma\gamma collisions in weakly radiative accretion flows around a black hole of mass M and accretion rate \dot{M}. Our flow model is drawn from general relativistic magnetohydrodynamic simulations, and our radiation field is computed by a Monte Carlo transport scheme assuming the electron distribution function is thermal. We argue that the pair production rate scales as r^{-6} M^{-1} \dot{M}^{6}. We confirm this numerically and calibrate the scaling relation. This relation is self-consistent in a wedge in M, \dot{M} parameter space. If \dot{M} is too low the implied pair density over the poles of the black hole is below the Goldreich-Julian density and \gamma\gamma pair production is relatively unimportant; if \dot{M} is too high the models are radiatively efficient. We also argue that for a power-law spectrum the pair production rate should scale with the observables L_X \equiv X-ray luminosity and M as L_X^2 M^{-4}. We confirm this numerically and argue that this relation likely holds even for radiatively efficient flows. The pair production rates are sensitive to black hole spin and to the ion-electron temperature ratio which are fixed in this exploratory calculation. We finish with a brief discussion of the implications for Sgr A* and M87.

We review results from general relativistic axisymmetric magnetohydrodynamic simulations of accretion in Sgr A*. We use general relativistic radiative transfer methods and to produce a broad band (from millimeter to gamma-rays) spectrum. Using a ray tracing scheme we also model images of Sgr A* and compare the size of image to the VLBI observations at 230 GHz. We perform a parameter survey and study radiative properties of the flow models for various black hole spins, ion to electron temperature ratios, and inclinations. We scale our models to reconstruct the flux and the spectral slope around 230 GHz. The combination of Monte Carlo spectral energy distribution calculations and 230 GHz image modeling constrains the parameter space of the numerical models. Our models suggest rather high black hole spin ($a_*\approx 0.9$), electron temperatures close to the ion temperature ($T_i/T_e \sim 3$) and high inclination angles ($i \approx 90 \deg$).

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 Ti/Te. We find that: (1) models with Ti/Te=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 Ti/Te, but the Ti/Te = 10 models are generally large and at most marginally consistent with the 230 GHz VLBI data; (5) for models with Ti/Te=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* \sim 0.9. We finish with a discussion of the limitations of our model and prospects for future improvements.

We describe a Monte Carlo radiative transport code intended for calculating spectra of hot, optically thin plasmas in full general relativity. The version we describe here is designed to model hot accretion flows in the Kerr metric and therefore incorporates synchrotron emission and absorption, and Compton scattering. The code can be readily generalized, however, to account for other radiative processes and an arbitrary spacetime. We describe a suite of test problems, and demonstrate the expected $N^{-1/2}$ convergence rate, where $N$ is the number of Monte Carlo samples. Finally we illustrate the capabilities of the code with a model calculation, a spectrum of the slowly accreting black hole Sgr A* based on data provided by a numerical general relativistic MHD model of the accreting plasma.

The ultracompact binary systems V407 Vul (RX J1914.4+2456) and HM Cnc (RX J0806.3+1527) - a two-member subclass of the AM CVn stars - continue to pique interest because they defy unambiguous classification. Three proposed models remain viable at this time, but none of the three is significantly more compelling than the remaining two, and all three can satisfy the observational constraints if parameters in the models are tuned. One of the three proposed models is the direct impact model of Marsh & Steeghs (2002), in which the accretion stream impacts the surface of a rapidly-rotating primary white dwarf directly but at a near-glancing angle. One requirement of this model is that the accretion stream have a high enough density to advect its specific kinetic energy below the photosphere for progressively more-thermalized emission downstream, a constraint that requires an accretion spot size of roughly 1.2x10^5 km^2 or smaller. Having at hand a smoothed particle hydrodynamics code optimized for cataclysmic variable accretion disk simulations, it was relatively straightforward for us to adapt it to calculate the footprint of the accretion stream at the nominal radius of the primary white dwarf, and thus to test this constraint of the direct impact model. We find that the mass flux at the impact spot can be approximated by a bivariate Gaussian with standard deviation \sigma_{\phi} = 164 km in the orbital plane and \sigma_{\theta} = 23 km in the perpendicular direction. The area of the the 2\sigma ellipse into which 86% of the mass flux occurs is roughly 47,400 km^2, or roughly half the size estimated by Marsh & Steeghs (2002). We discuss the necessary parameters of a simple model of the luminosity distribution in the post-impact emission region.

We have discovered a spectral transition from the low/hard (LH) state to the high/soft (HS) state when Aquila X-1 was {\em declining} in observations made with the {\it Rossi X-Ray Timing Explorer (RXTE)}. The 2--200 keV energy flux corresponding to the state transition is $1.1\times{10}^{-9} ergs cm^{-2} s^{-1}$, an order of magnitude lower than observed in the past. The 2--200 keV peak flux of the following HS state is $1.6\times{10}^{-9} ergs cm^{-2} s^{-1}$. The relation between the luminosity of the hard-to-soft state transition and the peak luminosity of the following HS state confirms the linear relation found previously. This implies that the luminosity of the hard-to-soft state transition is not determined solely by the mass accretion rate, but appears to be determined by the peak luminosity of the soft X-ray outburst. We also found that the LH-to-HS state transition occurred at a luminosity similar to that of the corresponding HS-to-LH state transition, i.e., there is no apparent hysteresis. These results provide additional evidence that the mass in the accretion disk affects the luminosity of the hard-to-soft state transition, and that the accretion flow that powers the LH state is related to the accretion flow that powers the HS state at a later time.

We present the software tool FITDisk, a precompiled-binary Windows GUI version of our smoothed particle hydrodynamics cataclysmic variable accretion disk research code. Cataclysmic variables are binary star systems in which a compact stellar remnant, typically a white dwarf star, is stripping mass from a lower-main-sequence companion star by way of an accretion disk. Typically the disk is the brightest component of the system, because the plasma is heated dramatically as it spirals down in the gravitational well of the primary white dwarf star. The shortest-period systems can display disk "superhump" oscillations driven by the rotating tidal field of the secondary star. FITDisk models these accretion disk phenomena using a fully three-dimensional hydrodynamics calculation, and data can be visualized as they are computed or stored to hard drive for later playback at a fast frame rate. Simulations are visualized using OpenGL graphics and the viewing angle can be changed interactively. Pseudo light curves of simulated systems can be plotted along with the associated Fourier amplitude spectrum. FITDisk is available for free download at www.astro.fit.edu/cv/fitdisk.html.