Kiến trúc hệ thống tính toán thiên văn ephemeris, mô hình quy chiếu tọa độ hành tinh và chuẩn kiểm định độ chính xác Level-3 Quizzman

Tóm tắt (Abstract)

Bối cảnh (Background / Context): Các hệ thống tính toán tọa độ thiên văn, niên giám hành tinh và lập lịch hiện đại đòi hỏi sự kết hợp chặt chẽ giữa các mô hình giải tích cơ học thiên thể cổ điển và các nghiệm số tích phân quỹ đạo độ chính xác cao. Trong các ứng dụng dịch vụ tính toán thời gian thực, sự phụ thuộc đơn nhất vào các tệp dữ liệu kernel số trị kích thước lớn (như JPL DE series) thường gây suy giảm hiệu năng khởi động và hạn chế khả năng ngoại suy ngoài khoảng thời gian bao phủ, trong khi các công thức rút gọn cổ điển (như Jean Meeus rút gọn) lại phát sinh sai số vị trí đáng kể tại các điểm phân chí và ranh giới tuần trăng.

Phương pháp (Methodology): Công trình này giới thiệu kiến trúc phân tầng của hệ thống tính toán thiên văn @quizzman/qm-astro phiên bản 0.1.0, kết hợp giải pháp kép giữa nhân số trị JPL DE440 SPK (định dạng DAF/Chebyshev) và chuỗi giải tích hành tinh VSOP87D (bản đầy đủ của Bretagnon và Francou) cùng lý thuyết chuyển động Mặt Trăng ELP/MPP02 (Chapront và Francou). Hệ thống triển khai cơ chế điều phối chuyển mạch tự động CompositeProvider, áp dụng chuẩn quy đổi thời gian động học Trái Đất (TT/TDB) liên tục qua mô hình $\Delta T$ tích hợp giữa đa thức Espenak-Meeus và bảng tham số định hướng Trái Đất IERS finals2000A. Vị trí biểu kiến được quy chuẩn qua mô hình chương động IAU 2000B, quang hành sai hàng năm, hiệu chỉnh thời gian truyền ánh sáng, mô hình khúc xạ khí quyển Bennett / Kasten–Young có xét hạ độ cao chân trời địa hình, và thuật toán phân định ranh giới chòm sao IAU 1930 qua hình học đa giác phẳng.

Kết quả (Results): Báo cáo kiểm định thực nghiệm và bộ kiểm chuẩn độc lập Level-3 ghi nhận 100% tỷ lệ vượt qua (237/237 ca kiểm thử tự động) trên 6 mô-đun kiểm chứng. Độ sai biệt vectơ vị trí Cartesian so với hệ thống JPL Horizons đạt dưới 1 km đối với các hành tinh vòng trong. Kinh độ hoàng đạo biểu kiến của Mặt Trời tại điểm phân chí năm 2024 đạt sai số tồn dư $|\Delta\lambda_\odot| \ll 0.02^\circ$, thời điểm Sóc thiên văn khớp với dữ liệu quan sát của NASA trong khoảng 1.0 phút, và sai số thời điểm mọc/lặn của Mặt Trời tại kinh vĩ độ Hà Nội đạt dưới 30 giây khi tính đủ khúc xạ khí quyển.

Kết luận (Conclusion): Hệ thống @quizzman/qm-astro thiết lập một nền tảng tính toán ephemeris có cấu trúc phân tầng khả kiểm, đạt chuẩn độ chính xác Level-3, đóng vai trò nguồn chân lý thực thi (executable source of truth) phục vụ công cụ hiển thị bầu trời Sky Viewer, các dịch vụ tra cứu lịch thiên văn, và hạ tầng API vị trí hành tinh.

1 Phạm vi áp dụng (Scope)

Tài liệu này đặc tả kiến trúc kỹ thuật, mô hình toán học và các quy tắc chuẩn hóa của hệ thống tính toán thiên văn ephemeris Quizzman (@quizzman/qm-astro), bao gồm:

1. Quy chuẩn thang đo thời gian thiên văn và liên kết liên tục giữa Giờ Phối hợp Quốc tế (UTC), Giờ Quốc tế thực nghiệm (UT1), Giờ Nguyên tử Quốc tế (TAI), Giờ Động học Trái Đất (TT), và Giờ Động lực học Tâm tỉ cự (TDB).
2. Chuỗi tính toán giải tích vị trí Mặt Trời và các hành tinh qua lý thuyết VSOP87D, cùng vị trí địa tâm của Mặt Trăng qua lý thuyết ELP/MPP02.
3. Cơ chế đọc và nội suy đa thức Chebyshev từ tệp kernel số trị JPL DE440 SPK (định dạng DAF nhị phân).
4. Kiến trúc điều phối bộ cung cấp dữ liệu kép (CompositeProvider) với cơ chế fallback tự động bảo toàn tính liên tục ngoài khoảng thời gian bao phủ của kernel.
5. Quy trình quy chiếu tọa độ từ hệ quy chiếu thiên thể quốc tế (ICRS / J2000.0) sang vị trí biểu kiến trên kinh vĩ tuyến chân trời địa phương, bao gồm tiến trình tuế sai, chương động IAU 2000B, quang hành sai, khúc xạ khí quyển và độ võng chân trời do độ cao quan sát.
6. Phương pháp số hóa ranh giới 88 chòm sao theo định nghĩa IAU 1930 và lan truyền chuyển động riêng sao từ danh mục Hipparcos.
7. Giao diện tích hợp giữa lõi thiên văn học và hệ thống lịch pháp: thời điểm Sóc (New Moon), 24 Tiết Khí, và Giờ Thái Dương Thực.
8. Hồ sơ đo kiểm và tiêu chuẩn đánh giá độ chính xác phân cấp Level-3.

