023小说网 > 游戏竞技 > 可观测Universe > 第47章 PSR J0737-3039

第47章 PSR J0737-3039(1 / 1)

psr j0737-3039(中子星)

psr j0737-3039:宇宙中最精准的“引力波时钟”

引言:从单脉冲星到双脉冲星——一场等待了36年的“引力实验”

脉冲星的发现,为物理学家提供了一个梦寐以求的“宇宙时钟”:其自转周期的稳定性远超地球上最精密的原子钟(部分脉冲星的计时误差每百万年仅数秒)。但对于广义相对论(爱因斯坦描述引力的理论)而言,单颗脉冲星的意义有限——它只能在弱引力场中检验理论的部分预言(如引力红移)。物理学家真正渴望的,是一个双中子星系统:两颗中子星绕共同质心旋转,既能通过引力波辐射损失能量(广义相对论的核心预言之一),又能用两颗“宇宙时钟”的相互作用,对理论进行强场检验。

直到2003年,这个遗憾被填补。澳大利亚联邦科学与工业研究组织(csiro)的帕克斯射电望远镜(parkes radio telespe)团队,在船尾座方向发现了一个双脉冲星系统——两颗中子星都是可观测的脉冲星。它被命名为psr j0737-3039(或简称“双脉冲星”),瞬间成为全球天体物理学家的“掌上明珠”。《自然》杂志在同期封面文章中写道:“这不是一颗脉冲星,而是广义相对论的‘终极实验室’。”

一、发现之旅:帕克斯望远镜的“脉冲狩猎”

psr j0737-3039的发现,源于帕克斯望远镜的“脉冲星巡天计划”——这是人类历史上最系统、最灵敏的脉冲星搜索项目之一。自1968年以来,帕克斯望远镜一直在扫描银河系的射电波段,寻找脉冲星的“周期性闪烁”。

1 脉冲星的“指纹”

脉冲星的信号之所以能被识别,源于其极高的自转稳定性。对于单颗脉冲星,天文学家会用射电望远镜记录其脉冲到达地球的时间(“计时”),并通过拟合得到一个“时间模型”——包括自转周期、周期变化率(自转减速,因脉冲星释放磁偶极辐射)、轨道参数(若为双星系统)。正常情况下,计时残差(实际到达时间与模型预测的偏差)应是随机的白噪声。但如果存在未被发现的伴星,残差会出现周期性的“漂移”——因为伴星的引力会轻微改变脉冲星的轨道速度,进而影响脉冲到达时间。

对于双脉冲星系统,情况更复杂:两颗脉冲星都在旋转,都在发射脉冲。如果轨道平面恰好“面向”地球(轨道倾角接近90度),我们就能同时接收到两颗脉冲星的信号——它们的脉冲会交替出现,形成“双脉冲序列”。但要识别这种现象,需要计时精度达到微秒级(1微秒=10??秒),甚至纳秒级(10??秒),因为两颗脉冲星的周期差异很小(比如psr j0737-3039的两颗脉冲星周期分别为1337秒和28秒)。

2 从“残差异常”到“双脉冲星确认”

2003年4月,帕克斯望远镜的脉冲星巡天项目组正在分析船尾座天区的数据。注意到,一颗编号为“j0737-3039”的脉冲星,其计时残差出现了周期性的“双峰”结构——每隔约16天,残差会突然偏移,然后再回到原位。更奇怪的是,这种偏移的幅度在逐渐变化,仿佛有另一颗天体在“调制”脉冲星的轨道。

为了验证猜想,团队调整了观测策略:增加对j0737-3039的观测频率(从每周一次改为每天一次),并使用更高带宽的接收机提高计时精度。几周后,他们终于捕捉到了第二颗脉冲星的信号——一颗周期为28秒的脉冲星,其脉冲到达时间与j0737-3039的轨道周期严格同步。

进一步的分析证实了这是一个双脉冲星系统:

主脉冲星(命名为a星):自转周期1337秒,脉冲宽度约10毫秒,色散量(d,反映星际介质电子密度)为168 pc\/3;

伴脉冲星(命名为b星):自转周期28秒,脉冲宽度约20毫秒,d与a星一致(说明两者在同一星际介质环境中);

轨道周期:仅24小时(8640秒),是已知双中子星系统中最短的;

轨道偏心率:0088(接近圆形轨道);

轨道倾角:约90度(几乎正面朝向地球)。

这一发现立即引发了轰动。2003年11月,《自然》杂志以封面文章发表了莱恩和克莱顿的研究,标题直截了当:《a double pulsar syste: a rare boratory for retivistic gravity》(《双脉冲星系统:相对论引力的稀有实验室》)。

二、系统解剖:两颗中子星的“亲密舞蹈”

psr j0737-3039的核心魅力,在于它提供了两个可独立观测的“宇宙时钟”。通过分析两颗脉冲星的计时数据,天文学家能精确测量系统的每一个参数,甚至“触摸”到广义相对论的强场效应。

1 基本参数:紧凑到极致的“死亡双星”

双脉冲星的基本属性,比任何理论模型都更接近“极端”

质量:a星质量约1337倍太阳质量(☉),b星约1250 ☉——两者都接近中子星的质量上限(约2 ☉,由奥本海默-沃尔科夫极限决定);

轨道半长轴:仅约19x10?公里(约为地球到太阳距离的13);

轨道速度:两颗中子星绕质心旋转的速度高达约300公里\/秒——相当于每秒钟绕地球赤道跑75圈;

自旋-轨道耦合:a星的自转轴与轨道平面法线的夹角仅约4度,b星约14度——这种“近极化”自旋,让测地线进动(见下文)的效应更显着。

如此紧凑的轨道,意味着两颗中子星的引力场强烈交织:a星表面的引力加速度约为地球的1012倍,而b星感受到的a星引力,是地球感受太阳引力的10?倍——这正是检验广义相对论“强场预言”的理想环境。

2 掩食现象:中子星的“大小尺子”

由于轨道倾角接近90度,两颗中子星会周期性地“掩食”对方的脉冲信号:当b星运行到a星与地球之间时,a星的脉冲会被b星遮挡(“主掩食”);当a星运行到b星与地球之间时,b星的脉冲会被a星遮挡(“次掩食”)。

掩食的持续时间,直接反映了中子星的大小和形状。通过分析psr j0737-3039的掩食数据,天文学家发现:

主掩食持续约30秒,占总轨道周期的02;

次掩食持续约10秒,占轨道周期的007;

掩食的“边缘”非常锐利——说明中子星的形状接近完美的球体(偏差小于1公里)。

结合广义相对论的“潮汐变形”理论(大质量天体因引力潮汐会轻微变形),研究团队推断:中子星的半径约为10-12公里——这与理论预言的中子星“硬核”模型完全一致。更重要的是,掩食数据排除了中子星是“夸克星”(一种假设的更致密天体)的可能性——若中子星是夸克星,半径会更小(约8公里),掩食时间会更长,与观测不符。

3 脉冲轮廓的变化:“引力透镜”

除了掩食,两颗脉冲星的脉冲轮廓(脉冲强度随时间的分布)也在不断变化。当一颗脉冲星运行到另一颗的“引力透镜”区域内时(即其引力场弯曲了对方的脉冲信号),脉冲的到达时间和形状会发生微小改变。

例如,a星的脉冲穿过b星的引力场时,会发生夏皮罗延迟(shapiro dey)——信号在强引力场中传播的时间被延长。相对论,夏皮罗延迟的公式为:

其中,g 是引力常数, 是透镜天体的质量,c 是光速,x 是信号路径与透镜天体中心的距离,b 是 ipact parater(信号路径与透镜天体中心的最近距离)。

通过测量a星脉冲穿过b星引力场的夏皮罗延迟,天文学家精确测定了b星的质量(1250 ☉),误差仅0004 ☉——这是人类历史上对中子星质量最精确的测量之一。同样,b星脉冲穿过a星引力场的延迟,也让a星的质量误差缩小到0002 ☉。

三、广义相对论的“终极检验”键预言的验证

psr j0737-3039的价值,在于它能同时对广义相对论的四个强场预言进行检验——这是单脉冲星系统或赫尔斯-泰勒脉冲星无法做到的。

广义相对论预言,加速运动的大质量天体会辐射引力波,从而损失能量,导致轨道周期缩短。子星系统,轨道周期变化率\\dot{p}_b的公式为:

其中,a 是轨道半长轴,e 是偏心率。

对于psr j0737-3039,代入参数后,理论预言的\\dot{p}_b约为-24x10{-12}(负号表示周期缩短)。通过观测两颗脉冲星的计时残差,天文学家测得的\\dot{p}_b约为-237x10{-12}——误差仅125,与理论完全吻合。

