Researcher profile

Jeremy S. Heyl

Jeremy S. Heyl contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
28works
0followers
12topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

28 published item(s)

preprint2016arXiv

Eötvös Experiments with Supermassive Black Holes

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.

preprint2016arXiv

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

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.

preprint2015arXiv

Casimir Interactions between Magnetic Flux Tubes in a Dense Lattice

We use the worldline numerics technique to study a cylindrically symmetric model of magnetic flux tubes in a dense lattice and the non-local Casimir forces acting between regions of magnetic flux. Within a superconductor the magnetic field is constrained within magnetic flux tubes and if the background magnetic field is on the order the quantum critical field strength, $B_k = \frac{m^2}{e} = 4.4 \times 10^{13}$ Gauss, the magnetic field is likely to vary rapidly on the scales where \acs{QED} effects are important. In this paper, we construct a cylindrically symmetric toy model of a flux tube lattice in which the non-local influence of \acs{QED} on neighbouring flux tubes is taken into account. We compute the effective action densities using the worldline numerics technique. The numerics predict a greater effective energy density in the region of the flux tube, but a smaller energy density in the regions between the flux tubes compared to a locally-constant-field approximation. We also compute the interaction energy between a flux tube and its neighbours as the lattice spacing is reduced from infinity. Because our flux tubes exhibit compact support, this energy is entirely non-local and predicted to be zero in local approximations such as the derivative expansion. This Casimir-Polder energy can take positive or negative values depending on the distance between the flux tubes, and it may cause the flux tubes in neutron stars to form bunches. In addition to the above results we also discuss two important subtleties of determining the statistical uncertainties within the worldline numerics technique and recommend a form of jackknife analysis.

preprint2014arXiv

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

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.

preprint2014arXiv

Energy levels of light atoms in strong magnetic fields

In this review article we provide an overview of the field of atomic structure of light atoms in strong magnetic fields. There is a very rich history of this field which dates back to the very birth of quantum mechanics. At various points in the past significant discoveries in science and technology have repeatedly served to rejuvenate interest in atomic structure in strong fields, broadly speaking, resulting in three eras in the development of this field; the historical, the classical and the modern eras. The motivations for studying atomic structure have also changed significantly as time progressed. The review presents a chronological summary of the major advances that occurred during these eras and discusses new insights and impetus gained. The review is concluded with a description of the latest findings and the future prospects for one of the most remarkably cutting-edge fields of research in science today.

preprint2014arXiv

Parallel Worldline Numerics: Implementation and Error Analysis

We give an overview of the worldline numerics technique, and discuss the parallel CUDA implementation of a worldline numerics algorithm. In the worldline numerics technique, we wish to generate an ensemble of representative closed-loop particle trajectories, and use these to compute an approximate average value for Wilson loops. We show how this can be done with a specific emphasis on cylindrically symmetric magnetic fields. The fine-grained, massive parallelism provided by the GPU architecture results in considerable speedup in computing Wilson loop averages. Furthermore, we give a brief overview of uncertainty analysis in the worldline numerics method. There are uncertainties from discretizing each loop, and from using a statistical ensemble of representative loops. The former can be minimized so that the latter dominates. However, determining the statistical uncertainties is complicated by two subtleties. Firstly, the distributions generated by the worldline ensembles are highly non-Gaussian, and so the standard error in the mean is not a good measure of the statistical uncertainty. Secondly, because the same ensemble of worldlines is used to compute the Wilson loops at different values of $T$ and $x_\mathrm{ cm}$, the uncertainties associated with each computed value of the integrand are strongly correlated. We recommend a form of jackknife analysis which deals with both of these problems.

preprint2013arXiv

Diffractive Microlensing: A New Probe of the Local Universe

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.

preprint2013arXiv

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

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.

preprint2013arXiv

Magnetic Reconnection Instabilities in Soft-Gamma Repeaters

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.

preprint2013arXiv

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

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&#39;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.

preprint2012arXiv

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

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.

preprint2012arXiv

Hydrogen and Helium atoms in strong magnetic fields

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.

preprint2012arXiv

Is Mira a magneto-dusty rotator?

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&#39;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.

preprint2012arXiv

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

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.

preprint2012arXiv

