统一 Hirota 双线性方法下爱因斯坦真空场方程各种引力孤子剖面的研究

《Chaos, Solitons & Fractals》：The investigation of various gravitational soliton profiles for Einstein’s vacuum field equations under the unified Hirota bilinear approach

　　

**摘要**

爱因斯坦真空场方程（EVFE）在广义相对论中占据核心地位，支配无物质源区域的时空几何。本文分析的非线性偏微分方程是完整真空场方程的一个特定（1+1）维约化，通过施加Guilfoyle度量试探解并令$R_{tt}=0$得到；所有结果均针对此约化模型。尽管具有物理重要性，但针对此约化EVFE的精确非线性波解目录尚不包含空间局域化的多分量结构和呼吸子型结构：在本文之前，尚无此类解被构造。我们通过应用Hirota双线性方法并结合策略性选取的多项式试探函数，在统一框架内获得了一大类精确解析解：孤束、孤束1条纹、孤束2条纹、孤束-周期、流形周期、巨浪、Ma呼吸子、Kuznetsov–Ma呼吸子、广义呼吸子、Akhmediev呼吸子以及对应的标准巨浪。我们进一步构造了孤束、周期和扭波分量之间的相互作用解，并通过等高线分析展示了其三尺度空间结构。对每个解类生成了三维曲面轮廓和等高线图，通过系统参数变化，揭示了背景水平、波数和速度参数如何控制局域化几何、峰值幅度和传播方向。这些解与引力波物理、非线性光学和海洋波动力学直接相关，并为数值相对论和波散射研究提供基准构型。

**引言**

爱因斯坦场方程（EFE）是广义相对论（GR）的数学基石，将时空曲率与物质和能量分布联系起来。在其真空形式中，场方程简化为Ricci张量为零的条件$R_{\mu\nu}=0$，编码引力扰动在真空时空中的传播。John Archibald Wheeler的经典总结概括了其物理内涵："时空告诉物质如何移动，物质告诉时空如何弯曲"。

EVFE的物理范围异常广泛。在太阳系尺度上，EVFE预测了水星近日点进动、太阳引起的光线偏折和Shapiro时间延迟效应，所有这些都已被航天器和脉冲星计时以亚百分比精度确认[1], [2], [3]。在宇宙学尺度上，Schwarzschild和Kerr度量等精确真空解描述了黑洞外部、引力透镜以及激光干涉仪可探测的并合事件几何。LIGO于2015年首次探测到引力波的里程碑式发现提供了直接证据，证明EVFE支配引力波涟漪在天文距离上的产生和传播[4], [5]。本文研究的控制方程源自Guilfoyle度量[6]，该度量描述爱因斯坦-麦克斯韦框架下静态带电球形物质分布。令Ricci张量的$(t,t)$分量为零，将其简化为两个自变量的单一非线性偏微分方程。

非线性波方程的精确解的意义超越数学优雅。在引力物理背景下，它们提供了：(i)验证数值相对论代码的基准时空构型；(ii)引力波散射和记忆效应微扰研究的种子度量；(iii)其稳定性特征约束早期宇宙和致密双星并合中物理场景的原型波形。孤束波是理性局域化且代数衰减的解，自然地从双线性形式中产生，代表平坦背景上空间紧凑的引力激发。呼吸子是空间上或时间上周期性的结构，在固定位置有限时间内集中能量，类似于在光纤和深水实验中观测到的由调制不稳定性驱动的极端波。在呼吸子族中，Ma呼吸子在空间上周期、在时间上局域，在单一时刻达到最大幅度后返回背景。Akhmediev呼吸子则在时间上周期、在空间上局域，作为调制不稳定性的非线性饱和出现。两者在适当极限情况下收敛到Peregrine孤子。巨浪突然从准均匀背景中出现，携带三倍于背景水平的幅度，模型化了与随机引力背景涨落相关的间歇性极端幅度事件。有限振幅的引力波本质上是非线性的：它们通过Weyl张量作为进一步曲率的源[4]，而线性化近似恰好在巨浪型解所描述的高振幅区域失效。本工作正是针对这一真正的非线性区域。

关于EVFE的解析工作具有显赫但狭窄的历史。Kaur和Wazwaz[7]应用Painlevé分析和Lie对称方法建立可积性性质并获得相似性约化。Friedrich和Nagy[8]研究了真空方程边值问题的适定性，而Rutz[9]将EVFE扩展到Finsler几何设置，Heinicke和Hehl[10]研究了Kerr型真空解。El-Nabulsi和Golmankhaneh[11]引入分数阶导数公式，Jyoti和Kumar[12]通过对称性约化构造了非静态精确解。对其他可积系统的多孤束、呼吸子和巨浪解的平行研究[13]-[20]展示了Hirota双线性方法从单一统一试探解产生丰富精确解族的强大能力。近期进展包括(3+1)维KdV型方程的巨浪层次[21], [22]以及通过直接符号方法获得的新中心控制巨浪解[23], [24], [25]。然而，对于EVFE，之前仅有孤束型解可用[26]，呼吸子层次、巨浪族和多分量相互作用解仍完全未被探索。

