N-Body Dynamics of Intermediate Mass-Ratio Inspirals in Star Clusters Haster, Carl-Johan; Antonini, Fabio; Kalogera, Vicky; Mandel, Ilya

Total Page:16

File Type:pdf, Size:1020Kb

N-Body Dynamics of Intermediate Mass-Ratio Inspirals in Star Clusters Haster, Carl-Johan; Antonini, Fabio; Kalogera, Vicky; Mandel, Ilya University of Birmingham n-body dynamics of intermediate mass-ratio inspirals in star clusters Haster, Carl-Johan; Antonini, Fabio; Kalogera, Vicky; Mandel, Ilya DOI: 10.3847/0004-637X/832/2/192 License: None: All rights reserved Document Version Peer reviewed version Citation for published version (Harvard): Haster, C-J, Antonini, F, Kalogera, V & Mandel, I 2016, 'n-body dynamics of intermediate mass-ratio inspirals in star clusters', The Astrophysical Journal, vol. 832, no. 2. https://doi.org/10.3847/0004-637X/832/2/192 Link to publication on Research at Birmingham portal General rights Unless a licence is specified above, all rights (including copyright and moral rights) in this document are retained by the authors and/or the copyright holders. The express permission of the copyright holder must be obtained for any use of this material other than for purposes permitted by law. •Users may freely distribute the URL that is used to identify this publication. •Users may download and/or print one copy of the publication from the University of Birmingham research portal for the purpose of private study or non-commercial research. •User may use extracts from the document in line with the concept of ‘fair dealing’ under the Copyright, Designs and Patents Act 1988 (?) •Users may not further distribute the material nor use it for the purposes of commercial gain. Where a licence is displayed above, please note the terms and conditions of the licence govern your use of this document. When citing, please reference the published version. Take down policy While the University of Birmingham exercises care and attention in making items available there are rare occasions when an item has been uploaded in error or has been deemed to be commercially or otherwise sensitive. If you believe that this is the case for this document, please contact [email protected] providing details and we will remove access to the work immediately and investigate. Download date: 02. Oct. 2021 Draft version June 24, 2016 Preprint typeset using LATEX style emulateapj v. 01/23/15 N−BODY DYNAMICS OF INTERMEDIATE MASS-RATIO INSPIRALS IN GLOBULAR CLUSTERS Carl-Johan Haster School of Physics and Astronomy, University of Birmingham, Edgbaston, Birmingham B15 2TT, United Kingdom and Center for Interdisciplinary Exploration and Research in Astrophysics (CIERA) & Dept. of Physics and Astronomy, 2145 Sheridan Rd, Evanston, IL 60208, USA Fabio Antonini Center for Interdisciplinary Exploration and Research in Astrophysics (CIERA) & Dept. of Physics and Astronomy, 2145 Sheridan Rd, Evanston, IL 60208, USA Vicky Kalogera Center for Interdisciplinary Exploration and Research in Astrophysics (CIERA) & Dept. of Physics and Astronomy, 2145 Sheridan Rd, Evanston, IL 60208, USA Ilya Mandel School of Physics and Astronomy, University of Birmingham, Edgbaston, Birmingham B15 2TT, United Kingdom Draft version June 24, 2016 ABSTRACT The intermediate mass-ratio inspiral of a stellar compact remnant into an intermediate mass black hole (IMBH) can produce a gravitational wave (GW) signal that is potentially detectable by current ground-based GW detectors (e.g., Advanced LIGO) as well as by planned space-based interferome- ters (e.g., eLISA). Here, we present results from a direct integration of the post-Newtonian N-body equations of motion describing stellar clusters containing an IMBH and a population of stellar-mass black holes (BHs) and solar mass stars. We take particular care to simulate the dynamics closest to the IMBH, including post-Newtonian effects up to order 2.5. Our simulations show that the IMBH readily forms a binary with a BH companion. This binary is gradually hardened by transient 3-body or 4-body encounters, leading to frequent substitutions of the BH companion, while the binary’s eccen- tricity experiences large amplitude oscillations due to the Lidov-Kozai resonance. We also demonstrate suppression of these resonances by the relativistic precession of the binary orbit. We find an inter- mediate mass-ratio inspiral in one of the 12 cluster models we evolved for ∼ 100 Myr. This cluster hosts a 100M IMBH embedded in a population of 32 10M BH and 32,000 1M stars. At the end of the simulation, after ∼ 100 Myr of evolution, the IMBH merges with a BH companion. The IMBH–BH binary inspiral starts in the eLISA frequency window (& 1mHz) when the binary reaches an eccentricity 1 − e ' 10−3. After ' 105 years the binary moves into the LIGO frequency band with a negligible eccentricity. We comment on the implications for GW searches, with a possible detection within the next decade. 1. INTRODUCTION rell et al. 2009; Davis et al. 2011; Godet et al. 2014), Intermediate mass black holes (IMBHs) are conjec- but these dynamical measurements alone can not pro- tured to occupy the mass range between stellar-mass vide conclusive proof for the existence of IMBHs. If these lower-mass IMBHs reside in globular clusters, black holes (BHs), with masses . 100M , and super- 6 they will play an important role in cluster dynamics (e.g. massive black holes with masses & 10 M (see Miller & Colbert 2004, for a review). While the existence Trenti 2006; Umbreit & Rasio 2013; Konstantinidis et al. arXiv:1606.07097v1 [astro-ph.HE] 22 Jun 2016 of some IMBH candidates in dwarf spheroidal galaxies 2013; Leigh et al. 2014; MacLeod et al. 2016). Of particu- has been conjectured by extending the M–σ relation lar interest to our study is the likely tendency of IMBHs (Graham & Scott 2013) (but see Maccarone & Servillat to dynamically form compact binaries with other com- 2008), dynamical measurements of IMBHs in the few- pact remnants (e.g. Taniguchi et al. 2000; Miller 2002; hundred solar-mass range are extremely challenging (e.g., Miller & Hamilton 2002; Amaro-Seoane & Freitag 2006; Pasquato et al. 2016). The best evidence for such lower Brown et al. 2007; Mandel et al. 2008; Amaro-Seoane & Santamaría 2010; Mapelli et al. 2010; Mapelli 2016). mass IMBHs (with mass ∼ 100M ) could come from ul- traluminous X-ray sources (but see Berghea et al. 2008); Generally, these analyses find that the IMBH readily cap- for example, (Pasham et al. 2014) have claimed a mass of tures a binary companion. The binary is subsequently hardened through a sequence of 3-body and 4-body in- ∼ 400 M for M82 X-1 from quasi-periodic oscillations, 4 teractions, occasionally with substitutions which make a while a mass around 10 M has been suggested for the brightest ultraluminous X-ray source HLX-1 (e.g., Far- black hole (BH) of a few tens of solar masses the most likely IMBH companion, and possible Lidov-Kozai (LK) resonances (Lidov 1962; Kozai 1962) if hierarchical triples [email protected] 2 are formed. Eventually, the IMBH–BH binary merges The simulations were performed with the number of par- through the radiation of gravitational waves, emitting ticles N ⊂ {32768, 65536} and the mass of the IMBH, a signal that is potentially detectable by the Advanced M ⊂ {50, 100, 200} M and a cluster virial radius rv of LIGO ground-based GW detectors (Aasi et al. 2015; 3.5 pc. For N = 32768 and all three IMBH masses, sim- Abadie et al. 2010; Smith et al. 2013; Haster et al. 2016). ulations with rv ⊂ {0.35, 1.0} pc were also performed. Previous simulations of globular clusters with IMBH The inclusion of high-order pN terms fixes the physical coalescences have generally simplified the interactions in scale of the cluster, thus removing the conventional free- order to avoid excessive computational cost. For exam- dom for rescaling simulations in cluster size and density. ple, Gültekin et al.(2004) considered a series of individ- We observe the IMBH forming a binary with a BH ual Newtonian interactions interspersed with orbital evo- within . 20Myr in every simulated cluster. Only in one lution through GW emission. Mandel et al.(2008) car- cluster (N = 32768,M = 100M , rv = 3.5pc) we observe ried out analytical estimates of the hardening sequence a merger within the simulated time (' 100Myr). While to obtain the intermediate mass-ratio merger timescale. the result of the entire set of simulations will be pre- Leigh et al.(2014) simulated the entire cluster with a sented in a future paper, in what follows we will focus on mixture of analytical and numerical N-body analytical describing the detailed dynamics of the one cluster pro- calculations, while MacLeod et al.(2016) focused their ducing the merger event. We note in passing that as the N-body investigation on tidal disruptions of stars by the main focus of this study is the dynamical formation and IMBH as well as merger events. We note that in the evolution of binaries, and higher order N-tuples, with previous literature effects of pN terms are either not ac- the IMBH as the primary companion, all cluster parti- counted for (Leigh et al. 2014), or included only at the cles are solely characterized by their mass and no stellar 2.5pN level (Samsing et al. 2014; MacLeod et al. 2016). evolution is included in these simulations. In this paper we show a clear example in which lower order pN terms play a fundamental role in the dynamics. 3. RESULTS More specifically, an essential element that differs be- The simulated globular cluster was initialised with the tween the relativistic and non-relativistic dynamics turns IMBH at rest at the center while the remaining stars and out to be the 1pN precession of the periapsis. BHs follows a King model. Figure 1 shows the position We introduce our numerical method and the simula- of the IMBH and a subset of BH particles, and their tion setup in section 2. We describe our simulation re- subsequent movement within the cluster, relative to the sults in section 3.
Recommended publications
  • Rocket Equation
    Rocket Equation & Multistaging PROFESSOR CHRIS CHATWIN LECTURE FOR SATELLITE AND SPACE SYSTEMS MSC UNIVERSITY OF SUSSEX SCHOOL OF ENGINEERING & INFORMATICS 25TH APRIL 2017 Orbital Mechanics & the Escape Velocity The motion of a space craft is that of a body with a certain momentum in a gravitational field The spacecraft moves under the combines effects of its momentum and the gravitational attraction towards the centre of the earth For a circular orbit V = wr (1) F = mrw2 = centripetal acceleration (2) And this is balanced by the gravitational attraction F = mMG / r2 (3) 퐺푀 1/2 퐺푀 1/2 푤 = 표푟 푉 = (4) 푟3 푟 Orbital Velocity V 퐺 = 6.67 × 10−11 N푚2푘푔2 Gravitational constant M = 5.97× 1024 kg Mass of the Earth 6 푟0 = 6.371× 10 푚 Hence V = 7900 m/s For a non-circular (eccentric) orbit 1 퐺푀푚2 = 1 + 휀 푐표푠휃 (5) 푟 ℎ2 ℎ = 푚푟푉 angular momentum (6) ℎ2 휀 = 2 − 1 eccentricity (7) 퐺푀푚 푟0 1 퐺푀 2 푉 = 1 + 휀 푐표푠휃 (8) 푟 Orbit type For ε = 0 circular orbit For ε = 1 parabolic orbit For ε > 1 hyperbolic orbit (interplanetary fly-by) Escape Velocity For ε = 1 a parabolic orbit, the orbit ceases to be closed and the space vehicle will not return In this case at r0 using equation (7) ℎ2 1 = 2 − 1 but ℎ = 푚푟푣 (9) 퐺푀푚 푟0 1 2 2 2 푚 푉0 푟0 2퐺푀 2 2 = 2 therefore 푉0 = (10) 퐺푀푚 푟0 푟0 퐺 = 6.67 × 10−11 N푚2푘푔2 Gravitational constant M = 5.975× 1024 kg Mass of the Earth 6 푟0 = 6.371× 10 푚 Mean earth radius 1 2퐺푀 2 푉0 = Escape velocity (11) 푟0 푉0= 11,185 m/s (~ 25,000 mph) 7,900 m/s to achieve a circular orbit (~18,000 mph) Newton’s 3rd Law & the Rocket Equation N3: “to every action there is an equal and opposite reaction” A rocket is a device that propels itself by emitting a jet of matter.
    [Show full text]
  • Rocket Basics 1: Thrust and Specific Impulse
    4. The Congreve Rocket ROCKET BASICS 1: THRUST AND SPECIFIC IMPULSE In a liquid-propellant rocket, a liquid fuel and oxidizer are injected into the com- bustion chamber where the chemical reaction of combustion occurs. The combustion process produces hot gases that expand through the nozzle, a duct with the varying cross section. , In solid rockets, both fuel and oxidizer are combined in a solid form, called the grain. The grain surface burns, producing hot gases that expand in the nozzle. The rate of the gas production and, correspondingly, pressure in the rocket is roughly proportional to the burning area of the grain. Thus, the introduction of a central bore in the grain allowed an increase in burning area (compared to “end burning”) and, consequently, resulted in an increase in chamber pressure and rocket thrust. Modern nozzles are usually converging-diverging (De Laval nozzle), with the flow accelerating to supersonic exhaust velocities in the diverging part. Early rockets had only a crudely formed converging part. If the pressure in the rocket is high enough, then the exhaust from a converging nozzle would occur at sonic velocity, that is, with the Mach number equal to unity. The rocket thrust T equals PPAeae TmUPPAmUeeae e mU eq m where m is the propellant mass flow (mass of the propellant leaving the rocket each second), Ue (exhaust velocity) is the velocity with which the propellant leaves the nozzle, Pe (exit pressure) is the propellant pressure at the nozzle exit, Pa (ambient 34 pressure) is the pressure in the surrounding air (about one atmosphere at sea level), and Ae is the area of the nozzle exit; Ueq is called the equivalent exhaust velocity.
    [Show full text]
  • MARYLAND Orbital Mechanics
    Orbital Mechanics • Orbital Mechanics, continued – Time in orbits – Velocity components in orbit – Deorbit maneuvers – Atmospheric density models – Orbital decay (introduction) • Fundamentals of Rocket Performance – The rocket equation – Mass ratio and performance – Structural and payload mass fractions – Multistaging – Optimal ΔV distribution between stages (introduction) U N I V E R S I T Y O F Orbital Mechanics © 2004 David L. Akin - All rights reserved Launch and Entry Vehicle Design MARYLAND http://spacecraft.ssl.umd.edu Calculating Time in Orbit U N I V E R S I T Y O F Orbital Mechanics MARYLAND Launch and Entry Vehicle Design Time in Orbit • Period of an orbit a3 P = 2π µ • Mean motion (average angular velocity) µ n = a3 • Time since pericenter passage M = nt = E − esin E ➥M=mean anomaly E=eccentric anomaly U N I V E R S I T Y O F Orbital Mechanics MARYLAND Launch and Entry Vehicle Design Dealing with the Eccentric Anomaly • Relationship to orbit r = a(1− e cosE) • Relationship to true anomaly θ 1+ e E tan = tan 2 1− e 2 • Calculating M from time interval: iterate Ei+1 = nt + esin Ei until it converges U N I V E R S I T Y O F Orbital Mechanics MARYLAND Launch and Entry Vehicle Design Example: Time in Orbit • Hohmann transfer from LEO to GEO – h1=300 km --> r1=6378+300=6678 km – r2=42240 km • Time of transit (1/2 orbital period) U N I V E R S I T Y O F Orbital Mechanics MARYLAND Launch and Entry Vehicle Design Example: Time-based Position Find the spacecraft position 3 hours after perigee E=0; 1.783; 2.494; 2.222; 2.361; 2.294; 2.328;
    [Show full text]
  • The Effects of Barycentric and Asymmetric Transverse Velocities on Eclipse and Transit Times
    The Astrophysical Journal, 854:163 (13pp), 2018 February 20 https://doi.org/10.3847/1538-4357/aaa3ea © 2018. The American Astronomical Society. All rights reserved. The Effects of Barycentric and Asymmetric Transverse Velocities on Eclipse and Transit Times Kyle E. Conroy1,2 , Andrej Prša2 , Martin Horvat2,3 , and Keivan G. Stassun1,4 1 Vanderbilt University, Departmentof Physics and Astronomy, 6301 Stevenson Center Lane, Nashville TN, 37235, USA 2 Villanova University, Department of Astrophysics and Planetary Sciences, 800 E. Lancaster Avenue, Villanova, PA 19085, USA 3 University of Ljubljana, Departmentof Physics, Jadranska 19, SI-1000 Ljubljana, Slovenia 4 Fisk University, Department of Physics, 1000 17th Avenue N., Nashville, TN 37208, USA Received 2017 October 20; revised 2017 December 19; accepted 2017 December 19; published 2018 February 21 Abstract It has long been recognized that the finite speed of light can affect the observed time of an event. For example, as a source moves radially toward or away from an observer, the path length and therefore the light travel time to the observer decreases or increases, causing the event to appear earlier or later than otherwise expected, respectively. This light travel time effect has been applied to transits and eclipses for a variety of purposes, including studies of eclipse timing variations and transit timing variations that reveal the presence of additional bodies in the system. Here we highlight another non-relativistic effect on eclipse or transit times arising from the finite speed of light— caused by an asymmetry in the transverse velocity of the two eclipsing objects, relative to the observer. This asymmetry can be due to a non-unity mass ratio or to the presence of external barycentric motion.
    [Show full text]
  • A Conceptual Analysis of Spacecraft Air Launch Methods
    A Conceptual Analysis of Spacecraft Air Launch Methods Rebecca A. Mitchell1 Department of Aerospace Engineering Sciences, University of Colorado, Boulder, CO 80303 Air launch spacecraft have numerous advantages over traditional vertical launch configurations. There are five categories of air launch configurations: captive on top, captive on bottom, towed, aerial refueled, and internally carried. Numerous vehicles have been designed within these five groups, although not all are feasible with current technology. An analysis of mass savings shows that air launch systems can significantly reduce required liftoff mass as compared to vertical launch systems. Nomenclature Δv = change in velocity (m/s) µ = gravitational parameter (km3/s2) CG = Center of Gravity CP = Center of Pressure 2 g0 = standard gravity (m/s ) h = altitude (m) Isp = specific impulse (s) ISS = International Space Station LEO = Low Earth Orbit mf = final vehicle mass (kg) mi = initial vehicle mass (kg) mprop = propellant mass (kg) MR = mass ratio NASA = National Aeronautics and Space Administration r = orbital radius (km) 1 M.S. Student in Bioastronautics, [email protected] 1 T/W = thrust-to-weight ratio v = velocity (m/s) vc = carrier aircraft velocity (m/s) I. Introduction T HE cost of launching into space is often measured by the change in velocity required to reach the destination orbit, known as delta-v or Δv. The change in velocity is related to the required propellant mass by the ideal rocket equation: 푚푖 훥푣 = 퐼푠푝 ∗ 0 ∗ ln ( ) (1) 푚푓 where Isp is the specific impulse, g0 is standard gravity, mi initial mass, and mf is final mass. Specific impulse, measured in seconds, is the amount of time that a unit weight of a propellant can produce a unit weight of thrust.
    [Show full text]
  • Revolutionary Propulsion Concepts for Small Satellites
    SSC01-IX-4 Revolutionary Propulsion Concepts for Small Satellites Steven R. Wassom, Ph.D., P.E. Senior Mechanical Engineer Space Dynamics Laboratory/Utah State University 1695 North Research Park Way North Logan, Utah 84341-1947 Phone: 435-797-4600 E-mail: [email protected] Abstract. This paper addresses the results of a trade study in which four novel propulsion approaches are applied to a 100-kg-class satellite designed for rendezvous, reconnaissance, and other on-orbit operations. The technologies, which are currently at a NASA technology readiness level of 4, are known as solar thermal propulsion, digital solid motor, water-based propulsion, and solid pulse motor. Sizing calculations are carried out using analytical and empirical parameters to determine the propellant and inert masses and volumes. The results are compared to an off- the-shelf hydrazine system using a trade matrix “scorecard.” Other factors considered besides mass and volume include safety, storability, mission time, accuracy, and refueling. Most of the concepts scored higher than the hydrazine system and warrant further development. The digital solid motor had the highest score by a small margin. Introduction propulsion concepts to a state-of-the-art (SOTA) hydrazine system, as applied to a notional 100-kg Aerospace America1 recently provided some interesting satellite designed for rendezvous with other objects. data on the number of payloads proposed for launch in the next 10 years. In the year 2000 there was a 68% Basic Description of Concepts increase in the number of proposed payloads weighing less than 100 kg, a 67% increase in the number of Solar Thermal Propulsion (STP) proposed military payloads, and a 92% increase in the number of proposed reconnaissance/surveillance Solar Thermal Propulsion (STP) uses the sun’s energy satellites.
    [Show full text]
  • Debris Assessment Software User’S Guide
    JSC 64047 Debris Assessment Software User’s Guide Version 2.0 Astromaterials Research and Exploration Science Directorate Orbital Debris Program Office November 2007 National Aeronautics and Space Administration Lyndon B. Johnson Space Center Houston, Texas 77058 Debris Assessment Software Version 2.0 User’s Guide NASA Lyndon B. Johnson Space Center Houston, TX November 2007 Prepared by: John N. Opiela Eric Hillary David O. Whitlock Marsha Hennigan ESCG Approved by: Kira Abercromby ESCG Project Manager Orbital Debris Program Office Approved by: Eugene Stansbery NASA/JSC Program Manager Orbital Debris Program Office Table of Contents List of Figures................................................................................................................................... iii List of Tables..................................................................................................................................... iii List of Plots ........................................................................................................................................iv 1. Introduction ................................................................................................................................1 1.1 Major Changes since DAS 1.5.3 ..........................................................................................1 1.2 Software Installation and Removal ......................................................................................2 2. DAS Main Window Features ....................................................................................................4
    [Show full text]
  • RWTH Aachen University Institute of Jet Propulsion and Turbomachinery
    RWTH Aachen University Institute of Jet Propulsion and Turbomachinery DLR - German Aerospace Center Institute of Space Propulsion Optimization of a Reusable Launch Vehicle using Genetic Algorithms A thesis submitted for the degree of M.Sc. Aerospace Engineering by Simon Jentzsch June 20, 2020 Advisor: Felix Wrede, M.Sc. External Advisors: Kai Dresia, M.Sc. Dr. Günther Waxenegger-Wilfing Examiner: Prof. Dr. Michael Oschwald Statutory Declaration in Lieu of an Oath I hereby declare in lieu of an oath that I have completed the present Master the- sis entitled ‘Optimization of a Reusable Launch Vehicle using Genetic Algorithms’ independently and without illegitimate assistance from third parties. I have used no other than the specified sources and aids. In case that the thesis is additionally submitted in an electronic format, I declare that the written and electronic versions are fully identical. The thesis has not been submitted to any examination body in this, or similar, form. City, Date, Signature Abstract SpaceX has demonstrated that reusing the first stage of a rocket implies a significant cost reduction potential. In order to maximize cost savings, the identification of optimum rocket configurations is of paramount importance. Yet, the complexity of launch systems, which is further increased by the requirement of a vertical landing reusable first stage, impedes the prediction of launch vehicle characteristics. Therefore, in this thesis, a multidisciplinary system design optimization approach is applied to develop an optimization platform which is able to model a reusable launch vehicle with a large variety of variables and to optimize it according to a predefined launch mission and optimization objective.
    [Show full text]
  • Lecture 7 Launch Trajectories
    AA 284a Advanced Rocket Propulsion Lecture 7 Launch Trajectories Prepared by Arif Karabeyoglu Department of Aeronautics and Astronautics Stanford University and Mechanical Engineering KOC University Fall 2019 Stanford University AA284a Advanced Rocket Propulsion Orbital Mechanics - Review • Newton’s law of gravitation: mMG Fg = r2 – M, m: Mass of the bodies – r: Distance between the center of masses of the two bodies – Fg: Gravitational attraction force between the two bodies – G: Universal gravitational constant • Assume that m is the mass of the spacecraft and M is the mass of the celestial body. Arrange the force expression as (Note that m << M) m 2 Fg = = g m = M G g = r r2 • Here the gravitational parameter has been introduced for convenience. It is a constant for a given celestial mass. For Earth 3 2 = 398,600 km /sec • For circular orbit: centrifugal force balancing the gravitational force acting on the satellite – Orbital Velocity: V = co r r3 – Orbital Period: P = 2 Stanford University 2 Karabeyoglu AA284a Advanced Rocket Propulsion Orbital Mechanics - Review • Fundamental Assumptions: – Two body assumption • Motion of the spacecraft is only affected my a single central body – The mass of the spacecraft is negligible compared to the mass of the celestial body – The bodies are spherically symmetric with the masses concentrated at the center of the sphere – No forces other than gravity (and inertial forces) r r = k r3 • Solution: • Sections of a cone Stanford University 3 Karabeyoglu AA284a Advanced Rocket Propulsion Orbital Mechanics - Review • Solution: – Orbits of any conic section, elliptic, parabolic, hyperbolic – Energy is conserved in the conservative force field.
    [Show full text]
  • A New Radial Velocity Curve, Apsidal Motion, and the Alignment of Rotation and Orbit Axes
    A&A 586, A104 (2016) Astronomy DOI: 10.1051/0004-6361/201526662 & c ESO 2016 Astrophysics The α CrB binary system: A new radial velocity curve, apsidal motion, and the alignment of rotation and orbit axes J. H. M. M. Schmitt1, K.-P. Schröder2,1,G.Rauw3, A. Hempelmann1, M. Mittag1, J. N. González-Pérez1, S. Czesla1,U.Wolter1, and D. Jack2 1 Hamburger Sternwarte, Universität Hamburg, 21029 Hambourg, Germany e-mail: [email protected] 2 Universidad de Guanajuato, Departamento de Astronomía, 36000 Guanajuato, Mexico 3 Institut d’Astrophysique et de Géophysique Université de Liège, 4000 Liège, Belgium Received 3 June 2015 / Accepted 12 November 2015 ABSTRACT We present a new radial velocity curve for the two components of the eclipsing spectroscopic binary α CrB. This binary consists of two main-sequence stars of types A and G in a 17.3599-day orbit, according to the data from our robotic TIGRE facility that is located in Guanajuato, Mexico. We used a high-resolution solar spectrum to determine the radial velocities of the weak secondary component by cross-correlation and wavelength referencing with telluric lines for the strongly rotationally broadened primary lines (v sin(i) = 138 km s−1) to obtain radial velocities with an accuracy of a few hundred m/s. We combined our new RV data with older measurements, dating back to 1908 in the case of the primary, to search for evidence of apsidal motion. We find an apsidal motion period between 6600 and 10 600 yr. This value is consistent with the available data for both the primary and secondary and is also consistent with the assumption that the system has aligned orbit and rotation axes.
    [Show full text]
  • AIAA 2001-3377 Solar Thermal Propulsion for an Interstellar Probe
    (c)2001 American Institute of Aeronautics & Astronautics or Published with Permission of Author(s) and/or Author(s)' Sponsoring Organization. AM A A A01-34143 AIAA 2001-3377 Solar Thermal Propulsion a r nfo Interstellar Probe R. W. Lyman, M. E. Ewing, R. S. Krishnan, D. M. Lester Thiokol Propulsion Brigham City, UT R. McNuttL . Jr , Applied Physics Laboratory Johns Hopkins University Laurel, MD 37th AIAA/ASME/SAE/ASEE Joint Propulsion Conference July 8-11,2001 Salt Lake City, Utah For permission to copy or to republish, contact the copyright owner named on the first page. For AlAA-held copyright, write to AIAA Permissions Department, 1801 Alexander Bell Drive, Suite 500, Reston, VA, 20191-4344. (c)2001 American Institute of Aeronautics & Astronautics or Published with Permission of Author(s) and/or Author(s)' Sponsoring Organization. SOLAR THERMAL PROPULSION FOR AN INTERSTELLAR PROBE Ronal . LymandW , Mar . EwingkE , Rames . KrishnanhS , Dea . LesternM , Thiokol Propulsion Brigham CityT U , Ralp . McNuthL . tJr Applied Physics Laboratory Johns Hopkins University 11100 Johns Hopkins Road Laurel, MD ABSTRACT readiness and required system improvements will be The conceptual design of an interstellar identified. probe powered by solar thermal propulsion (STP) INTRODUCTION engines will be described. This spacecraft will use Although traveling to the stars is still engineP ST perforo st m solar gravity assist witha beyon reace dth knowf ho n technology proba , e perihelion of 3 to 4 solar radii to achieve a solar capabl penetratinf eo interstellae gth r mediue b y mma system exit velocity (see Figure 1). Solar thermal within that reach.1"6 Usin ghigh-thrusta , high propulsion use sun'e sth s energ heayo t workina t g specific impulse solar thermal propulsion system this fluid wit moleculaw hlo r weight, suc hydrogens ha , probe could potentially achiev asymptotin ea c escape vero t y high temperatures (aroune d Th 300 .
    [Show full text]
  • Optimization of Low-Thrust Gravity Assist Trajectories Under Application of Tisserand’S Criterion
    OPTIMIZATION OF LOW-THRUST GRAVITY-ASSIST TRAJECTORIES UNDER APPLICATION OF TISSERAND’S CRITERION Vom Fachbereich Produktionstechnik der UNIVERSITÄT BREMEN zur Erlangung des Grades Doktor-Ingenieur genehmigte Dissertation von Dipl.-Ing. Volker Maiwald Gutachter: Prof. Dr.-Ing. Andreas Rittweger Fachbereich Produktionstechnik, Universität Bremen Prof. Dr.-Ing. Bernd Dachwald Fachbereich Luft- und Raumfahrttechnik, Fachhochschule Aachen Tag der mündlichen Prüfung: 23. März 2018 I think we are going… because it’s in the nature of the human being to face challenges. - Neil Armstrong Acknowledgements My gratitude goes towards my colleagues and friends from the System Analysis Space Segment Department for their support and fellowship. I thank especially Dominik Quantius, Andy Braukhane, Daniel Schubert and Conrad Zeidler for their discussions, support and proof reading concerning my dissertation. I also thank Dr. Marco Sharringhausen, who always helped me whenever I had questions about mathematics and Dr. Wolfgang Seboldt for stepping up on a short notice during the final stages of this thesis and providing much appreciated insight and advice. I thank my parents, Dr. Ursula and Werner Maiwald for enabling me to follow my dreams – I did, thank you! – and my brothers Rüdiger and Tobias for their encouragement. And of course, I thank my incredible wife Inga who patiently endured my struggles with my work and supported me so much, despite the bumpy road far off an optimal path. Revisions for Publication This is a revised edition of the original dissertation. The changes have been made to allow a more self-enclosed reading of the dissertation. Next to editorial changes to address language and typing errors, the following changes have been made to the originally submitted version: 1) The paper by Chen, Kloster and Longuski (2008) has been included in the discussion of Chapters 2.4.2, 6 and 7.11.4.
    [Show full text]