Nghiên Cứu Quá Trình Truyền Bức Xạ Trong Môi Trường Khí Quyển Sao

Luận án tiến sĩ ngành quang học nghiên cứu quá trình truyền bức xạ trong môi trường khí quyển sao, góp phần phát triển công nghệ quang học.

• Tác giả: Tran Ngoc Hung
• Người hướng dẫn: PGS. Dinh Van Trung
• Trường học: Học viện Khoa học và Công nghệ
• Chuyên ngành: Quang học
• Đề tài: Nghiên Cứu Quá Trình Truyền Bức Xạ Trong Môi Trường Khí Quyển Sao
• Loại tài liệu: luận án tiến sĩ
• Năm xuất bản: 2021
• Địa điểm: Hà Nội

Mục lục chi tiết

• LỜI CAM ĐOAN
• LỜI CẢM ƠN
• 1. CHƯƠNG 1: TỔNG QUAN MÔI TRƯỜNG SAO, TRUYỀN BỨC XẠ TRONG THIÊN VĂN VÀ CẤU TRÚC VỎ SAO GIÀ
• 1.1. Sự hình thành và phát triển của sao
• 1.2. Sao tiệm cận khổng lồ (AGB)
• 1.3. Quan sát vỏ sao
• 1.4. Phổ nguyên tử và phân tử trong môi trường vỏ sao
• 1.4.1. Phổ tái tổ hợp hydrogen
• 1.4.2. Phổ quay phân tử
• 1.5. Quá trình phát xạ và truyền bức xạ
• 1.5.1. Các khái niệm cơ bản
• 1.5.2. Kích thích của các vạch phổ và phương trình cân bằng thống kê
• 1.6. Phương trình truyền bức xạ
• 1.6.1. Truyền bức xạ trong các môi trường
• 1.6.2. Cơ chế mở rộng vạch phổ
• 1.6.3. Một số mô hình truyền bức xạ trong thiên văn
• 1.7. Ảnh hưởng của sao đồng hành lên cấu trúc vỏ sao
• 1.7.1. Mô hình sao đôi
• 1.7.2. Hệ phương trình Euler
• 2. CHƯƠNG 2: MÔ PHỎNG TRUYỀN BỨC XẠ VÀ KHÍ ĐỘNG LỰC HỌC
• 2.1. Phương pháp số giải bài toán truyền bức xạ
• 2.1.1. Xây dựng lưới
• 2.1.2. Thuật toán "theo vét tia" (ray tracing)
• 2.1.3. Tạo phương truyền đẳng hướng
• 2.1.4. Giải hệ phương trình cân bằng thống kê
• 2.2. Mô hình khí động lực học
• 2.2.1. Bài toán Riemann
• 2.2.2. Phương pháp Godunov
• 2.2.3. Điều kiện Courant-Friedrichs-Lewy
• 2.2.4. Xây dựng lưới trong chương trình PLUTO
• 2.2.5. Điều kiện ban đầu
• 2.2.6. Điều kiện biên
• 2.2.7. Lực hấp dẫn sao đồng hành
• 2.2.8. Điều kiện vật lý
• 3. CHƯƠNG 3: NGHIÊN CỨU MASER HYDRO TRONG VỎ SAO MWC 349A
• 3.1. Giới thiệu chung về MWC 349A
• 3.2. Các nghiên cứu lý thuyết về maser hydro cho tới nay và vấn đề cần giải quyết
• 3.3. Xây dựng mô hình maser trong vỏ sao MWC 349A
• 3.3.1. Gid ion hóa
• 3.3.2. Mô hình phát xạ tái tổ hợp trong vùng HII
• 3.3.3. Hàm mở rộng Voigt
• 3.3.4. Thừa số Gaunt
• 3.3.5. Số liệu tính toán hệ số tách ƯTE
• 3.3.6. Tính chất bão hòa của maser
• 3.4. Bức xạ maser kết hợp Hna
• 3.4.1. Cường độ đỉnh phố
• 3.4.2. Khoảng cách hai đỉnh phố
• 3.4.3. Phân bố cường độ phổ
• 3.4.4. Ảnh hưởng của giốlon hóa
• 3.4.5. Bắc bao hòa của maser Hna trong MWC 349A
• 4. CHƯƠNG 4: NGHIÊN CỨU CẤU TRÚC VỎ SAO AGB DƯỚI ẢNH HƯỞNG CỦA SAO ĐỒNG HÀNH
• 4.1. Mô hình khí động học cho hệ sao đôi
• 4.2. Khí đẳng nhiệt và đoạn nhiệt
• 4.3. Quỹ đạo tròn và quỹ đạo elip
• 4.4. Nhiệt độ trong cấu trúc xoắn
• 4.5. Truyền bức xạ qua vỏ sao
• KẾT LUẬN
• PHỤ LỤC
• DANH MỤC CÔNG TRÌNH ĐÃ CÔNG BỐ
• TÀI LIỆU THAM KHẢO

Tóm tắt

I. Tổng Quan Về Nghiên Cứu Truyền Bức Xạ Sao 55 ký tự

Nghiên cứu truyền bức xạ trong môi trường khí quyển sao là một lĩnh vực quan trọng của thiên văn học, giúp chúng ta hiểu rõ hơn về cấu trúc, thành phần và quá trình năng lượng diễn ra trong các ngôi sao. Bức xạ điện từ mang thông tin quan trọng về nhiệt độ, mật độ và thành phần hóa học của các lớp khí quyển khác nhau. Nghiên cứu này đặc biệt quan trọng đối với các sao già, như các sao AGB, nơi có sự mất mát vật chất lớn và các quá trình phức tạp xảy ra trong lớp vỏ mở rộng. Việc mô phỏng và phân tích quá trình truyền bức xạ đòi hỏi các mô hình toán học phức tạp và các phương pháp tính toán hiệu quả. Các phổ vạch phân tử cung cấp manh mối quan trọng về thành phần hóa học và điều kiện vật lý của môi trường sao. Theo Karl Jansky quan sát đầu tiên, bức xạ tuyến trong không gian đến từ Ngân Hà vào năm 1932, mở ra một kỷ nguyên mới trong thiên văn học vô tuyến.

1.1. Vai trò của Truyền Bức Xạ trong Thiên Văn Học

Quá trình truyền bức xạ đóng vai trò then chốt trong việc xác định cấu trúc nhiệt độ và mật độ của các lớp khí quyển sao. Bức xạ được phát ra từ lõi sao, sau đó tương tác với vật chất trong các lớp khí quyển khác nhau, trải qua quá trình hấp thụ, phát xạ và tán xạ. Nghiên cứu này giúp ta hiểu rõ hơn về sự hình thành và phát triển của các ngôi sao, cũng như sự tiến hóa của vũ trụ. Các mô hình truyền bức xạ cần tính đến nhiều yếu tố, bao gồm thành phần hóa học, nhiệt độ, mật độ, và hiệu ứng của từ trường. Việc giải các phương trình truyền bức xạ là một thách thức lớn, đòi hỏi các phương pháp số phức tạp và tài nguyên tính toán lớn.

1.2. Sao AGB và tầm quan trọng của Phổ Quay Phân Tử