EVFE的双线性形式与聚焦非线性Schr?dinger（NLS）方程共享结构联系。两个方程支持相同的拓扑解族：孤束波、Ma呼吸子、Kuznetsov–Ma呼吸子、Akhmediev呼吸子和Peregrine巨浪，因为两者都允许Hirota双线性形式。关键区别在于色散关系：EVFE双线性算子$D_xD_t+D_x^4$给出$\omega(k)=-k^4$，而NLS给出$\omega(k)=k^2$。尽管色散关系发生变化，拓扑解结构得以保留。这种与NLS层次的联系在第2节中明确阐述，解释了我们在此推导的解为何与非线性和海洋波物理直接相关，其中呼吸子和巨浪族已被实验确认[27]。

本文推导的解具有跨多个活跃领域的应用。在引力波物理中，孤束解提供紧凑的非辐射度规扰动波形；巨浪解模型化了与随机引力背景相关的极端瞬态度规放大；呼吸子解描述引力波介质中调制不稳定性的有限幅度饱和。这些是已在实验中观测到的结构的引力对应物：Akhmediev呼吸子和Peregrine孤子已在光纤[27]和深水波池[28], [29]中得到确认。其他研究小组通过符号双线性方法研究了非线性演化方程中相关孤子的传播[21]-[25]。

在本工作中，我们填补了这一空白。从EVFE的双线性形式出发，我们通过选择适当的多项式和指数试探函数系统性地构造了十二个不同的解类。全文呈现了三维曲面轮廓和等高线图，每个参数在物理动机范围内变化。论文组织如下。第3节推导基本孤束解$u_3$。第4-5节构造一和二条纹孤束相互作用。第6节呈现孤束-周期解。第7节给出流形周期解。第8节推导巨浪$u_{13}$。第9节构造三分量相互作用。第10-14节呈现完整的呼吸子层次和标准巨浪。第15节与先前工作进行比较，第16节总结结果。为提高可读性，原第3节（独立摘要表）已重新命名并融入解呈现流程，讨论已重新组织为聚焦小节。

**分节片段**

**控制方程与双线性约化**

在Guilfoyle时空中[6]，函数$u(x,t)$起引力势的作用：$g_{tt}=-au^2+2u_{xx}$编码度规的推迟函数，而非对角项$g_{xt}=-u$耦合时间和空间坐标。物理上，$u$大的区域对应更强的时空曲率；在线性化极限$u\to m+\epsilon\phi$，$|\epsilon|\ll 1$时，方程(2)简化为$\phi$的线性波动方程，恢复弱场引力波区域。完整非线性方程(2)支配

**孤束解**

应用试探解$u=2b(\ln f)_x+m$并使用二次形式$f=\eta_1^2+\eta_2^2+a_7$，$\eta_i=a_ix+a_{i+1}t$，集合III给出孤束解：
$$u_3=m+2b\frac{2a_4(a_5t+a_4x)-2a_4a_5\frac{(a_2t-a_4a_5x)}{a_2}}{a_4^2\left[a_2^2(12m+1)+(a_5t+a_4x)^2+\frac{a_2^2(a_2t-a_4a_5x)^2}{a_2^2}\right]}$$

孤束波是理性局域化且在$(x,t)$平面所有方向代数衰减的非线性波结构。它们携带比一般孤子更高的传播能量，在理解海洋和非线性光学中的极端波事件中起重要作用。

**孤束1条纹**

使用$f=\eta_1^2+\eta_2^2+a_7+b_1e^{k_1x+k_2t}$，集合I给出孤束1条纹解$u_4$[26]。

孤束1条纹解源于理性孤束分量与单个指数扭状条纹之间的相互作用[18], [20]。辅助函数$f=\eta_1^2+\eta_2^2+a_7+b_1e^{k_1x+k_2t}$结合了正定二次项（产生孤束）和沿由比值$k_2/k_1$设定的条纹方向衰减（或增长）的指数项。参数$b_1$控制条纹强度；当$b_1\to$

**孤束2条纹**

