Research Paper Notes on Time-Delay Interferometry

Research Paper Notes on Time-Delay Interferometry

本文档除了包括推导, 疑惑 外,做 读书重点 的记录 $$$$

文献列表

 * Time-Delay Interferometry, Living Rev. Relativity, 17, (2014) 6, by Massimo Tinto and Sanjeev V. Dhurandhar
 * Time-Delay interferometry, Living Rev. Relativity, 24, (2021) 1, by Massimo Tinto and Sanjeev V. Dhurandhar
 * Elimination of Clock Jitter Noise inSpaceborn Laser Interferometers, arXiv:gr-qc/0012013v2, by Ronald W. Hellings
 * Effect of filters on the time-delay interferometry residual laser noise for LISA, arXiv:1811.01575v8, by Jean-Baptiste Bayle et al
 * Clock-jitter reduction in LISA time-delay interferometry combinations, arXiv:2005.02430v2, by Olaf Hartwig and Jean-Baptiste Bayle
 * Geometric Time Delay Interferometry, arXiv:gr-qc/0504145v2, by Michele Vallisneri

Time-Delay Interferometry, Living Rev. Relativity, 17, (2014) 6, by Massimo Tinto and Sanjeev V. Dhurandhar
本文是综述arXiv:gr-qc/0409034v2的2014更新版本.讨论探测引力波的空间迈克尔逊干涉仪的激光器相位噪声的消除TDI方法.本质上该方法利用匹配光程的方法来抵消未知形式的噪声,除了表面上简单的代数加减外,抵消方案即方程(21)或(28)的解,对应了从四个多项式多项式环到一个多项式环的映射$$\mathcal{R}^4\to\mathcal{R}$$的核.这进而可以通过Grobner基,环理想的生成元,和syzygies模等数学概念来系统的得到问题的解.后者的数学基础是值得深入学习的.

(7-9)

这里的推导文中的叙述并不清晰,可以参考比较(25)附近的推导.(8)是对光路中的激光相位噪声进行两次延时并相减后剩余的非零部分,它对应了噪声一次延时后的贡献.相当于(25)中算符乘积中的一次项.

(21)

因为(20)中的$$p$$可以是任何噪声形式,故(20)并不对应矩阵的齐次方程从而得到(21)的行列式为零,而是(21)本身必须为零.

