Reference System Issues in Binary Star Calculations Poster for Division A Meeting DAp.1.05 George H. Kaplan Consultant to U.S. Naval Observatory Washington, D.C., U.S.A. [email protected] or [email protected] IAU General Assembly Honolulu, August 2015 1 Poster DAp.1.05 Division!A! !!!!!!Coordinate!System!Issues!in!Binary!Star!Computa3ons! George!H.!Kaplan! Consultant!to!U.S.!Naval!Observatory!!([email protected][email protected])! It!has!been!es3mated!that!half!of!all!stars!are!components!of!binary!or!mul3ple!systems.!!Yet!the!number!of!known!orbits!for!astrometric!and!spectroscopic!binary!systems!together!is! less!than!7,000!(including!redundancies),!almost!all!of!them!for!bright!stars.!!A!new!genera3on!of!deep!allLsky!surveys!such!as!PanLSTARRS,!Gaia,!and!LSST!are!expected!to!lead!to!the! discovery!of!millions!of!new!systems.!!Although!for!many!of!these!systems,!the!orbits!may!be!poorly!sampled!ini3ally,!it!is!to!be!expected!that!combina3ons!of!new!and!old!data!sources! will!eventually!lead!to!many!more!orbits!being!known.!!As!a!result,!a!revolu3on!in!the!scien3fic!understanding!of!these!systems!may!be!upon!us.! The!current!database!of!visual!(astrometric)!binary!orbits!represents!them!rela3ve!to!the!“plane!of!the!sky”,!that!is,!the!plane!orthogonal!to!the!line!of!sight.!!Although!the!line!of!sight! to!stars!constantly!changes!due!to!proper!mo3on,!aberra3on,!and!other!effects,!there!is!no!agreed!upon!standard!for!what!line!of!sight!defines!the!orbital!reference!plane.!! Furthermore,!the!computa3on!of!differen3al!coordinates!(component!B!rela3ve!to!A)!for!a!given!date!must!be!based!on!the!binary!system’s!direc3on!at!that!date.!!Thus,!a!different! “plane!of!the!sky”!is!appropriate!for!each!such!date,!i.e.,!each!observa3on.!!However,!projec3on!effects!between!the!reference!planes,!differen3al!aberra3on,!and!the!curvature!of!the! sky!are!generally!neglected!in!such!computa3ons.!!Usually!the!only!correc3on!applied!is!for!the!change!in!the!north!direc3on!(posi3on!angle!zero)!due!to!precession!(and!some3mes! also!proper!mo3on).!!This!paper!will!present!an!algorithm!for!a!more!complete!model!of!the!geometry!involved,!and!will!show!that!such!a!model!is!necessary!to!avoid!errors!in!the! computed!observables!that!are!significant!at!modern!astrometric!accuracy.!!The!paper!will!also!suggest!where!conven3ons!need!to!be!established!to!avoid!ambigui3es!in!how!quan33es! related!to!binary!star!orbits!are!interpreted.! ! Mo#va#on' Ambiguous'Aspects'of'the'Orbital'Parameters' Other'Geometry'to'Consider'(cont.)' This!inves3ga3on!began!as!a!rela3vely!simple!project!to!add!a! Because!ORB6!is!a!collec3on!of!orbits!determined!by!different! Another!effect!arises!from!the!3LD!nature!of!the!rela3ve!posi3ons!of! binaryLstar!so[ware!module!to!the!Naval!Observatory!Vector! inves3gators,!it!is!a!heterogeneous!source!of!informa3on.!!Several! the!two!stars.!!For!nearby,!wide!systems,!the!fact!that!B!is!closer!(or! Astrometry!So[ware!(NOVAS).!!NOVAS!it!is!a!free!sourceLcode! aspects!of!the!data!in!ORB6!are!not!well!defined!and!therefore! farther)!than!A!affects!the!separa3on!we!measure;!!we!should!not! library!for!highLprecision!astrometric!calcula3ons!available!at! present!an!obstacle!to!precise!computa3ons.!!About!1/3!of!the!orbits! assume!a!viewpoint!at!infinity!for!these!systems.!! h]p://aa.usno.navy.mil/so[ware/novas/!! do!not!have!the!epoch!E"specified!!(most!of!the!specified!E’s!are! To"Earth,"assumed"at"infinity," in!Fortran,!C,!and!Python.!!It!is!likely! 2000).!!!Orbital!elements!P!and!T,!which!are!3mes,!are!o[en! at"*me"t" that!as!observa3onal!accuracy!inL! expressed!as!a!year!+!frac3on,!in!units!of!Besselian!years,!!but!the! creases,!binary!stars!will!become! conversion!to!more!standard!3me!units!is!not!given.!!In!a!more!subtle! B" an!increasingly!important!issue!in!! issue,!it!is!unclear!what!line!of!sight!should!be!used!to!define!the! sep"2" defining!celes3al!reference!frames.!!For!example,!even!when! “plane!of!the!sky”.!!Stars!move!on!the!celes3al!sphere!(i.e.,!with! sep"2" sep"1" i’ observa3ons!span!only!a!few!years,!the!components!of!wide!pairs! respect!to!the!ICRS!axes)!due!to!proper!mo3on,!aberra3on,!and!other! δ Earth"at"*me"t" d" " with!orbital!periods!of!centuries!will!show!detectable!curvatures!of! effects.!!Aberra3on!shi[s!our!line!of!sight!to!every!star,!regardless!of! (distance"Earth"to"star"A)" A" their!paths!if!the!observa3onal!accuracy!is!great!enough!(Kaplan!&! distance,!by!41!arcseconds!over!the!course!of!every!year.!! sep 1 = (B−A) cos i’ = cos (i + Δϕ) = 70.7038 arcsec Makarov!2004,!AN!324,!419;!Makarov!&!Kaplan!2005,!AJ!129,!2420).!! ! (from previous figure) sep 1 (B−A) cos i’ NOVAS!should!deal!with!all!types!of!astrometric!binary!systems.! To!move!forward,!I!had!to!make!several!(arbitrary)!assump3ons:! sep 2 = arctan = arctan d − δ d − (B−A) sin i ! ’ Plane"of"the" ★ If!E!is!not!specified,!then!E=2000,!i.e.,!E!=!J2000.0!=!JD!2451545.0!(TT).! sky"at"*me"t" The!binary!star!module!that!is!planned!for!NOVAS!and!is!the!subject! Again suppose B−A = 100 arcsec, i=45°, i’=45°.0555°, ★ (="*me"of" of!this!paper!will!provide!both!the!absolute!and!rela3ve!posi3ons!of! The!conversion!from!Besselian!to!Julian!epochs!and!3me!units!is!done! and d = 100 pc. ephemeris"or" observa*on)" the!two!components!of!a!binary!system,!at!a!requested!date,!given! according!to!a!conven3onal!expression!(using!fixedLlength!years)! Then sep 1 − sep 2 = 70.7038 arcsec − 70.7281 arcsec recommended!by!the!IAU!in!1976.! = −24.2 mas known!orbital!elements!and!standard!astrometric!parameters.!!Thus,! it!covers!what!may!called!the!“forward!problem”!of!producing!binary! ★ The!plane!of!the!sky!for!the!orbit!is!orthogonal!to!the!direc3on!toward! the!A!component!of!the!system!defined!by!its!mean"place"at"epoch"E.! star!ephemerides!given!an!orbit!already!determined.!!It!does!not! One!more!effect!to!consider:!differen3al!aberra3on!between!the! deal!with!the!“backward!problem”!of!determining!an!orbit!from! ★ The!direc3on!of!the!orbital!node,!i.e.,!the!sign!of!the!inclina3on,!is!! two!components.! known!(even!though!in!many!cases!it!is!ambiguous).! observa3ons.!!However,!once!an!ini3al!orbit!has!been!determined! aberrated geometric ★ Units!of!arcseconds!used!for!the!semimajor!axis!a!should!be!considered! Aberration moves the apparent places of for!a!binary!system,!followLup!observa3ons!can!improve!the!orbit! positions positions the two stars slightly differently on the to!be!a!measure!of!length;!to!convert!to!astronomical!units!(au),!divide! sky, which can affect both separations through!differen3al!correc3on;!!and!a!correct!formula3on!of!the! B" B" by!the!parallax!of!the!system.! v/c" and position angles. Differential forward!problem!is!necessary!to!avoid!systema3c!errors!in!the! projected"at"B" aberration can amount to sep"2" 1 mas for 10 arcsec separations and improved!orbital!elements.!!!That!is,!the!OLC’s!should!be!as!correct! !!!!!!!!!'''It'would'be'quite'useful'if'the'appropriate'IAU'body'would' sep"3" 10 mas for 100 arcsec separations. In this case, with sep 2 = 70.2781 arcsec, as!possible.!!Given!perfect!orbital!elements,!the!output!of!the! A" A" establish'conven#ons'such'as'the'above'for'future'binary'star'work.'' v/c" sep 2 − sep 3 could be up to ≈ 7 mas. NOVAS!module!(the!C’s)!should!be!accurate!to!be]er!than!1!mas.!!!! projected"at"A" B" ! Binary!stars!are!perhaps!the!purest!example!of!the!twoLbody! Geometry'Usually'Considered' Earth’s"orbital"velocity" problem!in!astronomy,!and!the!computa3on!of!the!rela3ve!posi3ons! v"" In!compu3ng!observables!from!orbital!elements,!we!can!take!the! A" of!the!two!components!in!a!local!spaceLfixed!reference!system!is! solu3on!of!the!twoLbody!problem!for!binary!stars!as!a!given.!! nonLcontroversial.!!However,!taking!that!informa3on!and!predic3ng!! Standard!texts!such!as!Aitken!(1964),!The"Binary"Stars,!then!advise! the!values!of!astrometric!parameters!for!an!EarthLbased!observer! correc3ng!for!the!change!in!the!north!direc3on,!due!to!precession,! requires!a!more!sophis3cated!model!of!the!whole!geometry!than!is! between!the!epoch!of!the!orbit!(E)!and!the!epoch!of!observa3on,! usually!applied.!!In!par3cular,!at!high!levels!of!astrometric!accuracy,! and!similarly!for!the!change!in!north!direc3on!due!to!the!star’s! NOVAS'Algorithm' the!3LD!nature!of!the!problem!comes!into!play!in!surprising!ways.! proper!mo3on!(although!aberra3on!is!usually!a!larger!effect).!! The!algorithm!adopted!for!the!NOVAS!double!star!algorithm!can!be! Simple!formulas!are!provided.!!Some!texts!also!give!correc3ons!for! summarized!as!follows.!!Some!details!(such!as!change!of!units)!are! change!of!scale!and!change!of!lightL3me,!for!those!systems!with! omi]ed,!and!for!prac3cal!reasons!the!order!of!computa3ons!in!the! Available'Orbital'Data'for'Binary'Stars' significant!radial!veloci3es.! so[ware!is!somewhat!different.!!! The!data!set!used!is!the!Sixth!Catalog!of!Orbits!of!Visual!Binary!Stars,! 1) Using!standard!astrometric!parameters!(e.g.,!from!the!Hipparcos!catalog)! 2 available!on!CD!from!the!U.S.!Naval!Observatory!and!online!at! Other'Geometry'To'Consider' for!component!A,!or!the!system’s!center!of!mass,!compute:! h]p://www.usno.navy.mil/USNO/astrometry/op3calLIRLprod/wds/orb6!! ''''''''''''a ,!the!mean!place!direc3on!vector!of!A!at!epoch!E"(epoch!of!orbit),! (it!is!an!update!to!Hartkopf,!Mason,!&!Worley!2001,!AJ!122,!3472).1!! In!the!module!developed!for!NOVAS,!the!change!in!north!direc3on!is! E part!of!a!general!transforma3on!between!two!3LD!reference!systems! ''''''''''''aO,!the!mean!place!direc3on!vector!of!A!at!the!3me!of!observa3on,! The!catalog!provides!2518!orbits!of!2413!systems,!although!not!all!