Tài liệu áp dụng cho toàn bộ các thành phần dịch vụ tính toán vị trí thiên thể, công cụ đồ họa Sky Viewer, hệ thống tính lịch âm dương, và các điểm cuối API liên quan thuộc hạ tầng Quizzman.

2 Tài liệu viện dẫn (Normative references)

Các tài liệu viện dẫn sau đây chứa các quy định cấu thành một phần không thể tách rời của đặc tả kỹ thuật này:

• ISO 8601-1:2019 ^[1], Date and time — Representations for information interchange — Part 1: Basic rules.
• ISO 2145:1978 ^[2], Documentation — Numbering of divisions and subdivisions in written documents.
• ISO 690:2021 ^[3], Information and documentation — Guidelines for bibliographic references and citations to information resources.
• R. S. Park et al. (2021) ^[4], The JPL Planetary and Lunar Ephemerides DE440 and DE441, The Astronomical Journal, Vol. 161, No. 3.
• P. Bretagnon, G. Francou (1988) ^[5], Planetary theories in rectangular and spherical variables. VSOP87 solutions, Astronomy and Astrophysics, Vol. 202, pp. 309–315.
• J. Chapront, G. Francou (2003) ^[6], The lunar theory ELP/MPP02, Astronomy and Astrophysics, Vol. 404, pp. 735–742.
• Jean Meeus (1998) ^[7], Astronomical Algorithms, 2nd Edition, Willmann-Bell, Richmond, Virginia.
• E. M. Standish (1998) ^[8], JPL Planetary and Lunar Ephemerides, DE405/LE405, JPL Interoffice Memorandum 312.F-98-048.
• D. D. McCarthy (2000) ^[9], IERS Conventions (2000): Chapter 5, IAU 2000B Nutation Model, IERS Technical Note 32.
• G. Petit, B. Luzum (2010) ^[10], IERS Conventions (2010), IERS Technical Note 36.
• F. Espenak, J. Meeus (2006) ^[11], Five Millennium Catalog of Solar Eclipses: -1999 to +3000, NASA TP-2006-214141.
• E. Delporte (1930) ^[12], Délimitation scientifique des constellations (tables et cartes), Cambridge University Press, IAU Commission 3.
• M. A. C. Perryman et al. (1997) ^[13], The HIPPARCOS Catalogue, Astronomy and Astrophysics, Vol. 323, pp. L49–L52.
• G. G. Bennett (1982) ^[14], The Calculation of Astronomical Refraction in Marine Navigation, Journal of Navigation, Vol. 35, No. 2, pp. 255–259.
• F. Kasten, A. T. Young (1989) ^[15], Revised optical air mass tables and approximation formula, Applied Optics, Vol. 28, No. 22, pp. 4735–4738.
• Quizzman Astronomy Working Group (2026) ^[16], Quizzman Astro Technical Reference, Level-3 Parity Audit and Verification Report, Technical Dossier TR-ASTRO-2026-07.

3 Thuật ngữ, định nghĩa và ký hiệu (Terms, definitions, and symbols)

Trong tài liệu này, các thuật ngữ và ký hiệu sau đây được áp dụng thống nhất:

3.1 Thang đo thời gian thiên văn (Astronomical time scales)

• UTC (Coordinated Universal Time): Giờ Phối hợp Quốc tế, thang đo thời gian dân sự dựa trên tần số nguyên tử SI nhưng được đồng bộ theo gia số giây nhuận (leap seconds) để duy trì sai lệch $|\text{UTC} - \text{UT1}| < 0.9\text{ s}$.
• UT1 (Universal Time): Giờ Quốc tế góc quay Trái Đất, tỷ lệ trực tiếp với góc quay thực tế của Trái Đất so với các nguồn vô tuyến chuẩn thiên văn ngoài thiên hà.
• TAI (International Atomic Time): Giờ Nguyên tử Quốc tế, thang đo thời gian nguyên tử liên tục chuẩn không chèn giây nhuận, liên hệ với UTC tại thời điểm hiện hành qua công thức: $$\text{TAI} = \text{UTC} + \Delta\text{AT}$$
• TT (Terrestrial Time): Giờ Động học Trái Đất, thang đo thời gian lý thuyết chuẩn dùng trong biểu diễn chuyển động của các thiên thể nhìn từ tâm Trái Đất, xác định qua: $$\text{TT} = \text{TAI} + 32.184\text{ s}$$
• TDB (Barycentric Dynamical Time): Giờ Động lực học Tâm tỉ cự, thang đo thời gian tương đối tính được định nghĩa tại tâm tỉ cự hệ Mặt Trời, chỉ khác biệt với TT bởi các số hạng tuần hoàn tuần suất rất nhỏ ($< 0.0017\text{ s}$).
• $\Delta T$ (Delta T): Hiệu số thời gian giữa Giờ Động học và Giờ Quốc tế, phản ánh sự suy giảm vận tốc tự quay của Trái Đất: $$\Delta T = \text{TT} - \text{UT1}$$

3.2 Điểm Sóc và kinh độ hoàng đạo (New Moon and ecliptic longitude)

• Kinh độ hoàng đạo ($\lambda$): Tọa độ góc của một thiên thể dọc theo mặt phẳng hoàng đạo tính từ điểm Xuân Phân ($0^\circ$).
• Thời điểm Sóc thiên văn (Astronomical New Moon): Thời điểm mà kinh độ hoàng đạo biểu kiến của Mặt Trăng trùng khớp chính xác với kinh độ hoàng đạo biểu kiến của Mặt Trời nhìn từ tâm Trái Đất ($\lambda_{\text{moon}} = \lambda_\odot$).

3.3 Hệ quy chiếu tọa độ thiên thể (Celestial reference systems)

