Showing 1–50 of 79 results for author: Heyl, J S

Tests of Gravity Theories Using Supermassive Black Holes

Authors: Jeremy Sakstein, Bhuvnesh Jain, Jeremy S. Heyl, Lam Hui

Abstract: Scalar-tensor theories of gravity generally violate the strong equivalence principle, namely compact objects have a suppressed coupling to the scalar force, causing them to fall slower. A black hole is the extreme example where such a coupling vanishes, i.e. black hole has no scalar hair. Following earlier work, we explore observational scenarios for detecting strong equivalence principle violatio… ▽ More Scalar-tensor theories of gravity generally violate the strong equivalence principle, namely compact objects have a suppressed coupling to the scalar force, causing them to fall slower. A black hole is the extreme example where such a coupling vanishes, i.e. black hole has no scalar hair. Following earlier work, we explore observational scenarios for detecting strong equivalence principle violation, focusing on galileon gravity as an example. For galaxies in-falling towards galaxy clusters, the supermassive black hole can be offset from the galaxy center away from the direction of the cluster. Hence, well resolved images of galaxies around nearby clusters can be used to identify the displaced black hole via the star cluster bound to it. We show that this signal is accessible with imaging surveys, both ongoing ones such as the Dark Energy Survey, and future ground and space based surveys. Already, the observation of the central black hole in M~87 places new constraints on the galileon parameters, which we present here. $\mathcal{O}(1)$ matter couplings are disfavored for a large region of the parameter space. We also find a novel phenomenon whereby the black hole can escape the galaxy completely in less than one billion years. △ Less

Submitted 24 July, 2017; v1 submitted 7 April, 2017; originally announced April 2017.

Comments: Updated to match published version

Journal ref: ApJL, 2017, Vol. 844, L14

Polluting White Dwarfs with Perturbed Exo-Comets

Authors: Ilaria Caiazzo, Jeremy S. Heyl

Abstract: We present a model to account for the observed debris disks around young white dwarfs and the presence of metal-lines in their spectra. Stellar evolution models predict that the mass-loss on the AGB will be pulsed; furthermore, observations indicate that the bulk of the mass-loss occurs on the AGB. In this case, if the progenitors of the white dwarfs had remnants of planetary formation like the Su… ▽ More We present a model to account for the observed debris disks around young white dwarfs and the presence of metal-lines in their spectra. Stellar evolution models predict that the mass-loss on the AGB will be pulsed; furthermore, observations indicate that the bulk of the mass-loss occurs on the AGB. In this case, if the progenitors of the white dwarfs had remnants of planetary formation like the Sun's Oort cloud or the Kuiper Belt and a planet lying within that cloud or nearby, we find that up to 2% of the planetesimals will fall either into planet-crossing orbits or into chaotic regions after the mass-loss, depending on the location and mass of the planet (from Mars to Neptune). This yields a sufficient mass of comets that can be scattered toward the star, form a debris disk and pollute the atmosphere. △ Less

Submitted 27 April, 2017; v1 submitted 24 February, 2017; originally announced February 2017.

Comments: 10 pages, 5 figures, final version accepted by MNRAS

Globular cluster absolute ages from cooling brown dwarfs

Authors: Ilaria Caiazzo, Jeremy S. Heyl, Harvey Richer, Jason Kalirai

Abstract: Globular clusters are the oldest conglomerates of stars in our Galaxy and can be useful laboratories to test theories from stellar evolution to cosmology. In this paper, we present a new method to estimate the absolute age of a globular cluster from observations of its brown dwarfs. The transition region between the end of the main sequence and the brown dwarf regime is characterized by a dearth o… ▽ More Globular clusters are the oldest conglomerates of stars in our Galaxy and can be useful laboratories to test theories from stellar evolution to cosmology. In this paper, we present a new method to estimate the absolute age of a globular cluster from observations of its brown dwarfs. The transition region between the end of the main sequence and the brown dwarf regime is characterized by a dearth of objects as function of magnitude. The brightest of the cooling brown dwarfs is easily identified by an increase in density in the color-magnitude diagram as you go fainter in magnitudes, and these brightest brown dwarfs get fainter with age. By identifying the brightest brown dwarfs, it is thus possible to determine the age of a globular cluster within a 1 Gyr precision with four-sigma confidence. This new method, which is independent of current methods of age estimation and which does not rely on the knowledge of the cluster's distance from Earth, will become feasible thanks to the high spatial resolution and incredible infrared sensitivity of the James Webb Space Telescope. △ Less

Submitted 31 January, 2017; originally announced February 2017.

Using Galaxy Formation Simulations to optimise LIGO Follow-Up Observations

Authors: Elisa Antolini, Ilaria Caiazzo, Romeel Davé, Jeremy S. Heyl

Abstract: The recent discovery of gravitational radiation from merging black holes poses a challenge of how to organize the electromagnetic follow-up of gravitational-wave events as well as observed bursts of neutrinos. We propose a technique to select the galaxies that are most likely to host the event given some assumptions of whether the particular event is associated with recent star formation, low meta… ▽ More The recent discovery of gravitational radiation from merging black holes poses a challenge of how to organize the electromagnetic follow-up of gravitational-wave events as well as observed bursts of neutrinos. We propose a technique to select the galaxies that are most likely to host the event given some assumptions of whether the particular event is associated with recent star formation, low metallicity stars or simply proportional to the total stellar mass in the galaxy. We combine data from the 2-MASS Photometric Redshift Galaxy Catalogue with results from galaxy formation simulations to develop observing strategies that potentially reduce the area of sky to search by up to a factor of two relative to an unweighted search of galaxies, and a factor twenty to a search over the entire LIGO localization region. △ Less

Submitted 13 December, 2016; originally announced December 2016.

Comments: 6 pages, 4 figures, accepted by Monthly Notices

A Fast Algorithm for Finding Point Sources in the Fermi Data Stream: FermiFAST

Authors: Asha Asvathaman, Conor Omand, Alistair Barton, Jeremy S. Heyl

Abstract: We present a new and efficient algorithm for finding point sources in the photon event data stream from the Fermi Gamma-Ray Space Telescope, FermiFAST. The key advantage of FermiFAST is that it constructs a catalogue of potential sources very fast by arranging the photon data in a hierarchical data structure. Using this structure FermiFAST rapidly finds the photons that could have originated from… ▽ More We present a new and efficient algorithm for finding point sources in the photon event data stream from the Fermi Gamma-Ray Space Telescope, FermiFAST. The key advantage of FermiFAST is that it constructs a catalogue of potential sources very fast by arranging the photon data in a hierarchical data structure. Using this structure FermiFAST rapidly finds the photons that could have originated from a potential gamma-ray source. It calculates a likehihood ratio for the contribution of the potential source using the angular distribution of the photons within the region of interest. It can find within a few minutes the most significant half of the Fermi Third Point Source catalogue (3FGL) with nearly 80\% purity from the four years of data used to construct the catalogue. If a higher purity sample is desirable, one can achieve a sample that includes the most significant third of the Fermi 3FGL with only five percent of the sources unassociated with Fermi sources. Outside the galaxy plane, all but eight of the 580 FermiFAST detections are associated with 3FGL sources. And of these eight, six yield significant detections of greater than five sigma when a further binned likelihood analysis is performed. This software allows for rapid exploration of the Fermi data, simulation of the source detection to calculate the selection function of various sources and the errors in the obtained parameters of the sources detected. △ Less

Submitted 13 December, 2016; v1 submitted 10 June, 2016; originally announced June 2016.

Comments: 13 pages, 10 figures, accepted by Monthly Notices

Using the 2-MASS Photometric Redshift Survey to optimize LIGO Follow-Up Observations

Authors: Elisa Antolini, Jeremy S. Heyl

Abstract: The initial discovery of LIGO on 14 September 2015 was the inspiral merger and ring-down of the black hole binary at a distance of about 500~Mpc or a redshift of about 0.1. The search for electromagnetic counterparts for the inspiral of binary black holes is impeded by coarse initial source localisations and a lack of a compelling model for the counterpart; therefore, rapid electromagnetic follow-… ▽ More The initial discovery of LIGO on 14 September 2015 was the inspiral merger and ring-down of the black hole binary at a distance of about 500~Mpc or a redshift of about 0.1. The search for electromagnetic counterparts for the inspiral of binary black holes is impeded by coarse initial source localisations and a lack of a compelling model for the counterpart; therefore, rapid electromagnetic follow-up is required to understand the astrophysical context of these sources. Because astrophysical sources of gravitational radiation are likely to reside in galaxies, it would make sense to search first in regions where the LIGO-Virgo probability is large and where the density of galaxies is large as well. Under the assumption that the probability of a gravitational-wave event from a given region of space is proportional to the density of galaxies within the probed volume, one can calculate an improved localisation of the position of the source simply by multiplying the LIGO-Virgo skymap by the density of galaxies in the range of redshifts. We propose using the 2-MASS Photometric Redshift Galaxy Catalogue for this purpose and demonstrate that using it can dramatically reduce the search region for electromagnetic counterparts. △ Less

