CN116413832B - Time sequence data-based space gravitational wave detection sensitivity calculation method and system - Google Patents
Time sequence data-based space gravitational wave detection sensitivity calculation method and system Download PDFInfo
- Publication number
- CN116413832B CN116413832B CN202310306618.XA CN202310306618A CN116413832B CN 116413832 B CN116413832 B CN 116413832B CN 202310306618 A CN202310306618 A CN 202310306618A CN 116413832 B CN116413832 B CN 116413832B
- Authority
- CN
- China
- Prior art keywords
- noise
- time sequence
- sequence data
- sensitivity
- spectral density
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 230000035945 sensitivity Effects 0.000 title claims abstract description 70
- 238000001514 detection method Methods 0.000 title claims abstract description 40
- 238000004364 calculation method Methods 0.000 title claims abstract description 18
- 230000003595 spectral effect Effects 0.000 claims abstract description 41
- 238000000034 method Methods 0.000 claims abstract description 28
- 238000005516 engineering process Methods 0.000 claims abstract description 13
- 238000005305 interferometry Methods 0.000 claims abstract description 12
- 238000012545 processing Methods 0.000 claims abstract description 12
- 230000003287 optical effect Effects 0.000 claims description 26
- 238000004088 simulation Methods 0.000 claims description 20
- 238000012360 testing method Methods 0.000 claims description 16
- 238000005316 response function Methods 0.000 claims description 9
- 238000001228 spectrum Methods 0.000 claims description 7
- 238000012372 quality testing Methods 0.000 claims description 2
- 238000012546 transfer Methods 0.000 abstract description 7
- 230000008859 change Effects 0.000 description 7
- 238000005259 measurement Methods 0.000 description 7
- 230000015572 biosynthetic process Effects 0.000 description 6
- 238000011160 research Methods 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 5
- 238000006073 displacement reaction Methods 0.000 description 4
- 238000013461 design Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000012795 verification Methods 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000005452 bending Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 238000004880 explosion Methods 0.000 description 1
- 230000005486 microgravity Effects 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000010287 polarization Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000000241 respiratory effect Effects 0.000 description 1
- 238000010206 sensitivity analysis Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V13/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices covered by groups G01V1/00 – G01V11/00
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V7/00—Measuring gravitational fields or waves; Gravimetric prospecting or detecting
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/10—Noise analysis or noise optimisation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
The invention relates to the gravitational wave detection field, in particular to a method and a system for calculating the detection sensitivity of a space gravitational wave based on time sequence data. The invention calculates the sensitivity of a spatial gravitational wave detector using mixed noise time sequence data, comprising: acquiring multi-channel mixed noise phase time sequence data; processing the multi-path mixed noise phase time sequence data by using a time delay interferometry technology to generate one-path noise phase time sequence data; converting the noise phase timing data into noise strain timing data; performing Fourier transform on the noise strain time sequence data to calculate noise strain power spectral density; converting the noise strain power spectral density to a noise relative frequency power spectral density; the sensitivity is calculated using the noise versus frequency power spectral density and the transfer function R. The invention realizes the calculation of TDI1.0X combined sensitivity by utilizing the mixed noise time sequence data and the resolved transfer function formula.
Description
Technical Field
The invention relates to the gravitational wave detection field, in particular to a method and a system for calculating the detection sensitivity of a space gravitational wave based on time sequence data.
Background
In the beginning of 2016, the us land laser interference gravitational wave astronomy (LIGO) published results of direct detection of gravitational waves, confirming the prediction of einstein about the existence of gravitational waves one hundred years ago, and opening up a new research direction, gravitational wave astronomy. Unlike electromagnetic waves, gravitational waves will create a completely new window for us to observe the universe. From a large explosion in the universe to a black hole, from a merging of the intermediate stars to a double star of the white dwarf, from a topological defect in the universe to dark energy of a dark substance, and the like.
The ground gravitational wave measurement is limited by the influence of the earth gravitational force and the size of the scale, and the detected gravitational wave range is mainly concentrated at high frequency. Spatial (in space) gravitational wave detection is not an option, as it is understood that a wider range of gravitational waves and wave sources are of a nature. As a representative of the spatial gravitational wave detection items, LISA (Laser Interferometer Space Antenna) describes clear paths and platforms for the spatial gravitational wave detection subjects at the task concept level, provides references for the design of other spatial laser interference gravitational wave detection items such as tai chi plan, tian qin plan and the like in the world, particularly promotes the development of the spatial gravitational wave detection subjects in the aspects of wave source analysis, data processing, relevant wave source astronomy, relative theory itself and the like, and promotes the spatial gravitational wave detection to gradually form a subject research field.
The method is initiated by the institute of science of China in 2008, a plurality of units of the academy of sciences and a plurality of surgical research units of the academy of sciences participate together, so that a space gravitational wave detection demonstration group of the academy of sciences is established, and a development route pattern of space gravitational wave detection of China in the coming decades is planned. Through the pilot pre-research of two-stage Chinese academy, a space gravitational wave detection Taiji plan is provided. In 2014, the university of Zhongshan and the university of Huazhong science and technology jointly propose a Tianqin plan. In 2019, the Taiji team and the Tianqin team successfully launched the technical experimental star Taiji one number and Tianqin one number respectively.
Taking a Taiji plan space gravitational wave detector as an example, the geometry of the Taiji plan space gravitational wave detector is triangular by taking three spacecrafts as vertexes, each spacecraft is provided with two testing quality and two lasers, and each side is an arm. The method has high measurement accuracy requirement, strong coupling of each link and great technical challenge, and the technical scheme and key technology are still in the process of deep research and perfection at present, and comprehensive optimization and verification means of core indexes and the technical scheme are urgently needed in the system level. Because the detection system has strict requirements on microgravity and electromagnetic environment, and is difficult to establish a complete equivalent simulation loop model on the ground to develop system-level function and performance verification, the full-link numerical simulation analysis is an effective means for quickly finding potential problems and risks of tasks, comprehensively optimizing the overall design scheme of the system and assisting in the research, development and verification of key technologies. In order to achieve the purpose, the national center of space science of China academy of sciences constructs a preliminary framework of a space gravitational wave detection full-link simulation system (hereinafter referred to as a simulation system) and has the capabilities of top detection index sensitivity analysis, track and formation design, noise analysis, time sequence signal simulation and the like.
In order to determine whether a particular gravitational wave source can be detected by the gravitational wave detector, it is necessary to know the sensitivity limit of the instrument, the detector sensitivity reflecting the sensitivity of the detector to the target signal. The sensitivity of these instruments is typically described by plotting the square root (ASD, amplitude spectral density, amplitude spectral density) of the processed instrument total noise PSD (power spectral density ) versus transfer function ratio. This means that the instrument noise level affects the sensitivity, the smaller the noise level, the higher the sensitivity. The previous sensitivity calculation is generally based on two detector noise top-level indexes, namely a displacement noise ASD index and an inertial sensor noise (inertial transmission noise) ASD index. And for convenience, the displacement noise only considers the main mattersShot noise and laser phase noise. The laser phase noise of the detector is quite large, and for a ground interferometer, the fixed lengths of the two arms are the same, so that the laser experiences the same delay, and the laser phase noise can be well eliminated. However, for a spatial interferometer, it is not possible to keep the length of each arm constant. This results in lasers in different arms having different delays and residual laser phase noise greatly affects sensitivity. To solve this problem, tinto et al. First, time delay interferometry (TDI, time Delay Interferometry) techniques were proposed to eliminate noise with different measurement data combinations. The laser phase noise is eliminated by using the combination of the X link under TDI1.0, so that only inertial noise and shot noise are considered in calculation, and according to the noise index requirement of Taiji III, the shot noise ASD is 5 multiplied by 10 -12 m/≡Hz, inertial noise ASD of 3×10 -15 ms -2 V Hz, combined with a theoretical formula. The index-based approach obviously does not truly reflect the sensitivity of the detector.
Disclosure of Invention
The invention aims to provide a method for calculating the detection sensitivity of a space gravitational wave based on time sequence data, which can utilize the mixed noise time sequence data simulated by a simulation system to calculate the sensitivity, and can provide a method support for the analysis of the simulation precision of a mixed noise model of the simulation system compared with the sensitivity based on indexes; on the other hand, the method can calculate more real sensitivity by using real detection data obtained by a future space gravitational wave detection formation system.
In order to achieve the above purpose, the present invention is realized by the following technical scheme.
The invention provides a time sequence data-based spatial gravitational wave detection sensitivity calculation method, which is used for calculating the sensitivity of a spatial gravitational wave detector by using mixed noise time sequence data, and comprises the following steps:
s1, acquiring phase time sequence data of multiple paths of mixed noise;
s2, processing the multi-path mixed noise phase time sequence data by using a time delay interferometry technology to generate one-path noise phase time sequence data;
s3, converting the noise phase time sequence data into noise strain time sequence data;
s4, carrying out Fourier transform on the noise strain time sequence data to calculate noise strain power spectral density;
s5, converting the noise strain power spectral density into noise relative frequency power spectral density;
s6, calculating sensitivity by using the noise relative frequency power spectral density and the response function R.
As one of the improvements of the above technical solution, in the step S1, a simulation system is adopted to simulate and obtain multiple paths of mixed noise phase time sequence data, or a transmitted spatial gravitational wave detector is used to obtain multiple paths of mixed noise phase time sequence data.
As one of the improvements of the technical scheme, the space gravitational wave detector comprises three spacecrafts with two lines in a triangle shape;
each spacecraft comprises a master optical bench and a slave optical bench;
each of the master optical bench and the slave optical bench comprises a test quality, a laser and an interferometry device;
each interferometry device includes a scientific interferometer, a test mass interferometer, and a reference interferometer;
the scientific interferometer is used for measuring the distance between corresponding spacecrafts according to the interference of laser emitted by lasers between a main optical bench or a slave optical bench of one of the spacecrafts and a slave optical bench or a main optical bench of other spacecrafts;
the test quality interferometer is used for measuring the distance between the test quality in each spacecraft according to the interference of laser emitted by the laser between the master and slave tables of each spacecraft;
the reference interferometer is used for providing reference phases of laser interference for the scientific interferometer and the quality testing interferometer.
As one of the improvements of the above technical solution, at least 18 paths of mixed noise phase timing data are acquired in the step S1 using all of the scientific interferometer, the test quality interferometer and the reference interferometer.
As one of the improvements of the technical scheme, the noise strain power spectral density generated by S4 is multiplied by (2pi fL/c) 2 Obtaining the relative frequency power spectral density of noisef is the frequency, L is the distance between any two spacecrafts, and c is the speed of light.
As one of the improvements of the above technical solution, the formula for calculating the sensitivity S (u) in the step S6 is:
wherein T is the observation time.
The invention also provides a space gravitational wave detection sensitivity calculation system based on time sequence data, which comprises:
the noise phase time sequence data input module is used for inputting the multipath mixed noise phase time sequence data acquired by the space gravitational wave detector;
the noise relative frequency power spectrum density module is used for processing the multi-path mixed noise phase time sequence data by utilizing a time delay interferometry technology to generate one-path noise phase time sequence data; the method comprises the steps of converting noise phase time sequence data into noise strain time sequence data, performing Fourier transform on the noise strain time sequence data to calculate noise strain power spectral density, and finally converting the noise strain power spectral density into noise relative frequency power spectral density; and
and the sensitivity calculation module is used for calculating the sensitivity by utilizing the noise relative frequency power spectral density and the response function R.
As one of the improvements of the technical scheme, in the noise relative frequency power spectral density module, the noise strain power spectral density generated by S4 is multiplied by (2 pi fL/c) 2 The power spectrum density of the noise relative to the frequency is obtained, f is the frequency, L is the distance between any two spacecrafts, and c is the speed of light.
As one of the improvements of the above technical solution, in the sensitivity calculation module, the calculation formula of the sensitivity S (u) is:
wherein T is the observation time.
Compared with the prior art, the invention has the advantages that:
1. the sensitivity of TDI1.0X combination is calculated by combining the mixed noise time sequence data simulated by the simulation system with the analyzed transfer function formula, which is unprecedented;
2. the method can calculate more real sensitivity different from a theoretical value by utilizing real detection data obtained by a future space gravitational wave detection formation system.
Drawings
FIG. 1 is a schematic diagram of a spatial laser interference gravitational wave detection system formation;
FIG. 2 is a PSD of one path of noise timing data after TDI1.0-X combination processing;
FIG. 3 is a flow chart of an embodiment of the present invention;
FIG. 4 is a plot of theoretical sensitivity under TDI1.0X combinations based on noise top level indicators;
FIG. 5 is a plot of sensitivity under TDI1.0X combinations based on several primary noise timing data;
FIG. 6 is a plot of sensitivity under TDI1.0X combinations based on all noise timing data;
FIG. 7 is a plot of sensitivity in TDI1.0X combination with some gravitational wave sources based on all noise timing data;
fig. 8 is a diagram of gravitational wave sources.
Detailed Description
The technical scheme of the invention is described in detail below with reference to the accompanying drawings and examples.
FIG. 1 is a schematic diagram of a spatial laser interference gravitational wave detector formation; three spacecrafts (SC 1, SC2 and SC 3) of the space gravitational wave detector are triangular, each SC comprises a master optical table and a slave optical table, a testing quality, a laser and an interferometry device are arranged on the optical table, a scientific interferometer link i (i=1, 2 and 3) and a scientific interferometer link is in the figure represent six interference arms corresponding to six scientific interferometers, for example, a scientific interferometer link 1s represents that the laser of the SC1 local laser of the optical table interferes with the laser emitted by the SC3 main optical table laser; similarly, the test quality interferometer links i and is are six interferometer arms corresponding to six test quality interferometers, and the reference interferometer links i and is are six interferometer arms corresponding to six reference interferometers. Each optical platform has a metallic mass thereon, known as the proof mass. According to the generalized relativity theory, when gravitational waves pass, small changes generated by the space-time structure cause the space distance between two tested masses to change, and the laser interference principle is used for measuring the small changes. The scientific interferometer measures the displacement change between two satellites, the interference light intensity is 100pw and 10mw, and the displacement measurement precision is required to be lpm/Hz; the reference interferometer provides reference phases for other interferometers to eliminate noise caused by non-optical platforms, such as optical fibers, frequency modulation devices and the like, and provides phase-locked error signals for weak light phase locking; the test mass interferometer provides inertial sensing, scientific interferometers with positional information of the test mass relative to the optical platform.
The gravitational wave can cause interference arm length change when passing through the detector, and 18 interferometers can measure laser phase change caused by the gravitational wave. In addition, noise also causes laser phase change, and the simulation system simulates 18 paths of phase change time sequence data caused by the gravitational wave and the noise. When simulating, only various noises are simulated, and the simulation system outputs 18 paths of mixed noise time sequence data.
The method for calculating the sensitivity theoretical model based on the noise top-level index TDI1.0X combination (TDI1.0X combination specifically means that the TDI technology eliminates laser phase noise by carrying out phase delay and linear combination on phases measured by a plurality of interferometers, TDI1.0 means that the arm length is fixed without considering spacecraft constellation rotation and scientific interferometer interference arm telescopic bending, and does not change with time and does not distinguish the difference of interference arms when the same arm is shot on and shot off, X combination is also called Michelson combination, and is an important combination mode of the TDI technology), and the method for calculating the sensitivity theoretical model is as follows:
noise indicators such as Taiji threeThe shot noise ASD was found to be 5X 10 -12 m/≡Hz, inertial noise ASD of 3×10 -15 ms -2 V Hz, substituting equation (1) to calculate the noise PSD.
To the left of the equation is the noise PSD (this equation noise PSD is in relative frequency), subscript X 1 Representing TDI1.0X combinations. Right side s of equation a Is inertial noise ASD, s x For shot noise ASD, L is 3000000000m, c is the speed of light 300000000m/s, u is 2pi fL/c, f is the frequency, the frequency band is 0Hz to 1.5Hz, every interval is 1.5X10 -5 Hz takes a frequency bin.
The noise PSD is taken into equation (2) to calculate the sensitivity S.
Wherein T is the observation time 31536000s, R is the average response function of the space gravitational wave detector on the whole celestial sphere obtained by averaging the gravitational wave source position and the polarization angle. Conventionally, sensitivity is typically calculated using an all-celestial average response function. In previous studies, the response function was numerically simulated or semi-resolved, and in recent years, an average response function was derived for resolution, and the TDI1.0 generation X combined resolved response function was used.
The TDI1.0X combined sensitivity theoretical model calculation method based on time sequence data is as follows:
s1, acquiring multi-channel mixed noise phase time sequence data by using a spatial gravitational wave detector;
s2, processing the multi-path mixed noise phase time sequence data by utilizing a time delay interference TDI1.0X combined measurement technology to generate one-path noise phase time sequence data;
s3, converting the noise phase time sequence data into a relation between the noise strain time sequence data phase and strain as shown in a formula (3), wherein h is the strain,is phase, lambda is laser wavelength, 1064 nm;
s4, carrying out Fourier transform on the noise strain time sequence data to calculate noise strain power spectral density;
s5, multiplying the noise strain power spectral density by (2pi fL/c) 2 Converting into noise relative frequency power spectral density;
s6, calculating the sensitivity by combining the noise relative frequency power spectral density and the transfer function R with the formula (2).
Example 1
The method for calculating the spatial gravitational wave detection sensitivity curve based on time sequence data utilizes one path of noise time sequence data generated by TDI processing of 18 paths of noise time sequence data of an interferometer to calculate noise PSD, and then calculates sensitivity. The noise relative frequency PSD calculated by this method is shown in FIG. 2.
A specific implementation flowchart is shown in fig. 3.
The method specifically comprises the following steps:
s1, acquiring multi-channel mixed noise phase time sequence data by using a spatial gravitational wave detector;
s2, processing the multi-path mixed noise phase time sequence data by utilizing a time delay interference TDI1.0X combined measurement technology to generate one-path noise phase time sequence data;
s3, converting the noise phase time sequence data into a relation between the noise strain time sequence data phase and strain as shown in a formula (3), wherein h is the strain,is phase, lambda is laser wavelength, 1064 nm;
s4, carrying out Fourier transform on the noise strain time sequence data to calculate noise strain power spectral density;
s5, multiplying the noise strain power spectral density by (2pi fL/c) 2 Converting into noise relative frequency power spectral density;
s6, calculating the sensitivity by combining the noise relative frequency power spectral density and the transfer function R with the formula (2).
As fig. 5 shows a sensitivity curve calculated based on several main noise items mixed with noise simulation time sequence data, only shot noise, inertial noise and laser phase noise are considered, and it can be seen from comparison with fig. 4 that the method is basically correct.
Fig. 6 is a sensitivity curve calculated based on all noise term mixed noise simulation timing data, accounting for shot noise, laser frequency noise, laser pointing jitter noise, lead pointing path noise, respiratory angle noise, phase meter noise, inertial total noise, and other noise.
The bottom curve in fig. 7 is a sensitivity curve calculated based on the mixed noise simulation time sequence data of all noise items, and other lines are some wave sources. The gravitational wave source is as shown in fig. 8: wherein J0651+2844 and J0934+4411 are double-white dwarf systems. The wave sources are already marked in the figure.
Example 2
The embodiment 2 of the invention provides a space gravitational wave detection sensitivity calculation system based on time sequence data, and the sensitivity can be obtained after the acquired multipath mixed noise phase time sequence data are input into the system. The system comprises:
the noise phase time sequence data input module is used for inputting multi-path mixed noise phase time sequence data simulated by the simulation system;
the noise phase time sequence data processing module is used for processing the multi-path mixed noise phase time sequence data by utilizing a time delay interference TDI1.0X combination measurement technology to generate one-path noise phase time sequence data;
the noise relative frequency power spectrum density module is used for converting the noise phase time sequence data into noise strain time sequence data, then carrying out Fourier transform on the noise strain time sequence data to calculate noise strain power spectrum density, and finally converting the noise strain power spectrum density into noise relative frequency power spectrum density; and
and the sensitivity calculation module is used for calculating the sensitivity by utilizing the noise relative frequency power spectral density and the transfer function R.
According to the invention, the spatial gravitational wave detection sensitivity curve is calculated based on the noise time sequence data, on one hand, the sensitivity can be calculated by utilizing the mixed noise time sequence data simulated by the simulation system, and compared with the sensitivity based on indexes, the method support is provided for the analysis of the simulation precision of the mixed noise model of the simulation system; on the other hand, method support is provided for calculating more realistic sensitivity for real detection data obtained by utilizing a future space gravitational wave detection formation system.
Finally, it should be noted that the above embodiments are only for illustrating the technical solution of the present invention and are not limiting. Although the present invention has been described in detail with reference to the embodiments, it should be understood by those skilled in the art that modifications and equivalents may be made thereto without departing from the spirit and scope of the present invention, which is intended to be covered by the appended claims.
Claims (9)
1. A method for calculating sensitivity of a spatial gravitational wave detection based on time series data for calculating sensitivity of a spatial gravitational wave detector using mixed noise time series data, the method comprising:
s1, acquiring phase time sequence data of multiple paths of mixed noise;
s2, processing the multi-path mixed noise phase time sequence data by using a time delay interferometry technology to generate one-path noise phase time sequence data;
s3, converting the noise phase time sequence data into noise strain time sequence data;
s4, carrying out Fourier transform on the noise strain time sequence data to calculate noise strain power spectral density;
s5, converting the noise strain power spectral density into noise relative frequency power spectral density;
s6, calculating sensitivity by using the noise relative frequency power spectral density and the response function R.
2. The method according to claim 1, wherein in step S1, the simulation system is adopted to simulate and obtain multiple paths of mixed noise phase time sequence data, or the transmitted spatial gravitational wave detector is used to obtain multiple paths of mixed noise phase time sequence data.
3. The method for calculating the detection sensitivity of the spatial gravitational wave based on time series data according to claim 1, wherein the spatial gravitational wave detector comprises three spacecrafts with two lines in a triangle shape;
each spacecraft comprises a master optical bench and a slave optical bench;
each of the master optical bench and the slave optical bench comprises a test quality, a laser and an interferometry device;
each interferometry device includes a scientific interferometer, a test mass interferometer, and a reference interferometer;
the scientific interferometer is used for measuring the distance between corresponding spacecrafts according to the interference of laser emitted by lasers between a main optical bench or a slave optical bench of one of the spacecrafts and a slave optical bench or a main optical bench of other spacecrafts;
the test quality interferometer is used for measuring the distance between the test quality in each spacecraft according to the interference of laser emitted by the laser between the main optical bench and the auxiliary optical bench of each spacecraft;
the reference interferometer is used for providing reference phases of laser interference for the scientific interferometer and the quality testing interferometer.
4. The method for calculating the spatial gravitational wave detection sensitivity based on time series data according to claim 3, wherein said step S1 uses all of the scientific interferometer, the test quality interferometer and the reference interferometer to obtain at least 18 paths of mixed noise phase time series data.
5. The method of calculating sensitivity of spatial gravitational wave detection based on time series data as claimed in claim 3, wherein in said step S5, the noise strain power spectral density generated in S4 is multiplied by (2πfL/c) 2 Obtaining the relative frequency power spectral density of noisef is the frequency, L is the distance between any two spacecrafts, and c is the speed of light.
6. The method for calculating the sensitivity of spatial gravitational wave detection based on time series data as claimed in claim 5, wherein the formula for calculating the sensitivity S (u) in step S6 is:
wherein T is the observation time.
7. A spatial gravitational wave detection sensitivity calculation system based on time series data, the system comprising:
the noise phase time sequence data input module is used for inputting the multipath mixed noise phase time sequence data acquired by the space gravitational wave detector;
the noise relative frequency power spectrum density module is used for processing the multi-path mixed noise phase time sequence data by utilizing a time delay interferometry technology to generate one-path noise phase time sequence data; the method comprises the steps of converting noise phase time sequence data into noise strain time sequence data, performing Fourier transform on the noise strain time sequence data to calculate noise strain power spectral density, and finally converting the noise strain power spectral density into noise relative frequency power spectral density; and
and the sensitivity calculation module is used for calculating the sensitivity by utilizing the noise relative frequency power spectral density and the response function R.
8. The time series data based spatial gravitational wave detection sensitivity calculation system as claimed in claim 7 wherein in said noise versus frequency power spectral density module, the generated noise strain power spectral density is multiplied by (2pi fL/c) 2 Obtaining the relative frequency power spectral density of noisef is the frequency, L is the distance between any two spacecrafts, and c is the speed of light.
9. The time series data based spatial gravitational wave detection sensitivity calculation system of claim 8, wherein in said sensitivity calculation module, the calculation formula of the sensitivity S (u) is:
wherein T is the observation time.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310306618.XA CN116413832B (en) | 2023-03-27 | 2023-03-27 | Time sequence data-based space gravitational wave detection sensitivity calculation method and system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310306618.XA CN116413832B (en) | 2023-03-27 | 2023-03-27 | Time sequence data-based space gravitational wave detection sensitivity calculation method and system |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116413832A CN116413832A (en) | 2023-07-11 |
CN116413832B true CN116413832B (en) | 2023-11-10 |
Family
ID=87055839
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310306618.XA Active CN116413832B (en) | 2023-03-27 | 2023-03-27 | Time sequence data-based space gravitational wave detection sensitivity calculation method and system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116413832B (en) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1118866A1 (en) * | 2000-01-20 | 2001-07-25 | Tektronix, Inc. | Method of estimating phase noise spectral density and jitter in a periodic signal |
CN106706124A (en) * | 2016-11-29 | 2017-05-24 | 浙江大学 | Method of measuring light source relative intensity noise power spectrum density |
CN109471188A (en) * | 2018-10-15 | 2019-03-15 | 吉林大学 | A kind of gravitation wave detecting method based on pseudorandomcode technology |
CN115437026A (en) * | 2022-08-22 | 2022-12-06 | 中国科学院国家空间科学中心 | Method and system for making frequency planning scheme of space-based gravitational wave detector |
CN115561824A (en) * | 2022-10-19 | 2023-01-03 | 中国科学院国家空间科学中心 | Time delay-based spatial gravitational wave detection sensitivity curve calculation method |
-
2023
- 2023-03-27 CN CN202310306618.XA patent/CN116413832B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1118866A1 (en) * | 2000-01-20 | 2001-07-25 | Tektronix, Inc. | Method of estimating phase noise spectral density and jitter in a periodic signal |
CN106706124A (en) * | 2016-11-29 | 2017-05-24 | 浙江大学 | Method of measuring light source relative intensity noise power spectrum density |
CN109471188A (en) * | 2018-10-15 | 2019-03-15 | 吉林大学 | A kind of gravitation wave detecting method based on pseudorandomcode technology |
CN115437026A (en) * | 2022-08-22 | 2022-12-06 | 中国科学院国家空间科学中心 | Method and system for making frequency planning scheme of space-based gravitational wave detector |
CN115561824A (en) * | 2022-10-19 | 2023-01-03 | 中国科学院国家空间科学中心 | Time delay-based spatial gravitational wave detection sensitivity curve calculation method |
Non-Patent Citations (2)
Title |
---|
空间激光干涉引力波探测;罗子人;白姗;边星;陈葛瑞;董鹏;董玉辉;高伟;龚雪飞;贺建武;李洪银;李向前;李玉琼;刘河山;邵明学;宋同消;孙保三;唐文林;徐鹏;徐生年;杨然;靳刚;;力学进展(第04期);全文 * |
罗子人 ; 白姗 ; 边星 ; 陈葛瑞 ; 董鹏 ; 董玉辉 ; 高伟 ; 龚雪飞 ; 贺建武 ; 李洪银 ; 李向前 ; 李玉琼 ; 刘河山 ; 邵明学 ; 宋同消 ; 孙保三 ; 唐文林 ; 徐鹏 ; 徐生年 ; 杨然 ; 靳刚 ; .空间激光干涉引力波探测.力学进展.2013,(第04期),全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN116413832A (en) | 2023-07-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Bailes et al. | Gravitational-wave physics and astronomy in the 2020s and 2030s | |
Otto | Time-delay interferometry simulations for the laser interferometer space antenna | |
Rubbo et al. | Forward modeling of space-borne gravitational wave detectors | |
Hartig et al. | Geometric tilt-to-length coupling in precision interferometry: mechanisms and analytical descriptions | |
CN103308952A (en) | Gravitational wave detection device design and method thereof | |
Wang et al. | Algorithm for time-delay interferometry numerical simulation and sensitivity investigation | |
CN110146258A (en) | A kind of phase recovery method when under Poisson noise model to expansion target | |
CN116413832B (en) | Time sequence data-based space gravitational wave detection sensitivity calculation method and system | |
Vigier | Relativistic interpretation (with non-zero photon mass) of the small ether drift velocity detected by Michelson, Morley and Miller | |
CN115561824B (en) | Space gravitational wave detection sensitivity curve calculation method based on time delay | |
Sivanathan et al. | Development of a hardware for frequency scanning interferometry for long range measurement | |
Muratore et al. | Effectiveness of null time-delay interferometry channels as instrument noise monitors in LISA | |
CN103675935A (en) | Design method of multifunctional space gravitational wave detector based on TRIZ | |
CN112325808B (en) | Flatness real-time calibration compensation measurement method based on multiple PSDs | |
Gondoin et al. | Darwin-GENIE: a nulling instrument at the VLTI | |
Tinto et al. | LISA time-delay interferometry zero-signal solution: Geometrical properties | |
Wanner et al. | In-depth modeling of tilt-to-length coupling in LISA’s interferometers and TDI Michelson observables | |
Jiang et al. | A compact high-precision frequency swept interferometer (FSI) for autonomous ranging of satellite constellation | |
Passe et al. | A metrology system for large deployable space structures using optical frequency domain reflectometry | |
Quirrenbach et al. | The PRIMA astrometric planet search project | |
Schutz | Astrometric and timing effects of gravitational waves | |
Crosta et al. | On the principle of Astrometric Gravitational Wave Antenna | |
Blandford | Gravitational lensing and the extragalactic distance scale By ROGER D. BLANDFORD AND TOMISLAV KUNDIĆ | |
Larson | LISA: A modern astrophysical observatory | |
Massimo et al. | Time-delay interferometry |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |