skip to main content


Title: Building an Acceleration Ladder with Tidal Streams and Pulsar Timing
Abstract

We analyse stellar streams in action-angle coordinates combined with recent local direct acceleration measurements to provide joint constraints on the potential of our galaxy. Our stream analysis uses the Kullback–Leibler divergence with a likelihood analysis based on the two-point correlation function. We provide joint constraints from pulsar accelerations and stellar streams for local and global parameters that describe the potential of the Milky Way (MW). Our goal is to build an “acceleration ladder,” where direct acceleration measurements that are currently limited in dynamic range are combined with indirect techniques that can access a much larger volume of the MW. To constrain the MW potential with stellar streams, we consider the Palomar 5, Orphan, Nyx, Helmi, and GD1 streams. Of the potential models that we have considered here, the preferred potential for the streams is a two-component Staeckel potential. We also compare the vertical accelerations from stellar streams and pulsar timing, defining the functionf(z)=α1pulsarzΦz, where Φ is the MW potential determined from stellar streams andα1 pulsarzis the vertical acceleration determined from pulsar timing observations. Our analysis indicates that the Oort limit determined from streams is consistently (regardless of the choice of potential) lower than that determined from pulsar timing observations. The calibration we have derived here may be used to correct the estimate of the acceleration from stellar streams.

 
more » « less
Award ID(s):
2007232
NSF-PAR ID:
10401470
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal Letters
Volume:
945
Issue:
2
ISSN:
2041-8205
Format(s):
Medium: X Size: Article No. L32
Size(s):
["Article No. L32"]
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Supermassive black hole binary systems (SMBHBs) should be the most powerful sources of gravitational waves (GWs) in the universe. Once pulsar timing arrays (PTAs) detect the stochastic GW background from their cosmic merger history, searching for individually resolvable binaries will take on new importance. Since these individual SMBHBs are expected to be rare, here we explore how strong gravitational lensing can act as a tool for increasing their detection prospects by magnifying fainter sources and bringing them into view. Unlike for electromagnetic waves, when the geometric optics limit is nearly always valid, for GWs the wave-diffraction-interference effects can become important when the wavelength of the GWs is larger than the Schwarzchild radius of the lens, i.e.,Mlens108fmHz1M. For the GW frequency range explored in this work, the geometric optics limit holds. We investigate GW signals from SMBHBs that might be detectable with current and future PTAs under the assumption that quasars serve as bright beacons that signal a recent merger. Using the black hole mass function derived from quasars and a physically motivated magnification distribution, we expect to detect a few strongly lensed binary systems out toz≈ 2. Additionally, for a range of fixed magnifications 2 ≤μ≤ 100, strong lensing adds up to ∼30 more detectable binaries for PTAs. Finally, we investigate the possibility of observing both time-delayed electromagnetic signals and GW signals from these strongly lensed binary systems—that will provide us with unprecedented multi-messenger insights into their orbital evolution.

     
    more » « less
  2. Abstract

    The merger of two galaxies, each hosting a supermassive black hole (SMBH) of mass 106Mor more, could yield a bound SMBH binary. For the early-type galaxy NGC 4472, we study how astrometry with a next-generation Very Large Array could be used to monitor the reflex motion of the primary SMBH of massMpri, as it is tugged on by the secondary SMBH of massMsec. Casting the orbit of the putative SMBH binary in terms of its periodP, semimajor axisabin, and mass ratioq=Msec/Mpri1, we find the following: (1) Orbits with fiducial periods ofP= 4 yr and 40 yr could be spatially resolved and monitored. (2) For a 95% accuracy of 2μas per monitoring epoch, subparsec values ofabincould be accessed over a range of mass ratios notionally encompassing majorq>14and minorq<14galaxy mergers. (3) If no reflex motion is detected forMpriafter 1 (10) yr of monitoring, an SMBH binary with periodP= 4 (40) yr and mass ratioq> 0.01 (0.003) could be excluded. This would suggest no present-day evidence for a past major merger like that recently simulated, where scouring by aq∼ 1 SMBH binary formed a stellar core with kinematic traits like those of NGC 4472. (4) Astrometric monitoring could independently check the upper limits onqfrom searches for continuous gravitational waves from NGC 4472.

     
    more » « less
  3. Abstract

    The genericity of Arnold diffusion in the analytic category is an open problem. In this paper, we study this problem in the followinga prioriunstable Hamiltonian system with a time-periodic perturbationHε(p,q,I,φ,t)=h(I)+i=1n±12pi2+Vi(qi)+εH1(p,q,I,φ,t),where(p,q)Rn×Tn,(I,φ)Rd×Tdwithn,d⩾ 1,Viare Morse potentials, andɛis a small non-zero parameter. The unperturbed Hamiltonian is not necessarily convex, and the induced inner dynamics does not need to satisfy a twist condition. Using geometric methods we prove that Arnold diffusion occurs for generic analytic perturbationsH1. Indeed, the set of admissibleH1isCωdense andC3open (a fortiori,Cωopen). Our perturbative technique for the genericity is valid in theCktopology for allk∈ [3, ∞) ∪ {∞,ω}.

     
    more » « less
  4. Abstract

    The stellar initial mass function (IMF) is a fundamental property in the measurement of stellar masses and galaxy star formation histories. In this work, we focus on the most massive galaxies in the nearby universelog(M/M)>11.2. We obtain high-quality Magellan/LDSS-3 long-slit spectroscopy with a wide wavelength coverage of 0.4–1.01μm for 41 early-type galaxies (ETGs) in the MASSIVE survey and derive high signal-to-noise spectra within an aperture ofRe/8. Using detailed stellar synthesis models, we constrain the elemental abundances and stellar IMF of each galaxy through full spectral modeling. All the ETGs in our sample have an IMF that is steeper than a Milky Way (Kroupa) IMF. The best-fit IMF mismatch parameter,αIMF= (M/L)/(M/L)MW, ranges from 1.1 to 3.1, with an average of 〈αIMF〉 = 1.84, suggesting that on average, the IMF is more bottom heavy than Salpeter. Comparing the estimated stellar masses with the dynamical masses, we find that most galaxies have stellar masses that are smaller than their dynamical masses within the 1σuncertainty. We complement our sample with lower-mass galaxies from the literature and confirm thatlog(αIMF)is positively correlated withlog(σ),log(M), andlog(Mdyn). From the combined sample, we show that the IMF in the centers of more massive ETGs is more bottom heavy. In addition, we find thatlog(αIMF)is positively correlated with both [Mg/Fe] and the estimated total metallicity [Z/H]. We find suggestive evidence that the effective stellar surface density ΣKroupamight be responsible for the variation ofαIMF. We conclude thatσ, [Mg/Fe], and [Z/H] are the primary drivers of the global stellar IMF variation.

     
    more » « less
  5. Abstract

    A steady-state, semi-analytical model of energetic particle acceleration in radio-jet shear flows due to cosmic-ray viscosity obtained by Webb et al. is generalized to take into account more general cosmic-ray boundary spectra. This involves solving a mixed Dirichlet–Von Neumann boundary value problem at the edge of the jet. The energetic particle distribution functionf0(r,p) at cylindrical radiusrfrom the jet axis (assumed to lie along thez-axis) is given by convolving the particle momentum spectrumf0(,p)with the Green’s functionG(r,p;p), which describes the monoenergetic spectrum solution in whichf0δ(pp)asr→ ∞ . Previous work by Webb et al. studied only the Green’s function solution forG(r,p;p). In this paper, we explore for the first time, solutions for more general and realistic forms forf0(,p). The flow velocityu=u(r)ezis along the axis of the jet (thez-axis).uis independent ofz, andu(r) is a monotonic decreasing function ofr. The scattering timeτ(r,p)=τ0(p/p0)αin the shear flow region 0 <r<r2, andτ(r,p)=τ0(p/p0)α(r/r2)s, wheres> 0 in the regionr>r2is outside the jet. Other original aspects of the analysis are (i) the use of cosmic ray flow lines in (r,p) space to clarify the particle spatial transport and momentum changes and (ii) the determination of the probability distributionψp(r,p;p)that particles observed at (r,p) originated fromr→ ∞ with momentump. The acceleration of ultrahigh-energy cosmic rays in active galactic nuclei jet sources is discussed. Leaky box models for electron acceleration are described.

     
    more » « less