Submitted 14 July, 2016; v1 submitted 24 February, 2016; originally announced February 2016.

Comments: 8 pages, 6 figures, accepted for Monthly Notices

Eötvös Experiments with Supermassive Black Holes

Authors: Asha Asvathaman, Jeremy S. Heyl, Lam Hui

Abstract: By examining the locations of central black holes in two elliptical galaxies, M\,32 and M\,87, we derive constraints on the violation of the strong equivalence principle for purely gravitational objects, i.e. black holes, of less than about two-thirds, $η_N<0.68$ from the gravitational interaction of M\,87 with its neighbours in the Virgo cluster. Although M\,32 appears to be a good candidate for… ▽ More By examining the locations of central black holes in two elliptical galaxies, M\,32 and M\,87, we derive constraints on the violation of the strong equivalence principle for purely gravitational objects, i.e. black holes, of less than about two-thirds, $η_N<0.68$ from the gravitational interaction of M\,87 with its neighbours in the Virgo cluster. Although M\,32 appears to be a good candidate for this technique, the high concentration of stars near its centre substantially weakens the constraints. On the other hand, if a central black hole is found in NGC 205 or one of the other satellite ellipticals of M\,31, substantially better constraints could be obtained. In all cases the constraints could improve dramatically with better astrometry. △ Less

Submitted 13 December, 2016; v1 submitted 25 June, 2015; originally announced June 2015.

Comments: 7 pages, 3 figures, accepted by Monthly Notices

Statistical ages and the cooling rate of X-ray dim isolated neutron stars

Authors: Ramandeep Gill, Jeremy S. Heyl

Abstract: The cooling theory of neutron stars is corroborated by its comparison with observations of thermally emitting isolated neutron stars and accreting neutron stars in binary systems. An important ingredient for such an analysis is the age of the object, which, typically, is obtained from the spin-down history. This age is highly uncertain if the object's magnetic field varies appreciably over time. O… ▽ More The cooling theory of neutron stars is corroborated by its comparison with observations of thermally emitting isolated neutron stars and accreting neutron stars in binary systems. An important ingredient for such an analysis is the age of the object, which, typically, is obtained from the spin-down history. This age is highly uncertain if the object's magnetic field varies appreciably over time. Other age estimators, such as supernova remnant ages and kinematic ages, only apply to few handful of neutron stars. We conduct a population synthesis study of the nearby isolated thermal emitters and obtain their ages statistically from the observed luminosity function of these objects. We argue that a more sensitive blind scan of the galactic disk with the upcoming space telescopes can help to constrain the ages to higher accuracy. △ Less

Submitted 12 August, 2013; v1 submitted 4 May, 2013; originally announced May 2013.

Comments: 9 pages, 2 figures, 2 tables. Figures and tables updated. Added discussion on errors. Accepted for publication in MNRAS

k-d Match: A Fast Matching Algorithm for Sheared Stellar Samples

Authors: Jeremy S. Heyl

Abstract: This paper presents new and efficient algorithms for matching stellar catalogues where the transformation between the coordinate systems of the two catalagoues is unknown and may include shearing. Finding a given object whether a star or asterism from the first catalogue in the second is logarithmic in time rather than polynomial, yielding a dramatic speed up relative to a naive implementation. Bo… ▽ More This paper presents new and efficient algorithms for matching stellar catalogues where the transformation between the coordinate systems of the two catalagoues is unknown and may include shearing. Finding a given object whether a star or asterism from the first catalogue in the second is logarithmic in time rather than polynomial, yielding a dramatic speed up relative to a naive implementation. Both acceleration of the matching algorithm and the ability to solve for arbitrary affine transformations not only will allow the registration of stellar catalogues and images that are now impossible to use but also will find applications in machine vision and other imaging applications. △ Less

Submitted 2 May, 2013; v1 submitted 2 April, 2013; originally announced April 2013.

Comments: 5 pages, 4 figures (accepted for publication in Monthly Notices of the Royal Astronomical Society with minor revisions c.f. v1)

Diffractive Microlensing: A New Probe of the Local Universe

Authors: Jeremy S. Heyl

Abstract: Diffraction is important when nearby substellar objects gravitationally lens distant stars. If the wavelength of the observation is comparable to the Schwarzschild radius of lensing object, diffraction leaves an observable imprint on the lensing signature. The SKA may have sufficient sensitivity to detect the typical sources, giant stars in the bulge. The diffractive signatures in a lensing event… ▽ More Diffraction is important when nearby substellar objects gravitationally lens distant stars. If the wavelength of the observation is comparable to the Schwarzschild radius of lensing object, diffraction leaves an observable imprint on the lensing signature. The SKA may have sufficient sensitivity to detect the typical sources, giant stars in the bulge. The diffractive signatures in a lensing event break the degeneracies between the mass of the lens, its distance and proper motion. △ Less

Submitted 22 February, 2013; originally announced February 2013.

Comments: 3 pages, 2 figurs, submission for the Proceedings of the Thirteenth Marcel Grossman Meeting on General Relativity edited by Kjell Rosquist, Robert T Jantzen, Remo Ruffini World Scientific, Singapore, 2013

Magnetic Reconnection Instabilities in Soft-Gamma Repeaters

Authors: Jeremy S. Heyl, Ramandeep Gill

Abstract: We examine an external trigger mechanism that gives rise to the intense soft gamma-ray repeater (SGR) giant flares. Out of the three giant flares, two showcased the existence of a precursor, which we show to have had initiated the main flare. We develop a reconnection model based on the hypothesis that shearing motion of the footpoints causes the materialization of a Sweet-Parker current layer in… ▽ More We examine an external trigger mechanism that gives rise to the intense soft gamma-ray repeater (SGR) giant flares. Out of the three giant flares, two showcased the existence of a precursor, which we show to have had initiated the main flare. We develop a reconnection model based on the hypothesis that shearing motion of the footpoints causes the materialization of a Sweet-Parker current layer in the magnetosphere. The thinning of this macroscopic layer due to the development of an embedded super-hot turbulent current layer switches on the impulsive Hall reconnection, which powers the giant flare. We show that the thinning time is on the order of the pre-flare quiescent time. △ Less

Submitted 22 February, 2013; originally announced February 2013.

Comments: 4 pages, 2 figures, submission for the Proceedings of the Thirteenth Marcel Grossman Meeting on General Relativity edited by Kjell Rosquist, Robert T Jantzen, Remo Ruffini World Scientific, Singapore, 2013

Deep HST Imaging in NGC 6397: Stellar Dynamics

Authors: J. S. Heyl, H. Richer, J. Anderson, G. Fahlman, A. Dotter, J. Hurley, J. Kalirai, R. M. Rich, M. Shara, P. Stetson, K. A. Woodley, D. Zureko

Abstract: Multi-epoch observations with ACS on HST provide a unique and comprehensive probe of stellar dynamics within NGC 6397. We are able to confront analytic models of the globular cluster with the observed stellar proper motions. The measured proper motions probe well along the main sequence from 0.8 to below 0.1 M$_\odot$ as well as white dwarfs younger than one gigayear. The observed field lies just… ▽ More Multi-epoch observations with ACS on HST provide a unique and comprehensive probe of stellar dynamics within NGC 6397. We are able to confront analytic models of the globular cluster with the observed stellar proper motions. The measured proper motions probe well along the main sequence from 0.8 to below 0.1 M$_\odot$ as well as white dwarfs younger than one gigayear. The observed field lies just beyond the half-light radius where standard models of globular cluster dynamics (e.g. based on a lowered Maxwellian phase-space distribution) make very robust predictions for the stellar proper motions as a function of mass. The observed proper motions show no evidence for anisotropy in the velocity distribution; furthermore, the observations agree in detail with a straightforward model of the stellar distribution function. We do not find any evidence that the young white dwarfs have received a natal kick in contradiction with earlier results. Using the observed proper motions of the main-sequence stars, we obtain a kinematic estimate of the distance to NGC 6397 of $2.2^{+0.5}_{-0.7}$ kpc and a mass of the cluster of $1.1 \pm 0.1 \times 10^5 \mathrm{M}_\odot$ at the photometric distance of 2.53 kpc. One of the main-sequence stars appears to travel on a trajectory that will escape the cluster, yielding an estimate of the evaporation timescale, over which the number of stars in the cluster decreases by a factor of e, of about 3 Gyr. The proper motions of the youngest white dwarfs appear to resemble those of the most massive main-sequence stars, providing the first direct constraint on the relaxation time of the stars in a globular cluster of greater than or about 0.7 Gyr. △ Less