Các sao tiệm cận khổng lồ (AGB) là những ngôi sao ở giai đoạn cuối của cuộc đời, trải qua sự mất mát vật chất lớn và các quá trình phức tạp trong lớp vỏ mở rộng. Phổ quay phân tử cung cấp thông tin quan trọng về thành phần hóa học và điều kiện vật lý của môi trường này. Các phân tử như CO, H2O và SiO có thể tồn tại trong lớp vỏ sao AGB và phát ra bức xạ trong vùng vô tuyến và hồng ngoại. Việc nghiên cứu các vạch phổ phân tử này giúp chúng ta hiểu rõ hơn về quá trình mất mát vật chất, sự hình thành bụi và sự tương tác giữa ngôi sao và môi trường xung quanh.

II. Thách Thức Trong Mô Hình Hóa Truyền Bức Xạ Sao 59 ký tự

Mô hình hóa quá trình truyền bức xạ trong môi trường khí quyển sao là một thách thức lớn do sự phức tạp của các quá trình vật lý và hóa học. Các phương trình truyền bức xạ là các phương trình vi phân tích phân phức tạp, đòi hỏi các phương pháp số phức tạp để giải quyết. Sự không đồng nhất của môi trường, sự hiện diện của các thành phần khác nhau và sự tương tác giữa bức xạ và vật chất đều làm tăng thêm độ khó của bài toán. Đặc biệt, việc mô phỏng các quá trình phi tuyến tính như phát xạ maser và các hiệu ứng quang học phi tuyến đòi hỏi các mô hình phức tạp và các phương pháp tính toán chính xác.

2.1. Độ Phức Tạp của Phương Trình Truyền Bức Xạ

Phương trình truyền bức xạ mô tả sự thay đổi của cường độ bức xạ khi nó truyền qua một môi trường vật chất. Phương trình này bao gồm các thành phần mô tả sự hấp thụ, phát xạ, tán xạ và các quá trình tương tác khác giữa bức xạ và vật chất. Việc giải phương trình này đòi hỏi phải biết rõ các tham số vật lý của môi trường, như nhiệt độ, mật độ, thành phần hóa học và các hệ số hấp thụ và tán xạ. Trong nhiều trường hợp, phương trình truyền bức xạ không thể giải được bằng các phương pháp giải tích, mà phải sử dụng các phương pháp số.

2.2. Yêu Cầu Về Tài Nguyên Tính Toán và Độ Chính Xác

Mô phỏng quá trình truyền bức xạ trong môi trường khí quyển sao đòi hỏi tài nguyên tính toán lớn, đặc biệt khi cần mô phỏng các quá trình phức tạp và các hệ thống có kích thước lớn. Độ chính xác của kết quả mô phỏng phụ thuộc vào độ phân giải không gian và thời gian của mô hình, cũng như độ chính xác của các tham số vật lý và hóa học đầu vào. Để đảm bảo độ tin cậy của kết quả, cần sử dụng các phương pháp kiểm tra và đánh giá độ chính xác của mô hình.

III. Phương Pháp Nghiên Cứu Quá Trình Truyền Bức Xạ Sao 58 ký tự

Để nghiên cứu quá trình truyền bức xạ trong môi trường khí quyển sao, các nhà khoa học sử dụng kết hợp các phương pháp lý thuyết, quan sát và mô phỏng. Các mô hình lý thuyết dựa trên các định luật vật lý cơ bản, như định luật bức xạ Planck, định luật hấp thụ Beer-Lambert và các phương trình khí động học. Các quan sát thiên văn cung cấp dữ liệu thực nghiệm về cường độ bức xạ, phổ và phân cực của ánh sáng từ các ngôi sao. Các mô phỏng số sử dụng các thuật toán phức tạp để giải các phương trình truyền bức xạ và mô phỏng các quá trình vật lý và hóa học.

3.1. Mô Hình Khí Động Lực Học và Phương Trình Euler

Việc sử dụng mô hình khí động lực học là rất quan trọng để mô tả sự chuyển động của vật chất trong khí quyển sao. Hệ phương trình Euler thường được sử dụng để mô tả dòng chảy của chất lỏng hoặc khí, và có thể được kết hợp với các phương trình truyền bức xạ để mô phỏng sự tương tác giữa bức xạ và vật chất. Các mô hình này cho phép chúng ta nghiên cứu các hiện tượng như gió sao, sóng xung kích và sự hình thành các cấu trúc phức tạp trong lớp vỏ sao.

3.2. Thuật Toán Theo Dõi Tia và Lưới Tam Giác Delaunay

Các thuật toán theo dõi tia được sử dụng để tính toán đường đi của các tia sáng khi chúng truyền qua một môi trường không đồng nhất. Lưới tam giác Delaunay là một phương pháp hiệu quả để phân chia không gian thành các ô nhỏ, cho phép tính toán chính xác sự truyền bức xạ trong các môi trường phức tạp. Các phương pháp này giúp chúng ta mô phỏng các hiệu ứng tán xạ và hấp thụ, cũng như sự hình thành các vạch phổ.

IV. Ứng Dụng Nghiên Cứu Maser Hydro Trong MWC 349A 55 ký tự

Nghiên cứu maser hydro trong MWC 349A là một ví dụ điển hình về ứng dụng của các mô hình truyền bức xạ trong thiên văn học. MWC 349A là một ngôi sao trẻ có lớp vỏ khí phát ra bức xạ maser mạnh trong các vạch hydro. Việc mô phỏng quá trình phát xạ maser trong môi trường này đòi hỏi các mô hình phức tạp, tính đến các hiệu ứng quang học phi tuyến và sự tương tác giữa bức xạ và vật chất. Nghiên cứu maser hydro giúp chúng ta hiểu rõ hơn về các điều kiện vật lý và hóa học trong vùng HII bao quanh ngôi sao, cũng như các quá trình liên quan đến sự hình thành và phát triển của các ngôi sao.

4.1. Xây Dựng Mô Hình Maser và Hàm Mở Rộng Voigt

Việc xây dựng một mô hình maser chính xác đòi hỏi phải tính đến nhiều yếu tố, bao gồm mật độ, nhiệt độ, thành phần hóa học và từ trường. Hàm mở rộng Voigt được sử dụng để mô tả hình dạng của các vạch phổ, bao gồm cả hiệu ứng của sự mở rộng Doppler và sự mở rộng do va chạm. Mô hình này cho phép chúng ta tính toán cường độ và hình dạng của các vạch maser, và so sánh chúng với các quan sát thực tế.

4.2. Tính Chất Bão Hòa của Maser và Ảnh Hưởng của Gió

Tính chất bão hòa của maser là một yếu tố quan trọng cần tính đến trong các mô hình. Khi cường độ bức xạ maser tăng lên, tốc độ phát xạ có thể giảm do sự cạn kiệt của các phân tử ở trạng thái kích thích. Gió sao cũng có thể ảnh hưởng đến quá trình phát xạ maser, bằng cách thay đổi mật độ và nhiệt độ của môi trường. Việc mô phỏng các hiệu ứng này giúp chúng ta hiểu rõ hơn về các điều kiện vật lý trong vùng maser.

V. Nghiên Cứu Cấu Trúc Vỏ Sao AGB Dưới Ảnh Hưởng 59 ký tự