• ICRS (International Celestial Reference System): Hệ quy chiếu thiên thể quốc tế chuẩn, không quay so với các thiên hà xa, với gốc tọa độ đặt tại tâm tỉ cự hệ Mặt Trời. Trục cơ bản xấp xỉ hướng điểm Xuân Phân trung bình tại kỷ nguyên chuẩn J2000.0 ($JD = 2451545.0$).
• Vị trí biểu kiến (Apparent place): Vị trí của thiên thể trên thiên cầu nhìn từ vị trí người quan sát sau khi đã áp dụng đầy đủ các hiệu chỉnh thời gian truyền ánh sáng, quang hành sai, chương động và khúc xạ khí quyển.

3.4 Khung tọa độ khám phá (Explorer coordinate frame)

• Khung tọa độ khám phá (Explorer frame): Hệ tọa độ Cartesian hoàng đạo có gốc tọa độ đặt tại trọng tâm của một thiên thể được chỉ định tùy ý ($center \neq \text{Earth}$), với đơn vị khoảng cách đo bằng đơn vị thiên văn (AU), phục vụ mô phỏng không gian đa vật thể trong công cụ Sky Viewer.

3.5 Cấp độ kiểm chuẩn chính xác Level-3 (Level-3 accuracy classification)

• Level 1 (Rút gọn): Áp dụng công thức giải tích xấp xỉ tối giản (Kepler 3 số hạng hoặc Jean Meeus rút gọn), phù hợp cho tính toán dân dụng thông thường với sai số vài phút cung.
• Level 2 (Giải tích đầy đủ): Sử dụng lý thuyết giải tích hành tinh VSOP87D và Mặt Trăng ELP/MPP02 đầy đủ số hạng, chương động IAU, và $\Delta T$ nội suy bảng.
• Level 2+ (Kết hợp số trị): Tích hợp hạt nhân số trị JPL DE440 SPK vào đường ống tính toán trực tiếp bên cạnh tầng giải tích.
• Level 3 (Chuẩn kiểm định khoa học): Đạt sự tương thích sai số định lượng với JPL Horizons ($< 1\text{ km}$ vị trí không gian), ranh giới 88 chòm sao IAU 1930 dạng đa giác, thời điểm mọc/lặn dưới 30 giây có khúc xạ khí quyển, và vượt qua 100% các bài kiểm thử tự động trong quy trình tích hợp liên tục (CI).

4 Quy ước và giả định tính toán (Conventions and assumptions)

Hệ thống @quizzman/qm-astro tuân thủ các quy ước và giả định định lượng sau:

4.1 Hằng số thiên văn cơ bản (Fundamental astronomical constants)

Các hằng số cơ bản được thiết lập theo khuyến nghị của IAU 2012 và JPL DE440:

• Vận tốc ánh sáng trong chân không: $c = 299,792,458\text{ m/s}$.
• Đơn vị thiên văn: $1\text{ AU} = 149,597,870,700\text{ m}$.
• Kỷ nguyên chuẩn J2000.0: $JD = 2451545.0\text{ TT}$ (tương ứng ngày 01 tháng 01 năm 2000 lúc 12:00:00 TT).
• Độ nghiêng hoàng đạo trung bình tại J2000.0: $\epsilon_0 = 23^\circ 26' 21.406'' = 84,381.406''$.

4.2 Quy ước kinh tuyến chuẩn và thời gian dân sự (Standard meridian and civil time)

• Kinh tuyến chuẩn của Việt Nam (UTC+07:00): $\lambda_{\text{std}} = 105^\circ\text{E}$.
• Kinh tuyến chuẩn của Trung Quốc (UTC+08:00): $\lambda_{\text{std}} = 120^\circ\text{E}$.
• Ranh giới ngày dân sự được cố định tuyệt đối tại thời điểm 00:00:00 giờ chuẩn địa phương.

4.3 Giả định chiết quang khí quyển và hạ độ cao đường chân trời (Atmospheric refraction and horizon dip)

Trong các tính toán vị trí biểu kiến nằm ngang và thời điểm mặt trời mọc/lặn:

1. Áp suất khí quyển tiêu chuẩn tại mặt biển: $P_0 = 1010\text{ hPa}$.
2. Nhiệt độ khí quyển tiêu chuẩn: $T_0 = 10^\circ\text{C}$ ($283.15\text{ K}$).
3. Bán kính góc biểu kiến của Mặt Trời tại khoảng cách trung bình: $r_\odot \approx 16'$.
4. Chiết quang khí quyển tại chân trời tiêu chuẩn: $R_0 = 34'$.
5. Góc hạ độ cao đường chân trời (horizon dip angle) phụ thuộc độ cao quan sát $h$ (tính bằng mét): $$d = 0.0293^\circ \times \sqrt{h}$$
6. Góc thiên đỉnh biểu kiến của tâm đĩa Mặt Trời tại thời điểm mọc/lặn thực tế: $$z_0 = 90^\circ + r_\odot + R_0 + d = 90^\circ 50' + 0.0293^\circ \times \sqrt{h}$$

4.4 Phân định ranh giới kiến trúc với Engine Lịch pháp (Architectural boundary with Calendar Engine)

@quizzman/qm-astro là hạt nhân tính toán cơ học thiên thể thuần túy, chuyên trách:

1. Tính toán tọa độ và vận tốc không gian của các thiên thể trong hệ tọa độ quán tính ICRS/J2000;
2. Xác định các thời điểm giao hội Sóc thiên văn ($\lambda_{\text{Moon}} = \lambda_\odot$) và 24 Tiết khí ($\lambda_\odot = k \times 15^\circ$) theo các thang đo thời gian vật lý liên tục (UT1, TT, TDB);
3. Quy đổi vị trí biểu kiến địa tâm và chân trời quan sát (có khúc xạ và hạ độ cao đường chân trời).

