Cosmic rays (CRs) with ~GeV energies can contribute significantly to the energy and pressure budget in the interstellar, circumgalactic, and intergalactic medium (ISM, CGM, IGM). Recent cosmological simulations have begun to explore these effects, but almost all studies have been restricted to simplified models with constant CR diffusivity and/or streaming speeds. Physical models of CR propagation/scattering via extrinsic turbulence and self-excited waves predict transport coefficients which are complicated functions of local plasma properties. In a companion paper, we consider a wide range of observational constraints to identify proposed physically-motivated cosmic-ray propagation scalings which satisfy both detailed Milky Way (MW) and extra-galactic $\gamma$-ray constraints. Here, we compare the effects of these models relative to simpler 'diffusion+streaming' models on galaxy and CGM properties at dwarf through MW mass scales. The physical models predict large local variations in CR diffusivity, with median diffusivity increasing with galacto-centric radii and decreasing with galaxy mass and redshift. These effects lead to a more rapid dropoff of CR energy density in the CGM (compared to simpler models), in turn producing weaker effects of CRs on galaxy star formation rates (SFRs), CGM absorption profiles and galactic outflows. The predictions of the more physical CR models tend to lie 'in between' models which ignore CRs entirely and models which treat CRs with constant diffusivity.

Project AMIGA (Absorption Maps In the Gas of Andromeda) is a large ultraviolet Hubble Space Telescope program, which has assembled a sample of 43 QSOs that pierce the circumgalactic medium (CGM) of Andromeda (M31) from R=25 to 569 kpc (25 of them probing gas from 25 kpc to about the virial radius-Rvir = 300 kpc-of M31). Our large sample provides an unparalleled look at the physical conditions and distribution of metals in the CGM of a single galaxy using ions that probe a wide range of gas phases (Si II, Si III, Si IV, C II, C IV, and O VI, the latter being from the Far Ultraviolet Spectroscopic Explorer). We find that Si III and O VI have near unity covering factor maintained all the way out to 1.2Rvir and 1.9Rvir, respectively. We show that Si III is the dominant ion over Si II and Si IV at any R. While we do not find that the properties of the CGM of M31 depend strongly on the azimuth, we show that they change remarkably around 0.3-0.5Rvir, conveying that the inner regions of the CGM of M31 are more dynamic and have more complicated multi-phase gas-structures than at R>0.5Rvir. We estimate the metal mass of the CGM within Rvir as probed by Si II, Si III, and Si IV is 2x10^7 Msun and by O VI is >8x10^7 Msun, while the baryon mass of the 10^4-10^5.5 K gas is ~4x10^10 (Z/0.3 Zsun)^(-1) Msun within Rvir. We show that different zoom-in cosmological simulations of L* galaxies better reproduce the column density profile of O VI with R than Si III or the other studied ions. We find that observations of the M31 CGM and zoom-in simulations of L* galaxies have both lower ions showing higher column density dispersion and dependence on R than higher ions, indicating that the higher ionization structures are larger and/or more broadly distributed.

The microphysics of ~GeV cosmic ray (CR) transport on galactic scales remain deeply uncertain, with almost all studies adopting simple prescriptions (e.g. constant-diffusivity). We explore different physically-motivated, anisotropic, dynamical CR transport scalings in high-resolution cosmological FIRE simulations of dwarf and ~$L_{\ast}$ galaxies where scattering rates vary with local plasma properties motivated by extrinsic turbulence (ET) or self-confinement (SC) scenarios, with varying assumptions about e.g. turbulent power spectra on un-resolved scales, Alfven-wave damping, etc. We self-consistently predict observables including $\gamma$-rays ($L_{\gamma}$), grammage, residence times, and CR energy densities to constrain the models. We demonstrate many non-linear dynamical effects (not captured in simpler models) tend to enhance confinement. For example, in multi-phase media, even allowing arbitrary fast transport in neutral gas does not substantially reduce CR residence times (or $L_{\gamma}$), as transport is rate-limited by the ionized WIM and 'inner CGM' gaseous halo ($10^{4}-10^{6}$ K gas within 10-30 kpc), and $L_{\gamma}$ can be dominated by trapping in small 'patches.' Most physical ET models contribute negligible scattering of ~1-10 GeV CRs, but it is crucial to account for anisotropy and damping (especially of fast modes) or else scattering rates would violate observations. We show that the most widely-assumed scalings for SC models produce excessive confinement by factors >100 in the WIM and inner CGM, where turbulent and Landau damping dominate. This suggests either a breakdown of quasi-linear theory used to derive the CR transport parameters in SC, or that other novel damping mechanisms dominate in intermediate-density ionized gas.