Submitted 4 October, 2012; v1 submitted 2 October, 2012; originally announced October 2012.

Comments: 25 pages, 20 figures, accepted for publication in Astrophysical Journal

Is Mira a magneto-dusty rotator?

Authors: Anand Thirumalai, Jeremy S. Heyl

Abstract: We investigate the possibility that a magnetic field may be present in the star $o-$Ceti (hereafter, Mira) and that the field plays a role in the star's mass loss. The model presented here is an application of an earlier derived theory that has been successfully employed for intermediate and high-mass evolved stars, and is now extended to the low-mass end. The modelling shows that it is possible t… ▽ More We investigate the possibility that a magnetic field may be present in the star $o-$Ceti (hereafter, Mira) and that the field plays a role in the star's mass loss. The model presented here is an application of an earlier derived theory that has been successfully employed for intermediate and high-mass evolved stars, and is now extended to the low-mass end. The modelling shows that it is possible to obtain a hybrid magnetohydrodynamic-dust-driven wind scenario for Mira, in which the role of a magnetic field in the equatorial plane of the star is dynamically important for producing a stellar wind. The wind velocity and the temperatures obtained from the model appear consistent with findings elsewhere. △ Less

Submitted 31 May, 2012; originally announced June 2012.

Comments: 5 pages, 2 figures

Journal ref: MNRAS, Vol. 430, pp. 1359-1368 (2013)

Proper Motions and Internal Dynamics in the Core of the Globular Cluster M71

Authors: Raminder Singh Samra, Harvey B. Richer, Jeremy S. Heyl, Ryan Goldsbury, Karun G. Thanjavur, Gordon A. Walker, Kristin A. Woodley

Abstract: We have used Gemini North together with the NIRI-ALTAIR adaptive optics imager in the H and K bands to explore the core of the Galactic globular cluster M71 (NGC 6838). We obtained proper motions for 217 stars and have resolved its internal proper motion dispersion. Using a 3.8 year baseline, the proper motion dispersion in the core is found to be 179 +/- 17 microarcsec/yr. We find no evidence of… ▽ More We have used Gemini North together with the NIRI-ALTAIR adaptive optics imager in the H and K bands to explore the core of the Galactic globular cluster M71 (NGC 6838). We obtained proper motions for 217 stars and have resolved its internal proper motion dispersion. Using a 3.8 year baseline, the proper motion dispersion in the core is found to be 179 +/- 17 microarcsec/yr. We find no evidence of anisotropy in the motions and no radial variation in the proper motions with respect to distance from the cluster center. We also set an upper limit on any central black hole to be ~150 Msun at 90% confidence level. △ Less

Submitted 23 March, 2012; originally announced March 2012.

Comments: 5 pages, 4 figures, accepted for publication in Astrophysical Journal Letters

A two-dimensional pseudospectral Hartree-Fock method for low-Z atoms in intense magnetic fields

Authors: Anand Thirumalai, Jeremy S. Heyl

Abstract: The energy levels of the first few low-lying states of helium and lithium atoms in intense magnetic fields up to $\approx 10^8-10^9$~T are calculated in this study. A pseudospectral method is employed for the computational procedure. The methodology involves computing the eigenvalues and eigenvectors of the generalized two-dimensional Hartree-Fock partial differential equations for these two- and… ▽ More The energy levels of the first few low-lying states of helium and lithium atoms in intense magnetic fields up to $\approx 10^8-10^9$~T are calculated in this study. A pseudospectral method is employed for the computational procedure. The methodology involves computing the eigenvalues and eigenvectors of the generalized two-dimensional Hartree-Fock partial differential equations for these two- and three-electron systems in a self-consistent manner. The method exploits the natural symmetries of the problem without assumptions of any basis functions for expressing the wave functions of the electrons or the commonly employed adiabatic approximation. It is seen that the results obtained here for a few of the most tightly bound states of each of the atoms, helium and lithium, are in good agreement with findings elsewhere. In this regard, we report new data for two new states of lithium that have not been studied thus far in the literature. It is also seen that the pseudospectral method employed here is considerably more economical, from a computational point of view, than previously employed methods such as a finite-element based approach. The key enabling advantage of the method described here is the short computational times which are on the order of seconds for obtaining accurate results for heliumlike systems. △ Less

Submitted 23 April, 2014; v1 submitted 15 March, 2012; originally announced March 2012.

Comments: 24 pages, 8 figures, 7 tables

Journal ref: Phys. Rev. A 89, 052522 (2014)

The Magnetised Bellows of Betelgeuse

Authors: Anand Thirumalai, Jeremy S. Heyl

Abstract: We present calculations for a magnetised hybrid wind model for Betelegeuse ($α- $Orionis). The model is a direct application of our previously derived theory, combining a canonical Weber-Davis (WD) stellar wind with dust grains in the envelope of an AGB star \citep[see][]{Thirumalai2010}. The resulting hybrid picture provides a mechanism for solving the problem of lifting stellar material up from… ▽ More We present calculations for a magnetised hybrid wind model for Betelegeuse ($α- $Orionis). The model is a direct application of our previously derived theory, combining a canonical Weber-Davis (WD) stellar wind with dust grains in the envelope of an AGB star \citep[see][]{Thirumalai2010}. The resulting hybrid picture provides a mechanism for solving the problem of lifting stellar material up from the photosphere \citep[e.g.][] {Harper2009,Guandalini2006,Jura1984} and into the circumstellar envelope. It also predicts wind velocities in agreement with current estimates. Our approach reveals that magnetic fields in supergiant stars like Betelgeuse \citep[see][]{Auriere2010}, may play a vital role in determining the nature of the stellar outflow and consequently, opens a new avenue of investigation in the field of hybrid stellar winds. △ Less

Submitted 6 February, 2012; v1 submitted 23 September, 2011; originally announced September 2011.

Comments: 13 pages, 4 figures

Journal ref: MNRAS (2012) Vol. 422 1272-1282

Constraining the photon-axion coupling constant with magnetic white dwarfs

Authors: Ramandeep Gill, Jeremy S. Heyl

Abstract: The light pseudoscalar particle, dubbed the axion, borne out of the Peccei-Quinn solution to the strong CP problem in QCD remains elusive. One avenue of inferring its existence is through its coupling to electromagnetic radiation. So far, laboratory experiments have dedicated all efforts to detect the axion in the mass range $10^{-6} < m_a < 10^{-3}$ eV with a photon-axion coupling strength… ▽ More The light pseudoscalar particle, dubbed the axion, borne out of the Peccei-Quinn solution to the strong CP problem in QCD remains elusive. One avenue of inferring its existence is through its coupling to electromagnetic radiation. So far, laboratory experiments have dedicated all efforts to detect the axion in the mass range $10^{-6} < m_a < 10^{-3}$ eV with a photon-axion coupling strength $g_{aγγ} < 10^{-10}GeV^{-1}$, where the limits are derived from astrophysical considerations. In this study, we present a novel way of constraining $g_{aγγ}$ by looking at the level of linear polarization in the radiation emerging from magnetic white dwarfs (mWDs). We find that photon-axion oscillations in WD magnetospheres can enhance the degree of linear polarization. Observing that most mWDs show only 5% linear polarization, we derive upper limits on $g_{aγγ}$ for different axion masses. △ Less

Submitted 14 September, 2011; v1 submitted 10 May, 2011; originally announced May 2011.

Comments: Accepted for publication in PRD. 10 pages, 5 figs. Main results are changed, new references are added, Results section is extended to incorporate the changes. Outlook section is also extended with new references. Fig. 2 and 4 are modified

A hybrid steady-state magnetohydrodynamic dust-driven stellar wind model for AGB stars

Authors: Anand Thirumalai, Jeremy S. Heyl

Abstract: We present calculations for a magnetised hybrid wind model for Asymptotic Giant Branch (AGB) stars. The model incorporates a canonical Weber-Davis (WD) stellar wind with dust grains in the envelope of an AGB star. The resulting hybrid picture preserves traits of both types of winds. It is seen that this combination requires that the dust-parameter ($Γ_{d}$) be less than unity in order to achieve a… ▽ More We present calculations for a magnetised hybrid wind model for Asymptotic Giant Branch (AGB) stars. The model incorporates a canonical Weber-Davis (WD) stellar wind with dust grains in the envelope of an AGB star. The resulting hybrid picture preserves traits of both types of winds. It is seen that this combination requires that the dust-parameter ($Γ_{d}$) be less than unity in order to achieve an outflow. The emergence of critical points in the wind changes the nature of the dust-driven outflow, simultaneously, the presence of a dust condensation radius changes the morphology of the magnetohydrodynamic (MHD) solutions for the wind. In this context, we additionally investigate the effect of having magnetic-cold spots on the equator of an AGB star and its implications for dust formation; which are seen to be consistent with previous findings. △ Less

