Source author record

Philip Chang

Philip Chang appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

37works
13topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

37 published item(s)

preprint2022arXiv

Snowmass2021 Cosmic Frontier White Paper: Observational Facilities to Study Dark Matter

We present an overview of future observational facilities that will significantly enhance our understanding of the fundamental nature of dark matter. These facilities span a range of observational techniques including optical/near-infrared imaging and spectroscopy, measurements of the cosmic microwave background, pulsar timing, 21-cm observations of neutral hydrogen at high redshift, and the measurement of gravitational waves. Such facilities are a critical component of a multi-pronged experimental program to uncover the nature of dark matter, while often providing complementary measurements of dark energy, neutrino physics, and inflation.

preprint2021arXiv

Eclipse timing the Milky Way's gravitational potential

We show that a small, but \textit{measurable} shift in the eclipse mid-point time of eclipsing binary (EBs) stars of $\sim$ 0.1 seconds over a decade baseline can be used to directly measure the Galactic acceleration of stars in the Milky Way at $\sim$ kpc distances from the Sun. We consider contributions to the period drift rate from dynamical mechanisms other than the Galaxy's gravitational field, and show that the Galactic acceleration can be reliably measured using a sample of $\textit{Kepler}$ EBs with orbital and stellar parameters from the literature. Given the uncertainties on the formulation of tidal decay, our approach here is necessarily approximate, and the contribution from tidal decay is an upper limit assuming the stars are not tidally synchronized. We also use simple analytic relations to search for well-timed sources in the \textit{Kepler} field, and find $\sim$ 70 additional detached EBs with low eccentricities that have estimated timing precision better than 1 second. We illustrate the method with a prototypical, precisely timed EB using an archival \textit{Kepler} light curve and a modern synthetic \textit{HST} light curve (which provides a decade baseline). This novel method establishes a realistic possibility for obtaining fundamental Galactic parameters using eclipse timing to measure Galactic accelerations, along with other emerging new methods, including pulsar timing and extreme precision radial velocity observations. This acceleration signal grows quadratically with time. Therefore, given baselines established in the near-future for distant EBs, we can expect to measure the period drift in the future with space missions like \textit{JWST} and the \textit{Roman Space Telescope}.

preprint2020arXiv

General Relativistic Hydrodynamics on a Moving-mesh I: Static Spacetimes

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%

preprint2020arXiv

Simulation of a Compact Object with Outflows Moving Through a Gaseous 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.

preprint2020arXiv

The effect of diffusive nuclear burning in neutron star envelopes on cooling in accreting systems

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.

preprint2020arXiv

The growth of the longitudinal beam-plasma instability in the presence of an inhomogeneous background

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.

preprint2020arXiv

Time Dependent Radiation Hydrodynamics on a Moving Mesh

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.

preprint2019arXiv

Antlia2's role in driving the ripples in the outer gas disk of the Galaxy

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 $(μ_αcosδ, μ_δ) = (-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.

preprint2018arXiv

A Roadmap for HEP Software and Computing R&D for the 2020s

Particle physics has an ambitious and broad experimental programme for the coming decades. This programme requires large investments in detector hardware, either to build new facilities and experiments, or to upgrade existing ones. Similarly, it requires commensurate investment in the R&D of software to acquire, manage, process, and analyse the shear amounts of data to be recorded. In planning for the HL-LHC in particular, it is critical that all of the collaborating stakeholders agree on the software goals and priorities, and that the efforts complement each other. In this spirit, this white paper describes the R&D activities required to prepare for this software upgrade.

preprint2016arXiv

Bow Ties in the Sky I: The Angular Structure of Inverse Compton Gamma-ray Halos in the Fermi Sky

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.

preprint2016arXiv

Collapse in Self-gravitating Turbulent Fluids

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, $ρ(r,t ) \rightarrow ρ(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.

preprint2016arXiv

Galactoseismology: Discovery of a cluster of receding, variable halo stars

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.

preprint2016arXiv

The Linear Instability of Dilute Ultrarelativistic $e^{\pm}$ Pair Beams

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{ü}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.

preprint2015arXiv

Magnetized Moving Mesh Merger of a Carbon-Oxygen White Dwarf Binary

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.

preprint2015arXiv

Patchy blazar heating: diversifying the thermal history of the intergalactic medium

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.

preprint2015arXiv

Time Varying Dynamical Star Formation Rate

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.

preprint2014arXiv

An Eccentricity-Mass Relation for Galaxies from Tidally Disrupting Satellites

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.

preprint2014arXiv

Star Formation in Self-Gravitating Turbulent Fluids

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, $ρ(r,t)\toρ(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_*=ϕM_{\rm cl}(t/τ_{ff})^2$, with $ϕ$ a dimensionless number somewhat less than unity, $M_{\rm cl}$ the clump mass and $τ_{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.

preprint2014arXiv

The Effect of Nonlinear Landau Damping on Ultrarelativistic Beam Plasma Instabilities

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.

preprint2013arXiv

A Parameter-Space Study of Carbon-Oxygen White Dwarf Mergers

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.

preprint2013arXiv

Implications of Plasma Beam Instabilities for the Statistics of the Fermi Hard Gamma-ray Blazars and the Origin of the Extragalactic Gamma-Ray Background

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.

preprint2013arXiv

NGC1266 as a Local Candidate for Rapid Cessation of Star Formation

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.

preprint2012arXiv

The Cosmological Impact of Luminous TeV Blazars I: Implications of Plasma Instabilities for the Intergalactic Magnetic Field and Extragalactic Gamma-Ray Background

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.

preprint2012arXiv

The Cosmological Impact of Luminous TeV Blazars II: Rewriting the Thermal History of the Intergalactic Medium

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)

preprint2012arXiv

The Cosmological Impact of Luminous TeV Blazars III: Implications for Galaxy Clusters and the Formation of Dwarf Galaxies

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)

preprint2012arXiv

The Lyman-alpha forest in a blazar-heated Universe

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]

preprint2011arXiv

Dark Subhaloes and Disturbances in Extended HI Discs

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.

preprint2011arXiv

Finding Dwarf Galaxies From Their Tidal Imprints

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.

preprint2011arXiv

Properties of Carbon-Oxygen White Dwarf Merger Remnants

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.

preprint2011arXiv

Star Formation in Massive Clusters via Bondi Accretion

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 $τ_{\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_*=ε_{\rm ff} M_{\rm cl}/τ_{\rm cl}$, with $ε_{\rm ff} \approx 0.017$). However, after $\sim 2$ GMC free-fall times $τ_{\rm GMC}$, the clump accretion rate accelerates rapidly; formally, the clump can accrete the entire GMC in $\sim 3τ_{\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.

preprint2011arXiv

Warm Saturns: On the Nature of Rings around Extrasolar Planets that Reside Inside the Ice Line

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.

preprint2010arXiv

Diffusive Nuclear Burning of Helium on Neutron Stars

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 $α$ capturing material, accounting for thermal, mass defect, and Coulomb corrections on the stratification of material with the same zeroth order $μ_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.

preprint2010arXiv

Stellar Structure and Tests of Modified Gravity

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.

preprint2010arXiv

Sub-Chandrasekhar White Dwarf Mergers as the Progenitors of Type Ia Supernovae

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.

preprint2009arXiv

Migration of Star Clusters and Nuclear Rings

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.

preprint2009arXiv

Shock Breakout from Type Ia Supernova

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.

preprint2004arXiv

Magnetar Spindown, Hyper-Energetic Supernovae, and Gamma Ray Bursts

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.