We present a suite of high-resolution cosmological simulations, using the FIRE-2 feedback physics together with explicit treatment of magnetic fields, anisotropic conduction and viscosity, and cosmic rays (CRs) injected by supernovae (including anisotropic diffusion, streaming, adiabatic, hadronic and Coulomb losses). We survey systems from ultra-faint dwarf ($M_{\ast}\sim 10^{4}\,M_{\odot}$, $M_{\rm halo}\sim 10^{9}\,M_{\odot}$) through Milky Way masses, systematically vary CR parameters (e.g. the diffusion coefficient $\kappa$ and streaming velocity), and study an ensemble of galaxy properties (masses, star formation histories, mass profiles, phase structure, morphologies). We confirm previous conclusions that magnetic fields, conduction, and viscosity on resolved ($\gtrsim 1\,$pc) scales have small effects on bulk galaxy properties. CRs have relatively weak effects on all galaxy properties studied in dwarfs ($M_{\ast} \ll 10^{10}\,M_{\odot}$, $M_{\rm halo} \lesssim 10^{11}\,M_{\odot}$), or at high redshifts ($z\gtrsim 1-2$), for any physically-reasonable parameters. However at higher masses ($M_{\rm halo} \gtrsim 10^{11}\,M_{\odot}$) and $z\lesssim 1-2$, CRs can suppress star formation by factors $\sim 2-4$, given relatively high effective diffusion coefficients $\kappa \gtrsim 3\times10^{29}\,{\rm cm^{2}\,s^{-1}}$. At lower $\kappa$, CRs take too long to escape dense star-forming gas and lose energy to hadronic collisions, producing negligible effects on galaxies and violating empirical constraints from $\gamma$-ray emission. But around $\kappa\sim 3\times10^{29}\,{\rm cm^{2}\,s^{-1}}$, CRs escape the galaxy and build up a CR-pressure-dominated halo which supports dense, cool ($T\ll 10^{6}$ K) gas that would otherwise rain onto the galaxy. CR heating (from collisional and streaming losses) is never dominant.

We study the effects of cosmic rays (CRs) on outflows from star-forming galaxies in the circum and inter-galactic medium (CGM/IGM), in high-resolution, fully-cosmological FIRE-2 simulations (accounting for mechanical and radiative stellar feedback, magnetic fields, anisotropic conduction/viscosity/CR diffusion and streaming, and CR losses). We showed previously that massive ($M_{\rm halo}\gtrsim 10^{11}\,M_{\odot}$), low-redshift ($z\lesssim 1-2$) halos can have CR pressure dominate over thermal CGM pressure and balance gravity, giving rise to a cooler CGM with an equilibrium density profile. This dramatically alters outflows. Absent CRs, high gas thermal pressure in massive halos "traps" galactic outflows near the disk, so they recycle. With CRs injected in supernovae as modeled here, the low-pressure halo allows "escape" and CR pressure gradients continuously accelerate this material well into the IGM in "fast" outflows, while lower-density gas at large radii is accelerated in-situ into "slow" outflows that extend to $>$Mpc scales. CGM/IGM outflow morphologies are radically altered: they become mostly volume-filling (with inflow in a thin mid-plane layer) and coherently biconical from the disk to $>$Mpc. The CR-driven outflows are primarily cool ($T\sim10^{5}\,$K) and low-velocity. All of these effects weaken and eventually vanish at lower halo masses ($\lesssim 10^{11}\,M_{\odot}$) or higher redshifts ($z\gtrsim 1-2$), reflecting the ratio of CR to thermal+gravitational pressure in the outer halo. We present a simple analytic model which explains all of the above phenomena.

We use a well-known deep neural network framework, called Mask R-CNN, for identification of solar filaments in full-disk H-alpha images from Big Bear Solar Observatory (BBSO). The image data, collected from BBSO's archive, are integrated with the spatiotemporal metadata of filaments retrieved from the Heliophysics Events Knowledgebase (HEK) system. This integrated data is then treated as the ground-truth in the training process of the model. The available spatial metadata are the output of a currently running filament-detection module developed and maintained by the Feature Finding Team; an international consortium selected by NASA. Despite the known challenges in the identification and characterization of filaments by the existing module, which in turn are inherited into any other module that intends to learn from such outputs, Mask R-CNN shows promising results. Trained and validated on two years worth of BBSO data, this model is then tested on the three following years. Our case-by-case and overall analyses show that Mask R-CNN can clearly compete with the existing module and in some cases even perform better. Several cases of false positives and false negatives, that are correctly segmented by this model are also shown. The overall advantages of using the proposed model are two-fold: First, deep neural networks' performance generally improves as more annotated data, or better annotations are provided. Second, such a model can be scaled up to detect other solar events, as well as a single multi-purpose module. The results presented in this study introduce a proof of concept in benefits of employing deep neural networks for detection of solar events, and in particular, filaments.