使用$g=\eta_1^2+\eta_2^2+a_7+b_1e^{k_1x+k_2t}+b_2e^{k_3x+k_2t}$，集合II给出解$u_7$：
$$u_7=m+2b\frac{2a_1ia_5t^3+a_1x+2a_4(a_5t+a_4x)}{U(k_i,t)}$$
其中$U(k_i,t)=b_1e^{k_2t}+b_2e^{k_2t}+56a_1^2bk_2+(a_5t+a_4x)^2+i\frac{a_5t^3+a_1x^2}{2}$。孤束2条纹解向$g$引入第二条指数条带，创建支持三种不同分量相互作用的波形：理性孤束核心和两个扭状条纹。虚系数$a_1=(2i/3)a_4$和$a_2=(i/3)a_5$自然地出现在双线性形式的相容性条件中。

**孤束-周期解**

使用$f=\eta_1^2+\eta_2^2+a_7+b_1\cos(k_1x+p_1t)$，$p_1=(11k_1^3-4k_1m+2k_1)/2$，集合I给出$u_9$。

孤束-周期解源于将指数试探解$f$中的指数条纹替换为余弦调制[30]。时间频率$p_1=(11k_1^3-4k_1m+2k_1)/2$由双线性约束唯一确定，确保理性孤束轮廓与周期性背景之间的精确相容性。与条纹解不同，余弦项在两个方向上无限延伸，创建

**流形周期解**

使用$f=e^{-d_1\eta_1}+h_1e^{d_1\eta_1}+h_2\cos(d_2\eta_2)$，集合II给出：
$$u_{12}=-2a_4b\frac{d_2h_2\sin d_2}{a_4+(2a_4-11a_4^3d_2^2)t^2+a_4x}\cdot\frac{1}{e^{-d_1(a_3+a_2t)}+e^{d_1(a_3+a_2t)}+h_2\cos(\cdots)}$$

流形周期解采用混合试探解，在一个坐标方向上组合双曲（cosh型）函数，在另一方向上组合余弦[30]。变换$f=e^{-d_1\eta_1}+h_1e^{d_1\eta_1}+h_2\cos(d_2\eta_2)$创建非均匀周期晶格，其峰值高度由双曲包络调制。令$h_1=0$简化为单侧指数，而$h_1=1$

**巨浪解**

使用$f=\xi_1^2+\xi_2^2+a_7+b_1\cosh(k_1z+p_1t)$，集合I给出巨浪解$u_{13}$。巨浪是突然从准均匀背景中出现并同样迅速消失的极端波事件[19], [31]。它们在数学上以空间和时间的理性局域化为特征，与孤子（空间局域但无限传播于时间）和周期解（空间延伸）区分。双曲试探解$f=\xi_1^2+\xi_2^2+a_7+b_1\cosh(k_1z+p_1t)$

**孤束、周期和1扭波的相互作用**

使用$f=\eta_1^2+\eta_2^2+a_7+b_1\cosh(k_1z+p_1t)+b_2e^{k_3z+k_4t}$，集合I给出$u_{14}$。

相互作用解是本工作中推导的最复杂波形，同时编码三种不同的物理机制。复合试探解$f=\eta_1^2+\eta_2^2+a_7+b_1\cosh(k_1z+p_1t)+b_2e^{k_3z+k_4t}$包含：(i)生成理性孤束的二次部分；(ii)创建时间周期山脊的cosh部分；以及(iii)引入单调背景梯度的指数扭部分。三个分量相互作用

**Ma呼吸子及其巨浪**

使用$f=1+a_1+e^{ip_1x+k_1t+k_2}+e^{-ip_1x+k_1t+k_2}+b_1e^{2(k_1t+k_2)}$，集合I给出Ma呼吸子$u_{15}$：
$$u_{15}=2a_1b\frac{e^{k_2-k_1t}(-ie^{-ip_1x}p_1+ie^{ip_1x}p_1)}{1+\frac{9}{4}a_1^2e^{2(k_2-k_1t)}+a_1e^{k_2-k_1t}(e^{-ip_1x}+e^{ip_1x})-im}$$

Ma呼吸子（MB）是空间周期、时间局域的非线性波结构，以Ma命名[31]。它可以被视为调制不稳定性过程的有限幅度极限：平面波背景被扰动，在$t=0$时增长到最大幅度，然后返回背景。复指数

**Kuznetsov–Ma呼吸子**

使用$f=e^{-p_1(x-b_1t)}+a_1\cos(p(x-b_1t))$，$b_1=-1$，集合I给出$u_{18}$（方程(48)）。

