A Bayesian Analysis of Extrasolar Planet Data for HD 73526

Total Page:16

File Type:pdf, Size:1020Kb

A Bayesian Analysis of Extrasolar Planet Data for HD 73526 Astrophysical Journal, 22 Dec. 2004, revised May 2005 A Bayesian Analysis of Extrasolar Planet Data for HD 73526 P. C. Gregory Physics and Astronomy Department, University of British Columbia, Vancouver, BC V6T 1Z1 ABSTRACT A new Bayesian tool for nonlinear model fitting has been developed, that employs a parallel tempering Markov Chain Monte Carlo algorithm with a novel statistical control system. The algorithm has been used to re-analyze precision radial velocity data for HD 73526. For a single planet model, three possible +0.37 +1.8 +1.4 orbits were found with periods of 127.88−0.09, 190.4−2.1 and 376.2−4.3 days. The 128 day orbit, with an eccentricity of 0.56 0.03, has a maximum value of prior × likelihood which is 16 times larger than for the next highest solution at 376 +0.05 days. However, the 376 day orbit, with an eccentricity of 0.10−0.10, is formally more probable because for this sparse data set there is a much larger volume of parameter space with a significant probability density in the vicinity of the 376 day peak. The previously reported orbit (Tinney et al. 2003) of 190.5 3.0 days corresponds to our least probable orbit. The analysis highlights the need for measurements around phase 0.5 for the 376 day period. Subject headings: planetary systems — methods: data analysis — stars: individ- ual (HD 73526) 1. Introduction Since the discovery of the extra-solar planet orbiting 51 Peg (Mayor & Queloz 1995), high-precision radial-velocity measurements of the reflex radial velocity of the host star have lead to the discovery of over 130 planetary candidate companions to solar-type stars both in single and multiple star systems (e.g., see Butler et al. 2002, Eggenberger et al. 2003). Most of these planets have M sin i comparable to or greater than Jupiter, where M is the mass of the planet and i is the angle between the orbital plane and a plane perpendicular to the line of sight to the planetary system. However, as Doppler precision increases, lower This is an unedited preprint of an article accepted for publication in The Astrophysical Journal. The final published article may differ from this preprint. Copyright 2005 by The American Astronomical Society. Please cite as 'ApJ preprint doi:10.1086/'432594''. This is an unedited preprint of an article accepted for publication in The Astrophysical Journal. The final published article may differ from this preprint. Copyright 2005 by The American Astronomical Society. Please cite as 'ApJ preprint doi:10.1086/'432594''. –2– mass planets to the level of Neptune are also being discovered (McArthur et al. 2004, Fisher et al. 2003). Existing orbital solutions are based on the use of nonlinear-least squares methods which typically require a good initial guess of the parameter values. Frequently, the first indication of a periodic signal comes from a periodogram or epoch folding analysis of the data. These periodogram methods do not make use of relevant prior information about the known shape of the periodic signal. In the extrasolar planet Kepler problem, the mathematical form of the periodic signal is well known and is readily built into a Bayesian analysis. In this work, we re-analyze recently published high-precision radial-velocity measurements using a Bayesian Markov Chain Monte Carlo (MCMC) algorithm that is effective for both detecting and characterizing the orbits of extrasolar planets. This approach has the following advantages: a) Good initial guesses of the parameter values are not required. A MCMC is capable of efficiently exploring all regions of joint prior parameter space having significant probability 1. There is thus no need to carry out a separate search for the orbital period(s). b) The method is also capable of fitting a portion of an orbital period, so search periods longer than the duration of the data are possible. c) The built-in Occam’s razor in the Bayesian analysis can save considerable time in deciding whether a detection is believable. More complicated models contain larger numbers of parameters and thus incur a larger Occam penalty, which is automatically incorporated in a Bayesian model selection analysis in a quantitative fashion. The analysis yields the relative probability of each of the models explored. d) The analysis yields the full marginal posterior probability density function (PDF) for each model parameter, not just the maximum a posterior (MAP) values and a Gaussian approximation of their uncertainties. In 2003 Tinney et al. reported the detection of four new planets orbiting HD 2039, 30177, 73526, and 76700. Our results for HD 3029, together with some of the details of the Bayesian MCMC approach, were presented in Chapter 12 of Gregory (2005). In this paper we describe a modified form of the Bayesian MCMC approach and report the results of our re-analysis of the HD 73526 data which differ significantly from the conclusions of Tinney et al. (2003). For another application of MCMC methods to extrasolar planet detection see Ford (2004). 1It may fail to detect pathologically narrow probability peaks in a reasonable number of iterations. This is an unedited preprint of an article accepted for publication in The Astrophysical Journal. The final published article may differ from this preprint. Copyright 2005 by The American Astronomical Society. Please cite as 'ApJ preprint doi:10.1086/'432594''. –3– 2. Planetary Orbital Models There are two questions of interest. Initially, we are interested in whether the reflex motion of the host star provides evidence for 0, 1, 2, 3, ···planets. Secondly, we are interested in the orbital parameters of each planet detected. 2.1. Model Space The starting point for any Bayesian analysis is Bayes’ theorem given by equation (1). p(M|I)p(D|M,I) p(M |D, I)= i i (1) i p(D|I) where p(Mi|D, I) is the posterior probability that model Mi is true, given the new data D and prior background information I. The logical conjunction D, I is a proposition representing our current state of knowledge that is relevant for computing the probability of the truth of the models in question. A Bayesian analysis is an encoding of our current state of knowledge into a probability distribution concerning the hypothesis space of interest. Of course, as our state of knowledge changes (e.g., we acquire new data) so in general will the probabilities. The posterior is proportional to p(Mi|I), the prior probability that Mi is true, times the likelihood function, p(D|Mi,I). The likelihood is the probability that we would have gotten the data D that we did, if Mi is true. The likelihood is often written as L(Mi). In a well posed problem the prior (background) information (I) specifies the hypothesis space of interest (model space), and how to compute the likelihood function. In general, the models have parameters that are discussed below. Table 1 summarizes the hypothesis space of interest. For multiple planet models, there is no analytic expression for the exact radial velocity perturbation 2. In many cases, the radial velocity perturbation can be well modeled as the sum of multiple independent Keplerian orbits which is what has been assumed in this paper. 2The deviations from the simple sum of Keplerian models can be divided into two types (Brown 2005): short-period interactions and secular interactions. The magnitude of the short-period interactions is often small compared to the magnitude of the Keplerian perturbation and the observational uncertainties. The secular interactions are typically modeled as changes in the Keplerian orbital parameters. While the secular interactions can have large effects on the observed radial velocities, the timescales are typically much longer than the time span of observations. This is an unedited preprint of an article accepted for publication in The Astrophysical Journal. The final published article may differ from this preprint. Copyright 2005 by The American Astronomical Society. Please cite as 'ApJ preprint doi:10.1086/'432594''. –4– Table 1. Model space Symbol Model Number of parameters M0 constant velocity V 1 M0s constant velocity V + extra noise term 1 + s M1 V + elliptical orbit 6 M1s V + elliptical orbit + extra noise term 6 + s M2 V + 2 elliptical orbits 11 M2s V + 2 elliptical orbits + extra noise term 11 + s M3 V + 3 elliptical orbits 16 M3s V + 3 elliptical orbits + extra noise term 16 + s etc. This is an unedited preprint of an article accepted for publication in The Astrophysical Journal. The final published article may differ from this preprint. Copyright 2005 by The American Astronomical Society. Please cite as 'ApJ preprint doi:10.1086/'432594''. –5– In the case of model M0, the star is moving away or towards the observer with an unknown velocity. Note, the radial velocities published by Tinney et al. (2003) have an arbitrary zero point determined by a template observation, so the V inouranalysisisjust a constant velocity term. Model M0s assumes that there is an extra noise term in addition to the known measurement uncertainties. The nature of this extra noise term is discussed below. The simplest single planet model would have a circular orbit with four unknown parameters. Although the circular model is simpler, we will focus only on elliptical orbits because the prior probability of an orbit with eccentricity of precisely zero is considered to be extremely small. We illustrate the calculations for models M0, M0s, M1,andM1s. Much of the background material needed to appreciate the current analysis is presented in Chapter 12 of Gregory (2005). We can represent the measured velocities by the equation vi = fi + ei + e0, (2) where ei is the noise component arising from known but unequal measurement errors with an assumed Gaussian distribution, The term e0 accounts for any additional unknown mea- surement errors 3 plus any real signal in the data that cannot be explained by the model prediction fi.
Recommended publications
  • Naming the Extrasolar Planets
    Naming the extrasolar planets W. Lyra Max Planck Institute for Astronomy, K¨onigstuhl 17, 69177, Heidelberg, Germany [email protected] Abstract and OGLE-TR-182 b, which does not help educators convey the message that these planets are quite similar to Jupiter. Extrasolar planets are not named and are referred to only In stark contrast, the sentence“planet Apollo is a gas giant by their assigned scientific designation. The reason given like Jupiter” is heavily - yet invisibly - coated with Coper- by the IAU to not name the planets is that it is consid- nicanism. ered impractical as planets are expected to be common. I One reason given by the IAU for not considering naming advance some reasons as to why this logic is flawed, and sug- the extrasolar planets is that it is a task deemed impractical. gest names for the 403 extrasolar planet candidates known One source is quoted as having said “if planets are found to as of Oct 2009. The names follow a scheme of association occur very frequently in the Universe, a system of individual with the constellation that the host star pertains to, and names for planets might well rapidly be found equally im- therefore are mostly drawn from Roman-Greek mythology. practicable as it is for stars, as planet discoveries progress.” Other mythologies may also be used given that a suitable 1. This leads to a second argument. It is indeed impractical association is established. to name all stars. But some stars are named nonetheless. In fact, all other classes of astronomical bodies are named.
    [Show full text]
  • A Detailed Analysis of the Hd 73526 2:1 Resonant Planetary System∗
    The Astrophysical Journal, 780:140 (9pp), 2014 January 10 doi:10.1088/0004-637X/780/2/140 C 2014. The American Astronomical Society. All rights reserved. Printed in the U.S.A. A DETAILED ANALYSIS OF THE HD 73526 2:1 RESONANT PLANETARY SYSTEM∗ Robert A. Wittenmyer1,2,XianyuTan3,4, Man Hoi Lee3,5, Jonathan Horner1,2,C.G.Tinney1,2, R. P. Butler6, G. S. Salter1,2, B. D. Carter7,H.R.A.Jones8, S. J. O’Toole9, J. Bailey1,2, D. Wright1,2,J.D.Crane10, S. A. Schectman10, P. Arriagada6, I. Thompson10, D. Minniti11,12,andM.Diaz11 1 School of Physics, University of New South Wales, Sydney 2052, Australia; [email protected] 2 Australian Centre for Astrobiology, University of New South Wales, Sydney 2052, Australia 3 Department of Earth Sciences, The University of Hong Kong, Pokfulam Road, Hong Kong 4 Department of Planetary Sciences and Lunar and Planetary Laboratory, The University of Arizona, 1629 University Boulevard, Tucson, AZ 85721, USA 5 Department of Physics, The University of Hong Kong, Pokfulam Road, Hong Kong 6 Department of Terrestrial Magnetism, Carnegie Institution of Washington, 5241 Broad Branch Road, NW, Washington, DC 20015-1305, USA 7 Faculty of Sciences, University of Southern Queensland, Toowoomba, Queensland 4350, Australia 8 University of Hertfordshire, Centre for Astrophysics Research, Science and Technology Research Institute, College Lane, AL10 9AB, Hatfield, UK 9 Australian Astronomical Observatory, PO Box 915, North Ryde, NSW 1670, Australia 10 The Observatories of the Carnegie Institution of Washington, 813 Santa Barbara Street, Pasadena, CA 91101, USA 11 Institute of Astrophysics, Pontificia Universidad Catolica de Chile, Casilla 306, Santiago 22, Chile 12 Vatican Observatory, V00120 Vatican City State, Italy Received 2013 August 14; accepted 2013 November 26; published 2013 December 18 ABSTRACT We present six years of new radial velocity data from the Anglo-Australian and Magellan Telescopes on the HD 73526 2:1 resonant planetary system.
    [Show full text]
  • 2016 Publication Year 2021-04-23T14:32:39Z Acceptance in OA@INAF Age Consistency Between Exoplanet Hosts and Field Stars Title B
    Publication Year 2016 Acceptance in OA@INAF 2021-04-23T14:32:39Z Title Age consistency between exoplanet hosts and field stars Authors Bonfanti, A.; Ortolani, S.; NASCIMBENI, VALERIO DOI 10.1051/0004-6361/201527297 Handle http://hdl.handle.net/20.500.12386/30887 Journal ASTRONOMY & ASTROPHYSICS Number 585 A&A 585, A5 (2016) Astronomy DOI: 10.1051/0004-6361/201527297 & c ESO 2015 Astrophysics Age consistency between exoplanet hosts and field stars A. Bonfanti1;2, S. Ortolani1;2, and V. Nascimbeni2 1 Dipartimento di Fisica e Astronomia, Università degli Studi di Padova, Vicolo dell’Osservatorio 3, 35122 Padova, Italy e-mail: [email protected] 2 Osservatorio Astronomico di Padova, INAF, Vicolo dell’Osservatorio 5, 35122 Padova, Italy Received 2 September 2015 / Accepted 3 November 2015 ABSTRACT Context. Transiting planets around stars are discovered mostly through photometric surveys. Unlike radial velocity surveys, photo- metric surveys do not tend to target slow rotators, inactive or metal-rich stars. Nevertheless, we suspect that observational biases could also impact transiting-planet hosts. Aims. This paper aims to evaluate how selection effects reflect on the evolutionary stage of both a limited sample of transiting-planet host stars (TPH) and a wider sample of planet-hosting stars detected through radial velocity analysis. Then, thanks to uniform deriva- tion of stellar ages, a homogeneous comparison between exoplanet hosts and field star age distributions is developed. Methods. Stellar parameters have been computed through our custom-developed isochrone placement algorithm, according to Padova evolutionary models. The notable aspects of our algorithm include the treatment of element diffusion, activity checks in terms of 0 log RHK and v sin i, and the evaluation of the stellar evolutionary speed in the Hertzsprung-Russel diagram in order to better constrain age.
    [Show full text]
  • AMD-Stability and the Classification of Planetary Systems
    A&A 605, A72 (2017) DOI: 10.1051/0004-6361/201630022 Astronomy c ESO 2017 Astrophysics& AMD-stability and the classification of planetary systems? J. Laskar and A. C. Petit ASD/IMCCE, CNRS-UMR 8028, Observatoire de Paris, PSL, UPMC, 77 Avenue Denfert-Rochereau, 75014 Paris, France e-mail: [email protected] Received 7 November 2016 / Accepted 23 January 2017 ABSTRACT We present here in full detail the evolution of the angular momentum deficit (AMD) during collisions as it was described in Laskar (2000, Phys. Rev. Lett., 84, 3240). Since then, the AMD has been revealed to be a key parameter for the understanding of the outcome of planetary formation models. We define here the AMD-stability criterion that can be easily verified on a newly discovered planetary system. We show how AMD-stability can be used to establish a classification of the multiplanet systems in order to exhibit the planetary systems that are long-term stable because they are AMD-stable, and those that are AMD-unstable which then require some additional dynamical studies to conclude on their stability. The AMD-stability classification is applied to the 131 multiplanet systems from The Extrasolar Planet Encyclopaedia database for which the orbital elements are sufficiently well known. Key words. chaos – celestial mechanics – planets and satellites: dynamical evolution and stability – planets and satellites: formation – planets and satellites: general 1. Introduction motion resonances (MMR, Wisdom 1980; Deck et al. 2013; Ramos et al. 2015) could justify the Hill-type criteria, but the The increasing number of planetary systems has made it nec- results on the overlap of the MMR island are valid only for close essary to search for a possible classification of these planetary orbits and for short-term stability.
    [Show full text]
  • A Bayesian Analysis of Extrasolar Planet Data for HD 73526
    Astrophysical Journal, 631, 1198, 2005 A Bayesian Analysis of Extrasolar Planet Data for HD 73526 P. C. Gregory Physics and Astronomy Department, University of British Columbia, Vancouver, BC V6T 1Z1 ABSTRACT A new Bayesian tool for nonlinear model fitting has been developed, that employs a parallel tempering Markov Chain Monte Carlo algorithm with a novel statistical control system. The algorithm has been used to re-analyze precision radial velocity data for HD 73526. For a single planet model, three possible +0.37 +1.8 +1.4 orbits were found with periods of 127.88−0.09, 190.4−2.1 and 376.2−4.3 days. The 128 day orbit, with an eccentricity of 0.56 ± 0.03, has a maximum value of prior × likelihood which is 16 times larger than for the next highest solution at 376 +0.05 days. However, the 376 day orbit, with an eccentricity of 0.10−0.10, is formally more probable because for this sparse data set there is a much larger volume of parameter space with a significant probability density in the vicinity of the 376 day peak. The previously reported orbit (Tinney et al. 2003) of 190.5 ± 3.0 days corresponds to our least probable orbit. The analysis highlights the need for measurements around phase 0.5 for the 376 day period. Subject headings: Bayesian methods, time series analysis, extrasolar planets: individual( HD 73526) 1. Introduction Since the discovery of the extra-solar planet orbiting 51 Peg (Mayor & Queloz 1995), high-precision radial-velocity measurements of the reflex radial velocity of the host star have lead to the discovery of over 130 planetary candidate companions to solar-type stars both in single and multiple star systems (e.g., Eggenberger et al.
    [Show full text]
  • Emerging Trends in Metallicity and Lithium Properties of Debris Disc Stars
    MNRAS 000, 1–19 (2019) Preprint 30 May 2019 Compiled using MNRAS LATEX style file v3.0 Emerging trends in metallicity and lithium properties of debris disc stars ⋆ C. Chavero,1† R. de la Reza,2 L. Ghezzi,2,6 F. Llorente de Andr´es,3 C. B. Pereira,2 C. Giuppone,4 and G. Pinz´on5 1 Universidad Nacional de C´ordoba, Observatorio Astron´omico, Laprida 854, 5000 C´ordoba, CONICET, Argentina 2Observat´orio Nacional, Rua General Jos´eCristino, 77, 20921-400, S˜ao Crist´ov˜ao, Rio de Janeiro, RJ, Brazil 3 Ateneo de Almagro (Secc. Ciencia & Tecnolog´ıa)- 13270 Almagro, Spain 4 Universidad Nacional de C´ordoba, Observatorio Astron´omico, IATE, Laprida 854, 5000 C´ordoba, Argentina 5 Universidad Nacional de Colombia, Colombia 6 Observat´orio do Valongo, Universidade Federal do Rio de Janeiro, Ladeira do Pedro Antˆonio 43, Rio de Janeiro, RJ 20080-090 Accepted 2019 May 21; Received 2019 May 21; in original form 2018 October 4 ABSTRACT Dwarf stars with debris discs and planets appear to be excellent laboratories to study the core accretion theory of planets formation. These systems are however, insuffi- ciently studied. In this paper we present the main metallicity and lithium abundance properties of these stars together with stars with only debris discs and stars with only planets. Stars without detected planets nor discs are also considered. The anal- ysed sample is formed by main-sequence FGK field single stars. Apart from the basic stellar parameters, we include the use of dusty discs masses. The main results show for the first time that the dust mass of debris disc stars with planets correlate with metallicity.
    [Show full text]
  • Statistical Properties of Exoplanets
    A&A 398, 363–376 (2003) Astronomy DOI: 10.1051/0004-6361:20021637 & c ESO 2003 Astrophysics Statistical properties of exoplanets II. Metallicity, orbital parameters, and space velocities N. C. Santos1, G. Israelian2, M. Mayor1,R.Rebolo2,3,andS.Udry1 1 Observatoire de Gen`eve, 51 ch. des Maillettes, 1290 Sauverny, Switzerland 2 Instituto de Astrof´ısica de Canarias, 38200 La Laguna, Tenerife, Spain 3 Consejo Superior de Investigaciones Cient´ıficas, Spain Received 19 September 2002 / Accepted 7 November 2002 Abstract. In this article we present a detailed spectroscopic analysis of more than 50 extra-solar planet host stars. Stellar atmospheric parameters and metallicities are derived using high resolution and high S/N spectra. The spectroscopy results, added to the previous studies, imply that we have access to a large and uniform sample of metallicities for about 80 planet hosts stars. We make use of this sample to confirm the metal-rich nature of stars with planets, and to show that the planetary frequency is rising as a function of the [Fe/H]. Furthermore, the source of this high metallicity is shown to have most probably a “primordial” source, confirming previous results. The comparison of the orbital properties (period and eccentricity) and minimum masses of the planets with the stellar properties also reveal some emerging but still not significant trends. These are discussed and some explanations are proposed. Finally, we show that the planet host stars included in the CORALIE survey have similar kinematical properties as the whole CORALIE volume-limited planet search sample. Planet hosts simply seem to occupy the metal-rich envelope of this latter population.
    [Show full text]
  • Stability and Formation of the Resonant System HD 73526
    A&A 472, 981–992 (2007) Astronomy DOI: 10.1051/0004-6361:20077345 & c ESO 2007 Astrophysics Stability and formation of the resonant system HD 73526 Zs. Sándor1,2,W.Kley2, and P. Klagyivik1 1 Department of Astronomy, Eötvös Loránd University, Pázmány Péter sétány 1/A, 1117 Budapest, Hungary e-mail: [Zs.Sandor;P.Klagyivik]@astro.elte.hu 2 Institut für Astronomie und Astrophysik, Abt. Computational Physics, Universität Tübingen, Auf der Morgenstelle 10, 72076 Tübingen, Germany e-mail: [email protected] Received 23 February 2007 / Accepted 4 June 2007 ABSTRACT Context. Based on radial velocity measurements, it has been found that the two giant planets detected around the star HD 73526 are in 2:1 resonance. However, as our numerical integration shows, the derived orbital data for this system result in chaotic behavior of the giant planets, which is uncommon among the resonant extrasolar planetary systems. Aims. We present regular (non-chaotic) orbital solutions for the giant planets in the system HD 73526 and offer formation scenarios based on combining planetary migration and sudden perturbative effects such as planet-planet scattering or rapid dispersal of the protoplanetary disk. A comparison with the already-studied resonant system HD 128311, exhibiting similar behavior, is given. Methods. The new sets of orbital solutions were derived using the Systemic Console. The stability of these solutions was investigated using the Relative Lyapunov indicator, while the migration and scattering effects are studied by gravitational N-body simulations ap- plying non-conservative forces. Additionally, hydrodynamic simulations of embedded planets in protoplanetary disks were performed to follow the capture into resonance.
    [Show full text]
  • Quantization of Planetary Systems and Its Dependency on Stellar Rotation Jean-Paul A
    Quantization of Planetary Systems and its Dependency on Stellar Rotation Jean-Paul A. Zoghbi∗ ABSTRACT With the discovery of now more than 500 exoplanets, we present a statistical analysis of the planetary orbital periods and their relationship to the rotation periods of their parent stars. We test whether the structure of planetary orbits, i.e. planetary angular momentum and orbital periods are ‘quantized’ in integer or half-integer multiples with respect to the parent stars’ rotation period. The Solar System is first shown to exhibit quantized planetary orbits that correlate with the Sun’s rotation period. The analysis is then expanded over 443 exoplanets to statistically validate this quantization and its association with stellar rotation. The results imply that the exoplanetary orbital periods are highly correlated with the parent star’s rotation periods and follow a discrete half-integer relationship with orbital ranks n=0.5, 1.0, 1.5, 2.0, 2.5, etc. The probability of obtaining these results by pure chance is p<0.024. We discuss various mechanisms that could justify this planetary quantization, such as the hybrid gravitational instability models of planet formation, along with possible physical mechanisms such as inner discs magnetospheric truncation, tidal dissipation, and resonance trapping. In conclusion, we statistically demonstrate that a quantized orbital structure should emerge naturally from the formation processes of planetary systems and that this orbital quantization is highly dependent on the parent stars rotation periods. Key words: planetary systems: formation – star: rotation – solar system: formation 1. INTRODUCTION The discovery of now more than 500 exoplanets has provided the opportunity to study the various properties of planetary systems and has considerably advanced our understanding of planetary formation processes.
    [Show full text]
  • Radio Pulsars in Binary Systems
    Radio Pulsars in Binary Systems Ren´ePaul Breton Department of Physics McGill University Montr´eal,Qu´ebec Canada December 2008 A Thesis submitted to McGill University in partial fulfillment of the requirements of the degree of Doctor of Philosophy c Ren´eP. Breton, 2008 Pour mon p`ere,ma m`ereet mon fr`ere. Votre support depuis toujours est un cadeau inestimable. Contents Abstract xi Resum´ e´ xii Acknowledgments xiii Contributions of Authors xvi 1 Introduction 1 1.1 In a Nutshell... 1 1.2 Historical Background . 4 1.2.1 The Great Discovery . 4 1.2.2 Theoretical Precursor Work . 11 1.2.3 Parallel Discoveries . 13 1.3 Thesis Outline . 14 2 Neutron Stars and Pulsars 17 2.1 Neutron Stars: An Overview . 17 2.1.1 Internal Structure and Equation of State . 20 2.1.2 Magnetic Field . 24 2.2 Pulsars: An Overview . 25 2.2.1 Rotation-Powered Pulsars . 26 2.2.2 Non-Rotation-Powered Pulsars . 35 2.3 Measurements of Pulsar Properties . 37 2.3.1 Mass . 37 2.3.2 Radius . 43 2.3.3 Magnetic Field . 44 2.4 Radio Telescopes . 47 2.5 Pulsar Timing . 51 2.5.1 Extrinsic Timing Components . 54 2.5.2 Intrinsic Timing Components . 58 2.5.3 Binary Pulsar Timing . 60 2.6 Binary Pulsars . 64 2.6.1 Binary Evolution . 66 2.6.2 Binary Radio Pulsar Population . 73 2.6.3 Other Pulsar Binaries . 80 2.7 Binary Pulsars as Benchmarks for Gravity Theories . 86 2.7.1 “Classical” Tests of Gravity .
    [Show full text]
  • Stability and Formation of the Resonant System HD 73526 Oscillations of the Eccentricities and the Resonant Angles
    Astronomy & Astrophysics manuscript no. sandoretal c ESO 2018 June 1, 2018 Stability and Formation of the Resonant System HD 73526 Zs. S´andor,1,2 W. Kley,2 and P. Klagyivik1 1 Department of Astronomy, E¨otv¨os Lor´and University, P´azm´any P´eter s´et´any 1/A, H–1117 Budapest, Hungary e-mail: [email protected] e-mail: [email protected] 2 Institut f¨ur Astronomie und Astrophysik, Abt. Computational Physics, Universit¨at T¨ubingen, Auf der Morgenstelle 10, D–72076 T¨ubingen, Germany e-mail: [email protected] Received ; accepted ABSTRACT Context. Based on radial velocity measurements it has been found recently that the two giant planets detected around the star HD 73526 are in 2:1 resonance. However, as our numerical integration shows, the derived orbital data for this system result in chaotic behavior of the giant planets, which is uncommon among the resonant extrasolar planetary systems. Aims. We intend to present regular (non-chaotic) orbital solutions for the giant planets in the system HD 73526 and offer formation scenarios based on combining planetary migration and sudden perturbative effects such as planet-planet scattering or rapid dispersal of the protoplanetary disk. A comparison with the already studied resonant system HD 128311, exhibiting similar behavior, is also done. Methods. The new sets of orbital solutions have been derived by the Systemic Console (http://www.oklo.org). The stability of these solutions has been investigated by the Relative Lyapunov indicator, while the migration and scattering effects are studied by gravitational N-body simulations applying non-conservative forces as well.
    [Show full text]
  • Solar System Analogues Among Exoplanetary Systems
    Solar System analogues among exoplanetary systems Maria Lomaeva Lund Observatory Lund University ´´ 2016-EXA105 Degree project of 15 higher education credits June 2016 Supervisor: Piero Ranalli Lund Observatory Box 43 SE-221 00 Lund Sweden Populärvetenskaplig sammanfattning Människans intresse för rymden har alltid varit stort. Man har antagit att andra plan- etsystem, om de existerar, ser ut som vårt: med mindre stenplaneter i banor närmast stjärnan och gas- samt isjättar i de yttre banorna. Idag känner man till drygt 2 000 exoplaneter, d.v.s., planeter som kretsar kring andra stjärnor än solen. Man vet även att vissa av dem saknar motsvarighet i solsystemet, t. ex., heta jupitrar (gasjättar som har migrerat inåt och kretsar väldigt nära stjärnan) och superjordar (stenplaneter större än jorden). Därför blir frågan om hur unikt solsystemet är ännu mer intressant, vilket vi försöker ta reda på i det här projektet. Det finns olika sätt att detektera exoplaneter på men två av dem har gett flest resultat: transitmetoden och dopplerspektroskopin. Med transitmetoden mäter man minsknin- gen av en stjärnas ljus när en planet passerar framför den. Den metoden passar bäst för stora planeter med små omloppsbanor. Dopplerspektroskopin använder sig av Doppler effekten som innebär att ljuset utsänt från en stjärna verkar blåare respektive rödare när en stjärna förflyttar sig fram och tillbaka från observatören. Denna rörelse avslöjar att det finns en planet som kretsar kring stjärnan och påverkar den med sin gravita- tion. Dopplerspektroskopin är lämpligast för massiva planeter med små omloppsbanor. Under projektets gång har vi inte bara letat efter solsystemets motsvarigheter utan även studerat planetsystem som är annorlunda.
    [Show full text]