We investigate the impact of cosmic rays (CRs) on the circumgalactic medium (CGM) in FIRE-2 simulations, for ultra-faint dwarf through Milky Way (MW)-mass halos hosting star-forming (SF) galaxies. Our CR treatment includes injection by supernovae, anisotropic streaming and diffusion along magnetic field lines, collisional and streaming losses, with constant parallel diffusivity $\kappa\sim3\times10^{29}\,\mathrm{cm^2\ s^{-1}}$ chosen to match $\gamma$-ray observations. With this, CRs become more important at larger halo masses and lower redshifts, and dominate the pressure in the CGM in MW-mass halos at $z\lesssim 1-2$. The gas in these ``CR-dominated'' halos differs significantly from runs without CRs: the gas is primarily cool (a few $\sim10^{4}\,$K), and the cool phase is volume-filling and has a thermal pressure below that needed for virial or local thermal pressure balance. Ionization of the ``low'' and ``mid'' ions in this diffuse cool gas is dominated by photo-ionization, with O VI columns $\gtrsim 10^{14.5}\,\mathrm{cm^{-2}}$ at distances $\gtrsim 150\,\mathrm{kpc}$. CR and thermal gas pressure are locally anti-correlated, maintaining total pressure balance, and the CGM gas density profile is determined by the balance of CR pressure gradients and gravity. Neglecting CRs, the same halos are primarily warm/hot ($T\gtrsim 10^{5}\,$K) with thermal pressure balancing gravity, collisional ionization dominates, O VI columns are lower and Ne VIII higher, and the cool phase is confined to dense filaments in local thermal pressure equilibrium with the hot phase.

