We present the first-ever moving-mesh general relativistic hydrodynamics solver for static spacetimes as implemented in the code, MANGA. Our implementation builds on the architectures of MANGA and the numerical relativity Python package NRPy+. We review the general algorithm to solve these equations and, in particular, detail the time stepping; Riemann solution across moving faces; conversion between primitive and conservative variables; validation and correction of hydrodynamic variables; and mapping of the metric to a Voronoi moving-mesh grid. We present test results for the numerical integration of an unmagnetized Tolman-Oppenheimer-Volkoff star for 24 dynamical times. We demonstrate that at a resolution of $10^6$ mesh generating points, the star is stable and its central density drifts downward by 2% over this timescale. At a lower resolution the central density drift increases in a manner consistent with the adopted second order spatial reconstruction scheme. These results agree well with the exact solutions, and we find the error behavior to be similar to Eulerian codes with second-order spatial reconstruction. We also demonstrate that the new code recovers the fundamental mode frequency for the same TOV star but with its initial pressure depleted by 10%

It is known that the Higgs potential in the Standard Model can drive a successful inflation as long as the Higgs field couples non-minimally to gravity. It is then inevitable to take into account the loop corrections of the Standard Model particles to the Higgs potential in the Higgs inflation. In this paper, we discuss the one-loop corrections at finite temperature to the curvature perturbation generated during the Higgs inflation. We find that the thermal loop effects can suppress the power of the curvature perturbation at large scales, thus resulting in a low quadrupole of the temperature anisotropy in the cosmic microwave background.

A compact object moving relative to surrounding gas accretes material and perturbs the density of gas in its vicinity. In the classical picture of Bondi-Hoyle-Lyttleton accretion, the perturbation takes the form of an overdense wake behind the object, which exerts a dynamical friction drag. We use hydrodynamic simulations to investigate how the accretion rate and strength of dynamical friction are modified by the presence of outflow from the compact object. We show that the destruction of the wake by an outflow reduces dynamical friction, and reverses its sign when the outflow is strong enough, in good quantitative agreement with analytic calculations. For a strong isotropic outflow, the outcome on scales that we have simulated is a negative dynamical friction, i.e., net acceleration. For jet-like outflows driven by reprocessed accretion, both the rate of accretion and the magnitude of dynamical friction drop for more powerful jets. The accretion rate is strongly intermittent when the jet points to the same direction as the motion of the compact object. The dynamical effects of outflows may be important for the evolution of compact objects during the common envelope phase of binary systems, and for accreting compact objects and massive stars encountering AGN discs.

Valuable information about the neutron star interior can be obtained by comparing observations of thermal radiation from a cooling neutron star crust with theoretical models. Nuclear burning of lighter elements that diffuse to deeper layers of the envelope can alter the relation between surface and interior temperatures and can change the chemical composition over time. We calculate new temperature relations and consider two effects of diffusive nuclear burning (DNB) for H-C envelopes. First, we consider the effect of a changing envelope composition and find that hydrogen is consumed on short timescales and our temperature evolution simulations correspond to those of a hydrogen-poor envelope within ~100 days. The transition from a hydrogen-rich to a hydrogen-poor envelope is potentially observable in accreting NS systems as an additional initial decline in surface temperature at early times after the outburst. Second, we find that DNB can produce a non-negligible heat flux, such that the total luminosity can be dominated by DNB in the envelope rather than heat from the deep interior. However, without continual accretion, heating by DNB in H-C envelopes is only relevant for <1-80 days after the end of an accretion outburst, as the amount of light elements is rapidly depleted. Comparison to crust cooling data shows that DNB does not remove the need for an additional shallow heating source. We conclude that solving the time-dependent equations of the burning region in the envelope self-consistently in thermal evolution models instead of using static temperature relations would be valuable in future cooling studies.

We study the longitudinal stability of beam-plasma systems in the presence of a density inhomogeneity in the background plasma. Previous works have focused on the non-relativistic regime where hydrodynamical models are used to evolve pre-existing Langmuir waves within inhomogeneous background plasmas. Here, for the first time we study the problem with kinetic equations in a fully-relativistic way. We do not assume the existence of Langmuir waves, and we focus on the rate and the mechanism by which waves are excited in such systems from an initial perturbation. We derive the structure of the unstable modes and compute an analytical approximation for their growth rates. Our computation is limited to dilute and cold beams, and shows an excellent agreement with particle-in-cell simulations performed using the SHARP code. We show that, due to such an inhomogeneity, the virulent beam-plasma instabilities in the intergalactic medium are not suppressed but their counterparts in the solar wind can be suppressed as evidenced by propagating type-III solar radio bursts.

We describe the structure and implementation of a radiation hydrodynamic solver for MANGA, the moving-mesh hydrodynamics module of the large-scale parallel code, Charm N-body GrAvity solver (ChaNGa). We solve the equations of time dependent radiative transfer using a reduced speed of light approximation following the algorithm of Jiang et al (2014). By writing the radiative transfer equations as a generalized conservation equation, we solve the transport part of these equations on an unstructured Voronoi mesh. We then solve the source part of the radiative transfer equations following Jiang et al (2014) using an implicit solver, and couple this to the hydrodynamic equations. The use of an implicit solver ensure reliable convergence and preserves the conservation properties of these equations even in situations where the source terms are stiff due to the small coupling timescales between radiation and matter. We present the results of a limited number of test cases (energy conservation, momentum conservation, dynamic diffusion, linear waves, crossing beams, and multiple shadows) to show convergence with analytic results and numerical stability. We also show that it produces qualitatively the correct results in the presence of multiple sources in the optically thin case.

Pair creation on the cosmic infrared background and subsequent inverse-Compton scattering on the CMB potentially reprocesses the TeV emission of blazars into faint GeV halos with structures sensitive to intergalactic magnetic fields (IGMF). Previous work has shown that these halos are then highly-anisotropic and extended. If the coherence length of the IGMF is greater than the inverse-Compton cooling length of the pairs, then the orientation of the gamma-ray halo will be correlated with the direction of the magnetic field which is unknown and expected to change for each source. In order to constructively add each source we then use angular power spectra which are insensitive to the jet orientation. By looking at known GeV blazars detected by Fermi, we exclude the existence of an IGMF with coherence lengths $>100$ Mpc at greater than $3.9\sigma$ with current-day strengths in the range $10^{-16}$ to $10^{-15}$ G, and at 2$\sigma$ from $10^{-17}$ to $10^{-14}$ G. This provides a direct measurement of the non-existence of gamma-ray halos, providing an important check on previous results.

We employ the earlier published proper motions of the newly discovered Antlia 2 dwarf galaxy derived from Gaia data to calculate its orbital distribution in the cosmologically recent past. Using these observationally motivated orbits, we calculate the effect of the Antlia 2 dwarf galaxy on the outer HI disk of the Milky Way, using both test particle and Smoothed Particle Hydrodynamics simulations. We find that orbits with low pericenters, $\sim$ 10 kpc, produce disturbances that match the observed outer HI disk perturbations. We have independently recalculated the proper motion of the Antlia 2 dwarf from Gaia data and found a proper motion of $(\mu_{\alpha}cos\delta, \mu_{\delta}) = (-0.068,0.032) \pm (0.023,-0.031)~\rm mas/yr$, which agrees with results from Torrealba et al. (2019) within the errors, but gives lower mean pericenters, e.g., $\sim$ 15 kpc for our fiducial model of the Milky Way. We also show that the Sagittarius dwarf galaxy interaction does not match the observed perturbations in the outer gas disk. Thus, Antlia 2 may be the driver of the observed large perturbations in the outer gas disk of the Galaxy. The current location of the Antlia 2 dwarf galaxy closely matches that predicted by an earlier dynamical analysis (Chakrabarti \& Blitz 2009) of the dwarf that drove ripples in the outer Galaxy, and, in particular, its orbit is nearly coplanar to the Galactic disk. If the Antlia 2 dwarf galaxy is responsible for the perturbations in the outer Galactic disk, it would have a specific range of proper motions that we predict here; this can be tested soon with Gaia DR-3 and Gaia DR-4 data.

We outline the methodology of simulating common envelope evolution (CEE) with the moving-mesh code MANGA. We extend MANGA to include multiple time-steps. This provides substantial speedups for problems with large dynamic range. We describe the implementation of realistic equations of state relevant in stellar structure and the generation of suitable initial conditions. We then carry out two example simulations of a 2 $\rm{M}_{\odot}$ red giant with a 0.36 $\rm{M}_{\odot}$ core and a 1 $\rm{M}_{\odot}$ companion undergoing CEE for 240 days. In one simulation the red giant is set into corotation with the orbital motion and in the other it is non-rotating. We find that the separation between the companion and red giant core shrinks from 52 $\rm{R}_{\odot}$ to 3.6 $\rm{R}_{\odot}$ and 3.2 $\rm{R}_{\odot}$ respectively, ending with an eccentricity of 0.1. We also find that 66 and 63 per cent of the envelope mass is ejected. This is higher than in many previous works. Several reasons for this are discussed. These include our inclusion of recombination energy. Our simulations show that putting giants in corotation increases the fraction of mass ejected from the system and results in a larger final orbital separation. We conclude that the entire envelope of the red giant might be ejected during the plunge phase of CEE in this region of parameter space.

The study of how neutron stars cool over time can provide invaluable insights into fundamental physics such as the nuclear equation of state and superconductivity and superfluidity. A critical relation in neutron star cooling is the one between observed surface temperature and interior temperature. This relation is determined by the composition of the neutron star envelope and can be influenced by the process of diffusive nuclear burning (DNB). We calculate models of envelopes that include DNB and find that DNB can lead to a rapidly changing envelope composition which can be relevant for understanding the long-term cooling behavior of neutron stars. We also report on analysis of the latest temperature measurements of the young neutron star in the Cassiopeia A supernova remnant. The 13 Chandra observations over 18 years show that the neutron star's temperature is decreasing at a rate of 2-3 percent per decade, and this rapid cooling can be explained by the presence of a proton superconductor and neutron superfluid in the core of the star.

The identification of the electromagnetic counterpart of the gravitational wave event, GW170817, and discovery of neutrinos and gamma-rays from TXS 0506+056 heralded the new era of multi-messenger astrophysics. As the number of multi-messenger events rapidly grow over the next decade, the cyberinfrastructure requirements to handle the increase in data rates, data volume, need for event follow up, and analysis across the different messengers will also explosively grow. The cyberinfrastructure requirements to enhance multi-messenger astrophysics will both be a major challenge and opportunity for astronomers, physicists, computer scientists and cyberinfrastructure specialists. Here we outline some of these requirements and argue for a distributed cyberinfrastructure institute for multi-messenger astrophysics to meet these challenges.