Submitted 6 February, 2012; v1 submitted 10 June, 2010; originally announced June 2010.

Comments: 15 pages, 9 figures

Journal ref: MNRAS, Volume 409, Issue 4, pp. 1669-1681, (2010)

Cosmic Rays from Pulsars and Magnetars

Authors: Jeremy S. Heyl, Ramandeep Gill, Lars Hernquist

Abstract: We compare the expected abundance of cosmic-ray electrons and positrons from pulsars and magnetars. We assume that the distribution of infant pulsars and magnetars follows that of high-mass stars in the Milky Way and that the production rate of cosmic rays is proportional to the spin-down and magnetic-decay power of pulsars and magnetars, respectively. In combination with primary and secondary c… ▽ More We compare the expected abundance of cosmic-ray electrons and positrons from pulsars and magnetars. We assume that the distribution of infant pulsars and magnetars follows that of high-mass stars in the Milky Way and that the production rate of cosmic rays is proportional to the spin-down and magnetic-decay power of pulsars and magnetars, respectively. In combination with primary and secondary cosmic-ray leptons from other sources (especially supernova remnants), we find that both magnetars and pulsars can easily account for the observed cosmic-ray spectrum, in particular the dip seen by HESS at several TeV and the increase in positron fraction found by PAMELA. △ Less

Submitted 6 May, 2010; originally announced May 2010.

Comments: 5 pages, 3 figures, accepted to MNRAS

Diffractive Microlensing III: Astrometric Signatures

Authors: Jeremy S. Heyl

Abstract: Gravitational lensing is generally treated in the geometric optics limit; however, when the wavelength of the radiation approaches or exceeds the Schwarzschild radius of the lens, diffraction becomes important. Although the magnification generated by diffractive gravitational lensing is well understood, the astrometric signatures of diffractive microlensing are first derived in this paper along wi… ▽ More Gravitational lensing is generally treated in the geometric optics limit; however, when the wavelength of the radiation approaches or exceeds the Schwarzschild radius of the lens, diffraction becomes important. Although the magnification generated by diffractive gravitational lensing is well understood, the astrometric signatures of diffractive microlensing are first derived in this paper along with a simple closed-form bound for the astrometric shift. This simple bound yields the maximal shifts for substellar lenses in solar neighbourhood observed at 20~GHz, accessible to high sensitivity, high angular resolution radio telescopes such as the proposed Square Kilometre Array (SKA). △ Less

Submitted 4 October, 2010; v1 submitted 1 March, 2010; originally announced March 2010.

Comments: 6 pages, 3 figures, added calculations of expected rates, number of sources and lenses; version accepted to MNRAS

On the trigger mechanisms for SGR giant flares

Authors: Ramandeep Gill, Jeremy S. Heyl

Abstract: We examine two trigger mechanisms, one internal and the other external to the neutron star, that give rise to the intense soft gamma-ray repeater (SGR) giant flares. So far, three giant flares have been observed from the three out of the seven confirmed SGRs on March 5, 1979, August 27, 1998, and December 27, 2004. The last two events were found to be much more powerful than the first, and both sh… ▽ More We examine two trigger mechanisms, one internal and the other external to the neutron star, that give rise to the intense soft gamma-ray repeater (SGR) giant flares. So far, three giant flares have been observed from the three out of the seven confirmed SGRs on March 5, 1979, August 27, 1998, and December 27, 2004. The last two events were found to be much more powerful than the first, and both showcased the existence of a precursor, that we show to have had initiated the main flare. In the internal mechanism, we propose that the strongly wound up poloidal magnetic field develops tangential discontinuities and dissipates its torsional energy in heating the crust. The timescale for the instability to develop coincides with the duration of the quiescent state that followed the precursor. Alternatively, we develop a reconnection model based on the hypothesis that shearing motion of the footpoints causes the materialization of a Sweet-Parker current layer in the magnetosphere. The thinning of this macroscopic layer due to the development of an embedded super-hot turbulent current layer switches on the impulsive Hall reconnection, which powers the giant flare. Again, we show that the thinning time is on the order of the preflare quiescent time. This model naturally explains the origin of the observed nonthermal radiation during the flares, as well as the post flare radio afterglows. △ Less

Submitted 18 May, 2010; v1 submitted 18 February, 2010; originally announced February 2010.

Comments: 8 pages, 1 figure, version accepted by MNRAS. Changes: New references, added a section to introduction, added a paragraph to discussion. The results of this study remain unchanged.

Diffractive Microlensing II: Substellar Disk and Halo Objects

Authors: Jeremy S. Heyl

Abstract: Microlensing is generally studied in the geometric optics limit. However, diffraction may be important when nearby substellar objects lens occult distant stars. In particular the effects of diffraction become more important as the wavelength of the observation increases. Typically if the wavelength of the observation is comparable to the Schwarzschild radius of lensing object, diffraction leaves a… ▽ More Microlensing is generally studied in the geometric optics limit. However, diffraction may be important when nearby substellar objects lens occult distant stars. In particular the effects of diffraction become more important as the wavelength of the observation increases. Typically if the wavelength of the observation is comparable to the Schwarzschild radius of lensing object, diffraction leaves an observable imprint on the lensing signature. The commissioning of the Square Kilometre Array (SKA) over the next decade begs the question of whether it will become possible to follow up lensing events with radio observations because the SKA may have sufficient sensitivity to detect the typical sources, giant stars in the bulge. The detection of diffractive lensing in a lensing event would place unique constraints on the mass of the lens and its distance. In particular it would distinguish rapidly moving stellar mass lenses (e.g. neutron stars) from slowly moving substellar objects such freely floating planets. An analysis of the sensitivity of the SKA along with new simple closed-form estimates of the expected signal applied to local exemplars for stellar radio emission reveals that this effect can nearly be detected with the SKA. If the radio emission from bulge giants is stronger than expected, the SKA could detect the diffractive microlensing signature from Earth-like interstellar planets in the solar neighborhood. △ Less

Submitted 4 October, 2010; v1 submitted 16 February, 2010; originally announced February 2010.

Comments: 8 pages, 5 figures; added calculations of the event rate, number of sources and lenses; version accepted to MNRAS

Nonlinear Electromagnetic Waves in Magnetosphere of a Magnetar

Authors: Dan Mazur, Jeremy S. Heyl

Abstract: We compute electromagnetic wave propagation through the magnetosphere of a magnetar. The magnetosphere is modeled as the QED vacuum and a cold, strongly magnetized plasma. The background field and electromagnetic waves are treated nonperturbatively and can be arbitrarily strong. This technique is particularly useful for examining non-linear effects in propagating waves. Waves travelling through su… ▽ More We compute electromagnetic wave propagation through the magnetosphere of a magnetar. The magnetosphere is modeled as the QED vacuum and a cold, strongly magnetized plasma. The background field and electromagnetic waves are treated nonperturbatively and can be arbitrarily strong. This technique is particularly useful for examining non-linear effects in propagating waves. Waves travelling through such a medium typically form shocks; on the other hand we focus on the possible existence of waves that travel without evolving. Therefore, in order to examine the nonlinear effects, we make a travelling wave ansatz and numerically explore the resulting wave equations. We discover a class of solutions in a homogeneous plasma which are stabilized against forming shocks by exciting nonorthogonal components which exhibit strong nonlinear behaviour. These waves may be an important part of the energy transmission processes near pulsars and magnetars. △ Less

Submitted 10 November, 2010; v1 submitted 15 February, 2010; originally announced February 2010.

Comments: 8 pages, 6 figures, edited for clarity and references added, version accepted for publication by MNRAS

Dispersion Relations for Bernstein Waves in a Relativistic Pair Plasma

Authors: Ramandeep Gill, Jeremy S. Heyl

Abstract: A fully relativistic treatment of Bernstein waves in an electron-positron pair plasma has remained too formidable a task owing to the very complex nature of the problem. In this article, we perform contour integration of the dielectric response function and numerically compute the dispersion curves for a uniform, magnetized, relativistic electron-positron pair plasma. The behavior of the dispers… ▽ More A fully relativistic treatment of Bernstein waves in an electron-positron pair plasma has remained too formidable a task owing to the very complex nature of the problem. In this article, we perform contour integration of the dielectric response function and numerically compute the dispersion curves for a uniform, magnetized, relativistic electron-positron pair plasma. The behavior of the dispersion solution for several cases with different plasma temperatures is highlighted. In particular, we find two wave modes that exist only for large wavelengths and frequencies similar to the cyclotron frequency in a moderately relativistic pair plasma. The results presented here have important implications for the study of those objects where a hot magnetized electron-positron plasma plays a fundamental role in generating the observed radiation. △ Less

Submitted 21 August, 2009; v1 submitted 25 June, 2009; originally announced June 2009.