Radiative turbulent mixing layers should be ubiquitous in multi-phase gas with shear flow. They are a potentially attractive explanation for the high ions such as OVI seen in high velocity clouds and the circumgalactic medium (CGM) of galaxies. We perform 3D MHD simulations with non-equilibrium (NEI) and photoionization modeling, with an eye towards testing simple analytic models. Even purely hydrodynamic collisional ionization equilibrium (CIE) calculations have column densities much lower than observations. Characteristic inflow and turbulent velocities are much less than the shear velocity, and the layer width $h \propto t_\mathrm{cool}^{1/2}$ rather than $h \propto t_\mathrm{cool}$. Column densities are not independent of density or metallicity as analytic scalings predict, and show surprisingly weak dependence on shear velocity and density contrast. Radiative cooling, rather than Kelvin-Helmholtz instability, appears paramount in determining the saturated state. Low pressure due to fast cooling both seeds turbulence and sets the entrainment rate of hot gas, whose enthalpy flux, along with turbulent dissipation, energizes the layer. Regardless of initial geometry, magnetic fields are amplified and stabilize the mixing layer via magnetic tension, producing almost laminar flow and depressing column densities. NEI effects can boost column densities by factors of a few. Suppression of cooling by NEI or photoionization can in principle also increase OVI column densities, but in practice is unimportant for CGM conditions. To explain observations, sightlines must pierce hundreds or thousands of mixing layers, which may be plausible if the CGM exists as a `fog' of tiny cloudlets.

This study is motivated by recent observations on ubiquitous interstellar density filaments and guided by modern theories of compressible magnetohydrodynamic (MHD) turbulence. The interstellar turbulence shapes the observed density structures. As the fundamental dynamics of compressible MHD turbulence, perpendicular turbulent mixing of density fluctuations entails elongated density structures aligned with the local magnetic field, accounting for low-density parallel filaments seen in diffuse atomic and molecular gas. The elongation of low-density parallel filaments depends on the turbulence anisotropy. When taking into account the partial ionization, we find that the minimum width of parallel filaments in the cold neutral medium and molecular clouds is determined by the neutral-ion decoupling scale perpendicular to magnetic field. In highly supersonic MHD turbulence in molecular clouds, both low-density parallel filaments due to anisotropic turbulent mixing and high-density filaments due to shock compression exist.

Cold ($T\sim 10^{4} \ \mathrm{K}$) gas is very commonly found in both galactic and cluster halos. There is no clear consensus on its origin. Such gas could be uplifted from the central galaxy by galactic or AGN winds. Alternatively, it could form in situ by thermal instability. Fragmentation into a multi-phase medium has previously been shown in hydrodynamic simulations to take place once $t_\mathrm{cool}/t_\mathrm{ff}$, the ratio of the cooling time to the free-fall time, falls below a threshold value. Here, we use 3D plane-parallel MHD simulations to investigate the influence of magnetic fields. We find that because magnetic tension suppresses buoyant oscillations of condensing gas, it destabilizes all scales below $l_\mathrm{A}^\mathrm{cool} \sim v_\mathrm{A} t_\mathrm{cool}$, enhancing thermal instability. This effect is surprisingly independent of magnetic field orientation or cooling curve shape, and sets in even at very low magnetic field strengths. Magnetic fields critically modify both the amplitude and morphology of thermal instability, with $\delta \rho/\rho \propto \beta^{-1/2}$, where $\beta$ is the ratio of thermal to magnetic pressure. In galactic halos, magnetic fields can render gas throughout the entire halo thermally unstable, and may be an attractive explanation for the ubiquity of cold gas, even in the halos of passive, quenched galaxies.

We investigate the apparent power-law scaling of the pseudo phase space density (PPSD) in CDM halos. We study fluid collapse, using the close analogy between the gas entropy and the PPSD in the fluid approximation. Our hydrodynamic calculations allow for a precise evaluation of logarithmic derivatives. For scale-free initial conditions, entropy is a power law in Lagrangian (mass) coordinates, but not in Eulerian (radial) coordinates. The deviation from a radial power law arises from incomplete hydrostatic equilibrium (HSE), linked to bulk inflow and mass accretion, and the convergence to the asymptotic central power-law slope is very slow. For more realistic collapse, entropy is not a power law with either radius or mass due to deviations from HSE and scale-dependent initial conditions. Instead, it is a slowly rolling power law that appears approximately linear on a log-log plot. Our fluid calculations recover PPSD power-law slopes and residual amplitudes similar to N-body simulations, indicating that deviations from a power law are not numerical artefacts. In addition, we find that realistic collapse is not self-similar: scale lengths such as the shock radius and the turnaround radius are not power-law functions of time. We therefore argue that the apparent power-law PPSD cannot be used to make detailed dynamical inferences or extrapolate halo profiles inward, and that it does not indicate any hidden integrals of motion. We also suggest that the apparent agreement between the PPSD and the asymptotic Bertschinger slope is purely coincidental.

Increasing observational and theoretical evidence points to binary white dwarf mergers as the origin of some if not most normal Type Ia supernovae (SNe Ia). In this paper, we discuss the post-merger evolution of binary white dwarf (WD) mergers, and their relevance to the double-degenerate channel of SNe Ia. We present 3D simulations of carbon-oxygen (C/O) WD binary systems undergoing unstable mass transfer, varying both the total mass and the mass ratio. We demonstrate that these systems generally give rise to a one-armed gravitational spiral instability. The spiral density modes transport mass and angular momentum in the disk even in the absence of a magnetic field, and are most pronounced for secondary-to-primary mass ratios larger than $0.6$. We further analyze carbon burning in these systems to assess the possibility of detonation. Unlike the case of a $1.1 + 1.0 M_{\odot}$ C/O WD binary, we find that WD binary systems with lower mass and smaller mass ratios do not detonate as SNe Ia up to $\sim8-22$ outer dynamical times. Two additional models do however undergo net heating, and their secular increase in temperature could possibly result in a detonation on timescales longer than those considered here.

Turbulent dynamo field amplification has often been invoked to explain the strong field strengths in thin rims in supernova shocks ($\sim 100 \, \mu$G) and in radio relics in galaxy clusters ($\sim \mu$G). We present high resolution MHD simulations of the interaction between pre-shock turbulence, clumping and shocks, to quantify the conditions under which turbulent dynamo amplification can be significant. We demonstrate numerically converged field amplification which scales with Alfv\'en Mach number, $B/B_0 \propto {\mathcal M}_{\rm A}$, up to ${\mathcal M}_{\rm A} \sim 150$. This implies that the post-shock field strength is relatively independent of the seed field. Amplification is dominated by compression at low ${\mathcal M}_{\rm A}$, and stretching (turbulent amplification) at high ${\mathcal M}_{\rm A}$. For high $\mathcal{M}_{\rm A}$, the $B$-field grows exponentially and saturates at equipartition with turbulence, while the vorticity jumps sharply at the shock and subsequently decays; the resulting field is orientated predominately along the shock normal (an effect only apparent in 3D and not 2D). This agrees with the radial field bias seen in supernova remnants. By contrast, for low $\mathcal{M}_{\rm A}$, field amplification is mostly compressional, relatively modest, and results in a predominantly perpendicular field. The latter is consistent with the polarization seen in radio relics. Our results are relatively robust to the assumed level of gas clumping. Our results imply that the turbulent dynamo may be important for supernovae, but is only consistent with the field strength, and not geometry, for cluster radio relics. For the latter, this implies strong pre-existing $B$-fields in the ambient cluster outskirts.

Double-degenerate (DD) mergers of carbon-oxygen white dwarfs have recently emerged as a leading candidate for normal Type Ia supernovae (SNe Ia). However, many outstanding questions surround DD mergers, including the characteristics of their light curves and spectra. We have recently identified a spiral instability in the post-merger phase of DD mergers and demonstrated that this instability self-consistently leads to detonation in some cases. We call this the spiral merger SN Ia model. Here, we utilize the SuperNu radiative transfer software to calculate three-dimensional synthetic light curves and spectra of the spiral merger simulation with a system mass of 2.1 $M_\odot$ from Kashyap et al. Because of their large system masses, both violent and spiral merger light curves are slowly declining. The spiral merger resembles very slowly declining SNe Ia, including SN 2001ay, and provides a more natural explanation for its observed properties than other SN Ia explosion models. Previous synthetic light curves and spectra of violent DD mergers demonstrate a strong dependence on viewing angle, which is in conflict with observations. Here, we demonstrate that the light curves and spectra of the spiral merger are less sensitive to the viewing angle than violent mergers, in closer agreement with observation. We find that the spatial distribution of 56Ni and IMEs follows a characteristic hourglass shape. We discuss the implications of the asymmetric distribution of 56Ni for the early-time gamma-ray observations of 56Ni from SN 2014J. We suggest that DD mergers that agree with the light curves and spectra of normal SNe Ia will likely require a lower system mass.

Thermonuclear, or Type Ia supernovae (SNe Ia), originate from the explosion of carbon--oxygen white dwarfs, and serve as standardizable cosmological candles. However, despite their importance, the nature of the progenitor systems that give rise to SNe Ia has not been hitherto elucidated. Observational evidence favors the double-degenerate channel in which merging white dwarf binaries lead to SNe Ia. Furthermore, significant discrepancies exist between observations and theory, and to date, there has been no self-consistent merger model that yields a SNe Ia. Here we show that a spiral mode instability in the accretion disk formed during a binary white dwarf merger leads to a detonation on a dynamical timescale. This mechanism sheds light on how white dwarf mergers may frequently yield SNe Ia.

Type Ia supernovae (SNe Ia) play a crucial role as standardizable cosmological candles, though the nature of their progenitors is a subject of active investigation. Recent observational and theoretical work has pointed to merging white dwarf binaries, referred to as the double-degenerate channel, as the possible progenitor systems for some SNe Ia. Additionally, recent theoretical work suggests that mergers which fail to detonate may produce magnetized, rapidly-rotating white dwarfs. In this paper, we present the first multidimensional simulations of the post-merger evolution of white dwarf binaries to include the effect of the magnetic field. In these systems, the two white dwarfs complete a final merger on a dynamical timescale, and are tidally disrupted, producing a rapidly-rotating white dwarf merger surrounded by a hot corona and a thick, differentially-rotating disk. The disk is strongly susceptible to the magnetorotational instability (MRI), and we demonstrate that this leads to the rapid growth of an initially dynamically weak magnetic field in the disk, the spin-down of the white dwarf merger, and to the subsequent central ignition of the white dwarf merger. Additionally, these magnetized models exhibit new features not present in prior hydrodynamic studies of white dwarf mergers, including the development of MRI turbulence in the hot disk, magnetized outflows carrying a significant fraction of the disk mass, and the magnetization of the white dwarf merger to field strengths $\sim 2 \times 10^8$ G. We discuss the impact of our findings on the origins, circumstellar media, and observed properties of SNe Ia and magnetized white dwarfs.