Kuznetsov–Ma（KM）呼吸子是Ma呼吸子的时间对应物：它在$t$上周期、在空间上局域，在固定空间位置产生一系列规则脉冲[31], [32], [33]。Kuznetsov–Ma呼吸子首次由Kuznetsov（1977）和独立地由Ma（1979）推导，现已在光纤和深水波池中实验充分表征[32], [33]。试探解$f=e^{-p_1(x-b_1t)}+a_1$

**广义呼吸子**

使用[34]：
$$\Psi=e^{it(1-4c)}\frac{\cosh(\sigma t)+2c\cos(\rho x)}{2c\cos(\rho x)-\cosh(\sigma t)}+i\sigma\frac{\sinh(\sigma t)}{2c\cos(\rho x)-\cosh(\sigma t)}$$
其中$c=\frac{1}{16}+\frac{i3}{948}$，通过$V=2bc(6/\kappa\ln\Psi)_x+m$获得广义呼吸子解$u_{20}$。

广义呼吸子（GB）解从方程(11)的复试探解$\Psi$中获得，遵循Singh等人[34]的工作。与Ma和KM呼吸子在精确一个方向上局域不同，GB结合空间周期性（来自$\cos(\rho x)$）和时间周期性（来自$\cosh(\sigma t)$和$\sinh(\sigma t)$），创建空间和时间双重周期波。

**Akhmediev呼吸子**

使用[27]：
$$\psi=p_0\frac{(1-4a)\cosh(bx)+ib\sinh(bx)+2a\cos(\omega_{\text{mod}}T)}{2a\cos(\omega_{\text{mod}}T)-\cosh(bx)}$$
集合II给出$u_{22}$。

Akhmediev呼吸子（AB）是调制不稳定性循环的极限情况，其中调制频率$\omega_{\text{mod}}$亚临界，即$\omega_{\text{mod}}<\omega_c$[27]。在此区域中，扰动从连续波背景增长，在$t=0$时达到最大值，然后返回。

AB是空间周期性的，其波数为b=[8a(1?2a)]^(1/2)，由调制深度2a=1?(ω_mod/ω_c)2决定。

**标准rogue波解**

利用文献[31]中的ψ=?e^(it)(1?4(1+2it)/(1+4x2+4t2))，参数组I给出SRW解：

u??=?80(54i2+15)/[9a(58i2+3)]?384bc(1+2it)/[xk(1+4t2+4x2)2]?1+4(1+2it)/(1+4t2+4x2)。

标准rogue波（SRW）是Akhmediev呼吸子解和KM呼吸子解在调制频率与局域化长度均趋于零时的有理极限[31]。Peregrine型试探解ψ=?e^(it)[1?4(1+2it)/(1+4x2+4t2)]所给出的解在x和t方向均呈代数衰减，分母1+4x2+4t2确保了无奇异性。

**与引力物理的内在联系**

EVFE(2)控制Guilfoyle时空族中的度规势u(x,t)[6]。因此，每个解都描述了度规分量g_tt∝u2的一种不同形变。我们强调，下面给出的物理解释是定性和探索性的：它们是数学类比，旨在激发解的分类并将其与已知的引力现象联系起来。这些联系并未得到曲率不变量、测地线偏差、潮汐……的计算支持。（注：原文此处被截断）

**结束语**

我们通过Hirota双线性方法推导了爱因斯坦真空场方程的十二类精确解。其中十类对该方程而言是全新的：lump 2条带解u?、lump-周期解u?、流形周期解u??、rogue波u??、三分量相互作用解u??、Ma呼吸子u??、Kuznetsov–Ma呼吸子u??、广义呼吸子u??、Akhmediev呼吸子u??以及标准rogue波u??。

主要贡献包括：(i)从单一控制方程出发，通过统一的双线性推导方法获得全部十二类解；……（注：原文此处被截断）

**CRediT作者贡献声明**

潘朝红（Chaohong Pan）：形式分析。
Allah Ditta：软件、形式分析。
Irfan Mahmood：撰写——初稿。
Abaid ur Rehman Virk：撰写——审阅与编辑。
Nadia Cheemaa：撰写——初稿、方法论。
Farruh Atamurotov：软件。
Alisher Abduvokhidov：监督、概念化。
Rana Muhammad Zulqarnain：撰写——审阅与编辑、撰写——初稿。

**利益冲突声明**

作者声明没有已知的竞争性经济利益或个人关系可能影响本文所报告的工作。

**致谢**

湖南省大数据分析与应用省大学重点实验室，长沙，410205，中国。

作者：Chaohong Pan | Allah Ditta | Irfan Mahmood | Abaid ur Rehman Virk | Nadia Cheemaa | Farruh Atamurotov | Alisher Abduvokhidov | Rana Muhammad Zulqarnain