A critical relation in the study of neutron star cooling is the one between surface temperature and interior temperature. This relation is determined by the composition of the neutron star envelope and can be affected by the process of diffusive nuclear burning (DNB), which occurs when elements diffuse to depths where the density and temperature are sufficiently high to ignite nuclear burning. We calculate models of H-He and He-C envelopes that include DNB and obtain analytic temperature relations that can be used in neutron star cooling simulations. We find that DNB can lead to a rapidly changing envelope composition and prevents the build-up of thermally stable hydrogen columns y$_H$ > 10$^{7}$ g cm$^{-2}$, while DNB can make helium envelopes more transparent to heat flux for surface temperatures $T_s$ > 2 $\times 10^6$ K. We perform neutron star cooling simulations in which we evolve temperature and envelope composition, with the latter due to DNB and accretion from the interstellar medium. We find that a time-dependent envelope composition can be relevant for understanding the long-term cooling behaviour of isolated neutron stars. We also report on the latest Chandra observations of the young neutron star in the Cassiopeia A supernova remnant; the resulting 13 temperature measurements over more than 18 years yield a ten-year cooling rate of $\approx$ 2%. Finally, we fit the observed cooling trend of the Cassiopeia A neutron star with a model that includes DNB in the envelope.

An intergalactic magnetic field stronger than $3\times10^{-13}$~G would explain the lack of a bright, extended degree-scale, GeV-energy inverse Compton component in the gamma-ray spectra of TeV-blazars. A robustly predicted consequence of the presence of such a field is the existence of degree-scale GeV-energy gamma-ray halos -- gamma-ray bow ties -- about TeV-bright active galactic nuclei, corresponding to more than half of all radio galaxies. However, the emitting regions of these halos are confined to and aligned with the direction of the relativistic jets associated with gamma-ray sources. Based on the orientation of radio jets, we align and stack corresponding degree-scale gamma-ray images of isolated Fanaroff-Riley class I and II objects and exclude the existence of these halos at overwhelming confidence, limiting the intergalactic field strength to $<10^{-15}$~G for large-scale fields and progressively larger in the diffusive regime when the correlation length of the field becomes small in comparison to 1 Mpc. When combined with prior limits on the strength of the intergalactic magnetic field, this excludes a purely magnetic explanation for the absence of halos. Thus, it requires the existence of novel physical processes that preempt the creation of halos, e.g., the presence of beam-plasma instabilities in the intergalactic medium or a drastic cutoff of the very high energy spectrum of these sources.

We extend the work of Chang & Chakrabarti (2011) to find simple scaling relations between the density response of the gas disk of a spiral galaxy and the pericenter distance and mass ratio of a perturbing satellite. From the analysis of results from a test particle code, we obtained a simple scaling relation for the density response due to a single satellite interacting with a galactic disk, over a wide range of satellite masses and pericenter distances. We have also explored the effects of multiple satellites on the galactic disk, focusing on cases that are commonly found in cosmological simulations. Here, we use orbits for the satellites that are drawn from cosmological simulations. For these cases, we compare our approximate scaling relations to the density response generated by satellites, and find that for two satellite interactions, our scaling relations approximately recover the response of the galactic disk. We have also examined the observed HI data in the outskirts of several spiral galaxies from the THINGS sample and compared the observed perturbations to that of cosmological simulations and our own scaling relations. While small perturbations can be excited by satellites drawn from cosmological simulations, we find that large perturbations (such as those that are seen in some THINGS galaxies like M51) are not recovered by satellites drawn from cosmological simulations that are similar to Milky Way galaxies.

We explore how inhomogeneity in the background plasma number density alters the growth of electrostatic unstable wavemodes of beam plasma systems. This is particularly interesting for blazar-driven beam-plasma instabilities, which may be suppressed by inhomogeneities in the intergalactic medium as was recently claimed in the literature. Using high resolution Particle-In-Cell simulations with the SHARP code, we show that the growth of the instability is local, i.e., regions with almost homogeneous background density will support the growth of the Langmuir waves, even when they are separated by strongly inhomogeneous regions, resulting in an overall slower growth of the instability. We also show that if the background density is continuously varying, the growth rate of the instability is lower; though in all cases, the system remains within the linear regime longer and the instability is not extinguished. In all cases, the beam loses approximately the same fraction of its initial kinetic energy in comparison to the uniform case at non-linear saturation. Thus, inhomogeneities in the intergalactic medium are unlikely to suppress the growth of blazar-driven beam-plasma instabilities.

In principle, the angular anisotropy in the extragalactic gamma-ray background (EGRB) places severe constraints upon putative populations of unresolved gamma-ray point sources. Existing estimates of the EGRB anisotropy have been constructed by excising known point sources, e.g., taken from the First or 2 Year Fermi-LAT Source Catalog (1FGL or 2FGL, respectively) and statistically analyzing the residual gamma-ray sky maps. We perform an independent check of the EGRB anisotropy limits by comparing the values obtained from the 1FGL-masked sky maps to the signal implied by sources that lie below the 1FGL detection threshold in the more sensitive 2FGL and 1FHL (First Fermi-LAT catalog of >10 GeV sources). As such, our analysis provides an internal consistency check of implications for source counts and spectral index distributions of gamma-ray bright active galactic nuclei obtained from Fermi-LAT data. Based on this, we find evidence for substantially larger anisotropies than those previously reported at energies above 5 GeV, where BL Lac objects are likely to provide the bulk of their contribution to the EGRB. This uncertainty in the EGRB anisotropy cautions against using it as an independent constraint for the high-redshift gamma-ray universe. Moreover, this would suggest that contrary to previous claims, smooth extensions of the resolved point-source population may be able to simultaneously explain both the isotropic and anisotropic components of the EGRB.

We study the temporal and spatial distribution of star formation rates in four well-studied star-forming regions in local molecular clouds(MCs): Taurus, Perseus, $\rho$ Ophiuchi, and Orion A. Using published mass and age estimates for young stellar objects in each system, we show that the rate of star formation over the last 10 Myrs has been accelerating and is (roughly) consistent with a $t^2$ power law. This is in line with previous studies of the star formation history of molecular clouds and with recent theoretical studies. We further study the clustering of star formation in the Orion Nebula Cluster(ONC). We examine the distribution of young stellar objects as a function of their age by computing an effective half-light radius for these young stars subdivided into age bins. We show that the distribution of young stellar objects is broadly consistent with the star formation being entirely localized within the central region. We also find a slow radial expansion of the newly formed stars at a velocity of $v=0.17\,{\rm km\,s}^{-1}$, which is roughly the sound speed of the cold molecular gas. This strongly suggests the dense structures that form stars persist much longer than the local dynamical time. We argue that this structure is quasi-static in nature and is likely the result of the density profile approaching an attractor solution as suggested by recent analytic and numerical analysis.

We present results of hydrodynamic simulations of massive star forming regions with and without protostellar jets. We show that jets change the normalization of the stellar mass accretion rate, but do not strongly affect the dynamics of star formation. In particular, $M_*(t) \propto f^2 (t-t_*)^2$ where $f = 1 - f_{\rm jet}$ is the fraction of mass accreted onto the protostar, $f_{\rm jet}$ is the fraction ejected by the jet, and $(t-t_*)^2$ is the time elapsed since the formation of the first star. The star formation efficiency is nonlinear in time. We find that jets have only a small effect (of order 25\%) on the accretion rate onto the protostellar disk (the "raw" accretion rate). We show that the small scale structure -- the radial density, velocity, and mass accretion profiles are very similar in the jet and no-jet cases. Finally, we show that the inclusion of jets does drive turbulence but only on small (parsec) scales.

We investigate the contribution of outer HI disks to the observable population of merging black hole binaries. Like dwarf galaxies, the outer HI disks of spirals have low star formation rates and lower metallicities than the inner disks of spirals. Since low-metallicity star formation can produce more detectable compact binaries than typical star formation, the environments in the outskirts of spiral galaxies may be conducive to producing a rich population of massive binary black holes. We consider here both detailed controlled simulations of spirals and cosmological simulations, as well as the current range of observed values for metallicity and star formation in outer disks. We find that outer HI disks contribute at least as much as dwarf galaxies do to the observed LIGO/Virgo detection rates. Identifying the host galaxies of merging massive black holes should provide constraints on cosmological parameters and insights into the formation channels of binary mergers.

The light curve of GW170817 is surprisingly blue and bright. Assuming that the event is a binary neutron star merger, we argue that blueness and brightness of the light curve is the result of ejecta that contains an substantial amount of thermal energy. To achieve this, the ejecta must be reheated at a substantial distance (1 to 2000 solar radii) from the merger to avoid losing the energy to adiabatic cooling. We show that this reheating can occur if the merger occurs in a hierarchical triple system where the outer star has evolved and filled its Roche lobe. The outer star feeds mass to the inner binary, forming a circumbinary disc, driving the inner binary to merge. Because the outer star fills its Roche lobe, a substantial fraction of the dynamical ejecta collides with the evolved star, reheating the ejecta in the process. We suggest that the process of mass transfer in heirarchical triples tends to form coplanar triple systems such as PSR J0337+1715, and may provide electromagnetic counterparts to binary black hole mergers.

Many astrophysical plasmas are prone to beam-plasma instabilities. For relativistic and dilute beams, the {\it spectral} support of the beam-plasma instabilities is narrow, i.e., the linearly unstable modes that grow with rates comparable to the maximum growth rate occupy a narrow range of wave numbers. This places stringent requirements on the box-sizes when simulating the evolution of the instabilities. We identify the implied lower limits on the box size imposed by the longitudinal beam plasma instability, i.e., typically the most stringent condition required to correctly capture the linear evolution of the instabilities in multidimensional simulations. We find that sizes many orders of magnitude larger than the resonant wavelength are typically required. Using one-dimensional particle-in-cell simulations, we show that the failure to sufficiently resolve the spectral support of the longitudinal instability yields slower growth and lower levels of saturation, potentially leading to erroneous physical conclusion.