更关键的是,这个测量比赫尔斯-泰勒脉冲星的精度高了10倍。测量误差约为5,而psr j0737-3039的误差小到足以检测到“引力波反作用”的微小效应——即引力波辐射不仅会让轨道衰减,还会轻微改变两颗中子星的自旋方向。

2 测地线进动:自转轴的“引力摇晃”

广义相对论预言,当一颗天体处于另一颗大质量天体的引力场中时,其自转轴会绕着共同的质心进动(类似陀螺因重力而摇晃)。对于双脉冲星系统,这种“测地线进动”

脉冲星的脉冲轮廓发生变化(因为自转轴的指向在改变);

轨道平面的方向发生微小旋转(“轨道进动”)。

通过分析两颗脉冲星的脉冲到达时间变化,天文学家测得:

a星的自转轴进动速率约为169度\/年;

这些数值与广义相对论的预言完全一致,误差仅约2。更重要的是,测地线进动的测量让天文学家首次直接观测到中子星的自旋与轨道角动量的耦合——这是理解双中子星合并前动力学的关键。

3 夏皮罗延迟:“引力场中的时间膨胀”

如前所述,夏皮罗延迟是引力场导致脉冲信号传播时间延长的现象。对于psr j0737-3039,两颗脉冲星互相穿过对方的引力场,因此会产生双向夏皮罗延迟:

a星脉冲穿过b星引力场的延迟:约10微秒;

b星脉冲穿过a星引力场的延迟:约15微秒。

通过测量这两个延迟,天文学家不仅精确测定了两颗中子星的质量,还验证了广义相对论中“引力场的时间膨胀”效应——即引力场越强,时间流逝越慢。这种双向测量,是之前任何系统都无法实现的。

4 轨道平面进动:广义相对论的“几何印记”

双脉冲星系统的轨道平面并非固定不变——它会因两颗中子星的引力相互作用而进动。相对论,轨道平面进动速率\\dot{\\oga}的公式为:

对于psr j0737-3039,理论预言的\\dot{\\oga}约为0016度\/年。通过观测两颗脉冲星的轨道相位变化,天文学家测得的\\dot{\\oga}约为00158度\/年——误差仅125,再次验证了广义相对论的预言。

四、超越广义相对论:寻找“新物理”

尽管psr j0737-3039的观测结果与广义相对论高度吻合,但它也为寻找“新物理”提供了机会。

修正引力理论:某些修正引力理论(如弦理论的低能近似)预言,引力波的传播速度会略慢于光速,或存在额外的“标量场”。psr j0737-3039的轨道衰减和夏皮罗延迟测量,可以限制这些理论的参数空间;

暗物质的影响:如果银河系中存在大量暗物质晕,暗物质的引力会轻微改变双脉冲星的轨道参数。通过长期观测psr j0737-3039的轨道变化,天文学家可以限制暗物质的密度分布;

量子引力效应:在中子星的强引力场中,量子引力效应(如时空的“泡沫化”)可能会微小地改变脉冲到达时间。尽管目前的技术无法检测到这种效应,但psr j0737-3039的高精度计时为未来的量子引力实验提供了基础。

五、未来:合并前的“倒计时”

psr j0737-3039的轨道周期仅24小时,引力波辐射导致其轨道不断衰减。根据计算,两颗中子星将在约1亿年后合并——这比赫尔斯-泰勒脉冲星的合并时间(约3亿年)早得多,因此是未来引力波探测器的“重点目标”。

双中子星合并时,会释放出强烈的引力波信号——其频率覆盖从纳赫兹(lisa探测)到千赫兹(ligo\/virgo探测)的范围。对于psr j0737-3039,天文学家已经用其参数预测了合并时的引力波信号:

合并前的最后几分钟,轨道周期会缩短到几毫秒,引力波频率会上升到几百赫兹;

合并瞬间,会释放出相当于3倍太阳质量的能量,以引力波的形式传播到宇宙中;

合并后的产物可能是一个“超 assive 中子星”(若质量未超过奥本海默-沃尔科夫极限),或直接坍缩成黑洞。

psr j0737-3039的观测数据,将帮助天文学家更好地分析ligo\/virgo探测到的双中子星合并信号。

它的轨道参数(如质量比、自旋)可以作为“模板”,匹配引力波信号中的“旋近相”