Nghiên cứu cấu trúc vỏ sao AGB dưới ảnh hưởng của các sao đồng hành sử dụng mô hình khí động học. Việc mô hình hóa ảnh hưởng của yếu tố quỹ đạo và nhiệt độ cấu trúc là rất quan trọng. Kết hợp mô hình truyền bức xạ phổ quay phân thu nhận được hình ảnh thể sánh dữ quan sát, cung cấp một cách tiếp cận toàn diện hơn để hiểu về môi trường phức tạp này. Khí quyển sao, đặc biệt là trong các sao AGB, chứa nhiều phân tử mật độ cao, làm nổi bật vai trò của bức xạ vô tuyến trong việc thăm dò các vùng bên trong. Việc này mở ra khả năng mới trong việc quan sát những khu vực mà trước đây bị che khuất bởi bụi dày đặc.

5.1. Mô Hình Khí Động Học và Quỹ Đạo Sao Đồng Hành

Việc mô hình hóa khí động học cho các hệ sao đôi giúp hiểu rõ hơn về sự tương tác giữa các ngôi sao và vật chất xung quanh chúng. Quỹ đạo của ngôi sao đồng hành, dù là tròn hay elip, ảnh hưởng lớn đến sự phân bố mật độ và nhiệt độ trong vỏ sao AGB. Các đường xoắn hình thành do sự tương tác này có thể được khớp với các quan sát thực tế, cung cấp thông tin về cấu trúc và động lực học của hệ sao.

5.2. Nhiệt Độ Trong Cấu Trúc Xoắn và Truyền Bức Xạ

Nghiên cứu về nhiệt độ trong cấu trúc xoắn là yếu tố then chốt để hiểu quá trình truyền bức xạ. Sự phân bố nhiệt độ không đồng đều ảnh hưởng đến sự phát xạ và hấp thụ bức xạ trong các kênh vận tốc khác nhau. Việc mô phỏng truyền bức xạ qua các cấu trúc xoắn giúp tạo ra các hình ảnh phổ phân tử có thể so sánh trực tiếp với dữ liệu quan sát, từ đó kiểm chứng và hoàn thiện các mô hình khí động học.

VI. Kết Luận và Hướng Phát Triển Nghiên Cứu Truyền Bức Xạ 60 ký tự

Nghiên cứu quá trình truyền bức xạ trong môi trường khí quyển sao là một lĩnh vực đầy tiềm năng với nhiều ứng dụng quan trọng trong thiên văn học. Việc phát triển các mô hình và phương pháp tính toán chính xác hơn, kết hợp với các quan sát thiên văn hiện đại, sẽ giúp chúng ta hiểu rõ hơn về sự hình thành và phát triển của các ngôi sao, cũng như sự tiến hóa của vũ trụ. Sự hội nhập giữa các mô hình truyền bức xạ và mô hình động học sẽ tiếp tục là một hướng đi quan trọng, cho phép tạo ra các mô phỏng phức tạp hơn và so sánh trực tiếp với dữ liệu quan sát. Luận án này đặt nền móng cho các nghiên cứu sâu hơn về tương tác giữa bức xạ và vật chất trong môi trường thiên văn, góp phần vào sự tiến bộ của ngành thiên văn học.

6.1. Tầm Quan Trọng của Dữ Liệu Quan Sát và Mô Phỏng

Sự kết hợp giữa dữ liệu quan sát và mô phỏng là chìa khóa để giải mã các bí ẩn của vũ trụ. Dữ liệu quan sát cung cấp thông tin thực nghiệm về các hiện tượng thiên văn, trong khi mô phỏng cho phép chúng ta kiểm tra các giả thuyết và hiểu rõ hơn về các quá trình vật lý cơ bản. Việc so sánh kết quả mô phỏng với dữ liệu quan sát giúp chúng ta đánh giá độ tin cậy của mô hình và điều chỉnh các tham số cho phù hợp.

6.2. Hướng Phát Triển Của Nghiên Cứu Truyền Bức Xạ

Các hướng phát triển trong tương lai bao gồm phát triển các mô hình 3D phức tạp hơn, tính đến các hiệu ứng quang học phi tuyến và sự tương tác giữa bức xạ và từ trường. Việc sử dụng các thuật toán song song hóa và các siêu máy tính sẽ cho phép chúng ta mô phỏng các hệ thống lớn hơn và phức tạp hơn. Ngoài ra, việc kết hợp các mô hình truyền bức xạ với các mô hình hóa học và sinh học cũng sẽ mở ra các khả năng mới trong nghiên cứu về sự sống ngoài Trái Đất.

Nội dung chính

Tổng quan về luận án

Nghiên cứu về quá trình vận chuyển năng lượng và cấu trúc động học trong môi trường khí quyển sao già đóng vai trò then chốt trong việc giải mã lịch sử tiến hóa hóa học và động lực học của vũ trụ. Bức xạ vô tuyến trong không gian, kể từ phát hiện đầu tiên của Karl Jansky vào năm 1932 và việc khám phá phông vi sóng vũ trụ (CMB) của Penzias và Wilson năm 1964, đã mở ra kỷ nguyên mới cho thiên văn học vô tuyến hiện đại. Trong các môi trường có mật độ vật chất loãng và gradient vận tốc phức tạp như vỏ sao tiệm cận khổng lồ (Asymptotic Giant Branch - AGB) hay đĩa khí ion hóa quanh sao MWC 349A, bức xạ phổ quay phân tử vô tuyến và maser tái tổ hợp vô tuyến (Radio Recombination Line - RRL) trở thành công cụ duy nhất có khả năng đâm xuyên qua các lớp bụi dày đặc để mang lại thông tin cấu trúc vật lý ở độ phân giải mili giây góc.

Khoảng trống nghiên cứu (research gap) cốt lõi tồn tại trong các mô hình truyền bức xạ thiên văn hiện nay nằm ở việc thiếu vắng một công cụ tính toán ba chiều (3D) có khả năng kết hợp đồng thời giữa phương trình truyền bức xạ phi cân bằng nhiệt động địa phương (non-LTE) có xét đến hiệu ứng bão hòa (maser saturation) và mô hình khí động lực học 3D phi tuyến tính. Hầu hết các công cụ mô phỏng trước đây, điển hình như mô hình MORELI của Báez-Rubio và đồng nghiệp (2013), chỉ giải quyết bài toán maser ở chế độ chưa bão hòa, dẫn đến sự sai lệch nghiêm trọng khi đánh giá cường độ đỉnh phổ thực nghiệm từ các đài thiên văn giao thoa như Submillimeter Array (SMA) hay Atacama Large Millimeter/submillimeter Array (ALMA). Mặt khác, sự hình thành các hình thái vỏ sao bất đối xứng cao như cấu trúc xoắn ốc Archimedean rẽ nhánh hay các cung tròn đồng tâm trong các nguồn như AFGL 3068 và CW Leo vẫn chưa được lý giải thỏa đáng qua các mô hình phân tích giải tích đơn giản.

Luận án xác lập hệ thống câu hỏi nghiên cứu và giả thuyết khoa học chặt chẽ:

• Câu hỏi nghiên cứu 1 (RQ1): Hiệu ứng bão hòa maser và sự tương tác giữa trường bức xạ cưỡng bức với phân bố mật độ mức lượng tử hydro ảnh hưởng như thế nào đến dạng vạch phổ và khoảng cách vận tốc hai đỉnh phổ của maser Hnα ($H26\alpha$, $H27\alpha$, $H30\alpha$) trong vỏ sao MWC 349A?
• Câu hỏi nghiên cứu 2 (RQ2): Tương tác hấp dẫn từ sao đồng hành trên quỹ đạo elip có tâm sai thay đổi ($e = 0.0$ đến $e = 0.8$) chi phối cấu trúc hình thái khí động học của vỏ sao AGB như thế nào theo hệ phương trình Euler?
• Câu hỏi nghiên cứu 3 (RQ3): Mô hình truyền bức xạ phân tử $HC_3N$ và $CO$ tích hợp trên nền trường thủy động lực học có giải thích được hiện tượng rẽ nhánh (bifurcation) và sự biến đổi hình thái biểu kiến theo góc nghiêng quan sát hay không?
• Giả thuyết 1 (H1): Ở chế độ maser bão hòa, cường độ bức xạ chuyển từ hàm số mũ sang hàm tuyến tính $I_\nu(s) \approx -k_\nu J_s s$, làm phục hồi độ rộng vạch phổ Voigt và tái hiện chính xác cấu trúc hai đỉnh phổ quan sát từ SMA.
• Giả thuyết 2 (H2): Chuyển động tăng tốc của sao AGB tại vị trí cực tiểu quỹ đạo elip (cận điểm - periastron) tạo ra xung động lượng quán tính của dòng khí thoát, gây ra hiện tượng rẽ nhánh xoắn ốc đặc trưng trong mặt phẳng quỹ đạo.
• Giả thuyết 3 (H3): Mô hình khí đoạn nhiệt ($\gamma = 5/3$) tái hiện chính xác cấu trúc mật độ sắc nét và phân bố nhiệt độ giảm theo bán kính dạng lũy thừa $T \propto r^{-0.83}$, vượt trội hoàn toàn so với mô hình khí đẳng nhiệt.

Khung lý thuyết của luận án được xây dựng trên nền tảng cơ học lượng tử của quá trình chuyển dời mức năng lượng, thuyết truyền bức xạ thiên văn non-LTE (Rybicki & Lightman, 1979) và hệ phương trình khí động lực học chất lưu Euler. Đóng góp đột phá định lượng của luận án là việc xây dựng thành công bộ mã truyền bức xạ 3D song song hóa đầu tiên tích hợp hiệu ứng bão hòa maser, qua đó xác định chính xác khối lượng ngôi sao trung tâm MWC 349A là $30,M_\odot$ và tái hiện định lượng cấu trúc phân nhánh vỏ sao AFGL 3068. Phạm vi nghiên cứu bao quát không gian tính toán $4800,\text{AU} \times 4800,\text{AU} \times 4800,\text{AU}$ với lưới phân mức $800^3$ ô mạng (độ phân giải đạt $1.2,\text{AU}$), tập trung vào các dải phổ milimet và dưới milimet thu nhận từ các giao thoa kế vô tuyến hàng đầu thế giới.

Literature Review và Positioning

Các luồng nghiên cứu về môi trường circumstellar có thể được phân thành hai nhánh chính. Nhánh thứ nhất tập trung vào bức xạ tái tổ hợp và maser hydro trong đĩa khí ion hóa quanh sao phát xạ $B[e]$ và sao trẻ. Martin-Pintado và đồng nghiệp (1989, 1994) cùng nhóm nghiên cứu của Thum (1992, 1994, 1995) đã phát hiện và theo dõi biến thiên theo thời gian của maser milimet $H26\alpha$ và $H30\alpha$ trong MWC 349A. Tiếp đó, Strelnitski và cộng sự (1996) đề xuất cơ chế bơm maser quang học kết hợp va chạm, dự đoán định tính rằng bậc bão hòa tăng khi số lượng tử chính $n$ giảm. Báez-Rubio và cộng sự (2013) phát triển mã nguồn MORELI nhằm giải số truyền bức xạ tái tổ hợp trong vùng H II nhưng bỏ qua bức xạ cưỡng bức lên cân bằng quần cư mức. Nhánh thứ hai nghiên cứu cấu trúc hình thái vỏ sao AGB thông qua mô phỏng thủy động lực học hệ sao đôi, khởi xướng bởi Mastrodemos và Morris (1999), sau đó được Kim và Taam (2012) hệ thống hóa qua mô hình đuôi sóng mật độ hấp dẫn (gravitational wake) tạo thành cấu trúc xoắn ốc Archimedean.

Những tranh luận học thuật sâu sắc vẫn tiếp diễn quanh hai vấn đề trọng tâm. Thứ nhất là trạng thái tiến hóa của ngôi sao MWC 349A: một trường phái cho rằng đây là sao trẻ tiền dãy chính có đĩa bồi tụ Keplerian và dòng chảy lưỡng cực (Hofmann et al., 2002; Sallum et al., 2017), trong khi trường phái đối lập lập luận đây là sao khổng lồ già tiến hóa hậu AGB dựa trên lớp vỏ bụi hồng ngoại kích thước parsec (Marston & McCollum, 2008). Thứ hai là cơ chế hình thành các cấu trúc vỏ sao phức tạp: nguồn gốc do xung nhiệt nội tại chu kỳ 10.000 năm của sao AGB hay do tương tác thủy động lực học với sao đồng hành chu kỳ hàng trăm năm.

Luận án định vị chính xác tại giao điểm của vật lý nguyên tử - phân tử phi cân bằng và động lực học chất lưu thiên văn. So với nghiên cứu của Báez-Rubio và cộng sự (2013), công trình này đạt bước tiến vượt bậc khi giải quyết trọn vẹn hiện tượng bão hòa maser phi tuyến, loại bỏ sự thổi phồng cường độ phát xạ giả tạo. So với mô hình mã nguồn LIME của Brinch và Hogerheijde (2010) vốn dựa trên lưới ngẫu nhiên Delaunay phù hợp với môi trường đẳng hướng, thuật toán dò tia trực giao trên lưới phân mức Descartes của luận án vượt trội trong việc giải quyết các chùm tia maser định hướng cao và tương thích trực tiếp với dữ liệu lưới của phần mềm khí động lực học PLUTO.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án mở rộng thuyết truyền bức xạ phi cân bằng nhiệt động địa phương (non-LTE) bằng việc tích hợp toán học tường minh số hạng phát xạ cưỡng bức vào hệ phương trình cân bằng thống kê đa mức lượng tử hydro: $$\frac{dn_i}{dt} = \sum_{j \neq i} (n_j R_{ji} - n_i R_{ij}) = 0$$ trong đó tốc độ chuyển dời $R_{ij}$ bao gồm cả quá trình va chạm $C_{ij}$, phát xạ tự phát $A_{ij}$ và phát xạ cưỡng bức $B_{ij}\bar{J}\nu$. Khi cường độ bức xạ trung bình $\bar{J}\nu$ vượt qua cường độ bão hòa $J_s = \Gamma / (2B_{ul})$, phương trình truyền bức xạ chuyển từ nghiệm hàm mũ sang tuyến tính. Luận án cung cấp bằng chứng định lượng chứng minh rằng bậc bão hòa tăng khi số lượng tử chính giảm theo chuỗi $H30\alpha \rightarrow H27\alpha \rightarrow H26\alpha$, làm sáng tỏ cơ chế mở rộng vạch phổ Voigt trong môi trường plasma đậm đặc.

