الفلك

رسم بيانات السرعة الشعاعية GAIA barycentric باستخدام الثعبان

رسم بيانات السرعة الشعاعية GAIA barycentric باستخدام الثعبان


We are searching data for your request:

Forums and discussions:
Manuals and reference books:
Data from registers:
Wait the end of the search in all databases.
Upon completion, a link will appear to access the found materials.

هل يعرف أي شخص درسًا تعليميًا جيدًا حول كيفية استخدام Python و Astropy لإنشاء قطع من GAIA أو RAVE أو أي قاعدة بيانات أخرى لبيانات السرعة الشعاعية المركزية (مركزية الشمس) للنجوم في مجرتنا؟

أرغب في إنشاء مخططات مختلفة لسرعات مركزية الشمس للنجوم في مجرتنا ، مع وضع قيود مختلفة على المسافات من الشمس ، وأود أن أرى كيف تبدو البيانات مع تصحيح الحركة المناسبة للشمس وبدون . الآن ، أعلم أنه بإمكاني العثور على قطع أرض قام بها الآخرون ، لكنني أريد حقًا تعلم كيفية إنشاء قطع الأرض الخاصة بي.

أنا جديد تمامًا على بايثون ، لذا يرجى استهداف الإجابات وفقًا لذلك ...

(تمت إضافته لاحقًا :) من المحتمل أن تكون هناك وظيفة فلكية بسيطة إذا استدعيتها بالمتغيرات الصحيحة ، على سبيل المثال ، شريحة معينة من السماء في RA و Dec ، ونطاق مسافة معين في kpc ، فستعيد قائمة بكل شيء النجوم في تلك المنطقة مع سرعاتها الشعاعية. أنا فقط بحاجة لمعرفة اسم هذه الوظيفة وكيفية تسميتها.


يشار إلى المجموعات الفرعية ([أ] ، [ب]):

  • جون بارنز (OU) [A]
  • إيوانا سيوكا (MSSL ، UCL) (طالب دكتوراه) [ب]
  • توم كروسلاند (UCL) (طالب دكتوراه) [A]
  • جراهام فليمنج (UCLan) (طالب ماجستير) [A]
  • كيث إينايت ، ليف جي إم (ماجستير) [أ]
  • كاتريونا موراي (IoA ، كامبريدج) [B]
  • جون موريل (عالم فلك هواة) [ب]
  • هايدي ثيمان (OU) [A]
  • يوجين فاسيليف (IoA ، كامبريدج) [ب]
  • أندرو ويلسون (إكستر) (طالب دكتوراه) [ب]
  • سيلفيا زين (MSSL ، UCL) [B]
  • .

مدرسون (المجموعات أ ، ب)

  • جيورجيا بوسو (IoA ، كامبريدج) [B]
  • نايجل هامبلي (IfA ، إدنبرة) [ب]
  • فلور فان ليوين (IoA ، كامبريدج) [A]
  • جورج سيبروك (MSSL ، UCL) [الاتصال الهاتفي عن بعد]
  • مارك تايلور (بريستول) [أ]
  • نيكولاس والتون (IoA ، كامبريدج) [A / B]

الحاضرين

  • شداب علم (IfA ، إدنبرة)
  • ديفيد بيرندت (IfA ، إدنبرة) [طالب دكتوراه]
  • بيث بيلر (IfA ، إدنبرة)
  • آرون برادلي (IfA ، إدنبرة) [طالب دكتوراه]
  • أبي تشون (باث) [طالب دكتوراه]
  • رفائيل إراني (IfA ، إدنبرة) [طالب دكتوراه]
  • أنيت فيرجسون (IfA ، إدنبرة)
  • أوليفيا جونز (UKATC ، إدنبرة)
  • باميلا كلاسن (UKATC ، إدنبرة)
  • بيت كوزما (IfA ، إدنبرة)
  • ماركو لام (ARI ، LJMU)
  • بول روس ماكويرتر (IfA ، إدنبرة) [طالب دكتوراه]
  • كيث بالين (IfA ، إدنبرة)
  • تاتيانا بافليدو (سانت أندروز) [طالبة دكتوراه]
  • روبن سانشيز يانسن (UKATC ، إدنبرة)
  • هيرمين شنيتلر (IfA ، إدنبرة)
  • آنا ليزا فاري (IfA ، إدنبرة)
  • روزماري وايز (JHU)

مدرسون

  • نايجل هامبلي (IfA ، إدنبرة)
  • فلور فان ليوين (IoA ، كامبريدج)
  • نيك رويل (IfA ، إدنبرة)
  • جورج سيبروك (MSSL ، UCL)
  • مارك تايلور (بريستول)
  • نيكولاس والتون (IoA ، كامبريدج)

دعونا نراجع كيف وصلنا إلى هذه النقطة.

قمنا بعمل استعلام ADQL لخادم Gaia للحصول على بيانات عن النجوم بالقرب من GD-1.

قمنا بتحويل الإحداثيات إلى إطار GD1Koposov10 حتى نتمكن من تحديد النجوم على طول الخط المركزي لـ GD-1.

لقد رسمنا الحركة المناسبة لنجوم خط الوسط لتحديد حدود منطقة الكثافة الزائدة.

لقد صنعنا قناعًا يختار النجوم التي تكون حركتها الصحيحة في منطقة الكثافة الزائدة.

في هذه المرحلة ، قمنا بتنزيل البيانات لعدد كبير نسبيًا من النجوم (أكثر من 100000) واخترنا عددًا صغيرًا نسبيًا (حوالي 1000).

سيكون من الأفضل استخدام ADQL لاختيار النجوم التي نحتاجها فقط. ومن شأن ذلك أيضًا أن يجعل من الممكن تنزيل البيانات التي تغطي مساحة أكبر من السماء.

ومع ذلك ، فإن التحديد الذي قمنا به كان يعتمد على الحركة المناسبة في إطار GD1Koposov10. من أجل القيام بالاختيار نفسه في ADQL ، يتعين علينا العمل مع الحركات المناسبة في ICRS.