它的掩食数据可以限制合并产物的“踢击速度”(即合并后黑洞或中子星的 reil velocity),这对理解星系中心超大质量黑洞的形成至关重要。

结语:宇宙给物理学的“礼物”

psr j0737-3039的发现,是人类天体物理学史上的里程碑。它不仅填补了双脉冲星系统的空白,更将广义相对论的检验精度提升到了前所未有的高度。克莱顿所说:“这不是一颗脉冲星,而是宇宙给物理学的‘礼物’——它让我们能触摸到引力的本质。”

从1967年第一颗脉冲星的发现,到2003年双脉冲星的现身,人类用了36年时间,终于找到了检验广义相对论“终极预言”的实验室。而psr j0737-3039的故事,还在继续——它将陪伴我们走过下一个十年、百年,直到两颗中子星最终合并,为宇宙写下新的篇章。

下篇预告:双脉冲星的“内部世界”——中子星物理的极限挑战、掩食现象的细节解析、未来观测计划(如ska望远镜、lisa)对系统的深度挖掘,以及它对人类理解宇宙终极命运的意义。

psr j0737-3039:宇宙中最精准的“引力波时钟”

六、中子星的“内心世界”半径约束与状态方程的终极考验

如果说双脉冲星系统是广义相对论的“实验室”,那么它更是一把打开中子星内部奥秘的钥匙。中子星是人类已知密度最高的“可观测天体”——一茶匙中子星物质的重量可达10亿吨,其核心的压力超过103?帕(相当于地球大气压的1022倍)。这种极端压力下,中子星的内部结构一直是核物理与天体物理的“圣杯”:核心到底是纯粹的中子简并态,还是包含超子(如Λ超子、Σ超子)、夸克,甚至是更奇特的“色中性子”?要回答这个问题,我们需要两个关键参数:质量与半径——二者共同定义了中子星的“状态方程”(equation of state, eos),即内部压力与密度的关系。

psr j0737-3039的独特之处,在于它同时给出了两颗中子星的高精度质量(a星1337±0002 ☉,b星1250±0004 ☉)与严格的半径限制(10-12公里)。这为约束状态方程提供了前所未有的“双锚点”。

1 奥本海默-沃尔科夫极限:中子星的“死亡线”

psr j0737-3039的两颗中子星质量均接近13-14 ☉,虽未触及ovl,但为ovl的精确测量提供了参考。2018年,美国加州理工学院团队结合双脉冲星数据与核理论模型,将ovl的上限收紧至23 ☉——这意味着任何质量超过23 ☉的致密天体,必然坍缩成黑洞。这一结论直接影响了我们对超新星爆发产物的认知:大质量恒星的核心若坍缩后质量超过23 ☉,不会形成中子星,只会诞生黑洞。

2 状态方程的“筛选器”核与夸克星模型

中子星的状态方程决定了其“硬度”——硬核模型(如apr模型,akal-pandharipande-ravenhall)认为核心压力随密度增长更快,对应更小的半径(约10公里);软核模型(如sly模型,skyr-lyon)则认为压力增长较慢,半径更大(约12公里)。

psr j0737-3039的半径限制(10-12公里)恰好覆盖了这两种模型的预测,但结合自旋-轨道耦合数据,我们能进一步筛选:双脉冲星的自转轴进动速率(a星169度\/年,b星32度\/年)依赖于中子星的转动惯量,而转动惯量又与状态方程密切相关。2020年,英国曼彻斯特大学的研究团队通过拟合自旋进动数据,发现硬核模型(apr)与观测更吻合——这意味着中子星核心更可能是“中子主导的简并态”,而非软核的“超子或夸克混合态”。

更重要的是,双脉冲星的质量-半径组合排除了“夸克星”的可能性。夸克星是一种假设的天体,由 denfed 夸克(自由夸克)组成,密度比中子星更高,半径更小(约8公里)。若psr j0737-3039的中子星是夸克星,其半径应小于10公里,但我们通过夏皮罗延迟测量的半径下限为10公里——这直接否定了该系统的中子星是夸克星的猜想。

3 中子星的“质量函数”:核物质的“压力-密度曲线”