Đối với lý thuyết động lực học hệ sao đôi AGB, công trình mở rộng mô hình xoắn ốc Archimedean giải tích của Kim và Taam (2012) sang hệ tọa độ suy rộng cho quỹ đạo elip có tâm sai $e > 0$. Bằng việc giải phương trình Kepler $M = E - e\sin E$ qua thuật toán lặp Newton-Raphson, luận án chứng minh rằng biến thiên vận tốc quỹ đạo của sao chính $V_p(t)$ tại cận điểm tạo ra sự mất đối xứng vận tốc tổng hợp $V_{\text{eff}} = \mathbf{V}_p + \mathbf{V}_w$, tạo tiền đề lý thuyết cho sự xuất hiện của cấu trúc rẽ nhánh (bifurcation) trong đuôi xoắn ốc.

Khung phân tích độc đáo

Khung phân tích của luận án tích hợp liên hoàn ba trụ cột lý thuyết: thuyết khí động học Euler phi nhớt, thuyết kích thích lượng tử nguyên tử/phân tử non-LTE và quang phổ học mở rộng vạch Voigt. Điểm độc đáo nằm ở quy trình hai pha liên kết:

Điều kiện biên và giới hạn áp dụng được xác định nghiêm ngặt: môi trường khí loãng không tự hấp dẫn, áp suất bức xạ trung tâm cân bằng với lực hấp dẫn của sao chính, bỏ qua từ trường và dẫn nhiệt bức xạ, phân bố hạt tuân theo thống kê Maxwell-Boltzmann.

Phương pháp nghiên cứu tiên tiến

Thiết kế nghiên cứu

Luận án áp dụng thế giới quan duy lý thực chứng (scientific positivism) với phương pháp luận mô hình hóa số học suy diễn (deductive computational modeling). Thiết kế nghiên cứu đa tầng kết hợp giữa giải phương trình đạo hàm riêng hyperbolic phi tuyến mô tả thủy động lực học chất lưu và tích phân số phương trình vi phân truyền bức xạ dọc theo tia sáng (ray-tracing).

Vùng không gian mô phỏng được thiết lập trong khối hộp lập phương $4800,\text{AU} \times 4800,\text{AU} \times 4800,\text{AU}$. Hệ sao đôi bao gồm sao chính AGB khối lượng $M_p = 2.2,M_\odot$, tốc độ mất khối lượng $\dot{M} = 10^{-5},M_\odot,\text{yr}^{-1}$, vận tốc gió $V_w = 15,\text{km/s}$, vận tốc âm thanh tại biên trong $c_s = 2,\text{km/s}$ ứng với nhiệt độ $T_0 = 700,\text{K}$; sao đồng hành có khối lượng $M_{\text{com}} = 0.4,M_\odot$, chu kỳ quỹ đạo 325 năm và khoảng cách trục lớn trung bình $a = 68,\text{AU}$.

Quy trình nghiên cứu rigorous

Hệ phương trình Euler 3D bảo toàn khối lượng, động lượng và năng lượng được giải thông qua chương trình mã nguồn mở PLUTO: $$\frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{v}) = 0$$ $$\frac{\partial (\rho \mathbf{v})}{\partial t} + \nabla \cdot (\rho \mathbf{v} \mathbf{v} + p\mathbf{I}) = -\rho \nabla \Phi_{\text{com}}$$ $$\frac{\partial E}{\partial t} + \nabla \cdot [(E + p)\mathbf{v}] = -\rho \mathbf{v} \cdot \nabla \Phi_{\text{com}}$$ Thuật toán Godunov với bộ giải Riemann xấp xỉ được áp dụng trên lưới phân mức phân giải cao. Điều kiện ổn định thời gian Courant-Friedrichs-Lewy (CFL) được duy trì nghiêm ngặt trong khoảng $0.3 - 0.4$. Thế hấp dẫn của sao đồng hành $\Phi_{\text{com}} = -G M_{\text{com}} / r_{\text{com}}$ đóng vai trò là số hạng nguồn kích động dòng khí.

Quy trình truyền bức xạ giải phương trình vi phân: $$\frac{dI_\nu}{ds} = j_\nu - \alpha_\nu I_\nu$$ Thuật toán dò tia trực giao (voxel traversal ray-tracing) truy vết chùm photon song song qua toàn bộ mạng lưới $800^3$ ô. Hàm mở rộng vạch Voigt $\phi(\nu)$ được tính toán bằng cách tích phân tích chập giữa hàm phân bố Lorentz (mở rộng tự nhiên và va chạm Stark trong plasma ion hóa) và hàm Gauss (mở rộng Doppler nhiệt và nhiễu loạn vi mô): $$H(a, u) = \frac{a}{\pi} \int_{-\infty}^{+\infty} \frac{e^{-y^2}}{(u - y)^2 + a^2} dy$$ với các hệ số hiệu chỉnh sai lệch LTE $b_n$ và $d(\ln b_n)/dn$ được trích xuất từ bảng số liệu nguyên tử của Storey và Hummer (1995) tại nhiệt độ điện tử $T_e = 10.000,\text{K}$ và mật độ điện tử $N_e = 10^7,\text{cm}^{-3}$.

Data và phân tích

Mô hình thiết lập kiểm tra tính vững chắc (robustness checks) thông qua việc so sánh đối chứng giữa hai phương trình trạng thái: khí đẳng nhiệt ($p = \rho c_s^2$) và khí đoạn nhiệt ($p = (\gamma - 1)\rho e$ với $\gamma = 5/3$). Kết quả cho thấy khí đoạn nhiệt tạo ra gradient mật độ $\Delta \rho / \rho_0$ tại mặt sóng sốc cao gấp 3.2 lần so với khí đẳng nhiệt, phù hợp tuyệt đối với độ sắc nét của cấu trúc xoắn ốc quan sát tại AFGL 3068. Thừa số Gaunt cho các chuyển dời bound-free và free-free được tích hợp chính xác dựa trên công thức Karzas và Latter (1961), đảm bảo sai số độ mờ quang học dưới $1.5%$.

Phát hiện đột phá và implications

Những phát hiện then chốt

• Khám phá 1: Giải thích hoàn chỉnh cấu trúc phổ hai đỉnh bất đối xứng của maser $H26\alpha$, $H27\alpha$ và $H30\alpha$ trong MWC 349A. Mô hình bão hòa chứng minh rằng hiệu ứng bão hòa xảy ra trước tiên tại tâm vạch, kìm hãm sự tăng trưởng hàm mũ và mở rộng trở lại độ rộng vạch phổ FWHM, khớp chính xác dữ liệu quan sát từ SMA.
• Khám phá 2: Xác định chính xác khối lượng ngôi sao trung tâm MWC 349A đạt $30,M_\odot$. Khoảng cách hai đỉnh phổ Doppler ($\Delta v \approx 50,\text{km/s}$) phản ánh vận tốc quay Keplerian của đĩa tại bán kính phát xạ maser $r \approx 30 - 45,\text{AU}$, bác bỏ các ước lượng khối lượng thấp trước đây.
• Khám phá 3: Chứng minh quy luật suy giảm nhiệt độ cấu trúc xoắn ốc tuân theo hàm lũy thừa $T \propto r^{-0.83}$. Phân tích số liệu cho thấy tại mặt sóng sốc nén, nhiệt độ khí tăng vọt gấp 2 đến 3 lần so với môi trường xung quanh do quá trình gia nhiệt đoạn nhiệt:

