中子星是恒星演化末期形成的致密天体，是大质量恒星核心在引力坍缩和超新星爆发后的终极产物。中子星以极端物理条件著称：其密度超过原子核密度、引力场仅次于黑洞、磁场强度可达地球的万亿倍、自转速度可达每秒数百转。自1967年发现脉冲星（旋转的中子星）以来，中子星一直是天体物理学研究的前沿，为极端物质状态、引力理论和核物理提供了天然的实验室。
1932年，詹姆斯·查德威克（James Chadwick）发现中子。仅一年后，天文学家沃尔特·巴德（Walter Baade）和弗里茨·兹威基（Fritz Zwicky）在1933年美国物理学会会议上首次提出"中子星"的概念，认为超新星爆发可能将普通恒星压缩成由中子构成的天体。然而，由于当时缺乏观测手段，这一理论预言未能得到验证。此后数十年，中子星仅停留在理论物理中。
1939年，罗伯特·奥本海默（Robert Oppenheimer）和乔治·沃尔科夫（George Volkoff）建立第一个中子星结构模型——Tolman-Oppenheimer-Volkoff（TOV）方程，证明中子星的存在上限约为3个太阳质量。这一理论框架至今仍是中子星研究的基石。
| 年份 | 里程碑事件 |
|---|---|
| 1932 | 查德威克发现中子 |
| 1933 | 巴德和兹威基提出中子星概念 |
| 1939 | 奥本海默和沃尔科夫建立TOV方程 |
| 1967 | 贝尔和休伊什发现第一颗脉冲星（CP 1919） |
| 1968 | 确认脉冲星为旋转中子星 |
| 1974 | 赫尔斯-泰勒脉冲双星发现（间接验证引力波） |
| 1982 | 发现毫秒脉冲星（PSR B1937+21，自转1.56 ms） |
| 1998 | 发现磁星（SGR 1806-20） |
| 2004 | SGR 1806-20巨型耀斑（地球电离层受扰） |
| 2017 | LIGO/Virgo探测到双中子星并合引力波（GW170817） |
| 2019 | NICER精确测量脉冲星质量半径 |
| 2023 | NANOGrav探测到纳赫兹引力波背景（脉冲星计时阵列） |
1967年，剑桥大学的博士生乔斯林·贝尔（Jocelyn Bell）和导师安东尼·休伊什（Antony Hewish）在利用射电望远镜阵列观测行星际闪烁时，发现了一个周期为1.337秒的规律脉冲信号。信号极其精准，周期稳定到纳秒级，一度被称为"小绿人一号"（LGM-1）。后续观测排除了地外文明的解释，确认其为快速旋转的中子星——脉冲星。
这一发现震惊了天文学界。脉冲星如同宇宙中的灯塔，其射电波束随自转扫过地球，形成脉冲信号。贝尔的发现开创了脉冲星天文学，休伊什因"脉冲星发现"获1974年诺贝尔物理学奖（贝尔未获奖，成为科学史上著名的争议话题）。
并非所有恒星都会演化成中子星。演化路径取决于恒星初始质量：
| 初始质量范围 | 演化终点 | 遗留天体 |
|---|---|---|
| 褐矮星 | 无（不足以点燃氢聚变） | |
| 行星状星云 | 白矮星 | |
| 超新星（II型/Ib/Ic型） | 中子星 | |
| 超新星（有时含伽马暴） | 黑洞 | |
| 直接坍缩（无超新星） | 黑洞 |
其中 为太阳质量。注意质量边界因金属丰度不同而有几个太阳质量的浮动，以上为典型范围。
大质量恒星核心在耗尽核燃料后，铁核无法通过聚变释放能量，核心的辐射压不足以对抗引力，触发核心坍缩。过程可分为几个阶段：
这一过程释放大量中微子
数值案例：一颗 的恒星在超新星爆发后：
电子俘获过程中的关键反应：
中子星的核心参数范围之极端令人难以想象：
| 参数 | 典型值 | 与地球对比 | 与太阳对比 |
|---|---|---|---|
| 质量 | （约 ） | 倍 | 倍 |
| 半径 | |||
| 平均密度 | 倍 | 倍 | |
| 表面重力 | 倍 | 倍 | |
| 逃逸速度 | 倍 | 倍 | |
| 表面温度 | 倍 | 倍 | |
| 磁场强度 | 倍 | 倍 | |
| 自转周期 | 倍 | 更快 |
一茶匙（约 5 毫升）的中子星物质：
中子星的半径主要由其质量决定。不同物态方程（EOS）给出不同的质量-半径关系：
| 质量 | 半径（硬EOS，km） | 半径（中等EOS，km） | 半径（软EOS，km） |
|---|---|---|---|
| 0.5 | 13.5 | 12.8 | 11.8 |
| 1.0 | 13.8 | 13.0 | 11.9 |
| 1.4 | 13.2 | 12.4 | 11.2 |
| 1.8 | 12.6 | 11.7 | 10.3 |
| 2.0 | 12.0 | 11.0 | 9.5 |
| 2.2 | 11.2 | 10.0 | — |
| 2.3 | 10.6 | 9.3 | — |
注：硬EOS（如MS1）预测更大半径，软EOS（如APR）预测更小半径。GW170817和NICER观测倾向于支持中等EOS，典型半径约12 km。
中子星的内部结构从外到内可分为若干层，每层的物理性质截然不同：
┌──────────────────────────────┐
│ 大气层（~1 cm） │ ← 等离子体，温度 $10^5-10^7$ K
├──────────────────────────────┤
│ 外 壳（~1 km） │ ← 晶格结构，电子简并
│ - 外外层（~500 m） │ 原子核 + 相对论电子气
│ - 内外壳层（~500 m） │ 自由中子开始出现
├──────────────────────────────┤
│ 外 核（~9 km） │ ← 超流中子和质子
│ │ 密度 $2-8\times 10^{17}\ \mathrm{kg/m^3}$
│ - 均匀核物质区 │ 中子和质子比例约20:1
├──────────────────────────────┤
│ 内 核（~1-3 km） │ ← 未知物态
│ │ 密度 > $8\times 10^{17}\ \mathrm{kg/m^3}$
│ - 可能含π介子凝聚 │
│ - 可能含超子 │
│ - 可能含夸克物质 │
│ - 甚至可能是奇异夸克星 │
└──────────────────────────────┘
外层外壳（Outer Crust）：
内层外壳（Inner Crust）：
外核：
内核（最神秘的部分）：
物态方程在密度超过 时变得不确定，以下是几种可能的内核组成：
| 内核组成模型 | 关键特征 | 预测最大质量 | 观测检验可能性 |
|---|---|---|---|
| 核子物质（npeμ） | 中子和质子维持费米液体 | GW170817支持 | |
| 超子物质 | 出现Λ、Σ、Ξ超子 | 与PSR J0740+6620矛盾 | |
| 夸克物质（混合相） | 核心有夸克滴 | LIGO可检验 | |
| 奇异夸克物质（u、d、s） | 奇异夸克星，无核壳 | 可通过冷却观测区分 | |
| π介子凝聚 | π-π相互作用形成Bose-Einstein凝聚 | 冷却行为改变 | |
| K介子凝聚 | 奇异凝聚 | 与大多数观测矛盾 |
数值案例：典型中子星的结构参数
以一颗质量 、半径 的中子星为例：
| 层 | 半径范围 | 质量占比 | 密度范围 | 温度 |
|---|---|---|---|---|
| 大气层 | ||||
| 外壳 | ||||
| 外核 | ||||
| 内核 |
中子星诞生时温度极高（约 ），通过中微子辐射和光子辐射快速冷却：
其中 是中微子光度， 是表面光子光度， 是热容。
冷却曲线取决于核心是否存在快速冷却机制（如超流中子→质子转变的URCA过程）：
| 年龄 | 温度（标准冷却） | 温度（快速冷却） | 观测到的表面温度 |
|---|---|---|---|
| 10年 | — | ||
| 100年 | N/A | ||
| 1000年 | 蟹状星云脉冲星 | ||
| 10000年 | Vela脉冲星 | ||
| 100万年 | 典型中年脉冲星 | ||
| 10亿年 | 不可观测 |
通过X射线望远镜观测中子星的表面温度，可以约束其核心的物态和冷却机制。
中子星根据其观测特征可分为多种类型：
定义：具有强磁场（）且高速旋转的中子星，其磁极发出的辐射（射电、光学、X射线、伽马射线）随自转扫过地球，形成脉冲信号。
脉冲机制：
磁轴 磁轴
│ │
──━┿━━── 自转轴 ──━┿━━── 自转轴
│ ╲ │ ╲
│ ╲ 辐射束 │ ╲ 辐射束
│ ╲ │ ╲
││││ ╲ ││││ ╲
││││┌───┼─╲ ││││┌───┼─╲
│ ── │ ──
──┴── ──┴──
磁极 磁极
当磁轴与自转轴不平行时，射电波束如灯塔般扫过宇宙空间。若地球恰好在波束扫过的平面内，就能观测到脉冲信号。
周期-周期导数图（ 图）：
脉冲星的特征可以用自转周期 和周期变化率 来描述，这是脉冲星分类的基本工具：
log(˙P)
|
-10| ● ●●● ● = X射线脉冲星
| ● ●● ●●●● ▲ = 伽马射线脉冲星
-12| ●● ●● ○ = 正常射电脉冲星
| ● ●●●● ■ = 毫秒脉冲星
-14| ● ○○○
| ○○ ○○○○ ○
-16| ○○ ○○○○ ○○
| ○○ ○○ ○ ○ 恢复线（recycling line）
-18| ■■
| ■■■■
-20| ■■■■■■ ■■■
└─────┬──────┬──────┬──────┬──
-3 -1 1 3 log(P/s)
典型特征：
周期与年龄的关系：
脉冲星的年龄可近似估计为：
其中 称为特征年龄。
数值案例：
实际年龄（公元1054年超新星爆发）：约 970 年，两者相近
定义：磁场强度极高（）的中子星。磁星的能量来源不是自转能，而是磁能的释放。磁场能量密度可达：
这相当于核物质能量密度的千分之一，比太阳内部能量密度高10个数量级。
磁星类型：
| 类型 | 全称 | 观测特征 | 已知数量 |
|---|---|---|---|
| SGR | 软伽马射线重复源（Soft Gamma Repeater） | 间歇性软伽马暴、巨型耀斑 | ~7 |
| AXP | 反常X射线脉冲星（Anomalous X-ray Pulsar） | 稳定X射线辐射，周期2-12s | ~14 |
巨型耀斑：
2004年12月27日，SGR 1806-20 释放了有记录以来最强大的中子星巨耀斑：
| 参数 | 数值 | 对比 |
|---|---|---|
| 总能量 | 太阳30万年总辐射能 | |
| 峰值光度 | 银河系总光度的100倍 | |
| 持续时间 | 主脉冲 | — |
| 对地球影响 | 电离层扰动 | 卫星仪器饱和，堪比太阳耀斑 |
| 距离 | 银河系另一侧 |
许多中子星存在于双星系统中，通过与伴星的引力相互作用产生丰富现象：
| 系统类型 | 描述 | 质量来源 | 典型周期 | 示例 |
|---|---|---|---|---|
| X射线双星 | 吸积伴星物质产生X射线 | 恒星级伴星 | 数小时-数天 | 天鹅座X-1 |
| 毫秒脉冲星双星 | 从伴星吸积物质被"回收"再加速 | 白矮星或低质量恒星 | 数天 | PSR J0437-4715 |
| 双脉冲星 | 两颗中子星互相绕转 | N/A | 数小时 | PSR J0737-3039 |
| 脉冲星-黑洞 | 中子星绕黑洞运动 | N/A | — | 尚未确认 |
吸积过程的关键物理：
当物质从中子星伴星落入其中子星的引力势阱时，释放的引力能远超核聚变能：
对于每千克物质，释放的引力势能约为 ，而核聚变释放的约为 ——相差 250倍。
射电波段是脉冲星最常用的观测窗口。射电望远镜通过记录脉冲到达时间（TOA）来测量脉冲星的自转参数和位置：
主要射电望远镜：
| 望远镜 | 位置 | 有效口径 | 脉冲星发现数 |
|---|---|---|---|
| FAST（中国天眼） | 中国贵州 | 300 m | 800+ |
| Arecibo（已退役） | 波多黎各 | 305 m | 2000+ |
| Parkes | 澳大利亚 | 64 m | 1500+ |
| Green Bank | 美国西弗吉尼亚 | 100 m | 300+ |
| LOFAR | 荷兰（阵列） | 等效数千平米 | 200+ |
脉冲星计时：
通过精确测量脉冲的到达时间，可以实现极高的测量精度：
其中 是脉冲宽度，S/N 是信噪比。
对于蟹状星云脉冲星，经过20分钟积分，TOA精度可达 ，对应位置精度约 （在距离2200 pc处）！
中子星在不同波段辐射的物理机制不同：
| 波段 | 辐射机制 | 典型天体 | 关键望远镜 |
|---|---|---|---|
| 射电 | 相干曲率辐射 | 射电脉冲星 | FAST, Parkes, LOFAR |
| 光学 | 非热同步辐射（+热辐射） | 蟹状星云脉冲星 | HST, VLT |
| X射线 | 热表面辐射 + 非热磁层辐射 | 磁星、X射线双星 | Chandra, XMM-Newton, NICER |
| 伽马射线 | 磁层外间隙加速辐射 | 伽马脉冲星 | Fermi-LAT, MAGIC |
| 引力波 | 并合过程、非球对称形变 | 双中子星 | LIGO, Virgo, KAGRA |
NICER实验：
NASA的NICER（中子星内部组成探测器）于2017年安装在国际空间站，通过精确测量脉冲星的X射线轮廓，约束中子星的质量和半径。2019年，NICER测量了脉冲星PSR J0030+0451的质量和半径：
2021年，更高质量的PSR J0740+6620测量（已知最重中子星之一）：
这些测量对约束高密度物态方程至关重要。
GW170817事件（2017年8月17日）是首个被探测到的双中子星并合引力波信号，标志着多信使天文学新时代的开启：
| 参数 | 测量值 |
|---|---|
| 信号持续时间 | （在LIGO频带内） |
| 频率范围 | |
| 信噪比 | |
| 传播距离 | （约1.3亿光年） |
| 并合前总质量 | |
| 并合产物 | 可能是大质量中子星短暂存在后坍缩为黑洞 |
| 电磁对应体 | 短伽马暴GRB 170817A + 千新星AT 2017gfo |
| 引力波速度 | 等于光速，精度 |
GW170817提供了多个重要约束：
| 名称 | 质量 | 类型 | 距离 | 发现年份 | 特点 |
|---|---|---|---|---|---|
| PSR J0740+6620 | 毫秒脉冲星 | 约3800 ly | 2019 | 目前确认最重 | |
| PSR J2215+5135 | 黑洞掩星脉冲星 | 约10000 ly | 2018 | 争议较大 | |
| PSR J0348+0432 | 脉冲星-白矮星 | 约7000 ly | 2013 | 精确测量 | |
| PSR B1516+02B | 球状星团脉冲星 | 约21000 ly | 2023 | 新发现 | |
| Vela X-1 | 吸积X射线脉冲星 | 约6200 ly | 1970s | 高质量 |
| 名称 | 位置 | 周期 | 特点 | 发现年份 |
|---|---|---|---|---|
| 蟹状星云脉冲星 (PSR B0531+21) | 蟹状星云中心 | 33 ms | 最亮的脉冲星，可见光也可观测 | 1968 |
| Vela脉冲星 (PSR B0833-45) | 船帆座超新星遗迹 | 89 ms | 首个观测到星震的脉冲星 | 1968 |
| 赫尔斯-泰勒脉冲星 (PSR B1913+16) | 天鹰座 | 59 ms | 首个脉冲双星，引力波间接验证 | 1974 |
| 双脉冲星 (PSR J0737-3039) | 人马座 | 22.7 ms | 唯一已知的双脉冲星系统 | 2003 |
| 毫秒脉冲星 (PSR B1937+21) | 狐狸座 | 1.56 ms | 首个毫秒脉冲星 | 1982 |
| 第谷超新星遗迹 (PSR J0454+5543) | 仙后座 | — | 2018年发现 | 2018 |
| M15中子星 (PSR B2127+11C) | M15球状星团 | 30.5 ms | 球状星团中脉冲星 | 1990 |
作为研究最深入的中子星，蟹状星云脉冲星提供了丰富的信息：
中子星可以有多重？这是物态方程研究的核心问题之一。
已知约束：
理论极限：
由广义相对论的因果率和稳定条件，可推导出极限：
不同物态方程预测：
| EOS | 预测最大质量 | GW170817约束 | 是否被排除 |
|---|---|---|---|
| MS1 | 兼容 | 否 | |
| H4 | 兼容 | 否 | |
| APR | 兼容 | 否 | |
| ALF2 | 兼容 | 否 | |
| Sly | 兼容 | 否 | |
| GM1 | 兼容 | 否 | |
| BHF | 兼容 | 边缘 | |
| PAL6 | 不兼容 | 已排除 | |
| G240 | 不兼容 | 已排除 |
中子星内核可能存在退禁闭的夸克物质（u-d-s夸克的等离子体）。如果整个星体都由奇异夸克物质组成，则称为奇异夸克星。
夸克星与普通中子星的区别：
| 特征 | 普通中子星 | 奇异夸克星 |
|---|---|---|
| 主要成分 | 中子+质子+超子 | u,d,s夸克 |
| 表面 | 固态外壳（原子核晶格） | 夸克物质直接暴露（或极薄外壳） |
| 密度分布 | 中心高→表面低 | 更均匀 |
| 最大质量 | ||
| 半径-质量关系 | 质量↑半径↓ | 质量↑半径可能↑ |
| 冷却速率 | 中微子主导 | 夸克直URCA更快 |
| 观测线索 | 多种 | 快速冷却、小半径 |
利用数十颗毫秒脉冲星的精确计时，可以探测纳赫兹引力波：
基本原理：
引力波经过时，会改变地球和脉冲星之间空间的时间流逝，使所有脉冲星的TOA出现关联模式。
其中 是天线响应函数， 是引力波极化分量。
主要实验：
| 项目 | 望远镜 | 脉冲星数 | 观测时长 | 频率灵敏度 |
|---|---|---|---|---|
| NANOGrav | Arecibo / Green Bank | ~80 | 15年 | 1-100 nHz |
| PPTA | Parkes | ~50 | 18年 | 1-100 nHz |
| EPTA | 欧洲5台 | ~60 | 25年 | 1-100 nHz |
| CPTA | FAST | ~30 | 5年 | 1-100 nHz |
| IPTA | 联合 | 100+ | 20年 | 1-100 nHz |
2023年重大突破：NANOGrav等合作组宣布探测到纳赫兹引力波背景的强有力证据，可能源于超大质量双黑洞的并合过程。
GW170817及其电磁对应体证实了双中子星并合是宇宙中重元素（银、金、铂、铀等）的主要合成场所之一。
r-过程核合成：
在并合过程中，抛射出的中子富集物质经历快速中子俘获过程（r-process），产生约一半的比铁重的元素。
| 元素 | 原子序数 | 典型产物 | 衰变产物 | 观测已确认 |
|---|---|---|---|---|
| 锶 | 38 | 稳定 | ✅（光谱） | |
| 钇 | 39 | 稳定 | ✅ | |
| 碲 | 52 | 稳定 | 暂未确认 | |
| 铯 | 55 | 暂未确认 | ||
| 金 | 79 | 稳定 | 理论推测 | |
| 铂 | 78 | (等) | 稳定 | 理论推测 |
| 钍 | 90 | 长寿命放射性 | 理论推测 | |
| 铀 | 92 | 长寿命放射性 | 理论推测 |
千新星：
并合抛射物的放射性衰变加热产生光学/红外瞬态信号——千新星（Kilonova）：
| 问题 | 描述 | 可能解决途径 |
|---|---|---|
| 物态方程 | 超高密度下核物质的精确行为 | NICER + LIGO联合约束 + 重离子对撞实验 |
| 内核组成 | 是否存在夸克物质、超子、介子凝聚 | 质量-半径精确测量 + 冷却观测 |
| 最大质量 | 中子星不能超过的精确质量上限 | 更多高精度脉冲星质量测量 |
| 磁星起源 | 为什么某些中子星具有极端磁场 | 磁流体动力学模拟 + 磁星种群统计 |
| 脉冲辐射机制 | 射电到伽马的辐射产生精确位置和机制 | 多波段联合观测 + 等离子体模拟 |
| 星震机制 | 自转突变（glitch）的详细物理 | 毫秒分辨率计时 + 超流涡旋模型 |
| 并合产物 | 双中子星并合后是稳定中子星还是黑洞 | 更多引力波事件 + 电磁对应体 |
| 脉冲星计时引力波 | 纳赫兹引力波背景的详细来源 | 更长的计时基线 + 更多脉冲星 |
| 夸克星存在性 | 是否存在真正的奇异夸克星 | 半径-质量精确测量 + 冷却数据 |
| 中子星-黑洞双星 | 这类系统是否存在，比例如何 | 引力波探测 + 巡天观测 |
| 数据库 | 内容 | 网址 |
|---|---|---|
| ATNF脉冲星目录 | 最大脉冲星数据库（3300+颗） | https://www.atnf.csiro.au/research/pulsar/psrcat/ |
| 脉冲星计时数据 | NANOGrav/EPTA/PPTA公共数据 | https://data.nanograv.org/ |
| Chandra源目录 | X射线源数据 | https://cxc.harvard.edu/cda/ |
| Fermi伽马脉冲星目录 | 伽马射线脉冲星数据 | https://heasarc.gsfc.nasa.gov/W3Browse/fermi/fermipc |
| Gravitational Wave O4目录 | LIGO-Virgo-KAGRA引力波事件 | https://gracedb.ligo.org/ |
| 中子星质量数据库 | 中子星质量汇集 | https://stellarcollapse.org/nsmasses |