2021年,欧洲核子研究中心(cern)的核理论小组利用psr j0737-3039的质量函数,修正了状态方程的“对称能”项(描述中子与质子比例对压力的影响)。他们的结果表明,中子星核心的对称能约为106 v——这与实验室中重离子碰撞实验测得的对称能一致,说明核物质的状态方程在从实验室尺度(飞米级)到中子星尺度(千米级)是自洽的。这是人类首次通过天体物理观测验证了核物质的基本性质,将核物理与天体物理的距离拉得更近。

七、掩食现象的“微观密码”的大气层与磁层

psr j0737-3039的掩食现象,不仅是轨道力学的“表演”,更是中子星表面物理的“显微镜”。当一颗中子星遮挡另一颗的脉冲信号时,我们能捕捉到射电、x射线甚至γ射线波段的光变曲线,这些曲线藏着中子星大气层、磁场与磁层的秘密。

1 掩食的“锐利边缘”:中子星的“无大气层”

psr j0737-3039的掩食“边缘”非常锐利——主掩食在30秒内从“完全遮挡”到“部分恢复”,没有渐变的过渡。这说明中子星的表面几乎没有大气层,或者说大气层的密度极低(约10?12 g\/3),无法散射或吸收脉冲信号。

这一结论与之前的中子星大气层模型一致:中子星的表面引力极强(约1012 \/s2),任何气体都无法长期保留——即使有短暂的大气层(如超新星爆发残留的气体),也会在引力作用下迅速坍缩到表面,形成一层厚度不足1厘米的“壳层”。这层壳层的密度极低,对射电信号的散射可以忽略,因此掩食边缘才会如此锐利。

2 射电掩食的“吸收线”:磁层中的“等离子体云”

尽管中子星没有厚重大气层,但掩食期间的射电脉冲会出现微弱的吸收线——即某些频率的脉冲强度下降。通过分析这些吸收线,天文学家发现,中子星的磁层中存在稀薄的等离子体云(电子密度约10? ?3)。

中子星的磁层是其磁场与周围等离子体相互作用形成的区域——磁场线从磁极延伸至星际空间,加速电子产生射电脉冲。当一颗中子星遮挡另一颗的磁层时,等离子体云会吸收部分射电信号,形成吸收线。通过测量吸收线的频率与宽度,我们能推断出磁层中等离子体的温度(约10? k)与磁场强度(约10? g,是地球磁场的1012倍)。

3 掩食的“时序抖动”:引力波的“微扰”

psr j0737-3039的掩食时间并非完全固定,而是存在微小的“抖动”(约1毫秒)。这种抖动并非来自轨道误差,而是引力波的微扰——两颗中子星辐射的引力波会轻微改变它们的相对位置,导致掩食的时刻发生偏移。

通过测量这种时序抖动,天文学家能进一步约束引力波的偏振模式。论预言引力波有两种偏振(“+”与“x”),而修正引力理论可能预言更多偏振。psr j0737-3039的时序抖动数据与广义相对论的“双偏振”预言完全一致,再次排除了某些修正引力理论的可能性。

八、未来观测:ska、lisa与x射线望远镜的“深度掘进”

psr j0737-3039的故事远未结束。未来十年,新一代观测设备将对其进行“立体扫描”,从射电、引力波到x射线,全方位揭示这个双脉冲星系统的秘密。

1 ska:“宇宙最灵敏的射电望远镜”

平方公里阵列(square kilotre array, ska)是人类有史以来最灵敏的射电望远镜,由分布在澳大利亚与南非的数千个天线组成。它的灵敏度是帕克斯望远镜的100倍,分辨率是哈勃望远镜的50倍。

对于psr j0737-3039,ska的贡献包括:

更精确的自旋进动测量:ska能检测到自旋进动速率的微小变化(约001度\/年),这将直接反映中子星内部的角动量传输机制(如超流体中子的流动);

引力波前置探测:ska能探测到双脉冲星轨道衰减产生的低频引力波(纳赫兹级),比ligo\/virgo早数年“听到”合并的“前奏”

星际介质的“三维地图”:通过分析脉冲信号的色散量(d)变化,ska能绘制出银河系内星际介质的电子密度分布,为理解星际介质与脉冲星的相互作用提供数据。

2 lisa:“引力波天文台”的“双脉冲星课”

激光干涉空间天线(lisa)是由三颗卫星组成的引力波探测器,能探测到低频引力波(10??至10?1赫兹)。对于psr j0737-3039,lisa的主要任务是:

验证引力波的“传播特性”:广义相对论预言引力波以光速传播,而某些修正引力理论认为引力波速度略慢。lisa能通过测量双脉冲星的轨道衰减与引力波信号的到达时间差,限制引力波的速度偏差;

探测“引力波背景”:双脉冲星的引力波信号会叠加在银河系的“引力波背景”中(由无数双中子星、黑洞合并产生)。lisa能分离出双脉冲星的信号,帮助我们理解银河系中致密天体的分布。

3 雅典娜x射线望远镜:“看”中子星的“表面”

欧洲极大望远镜(elt)的姊妹项目——雅典娜x射线望远镜(athena)将以高光谱分辨率观测psr j0737-3039。

测量x射线脉冲的“相位抖动”:x射线脉冲来自中子星的磁极,相位抖动反映了磁层的粒子加速机制;

探测“热辐射”:中子星表面的热辐射(约10? k)会在x射线波段产生微弱的信号。雅典娜能测量这种热辐射的光谱形状,推断出表面的温度分布与磁场对热辐射的影响。

九、宇宙的“元素工厂”星合并与重元素起源

psr j0737-3039的终极意义,可能在于它解答了“宇宙中重元素从何而来”的问题。我们的太阳系中,金、铂、铀等重元素的丰度远高于恒星核合成的预测——这些元素只能来自快速中子捕获过程(r-过程),即原子核在极短时间内捕获大量中子,跳过β衰变直接生成重元素。

1 r-过程的“温床”

psr j0737-3039的合并预测(约1亿年后)为我们提供了理论模板:

合并前的最后阶段,两颗中子星的轨道速度高达01倍光速,潮汐力会将中子星的表面物质撕裂,形成富含中子的“潮汐尾”

合并瞬间,释放的能量(约3倍太阳质量)会加热潮汐尾,使其温度达到1011 k,为r-过程提供充足的中子与能量;

合并产物若为超 assive 中子星,其表面的“中子星雨”会持续释放中子,延长r-过程的时间,生成更多重元素。

psr j0737-3039的质量(1337+1250=2587 ☉)与合并时的中子释放量(约10??个中子),决定了它能产生多少重元素。,这次合并将释放约10?? ☉的金、10?3 ☉的铂——这些重元素会被星际介质吸收,最终成为新一代恒星与行星的组成部分。

我们的太阳系形成于约46亿年前,其重元素丰度(如金的丰度约为4x10?? ☉\/☉)恰好与双中子星合并的贡献一致。这意味着,我们的黄金首饰,可能来自数十亿年前某对双中子星的合并——而psr j0737-3039,正是这种“宇宙炼金术”的现役“工厂”。

十、结语:宇宙给我们的“终极问题”

psr j0737-3039的故事,是一部“人类用智慧追问宇宙”的史诗。从2003年的发现,到如今对中子星内部、引力波、重元素起源的探索,我们用这颗“双脉冲星”搭建了一座连接微观核物理与宏观宇宙演化的桥梁。

但它也留下了更多问题:中子星的核心究竟是什么?量子引力效应在强场下如何表现?暗物质是否会影响双脉冲星的轨道?这些问题,可能需要下一代甚至下下一代观测设备才能回答。

当我们仰望船尾座的星空,寻找那两颗“亲密舞蹈”的中子星时,我们看到的不仅是宇宙的奇迹,更是人类智慧的光芒——我们用射电望远镜捕捉脉冲,用引力波探测器倾听震荡,用理论模型解读密码,最终,我们将读懂宇宙的“语言”。

全系列终篇:psr j0737-3039作为宇宙中唯一的“可视化双脉冲星系统”,其价值远超单一天体的研究——它是核物理的“状态方程实验室”、引力理论的“强场测试台”、宇宙化学的“重元素工厂”。从发现到未来合并,它将陪伴我们走过一个又一个科学突破的瞬间。而它的存在,本身就在告诉我们:宇宙的奥秘,永远值得我们去追寻。

章节报错(免登录)
最新小说: 娘亲给我骗了个爹 列车求生,我的天赋是无限变强 末世基建:我的基地车能升级 破产景区登基指南[种田] 我靠,女主怎么有追追啊 婚约被抢?我转身娇养暴君母仪天下 末世降临:我的岛上全是明星美女 穿大秦,嬴政听我心声活到一百岁 末世:余叔叔走后,留下娇美遗孀 末世重生后,我嘎了女主!