We describe the structure and implementation of a moving-mesh hydrodynamics solver in the large-scale parallel code, Charm N-body GrAvity solver (ChaNGa). While largely based on the algorithm described by Springel (2010) that is implemented in AREPO, our algorithm differs a few aspects. We describe our use of the Voronoi tessellation library, VORO++, to compute the Voronoi tessellation directly. We also incorporate some recent advances in gradient estimation and reconstruction that gives better accuracy in hydrodynamic solutions at minimal computational cost. We validate this module with a small battery of test problems against the smooth particle hydrodynamics solver included in ChaNGa. Finally, we study one example of a scientific problem involving the mergers of two main sequence stars and highlight the small quantitative differences between smooth particle and moving-mesh hydrodynamics. We close with a discussion of anticipated future improvements and advancements.

Numerical heating in particle-in-cell (PIC) codes currently precludes the accurate simulation of cold, relativistic plasma over long periods, severely limiting their applications in astrophysical environments. We present a spatially higher-order accurate relativistic PIC algorithm in one spatial dimension, which conserves charge and momentum exactly. We utilize the smoothness implied by the usage of higher-order interpolation functions to achieve a spatially higher-order accurate algorithm (up to fifth order). We validate our algorithm against several test problems -- thermal stability of stationary plasma, stability of linear plasma waves, and two-stream instability in the relativistic and non-relativistic regimes. Comparing our simulations to exact solutions of the dispersion relations, we demonstrate that SHARP can quantitatively reproduce important kinetic features of the linear regime. Our simulations have a superior ability to control energy non-conservation and avoid numerical heating in comparison to common second-order schemes. We provide a natural definition for convergence of a general PIC algorithm: the complement of physical modes captured by the simulation, i.e., those that lie above the Poisson noise, must grow commensurately with the resolution. This implies that it is necessary to simultaneously increase the number of particles per cell and decrease the cell size. We demonstrate that traditional ways for testing for convergence fail, leading to plateauing of the energy error. This new PIC code enables us to faithfully study the long-term evolution of plasma problems that require absolute control of the energy and momentum conservation.

Many-degree-scale gamma-ray halos are expected to surround extragalactic high-energy gamma ray sources. These arise from the inverse Compton emission of an intergalactic population of relativistic electron/positron pairs generated by the annihilation of >100 GeV gamma rays on the extragalactic background light. These are typically anisotropic due to the jetted structure from which they originate or the presence of intergalactic magnetic fields. Here we propose a novel method for detecting these inverse-Compton gamma-ray halos based upon this anisotropic structure. Specifically, we show that by stacking suitably defined angular power spectra instead of images it is possible to robustly detect gamma-ray halos with existing Fermi Large Area Telescope (LAT) observations for a broad class of intergalactic magnetic fields. Importantly, these are largely insensitive to systematic uncertainties within the LAT instrumental response or associated with contaminating astronomical sources.

Despite over forty years of active research, the nature of the white dwarf progenitors of Type Ia supernovae remains unclear. However, in the last decade, various progenitor scenarios have highlighted the need for detonations to be the primary mechanism by which these white dwarfs are consumed, but it is unclear how these detonations are triggered. In this paper we study how detonations are spontaneously initiated due to temperature inhomogeneities, e.g., hotspots, in burning nuclear fuel in a simplified physical scenario. Following the earlier work by Zel'Dovich, we describe the physics of detonation initiation in terms of the comparison between the spontaneous wave speed and the Chapman-Jouguet speed. We develop an analytic expression for the spontaneous wave speed and utilize it to determine a semi-analytic criterion for the minimum size of a hotspot with a linear temperature gradient between a peak and base temperature for which detonations in burning carbon-oxygen material can occur. Our results suggest that spontaneous detonations may easily form under a diverse range of conditions, likely allowing a number of progenitor scenarios to initiate detonations that burn up the star.

Motivated by the nonlinear star formation efficiency found in recent numerical simulations by a number of workers, we perform high-resolution adaptive mesh refinement simulations of star formation in self-gravitating turbulently driven gas. As we follow the collapse of this gas, we find that the character of the flow changes at two radii, the disk radius $r_d$, and the radius $r_*$ where the enclosed gas mass exceeds the stellar mass. Accretion starts at large scales and works inwards. In line with recent analytical work, we find that the density evolves to a fixed attractor, $\rho(r,t ) \rightarrow \rho(r)$, for $r_d<r<r_*$; mass flows through this structure onto a sporadically gravitationally unstable disk, and from thence onto the star. In the bulk of the simulation box we find that the random motions $v_T \sim r^p$ with $p \sim 0.5$, in agreement with Larson's size-linewidth relation. In the vicinity of massive star forming regions we find $ p \sim 0.2-0.3$, as seen in observations. For $r<r_*$, $v_T$ increases inward, with $p=-1/2$. Finally, we find that the total stellar mass $M_*(t)\sim t^2$ in line with previous numerical and analytic work that suggests a nonlinear rate of star formation.