The Magnetised Bellows of Betelgeuse

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.

preprint2011arXiv

Constraining the photon-axion coupling constant with magnetic white dwarfs

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.

preprint2010arXiv

Cosmic Rays from Pulsars and Magnetars

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.

preprint2010arXiv

Diffractive Microlensing II: Substellar Disk and Halo Objects

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.

preprint2010arXiv

Diffractive Microlensing III: Astrometric Signatures

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).

preprint2010arXiv

Nonlinear Electromagnetic Waves in Magnetosphere of a Magnetar

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.

preprint2010arXiv

On the trigger mechanisms for SGR giant flares

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.

preprint2010arXiv

Pseudospectral methods for atoms in strong magnetic fields

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).

preprint2009arXiv

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

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.

preprint2004arXiv

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

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,&#39;&#39; 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).

preprint2002arXiv

QED and the High Polarization of the Thermal Radiation from Neutron Stars

The thermal emission of strongly magnetized neutron-star atmospheres is thought to be highly polarized. However, because of the different orientations of the magnetic field over the surface of the neutron star (NS), it is commonly assumed that the net observed polarization will be significantly reduced as the polarization from different regions will cancel each other. We show that the birefringence of the magnetized QED vacuum decouples the polarization modes in the magnetosphere; therefore, the direction of the polarization follows the direction of the magnetic field up to a large distance from the stellar surface. At this distance, the rays that leave the surface and are destined for our detectors pass through only a small solid angle; consequently, the polarization direction of the emission originating in different regions will tend to align together. The net observed polarization of the thermal radiation of NSs should therefore be very large. Measurement of this polarization will be the first direct evidence of the birefringence of the magnetized vacuum due to QED and a direct probe of behavior of the vacuum at magnetic fields of order of and above the critical QED field of 4.4 x 10 13 G. The large observable polarization will also help us learn more about the atmospheric properties of NSs.

preprint2000arXiv

Electron-Positron Jets from a Critically Magnetized Black Hole

The curved spacetime surrounding a rotating black hole dramatically alters the structure of nearby electromagnetic fields. The Wald field which is an asymptotically uniform magnetic field aligned with the angular momentum of the hole provides a convenient starting point to analyze the effects of radiative corrections on electrodynamics in curved spacetime. Since the curvature of the spacetime is small on the scale of the electron&#39;s Compton wavelength, the tools of quantum field theory in flat spacetime are reliable and show that a rotating black hole immersed in a magnetic field approaching the quantum critical value of $B_k=m^2 c^3/(e\hbar) \approx 4.4 \times 10^{13}$~G $\approx 1.3\times10^{-11}$ cm$^{-1}$ is unstable. Specifically, a maximally rotating three-solar-mass black hole immersed in a magnetic field of $2.3 \times 10^{12}$~G would be a copious producer of electron-positron pairs with a luminosity of $3 \times 10^{52}$ erg s$^{-1}$.

preprint1994arXiv

Galaxy Formation in a Variety of Hierarchical Models

We predict the observable properties of the galaxy population in popular hierarchical models of galaxy formation. We employ a detailed semianalytic procedure which incorporates the formation and merging of dark matter halos, the shock heating and radiative cooling of gas, self-regulated star formation, the merging of galaxies within dark matter halos, and the spectral evolution of the stellar populations. We contrast the standard CDM cosmogony with variants of the CDM model having either a low value of H_0, or a low value of Omega with or without a cosmological constant. In addition, we compare galaxy formation in these CDM universes with a CHDM model. We find that although the models have some success in remedying the shortcomings of the standard CDM cosmogony, none of these new models produce broad agreement with the whole range of observations. Although the low-Omega and Omega+Lambda=1 CDM models reduce the discrepancy between the predicted and observed Tully-Fisher relations (the main weakness of galaxy formation in standard CDM), these models predict an inverted colour-magnitude relation and do not produce an exponential cut-off at the bright end of the galaxy luminosity function. All of our models predict recent star formation and exhibit galaxy colours bluer than observed, but this problem is far more severe in the CHDM model which produces colours about two magnitudes too blue in B-K. Unlike in the variants of the CDM model in the CHDM case this result is not dependent on our model of stellar feedback, but is instead directly caused by the late epoch of structure formation in this model.