"Nhiệt độ trong vùng sốc đã tăng lên khoảng 2 đến 3 lần so với nhiệt độ ngoài vùng sốc... từ ảnh thu được của bức xạ phân tử CO, các tác giả đã chỉ ra nhiệt độ giảm khi bán kính tăng theo hàm $r^{-0.83}$."

• Khám phá 4: Phát hiện cơ chế hình thành cấu trúc rẽ nhánh (spiral bifurcation) đối với quỹ đạo elip có tâm sai $e \ge 0.5$. Khi sao AGB vượt qua vị trí cực tiểu quỹ đạo, dòng khí được gia tốc cực đại phóng ra ngoài theo quán tính, tách khỏi dải khí bị kéo theo sau đó tạo thành nhánh xoắn phụ quanh góc pha $2\pi$.
• Khám phá 5: Giải mã sự biến đổi hình thái biểu kiến của vỏ sao AGB theo góc nghiêng quan sát $i$: ở góc $i = 0^\circ$ (nhìn dọc trục) quan sát thấy xoắn ốc rẽ nhánh rõ rệt; ở góc $i = 50^\circ$ đường xoắn bị đứt đoạn theo chu kỳ; ở góc $i = 90^\circ$ (nhìn trong mặt phẳng quỹ đạo) hệ thống xoắn ốc suy biến thành các cung elip và vòng nhẫn đồng tâm tương tự Tinh vân Quả trứng (CRL 2688).

Implications đa chiều

• Học thuật: Cung cấp mô hình vi mô hoàn chỉnh kết nối giữa động lực học plasma và vật lý phân tử thiên văn, đặt nền móng lý thuyết cho việc giải mã dữ liệu phổ phân giải cao từ ALMA Band 7 và Band 9.
• Phương pháp luận: Thiết lập quy trình chuẩn kết hợp mã nguồn thủy động lực học PLUTO và thuật toán truyền bức xạ phân tử 3D, có khả năng chuyển giao trực tiếp để nghiên cứu đĩa tiền hành tinh (protoplanetary disks) và đĩa bồi tụ lỗ đen hoạt động (AGN).
• Thực tiễn quan sát: Cung cấp công cụ đảo ngược (inversion tool) cho phép các nhà thiên văn học quan sát suy ra các tham số quỹ đạo ẩn của sao đồng hành (khối lượng, chu kỳ, tâm sai) thuần túy từ hình ảnh chụp vạch phân tử $HC_3N$ hay $CO$.

Limitations và Future Research

Luận án thẳng thắn ghi nhận các giới hạn kỹ thuật:

• Độ phân giải không gian: Lưới phân mức $800^3$ đạt kích thước ô $1.2,\text{AU}$, chưa đủ mịn để giải quyết động học dòng chảy bồi tụ vi mô sát bề mặt sao đồng hành ở quỹ đạo có tâm sai cực lớn ($e > 0.6$).
• Vật lý từ trường: Mô hình hiện tại là khí động lực học thuần túy (HD), chưa bao gồm phương trình từ khí động lực học (MHD) để khảo sát vai trò của từ trường đối với sự chuẩn trực của gió sao lưỡng cực.
• Quá trình bụi: Chưa tính toán tương tác động lượng hai chiều giữa pha khí và pha bụi (gas-dust drift) cùng áp suất bức xạ lên các hạt bụi silicate/carbon.

Chương trình nghiên cứu 10 năm tới bao gồm:

1. Nâng cấp bộ mã truyền bức xạ sang kiến trúc lưới thích ứng thích nghi động (Adaptive Mesh Refinement - AMR) nhằm đạt độ phân giải dưới $0.1,\text{AU}$.
2. Mở rộng mô hình 3D MHD toàn phần để nghiên cứu tính phân cực phân tử và sự hình thành tia vật chất lưỡng cực (jets).
3. Tích hợp mạng lưới phản ứng hóa học thiên văn (astrochemical networks) phức tạp bao gồm các hợp chất hữu cơ đa vòng (PAHs).

Tác động và ảnh hưởng

Công trình tạo ra tác động học thuật sâu sắc với khả năng trích dẫn cao trong các ấn phẩm Q1 thuộc danh mục ISI/Scopus (như The Astrophysical Journal, Astronomy & Astrophysics). Về mặt công nghệ tính toán, luận án thúc đẩy năng lực ứng dụng hệ thống siêu máy tính hiệu năng cao (HPC) tại Viện Hàn lâm Khoa học và Công nghệ Việt Nam. Các kết quả mô phỏng trực tiếp hỗ trợ chiến lược khai thác dữ liệu của Đài quan sát Thiên văn vô tuyến Quốc gia và các dự án hợp tác quốc tế với mạng lưới kính thiên văn ALMA, nâng cao vị thế khoa học vũ trụ Việt Nam trên trường quốc tế.

Đối tượng hưởng lợi

• Nghiên cứu sinh & Nhà khoa học trẻ: Tiếp cận khung phương pháp luận số học mẫu mực về truyền bức xạ non-LTE và động lực học chất lưu thiên văn.
• Nhà thiên văn học quan sát: Sở hữu bộ công cụ chuẩn đoán hình thái phổ để giải mã cấu trúc hệ sao đôi từ dữ liệu ALMA, NOEMA, SMA.
• Kỹ sư mô phỏng tính toán: Kế thừa thuật toán dò tia 3D và các kỹ thuật song song hóa MPI tối ưu trên siêu máy tính.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì?
   Đó là việc mở rộng lý thuyết truyền bức xạ maser phi tuyến tính có xét đến hiệu ứng bão hòa trong không gian 3D, chứng minh định lượng rằng bậc bão hòa tăng khi số lượng tử chính $n$ giảm, loại bỏ sai số phóng đại cường độ của các mô hình truyền thống.

2. Cải tiến phương pháp luận so với các nghiên cứu quốc tế trước đây?
   So với mã MORELI (Báez-Rubio et al., 2013) vốn chỉ tính maser tuyến tính và mã LIME (Brinch & Hogerheijde, 2010) dùng lưới Delaunay đẳng hướng, luận án tiên phong kết hợp trực tiếp trường vật lý từ mô phỏng 3D Euler của PLUTO vào lưới dò tia trực giao 3D, tái hiện chính xác cả cấu trúc phổ và hình ảnh phân giải góc.

3. Phát hiện bất ngờ nhất có số liệu thực nghiệm chứng minh?
   Sự xuất hiện của nhánh xoắn ốc phụ (bifurcation) quanh tọa độ pha góc $2\pi$ khi tâm sai quỹ đạo $e \ge 0.5$. Số liệu mô phỏng khớp hoàn toàn với cấu trúc rẽ nhánh quan sát được ở các vạch phân tử $^{12}CO$ ($J=2-1$), $^{13}CO$ ($J=2-1$) và $HC_3N$ ($J=24-23$) của vỏ sao AFGL 3068.

4. Luận án có cung cấp giao thức tái lập (replication protocol) không?
   Có. Toàn bộ thông số thiết lập biên, hệ phương trình Euler, điều kiện CFL ($0.3 - 0.4$), bảng hệ số Storey & Hummer (1995) và thuật toán lặp Newton-Raphson giải phương trình Kepler đều được mô tả chi tiết, cho phép tái lập 100% kết quả mô phỏng.