Comments: 8 pages, 8 figures, Accepted for publication by Phys. Rev. E with minor changes

Journal ref: Phys.Rev.E80:036407,2009

Polarized X-rays from Magnetized Neutron Stars

Authors: Dong Lai, Wynn C. G. Ho, Matthew van Adelsberg, Chen Wang, Jeremy S. Heyl

Abstract: We review the polarization properties of X-ray emission from highly magnetized neutron stars, focusing on emission from the stellar surfaces. We discuss how x-ray polarization can be used to constrain neutron star magnetic field and emission geometry, and to probe strong-field quantum electrodynamics and possibly constrain the properties of axions. We review the polarization properties of X-ray emission from highly magnetized neutron stars, focusing on emission from the stellar surfaces. We discuss how x-ray polarization can be used to constrain neutron star magnetic field and emission geometry, and to probe strong-field quantum electrodynamics and possibly constrain the properties of axions. △ Less

Submitted 12 October, 2009; v1 submitted 25 June, 2009; originally announced June 2009.

Comments: to appear in "X-ray Polarimetry: A New Window in Astrophysics", edited by R. Bellazzini, E. Costa, G. Matt and G. Tagliaferri (Cambridge University Press)

Pseudospectral methods for atoms in strong magnetic fields

Authors: Jeremy S. Heyl, Anand Thirumalai

Abstract: We present a new pseudospectral algorithm for the calculation of the structure of atoms in strong magnetic fields. We have verified this technique for one, two and three-electron atoms in zero magnetic fields against laboratory results and find typically better than one-percent accuracy. We further verify this technique against the state-of-the-art calculations of hydrogen, helium and lithium in… ▽ More We present a new pseudospectral algorithm for the calculation of the structure of atoms in strong magnetic fields. We have verified this technique for one, two and three-electron atoms in zero magnetic fields against laboratory results and find typically better than one-percent accuracy. We further verify this technique against the state-of-the-art calculations of hydrogen, helium and lithium in strong magnetic fields (up to about $2\times 10^{6}$ T) and find a similar level of agreement. The key enabling advantages of the algorithm are its simplicity (about 130 lines of commented code) and its speed (about $10^2-10^5$ times faster than finite-element methods to achieve similar accuracy). △ Less

Submitted 28 April, 2010; v1 submitted 27 February, 2009; originally announced March 2009.

Comments: 10 pages, version accepted to MNRAS

Constraining white-dwarf kicks in globular clusters : IV. Retarding Core Collapse

Authors: Jeremy S. Heyl, Matthew Penrice

Abstract: Observations of white dwarfs in the globular clusters NGC 6397 and Omega Centauri indicate that these stars may get a velocity kick during their time as giants. If the mass loss while on the asymptotic giant branch is slightly asymmetric, the resulting white dwarf could be born with such a velocity kick. These energetic white dwarfs will impart their excess energy on other stars as they travel… ▽ More Observations of white dwarfs in the globular clusters NGC 6397 and Omega Centauri indicate that these stars may get a velocity kick during their time as giants. If the mass loss while on the asymptotic giant branch is slightly asymmetric, the resulting white dwarf could be born with such a velocity kick. These energetic white dwarfs will impart their excess energy on other stars as they travel through the cluster. A Monte-Carlo simulation of the white-dwarfs kicks combined with estimate of the phase-space diffusion of the white dwarfs reveals that as the white dwarfs equilibrate, they lose most of their energy in the central region of the cluster. They could possibly augment the effect of binaries, delaying core collapse or increasing the size of globular cluster cores. △ Less

Submitted 29 May, 2009; v1 submitted 13 January, 2009; originally announced January 2009.

Comments: 4 pages, 2 figures, changes to reflect version accepted to MNRAS

Hydrogen and Helium atoms in strong magnetic fields

Authors: Anand Thirumalai, Jeremy S. Heyl

Abstract: The energy levels of hydrogen and helium atoms in strong magnetic fields are calculated in this study. The current work contains estimates of the binding energies of the first few low-lying states of these systems that are improvements upon previous estimates. The methodology involves computing the eigenvalues and eigenvectors of the generalized two-dimensional Hartree-Fock partial differential eq… ▽ More The energy levels of hydrogen and helium atoms in strong magnetic fields are calculated in this study. The current work contains estimates of the binding energies of the first few low-lying states of these systems that are improvements upon previous estimates. The methodology involves computing the eigenvalues and eigenvectors of the generalized two-dimensional Hartree-Fock partial differential equations for these one- and two-electron systems in a self-consistent manner. The method described herein is applicable to calculations of atomic structure in magnetic fields of arbitrary strength as it exploits the natural symmetries of the problem without assumptions of any basis functions for expressing the wave functions of the electrons or the commonly employed adiabatic approximation. The method is found to be readily extendable to systems with more than two electrons. △ Less

Submitted 6 February, 2012; v1 submitted 18 June, 2008; originally announced June 2008.

Comments: 15 pages, 6 figures

Journal ref: Phys. Rev. A 79, 012514 (2009)

Constraining white-dwarf kicks in globular clusters : III. Cluster Heating

Authors: Jeremy S. Heyl

Abstract: Recent observations of white dwarfs in globular clusters indicate that these stars may get a velocity kick during their time as giants. This velocity kick could originate naturally if the mass loss while on the asymptotic giant branch is slightly asymmetric. The kicks may be large enough to dramatically change the radial distribution of young white dwarfs, giving them larger energies than other… ▽ More Recent observations of white dwarfs in globular clusters indicate that these stars may get a velocity kick during their time as giants. This velocity kick could originate naturally if the mass loss while on the asymptotic giant branch is slightly asymmetric. The kicks may be large enough to dramatically change the radial distribution of young white dwarfs, giving them larger energies than other stars in the cluster. As these energetic white dwarfs travel through the cluster they can impart their excess energy on the other stars in the cluster. This new heat source for globular clusters is expected to be largest during the clusters' youth. △ Less

Submitted 31 March, 2009; v1 submitted 18 March, 2008; originally announced March 2008.

Comments: 4 pages, 1 figure (changes to reflect published version)

Neutron Star Crustal Mass Fractions

Authors: Kelsey L. Hoffman, Jeremy S. Heyl

Abstract: We are investigating mass fractions on the crust of a neutron star which would remain after one year of cooling. We use cooling curves corresponding with various densities, or depths, of the neutron star just after its formation. We assume the modified Urca process dominates the energy budget of the outer layers of the star in order to calculate the temperature of the neutron star as a function… ▽ More We are investigating mass fractions on the crust of a neutron star which would remain after one year of cooling. We use cooling curves corresponding with various densities, or depths, of the neutron star just after its formation. We assume the modified Urca process dominates the energy budget of the outer layers of the star in order to calculate the temperature of the neutron star as a function of time. Using a nuclear reaction network up to technetium, we calculate how the distribution of nuclei quenches at various depths of the neutron star crust. The initial results indicate that $^{28}$Si is the lightest isotope to be optically thick on the surface after one year of cooling. △ Less

Submitted 1 October, 2007; originally announced October 2007.

Comments: Comments: 40 Years of Pulsars conference: Millisecond Pulsars, Magnetars and More. McGill University, Montreal, Canada, August 12-17, 2007. 3 pages 4 figures

Journal ref: AIP Conf.Proc.983:388-390,2008

Constraining white-dwarf kicks in globular clusters : II. Observational Significance

Authors: Jeremy S. Heyl

Abstract: If the winds of an asymptotic-giant-branch stars are sufficiently strong are slightly asymmetric, they can alter the star's trajectory through a globular cluster; therefore, if these winds are asymmetric, one would expect young white dwarfs to be less radially concentrated than either their progenitors or older white dwarfs in globular clusters. This latter effect has recently been observed. Add… ▽ More If the winds of an asymptotic-giant-branch stars are sufficiently strong are slightly asymmetric, they can alter the star's trajectory through a globular cluster; therefore, if these winds are asymmetric, one would expect young white dwarfs to be less radially concentrated than either their progenitors or older white dwarfs in globular clusters. This latter effect has recently been observed. Additionally the young white dwarfs should have larger typical velocities than their progenitors. After phase mixing this latter effect is vastly diminished relative to the changes in the spatial distribution of young white dwarfs with kicks, so it is more difficult to detect than the change in the spatial distribution. The most powerful kinematic signature is the change in the eccentricity of the orbits that is revealed through the distribution of the position angles of proper motion. △ Less

Submitted 19 September, 2007; originally announced September 2007.

Comments: 5 pages, 4 figures, submitted to MNRAS letters

Orbital evolution with white-dwarf kicks

Authors: Jeremy S. Heyl

