Arxiv:Astro-Ph/0608273 V1 13 Aug 2006 Ln Oteotclpt Ea Ocuetemanuscript the Conclude IV
Total Page:16
File Type:pdf, Size:1020Kb
CORE Metadata, citation and similar papers at core.ac.uk Provided by CERN Document Server Spherical Trigonometry of the Projected Baseline Angle Richard J. Mathar∗ Leiden Observatory, P.O. Box 9513, 2300 RA Leiden, The Netherlands (Dated: August 15, 2006) The basic geometry of a stellar interferometer with two telescopes consists of a baseline vector and a direction to a star. Two derived vectors are the delay vector, and the projected baseline vector in the plane of the wavefronts of the stellar light. The manuscript deals with the trigonometry of projecting the baseline further outwards onto the celestial sphere. The position angle of the projected baseline is defined, measured in a plane tangential to the celestial sphere, tangent point at the position of the star. This angle represents two orthogonal directions on the sky, differential star positions which are aligned with or orthogonal to the gradient of the delay recorded in the u−v plane. The North Celestial Pole is chosen as the reference direction of the projected baseline angle, adapted to the common definition of the “parallactic” angle. PACS numbers: 95.10.Jk, 95.75.Kk Keywords: projected baseline; parallactic angle; position angle; celestial sphere; spherical astronomy; stellar interferometer I. SCOPE II. SITE GEOMETRY A. Telescope Positions The paper describes a standard to define the plane that 1. Spherical Approximation contains the baseline of a stellar interferometer and the star, and its intersection with the Celestial Sphere. This In a geocentric coordinate system, the Cartesian coor- intersection is a great circle, and its section between the dinates of the two telescopes of a stellar interferometer star and the point where the stretched baseline meets the are related to the geographic longitude λi, latitude Φi Celestial Sphere defines an oriented projected baseline. and effective Earth radius ρ (sum of the Earth radius The incentive to use this projection onto the Celestial and altitude above sea level), Sphere is that this ties the baseline to sky coordinates, which defines directions in equatorial coordinates which cos λi cosΦi T are independent of the observatory’s location on Earth. i ≡ ρ sin λi cosΦi , i =1, 2. (1) sinΦ This is (i) the introduction of a polar coordinate system’s i g polar angle in the u − v coordinate system associated with the visibility tables of the Optical Interferometry We add a “g” to the geocentric coordinates to set them Exchange Format (OIFITS) [15, 16], detailing the or- apart from other Cartesian coordinates that will be used thogonal directions which are well and poorly resolved further down. A great circle of radius ρ ≈ 6380 km, cen- by the interferometer, and (ii) a definition of celestial tered at the Earth center, connects T1 and T2; a vector Jˆ 1 T T longitudes in a coordinate system where the pivot point perpendicular to this circle is ≡ ρ2 sin Z 1 × 2, where has been relocated from the North Celestial Pole (NCP) a baseline aperture angle Z has been defined under which arXiv:astro-ph/0608273 v1 13 Aug 2006 to the star. the baseline vector b = T2 − T1 (2) In overview, Section II introduces the standard set is seen from the Earth center: of variables in the familiar polar coordinate systems. Section III defines and computes the baseline position T1 · T2 = |T1| |T2| cos Z; (3) angle—which can be done purely in equatorial or hori- zontal coordinates or in a hybrid way if the parallactic cos Z = cosΦ1 cosΦ2 cos(λ2 − λ1)+sinΦ1 sinΦ2(4). angle is used to bridge these. Some comments on relating The wind rose direction from T1 to T2 is computed as distances away from the optical axis (star) in the tangent follows: define a tangent plane to the Earth at T1. Two plane to the optical path delay conclude the manuscript orthonormal vectors that span this plane are in Section IV. − cos λ1 sinΦ1 Nˆ 1 ≡ − sin λ1 sinΦ1 ∼ ∂T1/∂Φ1; |Nˆ 1| =1. ∗ cosΦ1 Electronic address: [email protected]; g URL: http://www.strw.leidenuniv.nl/~mathar (5) 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 defined 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 defines 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 <A<π<Ab and D < 0. t t We label this coordinate system t as “topocentric,” with (Note that the tables in [5] use a different definition.) the first component North, the second component West Straight forward algebra establishes τ from and the third component up. The value of A used in this script picks one of the (countably many) conventions: cosΦ1 sinΦ2 − sinΦ1 cosΦ2 cos(λ1 − λ2) cos τ = ;(7) South means A = 0 and West means A = +π/2. This sin Z conventional definition of the horizontal means the star cosΦ2 sin(λ2 − λ1) coordinates can be transformed to the equatorial param- sin τ = . (8) sin Z eters of hour angle h = l − α, declination δ, and right ascension α [9, (5.45)][7, (2.13)], 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.