The annihilation of TeV photons from extragalactic TeV sources and the extragalactic background light produces ultrarelativistic $e^{\pm}$ beams, which are subject to powerful plasma instabilities that sap their kinetic energy. Here we study the linear phase of the plasma instabilities that these pair beams drive. To this end, we calculate the linear growth rate of the beam plasma and oblique instability in the electrostatic approximation in both the reactive and kinetic regimes, assuming a Maxwell-J{\"u}ttner distribution for the pair beam. We reproduce the well-known reactive and kinetic growth rates for both the beam plasma and oblique mode. We demonstrate for the oblique instability that there is a broad spectrum of unstable modes that grow at the maximum rate for a wide range of beam temperatures and wave vector orientations relative to the beam. We also delineate the conditions for applicability for the reactive and kinetic regimes and find that the beam plasma mode transitions to the reactive regime at a lower Lorentz factor than the oblique mode due to a combination of their different scalings and the anisotropy of the velocity dispersions. Applying these results to the ultrarelativistic $e^{\pm}$ beams from TeV blazars, we confirm that these beams are unstable to both, the kinetic oblique mode and the reactive beam-plasma mode. These results are important in understanding how powerful plasma instabilities may sap the energy of the ultrarelativistic $e^{\pm}$ beams as they propagate through intergalactic space.

Extended inverse Compton halos are generally anticipated around extragalactic sources of gamma rays with energies above 100 GeV. These result from inverse Compton scattered cosmic microwave background photons by a population of high-energy electron/positron pairs produced by the annihilation of the high-energy gamma rays on the infrared background. Despite the observed attenuation of the high-energy gamma rays, the halo emission has yet to be directly detected. Here, we demonstrate that in most cases these halos are expected to be highly anisotropic, distributing the up-scattered gamma rays along axes defined either by the radio jets of the sources or oriented perpendicular to a global magnetic field. We present a pedagogical derivation of the angular structure in the inverse Compton halo and provide an analytic formalism that facilitates the generation of mock images. We discuss exploiting this fact for the purpose of detecting gamma-ray halos in a set of companion papers.

A dynamical characterization of dark matter dominated dwarf galaxies from their observed effects on galactic disks (i.e. Galactoseismology) has remained an elusive goal. Here, we present preliminary results from spectroscopic observations of three clustered Cepheid candidates identified from $K$-band light curves towards Norma. The average heliocentric radial velocity of these stars is $\sim$ 156 km/s, which is large and distinct from that of the Galaxy's stellar disk. These objects at $l \sim 333 ^\circ$ and $b \sim -1 ^\circ$ are therefore halo stars; using the $3.6~\micron$ period-luminosity relation of Type I Cepheids, they are at $\sim$ 73 kpc. Our ongoing $I$-band photometry indicates variability on the same time scale as the period determined from the $K_{s}$-band light curve. Distances determined from the $K$-band period-luminosity relation and the 3.6 $\micron$ period-luminosity relation are comparable. The observed radial velocity of these stars agrees roughly with predictions from dynamical models. If these stars are indeed members of the predicted dwarf galaxy that perturbed the outer HI disk of the Milky Way, this would mark the first application of Galactoseismology.

TeV-blazars potentially heat the intergalactic medium (IGM) as their gamma rays interact with photons of the extragalactic background light to produce electron-positron pairs, which lose their kinetic energy to the surrounding medium through plasma instabilities. This results in a heating mechanism that is only weakly sensitive to the local density, and therefore approximately spatially uniform, naturally producing an inverted temperature-density relation in underdense regions. In this paper we go beyond the approximation of uniform heating and quantify the heating rate fluctuations due to the clustered distribution of blazars and how this impacts on the thermal history of the IGM. We analytically compute a filtering function that relates the heating rate fluctuations to the underlying dark matter density field. We implement it in the cosmological code GADGET-3 and perform large scale simulations to determine the impact of inhomogeneous heating. We show that, because of blazar clustering, blazar heating is inhomogeneous for z>= 2. At high redshift, the temperature-density relation shows an important scatter and presents a low temperature envelope of unheated regions, in particular at low densities and within voids. However, the median temperature of the IGM is close to that in the uniform case, albeit slightly lower at low redshift. We find that blazar heating is more complex than initially assumed and that the temperature-density relation is not unique. Our analytic model for the heating rate fluctuations couples well with large scale simulations and provides a cost-effective alternative to subgrid models.

White dwarf (WD) binary mergers are possible progenitors to a number of unusual stars and transient phenomena, including type Ia supernovae. To date, simulations of mergers have not included magnetic fields, even though they are believed to play a significant role in the evolution of the merger remnant. We simulated a 0.625 - 0.65 $M_{\odot}$ carbon-oxygen WD binary merger in the magnetohydrodynamic moving mesh code Arepo. Each WD was given an initial dipole field with a surface value of $\sim10^3$ G. As in simulations of merging double neutron star binaries, we find exponential field growth within Kelvin-Helmholtz instability-generated vortices during the coalescence of the two stars. The final field has complex geometry, and a strength $>10^{10}$ G at the center of the merger remnant. Its energy is $\sim2\times10^{47}$ ergs, $\sim0.2$% of the remnant's total energy. The strong field likely influences further evolution of the merger remnant by providing a mechanism for angular momentum transfer and additional heating, potentially helping to ignite carbon fusion.

We present numerical evidence of dynamic star formation in which the accreted stellar mass grows superlinearly with time, roughly as $t^2$. We perform simulations of star formation in self-gravitating hydrodynamic and magneto-hydrodynamic turbulence that is continuously driven. By turning the self-gravity of the gas in the simulations on or off, we demonstrate that self-gravity is the dominant physical effect setting the mass accretion rate at early times before feedback effects take over, contrary to theories of turbulence-regulated star formation. We find that gravitational collapse steepens the density profile around stars, generating the power-law tail on what is otherwise a lognormal density probability distribution function. Furthermore, we find turbulent velocity profiles to flatten inside collapsing regions, altering the size-linewidth relation. This local flattening reflects enhancements of turbulent velocity on small scales, as verified by changes to the velocity power spectra. Our results indicate that gas self-gravity dynamically alters both density and velocity structures in clouds, giving rise to a time-varying star formation rate. We find that a substantial fraction of the gas that forms stars arrives via low density flows, as opposed to accreting through high density filaments.

NGC1266 is a nearby lenticular galaxy that harbors a massive outflow of molecular gas powered by the mechanical energy of an active galactic nucleus (AGN). It has been speculated that such outflows hinder star formation (SF) in their host galaxies, providing a form of feedback to the process of galaxy formation. Previous studies, however, indicated that only jets from extremely rare, high power quasars or radio galaxies could impart significant feedback on their hosts. Here we present detailed observations of the gas and dust continuum of NGC1266 at millimeter wavelengths. Our observations show that molecular gas is being driven out of the nuclear region at $\dot{M}_{\rm out} \approx 110 M_\odot$ yr$^{-1}$, of which the vast majority cannot escape the nucleus. Only 2 $M_\odot$ yr$^{-1}$ is actually capable of escaping the galaxy. Most of the molecular gas that remains is very inefficient at forming stars. The far-infrared emission is dominated by an ultra-compact ($\lesssim50$pc) source that could either be powered by an AGN or by an ultra-compact starburst. The ratio of the SF surface density ($\Sigma_{\rm SFR}$) to the gas surface density ($\Sigma_{\rm H_2}$) indicates that SF is suppressed by a factor of $\approx 50$ compared to normal star-forming galaxies if all gas is forming stars, and $\approx$150 for the outskirt (98%) dense molecular gas if the central region is is powered by an ultra-compact starburst. The AGN-driven bulk outflow could account for this extreme suppression by hindering the fragmentation and gravitational collapse necessary to form stars through a process of turbulent injection. This result suggests that even relatively common, low-power AGNs are able to alter the evolution of their host galaxies as their black holes grow onto the M-$\sigma$ relation.

Very-high energy gamma-rays from extragalactic sources pair-produce off of the extragalactic background light, yielding an electron-positron pair beam. This pair beam is unstable to various plasma instabilities, especially the "oblique" instability, which can be the dominant cooling mechanism for the beam. However, recently, it has been claimed that nonlinear Landau damping renders it physically irrelevant by reducing the effective damping rate to a low level. Here, we show with numerical calculations that the effective damping rate is $8\times 10^{-4}$ of the growth rate of the linear instability, which is sufficient for the "oblique" instability to be the dominant cooling mechanism of these pair beams. In particular, we show that previous estimates of this rate ignored the exponential cutoff in the scattering amplitude at large wavenumber and assumed that the damping of scattered waves entirely depends on collisions, ignoring collisionless processes. We find that the total wave energy eventually grows to approximate equipartition with the beam by increasingly depositing energy into long wavelength modes. As we have not included the effect of nonlinear wave-wave interactions on these long wavelength modes, this scenario represents the "worst-case" scenario for the oblique instability. As it continues to drain energy from the beam at a faster rate than other processes, we conclude that the "oblique" instability is sufficiently strong to make it the physically dominant cooling mechanism for high-energy pair beams in the intergalactic medium.

We present a model of star formation in self-gravitating turbulent gas. We treat the turbulent velocity $v_T$ as a dynamical variable, and assume that it is adiabatically heated by the collapse. The theory predicts the run of density, infall velocity, and turbulent velocity, and the rate of star formation in compact massive gas clouds. The turbulent pressure is dynamically important at all radii, a result of the adiabatic heating. The system evolves toward a coherent spatial structure with a fixed run of density, $\rho(r,t)\to\rho(r)$; mass flows through this structure onto the central star or star cluster. We define the sphere of influence of the accreted matter by $m_*=M_g(r_*)$, where $m_*$ is the stellar plus disk mass in the nascent star cluster and $M_g(r)$ is the gas mass inside radius $r$. The density is given by a broken power law with a slope $-1.5$ inside $r_*$ and $\sim -1.6$ to $-1.8$ outside $r_*$. Both $v_T$ and the infall velocity $|u_r|$ decrease with decreasing $r$ for $r>r_*$; $v_T(r)\sim r^p$, the size-linewidth relation, with $p\approx0.2-0.3$, explaining the observation that Larson's Law is altered in massive star forming regions. The infall velocity is generally smaller than the turbulent velocity at $r>r_*$. For $r<r_*$, the infall and turbulent velocities are again similar, and both increase with decreasing $r$ as $r^{-1/2}$, with a magnitude about half of the free-fall velocity. The accreted (stellar) mass grows super-linearly with time, $\dot M_*=\phi M_{\rm cl}(t/\tau_{ff})^2$, with $\phi$ a dimensionless number somewhat less than unity, $M_{\rm cl}$ the clump mass and $\tau_{ff}$ the free-fall time of the clump. We suggest that small values of p can be used as a tracer of convergent collapsing flows.

We infer the past orbit of the Sagittarius (Sgr) dwarf galaxy in the Milky Way halo by integrating backwards from its observed position and proper motions, including the effects of dynamical friction. Given measured proper motions, we show that there is a relation between the eccentricity ($e$) of Sgr's orbit and the mass of the Milky Way ($M_{T}$) in the limit of no dynamical friction. That relation can be fit by a power-law of the form: $e \approx 0.49 \left(M_{T}/10^{12} M_{\odot}\right)^{-0.88}$. At a fixed Milky Way mass, the dynamical friction term increases the mean eccentricity of the orbit and lowers the spread in eccentricities in proportion to the mass of the Sgr dwarf. We explore the implications of various observational constraints on Sgr's apocenter on the $e-M$ relation; Sgr masses outside the range $10^{9} M_{\odot} \la M_{\rm Sgr} \la 5 \times 10^{10} M_{\odot}$ are precluded, for Milky Way masses $\sim 1 - 2.5 \times 10^{12} M_{\odot}$. If Belokurov et al.'s (2014) observations represent the farthest point of Sgr's stream, then Milky Way masses in excess of $2 \times 10^{12} M_{\odot}$ are excluded for $M_{\rm Sgr} \la 10^{10} M_{\odot}$. Deeper observations of Sgr's tidal debris, from upcoming surveys such as GAIA, will allow better measurement of the Milky Way mass and of the Sgr dwarf.

We present new Spectrographic Areal Unit for Research on Optical Nebulae (SAURON) integral-field spectroscopy and Swift Ultraviolet Optical Telescope (UVOT) observations of molecular outflow host galaxy NGC 1266 that indicate NGC 1266 has experienced a rapid cessation of star formation. Both the SAURON maps of stellar population age and the Swift UVOT observations demonstrate the presence of young ($< 1$ Gyr) stellar populations within the central 1 kpc, while existing Combined Array for Research in Millimeter-wave Astronomy (CARMA) CO(1--0) maps indicate that the sites of current star formation are constrained to the inner few hundred parsecs of the galaxy only. The optical spectrum of NGC 1266 from Moustakas & Kennicutt (2006) reveal a characteristic post-starburst (K+A) stellar population and Davis et al. (2012) confirm that ionized gas emission in the system originate from a shock. Galaxies with K+A spectra and shock-like ionized gas line ratios may comprise an important, overlooked segment of the post-starburst population, containing exactly those objects in which the AGN is actively expelling the star-forming material. While AGN activity is not the likely driver of the post-starburst event that occurred 500 Myr ago, the faint spiral structure seen in the Hubble Space Telescope (HST) Wide-field Camera 3 (WFC3) Y-, J- and H-band imaging seems to point to the possibility of gravitational torques being the culprit. If the molecular gas were driven into the center at the same time as the larger scale galaxy disk underwent quenching, the AGN might be able to sustain the presence of molecular gas for $\gtrsim 1$ Gyr by cyclically injecting turbulent kinetic energy into the dense molecular gas via a radio jet, inhibiting star formation.

The extragalactic gamma-ray sky at TeV energies is dominated by blazars, a subclass of accreting super-massive black holes with powerful relativistic outflows directed at us. Only constituting a small fraction of the total power output of black holes, blazars were thought to have a minor impact on the universe at best. As we argue here, the opposite is true and the gamma-ray emission from TeV blazars can be thermalized via beam-plasma instabilities on cosmological scales with order unity efficiency, resulting in a potentially dramatic heating of the low-density intergalactic medium. Here, we review this novel heating mechanism and explore the consequences for the formation of structure in the universe. In particular, we show how it produces an inverted temperature-density relation of the intergalactic medium that is in agreement with observations of the Lyman-alpha forest. This suggests that blazar heating can potentially explain the paucity of dwarf galaxies in galactic halos and voids, and the bimodality of galaxy clusters. This also transforms our understanding of the evolution of blazars, their contribution to the extra-galactic gamma-ray background, and how their individual spectra can be used in constraining intergalactic magnetic fields.

Fermi has been instrumental in constraining the luminosity function and redshift evolution of gamma-ray bright blazars. This includes limits upon the spectrum and anisotropy of the extragalactic gamma-ray background (EGRB), redshift distribution of nearby Fermi active galactic nuclei (AGN), and the construction of a log(N)-log(S) relation. Based upon these, it has been argued that the evolution of the gamma-ray bright blazar population must be much less dramatic than that of other AGN. However, critical to such claims is the assumption that inverse Compton cascades reprocess emission above a TeV into the Fermi energy range, substantially enhancing the strength of the observed limits. Here we demonstrate that in the absence of such a process, due, e.g., to the presence of virulent plasma beam instabilities that preempt the cascade, a population of TeV-bright blazars that evolve similarly to quasars is consistent with the population of hard gamma-ray blazars observed by Fermi. Specifically, we show that a simple model for the properties and luminosity function is simultaneously able to reproduce their log(N)-log(S) relation, local redshift distribution, and contribution to the EGRB and its anisotropy without any free parameters. Insofar the naturalness of a picture in which the hard gamma-ray blazar population exhibits the strong redshift evolution observed in other tracers of the cosmological history of accretion onto halos is desirable, this lends support for the absence of the inverse Compton cascades and the existence of the beam plasma instabilities.

The merger of two carbon-oxygen white dwarfs can lead either to a spectacular transient, stable nuclear burning or a massive, rapidly rotating white dwarf. Simulations of mergers have shown that the outcome strongly depends on whether the white dwarfs are similar or dissimilar in mass. In the similar-mass case, both white dwarfs merge fully and the remnant is hot throughout, while in the dissimilar case, the more massive, denser white dwarf remains cold and essentially intact, with the disrupted lower mass one wrapped around it in a hot envelope and disk. In order to determine what constitutes "similar in mass" and more generally how the properties of the merger remnant depend on the input masses, we simulated unsynchronized carbon-oxygen white dwarf mergers for a large range of masses using smoothed-particle hydrodynamics. We find that the structure of the merger remnant varies smoothly as a function of the ratio of the central densities of the two white dwarfs. A density ratio of 0.6 approximately separates similar and dissimilar mass mergers. Confirming previous work, we find that the temperatures of most merger remnants are not high enough to immediately ignite carbon fusion. During subsequent viscous evolution, however, the interior will likely be compressed and heated as the disk accretes and the remnant spins down. We find from simple estimates that this evolution can lead to ignition for many remnants. For similar-mass mergers, this would likely occur under sufficiently degenerate conditions that a thermonuclear runaway would ensue.

Inverse-Compton cascades initiated by energetic gamma rays (E>100 GeV) enhance the GeV emission from bright, extragalactic TeV sources. The absence of this emission from bright TeV blazars has been used to constrain the intergalactic magnetic field (IGMF), and the stringent limits placed upon the unresolved extragalactic gamma-ray background (EGRB) by Fermi has been used to argue against a large number of such objects at high redshifts. However, these are predicated upon the assumption that inverse-Compton scattering is the primary energy-loss mechanism for the ultra-relativistic pairs produced by the annihilation of the energetic gamma rays on extragalactic background light photons. Here we show that for sufficiently bright TeV sources (isotropic-equivalent luminosities >10^{42} erg/s) plasma beam instabilities, specifically the "oblique" instability, present a plausible mechanism by which the energy of these pairs can be dissipated locally, heating the intergalactic medium. Since these instabilities typically grow on timescales short in comparison to the inverse-Compton cooling rate, they necessarily suppress the inverse-Compton cascades. As a consequence, this places a severe constraint upon efforts to limit the IGMF from the lack of a discernible GeV bump in TeV sources. Similarly, it considerably weakens the Fermi limits upon the evolution of blazar populations. Specifically, we construct a TeV-blazar luminosity function from those objects presently observed and find that it is very well described by the quasar luminosity function at z~0.1, shifted to lower luminosities and number densities, suggesting that both classes of sources are regulated by similar processes. Extending this relationship to higher redshifts, we show that the magnitude and shape of the EGRB above ~10 GeV is naturally reproduced with this particular example of a rapidly evolving TeV-blazar luminosity function.

The Universe is opaque to extragalactic very high-energy gamma rays (VHEGRs, E>100 GeV) because they annihilate and pair produce on the extragalactic background light. The resulting ultra-relativistic pairs are assumed to lose energy through inverse Compton scattering of CMB photons. In Broderick et al. (2011, Paper I of this three paper series), we argued that instead powerful plasma instabilities in the ultra-relativistic pair beam dissipate the kinetic energy of the TeV-generated pairs locally, heating the intergalactic medium (IGM). Here, we explore the effect of this heating upon the thermal history of the IGM. We collate the observed extragalactic VHEGR sources to determine a local VHEGR heating rate and correct for the pointed nature of VHEGR observations using Fermi observations of high and intermediate peaked BL Lacs. Because the local extragalactic VHEGR flux is dominated by TeV blazars, we tie the TeV blazar luminosity density to the quasar luminosity density, and produce a VHEGR heating rate as a function of redshift. This heating is relatively homogeneous for z<~4 with increasing spatial variation at higher redshift (order unity at z~6). This new heating process dominates photoheating at low redshift and the inclusion of TeV blazar heating qualitatively and quantitatively changes the structure and history of the IGM. TeV blazars produce a uniform volumetric heating rate that is sufficient to increase the temperature of the mean density IGM by nearly an order of magnitude, and at low densities by substantially more, naturally producing an inverted equation of state inferred by observations of the Ly-alpha forest, a feature that is difficult to reconcile with standard reionization models. Finally, we close with a discussion on the possibility of detecting this hot low-density IGM, but find that such measurements are currently not feasible. (abridged)

A subset of blazars emit TeV gamma rays which annihilate and pair produce on the extragalactic background light. We have argued in Broderick et al. (2011, Paper I) that plasma beam instabilities can dissipate the pairs' energy locally. This heats the intergalactic medium and dramatically increases its entropy after redshift z~2, with important implications for structure formation: (1) This suggests a scenario for the origin of the cool core (CC)/non-cool core (NCC) bimodality in galaxy clusters and groups. Early forming galaxy groups are unaffected because they can efficiently radiate the additional entropy, developing a CC. However, late forming groups do not have sufficient time to cool before the entropy is gravitationally reprocessed through successive mergers - counteracting cooling and raising the core entropy further. Hence blazar heating works different than feedback by active galactic nuclei, which balances radiative cooling but is unable to transform CC into NCC clusters due to the weak coupling to the cluster gas. (2) We predict a suppression of the Sunyaev-Zel'dovich power spectrum on angular scales smaller than 5' due to the globally reduced central pressure of groups and clusters forming after z~1. (3) Our redshift dependent entropy floor increases the characteristic halo mass below which dwarf galaxies cannot form by a factor of ~10 (50) at mean density (in voids) over that found in models that include photoionization alone. This prevents the formation of late forming dwarf galaxies (z<2) with masses ranging from 10^{10} to 10^{11} M_sun for redshifts z~2 to 0, respectively. This may help resolve the "missing satellite problem" in the Milky Way and the "void phenomenon" of the low observed abundances of dwarf satellites compared to cold dark matter simulations and may bring the observed early star formation histories into agreement with galaxy formation models. (abridged)

It has been realised only recently that TeV emission from blazars can significantly heat the intergalactic medium (IGM) by pair-producing high-energy electrons and positrons, which in turn excite vigorous plasma instabilities, leading to a local dissipation of the pairs' kinetic energy. In this work, we use cosmological hydrodynamical simulations to model the impact of this blazar heating on the Lyman-alpha forest at redshifts z~2-3. We find that blazar heating produces an inverted temperature-density relation in the IGM and naturally resolves many of the problems present in previous simulations of the forest that included photoheating alone. In particular, our simulations with blazar heating simultaneously reproduce the observed effective optical depth and temperature as a function of redshift, the observed probability distribution functions of the transmitted flux, and the observed flux power spectra, over the full redshift range 2<z<3 analysed here. Additionally, by deblending the Lyman-alpha forest into a sum of thermally broadened individual lines, we find superb agreement with the observed lower cutoff of the line-width distribution and abundances of neutral hydrogen column densities. Using the most recent constraints on the cosmic ultraviolet (UV) background, this excellent agreement with observations does not require rescaling the amplitude of the UV background; a procedure that was routinely used in the past to match the observed level of transmitted flux. We also show that our blazar-heated model matches the data better than standard simulations even when such a rescaling is allowed. This concordance between Lyman-alpha data and simulations, which are based on the most recent cosmological parameters, suggests that the inclusion of blazar heating alleviates previous tensions on constraints for sigma_8 derived from Lyman-alpha measurements and other cosmological data. [abridged]

Recent studies have shown that for suitable initial conditions both super- and sub-Chandrasekhar mass carbon-oxygen white dwarf mergers produce explosions similar to observed SNe Ia. The question remains, however, how much fine tuning is necessary to produce these conditions. We performed a large set of SPH merger simulations, sweeping the possible parameter space. We find trends for merger remnant properties, and discuss how our results affect the viability of our recently proposed sub-Chandrasekhar merger channel for SNe Ia.

We describe ongoing work on a new method that allows one to determine the mass and relative position (in galactocentric radius and azimuth) of galactic companions purely from analysis of observed disturbances in gas disks. Here, we demonstrate the validity of this method, which we call Tidal Analysis, by applying it to local spirals with known optical companions, namely M51 and NGC 1512. These galaxies span the range from having a very low mass companion ($\sim$ one-hundredth the mass of the primary galaxy) to a fairly massive companion ($\sim$ one-third the mass of the primary galaxy). This approach has broad implications for many areas of astrophysics -- for the indirect detection of dark matter (or dark-matter dominated dwarf galaxies), and for galaxy evolution in its use as a decipher for the dynamical impact of satellites on galactic disks. Here, we provide a proof of principle of the method by applying it to infer and quantitatively characterize optically visible galactic companions of local spirals, from the analysis of observed disturbances in outer gas disks.

Essentially all stars form in giant molecular clouds (GMCs). However, inside GMCs, most of the gas does not participate in star formation; rather, denser gas accumulates in clumps in the GMC, with the bulk of the stars in a given GMC forming in a few of the most massive clumps. In the Milky Way, these clumps have masses $M_{\rm cl}\lesssim 5\times 10^{-2}$ of the GMC, radii $r_{\rm cl} \sim 1$pc, and free-fall times $\tau_{\rm cl} \sim 2\times 10^5\yr$. We show that clumps inside giant molecular clouds should accrete at a modified Bondi accretion rate, which depends on clump mass as $\dot M_{\rm cl}\sim M_{\rm cl}^{5/4}$. This rate is initially rather slow, usually slower than the initial star formation rate inside the clump (we adopt the common assumption that inside the clump, $\dot M_*=\epsilon_{\rm ff} M_{\rm cl}/\tau_{\rm cl}$, with $\epsilon_{\rm ff} \approx 0.017$). However, after $\sim 2$ GMC free-fall times $\tau_{\rm GMC}$, the clump accretion rate accelerates rapidly; formally, the clump can accrete the entire GMC in $\sim 3\tau_{\rm GMC}$. At the same time, the star formation rate accelerates, tracking the Bondi accretion rate. If the GMC is disrupted by feedback from the largest clump, half the stars in that clump form in the final $\taug$ before the GMC is disrupted. The theory predicts that the distribution of effective star formation rates, measured per GMC free-fall time, is broad, ranging from $\sim 0.001$ up to 0.1 or larger and that the mass spectrum of star clusters is flatter than that of clumps, consistent with observations.

We discuss the nature of rings that may exist around extrasolar planets. Taking the general properties of rings around the gas giants in the Solar System, we infer the likely properties of rings around exoplanets that reside inside the ice line. Due to their proximity to their host star, rings around such exoplanets must primarily consist of rocky materials. However, we find that despite the higher densities of rock compared to ice, most of the observed extrasolar planets with reliable radii measurements have sufficiently large Roche radii to support rings. For the currently known transiting extrasolar planets, Poynting-Robertson drag is not effective in significantly altering the dynamics of individual ring particles over a time span of $10^8$ years provided that they exceed about 1 m in size. In addition, we show that significantly smaller ring particles can exist in optically thick rings, for which we find typical ring lifetimes ranging from a few times $10^6$ to a few times $10^9$ years. Most interestingly, we find that many of the rings could have nontrivial Laplacian planes due to the increased effects of the orbital quadrupole caused by the exoplanets' proximity to their host star, allowing a constraint on the $J_2$ of extrasolar planets from ring observations. This is particular exciting, since a planet's $J_2$ reveals information about its interior structure. Furthermore, measurements of an exoplanet's oblateness and of its $J_2$, from warped rings, would together place limits on its spin period. Based on the constraints that we have derived for extrasolar rings, we anticipate that the best candidates for ring detections will come from transit observations by the Kepler spacecraft of extrasolar planets with semi-major axes $\sim 0.1$ AU and larger.

We develop a perturbative approach to study the excitation of disturbances in the extended atomic hydrogen (HI) discs of galaxies produced by passing dark matter subhaloes. The shallow gravitational potential of the dark matter subhaloes (compared to the primary halo) allows us to use the epicyclic approximation, the equations of which we solve by modal analysis, i.e., assuming a disc is composed of N radial rings with M modes. We show that properties of dark matter subhaloes can be inferred from the profile and amplitude of the modal energy of the disc. Namely, we find that the overall amplitude of the response gives the mass of the dark sub-halo. Motivated by this modal analysis, we then show that the density response shows similar features. Finally, we show that our results agree with those from full hydrodynamic simulation. We find a simple scaling relation between the satellite mass and Fourier amplitudes of the resultant surface density of the gas disc where the effective Fourier amplitude (essentially a sum over the low order modes) scales as $m_{s}^{1/2}$, where $m_{s}$ is the satellite mass. The utility of this relation is that it can be readily applied to an observed HI map to deduce the satellite mass without recourse to full numerical simulation. This will greatly aid us in analyzing large samples of spiral galaxies to constrain the population of dwarf satellites in the Local Volume.

Theories that attempt to explain cosmic acceleration by modifying gravity typically introduces a long-range scalar force that needs to be screened on small scales. One common screening mechanism is the chameleon, where the scalar force is screened in environments with a sufficiently deep gravitational potential, but acts unimpeded in regions with a shallow gravitational potential. This leads to a variation in the overall gravitational G with environment. We show such a variation can occur within a star itself, significantly affecting its evolution and structure, provided that the host galaxy is unscreened. The effect is most pronounced for red giants, which would be smaller by a factor of tens of percent and thus hotter by 100's of K, depending on the parameters of the underlying scalar-tensor theory. Careful measurements of these stars in suitable environments (nearby dwarf galaxies not associated with groups or clusters) would provide constraints on the chameleon mechanism that are four orders of magnitude better than current large scale structure limits, and two orders of magnitude better than present solar system tests.

Type Ia supernovae are generally thought to be due to the thermonuclear explosions of carbon-oxygen white dwarfs with masses near the Chandrasekhar mass. This scenario, however, has two long-standing problems. First, the explosions do not naturally produce the correct mix of elements, but have to be finely tuned to proceed from sub-sonic deflagration to super-sonic detonation. Second, population models and observations give formation rates of near-Chandrasekhar white dwarfs that are far too small. Here, we suggest that type Ia supernovae instead result from mergers of roughly equal-mass carbon-oxygen white dwarfs, including those that produce sub-Chandrasekhar mass remnants. Numerical studies of such mergers have shown that the remnants consist of rapidly rotating cores that contain most of the mass and are hottest in the center, surrounded by dense, small disks. We argue that the disks accrete quickly, and that the resulting compressional heating likely leads to central carbon ignition. This ignition occurs at densities for which pure detonations lead to events similar to type Ia supernovae. With this merger scenario, we can understand the type Ia rates, and have plausible reasons for the observed range in luminosity and for the bias of more luminous supernovae towards younger populations. We speculate that explosions of white dwarfs slowly brought to the Chandrasekhar limit---which should also occur---are responsible for some of the "atypical" type Ia supernovae.

The formation of planetesimals via gravitational instability of the dust layer in a protoplanetary disks demands that there be local patches where dust is concentrated by a factor of $\sim$ a few $\times 10^3$ over the background value. Vortices in protoplanetary disks may concentrate dust to these values allowing them to be the nurseries of planetesimals. The concentration of dust in the cores of vortices increases the dust-gas ratio of the core compared to the background disk, creating a "heavy vortex." In this work, we show that these vortices are subject to an instability which we have called the heavy-core instability. Using Floquet theory, we show that this instability occurs in elliptical protoplanetary vortices when the gas-dust density of the core of the vortex is heavier than the ambient gas-dust density by a few tens of percent. The heavy-core instability grows very rapidly, with a growth timescale of a few vortex rotation periods. While the nonlinear evolution of this instability remains unknown, it will likely increase the velocity dispersion of the dust layer in the vortex because instability sets in well before sufficient dust can gather to form a protoplanetary seed. This instability may thus preclude vortices from being sites of planetesimal formation.

Diffusive nuclear burning of H by an underlying material capable of capturing protons can readily consume H from the surface of neutron stars (NSs) during their early cooling history. In the absence of subsequent accretion, it will be depleted from the photosphere. We now extend diffusive nuclear burning to He, motivated by the recent observation by Ho \& Heinke of a carbon atmosphere on the NS in the Cassiopeia A supernova remnant. We calculate the equilibrium structure of He on an underlying $\alpha$ capturing material, accounting for thermal, mass defect, and Coulomb corrections on the stratification of material with the same zeroth order $\mu_e = A/Z$. We show that Coulomb corrections dominate over thermal and mass defect corrections in the highly degenerate part of the envelope. We also show that the bulk of the He sits deep in the envelope rather than near the surface. Thus, even if the photospheric He abundance is low, the total He column could be substantially larger than the photospheric column, which may have implications for rapid surface evolution ($\approx 1$ yr timescales) of neutron stars. When nuclear reactions are taken into account, we find that for base temperatures $\gtrsim 1.6 \times 10^8$ K, He is readily captured onto C. As these high temperatures are present during the early stages of NS evolution, we expect that the primordial He is completely depleted from the NS surface like the case for primordial H. We also find that magnetic fields $\lesssim 10^{12}$ G do not affect our conclusions. Armed with the results of this work and our prior efforts, we expect that primordial H and He are depleted, and so any observed H or He on the surfaces of these NS must be due to subsequent accretion (with or without spallation). If this subsequent accretion can be prevented, the underlying mid-Z material would be exposed.

In low-collisionality plasmas, anisotropic heat conduction due to a magnetic field leads to buoyancy instabilities for any nonzero temperature gradient. We study analogous instabilities in degenerate {\it collisional} plasmas, i.e., when the electron collision frequency is large compared to the electron cyclotron frequency. Although heat conduction is nearly isotropic in this limit, the small residual anisotropy ensures that collisional degenerate plasmas are also convectively unstable independent of the sign of the temperature gradient. We show that the range of wavelengths that are unstable is independent of the magnetic field strength, while the growth time increases with decreasing magnetic field strength. We discuss the application of these collisional buoyancy instabilities to white dwarfs and neutron stars. Magnetic tension and the low specific heat of a degenerate plasma significantly limit their effectiveness; the most promising venues for growth are in the liquid oceans of young, weakly magnetized neutron stars ($B \lesssim 10^9$ G) and in the cores of young, high magnetic field white dwarfs ($B \sim 10^9$ G).

The mode of explosive burning in Type Ia SNe remains an outstanding problem. It is generally thought to begin as a subsonic deflagration, but this may transition into a supersonic detonation (the DDT). We argue that this transition leads to a breakout shock, which would provide the first unambiguous evidence that DDTs occur. Its main features are a hard X-ray flash (~20 keV) lasting ~0.01 s with a total radiated energy of ~10^{40} ergs, followed by a cooling tail. This creates a distinct feature in the visual light curve, which is separate from the nickel decay. This cooling tail has a maximum absolute visual magnitude of M_V = -9 to -10 at approximately 1 day, which depends most sensitively on the white dwarf radius at the time of the DDT. As the thermal diffusion wave moves in, the composition of these surface layers may be imprinted as spectral features, which would help to discern between SN Ia progenitor models. Since this feature should accompany every SNe Ia, future deep surveys (e.g., m=24) will see it out to a distance of approximately 80 Mpc, giving a maximum rate of ~60/yr. Archival data sets can also be used to study the early rise dictated by the shock heating (at about 20 days before maximum B-band light). A similar and slightly brighter event may also accompany core bounce during the accretion induced collapse to a neutron star, but with a lower occurrence rate.

The latest cosmological N-body simulations find two intriguing properties for dark matter haloes: (1) their radial density profile, rho, is better fit by a form that flattens to a constant at the halo center (the Einasto profile) than the widely-used NFW form; (2) the radial profile of the pseudo-phase-space density, rho/sigma3, on the other hand, continues to be well fit by a power law, as seen in earlier lower-resolution simulations. In this paper we use the Jeans equation to argue that (1) and (2) cannot both be true at all radii. We examine the implied radial dependence of rho/sigma3 over 12 orders of magnitude in radius by solving the Jeans equation for a broad range of input rho and velocity anisotropy beta. Independent of beta, we find that rho/sigma3 is approximately a power law only over the limited range of halo radius resolvable by current simulations (down to ~0.1% of the virial radius), and rho/sigma3 deviates significantly from a power-law below this scale for both the Einasto and NFW rho. The same conclusion also applies to a more general density-velocity relation rho/sigma_D^epsilon. Conversely, when we enforce rho/sigma^3 r^{-eta} as an input, none of the physically allowed rho (occurring for the narrow range 1.8<eta<1.9444) follows the Einasto form. We expect the next generation simulations with better spatial resolution to settle the debate: either the Einasto profile will continue to hold and rho/sigma3 will deviate from a power law, or rho/sigma3 will continue as a power law and rho will deviate from its current parameterizations.

Using a one-dimensional height integrated model, we calculate the evolution of an unequal mass binary black hole with a coplanar gas disk that contains a gap due to the presence of the secondary black hole. Viscous evolution of the outer circumbinary disk initially hardens the binary, while the inner disk drains onto the primary (central) black hole. As long as the inner disk remains cool and thin at low $\dot{M}_{\rm ext}$ (rather than becoming hot and geometrically thick), the mass of the inner disk reaches an asymptotic mass typically $\sim 10^{-3}-10^{-4}\Msun$. Once the semimajor axis shrinks below a critical value, angular momentum losses from gravitational waves dominate over viscous transport in hardening the binary. The inner disk then no longer responds viscously to the inspiraling black holes. Instead, tidal interactions with the secondary rapidly drive the inner disk into the primary. Tidal and viscous dissipation in the inner disk lead to a late time brightening in luminosity $L\propto t_{\rm minus}^{5/4}$, where $t_{\rm minus}$ is the time prior to the final merger. This late time brightening peaks $\sim 1$ day prior to the final merger at $\sim 0.1 L_{\rm Edd}$. This behavior is relatively robust because of self regulation in the coupled viscous-gravitational evolution of such binary systems. It constitutes a unique electromagnetic signature of a binary supermassive black hole merger and may allow the host galaxy to be identified if used in conjunction with the Laser Interferometric Space Antenna (LISA) localization.

Star clusters that form in nuclear rings appear to be at slightly larger radii than the gas. We argue that the star clusters move out from the gas in which they are formed because of satellite-disk tidal interactions. In calculating the dynamics of this star cluster and gas ring system, we include the effects of dynamical friction of the background stars in the host galaxy on the star cluster, and inflowing gas along the bar onto the nuclear ring at the two contact points. We show that the final separation is of the order of the Hill radius of the nuclear ring, which is typically 20-30% of its radius. Massive star clusters can reach half of this separation very quickly and produce a factor of a few enhancement in the gas surface density. If this leads to star formation in addition to the (ongoing) formation of star clusters near the contact points, a possible (initial) azimuthal age gradient may become diluted or even disappear. Finally, if the star cluster are massive and/or numerous enough, we expect the nuclear ring to migrate inward, away from the (possibly) associated (inner) Lindblad resonance. We discuss how these predictions may be tested observationally.

I examine the effectiveness of Kozai oscillations in the centres of galaxies and in particular the Galactic centre using standard techniques from celestial mechanics. In particular, I study the effects of a stellar bulge potential and general relativity on Kozai oscillations, which are induced by stellar discs. Lockmann et al (2008) recently suggested that Kozai oscillations induced by the two young massive stellar discs in the Galactic centre drives the orbits of the young stars to large eccentricity ($e\approx 1$). If some of these young eccentric stars are in binaries, they would be disrupted near pericentre, leaving one star in a tight orbit around the central SMBH and producing the S-star population. I find that the {\it spherical} stellar bulge suppresses Kozai oscillations, when its enclosed mass inside of a test body is of order the mass in the stellar disc(s). Since the stellar bulge in the Galactic centre is much larger than the stellar discs, Kozai oscillations {\it due to the stellar discs} are likely suppressed. Whether Kozai oscillations are induced from other nonspherical components to the potential (for instance, a flattened stellar bulge) is yet to be determined.

I calculate the action of a satellite, infalling through dynamical friction, on a coplanar gaseous disk of finite radial extent. The disk tides, raised by the infalling satellite, couple the satellite and disk. Dynamical friction acting on the satellite then shrinks the radius of the coupled satellite-disk system. Thus, the gas is ``shepherded'' to smaller radii. In addition, gas shepherding produces a large surface density enhancement at the disk edge. If the disk edge then becomes gravitationally unstable and fragments, it may give rise to enhanced star formation. On the other hand, if the satellite is sufficiently massive and dense, the gas may be transported from $\sim 100$ pc to inside of a 10 to 10s of parsecs before completely fragmenting into stars. I argue that gas shepherding may drive the fueling of active galaxies and central starbursts and I compare this scenario to competing scenarios. I argue that sufficiently large and dense super star clusters (acting as the shepherding satellites) can shepherd a gas disk down to ten to tens of parsecs. Inside of ten to tens of parsecs, another mechanism may operate, i.e., cloud-cloud collisions or a marginally (gravitationally) stable disk, that drives the gas $\lesssim 1$ pc, where it can be viscously accreted, feeding a central engine.

We study the long term evolution of magnetic fields generated by an initially unmagnetized collisionless relativistic $e^+e^-$ shock. Our 2D particle-in-cell numerical simulations show that downstream of such a Weibel-mediated shock, particle distributions are approximately isotropic, relativistic Maxwellians, and the magnetic turbulence is highly intermittent spatially, nonpropagating, and decaying. Using linear kinetic theory, we find a simple analytic form for these damping rates. Our theory predicts that overall magnetic energy decays like $(\omega_p t)^{-q}$ with $q \sim 1$, which compares favorably with simulations, but predicts overly rapid damping of short wavelength modes. Magnetic trapping of particles within the magnetic structures may be the origin of this discrepancy. We conclude that initially unmagnetized relativistic shocks in electron-positron plasmas are unable to form persistent downstream magnetic fields. These results put interesting constraints on synchrotron models for the prompt and afterglow emission from GRBs.

Following unstable ignition of carbon, but prior to explosion, a white dwarf (WD) in a Type Ia supernova (SN Ia) undergoes a simmering phase. During this time, a central convective region grows and encompasses ~1 Msun of the WD over a timescale of ~1000 yrs, which sets the thermal and turbulent profile for the subsequent explosion. We study this time-dependent convection and summarize some of the key features that differ from the traditional, steady-state case. We show that the long conductive timescale above the convective zone and the extraction of energy to heat the WD core leads to a decrease of the convective luminosity and characteristic velocities near the convective zone's top boundary. In addition, differences in the composition between the convective core and the conductive exterior will significantly alter the location of this boundary. In this respect, we find the biggest effect due to complete 22Ne sedimentation prior to carbon ignition. These effects add diversity to the possible WD models, which may alter the properties of the SN Ia explosion.

We study the long term evolution of magnetic fields generated by a collisionless relativistic $e^+e^-$ shock which is initially unmagnetized. Our 2D particle-in-cell numerical simulations show that downstream of such a Weibel-mediated shock, particle distributions are close to isotropic, relativistic Maxwellians, and the magnetic turbulence is highly intermittent spatially, with the non-propagating magnetic fields forming relatively isolated regions with transverse dimension $\sim 10-20$ skin depths. These structures decay in amplitude, with little sign of downstream merging. The fields start with magnetic energy density $\sim (0.1-0.2)$ of the upstream kinetic energy within the shock transition, but rapid downstream decay drives the fields to much smaller values, below $10^{-3}$ of equipartition after $10^3$ skin depths. In an attempt to construct a theory that follows field decay to these smaller values, we explore the hypothesis that the observed damping is a variant of Landau damping in an unmagnetized plasma. The model is based on the small value of the downstream magnetic energy density, which suggests that particle orbits are only weakly perturbed from straight line motion, if the turbulence is homogeneous. Using linear kinetic theory applied to electromagnetic fields in an isotropic, relativistic Maxwellian plasma, we find a simple analytic form for the damping rates, $\gamma_k$, in two and three dimensions for small amplitude, subluminous electromagnetic fields. We find that magnetic energy does damp due to phase mixing of current carrying particles as $(\omega_p t)^{-q}$ with $q \sim 1$. (abridged)

The triple nucleus of M31 consists of a population of old red stars in an eccentric disk (P1 and P2) and another population of younger A stars in a circular disk (P3) around M31's central supermassive black hole (SMBH). We argue that P1 and P2 determine the maximal radial extent of the younger A star population and provide the gas that fueled the starburst that generated P3. The eccentric stellar disk creates an $m=1$ non-axisymmetric perturbation to the potential. This perturbed potential drives gas into the inner parsec around the SMBH, if the pattern speed of the eccentric stellar disk is $\Omega_p \lesssim 3-10 {\rm km s^{-1} pc^{-1}}$. We show that stellar mass loss from P1 and P2 is sufficient to create a gravitationally unstable gaseous disk of $\sim 10^5\Msun$ every $0.1-1$ Gyrs, consistent with the 200 Myr age of P3. Similar processes may act in other systems to produce very compact nuclear starbursts.

We argue that the X-ray and UV flux illuminating the parsec-scale accretion disk around luminous active galactic nuclei (AGN) is super-Eddington with respect to the local far-infrared dust opacity. The far infrared opacity may be larger than in the interstellar medium of the Milky Way due to a combination of supersolar metallicity and the growth of dust grains in the dense accretion disk. Because of the irradiating flux, the outer accretion disk puffs up with a vertical thickness $h\sim R$. This provides a mechanism for generating a geometrically thick obscuring region from an intrinsically thin disk. We find obscuring columns $\sim 10^{22} - 10^{23} {\rm cm}^{-2}$, in reasonable agreement with observations.

RX J1856.5-3754 is one of the brightest nearby isolated neutron stars, and considerable observational resources have been devoted to it. However, current models are unable to satisfactorily explain the data. We show that our latest models of a thin, magnetic, partially ionized hydrogen atmosphere on top of a condensed surface can fit the entire spectrum, from X-rays to optical, of RX J1856.5-3754, within the uncertainties. In our simplest model, the best-fit parameters are an interstellar column density N_H \approx 1x10^20 cm^-2 and an emitting area with R^infty \approx 17 km (assuming a distance to RX J1856.5-3754 of 140 pc), temperature T^infty \approx 4.3x10^5 K, gravitational redshift z_g \sim 0.22, atmospheric hydrogen column y_H \approx 1 g cm^-2, and magnetic field B \approx (3-4)x10^12 G; the values for the temperature and magnetic field indicate an effective average over the surface. We also calculate a more realistic model, which accounts for magnetic field and temperature variations over the neutron star surface as well as general relativistic effects, to determine pulsations; we find there exist viewing geometries that produce pulsations near the currently observed limits. The origin of the thin atmospheres required to fit the data is an important question, and we briefly discuss mechanisms for producing these atmospheres. Our model thus represents the most self-consistent picture to date for explaining all the observations of RX J1856.5-3754.

The discovery of the first gravitationally redshifted spectral line from a neutron star (NS) by Cottam, Paerels and Mendez has triggered theoretical studies of the physics of atomic line formation in NS atmospheres. Chang, Bildsten and Wasserman showed that the hydrogenic Fe H$\alpha$ line formed above the photosphere of a bursting NS is intrinsically broad. We now include rotational broadening within general relativity and compare the resulting profile to that observed during Type I bursts from EXO 0748-676. We show that the fine structure splitting of the line precludes a meaningful constraint on the radius. Our fitting of the data show that the line forming Fe column is ${\rm log}_{10} (N_{\rm Fe, n=2}/{\rm cm^{-2}})=17.9_{-0.42}^{+0.27}$ and gravitational redshift $1+z =1.345_{-0.008}^{+0.005}$ with 95% confidence. We calculate the detectability of this spectral feature for a large range of spins and inclinations assuming that the emission comes from the entire surface. We find that at 300 (600) Hz only 10-20% (5-10%) of NSs would have spectral features as deep as that seen in EXO 0748-676.

Motivated by the measurement of redshifted Fe H$\alpha$ lines during type I X-ray bursts on EXO 0748-676 (Cottam, Paerels & Mendez), we study the formation of atomic Fe lines above the photosphere of a bursting neutron star ($k_BT_{\rm eff} \approx 1-2 {\rm keV}$). We discuss the effects of Stark broadening, resonant scattering and NLTE (level population) on the formation of hydrogenic Fe H$\alpha$, Ly$\alpha$ and P$\alpha$ lines. From the observed equivalent width of the Fe H$\alpha$ line, we find an implied Fe column of $1-3 \times 10^{20} {\rm cm}^{-2}$, which is 3-10 times larger than the Fe column calculated from the accretion/spallation model of Bildsten, Chang & Paerels. We also estimate that the implied Fe column is about a factor of 2-3 larger than a uniform solar metallicity atmosphere. We discuss the effects of rotational broadening and find that the rotation rate of \EXO must be slow, as confirmed by the recent measurement of a 45 Hz burst oscillation by Villarreal & Strohmayer. We also show that the Fe Ly$\alpha$ EW $\approx$ 15-20 eV (redshifted 11-15 eV) and the P$\alpha$ EW $\approx$ 4-7 eV (redshifted 3-5 eV) when the H$\alpha$ EW is 10 eV (redshifted 8 eV). The Ly$\alpha$ line is rotationally broadened to a depth of $\approx 10%$, making it difficult to observe with {\it Chandra}. We also show that radiative levitation can likely support the Fe column needed to explain the line.

We compute the rate of diffusive nuclear burning for hydrogen on the surface of a "magnetar" (Soft Gamma-Ray Repeater or Anomalous X-Ray Pulsar). We find that hydrogen at the photosphere will be burned on an extremely rapid timescale of hours to years, depending on composition of the underlying material. Improving on our previous studies, we explore the effect of a maximally thick "inert" helium layer, previously thought to slow down the burning rate. Since hydrogen diffuses faster in helium than through heavier elements, we find this helium buffer actually increases the burning rate for magnetars. We compute simple analytic scalings of the burning rate with temperature and magnetic field for a range of core temperature. We conclude that magnetar photospheres are very unlikely to contain hydrogen. This motivates theoretical work on heavy element atmospheres that are needed to measure effective temperature from the observed thermal emission and constrains models of AXPs that rely on magnetar cooling through thick light element envelopes.

The Kelvin-Helmholtz cooling epoch, lasting tens of seconds after the birth of a neutron star in a successful core-collapse supernova, is accompanied by a neutrino-driven wind. For magnetar-strength ($\sim10^{15}$ G) large scale surface magnetic fields, this outflow is magnetically-dominated during the entire cooling epoch.Because the strong magnetic field forces the wind to co-rotate with the protoneutron star,this outflow can significantly effect the neutron star's early angular momentum evolution, as in analogous models of stellar winds (e.g. Weber & Davis 1967). If the rotational energy is large in comparison with the supernova energy and the spindown timescale is short with respect to the time required for the supernova shockwave to traverse the stellar progenitor, the energy extracted may modify the supernova shock dynamics significantly. This effect is capable of producing hyper-energetic supernovae and, in some cases, provides conditions favorable for gamma ray bursts. We estimate spindown timescales for magnetized, rotating protoneutron stars and construct steady-state models of neutrino-magnetocentrifugally driven winds. We find that if magnetars are born rapidly rotating, with initial spin periods ($P$) of $\sim1$ millisecond, that of order $10^{51}-10^{52}$ erg of rotational energy can be extracted in $\sim10$ seconds. If magnetars are born slowly rotating ($P\gtrsim10$ ms) they can spin down to periods of $\sim1$ second on the Kelvin-Helmholtz timescale.

We extend our initial study of diffusive nuclear burning (DNB) for neutron stars (NSs) with Hydrogen atmospheres and an underlying layer of proton capturing nuclei. Our initial study showed that DNB can alter the photospheric abundance of Hydrogen on surprisingly short timescales ($10^{2-4}\yrs$). Significant composition evolution impacts the radiated thermal spectrum from the NS as well as its overall cooling rate. In this paper, we consider the case when the rate limiting step for the H consumption is diffusion to the burning layer, rather than the local nuclear timescale. This is relevant for NSs with surface temperatures in excess of $10^6 {\rm K}$, such as young ($<10^5$ yr) radio pulsars and accreting NSs in quiescence. When downward diffusion is the limiting rate in DNB, the rate of H consumption is suppressed by 1-2 orders of magnitude compared to a DNB estimate that assumes diffusive equilibrium. In order to apply our ongoing study to young neutron stars, we also include the important effects of strong magnetic fields ($B \sim 10^{12} {\rm G}$). In this initial study of magnetic modifications to DNB, we find that the H burning time is lengthened by 2-3 orders of magnitude for a $10^{12} {\rm G}$ field. However, even for NSs with dipole field strengths of $10^{12}$ G, we find that all of the H can be burned before the pulsar reaches an age of $\sim 10^5 \ {\rm yr}$, thus potentially revealing the underlying proton-capturing elements. Finally, we conclude by providing an overview of what can be learned about fallback and pulsar winds from measuring the surface composition of a young NS.

The gravitational redshift measured by Cottam, Paerels and Mendez for the neutron star (NS) in the low-mass X-ray binary EXO 0748-676 depends on the identification of an absorption line during a type I burst as the H$\alpha$ line from hydrogenic Fe. We show that Fe is present above the photosphere as long as $\dot M>4\times 10^{-13}M_\odot {\rm yr^{-1}}$ during the burst. In this limit, the total Fe column is $N_{\rm Fe}\approx 3\times 10^{19}{\rm cm^{-2}}$ for incident material of solar abundances and only depends on the nuclear physics of the proton spallation. The Fe destruction creates many heavy elements with $Z<26$ which may imprint photo-ionization edges on the NS spectra during a radius expansion event or in a burst cooling tail. Detecting these features in concert with those from Fe would confirm a redshift measurement. We also begin to address the radiative transfer problem, and find that a concentrated Fe layer with $kT=1.2-1.4 {\rm keV}$ and column $N_{\rm Fe}= 7-20 \times 10^{20} {\rm cm}^{-2}$ (depending on the line depth) above the hotter continuum photosphere is required to create the H$\alpha$ line of the observed strength. This estimate must be refined by considerations of non-LTE effects as well as resonant line transport. Until these are carried out, we cannot say whether the Fe column from accretion and spallation is in conflict with the observations. We also show that hydrogenic Fe might remain in the photosphere due to radiative levitation from the high burst flux.

We present a new mode of hydrogen burning on neutron stars (NSs) called diffusive nuclear burning (DNB). In DNB, the burning occurs in the exponentially suppressed tail of hydrogen that extends to the hotter regions of the envelope where protons are readily captured. Diffusive nuclear burning changes the compositional structure of the envelope on timescales $\sim 10^{2-4} {\rm yrs}$, much shorter than otherwise expected. This mechanism is applicable to the physics of young pulsars, millisecond radio pulsars (MSPs) and quiescent low mass X-ray binaries (LMXBs).

We calculate the rate of hydrogen burning for neutron stars (NSs) with hydrogen atmospheres and an underlying reservoir of nuclei capable of proton capture. This burning occurs in the exponentially suppressed diffusive tail of H that extends to the hotter depths of the envelope where protons are rapidly captured. This process, which we call diffusive nuclear burning (DNB), can change the H abundance at the NS photosphere on timescales as short as $10^{2-4}$ years. In the absence of diffusion, the hydrogen at the photosphere (where $T\approx 10^6 {\rm K}$ and $\rho\sim 0.1 {\rm g cm^{-2}}$) would last for far longer than a Hubble time. Our work impacts the understanding of the evolution of surface abundances of isolated NSs, which is important to their thermal spectrum and their effective temperature-core temperature relation. In this paper, we calculate the rate of H burning when the overall consumption rate is controlled by the nuclear timescales, rather than diffusion timescales. The immediate application is for H burning on millisecond radio pulsars and in quiescence for the accreting NS Cen X-4. We will apply this work to young radio pulsars and magnetars once we have incorporated the effects of strong $B>10^{12} {\rm G}$ magnetic fields.

The composition of the outer 100 m of a neutron star sets the heat flux that flows outwards from the core. For an accreting neutron star in an X-ray transient, the thermal quiescent flux depends sensitively on the amount of hydrogen and helium remaining on the surface after an accretion outburst and on the composition of the underlying ashes of previous H/He burning. Because H/He has a higher thermal conductivity, a larger mass of H/He implies a shallower thermal gradient through the low density envelope and hence a higher effective temperature for a given core temperature. The mass of residual H and He varies from outburst to outburst, so the thermal quiescent flux is variable even though the core temperature is constant for timescales < 10 000 yr. Heavy elements settle from a H/He envelope in a few hours; we therefore model the quiescent envelope as two distinct layers, H/He over heavier elements, and treat the mass of H/He as a free parameter. We find that the emergent thermal quiescent flux can vary by a factor of 2 to 3 between different quiescent epochs. The variation is more pronounced at lower interior temperatures, making systems with low quiescent luminosities and frequent outbursts, such as SAX J1808.4-3658, ideal candidates from which to observe this effect. We compute, for different ash compositions, the interior temperatures of Cen X-4, Aql X-1, and SAX J1808.4-3658. In the case of Aql X-1, the inferred high interior temperature suggests that neutrino cooling contributes to the neutron star's thermal balance.