No. 6, 2006KINEMATIC DISTANCE ERRORS 2379Fig. 4.—(a) Rotation curve given by <strong>the</strong> background potential (solid l<strong>in</strong>e) compared with <strong>the</strong> measured rotation curve. For <strong>the</strong> <strong>in</strong>ner <strong>galaxy</strong>, <strong>the</strong> rotations measuredfor both positive (dotted l<strong>in</strong>e) <strong>and</strong> negative longitudes (dashed l<strong>in</strong>e) are presented. The rotation curve correspond<strong>in</strong>g to <strong>the</strong> outer <strong>galaxy</strong> is also shown (dash-dotted l<strong>in</strong>e).(b) Rotation curve given by <strong>the</strong> background potential (solid l<strong>in</strong>e) compared with <strong>the</strong> average <strong>of</strong> <strong>the</strong> north <strong>and</strong> south rotation curves (dotted l<strong>in</strong>e) <strong>and</strong> <strong>the</strong> mean azimuthalvelocity <strong>of</strong> <strong>the</strong> gas <strong>in</strong> <strong>the</strong> simulation (dashed l<strong>in</strong>e). Note that both <strong>the</strong> mean velocity <strong>and</strong> <strong>the</strong> mean rotation curve are above <strong>the</strong> background rotation curve for most <strong>of</strong> <strong>the</strong>radial doma<strong>in</strong>.side <strong>of</strong> <strong>the</strong> tangent po<strong>in</strong>t by look<strong>in</strong>g at <strong>the</strong> galactic latitude extension<strong>of</strong> <strong>the</strong> s<strong>our</strong>ce (Fish et al. 2003), or by us<strong>in</strong>g observed <strong>in</strong>termediateabsorption features (Watson et al. 2003; Sewilo et al.2004). For this <strong>in</strong>vestigation I decided to cheat: I looked upwhich side <strong>of</strong> <strong>the</strong> tangent po<strong>in</strong>t <strong>the</strong> gas parcel fell on, <strong>and</strong> chose<strong>the</strong> measured distance accord<strong>in</strong>gly.Figure 5a shows <strong>the</strong> error <strong>in</strong> measured distance with respect to<strong>the</strong> real distance <strong>in</strong> <strong>the</strong> model. Recall<strong>in</strong>g Figure 4, <strong>the</strong> observerwould determ<strong>in</strong>e different rotation curves for <strong>the</strong> nor<strong>the</strong>rn <strong>and</strong>sou<strong>the</strong>rn sides <strong>of</strong> <strong>the</strong> galactic center. Accord<strong>in</strong>gly, <strong>in</strong> determ<strong>in</strong><strong>in</strong>g<strong>the</strong> k<strong>in</strong>ematic distance for Figure 5a, <strong>the</strong> rotation curve usedis that <strong>of</strong> <strong>the</strong> correspond<strong>in</strong>g side <strong>of</strong> <strong>the</strong> galactic center. It isFig. 5.—Error <strong>in</strong> <strong>the</strong> measured k<strong>in</strong>ematic distance (d ) obta<strong>in</strong>ed under <strong>the</strong> assumption <strong>of</strong> circular orbits follow<strong>in</strong>g (a) <strong>the</strong> measured rotation curve <strong>and</strong> (b) <strong>the</strong>rotation given by <strong>the</strong> background potential. Sections <strong>of</strong> 7 around <strong>the</strong> galactic longitudes l ¼ 0 <strong>and</strong> 180 were excluded. Although <strong>the</strong> error <strong>in</strong> most <strong>of</strong> <strong>the</strong> galactic diskis <strong>of</strong> <strong>the</strong> order <strong>of</strong> 0:5 kpc, it is significantly larger at <strong>the</strong> positions <strong>of</strong> <strong>the</strong> spiral arms. The sharp edges at <strong>the</strong> positions <strong>of</strong> <strong>the</strong> tangent po<strong>in</strong>ts are a consequence <strong>of</strong> <strong>the</strong> factthat <strong>the</strong> term<strong>in</strong>al velocities do not occur at those po<strong>in</strong>ts. The <strong>errors</strong> <strong>in</strong> measured k<strong>in</strong>ematic <strong>distances</strong> are larger when <strong>the</strong> real ( background) rotation curve is used. [See <strong>the</strong>electronic edition <strong>of</strong> <strong>the</strong> J<strong>our</strong>nal for a color version <strong>of</strong> this figure.]
2380 GÓMEZVol. 1321 In order to dim<strong>in</strong>ish spurious <strong>in</strong>terpolation effects, each gas parcel was spreadus<strong>in</strong>g a two-dimensional Gaussian weight function <strong>in</strong>to a 3 ; 3 grid-cell regionaround <strong>the</strong> position correspond<strong>in</strong>g to that parcel’s galactic longitude <strong>and</strong> measureddistance.Fig. 6.—Remapp<strong>in</strong>g <strong>of</strong> <strong>the</strong> gas distribution result<strong>in</strong>g from <strong>the</strong> k<strong>in</strong>ematic <strong>distances</strong>us<strong>in</strong>g <strong>the</strong> measured rotation curves <strong>in</strong> Fig. 4a <strong>and</strong> assum<strong>in</strong>g circular gasorbits. Note <strong>the</strong> regions near <strong>the</strong> tangent po<strong>in</strong>ts <strong>and</strong> <strong>the</strong> corotation circle, wherelittle or no gas is mapped to.noticeable that although <strong>the</strong> <strong>errors</strong> are on <strong>the</strong> order <strong>of</strong> 0:5 kpc<strong>in</strong>most <strong>of</strong> <strong>the</strong> galactic disk, <strong>the</strong>y are significantly larger at <strong>the</strong>positions <strong>of</strong> <strong>the</strong> spiral arms (as h<strong>in</strong>ted by Gómez & Cox 2004b).This fact has a special impact <strong>in</strong> studies <strong>of</strong> <strong>the</strong> spiral structure <strong>of</strong><strong>the</strong> Galaxy that rely on k<strong>in</strong>ematic <strong>distances</strong>, s<strong>in</strong>ce it distorts <strong>the</strong><strong>image</strong> <strong>the</strong> observer would generate (see x 4.1).There is ano<strong>the</strong>r significant feature <strong>in</strong> Figure 5. Although <strong>the</strong>term<strong>in</strong>al velocity does not really arise from <strong>the</strong> tangent po<strong>in</strong>t, <strong>the</strong>circular orbits assumption assigns gas observed near term<strong>in</strong>al velocityto that po<strong>in</strong>t. This generates a feature <strong>in</strong> <strong>the</strong> <strong>errors</strong> thatcorresponds to <strong>the</strong> locus <strong>of</strong> <strong>the</strong> tangent po<strong>in</strong>ts. Aga<strong>in</strong>, <strong>the</strong> error issignificant at <strong>the</strong> position <strong>of</strong> <strong>the</strong> spiral arms <strong>and</strong> would generatelarge <strong>errors</strong> <strong>in</strong> <strong>the</strong> determ<strong>in</strong>ation <strong>of</strong> <strong>distances</strong> to objects that trace<strong>the</strong> spiral structure.The assumptions <strong>of</strong> circular orbits <strong>and</strong> different rotation curvesfor positive <strong>and</strong> negative longitudes are, <strong>of</strong> c<strong>our</strong>se, <strong>in</strong>consistent.One solution is to fit a s<strong>in</strong>gle rotation curve to both sides <strong>of</strong> <strong>the</strong>Galaxy. In order to test this method <strong>the</strong> average <strong>of</strong> both rotationcurves was taken, <strong>and</strong> <strong>the</strong> equivalent <strong>of</strong> Figure 5a was calculated.The result was that <strong>the</strong> magnitude <strong>of</strong> <strong>the</strong> error <strong>in</strong> <strong>the</strong> k<strong>in</strong>ematic<strong>distances</strong> was approximately <strong>the</strong> same, but <strong>the</strong> area with error>0.5 kpc spanned a larger fraction <strong>of</strong> <strong>the</strong> disk.Suppose now that <strong>the</strong> imag<strong>in</strong>ary observer somehow managesto obta<strong>in</strong> <strong>the</strong> large-scale distribution <strong>of</strong> stellar mass <strong>in</strong> <strong>the</strong> model<strong>galaxy</strong>. This would allow <strong>the</strong> derivation <strong>of</strong> <strong>the</strong> real rotation curvefrom <strong>the</strong> background axisymmetric potential. If <strong>the</strong> observer nowuses that real rotation to determ<strong>in</strong>e k<strong>in</strong>ematic <strong>distances</strong>, evenlarger distance <strong>errors</strong> would be obta<strong>in</strong>ed, especially for <strong>the</strong> <strong>in</strong>ner<strong>galaxy</strong>, as shown <strong>in</strong> Figure 5b. This counter<strong>in</strong>tuitive result arisesbecause, at this po<strong>in</strong>t <strong>of</strong> <strong>the</strong> simulation, <strong>the</strong> gas has already adoptedorbits that are not only <strong>in</strong>fluenced by <strong>the</strong> background potentialbut also by <strong>the</strong> large-scale magnetic field (likely different from<strong>the</strong> field <strong>in</strong> <strong>the</strong> <strong>in</strong>itial conditions) <strong>and</strong> <strong>the</strong> torques <strong>and</strong> resonancesgenerated by <strong>the</strong> spiral perturbation. Although <strong>the</strong> real rotationcurve is consistent with <strong>the</strong> most important determ<strong>in</strong>ant <strong>of</strong> <strong>the</strong>gas rotation velocity (<strong>the</strong> background mass distribution), it doesnot <strong>in</strong>clude o<strong>the</strong>r <strong>in</strong>fluences <strong>in</strong> that velocity, while <strong>the</strong> ‘‘wrong’’rotation curve determ<strong>in</strong>ed from gaseous term<strong>in</strong>al velocities moreclosely reflects <strong>the</strong> real motion <strong>of</strong> <strong>the</strong> gas (recall Fig. 4, <strong>in</strong> which<strong>the</strong> measured rotation curve is systematically above <strong>the</strong> truerotation).Although <strong>in</strong>tr<strong>in</strong>sically <strong>in</strong>consistent, <strong>the</strong> two different measuredrotation curves are used <strong>in</strong> <strong>the</strong> rema<strong>in</strong>der <strong>of</strong> this <strong>in</strong>vestigation s<strong>in</strong>cethat procedure leads to smaller distance <strong>errors</strong>. The results presented<strong>in</strong> x 4.1 are even more noteworthy if <strong>the</strong> average or <strong>the</strong> realrotation curves are used.4.1. The Galaxy DistortedConsider now that <strong>the</strong> imag<strong>in</strong>ary observer is try<strong>in</strong>g to study<strong>the</strong> spiral structure <strong>of</strong> <strong>the</strong> <strong>galaxy</strong> he/she lives <strong>in</strong>. The procedurewould consist <strong>of</strong> translat<strong>in</strong>g <strong>the</strong> longitude-velocity data obta<strong>in</strong>edfrom a diffuse gas survey, for example, <strong>in</strong>to a spatial distributionus<strong>in</strong>g <strong>the</strong> k<strong>in</strong>ematic <strong>distances</strong> that result from <strong>the</strong> assumption <strong>of</strong>circular orbits that follow <strong>the</strong> measured rotation curve. 1 The result<strong>in</strong>gmap is shown <strong>in</strong> Figure 6. Note that <strong>the</strong> features describedfor Figure 1 (namely, <strong>the</strong> f<strong>our</strong> spiral arms, <strong>the</strong> 4 kpc high-densityr<strong>in</strong>g, <strong>and</strong> <strong>the</strong> corotation low-density r<strong>in</strong>g) all but disappear, whilenew fictitious features, such as <strong>the</strong> structure <strong>in</strong> <strong>the</strong> outer Galaxy,are formed as a consequence <strong>of</strong> <strong>the</strong> oscillations <strong>in</strong> <strong>the</strong> outer rotationcurve. Also significant <strong>in</strong> this figure are <strong>the</strong> regions wherelittle or no gas is assigned by <strong>the</strong> mapp<strong>in</strong>g, namely, <strong>the</strong> b<strong>and</strong>snear <strong>the</strong> corotation circle <strong>and</strong> <strong>the</strong> quasi-triangular regions near<strong>the</strong> tangent po<strong>in</strong>t locus. (These nearly empty regions are significantlylarger when <strong>the</strong> background or <strong>the</strong> mean rotation curvesare used to determ<strong>in</strong>e <strong>the</strong> distance to <strong>the</strong> observed gas parcel.)The imag<strong>in</strong>ary observer would likely conclude that his/herhome <strong>galaxy</strong> has two ill-def<strong>in</strong>ed spiral arms. If a logarithmic spiralmodel were forced, an 11 pitch angle <strong>and</strong> a density contrastmuch stronger than that <strong>in</strong> <strong>the</strong> numerical model would be found.Ano<strong>the</strong>r possibility for determ<strong>in</strong><strong>in</strong>g <strong>the</strong> distance to a gas parcelconsists <strong>of</strong> compar<strong>in</strong>g <strong>the</strong> l<strong>in</strong>e-<strong>of</strong>-sight velocity <strong>of</strong> <strong>the</strong> parcelwith <strong>the</strong> predicted velocity obta<strong>in</strong>ed from some model for <strong>the</strong>Galactic structure. For <strong>the</strong> numerical model described <strong>in</strong> x 2, givena Galactic longitude, Figure 3 is searched for <strong>the</strong> required velocity,<strong>and</strong> <strong>the</strong> correspond<strong>in</strong>g distance is read out. 2 Although<strong>the</strong> same procedure to solve <strong>the</strong> ambiguity with respect to <strong>the</strong>tangent po<strong>in</strong>t is used, <strong>the</strong> noncircular motions <strong>in</strong>troduce new distanceambiguities for certa<strong>in</strong> longitude-velocity values (up to 11,although 3 is a more typical number). When <strong>the</strong>se ambiguitiesappear, <strong>the</strong>y occur close to each o<strong>the</strong>r, mak<strong>in</strong>g <strong>the</strong>ir resolutiondifficult. So, when reconstruct<strong>in</strong>g <strong>the</strong> map <strong>of</strong> <strong>the</strong> <strong>galaxy</strong>, <strong>the</strong> gasdensity is equally split among <strong>the</strong>se positions.The result is shown <strong>in</strong> Figure 7. The new distance ambiguitiesstill <strong>in</strong>troduce spurious structure, such as <strong>the</strong> splitt<strong>in</strong>g <strong>of</strong> <strong>the</strong> spiralarms. Never<strong>the</strong>less, <strong>the</strong> number <strong>and</strong> position <strong>of</strong> <strong>the</strong> arms, <strong>the</strong>structure around <strong>the</strong> corotation radius, <strong>and</strong> <strong>the</strong> lack <strong>of</strong> features <strong>in</strong><strong>the</strong> outer <strong>galaxy</strong> are recovered. The imag<strong>in</strong>ary observer wouldlikely conclude that his/her home <strong>galaxy</strong> has f<strong>our</strong> arms with 9 <strong>and</strong> 12N5 pitch angles, although he/she would also f<strong>in</strong>d nonexistentbridges <strong>and</strong> spurs. On <strong>the</strong> o<strong>the</strong>r h<strong>and</strong>, it should be consideredthat <strong>the</strong> imag<strong>in</strong>ary observer would not see <strong>the</strong>rmal or turbulent2 A simple C language program that provides a distance given a Galactic longitude<strong>and</strong> l<strong>in</strong>e-<strong>of</strong>-sight velocity value <strong>and</strong> uncerta<strong>in</strong>ty is available at http://www.astrosmo.unam.mx/~g.gomez/publica/. In that program <strong>the</strong> result<strong>in</strong>g distanceis given as a range <strong>in</strong>stead <strong>of</strong> a central value <strong>and</strong> uncerta<strong>in</strong>ty, s<strong>in</strong>ce <strong>the</strong> velocitydistancemapp<strong>in</strong>g makes <strong>the</strong> distance probability distribution nei<strong>the</strong>r uniform <strong>in</strong><strong>the</strong> range nor peaked around a central value.