Engine thiên văn này hoàn toàn độc lập và không trực tiếp can thiệp vào các quy tắc lịch pháp dân sự mang tính pháp lý, lịch sử hoặc văn hóa xã hội (như lịch sử thay đổi múi giờ hành chính 1975, quy tắc tính tháng nhuận âm lịch, phân định ranh giới ngày dân sự nửa đêm, và hệ Can Chi). Toàn bộ các quy tắc lịch pháp đó được quy định độc lập và xử lý bởi engine @quizzman/qm-calendar (đặc tả tại [QP-CAL-001]).

5 Kiến trúc phân tầng của hệ thống ephemeris (Layered ephemeris system architecture)

Kiến trúc phần mềm của @quizzman/qm-astro được tổ chức thành 5 tầng độc lập, bảo đảm nguyên tắc phân tách trách nhiệm và kiểm thử cách ly:

┌────────────────────────────────────────────────────────┐
│ 5. Giao diện Ứng dụng & Dân sự (Civil & Almanac API)   │
│    Sóc, Tiết Khí, Giờ Thái Dương Thực, Sky Viewer JSON │
├────────────────────────────────────────────────────────┤
│ 4. Tầng Quy chiếu Biểu kiến (Apparent Reductions)      │
│    IAU 2000B, Quang hành sai, Khúc xạ khí quyển, IAU88 │
├────────────────────────────────────────────────────────┤
│ 3. Tầng Cung cấp Tọa độ Hợp nhất (CompositeProvider)   │
│    ├── Nhân số trị DE440 SPK (Khoảng: 1849 - 2150)    │
│    └── Nhân giải tích VSOP87D + ELP/MPP02 (Toàn vẹn)   │
├────────────────────────────────────────────────────────┤
│ 2. Tầng Chuyển động Thiên thể Nền tảng                 │
│    VSOP87D (Mặt Trời/Hành tinh), ELP/MPP02 (Mặt Trăng) │
├────────────────────────────────────────────────────────┤
│ 1. Tầng Cơ sở Thời gian & Hệ số Động học               │
│    JDN, UT ↔ TT, Espenak-Meeus ΔT, IERS finals2000A    │
└────────────────────────────────────────────────────────┘

5.1 Tầng thang đo thời gian và chuyển đổi động học (Time scale and dynamical time reduction)

Mọi phép tính toán vị trí bắt đầu bằng việc chuyển đổi thời điểm dân sự sang Julian Date (JD). Hiệu số $\Delta T$ được xác định theo quy tắc kết hợp:

1. Đối với các thời điểm thuộc giai đoạn hiện đại có dữ liệu quan trắc IERS: tra cứu trực tiếp giá trị thực nghiệm $\Delta T$ từ bảng finals2000A.
2. Đối với các thời điểm lịch sử và tương lai xa ngoài bảng: sử dụng đa thức xấp xỉ liên tục từng đoạn của Espenak và Meeus (2006) ^[11]. Ví dụ đối với giai đoạn 2005–2050: $$\Delta T = 62.92 + 0.32217,t + 0.005589,t^2$$ trong đó $t = y - 2000$.

5.2 Tầng lý thuyết giải tích hành tinh VSOP87D (Analytical planetary theory VSOP87D)

Lý thuyết VSOP87 phiên bản D cung cấp các tọa độ cầu nhật tâm của Trái Đất (kinh độ $L$, vĩ độ $B$, khoảng cách $R$) dưới dạng chuỗi lũy thừa thời gian $\tau$ (tính bằng thiên niên kỷ Julian từ J2000.0): $$L(\tau) = \sum_{k=0}^{5} \tau^k \sum_{i=1}^{N_k} A_i \cos(B_i + C_i \cdot \tau)$$ Hệ thống tải toàn bộ các hệ số của tệp VSOP87D.ear, bao gồm đầy đủ 6 bậc lũy thừa ($L_0$ đến $L_5$) với tổng cộng hàng trăm số hạng tuần hoàn, thay vì chỉ sử dụng 129 số hạng rút gọn.

Kinh độ nhật tâm của Trái Đất được chuyển dịch sang kinh độ địa tâm của Mặt Trời qua phép biến đổi: $$\lambda_\odot = L_{\text{helio}} + 180^\circ$$

5.3 Tầng lý thuyết chuyển động Mặt Trăng ELP/MPP02 (Analytical lunar theory ELP/MPP02)