للتذكير ، إليك المستطيل الذي اخترناه بناءً على الحركة المناسبة في إطار GD1Koposov10.


رسم بيانات السرعة الشعاعية GAIA barycentric باستخدام بيثون - علم الفلك

الوظائف المتعلقة بالقياس الفلكي مناسبة للاستخدام مع البيانات من مهمة قياس الفلك Gaia.

الأساليب هنا ليست خاصة بمهمة Gaia ، ولكن يتم تحديد معلمات الوظائف ووحداتها في شكل مناسب للاستخدام مع بيانات Gaia ، ولا سيما كتالوج gaia_source المتاح من http: //gea.esac.esa .int / أرشيف / ونسخ أو مرايا.

  • الموضع وحساب متجه السرعة والتلاعب به
  • تقدير المسافة من المنظر
  • انتشار الفلك إلى عصور مختلفة

متجهات الموقف والسرعة

يتم توفير وظائف لتحويل المعلمات الفلكية الموجودة في كتالوج Gaia إلى متجهات الموضع الديكارتي (XYZ) والسرعة (UVW) الخاصة بـ ICRS. يتم توفير وظائف أيضًا لتحويل هذه المتجهات بين إحداثيات ICRS و Galactic أو Ecliptic. الحسابات بسيطة إلى حد ما ، وتتبع المعادلات الموضحة في القسم 1.5.6 من كتالوجات Hipparcos و Tycho، ESA SP-1200 (1997) وكذلك القسم 3.1.7 من وثائق Gaia DR2 (2018).

غالبًا ما يتم دمج هذه الوظائف على سبيل المثال لحساب الموقع والسرعة في إحداثيات المجرة من قيم كتالوج Gaia ، قد تكون التعبيرات التالية مفيدة: على الرغم من ملاحظة أن هذه الأمثلة المعينة تعكس ببساطة اختلاف المنظر لتوفير تقديرات المسافة ، وهو أمر غير صحيح بشكل عام. لاحظ أيضًا أن هذه الوظائف لا تحاول تصحيح الحركة الشمسية. يجب إجراء هذه التعديلات يدويًا على نتائج هذه الوظائف إذا كانت مطلوبة.

لا تتوفر حاليًا وظائف لحساب الأخطاء على المكونات الديكارتية بناءً على الخطأ وكميات الارتباط من كتالوج Gaia. سوف تتطلب دعوات معقدة إلى حد ما. إذا كان هناك طلب فيمكن تنفيذها في المستقبل.

تقدير المسافة

  • سي أل بيلر جونز ، "تقدير المسافات من المنظر"، PASP 127، الصفحة 994 (2015) 2015PASP.127..994B
  • T.L. Astraatmadja و CL.Bailer-Jones ، "تقدير المسافات من المنظر. II. أداء مقدرات المسافة بايزي على كتالوج يشبه Gaia"، ApJ 832، a137 (2016) 2016ApJ. 832.137 أ
  • X.Luri et al. ، "إصدار بيانات Gaia 2: استخدام Gaia Parallaxes"، A & ampA في الصحافة (2018) arXiv: 1804.09376

تتوافق الوظائف المتوفرة هنا مع الحسابات من Astraatmadja & amp Bailer-Jones ، "تقدير المسافات من المنظر. ثالثا. مسافات مليوني نجمة في كتالوج Gaia DR1"، ApJ 833، a119 (2016) 2016ApJ. 833..119A على أساس تناقص كثافة الفضاء بشكل كبير محددة مسبقا فيه. تمت كتابة هذا التطبيق بالإشارة إلى تنفيذ Java بواسطة Enrique Utrilla (DPAC).

يتم تحديد معلمات هذه الوظائف من خلال مقياس الطول إل يحدد الانحلال الأسي (يكون وضع ملف PDF السابق في ص=2إل). يجب توفير بعض القيم لمقياس الطول هذا ، المحدد في الفرسخ ، للوظائف كمعامل lpc.

ملحوظة أن القيم التي توفرها هذه الوظائف تفعل ليس تطابق تلك الموجودة في الورقة Bailer-Jones et al. "تقدير المسافات من Parallaxes IV: المسافات إلى 1.33 مليار نجمة في Gaia Data Release 2 "، مقبول لـ AJ (2018) arXiv: 1804.10121. تختلف حسابات تلك الورقة عن تلك المعروضة هنا بعدة طرق: فهي تستخدم نموذجًا مجريًا لمقياس الطول المعتمد على الاتجاه غير المتاح حاليًا هنا ، وتطبق مسبقًا تصحيحًا للاختلاف المنظر يبلغ -0.029mas ، وتستخدم مقاييس عدم يقين مختلفة وفي بعض الحالات (ملف PDF ثنائي النسق) أفضل مقدر مختلف للمسافة.

انتشار العصر

يوفر كتالوج مصدر Gaia ، لبعض المصادر على الأقل ، الحل الفلكي المكون من ستة معلمات (الصعود الصحيح ، والانحراف ، والمنظر ، والحركة المناسبة في RA و Dec ، والسرعة الشعاعية) ، إلى جانب أخطاء في هذه القيم والارتباطات بين هذه الأخطاء. في حين أن تقديرًا أوليًا للموضع في فترة سابقة أو لاحقة من القياس يمكن إجراؤه بضرب مكونات الحركة المناسبة بفارق الحقبة والإضافة إلى الموضع المقاس ، فإن المعالجة الأكثر دقة مطلوبة للانتشار الدقيق بين فترات المعلمات الفلكية ، وإذا لزم الأمر أخطائها وارتباطاتها. تم توضيح التعبيرات الخاصة بذلك في القسم 1.5.5 (المجلد 1) من كتالوجات Hipparcos و Tycho، ESA SP-1200 (1997) (لكن انظر أدناه) ، ويستند الكود إلى تطبيق قام به Alexey Butkevich و Daniel Michalik (DPAC). يتم تطبيق التصحيح على معالجة SP-1200 لحالة عدم اليقين في السرعة الشعاعية التالية ميشاليك وآخرون 2014 2014 أمبير وأمبير. 571A.85M بسبب تعاملهم الأفضل مع السرعات الشعاعية الصغيرة أو المنظر.

