Models of long-period variables of the globular cluster 47 Tuc
Fadeyev
Stellar evolution computations were carried out for stars with a main sequence mass $M_\mathrm{ZAMS}=0.86M_\odot$ and initial metal abundance $Z=0.003$ and $Z=0.004$. Selected models of evolutionary sequences were used for calculation of radial pulsations in the RGB, eAGB and TP-AGB evolutionary stages. Not all pulsating red giants of the globular cluster 47 Tuc are shown to belong to the Mira variables because the lower limit of pulsation periods at the TP-AGB stage is $\approx 70$ day, whereas during the eAGB evolutionary stage the periods of radial oscillations range from $\approx 5$ to $\approx 40$ day. Periods and luminosities of hydrodynamic models of eAGB and TP--AGB pulsating stars locate along the common period-luminosity relation. Small masses of Mira variables in the globular cluster 47 Tuc ($0.54M_\odot\le M\le 0.70M_\odot$) is the main reason for irregular large-amplitude oscillations and the dynamical instability of outer stellar layers at pulsation periods $Î > 200$ day.
এই গবেষণায় প্রধান ক্রম ভর MZAMS=0.86M⊙, প্রাথমিক ধাতব প্রাচুর্য Z=0.003 এবং Z=0.004 সহ তারকাদের বিবর্তনীয় গণনা পরিচালিত হয়েছে। বিবর্তনীয় ক্রম থেকে নির্বাচিত মডেলগুলি RGB (লাল দৈত্য শাখা), eAGB (প্রাথমিক অ্যাসিম্পটোটিক দৈত্য শাখা) এবং TP-AGB (তাপীয় নাড়ি অ্যাসিম্পটোটিক দৈত্য শাখা) পর্যায়ে রেডিয়াল দোলন গণনার জন্য ব্যবহৃত হয়েছে। গবেষণা দেখায় যে গ্লোবুলার ক্লাস্টার ৪৭ টাকের দোলনশীল লাল দৈত্যগুলি সবই মিরা চলক নয়: TP-AGB পর্যায়ের দোলনশীল পর্যায়ের নিম্ন সীমা প্রায় ৭০ দিন, যখন eAGB পর্যায়ের রেডিয়াল দোলনের পর্যায় পরিসীমা প্রায় ৫ থেকে ৪০ দিন। eAGB এবং TP-AGB দোলনশীল তারকাদের তরল গতিশীল মডেলের পর্যায় এবং উজ্জ্বলতা একটি সাধারণ পর্যায়-উজ্জ্বলতা সম্পর্ক বরাবর বিতরণ করা হয়। ৪৭ টাকে মিরা চলকের ছোট ভর (0.54M⊙≤M≤0.70M⊙) পর্যায় Π>200 দিনে অনিয়মিত বৃহৎ-বিস্তৃতি দোলন এবং বাহ্যিক গতিশীল অস্থিতিশীলতার উপস্থিতির প্রধান কারণ।
এই গবেষণার লক্ষ্য গ্লোবুলার ক্লাস্টার ৪৭ টাকের লাল দৈত্যগুলির বিভিন্ন বিবর্তনীয় পর্যায়ে (RGB, eAGB, TP-AGB) রেডিয়াল দোলন উত্তেজনের শর্তগুলি নির্ধারণ করা এবং এই দীর্ঘ-পর্যায়ের চলক তারকাদের বিবর্তনীয় অবস্থা স্পষ্ট করা।
১. পর্যবেক্ষণগত বিশেষত্ব: ৪৭ টাক (NGC 104) অস্বাভাবিকভাবে উচ্চ ধাতব প্রাচুর্যের একটি গ্লোবুলার ক্লাস্টার, যা বেশিরভাগ গ্লোবুলার ক্লাস্টারের চেয়ে প্রায় এক মাত্রার বেশি, যা এটিকে বিপুল সংখ্যক দোলনশীল লাল দৈত্য রাখে (২০০১ সালের মধ্যে ১০২টি গ্লোবুলার ক্লাস্টারে আবিষ্কৃত ১১৭টি দীর্ঘ-পর্যায়ের চলকের মধ্যে ৪৭ টাক ১৪টি অন্তর্ভুক্ত করে)
२. ধুলো সংঘনন ঘটনা: অবলোহিত পর্যবেক্ষণ ৪৭ টাকের দীর্ঘ-পর্যায়ের চলকের বাহ্যিক বায়ুমণ্ডলে ধুলো কণার সংঘনন প্রদর্শন করে, যা বৃহৎ-বিস্তৃতি তারকা দোলন দ্বারা উত্পাদিত পর্যায়ক্রমিক শকওয়েভের সাথে ঘনিষ্ঠভাবে সম্পর্কিত
३. বিবর্তনীয় অবস্থা অস্পষ্ট: বিদ্যমান গবেষণা এই দোলনশীল লাল দৈত্যগুলির বিবর্তনীয় অবস্থার বিষয়ে বিরোধ রয়েছে, বিশেষ করে তারা সবই মিরা চলক কিনা তা সম্পর্কে
१. রৈখিক দোলন বিশ্লেষণ অপর্যাপ্ত: লেবজেল্টার এবং উড (২০০৫) Z=0.004, M=0.9M⊙ সহ তারকাদের অধ্যয়নের জন্য রৈখিক দোলন মডেল ব্যবহার করেছেন, কিন্তু প্রমাণ প্রদান করতে পারেননি
२. ধাতব প্রাচুর্য অনুমান পক্ষপাত: প্রাথমিক গবেষণা আধুনিক পর্যবেক্ষণ অনুমানের চেয়ে উল্লেখযোগ্যভাবে বেশি ধাতব প্রাচুর্য ব্যবহার করেছে
३. সামঞ্জস্যপূর্ণ গণনার অভাব: তারকা বিবর্তন এবং অ-রৈখিক দোলন একত্রিত করার সামঞ্জস্যপূর্ণ গণনার অভাব
তারকা বিবর্তন এবং অ-রৈখিক দোলন সামঞ্জস্যপূর্ণভাবে গণনা করে, বিবর্তনীয় ক্রম থেকে নির্বাচিত মডেলগুলি তরল গতিশীল সমীকরণের প্রাথমিক শর্ত হিসাবে ব্যবহার করে, বিভিন্ন বিবর্তনীয় পর্যায়ে ৪৭ টাক লাল দৈত্যের দোলন বৈশিষ্ট্য সিস্টেমেটিকভাবে অধ্যয়ন করা।
१. সম্পূর্ণ বিবর্তন-দোলন গণনা কাঠামো প্রতিষ্ঠা: MESA কোড (সংস্করণ r15140) ব্যবহার করে তারকা বিবর্তন গণনা, বিকিরণ তরল গতিশীলতা এবং সময়-পরিবর্তনশীল পরিচলন সমীকরণ সহ তারকা দোলন গণনা
२. বিভিন্ন বিবর্তনীয় পর্যায়ের দোলন বৈশিষ্ট্য স্পষ্ট করা:
RGB পর্যায়: দোলন অস্থিতিশীল অঞ্চল সংকীর্ণ বা অস্তিত্বহীন (Z=0.003 এ)
eAGB পর্যায়: পর্যায় পরিসীমা ৫-৪০ দিন, তুলনামূলকভাবে নিয়মিত বিস্তৃতি
TP-AGB পর্যায়: পর্যায় নিম্ন সীমা প্রায় ৭০ দিন, দীর্ঘ পর্যায় (>২০০ দিন) এ অনিয়মিত বৃহৎ-বিস্তৃতি দোলন
३. ৪৭ টাক দীর্ঘ-পর্যায়ের চলকের প্রধান বিবর্তনীয় অবস্থা নির্ধারণ: বেশিরভাগ দোলনশীল লাল দৈত্য TP-AGB পর্যায়ের পরিবর্তে eAGB পর্যায়ে রয়েছে
४. পর্যায়-উজ্জ্বলতা সম্পর্ক প্রতিষ্ঠা: eAGB এবং TP-AGB তারকাদের পর্যায়-উজ্জ্বলতা ডেটা একটি সাধারণ রিগ্রেশন লাইন বরাবর বিতরণ করা হয়
५. গতিশীল অস্থিতিশীলতা প্রক্রিয়া ব্যাখ্যা: ছোট ভর (0.54−0.70M⊙) দীর্ঘ পর্যায়ে বাহ্যিক অস্থিতিশীলতার প্রধান কারণ
তারকা বিবর্তন গণনা এবং অ-রৈখিক দোলন গণনা নির্বিঘ্নে সংযুক্ত করা, বিবর্তনীয় মডেল সরাসরি তরল গতিশীল প্রাথমিক শর্ত হিসাবে কাজ করে, পরামিতি অমিল সমস্যা এড়ায়।
নির্দিষ্ট যান্ত্রিক কাজের রেডিয়াল বিতরণ ∮PdV এর মাধ্যমে, দোলন উত্তেজনা/স্যুতি তে হাইড্রোজেন এবং হিলিয়াম আয়নীকরণ অঞ্চলের অবদান পরিমাণগতভাবে বিশ্লেষণ করা:
আংশিক আয়নীকৃত হাইড্রোজেন অঞ্চল (অ্যাডিয়াবেটিক সূচক Γ1<4/3) প্রধান চালিকা অঞ্চল
१. দোলন পর্যায়Π (দিন)
२. দোলন ধ্রুবকQ=Πρˉ/ρ⊙ (দিন, ρˉ গড় ঘনত্ব)
३. আপেক্ষিক রেডিয়াল স্থানচ্যুতি বিস্তৃতিΔR/R
४. দোলন অস্থিতিশীলতা স্থায়িত্বΔtev (মিলিয়ন বছর)
५. উজ্জ্বলতাlog(L/L⊙)
মূল সিদ্ধান্ত:
१. Z=0.003,ηR=0.3 মডেল দোলিত হয় না (দোলন অস্থিতিশীল অঞ্চল অস্তিত্বহীন)
२. ধাতব প্রাচুর্য এবং ভর ক্ষয় হার বৃদ্ধি বিবর্তনীয় ট্র্যাজেক্টরি আরও কম কার্যকর তাপমাত্রার দিকে স্থানান্তরিত করে
३. উচ্চতর Z এবং ছোট ভর তারকা রেডিয়াল দোলনের প্রতি আরও অস্থিতিশীল করে তোলে (নিম্ন ঘনত্ব বাহ্যিক পরিচলন অঞ্চল, বিস্তৃত আয়নীকরণ অঞ্চল)
४. RGB শীর্ষের কাছাকাছি দোলন স্থায়িত্ব খুবই সংক্ষিপ্ত (<१ Myr)
१. দোলন উপস্থিতির ক্রম: প্রথমে প্রথম ওভারটোন হিসাবে উপস্থিত হয়, তারকা ব্যাসার্ধ এবং উজ্জ্বলতা বৃদ্ধির সাথে, মৌলিক মোডে স্যুইচ করে
२. পর্যায় পরিসীমা: ५-४० দিন
३. বিস্তৃতি বৈশিষ্ট্য: ΔR/R≲0.3, দোলন তুলনামূলকভাবে নিয়মিত
४. উজ্জ্বলতা উপরের সীমা: যখন L≳103L⊙ দোলন বন্ধ হয়
५. স্থায়িত্ব: প্রায় 1.2×107 বছর (TP-AGB পর্যায়ের চেয়ে অনেক দীর্ঘ)
१. পর্যায় পরিসীমা: ७०-२४० দিন
२. অনিয়মিত দোলন থ্রেশহোল্ড: Π≈170 দিনে ΔR/R≈1 এর অনিয়মিত বৃহৎ-বিস্তৃতি দোলনে রূপান্তর
३. পর্যায় নিম্ন সীমা: প্রায় ७० দিন (আকাশগঙ্গা মিরা চলকের সাথে সামঞ্জস্যপূর্ণ)
१. বিবর্তনীয় অবস্থা নির্ধারণ: ४७ টাকের বেশিরভাগ দোলনশীল লাল দৈত্য TP-AGB পর্যায়ের পরিবর্তে eAGB পর্যায়ে রয়েছে
२. পর্যায় সীমানা:
পর্যায় < ७० দিন → eAGB তারকা
পর্যায় ≥ ७० দিন → সম্ভবত TP-AGB তারকা
३. ভর প্রভাব: ছোট ভর (0.54−0.70M⊙) দীর্ঘ পর্যায়ে বাহ্যিক গতিশীল অস্থিতিশীলতার প্রধান কারণ (>२०० দিন)
४. ধুলো সংঘনন প্রক্রিয়া: বৃহৎ-বিস্তৃতি দোলন দ্বারা উত্পাদিত পর্যায়ক্রমিক শকওয়েভ ধুলো কণা সংঘনের প্রধান প্রক্রিয়া
१. Clement et al. (२००१): १०२টি গ্লোবুলার ক্লাস্টারে ११७টি দীর্ঘ-পর্যায়ের চলক আবিষ্কার, ४७ টাক १४টি অন্তর্ভুক্ত
२. Lebzelter & Wood (२००५), Lebzelter et al. (२००५): আরও অনেক ४७ টাক দীর্ঘ-পর্যায়ের চলক আবিষ্কার, মোট সংখ্যা দ্বিগুণ
३. Percy & Gupta (२०२१): নতুন চলক আবিষ্কারের জন্য ক্রমাগত পর্যবেক্ষণ
४. Origlia et al. (२००२), McDonald et al. (२०११): অবলোহিত পর্যবেক্ষণ ধুলো সংঘনের প্রমাণ আবিষ্কার
१. Lebzelter & Wood (२००५): রৈখিক দোলন বিশ্লেষণ, Z=0.004,M=0.9M⊙ ব্যবহার, কিন্তু সিদ্ধান্ত অপর্যাপ্ত
२. Willson (२०००): পর্যায়ক্রমিক শকওয়েভ ধুলো সংঘনের দিকে পরিচালিত করার শারীরিক প্রক্রিয়া প্রস্তাব
३. Fadeyev (२०१३, २०२३, २०२४): লেখকের পূর্ববর্তী কাজ বিকিরণ তরল গতিশীলতা এবং সময়-পরিবর্তনশীল পরিচলন সমীকরণ কাঠামো প্রতিষ্ঠা
१. আরও নির্ভুল পরামিতি: আধুনিক ধাতব প্রাচুর্য অনুমান ব্যবহার (Z=0.003[Fe/H]=−0.67 এর সাথে সামঞ্জস্যপূর্ণ)
२. অ-রৈখিক গণনা: রৈখিক অনুমানের পরিবর্তে সম্পূর্ণ তরল গতিশীল সমীকরণ গ্রহণ
३. সামঞ্জস্যপূর্ণ কাঠামো: বিবর্তন এবং দোলন গণনা নির্বিঘ্নে সংযুক্ত
४. সিস্টেমেটিক বিশ্লেষণ: RGB, eAGB, TP-AGB তিনটি বিবর্তনীয় পর্যায় কভার করে
१. ভর ক্ষয় অনিশ্চয়তা: যদিও দুটি ηR মান বিবেচনা করা হয়েছে, প্রকৃত ভর ক্ষয় প্রক্রিয়া আরও জটিল হতে পারে
२. ধাতব প্রাচুর্য পরিসীমা: Z=0.004 ४७ টাক পর্যবেক্ষণ পরিসীমার উপরের সীমা অতিক্রম করে, প্রধানত পূর্ববর্তী কাজের সাথে তুলনার জন্য ব্যবহৃত
३. পরিচলন চিকিত্সা: সময়-পরিবর্তনশীল পরিচলন মডেল তাত্ত্বিক অনিশ্চয়তা রয়েছে
४. তাপীয় নাড়ির সংখ্যা: শুধুমাত্র ४টি তাপীয় নাড়ি অনুকরণ করা হয়েছে, TP-AGB পর্যায় সম্পূর্ণভাবে বর্ণনা করার জন্য অপর্যাপ্ত হতে পারে
१. আরও তাপীয় নাড়ি পর্যায়ে প্রসারিত করা, দেরী TP-AGB বিবর্তন অধ্যয়ন করা
२ ভর ক্ষয় হার সূত্র উন্নত করা, দোলন-বর্ধিত ভর ক্ষয় বিবেচনা করা
३ ত্রিমাত্রিক তরল গতিশীলতা অনুকরণ, পরিচলন এবং শকওয়েভ প্রচার আরও সঠিকভাবে বর্ণনা করা
४ পর্যবেক্ষণ ডেটার সাথে পরিসংখ্যানগত তুলনা, মডেল পূর্বাভাস যাচাই করা
१. উচ্চ ধাতব প্রাচুর্য গ্লোবুলার ক্লাস্টার (যেমন ω Cen) এর দীর্ঘ-পর্যায়ের চলক গবেষণা
२. নিম্ন ভর (∼0.86M⊙) তারকার বিবর্তন এবং দোলন গবেষণা
३. eAGB এবং TP-AGB তারকার পর্যবেক্ষণ পার্থক্য
१. বিভিন্ন ভর বা ধাতব প্রাচুর্যের তারকা: বিবর্তনীয় ক্রম পুনরায় গণনা প্রয়োজন
२. আরও তরুণ বা পুরানো ক্লাস্টার: প্রধান ক্রম ভর এবং বিবর্তনীয় সময়মাপ ভিন্ন
३. ক্ষেত্র তারকা (অ-ক্লাস্টার) চলক: ভর ক্ষয় হার ভিন্ন হতে পারে
१. নিম্ন ধাতব প্রাচুর্য গ্লোবুলার ক্লাস্টার (যেমন [Fe/H]<−1.5): দোলন প্রক্রিয়া ভিন্ন হতে পারে
२. বৃহৎ ভর তারকা (M>2M⊙): বিবর্তন পথ সম্পূর্ণ ভিন্ন
३. স্বল্প-পর্যায়ের চলক (যেমন RR Lyrae): বিভিন্ন দোলন অস্থিতিশীল অঞ্চল
१. Paxton et al. (२०१९): MESA কোড, Astrophys. J. Suppl. Ser. २४३, १०
२. Fadeyev (२०१३): তরল গতিশীলতা সমীকরণ, Astron. Lett. ३९, ३०६
३. Lebzelter & Wood (२००५): রৈখিক দোলন বিশ্লেষণ, Astron. Astrophys. ४४१, १११७
४. McDonald et al. (२०११): অবলোহিত পর্যবেক্ষণ এবং ধুলো, Astrophys. J. Suppl. Ser. १९३, २३
५. Willson (२०००): শকওয়েভ এবং ধুলো সংঘনন প্রক্রিয়া, Ann. Rev. Astron. Astrophys. ३८, ५७३
সামগ্রিক মূল্যায়ন: এটি একটি উচ্চ-মানের তারকা পদার্থবিজ্ঞান তাত্ত্বিক গবেষণা পেপার, কঠোর সংখ্যাসূচক গণনার মাধ্যমে ४७ টাকের দীর্ঘ-পর্যায়ের চলকের বিবর্তনীয় অবস্থার দীর্ঘমেয়াদী বিরোধ সমাধান করে। পদ্ধতি উদ্ভাবনী (বিবর্তন-দোলন সামঞ্জস্যপূর্ণ গণনা), সিদ্ধান্ত স্পষ্ট (বেশিরভাগ eAGB তারকা), শারীরিক চিত্র স্পষ্ট (ছোট ভর অস্থিতিশীলতা সৃষ্টি করে)। প্রধান অপূর্ণতা একক ভর গণনা এবং পর্যবেক্ষণের সাথে পরিসংখ্যানগত তুলনার অভাব। এই কাজ গ্লোবুলার ক্লাস্টার চলক গবেষণায় গুরুত্বপূর্ণ তাত্ত্বিক মূল্য রয়েছে, পদ্ধতি অন্যান্য জ্যোতির্বিজ্ঞান সমস্যায় প্রসারিত করা যায়।