Abstract: Recent observations of white dwarfs in globular clusters indicate that these stars may get a velocity kick during their time as giants. This velocity kick could originate naturally if the mass loss while on the asymptotic giant branch is slightly asymmetric. If white dwarfs get a kick comparable to the orbital velocity of the binary, the initial Runge-Lenz vector (eccentricity vector) of the orb… ▽ More Recent observations of white dwarfs in globular clusters indicate that these stars may get a velocity kick during their time as giants. This velocity kick could originate naturally if the mass loss while on the asymptotic giant branch is slightly asymmetric. If white dwarfs get a kick comparable to the orbital velocity of the binary, the initial Runge-Lenz vector (eccentricity vector) of the orbit is damped to be replaced by a component pointing toward the cross product of the initial angular momentum and the force. The final eccentricity may be of order unity and if the kick is sufficiently large, the system may be disrupted. These results may have important ramifications for the evolution of binary stars and planetary systems △ Less

Submitted 10 September, 2007; v1 submitted 6 June, 2007; originally announced June 2007.

Comments: 6 pages, 2 figures, added many intermediate equations, version accepted by MNRAS

Constraining white-dwarf kicks in globular clusters

Authors: Jeremy S. Heyl

Abstract: The wind of an asymptotic-giant-branch stars is sufficiently strong that if it is slightly asymmetric, it can propel the star outside of the open cluster of its birth or significantly alter its trajectory through a globular cluster; therefore, if these stellar winds are asymmetric, one would expect a deficit of white dwarfs of all ages in open clusters and for young white dwarfs to be less radia… ▽ More The wind of an asymptotic-giant-branch stars is sufficiently strong that if it is slightly asymmetric, it can propel the star outside of the open cluster of its birth or significantly alter its trajectory through a globular cluster; therefore, if these stellar winds are asymmetric, one would expect a deficit of white dwarfs of all ages in open clusters and for young white dwarfs to be less radially concentrated than either their progenitors or older white dwarfs in globular clusters. This latter effect has recently been observed. Hence, detailed studies of the radial distribution of young white dwarfs in globular clusters could provide a unique probe of mass loss on the asymptotic giant branch and during the formation of planetary nebulae both as a function of metallicity and a limited range of stellar mass. △ Less

Submitted 7 September, 2007; v1 submitted 6 June, 2007; originally announced June 2007.

Comments: 4 pages, 3 figures, changes to reflect version accepted for MNRAS Letters

Using long-term transit timing to detect terrestrial planets

Authors: Jeremy S. Heyl, Brett J. Gladman

Abstract: We propose that the presence of additional planets in extrasolar planetary systems can be detected by long-term transit timing studies. If a transiting planet is on an eccentric orbit then the presence of another planet causes a secular advance of the transiting planet's pericenter over and above the effect of general relativity. Although this secular effect is impractical to detect over a small… ▽ More We propose that the presence of additional planets in extrasolar planetary systems can be detected by long-term transit timing studies. If a transiting planet is on an eccentric orbit then the presence of another planet causes a secular advance of the transiting planet's pericenter over and above the effect of general relativity. Although this secular effect is impractical to detect over a small number of orbits, it causes long-term differences in when future transits occur, much like the long-term decay observed in pulsars. Measuring this transit-timing delay would thus allow the detection of either one or more additional planets in the system or the first measurements of non-zero oblateness ($J_2$) of the central stars. △ Less

Submitted 7 March, 2007; v1 submitted 10 October, 2006; originally announced October 2006.

Comments: 10 pages, 5 figures, accepted by Monthly Notices, updated to reflect accepted version

Journal ref: Mon.Not.Roy.Astron.Soc.377:1511-1519,2007

Ionization and dissociation equilibrium in strongly-magnetized helium atmosphere

Authors: Kaya Mori, Jeremy S. Heyl

Abstract: Recent observations and theoretical investigations of neutron stars indicate that their atmospheres consist not of hydrogen or iron but possibly other elements such as helium. We calculate the ionization and dissociation equilibrium of helium in the conditions found in the atmospheres of magnetized neutron stars. For the first time this investigation includes the internal degrees of freedom of t… ▽ More Recent observations and theoretical investigations of neutron stars indicate that their atmospheres consist not of hydrogen or iron but possibly other elements such as helium. We calculate the ionization and dissociation equilibrium of helium in the conditions found in the atmospheres of magnetized neutron stars. For the first time this investigation includes the internal degrees of freedom of the helium molecule. We found that at the temperatures and densities of neutron star atmospheres the rotovibrational excitations of helium molecules are populated. Including these excitations increases the expected abundance of molecules by up to two orders of magnitude relative to calculations that ignore the internal states of the molecule; therefore, if the atmospheres of neutron stars indeed consist of helium, helium molecules and possibly polymers will make the bulk of the atmosphere and leave signatures on the observed spectra from neutron stars. We applied our calculation to nearby radio-quiet neutron stars with B_{dipole}~10^{13}-10^{14} G. If helium comprises their atmospheres, our study indicates that isolated neutron stars with T_{BB}~10^6 K such as RXJ0720.4-3125 and RXJ1605.3+3249 will have He^+ ions predominantly, while isolated neutron stars with lower temperature (T_{BB}~5x10^5 K) such as RXJ1856.5-3754 and RXJ0420.0-5022 will have some fraction of helium molecules. If helium molecules are abundant, their spectroscopic signatures may be detected in the optical, UV and X-ray band. (Abridged) △ Less

Submitted 23 February, 2007; v1 submitted 9 October, 2006; originally announced October 2006.

Comments: 12 pages, 12 figures, accepted for publication in MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.376:895-906,2007

QED can explain the non-thermal emission from SGRs and AXPs : Variability

Authors: Jeremy S. Heyl

Abstract: Owing to effects arising from quantum electrodynamics (QED), magnetohydrodynamical fast modes of sufficient strength will break down to form electron-positron pairs while traversing the magnetospheres of strongly magnetised neutron stars. The bulk of the energy of the fast mode fuels the development of an electron-positron fireball. However, a small, but potentially observable, fraction of the e… ▽ More Owing to effects arising from quantum electrodynamics (QED), magnetohydrodynamical fast modes of sufficient strength will break down to form electron-positron pairs while traversing the magnetospheres of strongly magnetised neutron stars. The bulk of the energy of the fast mode fuels the development of an electron-positron fireball. However, a small, but potentially observable, fraction of the energy ($\sim 10^{33}$ ergs) can generate a non-thermal distribution of electrons and positrons far from the star. This paper examines the cooling and radiative output of these particles. Small-scale waves may produce only the non-thermal emission. The properties of this non-thermal emission in the absence of a fireball match those of the quiescent, non-thermal radiation recently observed non-thermal emission from several anomalous X-ray pulsars and soft-gamma repeaters. Initial estimates of the emission as a function of angle indicate that the non-thermal emission should be beamed and therefore one would expect this emission to be pulsed as well. According to this model the pulsation of the non-thermal emission should be between 90 and 180 degrees out of phase from the thermal emission from the stellar surface. △ Less

Submitted 25 August, 2006; originally announced August 2006.

Comments: 7 pages, 5 figures, to appear in the proceedings of the conference "Isolated Neutron Stars: from the Interior to the Surface" (April 2006, London), eds. D. Page, R. Turolla, & S. Zane, Astrophysics & Space Science

Journal ref: Astrophys.SpaceSci.308:101-107,2007

See a Black Hole on a Shoestring

Authors: Jeremy S. Heyl

Abstract: The modes of vibration of hanging and partially supported strings provide useful analogies to scalar fields travelling through spacetimes that admit conformally flat spatial sections. This wide class of spacetimes includes static, spherically symmetric spacetimes. The modes of a spacetime where the scale factor depends as a power-law on one of the coordinates provide a useful starting point and… ▽ More The modes of vibration of hanging and partially supported strings provide useful analogies to scalar fields travelling through spacetimes that admit conformally flat spatial sections. This wide class of spacetimes includes static, spherically symmetric spacetimes. The modes of a spacetime where the scale factor depends as a power-law on one of the coordinates provide a useful starting point and yield a new classification of these spacetimes on the basis of the shape of the string analogue. The family of corresponding strings follow a family of curves related to the cycloid, denoted here as hypercycloids (for reasons that will become apparent). Like the spacetimes that they emulate these strings exhibit horizons, typically at their bottommost points where the string tension vanishes; therefore, hanging strings may provide a new avenue for the exploration of the quantum mechanics of horizons. △ Less

Submitted 15 September, 2006; v1 submitted 16 February, 2006; originally announced February 2006.

Comments: 5 pages, 1 figure, extensive changes to refect version accepted to PRD

Journal ref: Phys.Rev. D74 (2006) 064029

The Long-Term Future of Space Travel

Authors: Jeremy S. Heyl