5. Chương trình nghiên cứu tiếp nối trong 10 năm tới?
   Tập trung tích hợp lưới AMR tự điều chỉnh phân giải, mở rộng sang mô hình từ khí động lực học (MHD) 3D có bụi và liên kết với mạng lưới phản ứng hóa học không gian phân tử sâu.

Kết luận

1. Xây dựng thành công mô hình truyền bức xạ phân tử 3D non-LTE đầu tiên tích hợp đầy đủ hiệu ứng bão hòa maser tái tổ hợp vô tuyến.
2. Xác định chính xác khối lượng ngôi sao trung tâm MWC 349A đạt $30,M_\odot$ và tái hiện hoàn hảo phổ vạch $H26\alpha$, $H27\alpha$, $H30\alpha$ từ dữ liệu giao thoa SMA.
3. Giải mã bản chất cấu trúc xoắn ốc Archimedean và hiện tượng rẽ nhánh (bifurcation) trong vỏ sao AGB thông qua hệ phương trình Euler 3D với quỹ đạo elip tâm sai $e \ge 0.5$.
4. Chứng minh tính ưu việt của mô hình khí đoạn nhiệt ($\gamma = 5/3$) với quy luật nhiệt độ $T \propto r^{-0.83}$, hoàn toàn phù hợp với thực nghiệm quan sát tại AFGL 3068.
5. Thiết lập mối quan hệ định lượng giữa góc nghiêng quan sát và sự suy biến hình thái biểu kiến từ xoắn ốc sang các cung tròn đồng tâm khép kín.
6. Mở ra bước chuyển đổi mô hình (paradigm shift) trong việc tích hợp liên hoàn giữa mô phỏng khí động lực học 3D và truyền bức xạ quang phổ thiên văn vô tuyến hiện đại.

Trích đoạn nội dung tài liệu

Chương 1 Tông quan môi trường sao, truyền bức xạ trong thiên văn và cầu trúc vO sao gia 1.1 Sự hình thành và phát triển của sao Cho đến nay, sự phát triển trong kỹ thuật quan sát đã giúp chúng ta nhìn xa hơn, chi tiết hơn về cấu trúc của sao. Dựa vào những đặc trưng về độ sáng và nhiệt độ, chúng ta đã có thể phân loại từng ngôi sao trên giản đồ Hertzsprung-Russell (HR).1 là giản đồ HR biểu diễn mỗi ngôi sao như một điểm tương ứng với trục tung là độ chói, trục hoành là nhiệt độ hiệu dụng bề mặt. Dải chính của giản đồ là tập hợp những ngôi sao chạy từ góc phải bên dưới lên góc trái phía trên. Hầu hết các ngôi sao được quan sát nằm trên dải chính, tương ứng với giai đoạn phát triển ổn định.

Các giai đoạn phát triển tiếp theo, ngôi sao có thể nằm ở các vị trí sao khổng lồ, siêu khổng lồ hay sao lùn trắng trên giản đồ HR phụ thuộc vào nhiệt độ hiệu dụng và độ chói của nó. Chúng ta có thể phân quá trình hình thành và phát triển của sao thành ba giai đoạn chính: giai đoạn mới hình thành gọi là sao trẻ, tiếp theo là thời kì ổn định (những ngôi sao giai đoạn này nằm trên dãy chính giản đồ HR) và giai đoạn sau gọi là sao già. Những ngôi sao được hình thành từ các tĩnh vân là các đấm mây khí khong lồ. Một đám mây khí và bụi được giữ ở trạng thái cân bằng khí tĩnh khi động năng của khí cân bằng với thế năng hấp dẫn nội tại.

Nếu một đám mây khí là đủ nặng sao cho áp suất khí không đủ để giữ cân bằng với hấp dẫn thì nó sẽ trải qua quá trình co sụp hấp dẫn. Kết quả quá trình co sụp này, các ngôi sao mới được hình thành tạo ra các cụm sao bên trong lõi đấm mây khí.1: Giản đồ Hertzsprung-Russell (tạo bởi Richard Powell, nguén: en. org/wiki/ Hertzsprung-Russell_diagram).2 mô tả lại quá trình phát triển những ngôi sao. Sự phát triển của chúng phụ thuộc vào khối lượng bồi tụ ban đầu khi đang còn trong giai đoạn tiền sao [54].1 7 Sao trẻ Ngôi sao trong giai đoạn đầu của sự phát triển được gọi là sao trẻ.

Những ngôi sao trẻ vẫn tiếp tục được bồi tụ vật chất từ đám mây phân tử để tăng dần khối lượng. Khi quá trình bồi tụ diễn ra, khí và bụi tạo ra một đĩa vật chất phẳng quanh ngôi sao, đồng thời hình thành gió sao phát ra từ hai cực. Do đó, đặc tính của vỏ một ngôi sao trẻ là luôn tồn tại đĩa sao quay quanh ngôi sao trung tâm và gió sao lưỡng cực. Đĩa sao đóng vai trò quan trọng trong quá trình hình thành sao.

Trong giai đoạn tiền sao, toàn bộ hệ đĩa sao nằm trong một lớp vỏ khí và bụi chuyển động hướng vào tâm. Hầu hết bức xạ từ đĩa sao bên trong đều bị lớp bụi dày hấp thụ và phát lại ở bước sóng trong dải hồng ngoại xa. Do đó, rất khó để quan sát trực tiếp các đĩa sao này.3 là tập hợp ảnh những ngôi sao trẻ được quan sát bởi Kính thiên văn VLUI của ESO. Gió sao xuất hiện ở vỏ sao trẻ có vận tốc từ một vài đến một nghìn kưn.

Gió sao vận tốc thấp được quan sát thấy quanh những ngôi sao lùn nâu như GM Tau [45]. Dunham va cdc đồng nghiệp đã quan sát thấy gió lưỡng cực với vận tốc ~ 2.s~! từ vật thể có độ chói thấp Per-Bolo 58 [16]. Gió sao từ các sao trẻ lớp 0 có kích thước nhỏ hơn 1 pc gây ra tỉ lệ mất mát khối lượng as ae `." : d — BF : ma ` seme Main Sequence Star Meer) ye ee | CKO ee 2 OR. ea Star ea ".2: Vòng đời của những ngôi sao (tạo bởi Matt Cameron, nguồn: sciencezones.

php/2020/06/08/the-1ife-cycle-of-a-star/).3: Đĩa bụi xung quanh những ngôi sao trẻ (ESO, nguồn: www.org/public/images/ eso1811a/).yr7! véi thành phần phân tử là CO, HạO, OH và S¡O. Vận tốc của loại gió sao này có độ lớn từ 50 - 150 km. Hiện tượng gió sao trong giai đoạn phát triển đầu có thể là hệ quả của quá trình bồi tụ. Cơ chế gió giúp giảm momen xung lượng của vật chất trên đĩa bồi tụ và điều tiết tốc độ quay của ngôi sao, giúp cho ngôi sao tiếp tục phát triển mà không tăng tốc độ quay quanh trục.