其中$$(q,q')$$是$$\mathcal{D}_i\ (i=1,2,3)$$的多项式,因为多项式本身构成环,而$$\mathcal{D}_i$$是相互对易的,故这些多项式也构成环.

(28)

这个表达式(28)左边通过四个多项式的线性组合得到一个多项式,故是一个从四个多项式多项式环到一个多项式环的映射$$\mathcal{R}^4\to\mathcal{R}$$.

如文中所述,方程(28)要求等式右边为零,故它的解就是上述映射的核.称为sizigies模(module).

这个问题的数学本质就是通过环论的工具来求解上述方程.

(29)

数学方案的第一步是讨论由集合(29)为生成元构造的理想.因为有四个生成元(即将这四个生成元乘以环的任何元素构成的理想,这区别于更简单情况下主理想通过一个生成元来产生的例子),故理想通过生成元来产生的形式是$$\sum v_i u_i$$,其中$$u_i$$就是由(29)决定的生成元(用符号$$u_i$$没有异议因为生成元本身必然属于通过其构造的理想$$u_i\in\mathcal{U}$$),而$$v_i\in\mathcal{R}$$是环中的任何元素,容易证明这样得到的集合满足理想的定义.上述定义的理想可以由某一组更为精简的生成元基来生成,这组基底的一个选组方案就是Grobner基.

文中指出Grobner基是由单一生成元产生的主理想情况的推广,因为这里涉及多个生成元.

对于单一生成元的简单情况,文中给出了这样的例子.考虑仅由一个(而非本文中涉及的三个)变量定义的多项式环.我们同样考虑由任意指定的$$m$$个环的元素$$(g_1,g_2,\cdots,g_m)$$为生成元定义的理想$$\mathcal{G}$$,理想的元素形式为$$\sum_{i=1}^m v_i g_i$$.其中$$v_i$$就是该多项式环的元素.现在要考虑某多项式$$f$$是否属于上述理想.我们注意到,按多项式环的定义,多项式系数本身也是某个环的元素,我们考虑一个简单的情况,即多项式系数是整数,即多项式系数属于整数环.这意味着,对生成元集合的线性组合可以得到其"最大公约数",它是某个确定的多项式$$g$$.所以,若$$f$$是$$g$$乘以其他换元素得到的,换言之,$$f$$可以被$$g$$整除(余数为零),我们就可以通过对生成元的线性组合得到$$g$$,然后得到$$f$$,即$$f$$属于环$$\mathcal{G}$$.反之,若$$f$$属于环$$\mathcal{G}$$,它必然可以表达为$$\sum_{i=1}^m v_i g_i$$的形式,但因为$$g$$是$$g_i$$的公约数,我们可以把它们各自都写成公约数乘以换元素的形式,这样提出公约数后$$g$$,显然形式上$$f$$能被$$g$$整除.上述论证意味着,一个变量的多项式环的任意理想必然是主理想,理想只有一个生产元,上述例子中理想$$\mathcal{G}$$的生产元就是$$g$$.而文中指出,多于多变量多项式环,上述简单结论不再成立,这样才有Groner基的概念.

文中指出,对多变量多项式环,虽然任意环元素都可以写成$$\sum_{i=1}^m v_i u_i$$的形式,但是把换元素用生产元$$u_i$$来除,并把得到的余数用下一个生成元来除,一直到穷尽所有生成元,最后的余数未必为零.这是因为保证余数最小的原则导致了,用第一个生成元$$u_1$$除以后的余数并不一定是$$\sum_{i=2}^m v_i u_i$$.而Broner基确保上述连续除法最后得到的余数为零.

(30)

具体计算结果得到这里的Grobner基的数目为3个(比之前(29)中给出系数数目少1).更为直观的细节,参见(115-116)附近的讨论.

(31)

进而,我们可以计算对应sizigies模的生成元,从而得到问题的解.因为模(module)的简单例子就是定义在标量场上的矢量线性空间,对应的模的每个生成元具有6个分量(即$$q,q'$$的分量总数).线性独立的生成元有4个.具体细节参见附录的讨论.

(118)

由(116)易知$$Af=(1-dc)f=0$$,取转置为$$f^T A^T=0$$.比较方程(114),知$$A^T$$的每一行都满足(114),它们的线性组合也必然满足.

我们指出,上述构造的基础是矩阵$$dc$$的存在性.这可以简单的由基的定义得到证明.参考(29)处采用的符号,因为Grobner基是理想的基,故对于任意一组环的元素$$v_i\in\mathcal{R}$$,我们得到一个理想中的元素$$r=fv\equiv \sum_i v_i u_i\in \mathcal{U}$$,必然存在另外一组对应的环的元素$$w_i\in\mathcal{R}$$满足$$r=gw\equiv\sum_i w_ig_i$$.反之,对任意一组换的元素$$w'_i\in\mathcal{R}$$,得到一个理想中的元素$$r'=gw'\equiv\sum_i w'_ig_i$$,必然能够找到另外一组对应的环的元素$$v'_i\in\mathcal{R}$$满足$$r'=fv'\equiv \sum_i v'_i u_i\in \mathcal{U}$$.由上述第一个条件,我们取$$v_i=(1,0,0,0)$$,对应的$$w_i$$就是$$d$$矩阵的第一行的元素,由此类推得到$$d$$矩阵.而类似的上述第二个条件可用于决定$$c$$矩阵.

这里给出的结果中$$a_2$$的第二个元素有打字错误,因为$$1-z^2$$而非$$z(1-z^2)$$.

(120)

这里具体描述了如何构造另外三个生成元.这里,矩阵$$b^*$$把$$f$$转换为Grobner基的S多项式,但并不为零.具体推导参见本综述2021版的(50)之前的讨论.

通过上述具体构造过程我们直观的看到为何Grobner基可以用于确定sizigies模的生成元.

本文档除了包括推导, 疑惑 外,做 读书重点 的记录 $$$$

Time-Delay interferometry, Living Rev. Relativity, 24, (2021) 1, by Massimo Tinto and Sanjeev V. Dhurandhar
这是综述的2021更新版本,对环论部分有较大更新.

(50)

这里给出了本文之前版本附录(120)的具体证明.比较(40),容易看到由(49)定义的$$b^*$$矩阵的行也是sizigies模的生成元.

(53)

这里讨论了syzygy的第二模的定义,它由sizigies第一模生成元的线性不独立性导致的.而对应的第二模的基是一维的.

而进一步推广,我们可以对希尔伯特的syzygy定理给出描述. 它形象的对应了文中给出的正合序列. 注意到中间部分的映射$$\mathcal{R}\to\mathcal{R}^4$$前后分别对应了syzygy第二模与第一模的生成元的维度,相应映射的核是零元素. 而中间的$$\mathcal{R}^4\to\mathcal{R}^4$$中的两个维度$$4$$分别是指由syzygy第一模式生成元的维度和(34)等式左边定义的理想的生成元的数目,相应的核正是(53). 而中后的$$\mathcal{R}^4\to\mathcal{R}$$中的两个维度$$4$$分别是指由(34)等式左边定义的理想的生成元的数目和等式右边的唯一多项式环,相应的核正是(34).

一般的,从syzygy第$$(k+1)$$模到第$$(k)$$模的映射,是一个理想.它是由前者的维度的环空间直积空间映射到后者维度的环空间直积空间.具体的,第$$(k+1)$$模的生成元的任一线性组合(决定的理想)按规则对应到第$$(k)$$模的生成元的某一个组合.显然,这个组合的组合系数的数目是$$(k)$$模生成元的数目.按定义,这个映射的相,就是由第$$(k)$$模的生成元定义的理想的核.由第$$(k)$$模的生成元定义的理想,是更低一阶的syzygy模的理想的核.而相应的,往另一个方向追溯,第$$(k+1)$$模的生成元定义了一个理想,如果这些生成元不是(resolved)线性独立的,那么这个理想有一个不平庸的核,它是比它更高一阶的syzygy模定义的映射的相.

上述讨论中的一个关键点是,对于多项式构成的对易环,对应某理想的完备的生成元并不一定是线性独立的.显然上述讨论中涉及的理想的核的生成元仍然满足理想的定义,且是上述理想的子理想(一般的情况下,模的核的生成元也构成模,且是子模).但是因为多项式环在除法运算下并不封闭(虽然可以引入整除和余数的定义,即两个环元素的除法并不一定仍是环的元素),容易理解环的生成元之间并不线性独立,但是因为不能保证整除,无法通过他们线性组合为零的关系把其中任何一个元素用其他元素表达(具体的,没有任何一个组合系数是其他系数的公因子)从而剔除出生成元集合.生成元线性相关的关系显然是用上述生成元表达理想的某种冗余.为了进一步消除这种冗余,我们可以试图找到生成元线性的关系系数(即生成元构成理想的映射的核)的生成元.上述做法显然可以不断重复,直至某阶的生成元是线性独立(映射的核为平庸)为止,这就是引入上述正和序列背后的动机.

本文档除了包括推导, 疑惑 外,做 读书重点 的记录 $$$$

Elimination of Clock Jitter Noise inSpaceborn Laser Interferometers, arXiv:gr-qc/0012013v2, by Ronald W. Hellings
本文是引入边带(sideband)以实现钟噪声消除方案的第一篇工作.

(7-10)

这是在时间域计算两束激光的相位差. 其中(10)等式右边的第一项,一方面它来自于将(8-9)直接代入(7),即非相对多普勒效应,另一方面,它与频率的相对论多普勒效应,即频率作为四矢的零分量在洛伦兹变换下的变化形式自洽.

(12ab)

这里加入了钟噪声的效应. 这里(12a)就是(10)的结果.而(12b)是在激光频率中加入钟频率,同时考虑了钟(频)噪声的结果.

我们注意到(12a)中等式右边的第一项$$(\nu_2-\nu_1-V_{12}\nu_2)t$$,其中$$(\nu_2-\nu_1)$$来自激光器频率差,文中指出,它的大小为300MHz左右,而$$V_{12}\nu_2$$是相对运动带来的多普勒效应的修正,数量级更小,仅为即MHz. 但是它们给所测量的拍频增加了一个很大的偏移,在下面,通过实验手段控制$$a_{21}$$的数值来消去它对测量的影响.

(13ab)

按之前讨论,这里在实验上的边带信号中通过相位锁定的办法,使得$$a_{21}f_1=\nu_2-\nu_1-V_{12}\nu_2$$,这样就能得到一般文献中考虑边带后的信号公式.

这里的问题是,实验上的相位锁定必然存在一些误差,这些误差乘以时间$$t$$故随着时间单调增加.我们需要证明它并不随着时间积累影响到引力波信号测量. 实际上,按arXiv:gr-qc/0112059v2的(4.20)的推导和结果,我们发现引力波对拍频的影响是对激光频率的相对值,故只要差值$$[a_{21}f_1-(\nu_2-\nu_1-V_{12}\nu_2)]$$足够小,它不会对引力波信号有影响.

本文档除了包括推导, 疑惑 外,做 读书重点 的记录 $$$$

Effect of filters on the time-delay interferometry residual laser noise for LISA, arXiv:1811.01575v8, by Jean-Baptiste Bayle et al
本文考虑反锯齿滤波,即高频滤波对TDI消除激光噪声方案的残余噪声的影响.残余噪声被表达为滤波与延时算符对易子的形式.

本文档除了包括推导, 疑惑 外,做 读书重点 的记录 $$$$

Clock-jitter reduction in LISA time-delay interferometry combinations, arXiv:2005.02430v2, by Olaf Hartwig and Jean-Baptiste Bayle
本文考虑TDI及边带消除激发及钟噪声方案的残余噪声.残余噪声被表达为延时算符对易子的形式.

本文档除了包括推导, 疑惑 外,做 读书重点 的记录 $$$$

Geometric Time Delay Interferometry, arXiv:gr-qc/0504145v2, by Michele Vallisneri
相对于"代数"方法,作者把本文提出的搜寻TDI组合的方法称为几何方法.

虽然物理上是很有启发性的,但是这个方法的数学基础应该并不清楚.

具体做法我们以Fig.4中unequal arm Michelson写出(9)为例. 这个方法要穷举所有从某给定飞船出发(TDI解显然对出发点具有轮换对称性)沿着实线和虚线按时间演化回到相同飞船的光信号路径,每两个飞船间的线段都带箭头.正确的TDI组合存在其中,显然,Fig.4中对应的图都是TDI解,且满足上述特性. 然后按以下规则得到该图的相应贡献. 每一截连接两个相邻飞船(link)的线段(臂)都给出其部分贡献,最后的贡献是和的形式. 每个线段的贡献就是其对应的相邻两飞船间光信号的干涉(差). 实线贡献为正,虚线贡献为负. 每个线段的贡献具有延时,其延时就是从最终飞船(也是出发飞船)位置逆着时间倒退直到出发位置对应的的延时算符的反复叠加. 顺着箭头和逆着箭头对应顺时针和逆时针的臂长(虽然是同一臂长).

(6)

这就是测量数据,它在图中对应相邻飞船间信号的干涉测量结果.其中包含延时不同的两个噪声和引力波信号.

Fig.2

这里没有用虚线,但是给出了最简单的测量信号(6)的差的表达式. 注意到(6)中等式左边的时间自变量是接收到光信号的非常的时间.

(7-8)

这是Fig.3对应不等臂第一代Michelson方案的具体计算.

(9)

这是Fig.4对应的不等臂第二代Michelson方案的具体计算.

本文档除了包括推导, 疑惑 外,做 读书重点 的记录 $$$$