Abstract: The fact that we apparently live in an accelerating universe places limitations on where humans might visit. If the current energy density of the universe is dominated by a cosmological constant, a rocket could reach a galaxy observed today at a redshift of 1.7 on a one-way journey or merely 0.65 on a round trip. Unfortunately these maximal trips are impractical as they require an infinite prope… ▽ More The fact that we apparently live in an accelerating universe places limitations on where humans might visit. If the current energy density of the universe is dominated by a cosmological constant, a rocket could reach a galaxy observed today at a redshift of 1.7 on a one-way journey or merely 0.65 on a round trip. Unfortunately these maximal trips are impractical as they require an infinite proper time to traverse. However, calculating the rocket trajectory in detail shows that a rocketeer could nearly reach such galaxies within a lifetime (a long lifetime admittedly -- about 100 years). For less negative values of $w$ the maximal redshift increases becoming infinite for $w\geq -1/3$. △ Less

Submitted 8 November, 2005; v1 submitted 9 September, 2005; originally announced September 2005.

Comments: 5 pages, 3 figures, minor changes to reflect version accepted to PRD

Journal ref: Phys.Rev.D72:107302,2005

Magnetars

Authors: Jeremy S. Heyl

Abstract: Ultramagnetized neutron stars or magnetars are magnetically powered neutron stars. Their strong magnetic fields dominate the physical processes in their crusts and their surroundings. The past few years have seen several advances in our theoretical and observational understanding of these objects. In spite of a surfeit of observations, their spectra are still poorly understood. I will discuss th… ▽ More Ultramagnetized neutron stars or magnetars are magnetically powered neutron stars. Their strong magnetic fields dominate the physical processes in their crusts and their surroundings. The past few years have seen several advances in our theoretical and observational understanding of these objects. In spite of a surfeit of observations, their spectra are still poorly understood. I will discuss the emission from strongly magnetized condensed matter surfaces of neutron stars, recent advances in our expectations of the surface composition of magnetars and a model for the non-thermal emission from these objects. △ Less

Submitted 4 April, 2005; originally announced April 2005.

Comments: will be appreciated. 12 pages, 5 figures, XXII Texas Symposium on Relatistic Astrophysics

Report number: 0036

R-Modes on Rapidly Rotating, Relativistic Stars : II. Blackbody Emission

Authors: Jeremy S. Heyl

Abstract: Rossby waves (or r-modes) on the surface of a neutron star have become a leading model for the oscillations observed during the tail of Type-I X-ray bursts. Their frequency evolution matches well with the observed frequency drifts of the oscillations in the bursts, and the burning appears to excite these waves quite naturally. This paper addresses the detailed shape of the expected flux profiles… ▽ More Rossby waves (or r-modes) on the surface of a neutron star have become a leading model for the oscillations observed during the tail of Type-I X-ray bursts. Their frequency evolution matches well with the observed frequency drifts of the oscillations in the bursts, and the burning appears to excite these waves quite naturally. This paper addresses the detailed shape of the expected flux profiles from r-modes on neutron stars as a function of energy. R-modes naturally account for both the small amplitude of the observed oscillations and their lack of harmonic content. However, the model predicts that the oscillation at higher energies leads the lower energy variation. The observed oscillations have the opposite trend which possibly indicates that the higher energy photons are upscattered in the plasma surrounding the neutron star and therefore delayed. △ Less

Submitted 24 February, 2005; originally announced February 2005.

Comments: 15 pages, 7 figures, submitted to MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc. 361 (2005) 504-510

The High-Energy Polarization-Limiting Radius of Neutron Star Magnetospheres II -- Magnetized Hydrogen Atmospheres

Authors: Jeremy S. Heyl, Don Lloyd, Nir J. Shaviv

Abstract: In the presence of strong magnetic fields, the vacuum becomes a birefringent medium. We show that this QED effect couples the direction of the polarization of photons leaving the NS surface, to the direction of the magnetic field along the ray's path. We analyze the consequences that this effect has on aligning the polarization vectors to generate large net polarizations, while considering therm… ▽ More In the presence of strong magnetic fields, the vacuum becomes a birefringent medium. We show that this QED effect couples the direction of the polarization of photons leaving the NS surface, to the direction of the magnetic field along the ray's path. We analyze the consequences that this effect has on aligning the polarization vectors to generate large net polarizations, while considering thermal radiation originating from a thermal hydrogen atmosphere. Counter to previous predictions, we show that the thermal radiation should be highly polarized even in the optical. When detected, this polarization will be the first demonstration of vacuum birefringence. It could be used as a tool to prove the high magnetic field nature of AXPs and it could also be used to constrain physical NS parameters, such as $R/M$, to which the net polarization is sensitive. △ Less

Submitted 17 February, 2005; originally announced February 2005.

Comments: 11 pages, 3 figures, submitted to MNRAS

A QED Model for Non-thermal Emission from SGRs and AXPs

Authors: Jeremy S. Heyl, Lars Hernquist

Abstract: Previously, we showed that, owing to effects arising from quantum electrodynamics (QED), magnetohydrodynamic fast modes of sufficient strength will break down to form electron-positron pairs while traversing the magnetospheres of strongly magnetised neutron stars. The bulk of the energy of the fast mode fuels the development of an electron-positron fireball. However, a small, but potentially obs… ▽ More Previously, we showed that, owing to effects arising from quantum electrodynamics (QED), magnetohydrodynamic fast modes of sufficient strength will break down to form electron-positron pairs while traversing the magnetospheres of strongly magnetised neutron stars. The bulk of the energy of the fast mode fuels the development of an electron-positron fireball. However, a small, but potentially observable, fraction of the energy ($\sim 10^{33}$ ergs) can generate a non-thermal distribution of electrons and positrons far from the star. In this paper, we examine the cooling and radiative output of these particles. We also investigate the properties of non-thermal emission in the absence of a fireball to understand the breakdown of fast modes that do not yield an optically thick pair plasma. This quiescent, non-thermal radiation associated with fast mode breakdown may account for the recently observed non-thermal emission from several anomalous X-ray pulsars and soft-gamma repeaters. △ Less

Submitted 17 February, 2005; originally announced February 2005.

Comments: 14 pages, 2 figures, submitted to MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc. 362 (2005) 777-783

Magnetospheric Birefringence Induces Unique Polarization Signatures in Neutron-Star Spectra

Authors: R. M. Shannon, Jeremy S. Heyl

Abstract: We study the propagation of polarization light through the magnetosphere of neutron stars. At intermediate frequencies (the optical through the infrared), both the birefringence induced by the plasma and by quantumelectrodynamics influence the observation polarization of radiation from the surface of the neutron star. Because these two processes compete in this regime, we find that polarization… ▽ More We study the propagation of polarization light through the magnetosphere of neutron stars. At intermediate frequencies (the optical through the infrared), both the birefringence induced by the plasma and by quantumelectrodynamics influence the observation polarization of radiation from the surface of the neutron star. Because these two processes compete in this regime, we find that polarization observations can constrain the properties of the neutron-star magnetosphere, specifically the total charge density. We calculate both the phase-resolved and the phase-averaged polarization signatures induced by magnetospheric birefringence. △ Less

Submitted 18 February, 2006; v1 submitted 22 October, 2004; originally announced October 2004.

Comments: 4 pages, 5 figures, changes to reflect accepted version

Journal ref: Mon.Not.Roy.Astron.Soc.368:1377-1380,2006

A QED Model for the Origin of Bursts from SGRs and AXPs

Authors: Jeremy S. Heyl, Lars Hernquist

Abstract: We propose a model to account for the bursts from soft gamma repeaters (SGRs) and anomalous X-ray pulsars (AXPs) in which quantum electrodynamics (QED) plays a vital role. In our theory, that we term ``fast-mode breakdown,'' magnetohydrodynamic (MHD) waves that are generated near the surface of a neutron star and propagate outward through the magnetosphere will be modified by the polarization of… ▽ More We propose a model to account for the bursts from soft gamma repeaters (SGRs) and anomalous X-ray pulsars (AXPs) in which quantum electrodynamics (QED) plays a vital role. In our theory, that we term ``fast-mode breakdown,'' magnetohydrodynamic (MHD) waves that are generated near the surface of a neutron star and propagate outward through the magnetosphere will be modified by the polarization of the vacuum. For neutron star magnetic fields $B_\rmscr{NS} \gtrsim B_\rmscr{QED} \approx 4.4 \times 10^{13}$ G, the interaction of the wave fields with the vacuum produces non-linearities in fast MHD waves that can steepen in a manner akin to the growth of hydrodynamic shocks. Under certain conditions, fast modes can develop field discontinuities on scales comparable to an electron Compton wavelength, at which point the wave energy will be dissipated through electron-positron pair production. We show that this process operates if the magnetic field of the neutron star is sufficiently strong and the ratio of the wavelength of the fast mode to its amplitude is sufficiently small, in which case the wave energy will be efficiently converted into an extended pair-plasma fireball. The radiative output from this fireball will consist of hard X-rays and soft $γ$-rays, with a spectrum similar to those seen in bursts from SGRs and AXPs (abridged). △ Less

