Arxiv:Astro-Ph/0608273 V1 13 Aug 2006 Ln Oteotclpt Ea Ocuetemanuscript the Conclude IV

Provided byCERNDocumentServer CORE arXiv:astro-ph/0608273 v1 13 Aug 2006 ln oteotclpt ea ocuetemanuscript the conclude IV. delay Section path tangent in optical the the in (star) to relating parallactic axis plane on optical the comments the Some if from away these. way distances bridge hori- hybrid to or used a is equatorial in angle position in or baseline purely systems. coordinates the done zontal coordinate be computes polar can and angle—which familiar defines III the Section in variables of (NCP) Pole Celestial star. North point celestial the the pivot to of from the relocated where definition been system a has coordinate resolved (ii) a in poorly and longitudes or- and interferometer, the well the are detailing by which Interferometry 16], directions Optical [15, thogonal the (OIFITS) of Format tables Exchange visibility the with oa nl nthe system’s in coordinate polar angle a Earth. of polar on introduction the location (i) observatory’s is coordinates, This the which sky of coordinates to independent equatorial Celestial baseline are in the the directions onto ties defines this which projection that this is baseline. use Sphere projected to oriented incentive an the The meets defines baseline the stretched Sphere between the section Celestial where its point the and This and circle, star Sphere. great Celestial the a the is and with intersection interferometer intersection stellar its a and star, of baseline the contains URL: ∗ lcrncades [email protected]; address: Electronic noeve,ScinI nrdcstesadr set standard the introduces II Section overview, In h ae ecie tnadt en h ln that plane the define to standard a describes paper The http://www.strw.leidenuniv.nl/~mathar trpstoswihaeaindwt rotooa otegr orthog the to two orthogonal represents or angle with This aligned are tangentia star. which manus plane the positions celest The star a of the in position light. onto measured the outwards stellar at defined, vecto is further the delay baseline baseline of the projected wavefronts the are the vectors projecting of derived of Two plane the star. in a to direction a ewrs rjce aeie aalci nl;positio angle; interferometer parallactic baseline; projected Keywords: 95.75.Kk 95.10.Jk, angle. numbers: “parallactic” PACS di the reference of the definition as common chosen the is to Pole adapted Celestial North The plane. h ai emtyo tla nefrmtrwt w tele two with interferometer stellar a of geometry basic The u peia rgnmtyo h rjce aeieAngle Baseline Projected the of Trigonometry Spherical − .SCOPE I. ednOsraoy ..Bx91,20 ALie,TeNethe The Leiden, RA 2300 9513, Box P.O. Observatory, Leiden v oriaesse associated system coordinate Dtd uut1,2006) 15, August (Dated: ihr .Mathar J. Richard nl;clsilshr;shrclatooy stellar astronomy; spherical sphere; celestial angle; n epniua oti iceis circle this to perpendicular n liueaoesalevel), sea above altitude and ee tteErhcne,connects center, Earth the at tered r eae otegorpi longitude interferometer geographic stellar the a to of related telescopes two are the of dinates ute on ra iceo radius of circle used great be A will that down. coordinates further Cartesian other from apart eada“ a add We n ffcieErhradius Earth effective and aeieaetr angle aperture baseline a h aeievector baseline the sse rmteErhcenter: Earth the from seen is h idrs ieto from direction rose wind The olw:dfieatnetpaet h at at Earth the to plane tangent a define follows: rhnra etr htsa hspaeare plane this span that vectors orthonormal nagoeti oriaesse,teCreincoor- Cartesian the system, coordinate geocentric a In N ˆ 1 T cos ≡ 1 deto h ea eoddi the in recorded delay the of adient ∗ a pee h oiinageo the of angle position The sphere. ial eto ftepoetdbsln angle, baseline projected the of rection oteclsilshr,tnetpoint tangent sphere, celestial the to l cpscnit fabsln etrand vector baseline a of consists scopes · Z   nldrcin ntesy differential sky, the on directions onal T T ,adtepoetdbsln vector baseline projected the and r, o Φ cos = − − 2 i g rp el ihtetrigonometry the with deals cript otegoeti oriae ostthem set to coordinates geocentric the to ” ≡ = cos sin .ShrclApproximation Spherical 1. o Φ cos I IEGEOMETRY SITE II. ρ .TlsoePositions Telescope A. | T λ λ   1 1 1 1 | | i Φ sin i Φ sin cos o Φ cos sin 1 T b 2 i Φ sin λ λ | = Metadata, citationandsimilarpapersatcore.ac.uk 1 i 1 i cos 2 Z o Φ cos o Φ cos   T cos( rlands a endfie ne which under defined been has i Z ρ 2 g (3) ; − ∼ J ˆ λ smo h at radius Earth the of (sum i i T 2 ≡ T   ∂ 1 − T 1 g to ρ λ 1 2 i , T /∂ sin 1 1 T i Φ sin + ) 1 ρ Φ Z 2 u and ≈ 1 = 1 − T scmue as computed is λ ; 30k,cen- km, 6380 i v 1 aiueΦ latitude , , T × 2 | N 2 ˆ . 1 T vector a ; T 1 i Φ sin 2 | where , 1 1 = Two . 2 (1) (2) (4) (5) . . i 2 to the North and for example

cos λ1 C    x  Eˆ 1 ≡ − sin λ1 ∼ ∂T1/∂λ1; |Eˆ 1| =1. (6) C = Cy , (9) 0 C  g  z g to the East. The unit vector from T1 to T2 along the plus local telescope coordinates 1 great circle is Jˆ × T1. The compass rose angle τ is ρ T ≡ C + t , i =1, 2. (10) Jˆ 1 T Nˆ i i computed by its decomposition × ρ 1 = cos τ 1 + Leaving aside geodetic variables, the array’s horizontal sin τEˆ 1 within the tangent plane. It becomes zero if T2 plane could be deﬁned as stretching from the center C is North of T1, and becomes π/2 if T2 is East of T1: Fig. C 1. orthogonal to to the directions called North and East, associating spherical coordinates North C ρ ≡ C2 + C2 + C2, λ = arctan y , (11) T2 q x y z Ab−A Cx C Φ = arcsin z . (12) 2 2 2 Cx + Cy + Cz Ab τ>0 q East with the array center by inversion of (1).

A>0 B. Sky coordinates T1 In a plane orthogonal to C this deﬁnes a star azimuth star A, a zenith angle z, and a star elevation a = π/2 − z, − cos A sin z − cos A cos a FIG. 1: In the horizontal coordinate system, the baseline di- s =  sin A sin z  =  sin A cos a  . (13) rection is characterized by the azimuth angle Ab. In the case cos z sin a shown, 0

2. Caveats cos a sin A = cos δ sin h; (14) cos a cos A = − sin δ cosΦ+ cos δ cos h sin Φ; (15) The angle τ of the previous section does not transform sin a = sin δ sinΦ + cos δ cos h cosΦ; (16) exactly into τ ± π if the roles of the two telescopes are cos δ cos h = sin a cosΦ+ cos a cos A sin Φ; (17) swapped, because the great circle, which has been used to define the direction, is not a loxodrome, This asym- sin δ = sin a sinΦ − cos a cos A cosΦ. (18) metry within the definition indicates that a more generic These convert (13) into framework to represent the geometry is useful. If we consider (i) geodetic rather than geocentric rep- cosΦsin δ − sinΦcos δ cos h resentations of geographic latitudes [14] or (ii) long base- s =  cos δ sin h  . (19) line interferometers build into rugged landscapes, the 5 sin Φ sin δ + cosΦcos δ cos h  t parameters above (ρ, λi, Φi) are too constrained to han- dle the 6 degrees of freedom of two earth-fixed telescope A star at hour angle h = 0 is on the Meridian, positions. Nevertheless, a telescope array “platform” is helpful to define a zenith, a horizon, and associated co- sin(δ − Φ)   ordinates like the zenith distance z or the star azimuth s = 0 , (h = 0), (20) C cos(δ − Φ) A. This leads to the OIFITS concept, an array center ,  t 3 north of the zenith at A = π if δ − Φ > 0, south at A =0 The standard definition of the delay vector D and pro- if δ − Φ < 0. A reference point on the Celestial Sphere jected baseline vector P is is the North Celestial Pole (NCP), given by insertion of δ = π/2 into (19), b = D + P; D k s; P ⊥ D. (30)

cosΦ P—and later its circle segment projected on the Celestial   b s+ = 0 . (21) Sphere—inherits its direction from such that heads and  sinΦ  tails of the vectors are associated with the same portion t of the wavefront: Fig. 2. We use the term “projected The angular distance between s and s+ is s s · s+ = sin δ. (22) D P The transformation from the local coordinate system to the geocentric coordinate system is [13] θ T1 b T2 T t 0  ix   ix    Tiy = Ugt · tiy + 0 , (23) FIG. 2: Sign conventions of the co-planar vectors b, D and T t ρ  iz g  iz t   P. via a rotation matrix baseline” both ways: for the straight vector (of length D b − sinΦcos λ sin λ cosΦcos λ P , units meter) that connects the tails of and , or   the line segment (length θ, units radian) on the Celestial Ugt(Φ, λ)= − sin Φ sin λ − cos λ cosΦsin λ . Sphere.  cosΦ 0 sinΦ  (24) The dot product of (13) by (28) is The inverse transformation proceeds by multiplying (23) D = s · b = b cos θ, (0 ≤ θ ≤ π), (31) from the left by the inverse matrix (which equals the transpose matrix) where the angular distance θ between the baseline and −1 † star directions is introduced as Utg ≡ Ugt = Ugt. (25) cos θ = cos ab cos a cos(Ab − A) + sin ab sin a. (32)

C. Baseline: Generic Position The star circles around the NCP in 24 hours, which changes the distance to the baseline in a centric periodic The baseline vector coordinates of the geocentric way (App. B). OIFITS system

B III. PROJECTED BASELINE ANGLES  x  b = By = T2 − T1 (26) B  z g A. Definition T T (defined to stretch from 1 to 2—the opposite sign is We define position angles of points on the Celestial also in use [17]) could be converted with Sphere as the bearing angle by which the NCP must be rotated around the star direction s (the axis) until the bx t2x t1x Bx         NCP and the particular point of the object are aligned b = by = t2y − t1y = U (Φ, λ)· By tg to the same direction (along a celestial circle) from the  bz   t2z   t1z   Bz  t t t g star. The sign convention is left-handed placing the head (27) of s at the center of the Celestial Sphere. This is equiv- to the topocentric system. Just like the star direction alent to drawing a line between the star and the NCP, bˆ (13), the unit vector along the baseline direction defines looking at it from inside the sphere, and defining posi- a baseline azimuth Ab and a baseline elevation ab tion angles of points as the angles in polar coordinates in the mathematical sign convention, using this line as − cos A cos a b b b the abscissa and the star as the center. Another, fully bˆ ≡ ≡  sin A cos a  . (28) b b equivalent definition uses a right-handed turn of the ob- b sin a  b t ject around the star until it is North of the star. And The baseline in most optical stellar interferometers is ap- finally, as an aid to memory, it is also the value of the proximately horizontal, nautical course while navigating a ship sailing the outer hull of the Celestial Sphere, which is currently poised at ab ≈ 0. (29) the star’s coordinates. 4

This defines values modulo 2π; whether these are fi- projection [3] we consider the two directions from the nally represented as numbers in the interval [0, 2π) or in origin toward the NCP, s+, on one hand and toward T2, the interval (−π, +π]—as the SLA BEAR and the SLA PA bˆ, on the other. The tangential plane is fixed by any two routines of the SLA library [22] do—is largely a matter unit vectors perpendicular to the vector s of (19). Rather of taste. arbitrarily they are chosen along ∂s/∂A and ∂s/∂a, ex- An overview on the relevant geometry is given in Fig- plicitly ure 3 which looks at the celestial sphere from outside—a sin A truly classical perspective. In Figure 3, the baseline has   eA = cos A , (34) 0  t zenith cos A sin a   ea = − sin A sin a . (35) cos a  t z (The equivalent exercise with a different choice of axes follows in Section III D.) These are orthonormal,

s·eA = s·ea = eA·ea = 0; |eA| = |ea| = 1; s×ea = eA. (36) NCP E T1 star Decomposition of the direction from s to the NCP deﬁnes three coeﬃcients c0,1,2,

S s+ = c0s + c1e + c2e . (37) θ a A a N Building the square on both sides yields the familiar formula for the sum of squares of the projected cosines, A -A 2 2 2 W T2 b c0 + c1 + c2 =1. (38) Two dot products of (37) using (21) and (34)–(35) solve FIG. 3: Celestial sphere, seen from the outside. The north for two coeﬃcients, direction through the object is given by the great circle pass- c1 = s+ · eA = cosΦsin A, (39) ing through the celestial poles and the object. The angle θ is s e the vector P = b − D projected on the sphere. c2 = + · a = cosΦcos A sin a + sinΦcos a. (40)

An angle angle ϕ between eA and s+ is defined via Fig. been infinitely extended straight outwards in both di- 4, rections which defines telescope coordinates T1 and T2 c1 also where the baseline penetrates the Celestial Sphere. cos ϕ = , (41) c2 + c2 Standing at the mid-point of the baseline, telescope 1 p 1 2 c2 then is at azimuth τ, telescope 2 at azimuth Ab = τ + π, sin ϕ = . (42) see Fig. 1. Let the object be at azimuth A and elevation c2 + c2 p 1 2 a. From (38), and since (22) equals c0, The projection of the baseline occurs in a plane in- cluding the object and the baseline and thus defining the c2 + c2 = cos δ. (43) great circle labeled θ in Figure 3. If θ is the angle on the q 1 2 sky between the object and the point on the horizon with The equivalent splitting of bˆ into components within and azimuth Ab, then the length P of the projected baseline perpendicular to the tangent plane is: is given by ˆ ′ ′ ′ b = c0s + c1eA + c2ea. (44) P = b sin θ, (33) Building the square on both sides yields where θ is calculated from the relation (32). With this ′2 ′2 ′2 auxiliary quantity, there are various ways of obtaining c0 + c1 + c2 =1. (45) the position angle pb of the projected baseline on the Dot products of (44) solve for the expansion coefficients, sky, some of which are detailed in Sections III B–III D. with (28) and (34)–(35): ′ ˆ c1 = b · eA B. In Horizontal Coordinates = − cos Ab cos ab sin A + sin Ab cos ab cos A

= cos ab sin(Ab − A), (46) We use a vector algebraic method to span a plane tan- ′ bˆ e gential to the celestial sphere; tangent point is the po- c2 = · a sition s of the star. Within this orthographic zenithal = − cos ab sin a cos(Ab − A) + sin ab cos a. (47) 5

ea and

′ ′ NCP c1c1 + c2c2 = cosΦ(cos A cos a cos θ − cos ab cos Ab) c 2 ′ +c2 sinΦcos a. (55) p The computational strategy is to build

′ ′ c1c2 − c1c2 tan pb = ′ ′ . (56) c1c1 + c2c2 p b ϕ The (positive) values of cos δ and sin θ in the denomi- nators of (52) and (53) need not to be calculated. The branch ambiguity of the arctan is typically handled by use of the atan2() functionality in the libraries of higher b c' programming languages (C/C++, Fortran, Java, tcl, 2 Maple, Python, IDL,. . . ). ϕ' pb C. Via Parallactic Angle star eA c 1 c'1 ψ If the position angle of the zenith p is known by any other sources—see App. A—a quicker approach to the calculation of pb employs an auxiliary angle ψ [10], FIG. 4: In the tangent plane spanned by the unit vectors ′ (34) and (35), the components of bˆ and s+ are c1, c1, c2 and pb ≡ p + π + ψ (mod 2π). (57) ′ c2 deﬁned by (37) and (44). The view is from the outside onto the plane, so eA, the axis vector into the direction of The calculation of ψ is delegated to the calculation of its increasing A, points to the left. sines and cosines

cos ψ = − cos pb cos p − sin pb sin p; (58) This defines an angle ϕ′ in the polar coordinates of the sin ψ = − sin pb cos p + cos pb sin p. (59) tangent plane, c′ In the right hand sides we insert (52)–(53) and (A1)– cos ϕ′ = 1 , (48) ′2 ′2 (A2), and after some standard manipulations c1 + c2 p ′ c cos ab sin(Ab − A) sin ϕ′ = 2 . (49) sin ψ = ; (60) ′2 ′2 c1 + c2 sin θ p sin a cos θ − sin ab ′ bˆ s cos ψ = Since c0 equals · = cos θ in (32), sin θ cos a cos ab sin a cos(Ab − A) − sin ab cos a ′2 ′2 = . (61) qc1 + c2 = sin θ. (50) sin θ The baseline position angle is the difference between the [These are c′ and −c′ of (46) and (47) divided by sin θ.] e 1 2 two angles, as to redefine the reference direction from A In Fig. 4, ψ can therefore be identified as the angle be- to the direction of the NCP (see Fig. 4): tween −ea and bˆ, and this is consistent with (57) which ′ pb = ϕ − ϕ (mod 2π). (51) decomposes the rotation into the angle p, a rotation by the angle π (which would join the end of p with the start Sine and cosine of this imply [1, 4.3.16, 4.3.17] of ψ in Fig. 4), and finally a rotation by ψ. ′ ′ In summary, the disadvantage of this approach is that sin pb = sin ϕ cos ϕ − cos ϕ sin ϕ ′ ′ p must be known by other means, and the advantage is c1c − c c2 = 2 1 ; (52) that computation of ψ via the arctan of (60) and (61) cos δ sin θ ′ ′ ′ ′ only needs c1 and c2, but not c1 or c2. cos pb = cos ϕ cos ϕ + sin ϕ sin ϕ ′ ′ c1c + c2c = 1 2 , (53) cos δ sin θ D. In Equatorial Coordinates where ′ ′ Section III B provides pb, given star coordinates (A, a). c1c2 − c1c2 = cos Φ(sin A sin ab cos a − cos ab sin a sin Ab) Subsequently we derive the same value in terms of δ and − sinΦcos a cos ab sin(Ab − A), (54) h. The straight forward approach is the substitution of 6

(14)–(18) in (54) and (55). A less exhaustive alternative (28) yields works as follows [8]: the axes eA and ea in the tangential ′ plane are replaced by two different orthonormal direc- d2 = − cos Ab cos ab sin Φ sin h + sin Ab cos ab cos h tions, better adapted to δ and h, namely the directions − sin ab cosΦsin h. (68) along along ∂s/∂h and ∂s/∂δ: Building the dot product of (67) with eδ yields sin Φ sin h   ′ eh = cos h ; (62) d1 = − cos Ab cos ab cosΦcos δ − cosΦsin h  t − cos Ab cos ab sin Φ sin δ cos h cosΦcos δ + sin Φ sin δ cos h − sin Ab cos ab sin δ sin h + sin ab sinΦcos δ e   δ = − sin δ sin h . (63) − sin ab cosΦsin δ cos h. (69) sinΦcos δ − cosΦsin δ cos h  t Some simplification in this formula is obtained by trading The further calculation follows the scheme of Section the hour angle h for the angle θ, which might be more IIIB: The axes are orthonormal, readily available since θ is measured via the delay: In the second term we use (15) to substitute s·eh = s·eδ = eh·eδ = 0; |eh| = |eδ| = 1; s×eδ = eh. (64) cosΦsin δ + cos A cos a sinΦcos h → . (70) The decomposition of s+ in this system defines three ex- cos δ pansion coefficients d0,1,2, In the third term we use (14) to substitute s+ = d0s + d1eδ + d2eh. (65) sin A cos a sin h → , (71) Multiplying this equation in turn with eh and eδ yields cos δ

d2 = 0; d1 = s+ · eδ = cos δ. (66) and in the last term we use (16) to substitute sin a − sin δ sinΦ The direction bˆ to T2 is projected into the same tangent cosΦcos h → . (72) ′ cos δ plane, deﬁning three expansion coeﬃcients d0,1,2, Further standard trigonometric identities [1, bˆ ′ s ′ e ′ e = d0 + d1 δ + d2 h. (67) 4.3.10,4.3.17] and (32) yield

Building the dot product of this equation with eh using sin a sinΦ − cos a cos A cosΦ − sin δ cos θ d′ = b b b . (73) 1 cos δ

eδ In the numerator we recognize a baseline declination δb,

s+ · bˆ ≡ sin δb = sin ab sinΦ − cos ab cos Ab cosΦ. (74)

pb is the angle that rotates the projected vector (d2, d1) within the tangent plane into the direction (d′ , d′ ) (Fig. pb NCP 2 1 d1 5),

d′ cos p = 1 , (75) d'1 b sin θ b d′ sin p = − 2 . (76) b sin θ Again, the sin θ does not need actually to be calculated, but only star d' −d′ 2 tan p = 2 , (77) eh b d′ pb 1 and again, selection of the correct branch of the arctan is easy with atan2() functions if the negative sign is kept ′ attached to d2. FIG. 5: Fig. 4 after switching from the (eA, ea) axes to Note also that the use of (73) is optional: d′ and d′ ˆ 2 1 (eh, eδ). Axis components of s+ and b are indicated. are completely deﬁned in terms of the baseline direction 7

(Ab,ab), the geographic latitude Φ and the star coordi- consistent with the decomposition in Fig. 4. The two nates (h,δ) via (68) and (69). With these one can imme- directions of constant delay are characterized by ∆θ = 0. diately proceed to (77). Equating the right hand side of the previous equation OIFITS [15, 16] deﬁnes no associated angle in the plane with zero, this means perpendicular to the direction of the phase center. ∆a The segment of the baseline in the tangent plane pri- = tan ψ (∆θ = ∆D = 0). (80) marily deﬁnes an orientation in the u−v plane. However, cos a∆A modal decomposition of the amplitudes in products of For the direction of maximum change in ∆D, which radial and angular basis functions (akin to Zernike poly- runs perpendicular to (80), nomials) means that the Fourier transform to the x − y plane preserves the angular basis functions (see the cal- ∆a 1 = − , (pb : max ∆D) (81) culation for the 3D case, Spherical Harmonics, in [12]). cos a ∆A tan ψ In this sense, the angle p is also “applicable” in the x−y b which leads to plane. ∆D = −b sin θ∆θ cos a ∆A ∆a = b sin θ = −b sin θ , (max ∆D(82)). IV. DIFFERENTIALS IN THE FIELD-OF-VIEW sin ψ cos ψ

′′ A. Sign Convention Example: a field of view of ∆θ = ±1 at a baseline of b = 100 m at a “random” position θ = 45◦ scans ∆D = ±340 µm. The details of the optics determine how The provisions of the previous section, namely far this contributes to the instrumental visibility (loss). • the OIFITS sign convention of the baseline vector, Fig. 6 sketches this geometry: the change ∆D when re- Eq. (2) and Fig. 2, pointing from one star by (∆A, ∆a) to another star along some segment of the Celestial Sphere can be split into a • the conventional formula (31) change associated with the radial direction toward/away from T2 along a great circle, and a zero change moving fix the sign of the delay D as follows: if the wavefront along a circle of radius cos θ centered on the baseline hits T2 prior to T1, the values of D and cos θ are positive; azimuth. if it hits T1 prior to T2, both values are negative. The transitional case D = 0 (θ = π/2) occurs if the zenith star passes through the “baseline meridian” plane perpendicular to the baseline; the two lines to the Northwest and Southeast in Fig. 1 show (projections of) these direc- star 1 tions. With (29) and (31) we see that this case is approx- ∆θ imately equivalent to cos(Ab − A) ≈ 0. Since cos a > 0, star 2 the sign of D coincides with the sign of cos(Ab −A). This is probably the fastest way to obtain the sign of the Op- tical Path Difference from FITS [6] header keywords; the advantage of this recipe is that the formula is independent of which azimuth convention is actually in use. The baseline length b follows immediately from the Eu- clidean distance of the two telescopes entries for column Ab STAXYZ of the OIFITS table OI ARRAY [16].

FIG. 6: Example of three concentric circles on the Celestial B. External Path Delay Sphere, centered at the baseline, which unite star directions of D = const. Two star positions are connected by the short di- agonal dash. A quarter of each of the two projected baselines, The total diﬀerential of (32) relates a direction which intersect at (Ab,ab) near the horizon, is also shown. (∆A, ∆a) away from the star to a change in the angular distance between the star and the baseline, We may repeat the calculation in the equatorial system and split ∆θ into ∆δ and ∆h, − sin θ∆θ = ∆a [− cos ab sin a cos(Ab − A) + sin ab cos a]

+∆A cos a cos ab sin(Ab − A). (78) ∆D = −b sin θ∆θ (83)

= b(∆δeδ · bˆ + ∆h cos δeh · bˆ) (84) (60) and (61) turn this into ′ ′ = b(∆δ d1 + ∆h cos δ d2) (85)

sin θ∆θ = sin θ(cos ψ∆a − sin ψ cos a∆A), (79) = b sin θ(∆δ cos pb − ∆h cos δ sin pb). (86) 8

This translates [3, App. C] to our variables. When The ratio of these two equations is ∆θ = 0, we have the directions of zero change in ∆D, perpendicular to the projected baseline, cosΦsin A tan p = . (A3) ∆δ cosΦcos A sin a + sinΦcos a = tan pb, (∆θ = ∆D = 0). (87) cos δ∆h This expression can be transformed from the horizontal The direction of maximum change is to equatorial coordinates as follows: Multiply numerator ∆δ 1 and denominator by cos a and use (14) in the numerator, = − , (pb : max ∆D), (88) cos δ∆h tan pb cosΦcos δ sin h tan p = . (A4) and along this gradient cosΦcos A sin a cos a + sinΦcos2 a cos δ∆h ∆δ ∆D = −b sin θ = b sin θ , (max ∆D). (89) Insert (18) in the denominator, sin pb cos pb cosΦcos δ sin h tan p = V. SUMMARY sin a(sin a sinΦ − sin δ)+sinΦcos2 a cosΦcos δ sin h = . (A5) To standardize the nomenclature, we propose to choose sinΦ − sin a sin δ the North Celestial Pole as the “reference” direction (direction of position angle zero) and a handedness to define Replace sin a with (16) in the denominator, eventually the sign of the position angles (mathematically positive divide numerator and denominator through cos Φ cos δ: if looking at the Celestial Sphere from the inside). cosΦcos δ sin h The set of position angles discussed here includes tan p = sinΦ − (sin δ sinΦ + cos δ cos h cosΦ)sin δ • the direction of the zenith, the “parallactic” angle, cosΦcos δ sin h = • the direction toward the point where the baseline cos2 δ sinΦ − cos δ sin δ cos h cosΦ pinpoints the Celestial Sphere, the “projected base- sin h line angle,” = . (A6) cos δ tanΦ − sin δ cos h • position angles of “secondary” stars in differential astrometry. These manipulations involved only multiplications with positive factors like cosΦ, cos δ and cos a; therefore one Computation of the angle proceeds via (56) if the object can use sin h and cos δ tanΦ − sin δ cos h separately as coordinates are given in the conventional system, or via the two arguments of the atan2 function to resolve the (77) if they are given in the equatorial system. branch ambiguity of the inverse tangent. The mathematics involved is applicable to observato- This definition of the parallactic angle coincides with ries at the Northern and the Southern hemisphere: posi- the dominant one in the literature [2, 18, 19, 20, 21]; tions (A, a) or (δ, h) are mapped onto a 2D zenithal co- others exist: Equation (1) in [4] generates the 90◦- ordinate system. The position angles play the role of the complement of our definition, for example. longitude (the North Celestial Pole the role of Greenwich If we insert (61)–(60) into (79) [1, 4.3.18] and compare or Aries). A distance between the object and the star (in with (86), we find the range from 0 to π along great circles through the Ce- lestial Pole in an ARC projection) might take the role of ∆a + tan p the polar distance—although an orthographic SIN pro- ∆δ cos a∆A = ∆a , (A7) jection had been used for the calculations in this script. cos δ∆h 1 − tan p cos a∆A

which means the parallactic angle mediates between di- APPENDIX A: PARALLACTIC ANGLE rections expressed as axis ratios in the two coordinate systems [11]. We deﬁne the parallactic angle p as the position angle of the zenith. Since the formulas in Section III dealt with the general case of coordinates on the Celestial Sphere, APPENDIX B: DIURNAL MOTION ′ one may just set ab = π/2 in (28), which induces c1 = 0, ′ c2 = cos a. This also turns (32) into cos θ = sin a. (52)– The diurnal motion of the external path diﬀerence is (53) become described by separating the terms ∝ sin h and ∝ cos h. cosΦsin A We rewrite (31) in equatorial coordinates by multiplying sin p = ; (A1) cos δ (19) and (28) and collect terms with the aid of (74): cosΦcos A sin a + sinΦcos a cos p = . (A2) s b cos δ D = · = b[sin δ sin δb + cos δ cos δb cos(h − hb)]. (B1) 9

A time-independent offset b sin δ sin δb and an amplitude therefore these closure relations of the projected baseline b cos δ cos δb, defined by products of the polar and equa- angles are correct for (i) both senses of defining their torial components of star and baseline, respectively, con- orientation, and (ii) any reference direction of the zero stitute the regular motion of the delay. The phase lag hb angle. is determined by the arctan of

cos δb sin hb = sin Ab cos ab; (B2)

cos δb cos hb = cos Ab cos ab sin Φ + sin ab cosΦ,(B3) APPENDIX D: NOTATIONS and plays the role of an interferometric hour angle.

α right ascension APPENDIX C: CLOSURE RELATIONS a, ∆a star elevation and its diﬀerence A, ∆A star azimuth angle, and its diﬀerence

An array T1, T2 and T3 of telescopes forms a baseline Ab baseline azimuth triangle b,b baseline vector and length

b12 + b23 + b31 =0. (C1) δ declination D signed external path diﬀerence s If all point to a common direction , some phase closure eh, eδ, eA, ea unit vector in tangent plane sum rules result: Φ geographic latitude h hour angle D12 + D23 + D31 =0; (C2) l local sidereal time D12 + D23 + D31 =0; (C3) λ geographic longitude b12 cos θ12 + b23 cos θ23 + b31 cos θ31 =0. (C4) p position angle of the zenith, “parallactic The triangle of projected baselines in the tangent plane angle” is pb position angle of the projected baseline P projected baseline length P P P 12 + 23 + 31 =0. (C5) ρ Earth radius s These vectors can be represented as numbers in a com- unit vector into star direction plex plane by introduction of three moduli P and three s+ unit vector into NCP direction orientation angles pb: θ angular separation between star and baseline vector ipb12 ipb23 ipb31 P12e + P23e + P31e ψ angle between projected baseline circle ipb12 ipb23 ipb31 = b12 sin θ12e + b23 sin θ23e + b31 sin θ31e and star meridian =0. (C6) z zenith angle Z baseline angle seen from the Earth cen- These equations remain correct if transformed by com- ter plex conjugation or multiplied by a complex phase factor;

[1] Abramowitz, M., and I. A. Stegun (eds.), 1972, Handbook Proc. SPIE, pp. 1540–1548. of Mathematical Functions (Dover Publications, New [6] Hanisch, R. J., A. Farris, E. W. Greisen, W. D. Pence, York), 9th edition, ISBN 0-486-61272-4. B. M. Schlesinger, P. J. Teuben, R. W. Thompson, and [2] Angeli, G. Z., B. C. Fitz-Patrick, and M. Lloyd-Hart, A. W. III, 2001, Astron. & Astrophys. 376(1), 359. 2000, in Telescope Structures, Enclosures, Controls, [7] Karttunen, H., P. Kröger, H. Oja, M. Poutanen, and K. J. Assembly/Integration/Validation, and Commissioning, Donner (eds.), 1987, Fundamental Astronomy (Springer, edited by T. A. Sebring and T. Andersen (Int. Soc. Op- Berlin, Heidelberg), ISBN 0-387-17264-5. tical Engineering), volume 4004 of Proc. SPIE, pp. 202– [8] Köhler, R., 2006, priv. commun. 211. [9] Lang, K. R., 1998, Astrophysical Formulae, volume II of [3] Calabretta, M. R., and E. W. Greisen, 2002, Astron. & Astronomy and Astrophysics Library (Springer Verlag, Astrophys. 395, 1077. Berlin, Heidelberg), 3rd edition, ISBN 3-540–64664-7, E: [4] Carbillet, M., S. Correia, P. Boccacci, and M. Bertero, In the second line of (5.45), sin A should read cos A. 2002, Astron. Astrophys. 387(3), 744. [10] Leinert, C., 2006, priv. commun. [5] Galliano, E., M. Schöller, M. Fischer, C. Hummel, [11] Makovoz, D., 2004, Publ. Astron. Soc. Pac. 116(824), S. Morel, F. Rantakyrö, and M. Vannier, 2004, in New 971. Frontiers in Stellar Interferometry, edited by W. A. [12] Mathar, R. J., 1999, arXiv:physics/9907051 . Traub (Int. Soc. Optical Engineering), volume 5491 of [13] Mathar, R. J., 2005, Baltic Astronomy 14(2), 277. 10

[14] National Imagery and Mapping Agency, 2000, partment, California Institute of Technology, URL Department of Defense World Geodetic System http://www.astro.caltech.edu/~tjp/citvlb/. 1984, Technical Report TR8350.2, NIMA, URL [18] Pettauer, T., 1994, Solar Physics 155(1), 195. http://earth-info.nga.mil/GandG/publications/tr8350.2/tr8350_2.html[19] Roe, H. G., 2002,. Publ. Astron. Soc. Pac. 114, 450. [15] Pauls, T. A., J. S. Young, W. D. Cotton, and J. D. Mon- [20] Smart, W. M., 1949, Text-book on Spherical Astronomy nier, 2004, in New Frontiers in Stellar Interferometry, (Cambridge University Press, Cambridge), 4th edition. edited by W. A. Traub (Int. Soc. Optical Engineering), [21] Vondr´ak, J., 2004, Serb. Astron. J. 168, 1. volume 5491 of Proc. SPIE, pp. 1231–1239. [22] Wallace, P. T., 2003, SLALIB—Positional As- [16] Pauls, T. A., J. S. Young, W. D. Cotton, and J. D. Mon- tronomy Library 2.4–12 Programmer’s Manual, nier, 2005, Publ. Astron. Soc. Pac. 117(837), 1255. http://www.starlink.rl.ac.uk/static www/soft further SLALIB.html. [17] Pearson, T. J., 1991, Introduction to the Caltech VLBI programs, Technical Report, Astronomy De-