留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

面向天基引力波探测的时间延迟干涉技术

王登峰 姚鑫 焦仲科 任帅 刘玄 钟兴旺

王登峰, 姚鑫, 焦仲科, 任帅, 刘玄, 钟兴旺. 面向天基引力波探测的时间延迟干涉技术[J]. 中国光学. doi: 10.37188/CO.2020-0098
引用本文: 王登峰, 姚鑫, 焦仲科, 任帅, 刘玄, 钟兴旺. 面向天基引力波探测的时间延迟干涉技术[J]. 中国光学. doi: 10.37188/CO.2020-0098
WANG Deng-feng, YAO Xin, JIAO Zhong-ke, REN Shuai, LIU Xuan, ZHONG Xing-wang. Time-delay interferometry for space-based gravitational wave detection[J]. Chinese Optics. doi: 10.37188/CO.2020-0098
Citation: WANG Deng-feng, YAO Xin, JIAO Zhong-ke, REN Shuai, LIU Xuan, ZHONG Xing-wang. Time-delay interferometry for space-based gravitational wave detection[J]. Chinese Optics. doi: 10.37188/CO.2020-0098

面向天基引力波探测的时间延迟干涉技术

doi: 10.37188/CO.2020-0098
基金项目: 中国航天科技集团公司自主研发项目(Y18-JTKJCX-01)
详细信息
    作者简介:

    王登峰(1979—),男,陕西西安人,硕士,高级工程师,2011年于西安电子科技大学获得硕士学位,主要从事星间精密测量与时间同步技术研究与产品开发。E-mail:wangdf@cast504.com

    姚 鑫(1990—),男,陕西渭南人,博士,工程师,2019年于清华大学获得博士学位,主要从事星间激光通信测距与星载相位计研究与产品开发。E-mail:yaox1@cast504.com

    焦仲科(1988—),男,甘肃白银人,博士,工程师,2017年于中国科学院光电技术研究所获得博士学位,主要从事星间精密测量研究。E-mail:jiaozk12@cast504.com

    任 帅(1988—),男,陕西榆林人,硕士,工程师,2016年于空间技术研究院获得硕士学位,主要从事星载高精度相位计研究。E-mail:rens@cast504.com

    刘 玄(1988—),男,陕西西安人,硕士,高级工程师,2013年于西安电子科技大学获得硕士学位,主要从事星载微波-激光混合测量链路研究。E-mail:liux73@cast504.com

    钟兴旺(1967—),男,山西忻州人,博士,研究员,1990年于西北工业大学获得学士学位,2010年于中国空间技术研究院获得博士学位。主要从事星间链路组网通信与基线测量设备研制。E-mail:zhongxw@cast504.com

    通讯作者: 姚鑫,中国空间技术研究院西安分院(陕西省西安市长安区东长安街504号,邮编:710100),联系电话:17791282559,029-89253514,E-mail:yaox1@cast504.com
  • 中图分类号: O436.1

Time-delay interferometry for space-based gravitational wave detection

