Arxiv:1707.08675V3 [Astro-Ph.EP] 18 Aug 2017
Total Page:16
File Type:pdf, Size:1020Kb
Draft version July 2, 2018 Typeset using LATEX twocolumn style in AASTeX61 PROSPECTS OF DYNAMICAL DETERMINATION OF GENERAL RELATIVITY PARAMETER β AND SOLAR QUADRUPOLE MOMENT J2 WITH ASTEROID RADAR ASTRONOMY Ashok K. Verma,1 Jean-Luc Margot,1, 2 and Adam H. Greenberg2 1Department of Earth, Planetary, and Space Sciences, University of California, Los Angeles, CA 90095, USA 2Department of Physics and Astronomy, University of California, Los Angeles, CA 90095, USA Submitted to ApJ ABSTRACT We evaluated the prospects of quantifying the parameterized post-Newtonian parameter β and solar quadrupole −1 moment J2 with observations of near-Earth asteroids with large orbital precession rates (9 to 27 arcsec century ). We considered existing optical and radar astrometry, as well as radar astrometry that can realistically be obtained with the Arecibo planetary radar in the next five years. Our sensitivity calculations relied on a traditional covariance analysis and Monte Carlo simulations. We found that independent estimates of β and J2 can be obtained with precisions of 6 × 10−4 and 3 × 10−8, respectively. Because we assumed rather conservative observational uncertainties, as is the usual practice when reporting radar astrometry, it is likely that the actual precision will be closer to 2 × 10−4 and 10−8, respectively. A purely dynamical determination of solar oblateness with asteroid radar astronomy may therefore rival the helioseismology determination. Keywords: astrometry | gravitation | minor planets, asteroids: general | relativistic processes | Sun: fundamental parameters | techniques: radar astronomy arXiv:1707.08675v3 [astro-ph.EP] 18 Aug 2017 Corresponding author: Ashok Kumar Verma [email protected] 2 Verma et al. 2017 1. INTRODUCTION (2:2±0:1)×10−7 (Will 2014). Dynamical estimates that The parameterized post-Newtonian (PPN) formal- do not rely on fits to helioseismology data yield similar −7 ism is a useful framework for testing metric theories values of J2 = 2:3 ± 0:25 × 10 (Fienga et al. 2015) −7 of gravity (Will 2014). It consists of 10 dimension- and J2 = 2:25 ± 0:09 × 10 (Park et al. 2017). High- less parameters that describe the general properties of precision dynamical estimates are important to validate the metric. In general relativity (GR), only 2 of the our understanding of the interior structure of the Sun. 10 parameters are non-zero. They are known as the Our simulations of the determination of β and J2 Eddington−Robertson−Schiff parameters γ and β. γ using a variety of asteroid orbits suggest that indepen- represents the amount of curvature produced by a unit dent values of β and J2 can be obtained with satisfac- mass, and β represents the amount of nonlinearity in tory precision: with the traditionally conservative as- the superposition law for gravity. signment of radar uncertainties, β can be constrained at −4 Several techniques have been used to place obser- the 6 × 10 level and J2 can be constrained at the −8 vational bounds on these parameters (Will 2014), in- 3 × 10 level. With uncertainties that more closely cluding observations of the bending and delay of light reflect measurement errors, this precision may be im- by spacecraft tracking (Bertotti et al. 2003, e.g.,) or proved by a factor of ∼3. (Section4). Very Long Baseline Interferometry (e.g., Lambert & The outline of this paper is as follows. In Section2, Le Poncin-Lafitte 2009), and fitting of ephemerides to we describe our choice of target asteroids. In Section3, observations of planetary positions (e.g., Folkner 2009; we discuss the estimation of asteroid orbits with opti- Fienga et al. 2011; Verma et al. 2014; Fienga et al. 2015). cal and radar measurements. Our dynamical model and In GR, γ and β are equal to one. Doppler track- data reduction procedures are described in Section 3.1 ing of the Cassini spacecraft has shown that γ does and 3.2, respectively. Orbit determination results are not differ from one by more than 2 × 10−5 (Bertotti presented in Section 3.3. Simulations of the determina- et al. 2003). Ephemeris-based studies prior to 2009 indi- tion of β and J2 are described in Section4. cated that β − 1 does not differ from zero by more than 2. TARGET ASTEROIDS 10−4 (Folkner 2009; Pitjeva & Pitjev 2014). More re- cently, the availability of precise ranging data from the The per-orbit secular advance in the angular position MESSENGER Mercury orbiter (Solomon et al. 2001) of the perihelion is given by (Misner et al. 1973) enabled improved estimates (Verma et al. 2014; Fienga 6πGM (2 − β + 2γ) 6π (1 − 3=2 sin2 i) δ! = + R2 J ; et al. 2015; Park et al. 2017). Here, we evaluate the a(1 − e2)c2 3 2 a2(1 − e2)2 2 prospect of asteroid orbit precession measurements to (1) place more stringent bounds on β. We consider Earth- where ! is the argument of perihelion, GM is the Sun's based radar observations of near-Earth asteroids with gravitational parameter, R is the radius of the Sun, c −1 perihelion shifts larger than 10 arcsec century . is the speed of light, and a, e, and i are the semi-major Orbital precession can also be caused by the nonuni- axis, eccentricity, and orbital inclination (with respect formity of the gravity field that results from the oblate to the solar equator) of a planetary body, respectively. shape of the Sun. The solar oblateness is character- Because both GR and solar oblateness affect perihelion ized by the solar quadrupole moment, J2 (e.g., Kaula precession, estimates of β and J2 are highly correlated 2000). Simultaneous estimation of β and J2 requires and it is desirable to track a variety of solar system that the precessional effects due to GR and to the Sun's bodies with a range of a, e, i values to disentangle the oblateness be disentangled. Fortunately, GR is a purely two effects. central effect, whereas the oblateness-induced precession Our selection of target asteroids follows the method of has an inclination dependence. The two effects also have Margot(2003). We select asteroids with both large per- a different distance dependence (Misner et al. 1973). As ihelion shift values and favorable observing conditions a result, observations of a small sample of near-Earth with radar (Table1 and Figure1). This sample of as- asteroids with a variety of semi-major axes and inclina- teroid orbits includes a wide range of semi-major axes, tions (Table1) can in principle be used to estimate β eccentricities, and inclinations, which are advantageous and J2 (Margot 2003; Margot & Giorgini 2009). when simultaneously solving for β and J2 . The pre- Current estimates of the solar quadrupole moment dicted rates of perihelion advance, δ!_ , shown in Figure are typically derived on the basis of interior models 1 and Table1 were computed assuming γ = β = 1 and of the Sun constrained by helioseismology data (e.g., −7 J2 = 2:2 × 10 . Mecheri et al. 2004; Antia et al. 2008). The current best value from the helioseismology literature is J2 = 3. METHODS Prospects of determination of GR parameter β with asteroid radar astronomy 3 Table 1. Selected asteroids and orbital elements: Semima- the value of the semi-major axis, but Yarkovsky drift jor Axis (a), Eccentricity (e), and Inclination with Respect does. This nongravitational effect has been shown to to the Ecliptic (iec) and Sun's equator (ieq). affect the semi-major axes of small bodies due to the 00 −1 anisotropic re-emission of absorbed sunlight (e.g., Bot- Target a (au) e iec (deg) ieq (deg) δ!_ ( cy ) tke et al. 2006). The change in semi-major axis with 1566 Icarus 1.078 0.827 22.9 15.8 10.1 time due to Yarkovsky orbital drift, hda=dti, was esti- 1998 TU3 0.787 0.484 5.41 3.41 9.11 mated for all target asteroids with the method of Green- 1999 KW4 0.642 0.688 38.9 46.0 22.1 berg et al.(2017). The values ranged in amplitude be- 1999 MN 0.674 0.665 2.02 5.25 18.5 tween 4 and 50 au/My, which is plausible for kilometer- 2000 BD19 0.876 0.895 25.7 28.0 26.9 sized bodies. Only one target (1566 Icarus) is common 2000 EE14 0.662 0.533 26.5 26.1 15.0 between our target list and the 42 Yarkovsky detections 2001 YE4 0.677 0.541 4.82 11.0 14.4 of Nugent et al.(2012), and only one target (1999 MN) 2004 KH17 0.712 0.499 22.1 14.9 12.0 is common between our target list and the 21 Yarkovsky 2006 CJ 0.676 0.755 10.3 16.1 23.7 detections of Farnocchia et al.(2013). In both cases, our Yarkovsky drift estimates are consistent with and better constrained than prior work. Note. The predicted rate of perihelion advance in arcsec To initialize the integration process, we used a priori century−1 (00 cy−1), δ!_ , was computed using Equation (1). state vectors extracted from the Minor Planet Center (MPC) database (Minor Planet Center 2017). We first determined nominal trajectories for asteroids 3.2. Existing optical and radar astrometry in our sample with astrometric (i.e., positional) data, both optical and radar (Table2). The process involved We used both optical and radar astrometry to de- three steps: (1) numerical integration of each asteroid's termine the nominal trajectory of each asteroid. Op- orbit and calculation of partial derivatives of the equa- tical measurements provide positional information on tions of motion with respect to the solve-for parameters the plane of the sky. They are typically expressed as (i.e., the six components of the state vectors), (2) eval- right ascension (R.A.) and declination (decl.) in the uation of simulated optical and radar observables and equatorial frame of epoch J2000.0.