Submitted 31 August, 2004; v1 submitted 23 December, 2003; originally announced December 2003.

Comments: 28 pages, 5 figures, accepted by ApJ, minor changes to reflect accepted version

Journal ref: Astrophys.J.618:463-473,2004

Rotational Evolution of Protoneutron Stars

Authors: Yefei Yuan, Jeremy S. Heyl

Abstract: We study the rotational evolution of a protoneutron star with hyperons and nucleons or solely nucleons in its core due to the escape of the trapped neutrinos. It is found that at the early stage of its evolution, the stellar crust contracts significantly, consequently the star spins up. At the late stage, a the protoneutron star with hyperons, it keeps shrinking and spinning up till all the trap… ▽ More We study the rotational evolution of a protoneutron star with hyperons and nucleons or solely nucleons in its core due to the escape of the trapped neutrinos. It is found that at the early stage of its evolution, the stellar crust contracts significantly, consequently the star spins up. At the late stage, a the protoneutron star with hyperons, it keeps shrinking and spinning up till all the trapped neutrinos escape. Consequently, the distribution of the stellar initial spin periods is skewed toward shorter periods. For a protoneutron star with only nucleons, the expansion of its core dominates, and the stellar rotation slows down. After the neutrinos escape, the range of the spin periods is narrower than the initial one, but the distribution is still nearly uniform. If the hyperonic star is metastable, its rotational frequency accelerates distinguishedly before it collapses to a black hole. △ Less

Submitted 6 May, 2003; originally announced May 2003.

Comments: 28 pages, 12 figures

Thermonuclear Stability of Material Accreting onto a Neutron Star

Authors: Ramesh Narayan, Jeremy S. Heyl

Abstract: We present a global linear stability analysis of nuclear fuel accumulating on the surface of an accreting neutron star and we identify the conditions under which thermonuclear bursts are triggered. The analysis reproduces all the recognized regimes of hydrogen and helium bursts, and in addition shows that at high accretion rates, near the limit of stable burning, there is a regime of ``delayed m… ▽ More We present a global linear stability analysis of nuclear fuel accumulating on the surface of an accreting neutron star and we identify the conditions under which thermonuclear bursts are triggered. The analysis reproduces all the recognized regimes of hydrogen and helium bursts, and in addition shows that at high accretion rates, near the limit of stable burning, there is a regime of ``delayed mixed bursts'' which is distinct from the more usual ``prompt mixed bursts.'' In delayed mixed bursts, a large fraction of the fuel is burned stably before the burst is triggered. Bursts thus have longer recurrence times, but at the same time have somewhat smaller fluences. Therefore, the parameter alpha, which measures the ratio of the energy released via accretion to that generated through nuclear reactions in the burst, is up to an order of magnitude larger than for prompt bursts. This increase in alpha near the threshold of stable burning has been seen in observations. We explore a wide range of mass accretion rates, neutron star radii and core temperatures, and calculate a variety of burst properties. From a preliminary comparison with data, we suggest that bursting neutron stars may have hot cores, with T_{core} >~ 10^{7.5} K, consistent with interior cooling via the modified URCA or similar low-efficiency process, rather than T_{core} ~ 10^7 K, as expected for the direct URCA process. There is also an indication that neutron star radii are somewhat small <~ 10 km. Both of these conclusions need to be confirmed by comparing more careful calculations with better data. △ Less

Submitted 22 August, 2003; v1 submitted 19 March, 2003; originally announced March 2003.

Comments: 67 pages, 19 figures, final version to appear in The Astrophysical Journal, vol. 599, no. 1, Dec. 10, 2003

Accretion limits the compactness of static stars

Authors: Jeremy S. Heyl

Abstract: General relativity limits the compactness of static stars. If the pressure of the fluid is positive and the density decreases with distance from the center, the value of the circumferential radius of the star must be greater than (9/4) G M/c^2, or equivalently the redshift of the surface must be less than two. If constraints on the equation of state of the material are relaxed, general relativit… ▽ More General relativity limits the compactness of static stars. If the pressure of the fluid is positive and the density decreases with distance from the center, the value of the circumferential radius of the star must be greater than (9/4) G M/c^2, or equivalently the redshift of the surface must be less than two. If constraints on the equation of state of the material are relaxed, general relativity alone does not restrict the redshift of a static stellar surface. However, because black hole candidates in the universe generally accrete material from their environs, the process of accretion provides upper limits on the redshift of a astrophysical black-hole candidates. △ Less

Submitted 2 March, 2003; originally announced March 2003.

Comments: 7 pages

The Synoptic Swift Synergy -- Catching Gamma-Ray Bursts Before They Fly

Authors: Jeremy S. Heyl

Abstract: The advent of large panoramic photometric surveys of the sky offers the possibly of exploring the association of gamma-ray bursts (GRBs) with supernovae. To date, a few gamma-ray bursts have been connected possibly with supernovae: GRB 980425 -- SN 1998bw, GRB 011121 -- SN 2001ke, GRB 970228 and GRB 980326. A combination of a large detection rate of GRBs and rapid coverage of a large portion of… ▽ More The advent of large panoramic photometric surveys of the sky offers the possibly of exploring the association of gamma-ray bursts (GRBs) with supernovae. To date, a few gamma-ray bursts have been connected possibly with supernovae: GRB 980425 -- SN 1998bw, GRB 011121 -- SN 2001ke, GRB 970228 and GRB 980326. A combination of a large detection rate of GRBs and rapid coverage of a large portion of the sky to faint magnitude limits offers the possibility of detecting a supernova preceding an associated GRB or at least placing limits on the rate of association between these two phenomena and the time delay between them. This would provide important constraints on theoretical models for gamma-ray bursts. △ Less

Submitted 3 April, 2003; v1 submitted 12 February, 2003; originally announced February 2003.

Comments: 8 pages, 1 figure, minor changes to reflect version accepted to Astrophys. J

Journal ref: Astrophys.J. 592 (2003) 401-403

The High-Energy Polarization-Limiting Radius of Neutron Star Magnetospheres I -- Slowly Rotating Neutron Stars

Authors: Jeremy S. Heyl, Nir J. Shaviv, Don Lloyd

Abstract: In the presence of strong magnetic fields, the vacuum becomes a birefringent medium. We show that this QED effect decouples the polarization modes of photons leaving the NS surface. Both the total intensity and the intensity in each of the two modes is preserved along a ray's path through the neutron-star magnetosphere. We analyze the consequences that this effect has on aligning the observed po… ▽ More In the presence of strong magnetic fields, the vacuum becomes a birefringent medium. We show that this QED effect decouples the polarization modes of photons leaving the NS surface. Both the total intensity and the intensity in each of the two modes is preserved along a ray's path through the neutron-star magnetosphere. We analyze the consequences that this effect has on aligning the observed polarization vectors across the image of the stellar surface to generate large net polarizations. Counter to previous predictions, we show that the thermal radiation of NSs should be highly polarized even in the optical. When detected, this polarization will be the first demonstration of vacuum birefringence. It could be used as a tool to prove the high magnetic field nature of AXPs and it could also be used to constrain physical NS parameters, such as $R/M$, to which the net polarization is sensitive. △ Less

Submitted 6 February, 2003; originally announced February 2003.

Comments: 23 pages, 9 figures

Journal ref: Mon.Not.Roy.Astron.Soc. 342 (2003) 134

Broadband Modeling of GRB 021004

Authors: Jeremy S. Heyl, Rosalba Perna

Abstract: We present a broadband modeling of the afterglow of GRB 021004. The optical transient of this burst has been detected very early and densely sampled in several bands. Its light curve shows significant deviations from a simple power law. We use the data from the X-ray to the $I-$band gathered in the first month of observations, and examine three models. Two involve variations in the energy of the… ▽ More We present a broadband modeling of the afterglow of GRB 021004. The optical transient of this burst has been detected very early and densely sampled in several bands. Its light curve shows significant deviations from a simple power law. We use the data from the X-ray to the $I-$band gathered in the first month of observations, and examine three models. Two involve variations in the energy of the shock. The first (energy injection) allows only increases to the shock energy, while the second (patchy shell) allows the energy to increase or decrease. In the final model (clumpy medium) the energy of the shock is constant while the density varies. While all three models reproduce well the optical bands, the variable density model can best account for the X-ray data, and the energy-injection model has the poorest fit. None of the models can account for the modest color variations observed during the first few days of the burst. △ Less

Submitted 5 February, 2003; v1 submitted 12 November, 2002; originally announced November 2002.

Comments: 13 pages, 5 figures, accepted to Astrophys. J. Letters, added analysis of X-ray lightcurve

Journal ref: Astrophys.J. 586 (2003) L13-L18