Để xác định vị trí Mặt Trăng khi không có kernel số trị, hệ thống tích hợp lý thuyết bán giải tích ELP/MPP02 (phiên bản khớp với DE405) của Chapront và Francou ^[6]. Tọa độ địa tâm của Mặt Trăng (kinh độ $\lambda_{\text{moon}}$, vĩ độ $\beta_{\text{moon}}$, khoảng cách $\Delta_{\text{moon}}$) được tổng hợp từ các biểu thức tuần hoàn phụ thuộc vào 4 góc Delaunay cơ bản ($D, l, l', F$) và các biến nhiễu loạn hành tinh: $$\lambda_{\text{moon}} = W_1 + \sum A_j \sin(\sum k_i \phi_i)$$

5.4 Tầng số trị chính xác cao JPL DE440 SPK (High-precision numerical ephemeris DE440 SPK)

Tầng số trị đọc trực tiếp tệp hạt nhân nhị phân de440s.bsp theo chuẩn Double Precision Array File (DAF) của NASA SPICE toolkit. Mỗi đoạn dữ liệu (segment) kiểu Type 2 chứa các hệ số đa thức Chebyshev biểu diễn tọa độ Cartesian $(X, Y, Z)$ và vận tốc $(\dot{X}, \dot{Y}, \dot{Z})$ của một thiên thể trong một khoảng thời gian con $[t_a, t_b]$.

Thuật toán nội suy Chebyshev bậc $N$ đánh giá tọa độ qua công thức hồi quy Clenshaw: $$P(x) = \sum_{j=0}^{N} c_j T_j(x), \quad x = \frac{2t - (t_a + t_b)}{t_b - t_a} \in [-1, 1]$$ trong đó $T_0(x) = 1, T_1(x) = x, T_{j+1}(x) = 2x T_j(x) - T_{j-1}(x)$.

5.5 Cơ chế CompositeProvider và chiến lược fallback an toàn (CompositeProvider and fail-safe fallback)

Hệ thống thiết lập lớp bao bọc hợp nhất CompositeProvider (Quy tắc định chuẩn ASTRO-R001):

• Nếu tệp kernel DE440 có mặt và thời điểm tính toán nằm trong khoảng bao phủ $[JD_{\min}, JD_{\max}]$ (khoảng năm 1849 đến 2150 đối với de440s): truy vấn vị trí từ nhân số trị DE440.
• Nếu kernel không tồn tại trên máy chủ hoặc thời điểm tính toán nằm ngoài khoảng hiệu lực: chuyển mạch tự động, êm thuận sang nhân giải tích VSOP87D + ELP/MPP02.
• Nhãn nguồn ephemeris (de440 hoặc analytical) được trả về công khai trong payload phản hồi để bảo đảm tính minh bạch khoa học.

6 Thuật toán quy chiếu vị trí biểu kiến và chuyển đổi hệ tọa độ (Apparent position reduction and coordinate transformations)

6.1 Quang hành sai và hiệu chỉnh thời gian truyền ánh sáng (Stellar aberration and light-time correction)

Đối với các thiên thể trong hệ Mặt Trời, vị trí của thiên thể tại thời điểm nhận tín hiệu $t$ được xác định tại thời điểm bức xạ ánh sáng $(t - \tau_{\text{light}})$, trong đó: $$\tau_{\text{light}} = \frac{|\mathbf{r}$${\text{target}}(t - \tau{\text{light}}) - \mathbf{r}_{\text{observer}}(t)|}{c} Nghiệm được giải lặp từ 2 đến 3 bước hội tụ.

Hiệu chỉnh quang hành sai thường niên đối với Mặt Trời được áp dụng theo công thức thực nghiệm Ron-Vondrák: $$\Delta\lambda_{\text{aber}} = -20.49552'' \times \frac{1\text{ AU}}{R_\odot} \approx -0.00569^\circ$$

6.2 Chương động thiên văn mô hình IAU 2000B (Astronomical nutation IAU 2000B model)

Mô hình chương động IAU 2000B (phiên bản rút gọn 77 số hạng của IERS, kế thừa từ ERFA) được áp dụng để xác định độ lệch kinh độ chương động $\Delta\psi$ và độ lệch góc nghiêng hoàng đạo $\Delta\epsilon$. Thành phần chính phụ thuộc vào kinh độ nút thăng của quỹ đạo Mặt Trăng $\Omega$: $$\Omega = 125.04452^\circ - 1934.136261^\circ,\tau$$ Kinh độ hoàng đạo biểu kiến của Mặt Trời sau khi áp dụng quang hành sai và chương động: $$\lambda_{\odot,\text{app}} = \lambda_\odot - 0.00569^\circ - 0.00478^\circ \sin(\Omega) + \Delta\psi \cos(\epsilon)$$

6.3 Chuyển đổi giữa hệ xích đạo, hoàng đạo, chân trời và ngân hà (Equatorial, ecliptic, horizontal, and galactic transforms)

Hệ thống cung cấp ma trận quay trực giao bảo đảm bảo toàn độ dài vectơ:

1. Hoàng đạo sang Xích đạo: $$\begin{pmatrix} X \ Y \ Z \end{pmatrix}$${\text{equ}} = \begin{pmatrix} 1 & 0 & 0 \ 0 & \cos\epsilon & -\sin\epsilon \ 0 & \sin\epsilon & \cos\epsilon \end{pmatrix} \begin{pmatrix} X \ Y \ Z \end{pmatrix}{\text{ecl}}
2. Xích đạo sang Chân trời (Góc phương vị $A$, Độ cao $h$): $$\sin(h) = \sin(\phi)\sin(\delta) + \cos(\phi)\cos(\delta)\cos(H)$$ $$\tan(A) = \frac{\sin(H)}{\cos(H)\sin(\phi) - \tan(\delta)\cos(\phi)}$$ trong đó $\phi$ là vĩ độ quan sát, $\delta$ là xích vĩ, và $H$ là góc giờ địa phương ($H = \text{LST} - \alpha$).

6.4 Lan truyền chuyển động riêng sao danh mục Hipparcos (Hipparcos proper-motion propagation)

Đối với các định tinh thuộc tập con danh mục Hipparcos, vị trí thiên cầu tại thời điểm $T$ được cập nhật từ tọa độ tại kỷ nguyên chuẩn J1991.25 qua các thành phần chuyển động riêng $\mu_\alpha \cos\delta$ và $\mu_\delta$: $$\alpha(T) = \alpha_0 + \frac{\mu_\alpha \cos\delta}{\cos\delta_0},(T - T_0), \quad \delta(T) = \delta_0 + \mu_\delta,(T - T_0)$$

6.5 Xác định ranh giới chòm sao IAU 1930 qua thuật toán đa giác (IAU 1930 constellation boundary polygon algorithm)

Để xác định một thiên thể thuộc chòm sao nào trong số 88 chòm sao thiên văn, hệ thống số hóa toàn bộ các đoạn cung biên giới do Eugène Delporte thiết lập năm 1930 ^[12]. Tọa độ thiên thể tại kỷ nguyên khảo sát được tiến động quy hồi về kỷ nguyên gốc B1875.0, sau đó thuật toán kiểm tra điểm nằm trong đa giác cầu (Spherical Point-in-Polygon, Quy tắc định chuẩn ASTRO-R006) xác định mã chòm sao chính thức (ví dụ: UMi, Ori, Sgr).

7 Giao diện liên kết lịch pháp và thuật số dân sự (Civil calendrics and almanac interface)

7.1 Giải thuật xác định điểm Sóc và liên kết lịch âm (New Moon solving and lunisolar integration)

Thời điểm Sóc thiên văn được xác định bằng phương pháp tìm nghiệm phi tuyến của phương trình hiệu kinh độ hoàng đạo: $$f(JDE) = \lambda_{\text{moon}}(JDE) - \lambda_\odot(JDE) = 0 \pmod{360^\circ}$$ Hệ thống sử dụng nghiệm xấp xỉ sơ bộ theo Jean Meeus (Chương 49) ^[7] làm điểm bắt đầu, sau đó áp dụng phép lặp Newton-Raphson: $$JDE_{n+1} = JDE_n - \frac{\lambda_{\text{moon}}(JDE_n) - \lambda_\odot(JDE_n)}{\dot{\lambda}$${\text{moon}} - \dot{\lambda}\odot} Nghiệm hội tụ khi sai số thời gian giữa hai bước lặp liên tiếp đạt $|\Delta t| < 0.1\text{ s}$.

7.2 Giải thuật xác định 24 Tiết Khí và điểm Cardinal (Solar term and cardinal point determination)

Mỗi tiết khí thứ $m$ ($m \in [0, 23]$) ứng với kinh độ hoàng đạo biểu kiến của Mặt Trời $\lambda_{\odot} = m \times 15^\circ$. Thời điểm chuyển dịch tiết khí được xác định chính xác theo thời gian thực UTC, độc lập hoàn toàn với việc làm tròn ngày dân sự.

Bảng 1 — 4 Điểm phân chí cơ bản (Cardinal Points)

7.3 Chuẩn hóa Giờ Thái Dương Thực và phương trình thời gian (True Solar Time and Equation of Time)

Giờ Thái Dương Thực (True Solar Time - TST) phản ánh trực tiếp vị trí đỉnh nhật của Mặt Trời qua kinh tuyến địa phương của người quan sát: $$\text{TST} = \text{CivilTime} + \Delta t_{\text{lon}} + \text{EOT}$$ trong đó:

• $\Delta t_{\text{lon}} = (\lambda_{\text{local}} - \lambda_{\text{std}}) \times 4\text{ min/deg}$ (tính bằng phút).
• $\text{EOT}$ là phương trình thời gian (Equation of Time), biểu thị hiệu số giữa Mặt Trời thực và Mặt Trời trung bình: $$\text{EOT} = 4 \times (\text{GHA}_\odot - \text{GMST} - 180^\circ) \text{ min}$$

8 Kết quả kiểm định thực nghiệm và đo kiểm Level-3 (Empirical validation and Level-3 benchmark results)

Quy trình kiểm định hệ thống @quizzman/qm-astro phiên bản 0.1.0 được thực thi trên môi trường Node.js v20 / Rust 1.80 trên hạ tầng máy chủ Linux x86_64, ghi nhận toàn bộ hồ sơ đo kiểm trong tệp báo cáo report.json công khai tại cổng thông tin https://astro.quizzman.com/paper/report.json.

8.1 Độ lệch vectơ so chuẩn JPL Horizons (JPL Horizons vector residual evaluation)

Khi tệp kernel số trị de440s.bsp được tải, các vectơ trạng thái ba chiều $(X, Y, Z)$ của Mặt Trời và các hành tinh vòng trong (Sao Thủy, Sao Kim, Trái Đất - Mặt Trăng barycenter, Sao Hỏa) được so chuẩn đối ứng với các vectơ phát hành bởi hệ thống JPL Horizons Web API tại kỷ nguyên J2000.0 và thời điểm ngẫu nhiên trong khoảng năm 1900–2100:

• Khoảng cách sai số vị trí không gian: $|\Delta\mathbf{r}| < 0.85\text{ km}$ (đáp ứng tiêu chuẩn Level-3 $< 1\text{ km}$).
• Sai số vận tốc quỹ đạo: $|\Delta\mathbf{v}| < 0.001\text{ m/s}$.

8.2 Sai số kinh độ hoàng đạo tại điểm phân chí 2024 (Ecliptic longitude residuals at 2024 equinoxes)

Bảng 2 — Độ sai lệch kinh độ hoàng đạo tại các điểm phân chí năm 2024

| Điểm phân chí | Thời điểm chuẩn UTC | Kinh độ lý thuyết | Kết quả @quizzman/qm-astro | Độ lệch tồn dư $|\Delta\lambda_\odot|$ | Đánh giá | |---|---|---|---|---|---| | Xuân Phân | 2024-03-20 03:06:21 | $0^\circ$ | $359.9982^\circ$ | $0.0018^\circ$ | Đạt ($< 0.02^\circ$) | | Hạ Chí | 2024-06-20 20:50:56 | $90^\circ$ | $89.9991^\circ$ | $0.0009^\circ$ | Đạt ($< 0.02^\circ$) | | Thu Phân | 2024-09-22 12:43:36 | $180^\circ$ | $179.9987^\circ$ | $0.0013^\circ$ | Đạt ($< 0.02^\circ$) | | Đông Chí | 2024-12-21 09:20:30 | $270^\circ$ | $269.9994^\circ$ | $0.0006^\circ$ | Đạt ($< 0.02^\circ$) |

Toàn bộ sai số tồn dư tại các điểm Cardinal đều nằm sâu dưới ngưỡng quy chuẩn $0.02^\circ$, xác nhận độ ổn định của chuỗi VSOP87D và bộ hiệu chỉnh chương động.

8.3 Đối chuẩn điểm Sóc thiên văn với dữ liệu NASA (New Moon validation against NASA catalog)

Thời điểm Sóc thiên văn ngày 11/01/2024 được so sánh đối chuẩn trực tiếp:

• Danh mục NASA Eclipse/Lunar Almanac ^[11]: $11:57\text{ UTC}$.
• Kết quả tính toán @quizzman/qm-astro: $11:58:02\text{ UTC}$.
• Độ lệch thời gian tuyệt đối: $\Delta t \approx 1.03\text{ phút}$.

Độ lệch này hoàn toàn nằm trong giới hạn biến động lý thuyết của các thuật toán ephemeris geocentric bậc cao khi không áp dụng vi hiệu chỉnh nhiễu loạn địa hình mặt trăng (topocentric lunar limb).

8.4 Đo kiểm thời điểm mặt trời mọc/lặn và khúc xạ khí quyển (Sunrise/sunset and atmospheric refraction verification)

Thời điểm mọc và lặn của Mặt Trời tại trạm quan trắc Hà Nội ($21.0285^\circ\text{N}, 105.8542^\circ\text{E}$, độ cao $12\text{ m}$) trong ngày Hạ chí 2024:

• Thời điểm mọc tính toán: 05:15:22 (Giờ địa phương UTC+07:00).
• Sai số so chuẩn mô hình tiêu chuẩn Hải quân Hoa Kỳ (USNO): $\Delta t_{\text{rise}} = 18\text{ giây}$.
• Sai số thời điểm lặn: $\Delta t_{\text{set}} = 24\text{ giây}$. Cả hai giá trị đều thỏa mãn tiêu chuẩn khắt khe Level-3 ($< 30\text{ giây}$).

8.5 Kiểm thử hồi quy 237 ca kiểm chuẩn tự động (Automated 237-case regression test matrix)

Bảng 3 — Tổng hợp kết quả đo kiểm Level-3 qua 6 mô-đun phần mềm

Tỷ lệ vượt qua đạt 100% (237/237 ca kiểm chuẩn), xác lập trạng thái sẵn sàng Level-3 của hệ thống.

9 Giới hạn kỹ thuật và phạm vi hiệu lực (Limitations)

Để duy trì tính nghiêm mật trong ứng dụng khoa học, các giới hạn kỹ thuật sau đây được công bố minh thị:

1. Phạm vi hiệu lực của kernel DE440 SPK: Tệp hạt nhân rút gọn de440s.bsp tích hợp trong hệ thống chỉ bao phủ khoảng thời gian từ năm 1849 đến năm 2150. Đối với các thời điểm ngoài khoảng này, hệ thống tự động suy thoái về chuỗi giải tích VSOP87D + ELP/MPP02. Độ chính xác ngoài khoảng này sẽ suy giảm dần theo thời gian thiên niên kỷ do sự tích lũy sai số thế kỷ của các số hạng tự do.
2. Danh mục sao: Hệ thống chỉ tích hợp tập con các định tinh sáng có chỉ số thị sai từ danh mục Hipparcos, không thay thế cho các danh mục hoàn chỉnh mật độ cao như Gaia DR3 (hơn 1.8 tỷ thiên thể).
3. Mô hình nhật thực toàn phần: Tính toán nhật thực hiện thời cung cấp pha cực đại địa tâm và tiếp xúc biểu kiến sơ bộ; thuật toán tính dải bóng tối đường đi chi tiết trên bề mặt Trái Đất (Besselian elements path track) hiện đang trong lộ trình hoàn thiện.
4. Tham số góc quay Trái Đất UT1-UTC: Dữ liệu dự báo EOP của IERS có độ tin cậy cao trong vòng 1 năm gần nhất; đối với các năm tương lai xa, hệ thống chuyển sang xấp xỉ thuần túy theo mô hình đa thức thế kỷ $\Delta T$.

10 Kết luận (Conclusion)

Hệ thống @quizzman/qm-astro thiết lập một kiến trúc tính toán ephemeris hai tầng (dual-backend) vững chắc, kết hợp thành công tính tiện dụng, liên tục của các lý thuyết giải tích hành tinh (VSOP87D, ELP/MPP02) với độ chính xác cao của nhân tích phân số trị hiện đại NASA JPL DE440. Toàn bộ chuỗi quy chiếu tọa độ biểu kiến được chuẩn hóa theo các nghị quyết của IAU và khuyến nghị của IERS, giải quyết triệt để sự thiếu đồng nhất giữa thời gian động học và thời gian dân sự.

Kết quả kiểm định khách quan với 237/237 ca kiểm chuẩn vượt qua đã chứng thực trạng thái độ chính xác Level-3 của hệ thống. Đây là nguồn chân lý thực thi chuẩn mực, bảo đảm tính tái lập khoa học và cung cấp nền tảng tính toán vững chắc cho hệ sinh thái thiên văn, niên giám và lịch pháp Quizzman.

Bibliography

• ^[1] International Organization for Standardization. (2019). ISO 8601-1:2019 Date and time — Representations for information interchange — Part 1: Basic rules. Geneva: ISO.
• ^[2] International Organization for Standardization. (1978). ISO 2145:1978 Documentation — Numbering of divisions and subdivisions in written documents. Geneva: ISO.
• ^[3] International Organization for Standardization. (2021). ISO 690:2021 Information and documentation — Guidelines for bibliographic references and citations to information resources. Geneva: ISO.
• ^[4] Park, R. S., Folkner, W. M., Williams, J. G., & Boggs, D. H. (2021). The JPL Planetary and Lunar Ephemerides DE440 and DE441. The Astronomical Journal, 161(3), 105.
• ^[5] Bretagnon, P., & Francou, G. (1988). Planetary theories in rectangular and spherical variables. VSOP87 solutions. Astronomy and Astrophysics, 202, 309–315.
• ^[6] Chapront, J., & Francou, G. (2003). The lunar theory ELP/MPP02. Astronomy and Astrophysics, 404, 735–742.
• ^[7] Meeus, J. (1998). Astronomical Algorithms (2nd ed.). Richmond, Virginia: Willmann-Bell.
• ^[8] Standish, E. M. (1998). JPL Planetary and Lunar Ephemerides, DE405/LE405 (JPL IOM 312.F-98-048). Pasadena: Jet Propulsion Laboratory.
• ^[9] McCarthy, D. D. (2000). IERS Conventions (2000): Chapter 5, IAU 2000B Nutation Model (IERS Technical Note 32). Frankfurt am Main: IERS.
• ^[10] Petit, G., & Luzum, B. (Eds.). (2010). IERS Conventions (2010) (IERS Technical Note 36). Frankfurt am Main: BKG.
• ^[11] Espenak, F., & Meeus, J. (2006). Five Millennium Catalog of Solar Eclipses: -1999 to +3000 (NASA TP-2006-214141). Greenbelt: NASA GSFC.
• ^[12] Delporte, E. (1930). Délimitation scientifique des constellations (tables et cartes). Cambridge: Cambridge University Press.
• ^[13] Perryman, M. A. C., Lindegren, L., Kovalevsky, J., & Hoeg, E. (1997). The HIPPARCOS Catalogue. Astronomy and Astrophysics, 323, L49–L52.
• ^[14] Bennett, G. G. (1982). The Calculation of Astronomical Refraction in Marine Navigation. Journal of Navigation, 35(2), 255–259.
• ^[15] Kasten, F., & Young, A. T. (1989). Revised optical air mass tables and approximation formula. Applied Optics, 28(22), 4735–4738.
• ^[16] Quizzman Astronomy Working Group. (2026). Quizzman Astro Technical Reference, Level-3 Parity Audit and Verification Report (Technical Report TR-ASTRO-2026-07). Hanoi: Quizzman Research Lab.

Annex A (normative) Danh mục quy tắc định chuẩn ASTRO-R (Normative Rule Catalog)

Phụ lục này quy định danh mục các quy tắc định chuẩn bắt buộc áp dụng trong phần mềm tính toán thiên văn của hệ sinh thái Quizzman:

A.1 ASTRO-R001: Lựa chọn bộ cung cấp ephemeris theo ranh giới thời gian

Phần mềm triển khai SHALL kiểm tra tính sẵn sàng của hạt nhân số trị DE440 SPK. Khi thời điểm tính toán $JD \in [2,396,400.5, 2,506,300.5]$ và tệp dữ liệu tồn tại, hệ thống SHALL sử dụng bộ nội suy đa thức Chebyshev từ kernel DE440. Khi nằm ngoài phạm vi này hoặc không tìm thấy kernel, hệ thống SHALL chuyển mạch tự động sang bộ giải tích VSOP87D và ELP/MPP02.

A.2 ASTRO-R002: Quy đổi thời gian động học TT

Phép quy đổi từ UTC sang TT SHALL tính đến số giây nhuận chính thức thông qua TAI đối với các niên đại từ năm 1972 đến nay. Đối với các thời điểm có ghi nhận trong bảng IERS, hệ thống SHALL áp dụng trực tiếp giá trị thực nghiệm $\Delta T$. Đối với các niên đại lịch sử hoặc tương lai ngoài bảng, hệ thống SHALL áp dụng hàm đa thức Espenak-Meeus (2006).

A.3 ASTRO-R003: Tính toán kinh độ hoàng đạo biểu kiến Mặt Trời

Kinh độ hoàng đạo biểu kiến của Mặt Trời SHALL được tính toán có áp dụng đầy đủ:

1. Quang hành sai thường niên $\Delta\lambda_{\text{aber}} = -20.49552'' / R$.
2. Hiệu chỉnh chương động kinh độ $\Delta\psi \cos(\epsilon)$ theo mô hình IAU 2000B hoặc ERFA.

A.4 ASTRO-R004: Giải nghiệm thời điểm Sóc thiên văn

Thời điểm Sóc thiên văn SHALL là nghiệm của phương trình $\lambda_{\text{moon}}(JDE) - \lambda_\odot(JDE) = 0 \pmod{360^\circ}$. Thuật toán lặp SHALL bảo đảm sai số hội tụ của thời điểm điểm Sóc nhỏ hơn hoặc bằng 0.1 giây.

A.5 ASTRO-R005: Quy đổi Giờ Thái Dương Thực

Giờ Thái Dương Thực SHALL là tổng của giờ đồng hồ dân sự, hiệu số kinh độ địa phương so với kinh tuyến chuẩn ($( \lambda_{\text{local}} - \lambda_{\text{std}} ) \times 4\text{ min/deg}$), và phương trình thời gian $\text{EOT}$ tại thời điểm khảo sát.

A.6 ASTRO-R006: Phân loại ranh giới chòm sao IAU 1930

Để phân loại một thiên thể vào 88 chòm sao thiên văn, tọa độ xích đạo của thiên thể SHALL được tuế sai ngược về kỷ nguyên chuẩn B1875.0 trước khi thực hiện thuật toán kiểm tra điểm trong đa giác đối với tập dữ liệu ranh giới Eugène Delporte (1930).

Annex B (informative) Bảng đối chuẩn phân chí và điểm Sóc thiên văn 2024–2026 (Astronomical Benchmark Tables 2024–2026)

Phụ lục này cung cấp dữ liệu đối chuẩn thực nghiệm để kiểm chứng chéo tính tương thích của các bản cài đặt phần mềm độc lập:

Bảng B.1 — Các mốc phân chí thiên văn giai đoạn 2024–2026 (Giờ UTC)

Bảng B.2 — Các điểm Sóc thiên văn năm 2024 (So chuẩn với NASA Catalog)