Funds: Supproted by China Aerospace Science and Technology Corporation’s Independent Research and Development Project (Y18-JTKJCX-01)
More Information
  • 摘要: 在天基引力波探测任务中,利用激光干涉仪对无拖曳的检验质量块间实现十皮米量级的位移测量精度。其中,激光源频率噪声和时钟频率噪声是两项主要噪声。在欧洲主导的LISA(Laser Interferometer Space Antenna)引力波探测项目中,利用时间延迟干涉技术(Time-delay Interferometry,TDI),通过对三星上的十二组相位测量值进行延迟和线性组合,构造出臂长相等的干涉仪,从而消除激光源噪声以及光学平台位移噪声。为了消除时钟噪声,将时钟信号倍频到GHz,通过相位调制的方式加载在星间激光链路上,最终从时钟边带拍频信号中提取出时钟噪声,并在TDI的数据组合中将时钟噪声项消除。为了实现TDI中的时间延迟处理,要求对星间绝对距离进行精确测量。因此,在TDI机制中,星间激光链路同时实现位移测量,时钟边带调制和绝对距离测量三个功能。其中,后两个功能分别大约消耗载波激光功率的10%和1%。LISA项目针对TDI技术的地面论证结果表明,TDI技术对激光源和时钟的噪声抑制分别达到了109和6×104倍。TDI技术对中国引力波探测项目及其它天基激光精密测量任务具有重要的参考价值。
  • 图  1  天基迈克尔逊干涉仪

    Figure  1.  Space-based Michelson interferometer

    图  2  TDI原理示意图

    Figure  2.  Diagram of the principles behind the TDI

    图  3  LISA光学系统示意图(已获参考文献[19]授权© American Physical Society)

    Figure  3.  The illustration of optical system in the LISA mission (Reprinted with permission from ref. [19] © American Physical Society)

    图  4  X型数据组合物理模型

    Figure  4.  The physics modal of the X type data combination

    图  5  当星座旋转时,${L'_i} \ne {L_i}$

    Figure  5.  In the case of the constellation rotation, ${L'_i} \ne {L_i}$

    图  6  时钟边带调制,EOM:电光调制器

    Figure  6.  Clock sideband modulation, EOM: electro-optic modulator

    图  7  探测端信号处理示意图,NCO:数控振荡器

    Figure  7.  TSignal processing illustration of photon detection, NCO: Numerically Controlled Oscillator

    图  8  在载波相位上调制数据码和伪随机码用于星间数据传输和绝对距离测量,EOM:电光调制器

    Figure  8.  The phase of the carrier is modulated by data codes and pseudo-random codes for the inter-satellite data transmission and absolute distance measurement, EOM: electro-optic modulator

    图  9  数据码和伪随机码参数设计(已获参考文献[24]授权© The Optical Society)

    Figure  9.  Parameter design of data codes and pseudo-random codes (Adapted with permission from ref. [24] © The Optical Society)

    图  10  光载波经过时钟信号和扩频调制后干涉测量结果(已获参考文献[24]授权© The Optical Society)

    Figure  10.  Interference measurement results when the optical carriers are modulated by the clock signal and spread spectrum (Reprinted with permission from ref. [24] © The Optical Society)

    图  11  LISA干涉仪测试床位移测量结果(已获参考文献[22]授权© American Physical Society)

    Figure  11.  The displacement measurment results in the LISA interferometry test bed (Reprinted with permission from ref. [22] © American Physical Society)

  • [1] ABBOTT B P, ABBOTT R, ABBOTT T D et al. Observation of gravitational waves from a binary black hole merger[J]. Physical Review Letters, 2016, 116(6): 061102. doi:  10.1103/PhysRevLett.116.061102
    [2] ABBOTT B P, ABBOTT R, ABBOTT T D et al. GW170817: Observation of gravitational waves from a binary neutron star inspiral[J]. Physical Review Letters, 2017, 119(16): 161101.
    [3] BARACK L, CARDOSO V, NISSANKE S, et al. Black holes, gravitational waves and fundamental physics: a roadmap[J]. Classical and Quantum Gravity, 2019, 36(14): 143001. doi:  10.1088/1361-6382/ab0587
    [4] DANZMANN K, RUDIGER A. LISA technology-concept, status, prospects[J]. Classical and Quantum Gravity, 2003, 20(10): S1-S9. doi:  10.1088/0264-9381/20/10/301
    [5] WANNER G. Space-based gravitational wave detection and how LISA Pathfinder successfully paved the way[J]. Nature Physics, 2019, 15(3): 200-202. doi:  10.1038/s41567-019-0462-3
    [6] ARMANO M, AUDLEY H, BAIRD J, et al. LISA pathfinder platform stability and drag-free performance[J]. Physical Review D, 2019, 99(8): 082001. doi:  10.1103/PhysRevD.99.082001
    [7] ABICH K, BRAXMAIER C, GOHLKE M, et al. In-orbit performance of the GRACE follow-on laser ranging interferometer[J]. Physical Review Letters, 2019, 123(3): 031101. doi:  10.1103/PhysRevLett.123.031101
    [8] 王智, 沙巍, 陈哲, 等. 空间引力波探测望远镜初步设计与分析[J]. 中国光学,2018,11(1):131-151. doi:  10.3788/co.20181101.0131

    WANG ZH, SHA W, CHEN ZH, et al. Preliminary design and analysis of telescope for space gravitational wave detection[J]. Chinese Optics, 2018, 11(1): 131-151. (in Chinese doi:  10.3788/co.20181101.0131
    [9] LUO Z R, LIU H SH, JIN G. The recent development of interferometer prototype for Chinese gravitational wave detection pathfinder mission[J]. Optics &Laser Technology, 2018, 105: 146-151.
    [10] 刘河山, 高瑞弘, 罗子人, 等. 空间引力波探测中的绝对距离测量及通信技术[J]. 中国光学,2019,12(3):486-492. doi:  10.3788/co.20191203.0486

    LIU H SH, GAO R H, LUO Z R, et al. Laser ranging and data communication for space gravitational wave detection[J]. Chinese Optics, 2019, 12(3): 486-492. (in Chinese doi:  10.3788/co.20191203.0486
    [11] LI Y P, LIU H SH, ZHAO Y, et al. Demonstration of an ultraprecise optical bench for the Taiji space gravitational wave detection pathfinder mission[J]. Applied Sciences, 2019, 9(10): 2087. doi:  10.3390/app9102087
    [12] 罗子人, 张敏, 靳刚, 等. 中国空间引力波探测“太极计划”及“太极1号”在轨测试[J]. 深空探测学报,2020,7(1):3-10.

    LUO Z R, ZHANG M, JIN G, et al. Introduction of Chinese space-borne gravitational wave detection program “Taiji” and “Taiji-1” satellite mission[J]. Journal of Deep Space Exploration, 2020, 7(1): 3-10. (in Chinese
    [13] LUO J, CHEN L SH, DUAN H Z, et al. TianQin: a space-borne gravitational wave detector[J]. Classical and Quantum Gravity, 2016, 33(3): 035010. doi:  10.1088/0264-9381/33/3/035010
    [14] SHI CH F, BAO J H, WANG H T, et al. Science with the TianQin observatory: preliminary results on testing the no-hair theorem with ringdown signals[J]. Physical Review D, 2019, 100(4): 044036. doi:  10.1103/PhysRevD.100.044036
    [15] ZHANG J Y, MING M, JIANG Y Z, et al. Inter-satellite laser link acquisition with dual-way scanning for Space Advanced Gravity Measurements mission[J]. Review of Scientific Instruments, 2018, 89(6): 064501. doi:  10.1063/1.5019433
    [16] FALLER J E, BENDER P L. A possible laser gravitational wave experiment in space[C]. Program and Abstracts of Second International Conference on Precision Measurement and Fundamental Constants (PMFC-II), National Bureau of Standards, 1981: 689-690.
    [17] TINTO M, DHURANDHAR S V. Time-delay interferometry[J]. Living Reviews in Relativity, 2014, 17: 6. doi:  10.12942/lrr-2014-6
    [18] ESTABROOK F B, TINTO M, ARMSTRONG J W. Time-delay analysis of LISA gravitational wave data: elimination of spacecraft motion effects[J]. Physical Review D, 2000, 62(4): 042002. doi:  10.1103/PhysRevD.62.042002
    [19] TINTO M, ESTABROOK F B, ARMSTRONG J W. Time-delay interferometry for LISA[J]. Physical Review D, 2002, 65(8): 082003. doi:  10.1103/PhysRevD.65.082003
    [20] TINTO M, HARTWIG O. Time-delay interferometry and clock-noise calibration[J]. Physical Review D, 2018, 98(4): 042003. doi:  10.1103/PhysRevD.98.042003
    [21] DE VINE G, RABELING D S, SLAGMOLEN B J J, et al. Picometer level displacement metrology with digitally enhanced heterodyne interferometry[J]. Optics Express, 2009, 17(2): 828-837. doi:  10.1364/OE.17.000828
    [22] DE VINE G, WARE B, MCKENZIE K, et al. Experimental demonstration of time-delay interferometry for the laser interferometer space antenna[J]. Physical Review Letters, 2010, 104(21): 211103. doi:  10.1103/PhysRevLett.104.211103
    [23] SPERO R, BACHMAN B, DE VINE G, et al. Progress in interferometry for LISA at JPL[J]. Classical and Quantum Gravity, 2011, 28(9): 094007. doi:  10.1088/0264-9381/28/9/094007
    [24] ESTEBAN J J, GARCÍA A F, BARKE S, et al. Experimental demonstration of weak-light laser ranging and data communication for LISA[J]. Optics Express, 2011, 19(17): 15937-15946. doi:  10.1364/OE.19.015937
    [25] HEINZEL G, ESTEBAN J J, BARKE S, et al. Auxiliary functions of the LISA laser link: ranging, clock noise transfer and data communication[J]. Classical and Quantum Gravity, 2011, 28(9): 094008. doi:  10.1088/0264-9381/28/9/094008
    [26] BRAUSE N C. Auxiliary function development for the LISA metrology system[D]. Hannover: Gottfried Wilhelm Leibniz Universität, 2018.
    [27] ISLEIF K S. Laser interferometry for LISA and satellite geodesy missions[D]. Hannover: Gottfried Wilhelm Leibniz Universität, 2018.
    [28] SCHWARZE T S, BARRANCO G F, PENKERT D, et al. Picometer-stable hexagonal optical bench to verify LISA phase extraction linearity and precision[J]. Physical Review Letters, 2019, 122(8): 081104. doi:  10.1103/PhysRevLett.122.081104
    [29] SHEARD B S, GRAY M B, MCCLELLAND D E, et al. Laser frequency stabilization by locking to a LISA arm[J]. Physics Letters A, 2003, 320(1): 9-21. doi:  10.1016/j.physleta.2003.10.076
    [30] SHADDOCK D, MCKENZIE K, SPERO R, et al.. LISA frequency control white paper[R]. 2009. (查阅所有网上资料, 未找到出版地出版者信息, 请联系作者补充)
    [31] KARLEN L, KUNDERMANN S, TORCHEBOEUF N, et al.. Laser system for the LISA mission[C]. 2019 Joint Conference of the IEEE International Frequency Control Symposium and European Frequency and Time Forum, IEEE, 2019: 1-2.
    [32] CORNISH N J, HELLINGS R W. The effects of orbital motion on LISA time delay interferometry[J]. Classical and Quantum Gravity, 2003, 20(22): 4851-4860. doi:  10.1088/0264-9381/20/22/009
    [33] KULLMANN J. Development of a digital phase measuring system with microradian precision for LISA[D]. Hannover: Gottfried Wilhelm Leibniz Universität, 2012.
    [34] GERBERDING O. Phase readout for satellite interferometry[D]. Hannover: Gottfried Wilhelm Leibniz Universität, 2014.
    [35] POLLACK S E, STEBBINS R T. Demonstration of the zero-crossing phasemeter with a LISA test-bed interferometer[J]. Classical and Quantum Gravity, 2006, 23(12): 4189-4200. doi:  10.1088/0264-9381/23/12/014
    [36] POLLACK S E, STEBBINS R T. A demonstration of LISA laser communication[J]. Classical and Quantum Gravity, 2006, 23(12): 4201-4213. doi:  10.1088/0264-9381/23/12/015
    [37] GRÜNING P, HALLOIN H, PRAT P, et al. Status of the eLISA on table (LOT) electro-optical simulator for space based, long arms interferometers[J]. Experimental Astronomy, 2015, 39(2): 281-302. doi:  10.1007/s10686-015-9448-z
    [38] LAPORTE B M, HALLOIN H, BRÉELLE E, et al. Status of the LISA On table experiment: a electro-optical simulator for LISA[J]. Journal of Physics:Conference Series, 2017, 840: 012014. doi:  10.1088/1742-6596/840/1/012014
    [39] FRANCIS S P, SHADDOCK D A, SUTTON A J, et al. Tone-assisted time delay interferometry on GRACE Follow-On[J]. Physical Review D, 2015, 92(1): 012005. doi:  10.1103/PhysRevD.92.012005
  • [1] 吕潇磊, 赵继广, 杜小平, 宋一铄, 张朋, 张建伟.  新型电光材料调制性能研究进展 . 中国光学, doi: 10.37188/CO.2020-0039
    [2] 赵家旺, 张运海, 王发民, 缪新, 施辛.  线扫描虚拟结构调制共聚焦显微成像 . 中国光学, doi: 10.37188/CO.2020-0120
    [3] 安玲坪, 王爽, 张耿, 李娟, 刘学斌.  大孔径静态干涉成像光谱仪径向畸变导致谱线偏移误差的校正方法 . 中国光学, doi: 10.37188/CO.2020-0084
    [4] 张石磊, 崔宇, 邢慕增, 闫斌斌.  光场成像目标测距技术 . 中国光学, doi: 10.37188/CO.2020-0043
    [5] 申帅, 张葆, 李贤涛, 张士涛.  基于伪微分和加速度反馈的航空光电稳定平台控制方法 . 中国光学, doi: 10.3788/CO.20171004.0491
    [6] 张玉莲, 储海荣, 张宏巍, 张明月, 陈阳, 李银海.  MEMS陀螺随机误差特性研究及补偿 . 中国光学, doi: 10.3788/CO.20160904.0501
    [7] 徐正平, 沈宏海, 许永森.  直接测距型激光主动成像系统发展现状 . 中国光学, doi: 10.3788/CO.20150801.0028
    [8] 梁静秋, 梁中翥, 吕金光, 秦余欣, 田超, 王维彪.  空间调制微型傅里叶变换红外光谱仪研究 . 中国光学, doi: 10.3788/CO.20150802.0277
    [9] 王竞, 李建中, 温伟峰, 陈光华, 彭其先, 陶世兴.  利用自相关方法实现光脉冲时间延迟精确测量 . 中国光学, doi: 10.3788/CO.20150802.0270
    [10] 代晓珂, 金春水, 于杰.  点衍射干涉仪波面参考源误差及公差分析 . 中国光学, doi: 10.3788/CO.20140705.0855
    [11] 孙鸽, 马冬梅, 张海涛.  干涉仪仪器传递函数的检测方法 . 中国光学, doi: 10.3788/CO.20140701.0137
    [12] 任玉, 蔡红星, 谭见瑶, 谭勇, 张喜和, 郑峰, 马文联.  声光调制光谱相机的成像漂移 . 中国光学, doi: 10.3788/CO.20130602.0179
    [13] 张绍练, 高世杰, 吴志勇.  RS码仿真与基于RiBM算法的硬件实现 . 中国光学, doi: 10.3788/CO.20130602.0171
    [14] 张春雷, 向阳, 于长淞.  高精度干涉仪针孔空间滤波器研制 . 中国光学, doi: 10.3788/CO.20130606.952
    [15] 苏志德, 史振广, 彭吉, 隋永新, 杨怀江.  菲索干涉仪中精确移相的实现 . 中国光学, doi: 10.3788/CO.20130602.0244
    [16] 李娜, 王红.  激光驾束制导仪中调制盘的设计 . 中国光学,
    [17] YANG T, LI C C, HO H P.  基于亚微米尺度干涉阵列的超小型高分辨光谱仪 . 中国光学,
    [18] 于杰.  用于相移点衍射干涉仪的加权最小二乘相位提取算法 . 中国光学,
    [19] 张海涛.  基于光学设计软件的相移点衍射干涉仪建模 . 中国光学,
    [20] 张未来, 宋克菲, 王云磊, 潘利华, 马庆军, 汪龙祺, 刘海波.  固相时间分辨荧光免疫分析仪数据采集系统设计 . 中国光学,
  • 加载中
图(11)
计量
  • 文章访问数:  101
  • HTML全文浏览量:  24
  • PDF下载量:  7
  • 被引次数: 0
出版历程

面向天基引力波探测的时间延迟干涉技术

doi: 10.37188/CO.2020-0098
    基金项目:  中国航天科技集团公司自主研发项目(Y18-JTKJCX-01)
    作者简介:

    王登峰(1979—),男,陕西西安人,硕士,高级工程师,2011年于西安电子科技大学获得硕士学位,主要从事星间精密测量与时间同步技术研究与产品开发。E-mail:wangdf@cast504.com

    姚 鑫(1990—),男,陕西渭南人,博士,工程师,2019年于清华大学获得博士学位,主要从事星间激光通信测距与星载相位计研究与产品开发。E-mail:yaox1@cast504.com

    焦仲科(1988—),男,甘肃白银人,博士,工程师,2017年于中国科学院光电技术研究所获得博士学位,主要从事星间精密测量研究。E-mail:jiaozk12@cast504.com

    任 帅(1988—),男,陕西榆林人,硕士,工程师,2016年于空间技术研究院获得硕士学位,主要从事星载高精度相位计研究。E-mail:rens@cast504.com

    刘 玄(1988—),男,陕西西安人,硕士,高级工程师,2013年于西安电子科技大学获得硕士学位,主要从事星载微波-激光混合测量链路研究。E-mail:liux73@cast504.com

    钟兴旺(1967—),男,山西忻州人,博士,研究员,1990年于西北工业大学获得学士学位,2010年于中国空间技术研究院获得博士学位。主要从事星间链路组网通信与基线测量设备研制。E-mail:zhongxw@cast504.com

    通讯作者: 姚鑫,中国空间技术研究院西安分院(陕西省西安市长安区东长安街504号,邮编:710100),联系电话:17791282559,029-89253514,E-mail:yaox1@cast504.com
  • 中图分类号: O436.1

摘要: 在天基引力波探测任务中,利用激光干涉仪对无拖曳的检验质量块间实现十皮米量级的位移测量精度。其中,激光源频率噪声和时钟频率噪声是两项主要噪声。在欧洲主导的LISA(Laser Interferometer Space Antenna)引力波探测项目中,利用时间延迟干涉技术(Time-delay Interferometry,TDI),通过对三星上的十二组相位测量值进行延迟和线性组合,构造出臂长相等的干涉仪,从而消除激光源噪声以及光学平台位移噪声。为了消除时钟噪声,将时钟信号倍频到GHz,通过相位调制的方式加载在星间激光链路上,最终从时钟边带拍频信号中提取出时钟噪声,并在TDI的数据组合中将时钟噪声项消除。为了实现TDI中的时间延迟处理,要求对星间绝对距离进行精确测量。因此,在TDI机制中,星间激光链路同时实现位移测量,时钟边带调制和绝对距离测量三个功能。其中,后两个功能分别大约消耗载波激光功率的10%和1%。LISA项目针对TDI技术的地面论证结果表明,TDI技术对激光源和时钟的噪声抑制分别达到了109和6×104倍。TDI技术对中国引力波探测项目及其它天基激光精密测量任务具有重要的参考价值。

English Abstract

王登峰, 姚鑫, 焦仲科, 任帅, 刘玄, 钟兴旺. 面向天基引力波探测的时间延迟干涉技术[J]. 中国光学. doi: 10.37188/CO.2020-0098
引用本文: 王登峰, 姚鑫, 焦仲科, 任帅, 刘玄, 钟兴旺. 面向天基引力波探测的时间延迟干涉技术[J]. 中国光学. doi: 10.37188/CO.2020-0098
WANG Deng-feng, YAO Xin, JIAO Zhong-ke, REN Shuai, LIU Xuan, ZHONG Xing-wang. Time-delay interferometry for space-based gravitational wave detection[J]. Chinese Optics. doi: 10.37188/CO.2020-0098
Citation: WANG Deng-feng, YAO Xin, JIAO Zhong-ke, REN Shuai, LIU Xuan, ZHONG Xing-wang. Time-delay interferometry for space-based gravitational wave detection[J]. Chinese Optics. doi: 10.37188/CO.2020-0098
    • 基于迈克尔逊激光干涉仪的引力波探测器开启了人类观测宇宙的新纪元。2016年以来,位于美国的臂长4 km的LIGO探测器和位于意大利的臂长3 km的VIRGO探测器联合探测到数十次致密双星的合并事件,标志着引力波天文学的迅猛发展。目前,地基引力波探测器已经持续运行数年,并且在不断改造升级以达到更高的灵敏度[1-3]。然而,由于臂长限制,地基探测器的观测窗口在1 Hz到KHz的高频段。由于在地面大幅度增加干涉仪臂长是不实际的,并且存在难以消除的低频地表震动噪声,因此LIGO等探测器无法有效观测到1 Hz以下的低频引力波信号。另一方面,在mHz波段宇宙存在着最丰富的引力波辐射源。在太空中构建具有十万公里甚至百万公里臂长的干涉仪是观测低频引力波的理想解决方案。欧洲主导的“激光干涉仪空间天线”(Laser Interferometer Space Antenna,LISA)项目便是瞄准探测和研究毫赫兹波段引力波辐射源的典型代表。在LISA任务中,三颗全同卫星组成边长为5×106 km(约16.67光秒)的等边三角形,依靠三颗卫星间的激光链路实现低频引力波的观测[4, 5]。2015年12月,欧洲航天局成功发射了一颗技术验证星LISA Pathfinder,验证了对检验质量块的无拖曳控制技术[5, 6]。2018年,基于激光干涉测距的GRACE Follow-On高精度地球重力场测量卫星发射成功,并在约220公里的低低跟踪双星编队间实现了纳米级的星间相对距离测量目标,验证了LISA项目的部分关键技术[7]。最近几年,为了抢占引力波探测科学高地,中国科学院[8-12]和中山大学[13-15]也相应开启了天基引力波探测研究项目,分别命名为“太极”和“天琴”计划,并在2019年相继成功发射了一颗技术验证星。

      在天基引力波探测中,需要对自由悬浮的检验质量块间实现十皮米量级的位移测量精度。然而,由于受轨道动力学影响,三颗卫星编队之间的基线长度最大相差2%。因此,注入到迈克尔逊干涉仪两臂的激光源频率噪声无法互相消除,从而无法有效探测到引力波导致的质量块间皮米级位移信号。另一方面,目前星载超稳晶振(Ultra-stable oscillator, USO)的频率稳定度无法支持皮米级的甚高测相精度,USO引入到相位计的噪声同样会导致无法有效探测到引力波信号。因此,激光源频率噪声和星载时钟频率噪声成为天基引力波探测任务中的两个主要问题。在2004年发布的LISA任务白皮书中,激光源噪声和时钟噪声消除问题位列69项任务风险中的前两项。为了解决天基引力波干涉仪中激光源噪声无法直接消除的问题,J. E. Faller等人首次提出了时间延迟干涉技术(Time Delay Interferometry, TDI)的概念[16]。在TDI中,对测相数据进行适当的延迟处理和线性组合,通过数据后处理的方式构造出臂长相等的干涉仪,从而消除激光源频率噪声[17, 18]。为了抑制相位计中的时钟频率噪声,M. Tinto等人提出将时钟信号调制到激光相位上,利用新的拍频信号和测相值提取出时钟噪声,从而在TDI数据组合中消除时钟噪声项[19, 20]。在后续的实验论证阶段,加州理工大学喷气推进实验室(JPL)[21-23]和德国马普引力物理研究所下属的阿尔伯特-爱因斯坦研究所(AEI)[24-28]等研究机构结合LISA的系统设计对TDI技术做了一系列地面测试,初步论证了在激光链路上完成相对位移测量,绝对距离测量和时钟噪声转移的三个功能,为天基引力波探测奠定了坚实的技术基础。

      本文的组织结构如下:第2节讨论激光源频率噪声对引力波探测的影响,分析了TDI技术抑制激光源噪声的物理图景,并结合LISA项目光路系统布局,给出了简要分析和相应的数据组合类型;在第3节,针对测相处理中的时钟噪声问题,LISA项目利用时钟边带调制的方案,消除了TDI数据组合中的时钟噪声项;第4节针对TDI数据处理过程中对星间绝对距离的测量需求,讨论了星间激光链路的附加功能:伪随机码测距和数据传输;在文章的最后一部分,概述了针对TDI技术的地面测试论证工作。希望对我国天基引力波探测任务以及其它激光干涉测量项目提供借鉴。

    • 图1所示,在天基激光干涉仪中,卫星1上的激光经过分束镜分束后分别注入到干涉仪两臂L2(卫星1和3之间)和L3(卫星1和2之间)中。激光经过长距离传输后到达远星处的检验质量块处,被锁相放大后传送回来,在卫星1上的分束镜合束后进入探测器。引力波会导致检验质量块产生微小位移,探测器将测量到激光干涉信号发生相位涨落。

      图  1  天基迈克尔逊干涉仪

      Figure 1.  Space-based Michelson interferometer

      然而,由于受轨道动力学影响,三星编队之间的基线会有一定的差别。如在LISA任务中,三星编队的中心位于日地轨道,编队基线的差别最大可达2%。因此,激光源注入到两臂的频率噪声无法在探测器处相互抵消。这样,激光源引起的频率噪声会淹没掉引力波信号引起的相位涨落。

      具体对于迈克尔逊干涉仪而言,星间位移测量精度为:

      $$\delta L = \Delta L\frac{{\delta \nu }}{{{\nu _0}}},$$ (1)

      其中,$\delta \nu $为引力波探测所关注的mHz波段的频率噪声,${\nu _0}$为激光源中心频率,$\Delta L$为干涉仪臂长差。

      假设星载激光源在mHz波段的频率噪声$\delta \nu $$300\;{\rm{ Hz}}/\sqrt {{\rm{Hz}}} $水平,激光源中心频率${\nu _0}$约为300 THz。以LISA任务为例,星间距离$L$,即干涉仪臂长,为5×109 m,臂长差$\Delta L$在108 m量级。那么,激光源噪声耦合到星间位移测量噪声在mHz波段高达$\delta L = {\rm{0}}.1 \;{\rm{mm}}/\sqrt {{\rm{Hz}}} $。另一方面,基于激光干涉仪的引力波探测器的灵敏度为:

      $$h = \frac{{\delta L}}{L},$$ (2)

      结合上述位移测量精度,引力波探测器的灵敏度仅为$2 \times {10^{ - 17}}/\sqrt {{\rm{Hz}}} $,不能达到天基引力波探测所要求的${10^{ - 21}}/\sqrt {{\rm{Hz}}} $灵敏度目标。因此,除了在硬件上通过稳频腔和锁臂等技术路线[29-31]降低激光源频率噪声,还有必要针对天基迈克尔逊干涉仪臂长不匹配的情况,从算法上消除激光源噪声。

    • 时间延迟干涉技术采用数据后处理的方式构造出虚拟的等臂长干涉仪。例如,在图1中,激光先经过L2臂的往返,再经过L3臂的往返,相当于激光在虚拟干涉仪的一臂上传输。类似地,同一激光源产生的激光先经过L2臂的往返,再经过L3臂的往返,相当于激光在另一臂上的传输。这样,虚拟干涉仪的两臂臂长均为L2+L3,因此克服了卫星基线长度不相等的问题,实现了消除激光源频率噪声的目标。

      具体而言,在卫星1处,有两个探测器,分别记录本地激光与经过L2L3)往返的激光干涉后的相位测量值:

      $${\varphi _{21}}(t) = p(t - 2{L_2}) - p(t) + h_{21}^{{\rm{gw}}}(t),$$ (3)
      $${\varphi _{31}}(t) = p(t - 2{L_3}) - p(t) + h_{31}^{{\rm{gw}}}(t),$$ (4)

      其中,$p(t)$是激光源相位噪声,$p(t - 2{L_2})$$p(t - 2{L_3})$分别是激光经过L2臂和L3臂后在t时刻返回卫星1处的相位噪声。这里为了书写简便,省略了光速c。值得注意的是,本文对下标的命名规则沿用TDI系列文献的惯例。例如,卫星1与3之间的链路时延标为L2${\varphi _{21}}(t)$则表示本地卫星1处的激光与经过L2臂传输来的激光的干涉测量值。$h_{21}^{{\rm{gw}}}(t)$$h_{31}^{{\rm{gw}}}(t)$分别是L2臂和L3臂上的引力波导致的相位涨落,即引力波信号。为简单起见,忽略其它噪声项。两组相位测量值的相减,即${\varphi _{21}}(t) - {\varphi _{31}}(t)$,消去相同的本振光相位因子,实际上就是图1中迈克尔逊干涉仪中探测器的测量值:

      $${\varphi _{21}}(t) - {\varphi _{31}}(t) = p(t - 2{L_2}) - p(t - 2{L_3}) + h_{21}^{{\rm{gw}}}(t) - h_{31}^{{\rm{gw}}}(t),$$ (5)

      如前所述,L2L3最大相差2%,因此两臂上的激光相位噪声无法相消,参考目前的星载稳频激光源的噪声水平和天基引力波探测的灵敏度要求,(5)式中的相位噪声项将会淹没掉引力波信号。

      在TDI技术中,为了消去激光源相位噪声,对获得的测相值进行时间延迟处理,即在时间轴上进行平移。具体而言,对${\varphi _{21}}(t)$延迟2L3秒,对${\varphi _{31}}(t)$延迟2L2秒,有:

      $$\begin{split} {\varphi _{21}}(t - 2{L_3}) = &p(t - 2{L_2} - 2{L_3}) - p(t - 2{L_3}) +\\ & h_{21}^{{\rm{gw}}}(t - 2{L_3}), \end{split} $$ (6)
      $$\begin{split} {\varphi _{31}}(t - 2{L_2}) =& p(t - 2{L_2} - 2{L_3}) - p(t - 2{L_2})+ \\ & h_{31}^{{\rm{gw}}}(t - 2{L_2}), \end{split} $$ (7)

      对上述两组延时处理的测相值作差有:

      $$\begin{split} {\varphi _{21}}(t - 2{L_3}) -& {\varphi _{31}}(t - 2{L_2}) = p(t - 2{L_2}) - p(t - 2{L_3}) +\\ & h_{21}^{{\rm{gw}}}(t - 2{L_3}) - h_{31}^{{\rm{gw}}}(t - 2{L_2}), \end{split} $$ (8)

      进一步地,比较(5)式和(8)式,均包含相同的激光相位噪声项。因此,通过(5)式与(8)式相减,可以得到一组含有引力波信号而不含有激光源噪声的数据组合,即:

      $$X \equiv \left[ {{\varphi _{21}}(t) - {\varphi _{31}}(t)} \right] - \left[ {{\varphi _{21}}(t - 2{L_3}) - {\varphi _{31}}(t - 2{L_2})} \right],$$ (9)

      为了更好地理解TDI数据组合背后的物理思想,将(9)式重新写作:

      $$X = \left[ {{\varphi _{21}}(t) + {\varphi _{31}}(t - 2{L_2})} \right] - \left[ {{\varphi _{31}}(t) + {\varphi _{21}}(t - 2{L_3})} \right],$$ (10)

      (10)式中$\left[ {{\varphi _{21}}(t) + {\varphi _{31}}(t - 2{L_2})} \right]$表示一束激光从卫星1处出发,先经过L2臂的往返后与本地激光进行干涉,随后再经过L3臂的往返后与本地激光进行干涉,如图2左图所示。其中,${\varphi _{31}}(t - 2{L_2})$指经过2L2秒延时后的相位测量值。类似地,$\left[ {{\varphi _{31}}(t) + {\varphi _{21}}(t - 2{L_3})} \right]$表示一束激光从卫星1处出发,先经过L3臂的往返后与本地激光进行干涉,随后再经过L2臂的往返后与本地激光进行干涉,如图2右图所示。这两种情况对应激光均经过了2L2+2L3的传输,相当于构造出臂长均为2L2+2L3的迈克尔逊干涉仪,从而实现了激光源相位噪声的消除。

      图  2  TDI原理示意图

      Figure 2.  Diagram of the principles behind the TDI

      可以看到,TDI技术实际上是对星间相位测量值进行适当的延时处理,再进行数据组合,从而构造出臂长相等的干涉仪,从而消除激光源相位噪声。

    • 为了具体给出TDI数据组合方式,以及方便下一节讨论消除时钟噪声的原理,这里以LISA项目为例,简要讨论天基引力波探测的光学系统布局。如图3所示,LISA星座中每一颗卫星均包含有两个光学平台。如图3插图所示,在卫星i处,当面对三星编队中心时,左手处的光学平台定义为i号平台,右手处的光学平台定义为i*号平台。卫星i与卫星j之间的距离定义为Lk。在卫星i处,定义本地激光与远星处经过Lj传输的激光干涉后的相位测量值为${\varphi _{ji}}$。另外,在同一光学平台上,还存在本地激光与相邻光学平台之间的激光干涉测量${\beta _{ji}}$。这样的下标命名规则是为了在得到一个卫星上的所有物理量表达式后,通过下标变换的方式,将表达式中下标“1”变为“2”,“2”变为“3”,“3”变为“1”,即“$1 \to 2 \to 3 \to 1$”式的下标循环变换,就能够直接得到其它卫星上相应物理量的表达式。

      图  3  LISA光学系统示意图(已获参考文献[19]授权© American Physical Society)

      Figure 3.  The illustration of optical system in the LISA mission (Reprinted with permission from ref. [19] © American Physical Society)

      图3给出了卫星1的两个光学平台上的四组测相值(${\varphi _{31}}$$\;{\beta _{31}}$${\varphi _{21}}$$\;{\beta _{21}}$)。${\varphi _{31}}$$\;{\beta _{31}}$是卫星1的1号光学平台上输出的测相值。${\varphi _{31}}$对应的干涉测量为:卫星2的2*平台上发射过来一束激光(相位噪声为$p_2^*$),经过L3星间链路,到达卫星1的1号平台,在被该平台上悬浮的检验质量块反射后,与本地激光${p_1}$进行干涉后输出测相值${\varphi _{31}}$$\;{\beta _{31}}$对应的干涉测量为:右侧的1*号平台上的激光$p_1^ * $经过本平台上的检验质量块的镜面反射后,通过光纤链路到达左侧的1号平台,与该平台上的激光${p_1}$发生干涉,输出测相值$\;{\beta _{31}}$。同理,${\varphi _{21}}$是卫星3上的激光经过星间链路L2后与卫星1的1*号平台激光的干涉测量结果,而$\;{\beta _{21}}$则是1*号平台上输出的相邻平台间激光干涉测相值。

      对于$ {\varphi _{31}}$,有:

      $${\varphi _{31}}(t) = p_2^*(t - {L_3}) - {p_1}(t) + h_{31}^{{\rm{gw}}}(t),$$ (11)

      $h_{31}^{{\rm{gw}}}(t)$为星间引力波引起的相位涨落。这里忽略了其它噪声项,它们主要包括光路指向误差,弱光散粒噪声,光学平台位移噪声和检验质量块加速度噪声。为了方便后续推导,定义时间延迟算符:${D_i}x(t) = x(t - {L_i})$。施加两次延迟算符,等于对该物理量做了两次时间延迟,即${D_j}{D_i}x(t) = x(t - {L_i} - {L_j})$。另外,当链路延迟Li不随时间变化时,时间延迟算符具有对易性,即${D_j}{D_i}x(t) = {D_i}{D_j}x(t)$。引入时间延迟算符后,(11)式可写作:

      $${\varphi _{31}} = {D_3}p_2^* - {p_1} + h_{31}^{{\rm{gw}}},$$ (12)

      对于卫星1内部相邻光学平台间的干涉测量${\;\beta _{31}}$,有:

      $${\beta _{31}} = p_1^* - {p_1},$$ (13)

      这里相邻光学平台之间的链路延迟较小,不考虑。同时,忽略了光学平台噪声和检验质量块噪声,以及光纤链路噪声。

      同理,对于${\varphi _{21}}$而言,有:

      $${\varphi _{21}} = {D_2}{p_3} - p_1^* + h_{21}^{{\rm{gw}}},$$ (14)

      对于$\;{\beta _{21}}$,有:

      $${\beta _{21}} = {p_1} - p_1^*,$$ (15)

      在得到了卫星1的四组相位测量的表达式后,通过“$1 \to 2 \to 3 \to 1$”下标循环变换,就可以得到其它两颗卫星上的八组测相表达式。

    • TDI数据组合的本质是通过对六组测相值${\varphi _{ij}}$进行链路延迟和线性组合,从而消除三星上六台激光源的相位噪声项${p_i}$$p_i^*$,保留引力波信号$h_{ij}^{{\rm{gw}}}$。本小节简要描述TDI数据组合的推导原理。

      同时,由于数据组合的目标是消除激光相位噪声项,因此忽略引力波信号和其它噪声项,六组测相值${\varphi _{ij}}$可以写作:

      $$\begin{split} &{\varphi _{31}} = {D_3}p_2^* - {p_1}, \\ &{\varphi _{12}} = {D_1}p_3^* - {p_2}, \\ &{\varphi _{23}} = {D_2}p_1^* - {p_3}, \end{split} $$ (16)
      $$\begin{split} & {\varphi _{21}} = {D_2}{p_3} - p_1^*, \\ & {\varphi _{32}} = {D_3}{p_1} - p_2^*, \\ &{\varphi _{13}} = {D_1}{p_2} - p_3^*, \end{split} $$ (17)

      (16)式描述三颗卫星左侧光学平台的测相值,(17)式描述右侧光学平台的测相值。因此,TDI数据组合的目标便是给出六个系数,使得$\displaystyle\sum\limits_{i,j} {{x_{ij}}{\varphi _{ij}}} = 0$,其中系数${x_{ij}}$包含时间延迟算符$D_1^{{k_1}}D_2^{{k_1}}D_3^{{k_2}}$,使得测相数据在时间上延迟${k_1}{L_1} + {k_2}{L_2} + {k_3}{L_3}$

      考虑一种简单的情况,同一颗卫星相邻光学平台上的激光源具有相同的相位噪声,即${p_i} = p_i^*$。对于$\displaystyle\sum\limits_{i,j} {{x_{ij}}{\varphi _{ij}}} = 0$,存在两类最简单的解,它们包含的时间延迟算符最少。第一类解称为萨格纳克(Sagnac)解,按惯例也称为$\zeta $型数据组合:

      $$\zeta = {\varphi _{31,1}} - {\varphi _{12,2}} - {\varphi _{23,3}} + {\varphi _{21,1}} + {\varphi _{32,2}} + {\varphi _{13,3}},$$ (18)

      这里,为了简便起见将延迟算符写入下标,例如将${D_1}{\varphi _{31}}$项写作${\varphi _{31,1}}$

      第二类解称为$\alpha $型数据组合:

      $$\alpha = - {\varphi _{31}} - {\varphi _{12,3}} - {\varphi _{23,13}} + {\varphi _{21}} + {\varphi _{32,12}} + {\varphi _{13,2}},$$ (19)

      通过“$1 \to 2 \to 3 \to 1$”下标变换,还可以得到两个新的数据组合,它们实质是等价的。

      上述两类解累计施加了6次时间延迟算符。除了这两类最简单的解,还存在一类常见的解,称为不等臂迈克尔逊数据组合,习惯称为X型数据组合:

      $$X = {\varphi _{21}} + {\varphi _{23,2}} + {\varphi _{31,22}} + {\varphi _{32,322}} - {\varphi _{31}} - {\varphi _{32,3}} - {\varphi _{21,33}} - {\varphi _{23,233}},$$ (20)

      同样地,对X型数据组合进行下标循环变换,可以得到等价的Y型和Z型数据组合。

      虽然对比前面的$\zeta $型和$\alpha $型数据组合,X型数据组合累积施加了12次时间延迟算符,然而,X型数据中不包含卫星2和卫星3间的干涉测量值,即${\varphi _{12}}$${\varphi _{13}}$,也不需要施加这两颗卫星间的链路延迟算符D1。因此,三星编队中某一条星间激光链路建立失效的情况下,X型数据组合仍然可以正常工作。因此,相比$\zeta $型和$\alpha $型数据组合。X型数据组合提升了系统冗余性,并且减少了对星间链路L1的测距需求。实际上,X型数据组合背后的物理模型与上一节中的图2所示是类似的。

      具体而言,(20)式中X型数据组合的前四项对应图4左图,即激光沿①-②-③-④路径传输的过程。例如,${\varphi _{21}}$对应激光从卫星3出发,经过星间链路L2后,到达卫星1与本地激光干涉。而第二项${\varphi _{23,2}}$对应经过L2秒时间延迟后的干涉测量,即${\varphi _{23,2}}$是紧接着${\varphi _{21}}$进行的。以此类推,${\varphi _{31,22}}$${\varphi _{32,322}}$分别对应接下来的图4左图的③和④。类似地,(20)式后4项对应图4右图的传输路径。可以看到,X型数据组合对应两臂臂长均为2L2+2L3的迈克尔逊干涉仪测量结果,因此激光源的相位噪声被消除。

      图  4  X型数据组合物理模型

      Figure 4.  The physics modal of the X type data combination

      值得注意的是,在上述分析中,存在相邻光学平台间的激光源噪声相等的假设,即${p_i} = p_i^*$。对于${p_i} \ne p_i^*$的情况,利用相邻光学平台间的干涉测量值${\beta _{ij}}$,再通过线性代换就可以消去$p_i^*$项。例如,通过${\varphi _{31}} - {D_3}{\beta _{12}}$就可以将${\varphi _{31}}$${D_3}p_2^* - {p_1}$变换为${D_3}{p_2} - {p_1}$,而待求解的数学方程形式仍然不变。因此,数据组合的类型仍然与前面的分析一致。以X型数据组合为例,考虑${p_i} \ne p_i^*$,同时考虑消除各个光学平台的位移噪声[18, 19],得到的数据组合形式为:

      $$\begin{split} X =& {\varphi _{21}} + {\varphi _{23,2}} + {\varphi _{31,22}} + {\varphi _{32,322}} - {\varphi _{31}} - {\varphi _{32,3}} -\\& {\varphi _{21,33}} - {\varphi _{23,233}} + \dfrac{1}{2}( - {\beta _{21,2233}} + {\beta _{21,33}} + {\beta _{21,22}} - {\beta _{21}} +\\& {\beta _{31,2233}} - {\beta _{31,33}} - {\beta _{31,22}} + {\beta _{31}}) \end{split} $$ (21)

      另外,相比$\zeta $型和$\alpha $型数据组合,考虑${p_i} \ne p_i^*$和消除光学平台噪声后,X型数据组合仅增加了卫星1上的测相值${\beta _{21}}$${\beta _{31}}$,而$\zeta $型和$\alpha $型数据组合需要增加全部的六组相邻平台测相值,这也是X型数据组合的优势。

      可以看到,TDI通过对各个卫星平台上的单程测相值进行时间延迟和线性组合,构造出等臂长的迈克尔逊型或萨格纳克型干涉仪,从而消除激光源频率噪声。

    • 前面关于TDI的讨论主要基于星座和臂长均处于静态的情况。实际上需要考虑星座旋转和臂长变化的情况,从而形成第二代TDI数据组合方式。

      当考虑星座旋转时,任意两颗卫星间的双向激光链路长度是不相等的。如图5所示,当星座顺时针刚性旋转时,1星与2星间的实际星间距离${l_2}$保持不变。然而,从1星到2星的激光链路${L_2}$会略小于${l_2}$,而从3星到1星的激光链路${L'_2}$会略大于${l_2}$。因此,需要对不同方向的激光链路进行区分。例如,将顺时针传输的激光链路标记为${L'_i}$,将逆时针传输的激光链路标记为${L_i}$。这样,对于如(20)式的X型数据组合,需要结合图4所示的具体传输路径进行修正,以区分不同方向的激光链路。进一步地,考虑星间距离变化时,时间延迟算符具有不对易性,即$\left[ {{D_i},{D_j}} \right] \ne 0$。因此,直接应用第一代TDI数据组合时,会存在若干激光噪声项无法相互抵消的情况。例如,在第一代X型数据组合中,会残余${D_{2'}}{D_2}{D_3}{D_{3'}}{p_1}$${D_3}{D_{3'}}{D_{2'}}{D_2}{p_1}$两项无法相互抵消。如果星间相对运动速度之差,即${\dot L_2}$${\dot L_3}$之差,在10 m/s量级,那么需要将激光频率噪声压低大约一个数量级,才能抵消由于延迟算符不对易性引入的激光残余噪声[32]。在第二代TDI数据组合中,对于${D_i}p(t) = p(t - {L_i}(t))$,仅考虑一阶线性项,忽略高阶项,可近似有${D_i}p(t) \approx p(t) + \dot p(t){L_i}(t)$。当施加两次时间延迟算符时,有${D_j}{D_i}p(t) \!=\! p(t \!-\! {L_j}(t) \!-\! {L_i}(t \!-\! {L_j}(t)))$。其中,${L_i}(t - {L_j}(t))$项可近似为${L_i}(t) - {\dot L_i}(t)L$,而L取典型的星间距离常量,如LISA中的16.67光秒。另外,${\dot L'_i}(t) \approx {\dot L_i}(t)$,即近似认为从j星到k星的星间距离变化率约等于从k星到j星的距离变化率。不难得出,对于第一代TDI数据组合中的$\left[ {{D_{2'}}{D_2},{D_3}{D_{3'}}} \right]{p_1}$项,会有激光相位噪声残差${\rm{4}}L({\dot L_3}(t) - {\dot L_2}(t)){\dot p_1}$。然而,在前述一阶近似的考虑下,可以得到$\left[ {D_{2'}}{D_2}{D_3}{D_{3'}}, {D_3}{D_{3'}}{D_{2'}}{D_2} \right]{p_1} \approx {\rm{0}}$。第二代X型TDI数据组合正是利用了这两个延迟算符项的对易性实现了激光相位噪声的消除[17]

      图  5  当星座旋转时,${L'_i} \ne {L_i}$

      Figure 5.  In the case of the constellation rotation, ${L'_i} \ne {L_i}$

    • 在相位测量系统中,时钟频率(相位)抖动会引入到测相误差中。因此,需要评估在甚高测相精度要求下,时钟频率稳定度需要达到的水准。

      天基引力波探测的位移测量精度为十皮米量级,分配到探测器后端的相位读出部分的噪声预算为$\varphi (f) = {10^{ - 6}}{\rm{cycles/}}\sqrt {{\rm{Hz}}} $,关心的频点为$f = 3\;{\rm{ MHz}}$。光电探测器拍频信号取星间最大的多普勒频移${f_0} = 20\;{\rm{ MHz}}$。因此,时钟频率稳定度的要求为[22]

      $$y(f) = 2\pi f\varphi (f)/{f_0},$$ (22)

      $y(f)$需要优于$9 \times {10^{ - 16}}/\sqrt {{\rm{Hz}}} $。目前的星载超稳晶振(Ultra-stable Oscillator,USO)的稳定度在${10^{ - 13}}/\sqrt {{\rm{Hz}}} $水平。因此,从引力波探测的灵敏度目标来看,消除USO噪声变得非常关键。只关注激光源噪声以及光学平台噪声的消除,不考虑USO噪声的抑制问题,会使得引力波信号隐没在USO噪声中,从而导致探测任务的失败。然而,从硬件上将USO稳定度提升3个数量级是非常困难的,因此,加州理工大学M. Tinto等人提出通过时钟边带调制的方法消除时钟噪声[19, 20]。目前,相关论证工作仍在进行中。

    • 在时钟边带调制方案中,时钟信号倍频到GHz,随后通过电光调制器加载到激光相位上。相位调制深度较小,因此调制到边带上的激光功率较低,约占总功率的10%。这里以卫星1和卫星3的干涉测量为例,如图6所示。卫星1的1*号光学平台上激光源中心频率为$\nu _1^*$,卫星3的3号光学平台上激光源中心频率为${\nu _3}$。两颗卫星上的激光源相位分别被GHz的时钟信号${f_1}$${f_3}$调制,分别在载波中心频率附近产生$\nu _1^* \pm {f_1}$${\nu _3} \pm {f_3}$的一阶边带。因此,在光电探测端,主要存在两个拍频信号:(1)载波之间的拍频信号,频率为${\nu _3} - \nu _1^*$,输出相位测量值${\varphi _{21}}$;(2)边带之间的拍频信号$({\nu _3} - \nu _1^*) \pm ({f_3} - {f_1})$,输出相位测量值${\varphi '_{21}}$。为了从载波拍频信号中区分出边带拍频信号,两颗卫星的GHz时钟信号通过要有MHz量级的差别。例如${f_3} - {f_1} = 5\;{\rm{MHz}}$,这样,在主要的载波拍频峰两侧的5 MHz处,会有边带拍频峰。通过后端滤波就可以分别处理两个信号。

      图  6  时钟边带调制,EOM:电光调制器

      Figure 6.  Clock sideband modulation, EOM: electro-optic modulator

      时钟边带调制消除时钟噪声的基本原理如下:载波拍频信号对应的测相值${\varphi _{21}}$中包含本地卫星1的时钟噪声,而加载时钟相位的边带拍频信号${\varphi '_{21}}$包含卫星1和3的时钟噪声。另一方面,${\varphi _{21}}$${\varphi '_{21}}$中的其他项,如引力波信号和激光源噪声项均相等。因此,对两组测相值作差就能提取出两颗卫星上的时钟噪声项。进一步地,在TDI数据组合中减去时钟噪声项,抑制时钟噪声问题对引力波探测的影响。

    • 本小节以X型数据组合为例,具体阐述时钟噪声消除的原理。图7图6探测端的信号处理示意图。这里进一步考虑了星间多普勒频移。

      图  7  探测端信号处理示意图,NCO:数控振荡器

      Figure 7.  TSignal processing illustration of photon detection, NCO: Numerically Controlled Oscillator

      图7所示,${\nu _3}(1 - {\dot L_2})$是3星激光到达1星处的载波中心频率。这里为了书写简便,省略了光速c。星间多普勒频移$\left| {{\nu _3}{{\dot L}_2}} \right| \leqslant 20\;{\rm{MHz}}$$({\nu _3} + {f_3}) (1 - {\dot L_2})$是3星激光到达1星处的一阶边带信号。探测器分别输出如图7所示的两个拍频信号,随后与数字信号处理单元的数控振荡器(Numerically Controlled Oscillator,NCO)混频,分别输出测相值${\varphi _{21}}$${\varphi '_{21}}$

      具体地,载波相位输出值${\varphi _{21}}$的表达式为:

      $${\varphi _{21}} = [{\nu _3}(1 - {\dot L_2}) - \nu _1^* - {a_{21}}{f_1}]t + {p_{3,2}} - p_1^* - {a_{21}}{q_1} + h_{{\rm{21}}}^{{\rm{gw}}},$$ (23)

      其中,${f_1}$是卫星1处倍频到GHz的时钟频率,${a_{21}}$是NCO输出系数,${p_{3,2}} = {D_2}{p_3} = {p_3}(t - {L_2})$是卫星3处的激光相位噪声,$p_1^*$是卫星1激光相位噪声,${a_{21}}{q_1}$是卫星1处的时钟通过NCO引入到相位测量中的噪声项。对于NCO系数${a_{21}}$,有:

      $${a_{21}} = \frac{{{\nu _3}(1 - {{\dot L}_2}) - \nu _1^*}}{{{f_1}}},$$ (24)

      类似地,时钟边带拍频信号${\varphi '_{21}}$的表达式为:

      $$\begin{split} {{\varphi '}_{21}} =& [({\nu _3} + {f_3})(1 - {{\dot L}_2}) - \nu _1^* - {f_1} - {b_{21}}{f_1}]t + {p_{3,2}} + {q_{3,2}} -\\ & p_1^* - {q_1} - {b_{21}}{q_1} + h_{{\rm{21}}}^{{\rm{gw}}}, \end{split} $$ (25)

      相比${\varphi _{21}}$${\varphi '_{21}}$表达式增加了时钟噪声项${q_{3,2}}$${q_1}$项,这正是因为两地激光源均被时钟信号进行了相位调制,时钟噪声因而混入了边带信号中。另外,与边带拍频信号混频的NCO系数${b_{21}}$为:

      $${b_{21}} = \frac{{({\nu _3} + {f_3})(1 - {{\dot L}_2}) - \nu _1^* - {f_1}}}{{{f_1}}},$$ (26)

      显而易见,由于${\varphi _{21}}$${\varphi '_{21}}$之间的区别仅在于时钟边带调制。因此,除了时钟噪声项外,${\varphi _{21}}$${\varphi '_{21}}$有着相同的激光源噪声项和引力波信号项,以及其它未写出的噪声项。通过对两组测相值作差,就能够提取出时钟噪声项:

      $${\varphi '_{21}} - {\varphi _{21}} = {q_{3,2}} - (1 - {\dot L_2}){f_3}\dfrac{{{q_1}}}{{{f_1}}},$$ (27)

      星间相对运动速度不超过20 m/s,远远小于光速,因此$(1 - {\dot L_2}) \approx 1$。定义变量${r_{21}} \equiv \dfrac{{{\varphi _{21}} - {{\varphi '}_{21}}}}{{{f_3}}}$,有:

      $${r_{21}} = \dfrac{{{q_1}}}{{{f_1}}} - \dfrac{{{q_{3,2}}}}{{{f_3}}},$$ (28)

      类似地,对于卫星1左侧的1号光学平台${\varphi _{31}}$${\varphi '_{31}}$,引入变量${r_{31}} \equiv \dfrac{{{\varphi _{31}} - {{\varphi '}_{31}}}}{{{f_2}}}$,同样有:

      $${r_{31}} = \frac{{{q_1}}}{{{f_1}}} - \frac{{{q_{2,3}}}}{{{f_2}}},$$ (29)

      通过下标循环变化“$1 \to 2 \to 3 \to 1$”,就可以得到其它两颗卫星上的四组表达式。

      对于(21)式所示的X型数据组合,当考虑如(23)式所示的时钟噪声项${q_i}$后,包含时钟噪声的Xq表达式为:

      $$\begin{split} {X_q} =& {X^{{\rm{gw}}}} + {a_{32}}[{q_{2,3}} - {q_{2,322}}] - {a_{23}}[{q_{3,2}} - {q_{3,233}}] + {a_{31}}[{q_1} -\\ & {q_{1,22}}] - {a_{21}}[{q_1} - {q_{1,33}}] + {c_{21}}[({q_1} - {q_{1,22}}) - {({q_1} - {q_{1,22}})_{,33}}], \end{split} $$ (30)

      其中,${c_{21}}$为卫星1的1*号平台处${\beta _{21}}$探测器(图3)后端NCO的输出系数:

      $${c_{21}} = - {c_{12}} = \frac{{{\nu _1} - v_1^*}}{{{f_1}}},$$ (31)

      通过${r_{ij}}$可以得到Xq式中$[{q_{2,3}} - {q_{2,322}}]$等项表达式:

      $$\begin{split} &[{q_1} - {q_{1,22}}] = {f_1}[{r_{21}} + {r_{23,2}}], \\ &[{q_1} - {q_{1,33}}] = {f_1}[{r_{31}} + {r_{32,3}}], \\ &[{q_{3,2}} - {q_{3,233}}] = {f_3}[{r_{21,33}} - {r_{21}} + {r_{31}} + {r_{32,3}}], \\ &[{q_{2,3}} - {q_{2,322}}] = {f_2}[{r_{31,22}} - {r_{31}} + {r_{21}} + {r_{23,2}}], \end{split} $$ (32)

      得到上述关系式后,就可以减去时钟噪声项,得到修正时钟噪声项的X型数据组合:

      $$\begin{split} X \equiv &{X_q} - {a_{32}}{f_2}[{r_{31,22}} - {r_{31}} + {r_{21}} + {r_{23,2}}] +\\ &{a_{23}}{f_3}[{r_{21,33}} - {r_{21}} + {r_{31}} + {r_{32,3}}] - {a_{31}}{f_1}[{r_{21}} + {r_{23,2}}] +\\ &{a_{21}}{f_1}[{r_{31}} + {r_{32,3}}] + {c_{21}}{f_1}[{r_{21,33}} + {r_{23,233}} - {r_{21}} - {r_{23,2}}], \end{split} $$ (33)

      可以看到,时钟调制产生了边带拍频信号,通过与载波拍频信号作差得到${r_{ij}}$,提取出时钟噪声因子,从而能够减去TDI数据组合中的时钟噪声项,最终实现消除时钟噪声的目的。值得注意的是,上述时钟噪声消除的讨论仅仅针对静态的第一代TDI数据处理方式。由于时间延迟算符的非对易性引入的时钟相位噪声可以忽略不计,因此在第二代TDI数据组合中,可以采用类似的处理方法,扣除相应时钟相位噪声项[20]

    • 在天基引力波探测中,为了实现${10^{ - 6}}$周的甚高测相精度,除了时钟相位噪声问题,还需要重点考虑模数转换(Analogue to Digital Conversion, ADC)过程引入的相位噪声。ADC噪声主要有量化噪声和采样时间抖动噪声。其中,量化噪声可以忽略不计,而采样时间抖动噪声由ADC器件本身的孔径抖动和外部的采样时钟抖动构成。为了消除采样时间抖动噪声,一种方案是从时钟处引出一路导频(Pilot Tone)信号,令导频信号经过ADC器件[33, 34]。导频信号通过数字锁相环路处理后读出相应的相位增量。假设导频信号的频率${f_{{\rm{pilot}}}}$非常稳定,那么理想的相位增量${f_{{\rm{pilot}}}}\Delta T$是恒定的。这里,$\Delta T$是相位读取周期。然而,测量得到的导频信号相位增量会由于ADC采样时间抖动而在${f_{{\rm{pilot}}}}\Delta T$附近抖动,相应的偏差$\delta {\varphi _{\det }}$即为该ADC器件由于采样时间抖动引入的测相噪声。这样,对于同样经过该ADC器件采样处理的拍频信号${f_{{\rm{beat}}}}$,添加相位纠正项${{\delta {\varphi _{\det }}{f_{{\rm{beat}}}}} / {{f_{{\rm{pilot}}}}}}$,即可消除ADC采样时间抖动引入的测相噪声。然而,导频方案的问题在于,USO有限的频率稳定度会导致${f_{{\rm{pilot}}}}\Delta T$本身不是恒定值,因此会导致$\delta {\varphi _{\det }}$不仅包含ADC采样时间抖动因素,还包含了时钟相位噪声${q_{{\rm{pilot}}}}$。由于时钟相位噪声与时钟信号的频率成正比,因此可以将${q_{{\rm{pilot}}}}$添加到前述测相表达式中。例如,在3.3小节中,1星上频率为${f_{\rm{1}}}$的时钟信号的相位噪声为${q_1}$,因此通过导频纠正方案引入的时钟相位噪声项${q_{{\rm{pilot}}}}{\rm{ = }}{{{q_1}{f_{{\rm{pilot}}}}} / {{f_{\rm{1}}}}}$。将其添加到测相表达式中,随后再进行相应的时钟相位噪声消除处理。可以看到,用来抑制模数转换噪声的导频方案需要建立在时钟噪声消除的基础上来实现,凸显了时钟噪声消除方案的重要性。

    • 在TDI的数据处理和重组中,为了消除激光源和时钟频率(相位)噪声,关键在于对相位测量值进行精确的时间延迟处理,这就要求对天基干涉仪的臂长,即星间距离有精确的认识。从这个角度看,TDI技术抑制激光源和时钟噪声的代价是需要对星间绝对距离进行一定精度的测量。

      当对星间绝度距离的测量精度为$\Delta L$时,那么,通过TDI数据处理构造出的干涉仪臂长差就是$\Delta L$。因此,激光源频率噪声耦合到星间相对距离的测量误差$\delta L$为:

      $$\delta L = \Delta L\frac{{\delta \nu }}{{{\nu _0}}},$$ (34)

      在mHz波段,假如激光源频率噪声$\delta \nu = 300\;{\rm{ Hz}}/\sqrt {{\rm{Hz}}} $,激光源中心频率${\nu _0} \approx 300\;{\rm{ THz}}$,激光源噪声耦合到位移测量噪声的误差分配在$1\;{\rm{ pm}}/\sqrt {{\rm{Hz}}} $时,因此星间绝对距离测量精度$\Delta L$应达到1 m量级。

      值得注意的是,(1)式与(34)式虽然形式相等,但$\Delta L$的物理含义不同。(1)式中的$\Delta L$是利用三星编队直接形成的迈克尔逊干涉仪的臂长差,高达108 m。(34)式是经过TDI算法构建的干涉仪臂长差。实际上(1)式与(34)式的$\Delta L$之比就是TDI对激光源噪声的抑制倍数,即达到了108倍。

    • 在LISA项目中,采用直接序列扩频进行数据传输和测距编码。数据码采用BPSK编码,伪随机码进行扩频后进行绝对距离测量。如图8所示,数据码和伪随机码构成的扩频码与时钟信号一同加载到载波激光相位上。类似于时钟信号的浅相位调制,扩频码的相位调制系数也很低,将消耗载波功率的1%左右。

      图  8  在载波相位上调制数据码和伪随机码用于星间数据传输和绝对距离测量,EOM:电光调制器

      Figure 8.  The phase of the carrier is modulated by data codes and pseudo-random codes for the inter-satellite data transmission and absolute distance measurement, EOM: electro-optic modulator

      图9是LISA项目扩频码的参数设计示意图。上层的数据码调制速率为24 kbps,对应的码周期为40.96 μs。底层的伪随机码调制速率为1.5 Mbps,码片周期为640 ns。伪随机码序列的重复频率为1.5 kHz,每段序列包含1024个伪随机码码片。因此,伪随机码序列覆盖的测距长度为200 km,测距信息刷新频率为1.5 kHz。对于LISA星座500万公里的星间距离,需要一个精度优于200 km量级的初始定位系统。深空网络结合每颗卫星上的恒星追踪器可以提供精度约为25 km的星间距离粗测量。经过初始定位后,就可以利用伪随机码来提供更精准的距离测量。

      图  9  数据码和伪随机码参数设计(已获参考文献[24]授权© The Optical Society)

      Figure 9.  Parameter design of data codes and pseudo-random codes (Adapted with permission from ref. [24] © The Optical Society)

      德国阿尔伯特-爱因斯坦研究所(AEI)针对时钟噪声转移和伪码测距等技术进行了一系列地面论证实验。图10是在光载波上调制GHz时钟频率和MHz扩频码的干涉测量结果。扩频码调制深度为0.1 rad,时钟边带调制深度为0.45 rad,分别约消耗光载波功率的1%和10%。

      图  10  光载波经过时钟信号和扩频调制后干涉测量结果(已获参考文献[24]授权© The Optical Society)

      Figure 10.  Interference measurement results when the optical carriers are modulated by the clock signal and spread spectrum (Reprinted with permission from ref. [24] © The Optical Society)

      图10左图是强光条件下的测试结果,中间的信号峰为载波拍频信号,附近则是扩频频谱。由于两台激光源的时钟调制频率相差5 MHz,因此在距离载波拍频峰±5 MHz处会出现时钟边带的拍频信号。图10右图是模拟天基引力波环境下的100 pW弱光探测结果,对应的扩频码所占功率仅为1 pW。在散粒噪声较大的情况下,模拟星上3 Hz伪码延时锁定环测量带宽,系统的测距均方根误差为42 cm,优于1 m的LISA任务绝对距离测量要求。另外,数据传输的误码率为26×10−3,应用前向纠错编码可使误码率进一步降低到可实际应用的10−9水平。

      值得注意的是,目前的伪距测量论证实验是在静态条件下完成的。在实际的在轨测距中,由于星间链路延时达数十秒,而星间相对运动速度在10 m/s量级,因此需要考虑在激光传输过程中卫星相对运动造成的测距误差。

    • 加州理工大学喷气推进实验室(JPL)针对TDI技术的激光源噪声和时钟噪声消除进行了地面验证实验。图11是在0.2 mHz到3 Hz频段几种测量方式下位移测量噪声线性谱密度的对比情况[22]。由于激光源中心波长为1064 nm,因此1 μcycle≈1 pm。

      图  11  LISA干涉仪测试床位移测量结果(已获参考文献[22]授权© American Physical Society)

      Figure 11.  The displacement measurment results in the LISA interferometry test bed (Reprinted with permission from ref. [22] © American Physical Society)

      图11曲线i)是原始的测量结果,由激光源相位噪声占主导。曲线ii)对应于采用α型数据组合后消除了激光源噪声的情况。然而此时受到曲线iii)所示的时钟相位噪声的限制,测量噪声仍然比较大。在典型的3 mHz波段,位移噪声和时钟噪声曲线是重合的。可以看到,尽管TDI数据处理令位移测量噪声下降了约5个数量级,然而由于在数据组合中未考虑时钟噪声项的修正问题,时钟噪声已经成为位移测量噪声的门限。实际上通过对比(1)式和(24)式可以发现,激光源噪声(${\rm{1}}{{\rm{0}}^{{\rm{ - 12}}}}/\sqrt {{\rm{Hz}}} $)和时钟噪声(${10^{ - 13}}/\sqrt {{\rm{Hz}}} $)耦合到星间位移测量噪声中的误差相差了约5个数量级。因此,消除了占主导作用的激光源噪声,还需要考虑消除次级的时钟噪声项。进一步地,通过第3节所示的时钟边带调制方案,消除α型数据组合中的时钟噪声项,得到了曲线iv)所示的位移噪声曲线,已经基本达到曲线v)所示的干涉仪极限水平,满足曲线vi)所示的LISA测距精度要求。可以看到,在3 mHz处,经过时钟噪声纠正后的TDI技术抑制了激光频率噪声109倍,抑制了时钟噪声5.8×104倍。

      与JPL实验室放置在真空腔内的短臂长干涉仪不同,科罗拉多大学博尔德分校的S. E. Pollack等人[35, 36]和法国巴黎狄德罗大学的P. Prat等人[37, 38]分别开发了小型的桌上电光系统,用来模拟LISA星座中的大范围光学传播时延和星间多普勒频移,以验证相位计与探测器等硬件部分和TDI算法的可靠性。S. E. Pollack等人搭建了一套LISA地面测试平台[35],用来生成接近实际在轨情况的激光干涉条纹,例如拍频信号频率从50 kHz变化到15 MHz,并加入了星间多普勒频移的变化速率以及加速度。干涉信号由100 pW量级的弱光与0.5 mW激光干涉产生,模拟了星间远距离传输来的激光与本地激光的干涉测量。实验论证了过零点计数相位计的有效性,验证了TDI算法对模拟产生的激光相位噪声的有效消除。进一步地,他们在测试平台上通过相位浅调制在激光链路上附加了数据传输和时钟边带调制功能,验证了包含科学测量信号的主拍频信号与链路附加信号的互不干扰[36]。P. Prat等人搭建了了一套光电模拟系统[37, 38]。该模拟系统具有光学模拟信道和电子模拟信道,产生激光拍频或射频间的混频信号,并最终进入相位计模块进行测相和数据处理。由于采用TDI中的X型数据组合方式,因此只需要模拟两条激光链路上的干涉测量即可。在电学模拟实验中,由直接数字信号合成器(DDS)产生3个100 MHz左右的射频信号。通过射频信号之间的混频来模拟星间激光干涉拍频信号。通过设置射频信号的相位噪声和时间延迟来模拟LISA中的噪声和链路时延。通过TDI技术使得信号噪声降低了109倍。在光学实验系统中,采样声光调制器(AOM)对激光链路进行移频的方式模拟星间多普勒频移。通过TDI技术使得光学信号的噪声降低了5×107倍。光学信号的降噪倍数较低,原因在于激光传输路径上存在光程抖动因素。另外,APC实验室还在相位计的ADC端口添加了如3.4小节所述的导频纠正方案。在后续的工作中APC实验室将考虑链路时延和多普勒频移随时间变化的情况,以使得光电模拟信道更加逼近真实的在轨情况。

    • 基于天基激光干涉仪的引力波探测项目由于高灵敏度探测要求,将极大推动星载激光收发终端和星间激光链路指向控制等研究方向的技术进步,并孵化出了面向高精度地球重力场测量的激光干涉测距载荷。针对十皮米量级的高精度星间位移测量要求,需要首先解决激光源频率噪声和时钟频率噪声的问题。时间延迟干涉技术(TDI)利用对三星编队上的多组测相值进行延时处理和线性组合,构造出了等臂长干涉仪,从而消除激光源噪声。结合星间绝对距离1米的测量目标,将抑制激光源噪声8个数量级。为了消除时钟噪声,通过时钟边带调制的方案产生新的拍频信号和测相值,并通过对测相值作差提取出时钟噪声项,从而修正TDI数据组合,将目前的星载时钟噪声抑制3个数量级。由于在降低测距误差中消除激光源噪声和时钟噪声具有最高的优先级,因此TDI数据处理方式将成为未来引力波探测的核心机制。另外,TDI技术对我国下一代基于激光干涉测距的重力场测量卫星载荷的研发也具有重要的参考价值。实际上双星激光干涉测距采用的锁相转发机制能够等价为两组单向测相值的TDI数据组合[39],进而可以降低光学锁相精度要求。可以看到,开展TDI方面的深入研究对我国未来开展包括引力波探测在内的天基激光精密测量任务具有重要的价值。

参考文献 (39)

目录

    /

    返回文章
    返回