تعطي الحسابات نفس النتائج ، وإن لم تكن بنفس الشكل بالضبط ، مثل وظائف انتشار العصر المتاحة في خدمة أرشيف Gaia.

بولاركسيز (فاي ، ثيتا ، ص) يحول من الإحداثيات القطبية الكروية إلى الإحداثيات الديكارتية.

  • المعلمات:
    • فاي (النقطة العائمة): خط الطول بالدرجات
    • ثيتا (النقطة العائمة): خط العرض بالدرجات
    • ص (النقطة العائمة): مسافة شعاعية
    • (صفيف النقطة العائمة): متجه ثلاثي العناصر يعطي الإحداثيات الديكارتية
    • polarXYZ (ra ، dec ، المسافة_estimate)
    • polarXYZ (l، b، 3262./parallax) - يحسب مكونات المتجه بوحدات السنة الضوئية في نظام المجرة ، على افتراض أن المسافة هي معكوس المنظر

    astromXYZ (ra، dec، parallax) تحسب المكونات الديكارتية للموضع من RA ، الانحراف والمنظر. هذه وظيفة ملائمة ، تعادل:

    لاحظ أن هذا يؤدي إلى قياس المسافة باستخدام انعكاس بسيط في اختلاف المنظر ، والذي لا يمكن الاعتماد عليه بشكل عام في حالات اختلاف المنظر مع أخطاء غير قابلة للإهمال. استخدام على مسؤوليتك الخاصة.

    • المعلمات:
      • را (النقطة العائمة): حق الصعود بالدرجات
      • ديسمبر (النقطة العائمة): الانحراف بالدرجات
      • المنظر (النقطة العائمة): اختلاف المنظر في ماس
      • (صفيف النقطة العائمة): متجه ثلاثي العناصر يعطي إحداثيات الفراغ الاستوائي بالفرسخ
      • astromXYZ (ra، dec، parallax)
      • icrsToGal (astromXYZ (ra، dec، parallax))

      icrsToGal (xyz) يحول متجهًا ثلاثي العناصر يمثل إحداثيات ICRS (الاستوائية) إلى إحداثيات مجرية. يمكن استخدام هذا مع متجهات الموضع أو السرعة.

      يتم ضرب متجه الإدخال بالمصفوفة أجي'، في المعادلة. 3.61 من وثائق Gaia DR2 ، بعد Eq. 1.5.13 من كتالوج Hipparcos.

      نظام إحداثيات الخرج يمين ، مع وجود ثلاثة مكونات موجبة في اتجاهات مركز المجرة ، ودوران المجرة ، وقطب المجرة الشمالي على التوالي.

      • المعلمات:
        • xyz (صفيف النقطة العائمة): متجه ثلاثي العناصر يعطي مكونات ICRS الديكارتية
        • (مجموعة من الفاصلة العائمة): ناقل ثلاثي العناصر يعطي مكونات مجرة ​​ديكارتية
        • icrsToGal (بولاركسيز (را ، ديسمبر ، مسافة))

        galToIcrs (xyz) يحول متجهًا ثلاثي العناصر يمثل إحداثيات المجرة إلى إحداثيات ICRS (الاستوائية). يمكن استخدام هذا مع متجهات الموضع أو السرعة.

        يتم ضرب متجه الإدخال بالمصفوفة أجي، في المعادلة. 3.61 من وثائق Gaia DR2 ، بعد Eq. 1.5.13 من كتالوج Hipparcos.

        نظام إحداثيات الإدخال يمين ، مع وجود ثلاثة مكونات موجبة في اتجاهات مركز المجرة ، ودوران المجرة ، والقطب الشمالي المجري على التوالي.

        • المعلمات:
          • xyz (مجموعة من الفاصلة العائمة): ناقل ثلاثي العناصر يعطي مكونات مجرة ​​ديكارتية
          • (صفيف النقطة العائمة): متجه ثلاثي العناصر يعطي مكونات ICRS الديكارتية
          • galToIcrs (polarXYZ (l ، b ، المسافة))

          icrsToEcl (xyz) يحول متجهًا مكونًا من 3 عناصر يمثل إحداثيات ICRS (الاستوائية) إلى إحداثيات مسير الشمس. يمكن استخدام هذا مع متجهات الموضع أو السرعة.

          التحويل يتوافق مع ذلك بين الإحداثيات (ra ، dec) و (ecl_lon ، ecl_lat) في كتالوج مصدر Gaia (DR2).

          • المعلمات:
            • xyz (مجموعة من الفاصلة العائمة): متجه ثلاثي العناصر يعطي مكونات ICRS الديكارتية
            • (صفيف النقطة العائمة): ناقل ثلاثي العناصر يعطي مكونات ديكارتية مسير الشمس
            • icrsToEcl (بولاركسيز (را ، ديسمبر ، مسافة))

            eclToIcrs (xyz) يحول متجهًا مكونًا من 3 عناصر يمثل إحداثيات مسير الشمس إلى إحداثيات ICRS (الاستوائية). يمكن استخدام هذا مع متجهات الموضع أو السرعة.

            التحويل يتوافق مع ذلك بين الإحداثيات (ecl_lon ، ecl_lat) و (ra ، dec) في كتالوج مصدر Gaia (DR2).

            • المعلمات:
              • xyz (مجموعة من الفاصلة العائمة): متجه ثلاثي العناصر يعطي إحداثيات ديكارتية مسير الشمس
              • (صفيف النقطة العائمة): متجه ثلاثي العناصر يعطي إحداثيات ICRS الديكارتية
              • eclToIcrs (polarXYZ (ecl_lon ، ecl_lat ، المسافة))

              astromUVW (astrom6) لحساب مكونات السرعة الديكارتية من الكميات المتوفرة في كتالوج مصدر Gaia. الإخراج في نفس نظام الإحداثيات مثل المدخلات ، أي ICRS لكميات Gaia المسماة بالمقابل.

              يتم تمثيل معلمات قياس الإدخال بمصفوفة مكونة من 6 عناصر ، مع العناصر التالية: الوحدات المستخدمة بواسطة هذه الوظيفة هي الوحدات المستخدمة في جدول gaia_source.

              تستدعي وظيفة الراحة هذه فقط وظيفة astromUVW ذات 7 حجة باستخدام المنظر المقلوب للمسافة الشعاعية ، وبدون استدعاء تصحيح دوبلر. إنه مكافئ تمامًا لـ: لاحظ أن هذا الانعكاس الساذج للاختلاف المنظر لتقدير المسافة لا يمكن الاعتماد عليه بشكل عام في حالات اختلاف المنظر مع أخطاء غير قابلة للإهمال.

              • المعلمات:
                • نجم 6 (مجموعة من الفاصلة العائمة): متجه لـ 6 معلمات فلكية على النحو المنصوص عليه في كتالوج مصدر Gaia
                • (مجموعة من الفاصلة العائمة): متجه ثلاثي العناصر يعطي مكونات السرعة الاستوائية بالكيلومتر / ثانية
                • astromUVW (مجموعة (ra، dec، parallax، pmra، pmdec، radial_velocity))
                • icrsToGal (astromUVW (مصفوفة (ra، dec، parallax، pmra، pmdec، radial_velocity)))

                astromUVW (ra، dec، pmra، pmdec، radial_velocity، r_parsec، useDoppler) تحسب المكونات الديكارتية للسرعة من الموضع المرصود والحركة المناسبة والسرعة الشعاعية والمسافة الشعاعية ، مع تصحيح اختياري لوقت الضوء. الإخراج في نفس نظام الإحداثيات مثل المدخلات ، أي ICRS لكميات Gaia المسماة بالمقابل.

                يجب توفير المسافة الشعاعية باستخدام المعلمة r_parsec. قد يتم إجراء تقدير ساذج من الكميات الموجودة في كتالوج مصدر Gaia باستخدام التعبير 1000./ Parallax ، على الرغم من ملاحظة أن هذا الانعكاس البسيط للاختلاف المنظر لا يمكن الاعتماد عليه بشكل عام لأوجه الاختلاف مع أخطاء غير قابلة للإهمال.

                الحسابات واضحة إلى حد ما ، بعد المعادلة. 1.5.74 من كتالوج Hipparcos. يمكن أيضًا تطبيق عامل دوبلر (صغير عادةً) الذي يمثل تأثيرات وقت الضوء بشكل اختياري. تأثير هذا هو ضرب المتجه العائد بعامل 1 / (1-radial_velocity / c) ، كما تمت مناقشته في المعادلة. 1.2.21 من كتالوج Hipparcos.

                لاحظ أنه لم يتم إجراء أي محاولة لضبط الحركة الشمسية.

                • المعلمات:
                  • را (النقطة العائمة): حق الصعود بالدرجات
                  • ديسمبر (النقطة العائمة): الانحراف بالدرجات
                  • pmra (النقطة العائمة): الحركة المناسبة في RA * cos (ديسمبر) في mas / yr
                  • pmdec (النقطة العائمة): الحركة الصحيحة في الانحراف في ماس / سنة
                  • radial_velocity (النقطة العائمة): السرعة الشعاعية بالكيلومتر / ثانية
                  • r_parsec (النقطة العائمة): المسافة الشعاعية في الفرسخ
                  • استخدم دوبلر (قيمة منطقية): ما إذا كان يجب تطبيق عامل دوبلر لحساب تأثيرات وقت الضوء
                  • (صفيف النقطة العائمة): متجه ثلاثي العناصر يعطي مكونات السرعة الاستوائية بالكيلومتر / ثانية
                  • astromUVW (ra، dec، pmra، pmdec، radial_velocity، dist، true)
                  • icrsToGal (astromUVW (ra، dec، pmra، pmdec، radial_velocity، 1000./parallax، false))

                  epochProp (tYr، astrom6) ينشر معلمات القياس الفلكي ، التي يتم توفيرها كمصفوفة مكونة من 6 عناصر ، إلى حقبة مختلفة.

                  يتم تمثيل كل من معلمات قياس المدخلات والمخرجات بمصفوفة مكونة من 6 عناصر ، مع العناصر التالية: الوحدات المستخدمة بواسطة هذه الوظيفة هي الوحدات المستخدمة في جدول gaia_source.

                  • المعلمات:
                    • tYr (النقطة العائمة): فرق الزمن بالسنوات
                    • نجم 6 (صفيف النقطة العائمة): القياس الفلكي في الوقت t0 ، ممثلة بمصفوفة من 6 عناصر على النحو الوارد أعلاه (يُسمح أيضًا بمصفوفة من 5 عناصر حيث تكون السرعة الشعاعية صفرًا أو غير معروفة)
                    • (مجموعة من الفاصلة العائمة): قياس الفلك في الوقت t0 + tYr ، ممثلة بمصفوفة مكونة من 6 عناصر على النحو الوارد أعلاه
                    • epochProp (-15.5، array (ra، dec، parallax، pmra، pmdec، radial_velocity)) - يحسب القياس الفلكي عند 2000.0 من قيم gaia_source التي تمت ملاحظتها في 2015.

                    epochPropErr (tYr، astrom22) ينشر معلمات القياس الفلكي وما يرتبط بها من أخطاء وارتباطات ، يتم توفيرها كمصفوفة مكونة من 22 عنصرًا ، إلى حقبة مختلفة.

                    يتم تمثيل كل من معلمات قياس المدخلات والمخرجات مع الخطأ المرتبط ومعلومات الارتباط بمصفوفة مكونة من 22 عنصرًا ، مع العناصر التالية: لاحظ أن معاملات الارتباط ، دائمًا في النطاق -1..1 ، بلا أبعاد.

                    من الواضح أن هذه وظيفة غير عملية للاستدعاء ، ولكن إذا كنت تستخدمها مع كتالوج gaia_source نفسه ، أو مع كتالوجات أخرى مشابهة لها نفس أسماء الأعمدة والوحدات ، فيمكنك استدعاؤها بمجرد نسخ ولصق المثال الموضح في هذه الوثائق.

                    لا ينطبق هذا التحول إلا على السرعات الشعاعية المحددة بشكل مستقل عن القياس الفلكي ، مثل تلك التي تم الحصول عليها باستخدام مقياس الطيف. لا ينطبق على التحويل الخلفي للبيانات التي تم نشرها بالفعل إلى حقبة أخرى.

                    • المعلمات:
                      • tYr (النقطة العائمة): فرق الزمن بالسنوات
                      • نجم 22 (مجموعة من الفاصلة العائمة): قياس الفلك في الوقت t0 ، ممثلة بمصفوفة مكونة من 22 عنصرًا على النحو الوارد أعلاه
                      • (مجموعة من الفاصلة العائمة): قياس الفلك في الوقت t0 + tYr ، ممثلة بمصفوفة مكونة من 22 عنصرًا على النحو الوارد أعلاه
                      • epochPropErr (-15.5، مجموعة (رضي، ديسمبر، المنظر، الخاصة بالوكالة، pmdec، radial_velocity، ra_error، dec_error، parallax_error، pmra_error، pmdec_error، radial_velocity_error، ra_dec_corr، ra_parallax_corr، ra_pmra_corr، ra_pmdec_corr، dec_parallax_corr، dec_pmra_corr، dec_pmdec_corr، parallax_pmra_corr، parallax_pmdec_corr، pmra_pmdec_corr )) - بحساب القياس الفلكي مع جميع الأخطاء والارتباطات عند 2000.0 لقيم gaia_source التي تمت ملاحظتها في 2015.5.

                      rvMasyrToKms (rvMasyr ، plxMas) يتحول من السرعة الشعاعية الطبيعية في ماس / سنة إلى سرعة شعاعية غير طبيعية بالكيلومتر / ثانية.

                      يتم احتساب الناتج كـ AU_YRKMS * rvMasyr / plxMas ، حيث AU_YRKMS = 4.740470446 هي وحدة فلكية واحدة بالكيلومتر.سنة / ثانية.

                      • المعلمات:
                        • rv مصر (النقطة العائمة): السرعة الشعاعية الطبيعية ، بالماس / السنة
                        • plxMas (النقطة العائمة): اختلاف المنظر في ماس
                        • (النقطة العائمة): السرعة الشعاعية بالكيلومتر / ثانية

                        rvKmsToMasyr (rvKms ، plxMas) يتحول من سرعة شعاعية غير طبيعية بالكيلومتر / ثانية إلى سرعة شعاعية طبيعية في ماس / سنة.

                        يتم احتساب الناتج كـ rvKms * plxMas / AU_YRKMS ، حيث AU_YRKMS = 4.740470446 هي وحدة فلكية واحدة بالكيلومتر.سنة / ثانية.

                        • المعلمات:
                          • rvKms (النقطة العائمة): سرعة شعاعية غير طبيعية ، بالماس / السنة
                          • plxMas (النقطة العائمة): اختلاف المنظر في ماس
                          • (النقطة العائمة): السرعة الشعاعية في ماس / السنة

                          DistanceEstimateEdsd (plxMas ، plxErrorMas ، lPc) أفضل تقدير للمسافة باستخدام كثافة الفضاء المتناقصة بشكل كبير مسبقًا. يتم توفير هذا التقدير من خلال صيغة PDF.

                          • المعلمات:
                            • plxMas (النقطة العائمة): اختلاف المنظر في ماس
                            • plxErrorMas (النقطة العائمة): خطأ اختلاف المنظر في ماس
                            • lPc (النقطة العائمة): مقياس الطول بالفرسخ
                            • (النقطة العائمة): أفضل تقدير للمسافة في الفرسخ

                            DistanceBoundsEdsd (plxMas ، plxErrorMas ، lPc) لحساب فترتي الثقة المئوية الخامس والتسعين على تقدير المسافة باستخدام القيمة السابقة لكثافة المساحة المتناقصة بشكل أسي.

                            لاحظ أن هذه الوظيفة يجب أن تدمج PDF عدديًا لتحديد القيم الكمية ، لذا فهي بطيئة نسبيًا.

                            • المعلمات:
                              • plxMas (النقطة العائمة): اختلاف المنظر في ماس
                              • plxErrorMas (النقطة العائمة): خطأ اختلاف المنظر في ماس
                              • lPc (النقطة العائمة): مقياس الطول بالفرسخ
                              • (مجموعة من الفاصلة العائمة): مصفوفة مكونة من عنصرين تعطي المئين الخامس والتسعين 95 بالفرسخ من مسافة EDSD PDF

                              DistanceQuantilesEdsd (plxMas ، plxErrorMas ، lPc ، qpoints ،.) لحساب الكميات التعسفية لتقدير المسافة باستخدام التناقص الأسي لكثافة الفضاء السابقة.

                              لاحظ أن هذه الوظيفة يجب أن تدمج PDF عدديًا لتحديد القيم الكمية ، لذا فهي بطيئة نسبيًا.

                              • المعلمات:
                                • plxMas (النقطة العائمة): اختلاف المنظر في ماس
                                • plxErrorMas (النقطة العائمة): خطأ اختلاف المنظر في ماس
                                • lPc (النقطة العائمة): مقياس الطول بالفرسخ
                                • كيوبوينتس (نقطة عائمة ، واحد أو أكثر): نقطة أو أكثر من نقاط القطع الكمي المطلوبة ، كل منها في النطاق 0..1
                                • (صفيف النقطة العائمة): صفيف يحتوي على عنصر واحد لكل من الكيوبوينتس المزودة مع إعطاء المسافة المقابلة بالفرسخ
                                • DistanceQuantilesEdsd (parallax، parallax_error، 1350، 0.5) [0] تحسب متوسط ​​مسافة EDSD PDF باستخدام مقياس طول 1.35 كيلو قطعة
                                • تُرجع الدالة DistanceQuantilesEdsd (اختلاف المنظر ، خطأ المنظر ، 3000 ، 0.01 ، 0.99) مصفوفة مكونة من عنصرين تعطي المئين الأول و 99 من تقدير المسافة باستخدام مقياس طول 3 كيلو قطعة

                                DistanceToModulus (distPc) تحويل المسافة بالفرسخ إلى معامل مسافة. الصيغة هي 5 * log10 (distPc) -5.

                                • المعلمات:
                                  • توزيعات (النقطة العائمة): المسافة في الفرسخ
                                  • (النقطة العائمة): معامل المسافة بالأحجام

                                  المعامل إلى المسافة (distmod) تحويل معامل المسافة إلى مسافة بالفرسخ. الصيغة هي 10 ^ (1 + distmod / 5).

                                  • المعلمات:
                                    • distmod (النقطة العائمة): معامل المسافة بالأحجام
                                    • (النقطة العائمة): المسافة في الفرسخ

                                    AU_YRKMS هذه الكمية هي A_v ، الوحدة الفلكية معبرًا عنها بالكيلومتر.سنة / ثانية. راجع كتالوج Hipparcos (ESA SP-1200) جدول 1.2.2 و Eq. 1.5.24.

                                    PC_AU فرسخ فلكي في الوحدات الفلكية ، يساوي 648000 / PI.


                                    معلومات الكاتب

                                    الانتماءات

                                    قسم الفيزياء والفلك ، جامعة شيفيلد ، شيفيلد ، المملكة المتحدة

                                    ستيفن جي بارسونز ، ألكسندر ج.براون ، ستيوارت بي ليتلفير ، فيكرام إس ديلون ، مارتن جيه داير وأمبير ديفيد آي ساهمان

                                    Instituto de Astrofísica de Canarias ، لا لاغونا ، تينيريفي ، إسبانيا

                                    قسم الفيزياء ، جامعة وارويك ، كوفنتري ، المملكة المتحدة

                                    توماس آر مارش وأمبير ماثيو ج. جرين

                                    قسم علم الفلك ، جامعة بوسطن ، بوسطن ، ماساتشوستس ، الولايات المتحدة الأمريكية

                                    قسم الفيزياء الفلكية ، جامعة رادبود نيميغن ، نيميغن ، هولندا

                                    معهد الفلك ، جامعة كامبريدج ، كامبريدج ، المملكة المتحدة

                                    يمكنك أيضًا البحث عن هذا المؤلف في PubMed Google Scholar

                                    يمكنك أيضًا البحث عن هذا المؤلف في PubMed Google Scholar

                                    يمكنك أيضًا البحث عن هذا المؤلف في PubMed Google Scholar

                                    يمكنك أيضًا البحث عن هذا المؤلف في PubMed Google Scholar

                                    يمكنك أيضًا البحث عن هذا المؤلف في PubMed Google Scholar

                                    يمكنك أيضًا البحث عن هذا المؤلف في PubMed Google Scholar

                                    يمكنك أيضًا البحث عن هذا المؤلف في PubMed Google Scholar

                                    يمكنك أيضًا البحث عن هذا المؤلف في PubMed Google Scholar

                                    يمكنك أيضًا البحث عن هذا المؤلف في PubMed Google Scholar

                                    يمكنك أيضًا البحث عن هذا المؤلف في PubMed Google Scholar

                                    يمكنك أيضًا البحث عن هذا المؤلف في PubMed Google Scholar

                                    مساهمات

                                    ساهم جميع المؤلفين في العمل المقدم في هذه الورقة. S.G.P. قلل من جميع البيانات الطيفية والقياسية الضوئية ونفذ تركيب SED. S.G.P. و V.S.D. أجرى ملاحظات Gran Telescopio Canarias. أ. نفذت السرعة الشعاعية وتركيب منحنى الضوء. م. كتب كود Python الذي نفذ العمليات Gaussian في تركيب المنحنى الخفيف. ، S.G.P. ، E.B. ، M.J.D. ، M.J.G. و د. كل ذلك ساهم في تطوير ودعم HiPERCAM. ج. ج. حلل النبضات في المنحنى الضوئي HiPERCAM. أ. بحث في التركيب الداخلي للأقزام البيضاء وتطور الثنائي. راجع جميع المؤلفين المخطوطة.

                                    المؤلف المراسل


                                    • ESA / ESAC: http://gea.esac.esa.int/archive/
                                    • CDS: http://cdsweb.u-strasbg.fr/gaia
                                    • ARI: http://gaia.ari.uni-heidelberg.de/index.html
                                    • AIP: https://gaia.aip.de/
                                    • ASI / ASCD: http://gaia.asdc.asi.it/
                                    • كتاب طبخ ADQL Gaia DR1: https://gaia.ac.uk/data/gaia-data-release-1/adql-cookbook
                                    • ورقة المعلومات السريعة لـ CDS ADQL: http://tapvizier.u-strasbg.fr/adql/help.html
                                    • ورقة تعليمات GAVO ADQL: http://docs.g-vo.org/adqlref/adqlref.pdf
                                    • صفحة البرنامج التعليمي لـ ARI ADQL: http://docs.g-vo.org/adql-gaia مع الشرائح المرتبطة http://docs.g-vo.org/adql-gaia/html/pop.pdf

                                    تثبيت العملاء "المساعدون":

                                    Topcat (الإصدار الأخير ، ويفضل v4.6 +): انظر http://www.star.bris.ac.uk/

                                    mbt / topcat / # تثبيت
                                    Aladin (الإصدار الأخير ، ويفضل v9 +): انظر http://aladin.u-strasbg.fr/java/nph-aladin.pl؟frame=downloading
                                    Java مطلوب من قبل Topcat و Aladin: أي إصدار حديث من Java لا بأس به ، لكن Oracle / Hotspot Java SE (من https://www.oracle.com/technetwork/java/index.html) موصى به بدلاً من OpenJDK. في نظام التشغيل Linux ، يمكنك تشغيل "إصدار جافا" لمعرفة ما لديك.

                                    تعليمات لتثبيت بايثون

                                    إذا لم يكن لديك بالفعل مكدس Python (علمي) على جهاز الكمبيوتر الخاص بك ، فيمكنك بسهولة تثبيت مكدس Anaconda. يأتي هذا مع "بطاريات مضمنة" ، وتأتي مليئة بالنومبي (الثعبان العددي) ، و scipy (الثعبان العلمي) ، و Astropy (مكتبة Python الفلكية) ، و matplotlib (مكتبة التخطيط).

                                    يمكنك الحصول على إصدار Python 2.7 أو Python 3.5 Anaconda من:

                                    إذا كنت لا تعرف إصدار Python الذي تختاره ، فحدد Python 3.5 (مثبت رسومي أو سطر أوامر - إما أنه موافق).


                                    رسم بيانات السرعة الشعاعية GAIA barycentric باستخدام بيثون - علم الفلك

                                    جو بوفي - بوفي في أسترو دوت أوتورونتو دوت كاليفورنيا

                                    راجع AUTHORS.txt للحصول على قائمة كاملة بالمساهمين.

                                    إذا وجدت هذا الرمز مفيدًا في بحثك ، فيرجى إبلاغي بذلك. إذا كنت تستخدم galpy في منشور ، فيرجى الاستشهاد به بوفي (2015) ورابط إلى http://github.com/jobovy/galpy. راجع قسم توثيق الإقرار للحصول على دليل أكثر تفصيلاً للاستشهاد بأجزاء من الكود. شكرا!

                                    يمكن العثور على أحدث الوثائق هنا. يمكنك أيضًا الانضمام إلى مجتمع galpy slack لأية أسئلة متعلقة بـ galpy.

                                    إذا وجدت أي خطأ في الكود ، يرجى الإبلاغ عن ذلك باستخدام متعقب المشكلات أو من خلال الانضمام إلى مجتمع galpy slack.

                                    إذا كنت تواجه مشكلات في تثبيت galpy ، فيرجى الرجوع أولاً إلى الأسئلة الشائعة حول التثبيت.

                                    إصدارات PYTHON والاعتمادات

                                    يدعم galpy لغة Python 3. على وجه التحديد ، يدعم galpy بايثون 3.7 و 3.8 و 3.9. يجب أن يعمل أيضًا على إصدارات Python 3. * السابقة ، ولكن هذا لم يتم اختباره على نطاق واسع على أساس مستمر ولأن المكتبات التي يعتمد عليها galpy تتخلى عن إصدارات Python 3. * السابقة ، فمن المحتمل ألا يعمل galpy نفسه بشكل كامل عليها. يبني GitHub Actions CI للتحقق بانتظام من دعم Python 3.9 (و 3.7 و 3.8 باستخدام مجموعة أكثر محدودية من الاختبارات الأساسية) على أنظمة Linux و Windows Appveyor تحقق بانتظام من دعم Python 3.9 على Windows. في حين أن معظم galpy يعمل على الأرجح في Python 2.7 ، إلا أنه لم يعد يتم اختباره ولم يعد Python 2.7 مدعومًا رسميًا.

                                    تتطلب هذه الحزمة Numpy و Scipy و Matplotlib. تتطلب بعض الميزات المتقدمة مكتبة GNU العلمية (GSL) ، مع حسابات الإجراءات التي تتطلب الإصدار 1.14 أو أعلى. تشمل التبعيات الاختيارية الأخرى:

                                    • يتم توفير الدعم لتوفير المدخلات والحصول على المخرجات مثل الكميات مع الوحدات من خلال Astropy.
                                    • يتطلب الاستعلام عن SIMBAD عن إحداثيات كائن في طريقة التهيئة Orbit.from_name استعلامًا استراتيجيًا.
                                    • يتطلب رسم الوظائف التعسفية لسمات Orbit numexpr.
                                    • تتطلب مدافع التباين الثابت في galpy.df.constantbetadf JAX.
                                    • يتطلب استخدام SnapshotRZPotential و InterpSnapshotRZPotential استخدام شخص واحد.

                                    قد تتطلب أجزاء أخرى من الكود حزمًا إضافية وسيتم تنبيهك بواسطة الكود إذا لم يتم تثبيتها.

                                    إذا كنت مهتمًا بالمساهمة في تطوير galpy ، فقم بإلقاء نظرة على هذا الدليل الموجز على الويكي. نأمل أن يساعدك هذا على البدء!

                                    يمكن العثور على بعض ملاحظات التطوير الإضافية على الويكي. يتضمن هذا قائمة بالامتدادات الصغيرة والأكبر من galpy والتي ستكون مفيدة هنا بالإضافة إلى خارطة طريق طويلة المدى هنا. يرجى إخبار المطور الرئيسي إذا كنت بحاجة إلى أي مساعدة للمساهمة!


                                    الشكل 8: التوزيعات الخلفية ثنائية الأبعاد لجميع معلمات النموذج. المعلمات الخمس الأولى للإشارة Keplerian ، b متخالفة ، و ln τ هما على التوالي مقياس الاتساع والوقت اللوغاريتمي لنموذج MA. يتم الإشارة إلى أسماء مجموعات البيانات بالأحرف العلوية للمعلمات ذات الصلة. تشير الأحرف المرتفعة "pre" و "post" إلى مجموعتي HARPSpre و HARPSpost ، على التوالي. تُظهر المخططات مستويات الثقة 68٪ و 90٪. تشير الخطوط المتقطعة إلى قيم معلمات MAP في التوزيع اللاحق بينما تُظهر الخطوط الصلبة قيم التوزيع بنسبة 10٪ و 90٪.

                                    الشكل 9: التوزيع الخلفي للكتلة والفترة المدارية لـ ϵ Indi A b. يتم حساب الكفاف من عينة فرعية مرسومة بواسطة سلاسل MCMC من التوزيع اللاحق. ثلاثة مستويات كفاف مقابلة لمستويات واحد ، اثنين وثلاثة سيجما موضحة بخطوط سوداء. يتم عرض قيمة MAP (خطوط صلبة بيضاء) والكميات 10٪ و 90٪ (خطوط بيضاء متقطعة) للفترة المدارية وكتلة الكواكب في الرسوم البيانية للكتلة والفترة المدارية. يشير الصليب الأسود في المحيط إلى حل MAP. يختلف حل MAP هذا عن الحل المتوسط ​​لأن توزيعات الكتلة والمعلمات المدارية ليست غاوسية وغير متماثلة.

                                    السمات مقابل الوظائف

                                    الشكل هو سمة ، لذلك نعرض قيمته دون تسميته كدالة.

                                    الرأس دالة ، لذلك نحتاج إلى الأقواس.

                                    الآن يمكننا استخراج الأعمدة التي نريدها من skycoord_gd1 وإضافتها كأعمدة في DataFrame. phi1 و phi2 يحتويان على الإحداثيات المحولة.

                                    تحتوي pm_phi1_cosphi2 و pm_phi2 على مكونات الحركة المناسبة في الرتل المحول.

                                    التفاصيل إذا لاحظت أن لدى SkyCoord سمة تسمى الحركة المناسبة ، فقد تتساءل عن سبب عدم استخدامها.

                                    يمكن أن يكون لدينا: Prop_motion يحتوي على نفس البيانات مثل pm_phi1_cosphi2 و pm_phi2 ، ولكن بتنسيق مختلف.

                                    قبل أن نذهب إلى أبعد من ذلك ، سنتخذ جميع الخطوات التي قمنا بها ودمجها في وظيفة واحدة يمكننا استخدامها لأخذ الإحداثيات والحركة المناسبة التي نحصل عليها كجدول Astropy من استعلام Gaia ، وقم بتحويلها إلى a Pandas DataFrame ، وإضافة أعمدة تمثل إحداثيات GD-1 المصححة منعكسة والحركات المناسبة. هذه وظيفة عامة سنستخدمها عدة مرات أثناء قيامنا ببناء استعلامات مختلفة لذلك نريد كتابتها مرة واحدة ثم استدعاء الوظيفة بدلاً من الاضطرار إلى نسخ الكود ولصقه مرارًا وتكرارًا.

                                    إليك كيف نستخدم الوظيفة:


                                    العثور على رفقاء ممتازين لم يتم حلهم في مضيفات كوكب كبلر باستخدام Gaia

                                    في اجتماع هذا الأسبوع ، قمت بتحديث المجموعة حول مشروعي مع D F-M ، وآندي كيسي (جامعة موناش) ، وتيم مورتون (USC- الساحل الغربي الأول). نحن نبحث في كيفية تحديد كيفية تأثير ثنائية مضيفي كوكب كبلر على تكوين الكواكب. هذا يبدو مثيرًا حقًا ولكن كلمة "تحديد" تقوم بالكثير من العمل في تلك الجملة. أولاً ، علينا معرفة النجوم التي (لم يتم حلها) ثنائيات. هناك طريقتان مختلفتان سنقوم بهما للقيام بذلك ، ولكن الطريقة التي تم تسليط الضوء عليها هذا الأسبوع هي ثمرة عمل آندي.

                                    إلى جانب القياسات الفلكية وقياسات RV ، يسرد Gaia أيضًا الأخطاء المرتبطة بالنوبات. بالنسبة للثنائيات التي لم يتم حلها ، فإن محاولة ملاءمة المنحنى الفلكي أو منحنى RV لنجم واحد سيؤدي إلى قدر كبير من الخطأ بشكل غير معهود. من خلال عمل آندي (الذي لا يمكنني الإدعاء بفهمه لما هو أبعد من مستوى السطح) ، يمكننا تحويل هذا الخطأ إلى احتمال ألا يكون النجم نجمًا واحدًا.

                                    بشكل لا يصدق ، يبدو أن هذا يعمل بشكل جيد حقًا! في الرسم التخطيطي لوني الحجم أعلاه لنجوم KIC بقياسات فلكية ، يمكننا اختيار النجوم التي لديها احتمالية منخفضة لكونها مفردة فوق التسلسل الرئيسي حيث نتوقع أن تكون الثنائيات.

                                    ومن المثير للاهتمام أن نفس الحبكة الخاصة بنجوم KIC ذات قياسات RV تبدو مختلفة تمامًا. هناك عدد قليل من تأثيرات الاختيار هنا سيتعين علينا اكتشافها: ليست كل النجوم لها قياسات فلكية وقياسات RV (كما أوضح تريفور) ، يختلف عدد الثنائيات التي اكتشفها القياس الفلكي عن ثنائيات RV. يتم التقاط الثنائيات الطيفية والكسوف المعروفة من خلال اهتزاز RV ولكن ليس الارتعاش الفلكي!

                                    سيكون اكتشاف مدى دقة هذا أمرًا أساسيًا. نظرًا لأن أيا من الطريقتين لا تضرب 1.000 ، فإن كيفية تحديد "أنا متأكد تمامًا من أن هذا النجم ليس وحيدًا" يمثل لغزًا مثيرًا للاهتمام لأنفسنا في المستقبل القريب.


                                    شاهد الفيديو: Man Catches The Worlds Most Dangerous Snake (قد 2022).