system, and Ðnally, in 7, we perform numerical tests which show how it can be used in practice.

Size: px
Start display at page:

Download "system, and Ðnally, in 7, we perform numerical tests which show how it can be used in practice."

Transcription

1 THE ASTROPHYSICAL JOURNAL 544:921È December 1 ( The American Astronomical Society. All rights reserved. Printed in U.S.A. IMPROVED TIMING FORMULA FOR THE PSR B1257]12 PLANETARY SYSTEM MACIEJ KONACKI1 AND ANDRZEJ J. MACIEJEWSKI2 Torun Centre for Astronomy Nicolaus Copernicus University Gagarina Torun Poland AND ALEX WOLSZCZAN3 Department of Astronomy and Astrophysics Pennsylvania State University University Park PA 16802; and Torun Centre for Astronomy Nicolaus Copernicus University ul. Gagarina Torun Poland Received 2000 February 8; accepted 2000 June 27 ABSTRACT We present a new analysis of the dynamics of the planetary system around the pulsar B1257]12. A semianalytical theory of perturbation between terrestrial-mass planets B and C is developed and applied to improve the multiorbit timing formula for this obect. We use numerical simulations of the pulse arrival times for PSR B1257]12 to demonstrate that our new timing model can serve as a tool to determine the masses of the two planets. Subect headings: celestial mechanics stellar dynamics È planetary systems È pulsars: individual (PSR B1257]12) 1. INTRODUCTION The Ðrst extrasolar planetary system has been discovered around the millisecond pulsar B1257]12 (Wolszczan & Frail 1992). The system consists of three planets named A B and C with planets B and C having the orbits close to a 3: 2 commensurability. This circumstance allows us to analyze the dynamics of the system beyond the classical Keplerian approximation. Namely in such a conðguration the gravitational interactions of planets B and C give rise to observable time variations of the B and C orbital elements. It was thoroughly discussed by Rasio et al. (1992) and Malhotra (1993a) (see also Malhotra et al. 1992; Malhotra 1993b; Rasio et al. 1993; Peale 1993). Subsequently these studies were used to conðrm the existence of the PSR B1257]12 planetary system through detecting in the timing observations of the pulsar the presence of those non- Keplerian variations (Wolszczan 1994). However the application of non-keplerian dynamics goes further than the conðrmation of the discovery. It can be used to derive some interesting information about the system which is not otherwise accessible. The aim of this paper is to apply the theory of perturbed planetary motion and derive an improved model for the timing observations of PSR B1257]12. Such a model as we show on simulations should lead to a determination of the masses of planets B and C as well as the inclinations of their orbits. To this end in 2 we analyze the equations of motion in the barycentric and Jacobi coordinate systems which we use in the paper. In 3 we show how these two slightly di erent representations are related to the commonly used timing model for pulsars with companions. In 4 we demonstrate how to express such a timing formula in terms of the osculating orbital elements. In 5 we show how to obtain the osculating elements of B and C. In 6 we present an improved timing model describing the motion of this 1 kmc=astri.uni.torun.pl. 2 macieka=astri.uni.torun.pl. 3 alex=astro.psu.edu. system and Ðnally in 7 we perform numerical tests which show how it can be used in practice. 2. EQUATIONS OF MOTION Let us consider a system consisting of a neutron star P 0 with mass m and N planets P with masses m. In an arbi- 0 i i trary inertial reference frame equations of motion of this system have the form m i R i \[G ; N m m i R3 (R i [ R ) (1) /0 i Ei R i \ p R i [ R p i \ 0... N and G is the gravitational constant. The Hamiltonian function for equation (1) has the form H \ 1 2 ; N 1 Gm m P2[ ; i (2) m i R i/0 i 0yi:yN i and equation (1) can be written as HamiltonÏs equations d dt R i \ LH d LP dt P i \[LH i \ 0... N. (3) LR i i Analytical perturbation theory for a planetary system is usually formulated in the so-called Jacobi coordinates r i \ 0... N which are deðned in the following way: i r \ R [ 1 k~1 ; mi R for k \ 1... N (4) k k k i k~1 i/0 r \ 1 N k ; mi R k \ ; mi. (5) 0 k i k N i/0 i/0 In other words r is the radius vector from the center of mass of bodies P i... P to body P and r is the center 0 i~1 i 0 of mass of the system. The above formulae deðne a one-to- 921

2 922 KONACKI MACIEJEWSKI & WOLSZCZAN Vol. 544 one relationship between Cartesian inertial coordinates and Jacobi coordinates. The inverse relationship has the form R \ r ] k N m k~1 r [ ; i r k \ 0... N (6) k 0 k k k i k i/k`1 i we assume k \ 0. In terms of canonical Jacobi coordinates the Hamiltonian ~1 (eq. [2]) reads H \ 1 N k N m m p2] ; i p2[g ; 0 i k 0 k m i r N i/1 i~1 i i/1 i ]H3 1 (r 1... r N ) (7) H3 is its perturbative part. If we 1 \ H3 1 (r 1... r N ) assume that masses of planets are small and all are of the same order v then H3 is of order v2 i.e. H3 1 1 \ H 1 ] O(v3) A 1 B H \[G ; m m [ r i Æ r. (8) 1 i r r 1yi:yN i i From the form of the Hamiltonian (eq. [7]) it follows that p is a Ðrst integral and that the equations for variables 0 Mr... r p... p N do not depend on its particular 1 N 1 N value. Thus we can assume that p \ 0 implying that r is 0 0 constant and can be set to r \ 0. This is equivalent to the 0 assumption that our inertial frame is a certain barycentric reference frame. The dynamics of the system are governed by HamiltonÏs equations with the Hamiltonian of the form H \ H ] H ] ÉÉÉ (9) 0 1 N k N m m H \ ; i p2[g ; 0 i (10) 0 k m i r i/1 i~1 i i/1 i is the unperturbed part of the Hamiltonian describing a system of N independent planets. It follows from the form of equation (10) that each planet moves in a Keplerian orbit in the same way as a body with the mass m k /k around a i i~1 i Ðxed gravitational center m m. Each of these Keplerian 0 i motions can be parameterized by Keplerian elements MT a e i u )NÈthe time of pericenter semimaor axis p eccentricity inclination argument of pericenter and longitude of ascending node respectively. 3. TIMING MODEL AND COORDINATE SYSTEMS Timing observations of pulsars represent measurements of the times of arrival (TOAs) of pulsars pulses. An extraordinary precision of timing measurements allows a detection of very low level e ects in timing residuals (see for review Lyne & Graham-Smith 1998). In the case of a binary pulsar the observed TOAs exhibit periodic variations resulting from the motion of the pulsar around the center of mass. Such variations are modeled with the formula *t \ x[(cos E [ e) sin u ] J1 [ e2 sin E cos u] x \ a sin i/c E [ e sin E \ n(t [ T p ) n \ 2n P (11) and *t is the additional delay/advance in TOAs E is the eccentric anomaly a e u sin i P T are Keplerian elements of the orbit and c is the speed of p light. When a pulsar has N planets the TOA variations become N *t \ ; x C (cos E [ e ) sin u /1 D ]J1 [ e2 (sin E cos u ). (12) Note that in practice the parameters of the model which we determine by means of the least-squares Ðt to the data are x e u P T i.e. the proection of semimaor axis of the pulsarïs orbit p eccentricity argument of pericenter orbital period and time of pericenter. In order to precisely understand and interpret these parameters we describe the pulsarïs motion in the barycentric reference frame with the z-axis of the system directed toward the barycenter of the solar system and x-y plane of the reference frame in the plane of the sky. This way we have *t \[ 1 c R 0 Æ ZŒ R 0 \[ 1 N ; mi R (13) m i 0 i/1 ZŒ is the unit vector along the z-axis of the pulsar system barycentric frame and R are barycentric positions of planets. Assuming that m a sin i \ x c G (1 ] m / )2 \ n 2 a 3 (14) with planetsï orbital parameters a n e u T we can p most naturally interpret the motion of the pulsar as a superposition of the elliptic motions of its planets around the barycenter of the system. However as it was mentioned in the previous section the analytical perturbation theory is usually formulated in Jacobi coordinates in which the TOA variations become *t \[ 1 c R 0 Æ ZŒ R 0 \[; N i r i \ m (15) k /1 r are positions of planets. Furthermore we have the following relations: Gm k i a sin i \ x c 0 \ G \ n2 a3. (16) k 1 [ i ~1 Thus from the timing formula for TOA variations of a pulsar with planets it is possible to obtain two somewhat di erent descriptions of the pulsar motion. Although they both represent a sum of a certain number of elliptic motions the interpretation of some of their parameters is slightly di erent. Throughout the rest of this paper we will use Jacobi coordinates as they are more convenient in the formulation of the theory of perturbed motion. 4. OSCULATING ORBITAL ELEMENTS The non-keplerian motion of the PSR B1257]12 system can be described by means of the osculating ellipses (i.e. by means of ellipses whose parameters change with time). The time evolution of orbital elements can be determined by solving the classical LagrangeÏs perturbation

3 No PSR B1257]12 PLANETARY SYSTEM 923 equations (see for example Brouwer & Clemence 1961) with the perturbation Hamiltonian given by equation (8). Such approach leads to the solution in the form x \ x0]*x(t [ t 0 ) (17) x stands for a speciðc orbital element x0 for its initial value and *x for its time-dependent part of small magnitude *x(t [ t )/x0 > 1. In the case of the PSR B1257]12 planetary systems 0 the most signiðcant part of the perturbations comes from planets B and C. Therefore we can assume that the orbital elements of planet A are approximately constant while the elements of planets B and C change with time. Thus using the formulae for *t from the previous section (eq. [12]) and assuming the time evolution of orbital elements in the form of equation (17) the additional TOA variations dt \ MB CN due to the inter- actions between planets B and C can be expressed as follows: cdt i sin i0 \[*h C3 0)D 2 a 0]1 2 a 0 cos (2 ] 1 2 *k a 0 sin (2 0) ] *a sin 0]* a 0 cos 0 (18) h \ e sin u k \ e cos u \ n (t [ T p ) ] u (19) and equation (18) is given to Ðrst order in *a * *h *k and the lowest order in e (or h and k since the eccentric- ities of planets B and C are very small). The above equations can be obtained from the following well-known expansions: cos E \[ e 2 ] 2 ; 1 k J (ke) cos (km) k~1 k Z0 sin E \ 2 e ; 1 k J (ke) sin (km) M \ [ u (20) k~1 k Z0 J (x) is a Bessel function of the Ðrst kind of order k and argument k x Z denotes the set of all positive and negative integers excluding 0 zero and an approximate relation for Bessel functions of small arguments the case of PSR B1257]12 periodic variations can be related to conunctions ( close encounters ÏÏ) of planets B and C with the frequency n \ n [ n and to the 3:2 near commensurability with the ce frequency B C n \ 3n [ 2n. The dynamics of this system was studied by r Rasio C et al. B (1992) and Malhotra (1993a) (see also Malhotra 1993b; Rasio et al. 1993; Peale 1993). Malhotra (1993a) solved the LagrangeÏs perturbation equations for the orbital elements assuming a nonresonant system with coplanar orbits. This Ðrst-order perturbation theory is valid for (sin i)~1 Z 10 (orbital inclinations larger than about 6 ) which corresponds to mass ratios m /m Z 6 ] 10~5 (and the nonresonant case). When planets B and 0 C are locked in the exact resonance it is necessary to develop a di erent theory of motion (Malhotra et al. 1992; Rasio et al. 1993). It turns out that the Ðrst-order perturbation theory for this system developed by Malhotra (1993a) is not accurate enough for the purpose of determination of the planetsï masses. Therefore in this paper we present a new approach to this problem which precisely addresses the issue. Namely Ðrst we assume that the orbits of B and C are not coplanar however their relative inclination I is small. The geometry of the system can be represented as in Figure 1 and the perturbative part of Hamiltonian expanded to the Ðrst order in I and the product of the mass of planets B and C m m has the following form: 1 2 H 1 \[ Gm 1 m 2 r 2 GC 1 [ 2 r 1 r2 cos t ] Ar 1 r2 B2D [ r 1 r2 cos t H Π 1= ω1 τ1 Π = ω τ z Π 1 Π2 n 12 (23) t \ f ] u [ f [ u [ q q \ q [ q r \ a (1 [ e 2) \ 12 1 ] e cos f and f is the true anomaly. Note that in the above expan- sion H does not depend on I. The Ðrst term in which I 1 J k (x) B xk 2kk! x > 1. (21) Now our next step is to Ðnd the explicit form of the functions O Ω 2 I τ 1 τ 2 i 2 y a (t) \ a 0]*a (t) (t) \ 0]* (t) h (t) \ h 0]*h (t) k (t) \ k 0]*k (t). (22) x Ω 1 n 1 i 1 n2 5. SEMIANALYTICAL PERTURBATION THEORY Mutual gravitational interactions between planets cause periodic and secular changes of their orbital elements. In FIG. 1.ÈGeometry of the system. The angles u u are the arguments of pericenter ) ) longitudes of ascending node 1 and 2 i i inclinations of the orbit of the 1 planets 2 B and C; q and q are the angles 1 2 n On and n On respectively. The angles I 1 q q can 2 be found by solving 1 12 the spherical 2 12 triangle n n n

4 924 KONACKI MACIEJEWSKI & WOLSZCZAN Vol. 544 appears is of the second order; precisely it is sin2 (I/2). Thus as long as sin2 (I/2) is negligible the above expansion is valid. Next from the classical form of LagrangeÏs perturbation equations (see for example Brouwer & Clemence 1961) we can obtain the following set of the Ðrst-order di erential equations for the elements a e u : a5 \ [2 k n a LH 1 Lp e5 \ 1 k n a 2 e C [ A 1 [ e 2 B LH 1 Lp ] J1 [ e 2 LH 1 Lu D u5 \ [J1 [ e 2 k n a 2 e LH 1 Le 5 \ n ] 2 LH e J1 [ e2 1 [ k n a La k n a2(1 ] J1 [ e2) (24) n is the mean motion and p is related to time of pericenter through p \[n T so the Kepler equation p reads E [ e sin E \ n t ] p \ [ u (25) and the true anomaly f is related to the eccentric anomaly E through tan f 2 \S1 ] e 1 [ e tan E 2. (26) And it should be noted that the elements ) i remain approximately constant in the case of small relative inclinations. In principle such a set of equations could be solved for a e u. In practice however it is extremely compli- cated. On the other hand in fact we do not need analytical formulae such as those presented in Malhotra (1993a). Thus the most suitable approach is to solve this problem numerically. In order to do so we ust need the explicit form of the right-hand-side functions of equation (24) which can be easily obtained using the perturbing Hamiltonian H. Subsequently the equations can be solved numerically 1 for a e u. As we show in 7 this approach gives very accurate results much more accurate than any reasonable analytical treatment. From the form of equation (24) it follows that the change of orbital parameters *x (strictly speaking here *x stands for * *h *k and *a/a ) is proportional to m /m and m /m. Precisely corrections 0 *x are proportional 1 to 0 m /m 2 and 0 *x are proportional to 1 m /m. Let us Ðnally note 2 0 that in this 2 model there is one additional 1 0 parameter q which in general is not known a priori. However it can be determined through a leastsquares analysis of the data. 6. TIMING FORMULA We now have all of the components necessary to derive a useful timing formula which will describe the motion of the PSR B1257]12 system. First of all let us note that due to the small mass of planet A we will have i A \ i 1 \ m A m i \ i \ B B m B B 2 m ] m ] m m ] A B 0 B m m i \ i \ C B C. (27) C 3 m ] m ] m ] m m ] m ] A B C 0 B C The problem of Ðnding the masses and inclinations of planets B and C can now be formulated as a least-squares problem in which we Ðt to the data the following function: *t \ *t K (() ] dt B (( c B c C ) ] dt C (( c B c C ) (28) ( \ Mx0 n0 e0 u0 T 0 x0 n0 e0 u0 T 0 N B B B B pb C C C C pc c B \ m B c C \ m C (29) *t (() describes the Keplerian part of the motion given by K equation (12) and dt B \ x B 0 G [*h B C3 2 ] 1 2 cos (2 B 0)D ] 1 2 *k B sin (2 B 0) ] *a B a B 0 sin B 0]* B cos B 0H dt C \ x C 0 G [*h C C3 2 ] 1 2 cos (2 C 0)D ] 1 2 *k C sin (2 C 0) ] *a C a C 0 sin C 0]* C cos C 0H (30) describe its non-keplerian part. In this model the initial values of osculating orbital elements replace the Keplerian elements as parameters of the model and we additionally have the parameters c and c related to the masses of B B C and C. We also need the derivatives of *t with respect to the parameters M( c c N. With sufficient accuracy the derivatives with respect B C to ( can be computed from the Keplerian part *t of the model only and the derivatives with respect to c K and c can be easily obtained as dt and dt are proportional B to them: C B C L*t Lc B \ dt C c B L*t Lc C \ dt B c C. (31) Now determination of the masses of B and C and inclinations i and i can be carried out in the following way. First we B assume C that the mass of the pulsar m is canonical (i.e. 1.4 M ) and derive m and m directly fro c and c. Subsequently _ inclinations B can be C computed from B the following C equations: \ MB CN. cx 0\i a 0 sin i 0 G 1 [ i \ n 02 a 03 (32)

5 No PSR B1257]12 PLANETARY SYSTEM NUMERICAL TESTS AND CONCLUSIONS We start the tests with the comparison of our approach for computing osculating elements with that of Malhotra (1993a). It turns out that the most signiðcant component of the non-keplerian motion comes from the changes of mean longitude. Therefore in Figure 2 we present * for the case of coplanar edge-on orbits. As we can see the approach developed in this paper is in fact as good as the integration of full equations of motion. For other elements we obtain very similar results. It should be mentioned that because the model of Malhotra (1993a) is not accurate enough in predicting the secular change of * and because * are proportional to c and c using such a model would result in a signiðcant B error in C determination of the planetsï masses. Next we compared the TOA residuals due to the non-keplerian part of the motion calculated from our model with those obtained by means of numerical integration of full equations of motion. In Figure 3 as one can see for small relative inclinations (Figs. 3a and 3b) our model is very accurate. For larger I (Figs. 3c and 3d) one can see small di erences but this is consistent with the assumptions made to obtain the model. Finally we performed the tests which show how the derived timing model can be used to obtain the masses of planets. To this end we used the TEMPO software package4 modiðed to include our model of the non- Keplerian motion of PSR B1257]12. We simulated two di erent sets of artiðcial TOAs assuming the parameters of PSR B1257]12 (Wolszczan 1994). The Ðrst set of TOAs sampled every day covered a period of 10 years. We also added Gaussian noise with p \ 0.1 ks. The second set was prepared to resemble the real timing observations of PSR B1257]12. Under such conditions we simulated TOAs for the cases described in Figure 3. Subsequently we applied the modiðed TEMPO to obtain the masses of planets B and C. The results are presented in Table 1. The tests performed on Set 1 were used to establish a limit on applicability of our model. As one can see for the relative inclinations I of about 10 the model becomes inaccurate and the relative error of determination of masses is at the level of 20%. Because the relative inclination weakens the interactions of planets using the model which assumes a small I results in such situations in the determined masses that are smaller than the real ones. From the tests performed on set 2 we also learn that this error is bigger than uncertainty originating in TOA measurement errors. On the other hand when I is relatively small we obtain a very 4 See FIG. 2.ÈTime evolution of the mean longitude *(t) \ (t) [ 0(t) for the planets B and C in the case of coplanar edge-on orbits. The solid line indicates the solution by means of the numerical integration of full equations of motion the open circles indicate the solution obtained with the approach presented in this paper and the dash-dotted line with the model by Malhotra (1993a). accurate determination of masses thus proving that our model can be successfully applied as long as the assumptions made to derive it are satisðed. To sum up our model can be applied in all the cases when the observational ÏÏ uncertainties of the planetsï masses are bigger than the uncertainties resulting from the assumptions made to derive the timing model. This in general depends on masses of planets and the relative inclination of the orbits I as well as the characteristic of the observations but one can estimate that in the case of the PSR B1257]12 data I should be smaller than 10. We should also mention that for non-coplanar orbits the angle q will be in general di erent from zero therefore we have to Ðnd it in order to get the proper values of the planetsï masses. It can be done by computing the least-squares value of s2 for a range of q and then choosing the q that corresponds to the minimum of s2. And eventually one has to remember that the PSR B1257]12 planetary system could be in such conðguration that a more signiðcant alteration must be done to our model. First of all the relative inclination of the orbits could be large. In such case one would have to use the full form of the perturbative Hamiltonian H which would lead to a much more complicated model with 1 four additional parameters ) ) i i (instead of one q). Second the presence of a massive B C distant B C planet if signiðcant for the theory of motion would have to be taken into account. TABLE 1 ASSUMED AND DERIVED PARAMETERS OF SIMULATIONS ASSUMED SET 1 SET 2 q I m B m C m B m C m B m C SIMULATION (deg) (deg) (M^) (M^) (M^) (M^) (M^) (M^) a (1) 2.82(1) 3.53(51) 2.63(51) b (1) 4.26(1) 5.07(51) 4.06(48) c (3) 4.11(3) 4.06(51) 3.79(48) d (12) 13.39(12) 8.25(45) 13.16(42) NOTE.ÈValues in parentheses are 3 p uncertainties in the last digits quoted.

6 926 KONACKI MACIEJEWSKI & WOLSZCZAN FIG. 3.ÈTOA residuals due to the non-keplerian part of motion of the PSR B1257]12 in four di erent conðgurations. The solution obtained by means of the numerical integration of full equations of motion is indicated with the dots and the one obtained using our model with the solid line. Brouwer D. & Clemence G. M Methods of Celestial Mechanics (New York: Academic Press) Lyne A. G. & Graham-Smith F Pulsar Astronomy (New York: Cambridge Univ. Press) Malhotra R. 1993a ApJ ÈÈÈ. 1993b in ASP Conf. Ser. 36 Planets around Pulsars ed. J. A. Phillips S. E. Thorsett & S. R. Kulkarni (San Francisco: ASP) 89 Malhotra R. Black D. Eck A. & Jackson A Nature REFERENCES Peale S. J AJ Rasio F. A. Nicholson P. D. Shapiro S. L. & Teukolsky S. A Nature ÈÈÈ in ASP Conf. Ser. 36 Planets around Pulsars ed. J. A. Phillips S. E. Thorsett & S. R. Kulkarni (San Francisco: ASP) 107 Wolszczan A Science Wolszczan A. & Frail D. A Nature

arxiv:astro-ph/ v1 28 May 2003

arxiv:astro-ph/ v1 28 May 2003 Masses and orbital inclinations of planets in the PSR B1257+12 system Maciej Konacki 1,2 and arxiv:astro-ph/0305536v1 28 May 2003 Alex Wolszczan 3,4 ABSTRACT We present measurements of the true masses

More information

imin...

imin... Pulsar Timing For a detailed look at pulsar timing and other pulsar observing techniques, see the Handbook of Pulsar Astronomy by Duncan Lorimer and Michael Kramer. Pulsars are intrinsically interesting

More information

This research was done while the author was a Sta Scientist at the Lunar. and Planetary Institute which is operated by the Universities Space Research

This research was done while the author was a Sta Scientist at the Lunar. and Planetary Institute which is operated by the Universities Space Research ACKNOWLEDGMENTS This research was done while the author was a Sta Scientist at the Lunar and Planetary Institute which is operated by the Universities Space Research Association under contract no. NASW-4574

More information

arxiv: v1 [astro-ph.im] 3 Oct 2017

arxiv: v1 [astro-ph.im] 3 Oct 2017 Starovoit E.D. Rodin A.E. On The Existence of Planets Around the Pulsar PSR B0329+54 E. D. Starovoit 1, 2, * 1, ** and A. E. Rodin 1 Pushchino Radio Astronomy Observatory, arxiv:1710.01153v1 [astro-ph.im]

More information

Orbital Evolution in Extra-solar systems

Orbital Evolution in Extra-solar systems Orbital Evolution in Extra-solar systems George Voyatzis Section of Astronomy, Astrophysics and Mechanics, Department of Physics, Aristotle University of Thessaloniki, Greece. Abstract Nowadays, extra-solar

More information

Dynamics of two planets in the 3/2 mean-motion resonance: application to the planetary system of the pulsar PSR B

Dynamics of two planets in the 3/2 mean-motion resonance: application to the planetary system of the pulsar PSR B Celestial Mech Dyn Astr (2006) 94:381 397 DOI 10.1007/s10569-006-9002-4 ORIGINAL ARTICLE Dynamics of two planets in the 3/2 mean-motion resonance: application to the planetary system of the pulsar PSR

More information

Analysis of Radial Velocity Measurements

Analysis of Radial Velocity Measurements Analysis of Radial Velocity Measurements Sistemas Planetários - 2nd project - 2011/2012 João Faria 1. Introduction The existence of solar systems other than our own has been speculated for a long time,

More information

CALCULATION OF POSITION AND VELOCITY OF GLONASS SATELLITE BASED ON ANALYTICAL THEORY OF MOTION

CALCULATION OF POSITION AND VELOCITY OF GLONASS SATELLITE BASED ON ANALYTICAL THEORY OF MOTION ARTIFICIAL SATELLITES, Vol. 50, No. 3 2015 DOI: 10.1515/arsa-2015-0008 CALCULATION OF POSITION AND VELOCITY OF GLONASS SATELLITE BASED ON ANALYTICAL THEORY OF MOTION W. Góral, B. Skorupa AGH University

More information

Probing Relativistic Gravity with the Double Pulsar

Probing Relativistic Gravity with the Double Pulsar Probing Relativistic Gravity with the Double Pulsar Marta Burgay INAF Osservatorio Astronomico di Cagliari The spin period of the original millisecond pulsar PSR B1937+21: P = 0.0015578064924327 ± 0.0000000000000004

More information

Eccentricity pumping of a planet on an inclined orbit by a disc

Eccentricity pumping of a planet on an inclined orbit by a disc Mon. Not. R. Astron. Soc. 44, 49 414 21) doi:1.1111/j.1365-2966.21.16295.x Eccentricity pumping of a planet on an inclined orbit by a disc Caroline Terquem 1,2 and Aikel Ajmia 1 1 Institut d Astrophysique

More information

Habitability in the Upsilon Andromedae System

Habitability in the Upsilon Andromedae System Habitability in the Upsilon Andromedae System Adrienne Dove University of Missouri Columbia Institute for Astronomy University of Hawaii Mentor: Nader Haghighipour ABSTRACT We investigate the habitability

More information

GEODETIC PRECESSION AND TIMING OF THE RELATIVISTIC BINARY PULSARS PSR B AND PSR B Maciej Konacki. Alex Wolszczan. and Ingrid H.

GEODETIC PRECESSION AND TIMING OF THE RELATIVISTIC BINARY PULSARS PSR B AND PSR B Maciej Konacki. Alex Wolszczan. and Ingrid H. The Astrophysical Journal, 589:495 502, 2003 May 20 # 2003. The American Astronomical Society. All rights reserved. Printed in U.S.A. GEODETIC PRECESSION AND TIMING OF THE RELATIVISTIC BINARY PULSARS PSR

More information

arxiv:astro-ph/ v1 7 Jul 2004

arxiv:astro-ph/ v1 7 Jul 2004 Binary Radio Pulsars ASP Conference Series, Vol. TBD, 2004 eds. F.A. Rasio & I.H. Stairs Relativistic Binary Pulsar B1913+16: Thirty Years of Observations and Analysis Joel M. Weisberg Dept. of Physics

More information

Joseph Castro Mentor: Nader Haghighipour

Joseph Castro Mentor: Nader Haghighipour ON THE POSSIBILITY OF ADDITIONAL PLANETS IN THE γ CEPHEI BINARY-PLANETARY SYSTEM Joseph Castro Mentor: Nader Haghighipour ABSTRACT Results of the simulations of the dynamical stability of additional hypothetical

More information

Dynamical properties of the Solar System. Second Kepler s Law. Dynamics of planetary orbits. ν: true anomaly

Dynamical properties of the Solar System. Second Kepler s Law. Dynamics of planetary orbits. ν: true anomaly First Kepler s Law The secondary body moves in an elliptical orbit, with the primary body at the focus Valid for bound orbits with E < 0 The conservation of the total energy E yields a constant semi-major

More information

Triple in M Analysis of the Pulsar Timing Data

Triple in M Analysis of the Pulsar Timing Data Triple in M4 2 sponding orbital periods in the range 10 2 10 3 yr (Michel 1994; Rasio 1994; Sigurdsson 1995). More recent calculations using four frequency derivatives and preliminary measurements of the

More information

2 Ford, Rasio, & Yu. 2. Two Planets, Unequal Masses

2 Ford, Rasio, & Yu. 2. Two Planets, Unequal Masses 2 Ford, Rasio, & Yu unlikely to have developed such a large eccentricity, since dissipation in the disk tends to circularize orbits. Dynamical instabilities leading to the ejection of one planet while

More information

h librate about 180. By considering the condition that -5 for stable simultaneous

h librate about 180. By considering the condition that -5 for stable simultaneous THE ASTROPHYSICAL JOURNAL, 567:596È609, 00 March ( 00. The American Astronomical Society. All rights reserved. Printed in U.S.A. DYNAMICS AND ORIGIN OF THE : ORBITAL RESONANCES OF THE GJ 876 PLANETS MAN

More information

Radio timing observations of the pulsar by Kaspi et al. (1994) have

Radio timing observations of the pulsar by Kaspi et al. (1994) have TIMING OBSERVATIONS OF THE SMC BINARY PSR J0045?7319 V. M. KASPI IPAC/Caltech/Jet Propulsion Laboratory Pasadena, CA, USA 91125 R. N. MANCHESTER AND M. BAILES ATNF/CSIRO Epping, Australia AND J. F. BELL

More information

orbits Moon, Planets Spacecrafts Calculating the and by Dr. Shiu-Sing TONG

orbits Moon, Planets Spacecrafts Calculating the and by Dr. Shiu-Sing TONG A Science Enrichment Programme for Secondary 3-4 Students : Teaching and Learning Resources the and Spacecrafts orbits Moon, Planets Calculating the 171 of by Dr. Shiu-Sing TONG 172 Calculating the orbits

More information

Binary Pulsars and Evidence for Gravitational Radiation

Binary Pulsars and Evidence for Gravitational Radiation Binary Pulsars and Evidence for Gravitational Radiation Matthew S. Paoletti Institute for Research in Electronics and Applied Physics Energy Research Facility, Bldg. #3 University of Maryland College Park,

More information

arxiv:astro-ph/ v1 16 Apr 2002

arxiv:astro-ph/ v1 16 Apr 2002 Interstellar Scintillation Velocities of a Relativistic Binary PSR B1534+12 and Three Other Millisecond Pulsars Slavko Bogdanov arxiv:astro-ph/0204270v1 16 Apr 2002 Department of Astronomy and Astrophysics,

More information

Binary Pulsars. By: Kristo er Hultgren. Karlstad University Faculty of Technology and Science Department of Physics

Binary Pulsars. By: Kristo er Hultgren. Karlstad University Faculty of Technology and Science Department of Physics Binary Pulsars By: Kristo er Hultgren Karlstad University Faculty of Technology and Science Department of Physics Course: Analytical Mechanics (FYSCD2) Examiner: Prof. Jürgen Fuchs January 5, 2007 Abstract

More information

LONG-TERM STABILITY AND DYNAMICAL ENVIRONMENT OF THE PSR PLANETARY SYSTEM

LONG-TERM STABILITY AND DYNAMICAL ENVIRONMENT OF THE PSR PLANETARY SYSTEM The Astrophysical Journal, 619:1084 1097, 2005 February 1 # 2005. The American Astronomical Society. All rights reserved. Printed in U.S.A. A LONG-TERM STABILITY AND DYNAMICAL ENVIRONMENT OF THE PSR 1257+12

More information

A study upon Eris. I. Describing and characterizing the orbit of Eris around the Sun. I. Breda 1

A study upon Eris. I. Describing and characterizing the orbit of Eris around the Sun. I. Breda 1 Astronomy & Astrophysics manuscript no. Eris c ESO 2013 March 27, 2013 A study upon Eris I. Describing and characterizing the orbit of Eris around the Sun I. Breda 1 Faculty of Sciences (FCUP), University

More information

The Use of Transit Timing to Detect Extrasolar Planets with Masses as Small as Earth

The Use of Transit Timing to Detect Extrasolar Planets with Masses as Small as Earth arxiv:astro-ph/41228v1 1 Dec 24 The Use of Transit Timing to Detect Extrasolar Planets with Masses as Small as Earth Matthew J. Holman, 1 Norman W. Murray 2,3 1 Harvard-Smithsonian Center for Astrophysics,

More information

arxiv: v1 [astro-ph] 28 May 2008

arxiv: v1 [astro-ph] 28 May 2008 Submitted to the The Astrophysical Journal Letters Possibility of Detecting Moons of Pulsar Planets Through Time-of-Arrival Analysis arxiv:0805.4263v1 [astro-ph 28 May 2008 Karen M. Lewis School of Mathematical

More information

AST111, Lecture 1b. Measurements of bodies in the solar system (overview continued) Orbital elements

AST111, Lecture 1b. Measurements of bodies in the solar system (overview continued) Orbital elements AST111, Lecture 1b Measurements of bodies in the solar system (overview continued) Orbital elements Planetary properties (continued): Measuring Mass The orbital period of a moon about a planet depends

More information

arxiv:astro-ph/ v1 15 Aug 2001

arxiv:astro-ph/ v1 15 Aug 2001 A test of general relativity from the three-dimensional orbital geometry of a binary pulsar arxiv:astro-ph/0108254 v1 15 Aug 2001 W. van Straten, M. Bailes, M.C. Britton, S.R. Kulkarni, S.B. Anderson,

More information

Planetenbewegung in Sternsystemen. The Effect of Resonances. Part 2

Planetenbewegung in Sternsystemen. The Effect of Resonances. Part 2 Planetenbewegung in Sternsystemen The Effect of Resonances Part 2 Topics overview 1. Definition and examples of resonances 2. Disturbing function 3. Mean-motion resonance (MMR) 4. Secular resonance (SR)

More information

discovers a radio-quiet gamma-ray millisecond Journal Group

discovers a radio-quiet gamma-ray millisecond Journal Group Einstein@Home discovers a radio-quiet gamma-ray millisecond pulsar @CHEA Journal Group Contribution of the paper Contribution of the paper Millisecond Pulsars (MSPs) Ver y rapid rotating neutron star

More information

Restricted three body problems in the Solar System: simulations

Restricted three body problems in the Solar System: simulations Author:. Facultat de Física, Universitat de Barcelona, Diagonal 645, 0808 Barcelona, Spain. Advisor: Antoni Benseny i Ardiaca. Facultat de Matemàtiques, Universitat de Barcelona, Gran Via de les Corts

More information

Planetary system dynamics Part III Mathematics / Part III Astrophysics

Planetary system dynamics Part III Mathematics / Part III Astrophysics Planetary system dynamics Part III Mathematics / Part III Astrophysics Lecturer: Prof. Mark Wyatt (Dr. Amy Bonsor on 9,11 Oct) Schedule: Michaelmas 2017 Mon, Wed, Fri at 10am MR11, 24 lectures, start Fri

More information

Importance of the study of extrasolar planets. Exoplanets Introduction. Importance of the study of extrasolar planets

Importance of the study of extrasolar planets. Exoplanets Introduction. Importance of the study of extrasolar planets Importance of the study of extrasolar planets Exoplanets Introduction Planets and Astrobiology (2017-2018) G. Vladilo Technological and scientific spin-offs Exoplanet observations are driving huge technological

More information

arxiv:astro-ph/ v1 23 Oct 2000

arxiv:astro-ph/ v1 23 Oct 2000 Mon. Not. R. Astron. Soc. 000, 000 000 (0000) Printed 27 October 2018 (MN LATEX style file v1.4) Determination of the orbital parameters of binary pulsars P. C. Freire, M. Kramer, A. G. Lyne arxiv:astro-ph/0010463v1

More information

{ 2{ as well as the structure and dynamics of globular clusters. In addition, the study of individual exotic binaries formed in the dense central core

{ 2{ as well as the structure and dynamics of globular clusters. In addition, the study of individual exotic binaries formed in the dense central core Observations of two Millisecond Pulsars in the Globular Cluster NGC5904 S. B. Anderson 1 and A. Wolszczan Department of Astronomy and Astrophysics, Pennsylvania State University, University Park, PA 16802

More information

Regular n-gon as a model of discrete gravitational system. Rosaev A.E. OAO NPC NEDRA, Jaroslavl Russia,

Regular n-gon as a model of discrete gravitational system. Rosaev A.E. OAO NPC NEDRA, Jaroslavl Russia, Regular n-gon as a model of discrete gravitational system Rosaev A.E. OAO NPC NEDRA, Jaroslavl Russia, E-mail: hegem@mail.ru Introduction A system of N points, each having mass m, forming a planar regular

More information

Corralling a distant planet with extreme resonant Kuiper belt objects

Corralling a distant planet with extreme resonant Kuiper belt objects Corralling a distant planet with extreme resonant Kuiper belt objects Renu Malhotra 1 Kathryn Volk Xianyu Wang 2 arxiv:1603.02196v3 [astro-ph.ep] 2 Jun 2016 Lunar and Planetary Laboratory, The University

More information

Extrasolar Planets. Methods of detection Characterization Theoretical ideas Future prospects

Extrasolar Planets. Methods of detection Characterization Theoretical ideas Future prospects Extrasolar Planets Methods of detection Characterization Theoretical ideas Future prospects Methods of detection Methods of detection Methods of detection Pulsar timing Planetary motion around pulsar

More information

arxiv: v1 [astro-ph.ep] 3 Apr 2018

arxiv: v1 [astro-ph.ep] 3 Apr 2018 Astronomy& Astrophysics manuscript no. zanardi_ c ESO 28 September 2, 28 The role of the general relativity on icy body reservoirs under the effects of an inner eccentric Jupiter M. Zanardi, 2, G. C. de

More information

SECULAR EVOLUTION OF HIERARCHICAL PLANETARY SYSTEMS Man Hoi Lee and S. J. Peale

SECULAR EVOLUTION OF HIERARCHICAL PLANETARY SYSTEMS Man Hoi Lee and S. J. Peale The Astrophysical Journal, 59:0 6, 003 August # 003. The American Astronomical Society. All rights reserved. Printed in U.S.A. SECULAR EVOLUTION OF HIERARCHICAL PLANETARY SYSTEMS Man Hoi Lee and S. J.

More information

The Transit Timing method: How to detect extrasolar planets? Jean Coupon. (Supervisor: Norman W. Murray)

The Transit Timing method: How to detect extrasolar planets? Jean Coupon. (Supervisor: Norman W. Murray) The Transit Timing method: How to detect extrasolar planets? Jean Coupon (Supervisor: Norman W. Murray) Canadian Institute for Theoretical Astrophysics Institut d Astrophysique de Paris June 26th, 2006

More information

Measurements of Neutron Star Masses with a strong emphasis on millisecond binary radio pulsar timing

Measurements of Neutron Star Masses with a strong emphasis on millisecond binary radio pulsar timing Measurements of Neutron Star Masses with a strong emphasis on millisecond binary radio pulsar timing David Nice, Lafayette College Physics of Neutron Stars 2014, Saint Petersburg 1. Motivation 2. How to

More information

On the Secular Evolution of Extrasolar Planetary Systems

On the Secular Evolution of Extrasolar Planetary Systems On the Secular Evolution of Extrasolar Planetary Systems Marco Sansottera [a] [a] Namur Center for Complex Systems (naxys) Based on a research work in collaboration with Anne-Sophie Libert [a] JISD 2012

More information

arxiv: v1 [astro-ph.sr] 30 Nov 2009

arxiv: v1 [astro-ph.sr] 30 Nov 2009 Compact stars in the QCD phase diagram II (CSQCD II) May 20-24, 2009, KIAA at Peking University, Beijing - P. R. China http://vega.bac.pku.edu.cn/rxxu/csqcd.htm The timing behaviour of radio pulsars arxiv:0911.5534v1

More information

THE PLANE OF THE KUIPER BELT

THE PLANE OF THE KUIPER BELT The Astronomical Journal, 127:2418 2423, 2004 April # 2004. The American Astronomical Society. All rights reserved. Printed in U.S.A. THE PLANE OF THE KUIPER BELT Michael E. Brown Division of Geological

More information

Spacecraft Dynamics and Control

Spacecraft Dynamics and Control Spacecraft Dynamics and Control Matthew M. Peet Arizona State University Lecture 5: Hyperbolic Orbits Introduction In this Lecture, you will learn: Hyperbolic orbits Hyperbolic Anomaly Kepler s Equation,

More information

The Restricted 3-Body Problem

The Restricted 3-Body Problem The Restricted 3-Body Problem John Bremseth and John Grasel 12/10/2010 Abstract Though the 3-body problem is difficult to solve, it can be modeled if one mass is so small that its effect on the other two

More information

5.1. Accelerated Coordinate Systems:

5.1. Accelerated Coordinate Systems: 5.1. Accelerated Coordinate Systems: Recall: Uniformly moving reference frames (e.g. those considered at 'rest' or moving with constant velocity in a straight line) are called inertial reference frames.

More information

Observational Astronomy - Lecture 4 Orbits, Motions, Kepler s and Newton s Laws

Observational Astronomy - Lecture 4 Orbits, Motions, Kepler s and Newton s Laws Observational Astronomy - Lecture 4 Orbits, Motions, Kepler s and Newton s Laws Craig Lage New York University - Department of Physics craig.lage@nyu.edu February 24, 2014 1 / 21 Tycho Brahe s Equatorial

More information

THE exploration of planetary satellites is currently an active area

THE exploration of planetary satellites is currently an active area JOURNAL OF GUIDANCE, CONTROL, AND DYNAMICS Vol. 9, No. 5, September October 6 Design of Science Orbits About Planetary Satellites: Application to Europa Marci E. Paskowitz and Daniel J. Scheeres University

More information

Chapter 2 Introduction to Binary Systems

Chapter 2 Introduction to Binary Systems Chapter 2 Introduction to Binary Systems In order to model stars, we must first have a knowledge of their physical properties. In this chapter, we describe how we know the stellar properties that stellar

More information

ANNEX 1. DEFINITION OF ORBITAL PARAMETERS AND IMPORTANT CONCEPTS OF CELESTIAL MECHANICS

ANNEX 1. DEFINITION OF ORBITAL PARAMETERS AND IMPORTANT CONCEPTS OF CELESTIAL MECHANICS ANNEX 1. DEFINITION OF ORBITAL PARAMETERS AND IMPORTANT CONCEPTS OF CELESTIAL MECHANICS A1.1. Kepler s laws Johannes Kepler (1571-1630) discovered the laws of orbital motion, now called Kepler's laws.

More information

Remarks on Quadratic Hamiltonians in Spaceflight Mechanics

Remarks on Quadratic Hamiltonians in Spaceflight Mechanics Remarks on Quadratic Hamiltonians in Spaceflight Mechanics Bernard Bonnard 1, Jean-Baptiste Caillau 2, and Romain Dujol 2 1 Institut de mathématiques de Bourgogne, CNRS, Dijon, France, bernard.bonnard@u-bourgogne.fr

More information

Pulsars. in this talk. Pulsar timing. Pulsar timing. Pulsar timing. Pulsar timing. How to listen to what exotic. are telling us! Paulo César C.

Pulsars. in this talk. Pulsar timing. Pulsar timing. Pulsar timing. Pulsar timing. How to listen to what exotic. are telling us! Paulo César C. How to listen to what exotic Pulsars are telling us! in this talk 1. 2. 3. Test of gravitational theories using binary pulsars 4. Probing the equation of state of super-dense matter Paulo César C. Freire

More information

10/16/ Detecting Planets Around Other Stars. Chapter 10: Other Planetary Systems The New Science of Distant Worlds

10/16/ Detecting Planets Around Other Stars. Chapter 10: Other Planetary Systems The New Science of Distant Worlds 10/16/17 Lecture Outline 10.1 Detecting Planets Around Other Stars Chapter 10: Other Planetary Systems The New Science of Distant Worlds Our goals for learning: How do we detect planets around other stars?

More information

Planetary Perturbations on the 2 : 3 Mean Motion Resonance with Neptune

Planetary Perturbations on the 2 : 3 Mean Motion Resonance with Neptune PASJ: Publ. Astron. Soc. Japan 54, 493 499, 2002 June 25 c 2002. Astronomical Society of Japan. Planetary Perturbations on the 2 : 3 Mean Motion Resonance with Neptune Tetsuharu FUSE Subaru Telescope,

More information

Observations of extrasolar planets

Observations of extrasolar planets Observations of extrasolar planets 1 Mercury 2 Venus radar image from Magellan (vertical scale exaggerated 10 X) 3 Mars 4 Jupiter 5 Saturn 6 Saturn 7 Uranus and Neptune 8 we need to look out about 10 parsecs

More information

Lecture 20: Planet formation II. Clues from Exoplanets

Lecture 20: Planet formation II. Clues from Exoplanets Lecture 20: Planet formation II. Clues from Exoplanets 1 Outline Definition of a planet Properties of exoplanets Formation models for exoplanets gravitational instability model core accretion scenario

More information

1. INTRODUCTION. THE ASTRONOMICAL JOURNAL, 121:1730È1735, 2001 March ( The American Astronomical Society. All rights reserved. Printed in U.S.A.

1. INTRODUCTION. THE ASTRONOMICAL JOURNAL, 121:1730È1735, 2001 March ( The American Astronomical Society. All rights reserved. Printed in U.S.A. THE ASTRONOMICAL JOURNAL, 121:1730È1735, 2001 March ( 2001. The American Astronomical Society. All rights reserved. Printed in U.S.A. ON THE SIZE DEPENDENCE OF THE INCLINATION DISTRIBUTION OF THE MAIN

More information

The Same Physics Underlying SGRs, AXPs and Radio Pulsars

The Same Physics Underlying SGRs, AXPs and Radio Pulsars Chin. J. Astron. Astrophys. Vol. 6 (2006), Suppl. 2, 273 278 (http://www.chjaa.org) Chinese Journal of Astronomy and Astrophysics The Same Physics Underlying SGRs, AXPs and Radio Pulsars Biping Gong National

More information

Indirect Methods: gravitational perturbation of the stellar motion. Exoplanets Doppler method

Indirect Methods: gravitational perturbation of the stellar motion. Exoplanets Doppler method Indirect Methods: gravitational perturbation of the stellar motion Exoplanets The reflex motion of the star is proportional to M p /M * This introduces an observational bias that favours the detection

More information

The two-body Kepler problem

The two-body Kepler problem The two-body Kepler problem set center of mass at the origin (X = 0) ignore all multipole moments (spherical bodies or point masses) define r := r 1 r 2,r:= r,m:= m 1 + m 2,µ:= m 1 m 2 /m reduces to effective

More information

Post-Keplerian effects in binary systems

Post-Keplerian effects in binary systems Post-Keplerian effects in binary systems Laboratoire Univers et Théories Observatoire de Paris / CNRS The problem of binary pulsar timing (Credit: N. Wex) Some classical tests of General Relativity Gravitational

More information

Exoplanet Search Techniques: Overview. PHY 688, Lecture 28 April 3, 2009

Exoplanet Search Techniques: Overview. PHY 688, Lecture 28 April 3, 2009 Exoplanet Search Techniques: Overview PHY 688, Lecture 28 April 3, 2009 Course administration final presentations Outline see me for paper recommendations 2 3 weeks before talk see me with draft of presentation

More information

Phys 2101 Gabriela González

Phys 2101 Gabriela González Phys 2101 Gabriela González Newton s law : F = Gm 1 m 2 /r 2 Explains why apples fall, why the planets move around the Sun, sciencebulletins.amnh.org And in YouTube! Explains just as well as Newtons why

More information

Orbital and Celestial Mechanics

Orbital and Celestial Mechanics Orbital and Celestial Mechanics John P. Vinti Edited by Gim J. Der TRW Los Angeles, California Nino L. Bonavito NASA Goddard Space Flight Center Greenbelt, Maryland Volume 177 PROGRESS IN ASTRONAUTICS

More information

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Physics. Physics 8.901: Astrophysics I Spring Term 2006 PROBLEM SET 1

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Physics. Physics 8.901: Astrophysics I Spring Term 2006 PROBLEM SET 1 MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Physics Physics 8.901: Astrophysics I Spring Term 2006 PROBLEM SET 1 Due: Thursday, February 16 in class Reading: Hansen, Kawaler, & Trimble, Chapter

More information

arxiv: v1 [astro-ph] 8 Dec 2008

arxiv: v1 [astro-ph] 8 Dec 2008 Warping the young stellar disc in the Galactic Centre arxiv:0812.1561v1 [astro-ph] 8 Dec 2008 Ladislav Šubr 1,2 and Jaroslava Schovancová 1 1 Faculty of Mathematics and Physics, Charles University, V Holešovičkách

More information

UNIVERSITY of CALIFORNIA SANTA CRUZ

UNIVERSITY of CALIFORNIA SANTA CRUZ UNIVERSITY of CALIFORNIA SANTA CRUZ NUMERICAL AND ANALYTIC RESEARCH INTO THE FORMATION OF THE HD 80606B PLANETARY SYSTEM A thesis submitted in partial satisfaction of the requirements for the degree of

More information

22 Years of a Pulsar-Be Binary System: From Parkes to the Heavens (Fermi) Ryan Shannon Postdoctoral Fellow, CSIRO Astronomy and Space Science

22 Years of a Pulsar-Be Binary System: From Parkes to the Heavens (Fermi) Ryan Shannon Postdoctoral Fellow, CSIRO Astronomy and Space Science Credit: CSIRO/John Sarkissian 22 Years of a Pulsar-Be Binary System: From Parkes to the Heavens (Fermi) Ryan Shannon Postdoctoral Fellow, CSIRO Astronomy and Space Science Big Picture and Talk Outline

More information

Chapter 13 Other Planetary Systems. The New Science of Distant Worlds

Chapter 13 Other Planetary Systems. The New Science of Distant Worlds Chapter 13 Other Planetary Systems The New Science of Distant Worlds 13.1 Detecting Extrasolar Planets Our goals for learning Why is it so difficult to detect planets around other stars? How do we detect

More information

Secular dynamics of extrasolar-systems

Secular dynamics of extrasolar-systems Secular dynamics of extrasolar-systems Marco Sansottera [1] [1] Università degli Studi di Milano [2] Namur Center for Complex Systems (naxys) Based on a research project in collaboration with A.-S. Libert

More information

Velocity Curve Analysis of the Spectroscopic Binary Stars PV Pup, HD , EE Cet and V921 Her by Nonlinear Regression

Velocity Curve Analysis of the Spectroscopic Binary Stars PV Pup, HD , EE Cet and V921 Her by Nonlinear Regression J. Astrophys. Astr. (2007) 28, 217 230 Velocity Curve Analysis of the Spectroscopic Binary Stars PV Pup, HD 141929, EE Cet and V921 Her by Nonlinear Regression K. Karami 1,2,3, & R. Mohebi 1, 1 Department

More information

arxiv:astro-ph/ v3 2 Mar 2006

arxiv:astro-ph/ v3 2 Mar 2006 Mon. Not. R. Astron. Soc. 000, 1?? (2006) Printed 5 February 2008 (MN LATEX style file v2.2) Pluto s moon system: survey of the phase space I. I. Nagy, Á. Süli and B. Érdi Department of Astronomy, Eötvös

More information

ISIMA lectures on celestial mechanics. 1

ISIMA lectures on celestial mechanics. 1 ISIMA lectures on celestial mechanics. 1 Scott Tremaine, Institute for Advanced Study July 2014 The roots of solar system dynamics can be traced to two fundamental discoveries by Isaac Newton: first, that

More information

ISIMA lectures on celestial mechanics. 3

ISIMA lectures on celestial mechanics. 3 ISIMA lectures on celestial mechanics. 3 Scott Tremaine, Institute for Advanced Study July 2014 1. The stability of planetary systems To understand the formation and evolution of exoplanet systems, we

More information

Chapter 13: universal gravitation

Chapter 13: universal gravitation Chapter 13: universal gravitation Newton s Law of Gravitation Weight Gravitational Potential Energy The Motion of Satellites Kepler s Laws and the Motion of Planets Spherical Mass Distributions Apparent

More information

Advanced Newtonian gravity

Advanced Newtonian gravity Foundations of Newtonian gravity Solutions Motion of extended bodies, University of Guelph h treatment of Newtonian gravity, the book develops approximation methods to obtain weak-field solutions es the

More information

Targeting Mars. 19th June University of Ljubljana. Faculty of Mathematic and Physics. Seminar: 4. year-university program

Targeting Mars. 19th June University of Ljubljana. Faculty of Mathematic and Physics. Seminar: 4. year-university program University of Ljubljana Faculty of Mathematic and Physics Targeting Mars Seminar: 4. year-university program 19th June 016 Author: Jo²t Vrabi Koren Advisor: Prof. Dr. Tomaº Zwitter Ljubljana, June 016

More information

Chapter 13 Lecture. The Cosmic Perspective. Seventh Edition. Other Planetary Systems: The New Science of Distant Worlds Pearson Education, Inc.

Chapter 13 Lecture. The Cosmic Perspective. Seventh Edition. Other Planetary Systems: The New Science of Distant Worlds Pearson Education, Inc. Chapter 13 Lecture The Cosmic Perspective Seventh Edition Other Planetary Systems: The New Science of Distant Worlds 13.1 Detecting Planets Around Other Stars Our goals for learning: Why is it so challenging

More information

Dynamically Unstable Planetary Systems Emerging Out of Gas Disks

Dynamically Unstable Planetary Systems Emerging Out of Gas Disks EXTREME SOLAR SYSTEMS ASP Conference Series, Vol. 398, 2008 D. Fischer, F. A. Rasio, S. E. Thorsett, and A. Wolszczan Dynamically Unstable Planetary Systems Emerging Out of Gas Disks Soko Matsumura, Edward

More information

Secular Evolution of Extrasolar Planetary Systems:

Secular Evolution of Extrasolar Planetary Systems: Secular Evolution of Extrasolar Planetary Systems: an Extension of the Laplace-Lagrange Secular Theory Marco Sansottera [a] [a] Namur Center for Complex Systems (naxys) Based on a research work in collaboration

More information

Lecture 22: Gravitational Orbits

Lecture 22: Gravitational Orbits Lecture : Gravitational Orbits Astronomers were observing the motion of planets long before Newton s time Some even developed heliocentric models, in which the planets moved around the sun Analysis of

More information

THE PARKES MULTIBEAM PULSAR SURVEY

THE PARKES MULTIBEAM PULSAR SURVEY THE PARKES MULTIBEAM PULSAR SURVEY I. H. STAIRS, A. G. LYNE, F. CAMILO, N. P. F. MCKAY, D. C. SHEPPARD, D. J. MORRIS University of Manchester, Nuffield Radio Astronomy Laboratories, Jodrell Bank, Macclesfield,

More information

Planets are plentiful

Planets are plentiful Extra-Solar Planets Planets are plentiful The first planet orbiting another Sun-like star was discovered in 1995. We now know of 209 (Feb 07). Including several stars with more than one planet - true planetary

More information

How Kepler discovered his laws

How Kepler discovered his laws How Kepler discovered his laws A. Eremenko December 17, 2016 This question is frequently asked on History of Science and Math, by people who know little about astronomy and its history, so I decided to

More information

An Introduction to Celestial Mechanics

An Introduction to Celestial Mechanics An Introduction to Celestial Mechanics This accessible text on classical celestial mechanics the principles governing the motions of bodies in the solar system provides a clear and concise treatment of

More information

Dynamical Stability of Terrestrial and Giant Planets in the HD Planetary System

Dynamical Stability of Terrestrial and Giant Planets in the HD Planetary System Dynamical Stability of Terrestrial and Giant Planets in the HD 155358 Planetary System James Haynes Advisor: Nader Haghighipour ABSTRACT The results of a study of the dynamical evolution and the habitability

More information

arxiv:astro-ph/ v1 17 Jul 2003

arxiv:astro-ph/ v1 17 Jul 2003 Astronomy & Astrophysics manuscript no. RB-letter September 25, 217 (DOI: will be inserted by hand later) Mercury s spin-orbit model and signature of its dynamical parameters arxiv:astro-ph/37328v1 17

More information

Dynamic Exoplanets. Alexander James Mustill

Dynamic Exoplanets. Alexander James Mustill Dynamic Exoplanets Alexander James Mustill Exoplanets: not (all) like the Solar System Exoplanets: not (all) like the Solar System Solar System Lissauer et al 14 Key questions to bear in mind What is role

More information

arxiv:astro-ph/ v1 14 Sep 2005

arxiv:astro-ph/ v1 14 Sep 2005 For publication in Bayesian Inference and Maximum Entropy Methods, San Jose 25, K. H. Knuth, A. E. Abbas, R. D. Morris, J. P. Castle (eds.), AIP Conference Proceeding A Bayesian Analysis of Extrasolar

More information

Chapter 13 Other Planetary Systems. Why is it so difficult to detect planets around other stars? Size Difference. Brightness Difference

Chapter 13 Other Planetary Systems. Why is it so difficult to detect planets around other stars? Size Difference. Brightness Difference Chapter 13 Other Planetary Systems Why is it so difficult to detect planets around other stars? Size Difference Planets are small compared to interstellar distances 10 billion to 1 scale Sun is size of

More information

Chapter 13 Other Planetary Systems. Why is it so difficult to detect planets around other stars? Brightness Difference

Chapter 13 Other Planetary Systems. Why is it so difficult to detect planets around other stars? Brightness Difference Chapter 13 Other Planetary Systems The New Science of Distant Worlds 13.1 Detecting Extrasolar Planets Our goals for learning:! Why is it so difficult to detect planets around other stars?! How do we detect

More information

Astrodynamics (AERO0024)

Astrodynamics (AERO0024) Astrodynamics (AERO0024) 3B. The Orbit in Space and Time Gaëtan Kerschen Space Structures & Systems Lab (S3L) Previous Lecture: The Orbit in Time 3.1 ORBITAL POSITION AS A FUNCTION OF TIME 3.1.1 Kepler

More information

On the uncertainty of the minimal distance between two confocal Keplerian orbits

On the uncertainty of the minimal distance between two confocal Keplerian orbits On the uncertainty of the minimal distance between two confocal Keplerian orbits Giovanni F. Gronchi and Giacomo Tommei Dipartimento di Matematica, Università di Pisa Largo B. Pontecorvo 5, 56127 Pisa,

More information

Gravity with the SKA

Gravity with the SKA Gravity with the SKA Strong-field tests of gravity using Pulsars and Black Holes Michael Kramer Jodrell Bank Observatory University of Manchester With Don Backer, Jim Cordes, Simon Johnston, Joe Lazio

More information

General Relativity Tests with Pulsars

General Relativity Tests with Pulsars General Relativity Tests with Pulsars Ingrid Stairs UBC SLAC Summer Institute July 27, 2005 Much of this material is in Living Reviews in Relativity 2003 5. Pulsars: rotating, magnetized neutron stars.

More information

TROJAN PAIRS IN THE HD AND HD PLANETARY SYSTEMS?

TROJAN PAIRS IN THE HD AND HD PLANETARY SYSTEMS? The Astrophysical Journal, 647:573Y586, 2006 August 10 # 2006. The American Astronomical Society. All rights reserved. Printed in U.S.A. A TROJAN PAIRS IN THE HD 128311 AND HD 82943 PLANETARY SYSTEMS?

More information

Resonance Capture. Alice Quillen. University of Rochester. Isaac Newton Institute Dec 2009

Resonance Capture. Alice Quillen. University of Rochester. Isaac Newton Institute Dec 2009 Resonance Capture Alice Quillen University of Rochester Isaac Newton Institute Dec 2009 This Talk Resonance Capture extension to nonadiabatic regime Astrophysical settings: Dust spiraling inward via radiation

More information