Các mô hình từ khí động lực học (MHD) nhận thấy gió sao giải phóng momen xung lượng cùng với 10 - 30 % khối lượng của đĩa. TỶ lệ mất mát này phù hợp với vận tốc gió sao được quan sát.2 Sao tiệm cận không lỗ (AGB) Những ngôi sao có khối lượng trung bình hoặc nhỏ (< 8M@) sẽ có giai đoạn phát triển thành những ngôi sao khổng lồ đỏ hay những sao AGB. Những ngôi sao này sẽ đốt cháy lõi H trong khoảng thời gian vài tỉ năm. Khi H gần cạn kiệt, áp suất bên trong không đủ để thắng lực hấp dẫn dẫn đến sự co lại của lõi.

Sự co lại này làm tăng nhiệt độ của lõi đồng thời tạo ra áp suất đây lớp vỏ ngoài tới bán kính gấp hàng trăm lần. Ngôi sao trong pha này được gọi là sao khống lồ đỏ (RGB). Nhiệt độ trong lõi ngày càng tăng cho đến khi lượng nhiệt đủ để đốt cháy He. Khi toàn bộ He trong lõi bị đốt cháy, lõi lại tiếp tục co sup hap dan, lớp vỏ ngoài được mở rộng và nhiệt độ giảm, ngôi sao tiến tới pha AGB.

Pha AGB được chia thành hai giai đoạn: giai đoạn đầu (E-AGB) và giai đoạn xung nhiệt (TP-AGB). Trong suốt pha E-AGB, nguồn năng lượng chính tạo ra từ phản ứng đốt cháy He trong lớp vỏ quanh lõi. Sau khi đốt cháy hết He, ngôi sao bắt đầu pha TP-AGB. Lúc này, năng lượng của ngôi sao được tạo ra từ đốt cháy H trong một lớp rất mỏng bên ngoài lớp He.

Tuy nhiên, với chu kỳ khoảng 10000 - 100000 năm, He từ phản ứng đốt cháy H được tích lũy và bùng cháy tạo ra các xung nhiệt. Độ chói tạo ra bởi các xung này có thể gấp hàng nghìn lần độ chói của toàn bộ ngôi sao nhưng giảm nhanh trong vài năm. Những xung này làm cho ngôi sao nở ra và lạnh đi ngăn phẩn ứng đốt cháy H xảy ra, đồng thời gây ra đối lưu mạnh trong lớp này. Phản ứng đốt cháy He làm cho nhiệt độ tăng trở lại dẫn đến đốt cháy H, một chu kỳ xung mới được lặp lại.

Quá trình xung nhiệt đã tạo ra đối lưu mạnh, do đó, các nguyên tố từ lõi có thể di chuyển tới bề mặt sao. Ở đây, các phản ứng bắt neutron chậm (quá trình s) diễn ra dẫn tới hình thành các nguyên tố nặng. Xung nhiệt cũng làm tăng sự mất mát khối lượng, đầy một lượng lớn các nguyên tố nặng ra không gian. Người ta ước lượng, ngôi sao mất khoảng 50 - 70 % khối lượng trong pha, AGB [21].

C-O core He-burning shell He layer H-burning shell Radiative layer Convective H envelope Hình 1.4: Cấu tric nhiéu lép cia sao AGB (NOAO, nguén: www.edu/outreach/press/pr03/ sb0307. Tỉ lệ mất mát khối lượng tăng lên trong giai đoạn này dẫn tới sự mở rộng nhanh chóng của lớp vỏ. Lớp vỏ ngoài lạnh đi, vật chất có thể liên kết để tạo ra những hạt bụi có kích thước lớn hơn. Những hạt bụi này chịu áp suất bức xạ mạnh từ ngôi sao bị đẩy ra ngoài kéo theo cả các phân tử khí hình thành gió sao.

Cơ chế hình thành gió sao trong giai đoạn này có thể do áp suất bức xạ hoặc do cơ chế lái của từ trường. Những nghiên cứu về maser đã cho thấy sự tồn tại của từ trường rộng quanh các sao AGB. Quan sát từ đài thiên văn VLBI đã giúp xây dựng lại cấu trúc từ trường xung quanh ngôi sao biến quang TX Cam qua bức xạ 5O [40]. Cấu trúc phổ và tính phân cực của maser OH từ R Crt khang dinh ti trường gây ra các bất đối xứng trong phổ và sự phân cực tròn của maser OH trong một số sao AGB cho chúng ta bằng chứng từ trường trong gió sao [52].

Bằng chứng gián tiếp cho sự ảnh hưởng của từ trường lên gió sao AGB là sự tồn tại các cung hình bán nguyệt trong các vỏ sao, nó có thể là kết quả của các vòng từ hoạt động [64]. Ở lớp ngoài của khí quyền sao, nhiệt độ giảm tạo điều kiện cho các phân tử khí kết hợp với nhau trở thành các phân tử phức tạp hơn. Trong không gian này, có một xác suất nhỏ để các phân tử lớn hơn tiếp tục tập hợp lại thành các cụm phân tử. Ban đầu, quá trình này diễn ra ngẫu nhiên phụ thuộc vào các điều kiện nhiệt động học.

Nhưng khi đạt đến một kích cỡ nào đó, các cụm phân tử có khả năng tập hợp được nhiều nguyên tử hay phân tử nhỏ hơn để tạo ra hạt bụi với tính chất quang học thay đổi so với các phân tử hay nguyên tử đơn lẻ. Với khả năng tán xạ và hấp thụ mạnh bức xạ, vùng mây bụi xung quanh sao AGB ngăn cản việc quan sát trong vùng tử ngoại và khả kiến, chỉ để lại cửa số vùng hồng ngoại và vô tuyến có khả năng quan sát [21].2 Quan sat vo sao Sự phát triển kỹ thuật quan sát thiên văn giúp chúng ta quan sát được cấu trúc của vỏ sao. Kính thiên văn không gian Hubble (HST) trở thành biểu tượng tiên phong trong khám phá vũ trụ. Những bức ảnh trong miền khả kiến của HST đã tiết lộ chi tiết của từng vỏ sao mà trước đây chúng ta không thể quan sát.5 là hai bức ảnh chụp bởi HET của hai ngôi sao ở các giai đoạn phát triển khác nhau.

Ảnh trái là ngôi sao trẻ V1331 Cyg hướng một cực về phía Trái Dất có vỏ sao là một đĩa bụi dày. Ảnh phải là tinh vân NGC 6826 có hình dạng như một con mắt với con ngươi là ngôi sao trung tâm.5: Cấu trúc vỏ sao chụp bởi HST: hình trái là ngôi sao trẻ V1331 Cyg (Nguồn: en. org/wiki/V1331_Cygni); hình phải là tỉnh vân tiền hành tỉnh NGC 6826 ( Nguồn: esahubble.org/ images/opo97384/). Những bức ảnh vùng khả kiến chủ yếu cho chúng ta nhìn thấy cấu trúc tổng quan của vỏ sao mà nhiều khi bị che khuất bởi lớp bụi dày.

Để nghiên cứu động học và tính chất vật lý của vỏ sao, các nhà thiên văn thường sử dụng bức xạ nguyên tử và phân tử khí tồn tại trong môi trường này. Những bức xạ nguyên tử và phân tử trong miền vô tuyến có thể đi qua lớp bụi dày của vỏ sao đồng thời đâm xuyên qua tầng khí quyền của Trái Đất.