Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Studying an off-axis optical bench for future gravity missions from the perspective of carrier-to-noise ratio

Open Access Open Access

Abstract

The inter-satellite laser ranging heterodyne interferometer is vital for future gravity missions to achieve high ranging accuracy. This paper proposes a novel off-axis optical bench design which integrates merits of the off-axis optical bench design of GRACE Follow-On mission and other on-axis designs. This design makes use of lens systems subtly to restrict the tilt-to-length coupling noise and takes advantage of the DWS feedback loop to maintain the transmitting beam and receiving beam anti-parallel. The critical parameters of the optical components are determined and the carrier-to-noise ratio for a single channel of the photoreceiver is calculated to be more than 100 dB-Hz for the high case. The off-axis optical bench design is a potential candidate for China’s future gravity missions.

© 2023 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement

1. Introduction

The Earth gravity field measurements, especially the time-variable gravity field measurements, are able to give a wealth of knowledge about the hydrologic cycle, sea-level changes, glacier mass changes, as well as mass movements inside Earth from seasonal to decadal and provide important help for solving the increasingly serious problems of resources, environment and disasters facing mankind. The global and homogeneous gravity field can be obtained by means of the inter-satellite ranging. In 2002, the Gravity Recovery and Climate Experiment (GRACE) mission [1], which tracked the distance between two satellites and its rate of change by using inter-satellite microwave system in K/Ka band to derive the model of the Earth gravity field, was launched and worked until 2017. Its successor, GRACE Follow-On mission [2], launched successfully in 2018, was loaded with an inter-satellite laser ranging interferometer(LRI) as a technology demonstrator for the first time and achieved the ranging precision of $10 \ nm/\sqrt {Hz}$ at a frequency of $40\ mHz$.

With the success of the GRACE and GRACE Follow-On missions, the researches on future gravity missions with higher ranging precision and longer lifetime are ongoing. Different to the off-axis optical bench design of GRACE Follow-On LRI project, which produces the lateral offset between the transmitting(TX) beam and receiving(RX) beam and keeps the beams anti-parallel with the method of triple mirrors assembly(TMA), several kinds of on-axis optical bench designs that use the polarized light and multiple sets of lens system have been proposed in recent years [35]. By combining advantages of less optical components of GRACE Follow-On LRI project and higher carrier-to-noise ratio(CNR) of the on-axis designs with the help of telescopes, the new off-axis topology of the optical bench for future gravity missions is put forward in this paper. In this design, five sets of lens system are used to restrain the tilt-to-length(TTL) coupling noise and the differential wavefront sensing(DWS) feedback loop is applied to make sure the TX beam and RX beam anti-parallel to each other. Furthermore, the higher CNR for a single channel of the photoreceiver will be realized, as the following content will demonstrate. Hence, the off-axis optical bench design will be a competitive choice for China’s future gravity missions and serves as a technology demonstrator for the ESA-NASA Laser Interferometer Space Antenna (LISA) mission [6], as well as the Chinese space-borne gravitational wave observatories Taiji [7] and TianQin [8] missions.

The structure of this paper is organized as follows. In Section 2, we introduce the design principle of the off-axis optical bench. In Section 3, the calculating method of the CNR for one segment of the photodiode is discussed in detail. The key parameters of the optical bench design are determined and results of the CNR are given in Section 4. Finally, we draw our conclusions in Section 5.

2. Design concept of the off-axis optical bench

As its name implies, the off-axis optical bench design means that the TX beam and RX beam are anti-parallel to each other but not on the same axis, as shown in Fig. 1. Clipped by the receiver aperture at first, the RX beam from the distant satellite passes through the lens L1, then the most part of the RX beam transmits through the beamsplitter BS1. Next, about one half of the beam penetrates the beamsplitter BS2, then passes through the lens L4 and finally arrives at the quadrant photodiode QPD1. The other half is reflected by the BS2, then passes through the lens L5 and is detected by the QPD2. On the other hand, the local(LO) Gaussian beam, whose waist is set to be at the center point of the reflecting surface of the fast steering mirror(FSM), is injected from the fiber coupler to the optical bench and is reflected by the FSM, then transmits through the lens L2. The majority of the LO beam penetrates the BS1, lens L3, then is reflected by the mirrors M1, M2 in turn and is transmitted to the remote satellite. It needs to be pointed that the use of two mirrors instead of a single mirror is for maintaining the direction of the TX beam relative to the optical axis. At the same time, the small part of the LO beam is reflected by the BS1 and goes along with the RX beam passing through the BS1. At the end of the light path, the heterodyne interference happens between the LO and RX beam at the QPDs which convert an optical power into an electric current by exploiting the inner photoelectric effect. The photocurrent is then converted to a voltage by a transimpedance amplifier in the photoreceiver. At last, the phasemeter track the phase of voltage signals provided by the photoreceiver and this phase carries information about changes in distance between the satellites.

 figure: Fig. 1.

Fig. 1. The diagrammatic drawing of the optical bench design for future gravity missions.

Download Full Size | PDF

In this design, there are five sets of Keplerian telescopes: L1-L4(L1-L5) as the receiving telescopes, L2-L4(L2-L5) as the local telescopes and L2-L3 as the transmitting telescope. For the receiving telescopes, the front focus of L1 coincides with the center of the receiver aperture and the back focus of L4(L5) is placed at the center of the active area of the QPD1(QPD2), in order to keep the optical path of the RX beam from the aperture to photodiodes unchanged when the RX beam rotates around the center point of the aperture, which is called the RX reference point(RX RP). Similarly, the front focus of the L2 is positioned at the center point of the reflecting surface of the FSM, and the back focus of the L3 is named as the TX reference point(TX RP). Therefore, the optical paths of both the LO beam from the FSM to photodiodes and the TX beam from the FSM to the TX RP remain constant even though the FSM rotates around the center point of its reflecting surface. The midpoint of the TX RP and RX RP is defined as the reference point(RP) of the optical bench, which coincides with the center of the test mass in the accelerometer and the mass center of the whole satellite. While the local satellite rotates around the RP, the DWS signal of the QPD1(or QPD2) is used to control the FSM attitude. As the DWS signal goes close to 0, the LO beam overlaps gradually with the RX beam and the TX and RX beam are nearly anti-parallel to each other in the meantime. At the same time, the contributions of displacements of the TX RP and RX RP to the ranging between two satellites cancel out so that the TTL coupling of the total measurement system is restricted significantly. Besides, the arrangement of double QPDs is for hot/cold redundancy as needed.

Next, the relations between transversal magnifications of the receiving, transmitting and local telescopes will be derived. $m_{RX}, m_{LO}$ and $m_{TX}$ denote the absolute value of transversal magnifications of the receiving, local and transmitting telescopes respectively. It needs to be noted that the transmitting telescope expands the TX beam while the receiving and local telescopes shrink the RX beam and LO beam respectively. The incident angles into telescopes and the emergent angles out of telescopes of the RX, TX and LO beam satisfy the following equations:

$$|\theta_{RX,out}|= |\theta_{RX,in}|/ m_{RX},$$
$$|\theta_{TX,out}|= |\theta_{TX,in}|/ m_{TX} ,$$
$$|\theta_{LO,out}|= |\theta_{LO,in}|/ m_{LO} .$$

In order to make the TX beam anti-parallel to the RX beam, the equation

$$\theta_{TX,out}=\theta_{RX,in}$$
needs to be satisfied. And it’s known from the optical bench design that the incident angles of the TX beam and the LO beam are equal, e.g.
$$\theta_{TX,in}=\theta_{LO,in}.$$

Combining all above equations, we have

$$m_{TX}=\frac{|\theta_{TX,in}|}{|\theta_{TX,out}|}=\frac{|\theta_{LO,in}|}{|\theta_{RX,in}|}=\frac{|\theta_{RX,out}|}{|\theta_{RX,in}|}\cdot \frac{|\theta_{LO,in}|}{|\theta_{LO,out}|}=\frac{m_{LO}}{m_{RX}}.$$

The relation will be helpful for the derivation of the CNR in the next section.

3. Calculations of carrier-to-noise ratio

A major function of the optical bench is to provide the heterodyne interference signal with the high CNR because the sensitivity of the phase readout of the phasemeter can be expressed as the inverse of the CNR $C/N_0$,

$$PSD[\delta \phi_{PM}]=\frac{1\ rad^2}{C/N_0}.$$

Here the unit of the CNR is Hz and it’s usually expressed by dB-Hz. The criteria is generally used for the phase tracking within GPS receivers [9], and has been applied in GRACE and GRACE Follow-On missions to evaluate system noise in the microwave system [10]. In order to ensure that the phasemeter continuously tracks the phase of the input signal without the loss of the phase lock, the minimum requirement for the CNR in GRACE Follow-On LRI project is 70 dB-Hz corresponding to $10^7$ Hz, which means a phase readout accuracy of $0.3\ mrad/\sqrt {Hz}$ according to eq. (7). This requirement will be applied in the following discussions about the CNR of the optical design. The calculations of the CNR below follow the approach presented in the reference [3], whereby the corresponding equations have been modified to account for our optical bench design.

The CNR is generally defined as the ratio of carrier power to noise power spectral density(PSD), which is used to estimate the quality of the beatnote signal. Therefore the CNR of a photoreceiver can be expressed as

$$C/N_0=\frac{carrier\ RMS\ power}{noise\ PSD}=\frac{I^2_{carrier}}{PSD[I_{noise}]},$$
where the impedance is cancelled because of the same impedance for both the carrier power and the noise PSD. $I_{carrier}$ is the effective value of the photocurrent. Based on the principle of laser interference, $I_{carrier}^2$ can be written in the form of
$$I^2_{carrier}=2\cdot \eta \cdot P_{RX} \cdot P_{LO} \cdot \eta^2_{PD},$$
where $P_{LO}$ and $P_{RX}$ denotes the local and receiving laser power respectively, $\eta$ is the heterodyne efficiency and $\eta _{PD}$ is the photon-to-electron conversion efficiency of the photodiode. The typical value of $\eta _{PD}$ is $0.8\ A/W$. In eq. (8), the PSD of the current noise $PSD[I_{noise}]$ can be divided into three dominating contributors as
$$PSD[I_{noise}]=PSD[I_{PR}]+PSD[I_{SN}]+PSD[I_{RIN}].$$

Here $I_{PR}$ is the equivalent input current noise of the photoreceiver [11] and its typical value of

$$PSD[I_{PR}]=25 pA^2/Hz, \ 4MHz<f<20MHz$$
has been shown in GRACE Follow-On LRI on-ground testing [3]. The shot noise originates from the energy quantization of light field and the statistics of arriving photons obeys Poisson distribution. Besides, the received laser power in the inter-satellite ranging is so weak that the strong local laser dominates the average number of photons. Hence, the single-sided PSD of the shot noise is
$$PSD[I_{SN}]=2\cdot e \cdot \eta_{PD} \cdot P_{LO},$$
where $e=1.609\times 10^{19}\ C$, is the electric charge quantity of an electron. At last, the RIN contribution comes from the relative intensity fluctuation of the laser and it can be expressed as
$$PSD[I_{RIN}]=PSD[\delta P_{LO}(f)/P_{LO}]\cdot (\eta_{PD}\cdot P_{LO})^2.$$

The $PSD[\delta P_{LO}(f)/P_{LO}]$, which denotes relative fluctuations of the local laser power, is required to be within $1.0\times 10^{-15}/Hz(f>4MHz)$ for gravity field measurements [12].

Consequently, the $C/N_0$ for a single channel of the photoreceiver can be expressed approximately as

$$C/N_0=\frac{2\cdot \eta \cdot P_{RX} \cdot P_{LO} \cdot \eta^2_{PD}}{25 pA^2/Hz+2\cdot e \cdot \eta_{PD} \cdot P_{LO}+1.0\times10^{{-}15}/Hz \cdot (\eta_{PD}\cdot P_{LO})^2}.$$

3.1 Optimal local beam power

In order to calculate the $C/N_0$ in eq. (14), we need to know these parameters: the LO beam power $P_{LO}$, RX beam power $P_{RX}$ and heterodyne efficiency $\eta$. For the sake of convenience, eq. (14) can be rewritten in the form of

$$C/N_0 =f(P_{LO})\cdot \frac{P_{RX}\cdot \eta}{1pW},$$
$$f(P_{LO}) =\frac{2 \cdot P_{LO} \cdot 1pW \cdot \eta^2_{PD}}{25 pA^2/Hz+2\cdot e \cdot \eta_{PD} \cdot P_{LO}+1.0\times10^{{-}15}/Hz \cdot (\eta_{PD}\cdot P_{LO})^2}.$$

It is evident that $f(P_{LO})$ is just the function of the LO beam power and it has a maximum at some value of $P_{LO}$. The function $f(P_{LO})$ is displayed in Fig. 2. After simple calculations, it comes to a conclusion that $f(P_{LO})$ reaches its maximum value 63.5 dB-Hz when $P_{LO}=0.198\ mW$. Therefore, when the LO beam power is set to be $0.198\ mW$, eq. (15) can also be changed to the form of

$$C/N_0=63.5 + 10\cdot log10\left(\frac{P_{RX}\cdot \eta}{1pW}\right)\ \text{dB-Hz}.$$

For meeting the preceding requirement of 70 dB-Hz, the second term on the right-hand side of the equal sign in the above equation needs to excess 6.5 dB-Hz. Next, the expression of the RX beam power $P_{RX}$ and the heterodyne efficiency $\eta$ is going to be determined.

 figure: Fig. 2.

Fig. 2. $C/N_0$ in dB-Hz as a function of the LO beam power for a single channel of the photoreceiver.

Download Full Size | PDF

3.2 Receiving beam power

It is obvious that the $C/N_0$ increases with the RX beam power from eq. (15). The laser beam intensity and the aperture size decide the power within the receiver aperture of the local satellite. The far-field intensity of the Gaussian beam for small misalignments of the TX beam relative to the line-of-sight is

$$I_{RX}\approx \frac{2\cdot P_{TX}}{\pi \omega (L)^2} \cdot e^{{-}2\alpha _{TX}^2\cdot L^2/\omega(L)^2},$$
where $P_{TX}$ denotes the TX beam power from the remote satellite, $\omega (L)$ is the spot size of the Gaussian beam that travels $L$ distance and $\alpha _{TX}$ represents the TX beam misalignment with respect to the line-of-sight. In addition, the effective area of the receiver aperture is
$$S=\pi \cdot r_{RAP}^2 \cdot \cos \beta_{RX}\approx \pi \cdot r_{RAP}^2 ,$$
where $r_{RAP}$ is the radius of the receiver aperture and $\beta _{RX}$, which denotes the misalignment angle of the receiver, is rather small so that $\cos \beta _{RX} \approx 1$. Combining eq. (18) and eq. (19), the laser beam power within the receiver aperture can be derived as
$$P_{RAP}=\frac{2\cdot P_{TX} \cdot r_{RAP}^2}{L^2 \theta_{TX}^2}\cdot e^{{-}2\alpha_{TX}^2/\theta_{TX}^2},$$
where the condition that $\omega _{L}\approx L\cdot \theta _{TX}$ is used and $\theta _{TX}$ is the divergence angle of the TX beam. however, the Gaussian beam tends to be clipped by the transmitting telescope when it’s transmitted from the satellite. The power in the receiver aperture for a clipped Gaussian beam can be expressed approximately as [3]
$$P_{RAP}\approx \frac{2\cdot P_{TX} \cdot r_{RAP}^2}{L^2 \theta_{TX}^2}\cdot (1-e^{{-}r_{TAP}^2/(\omega_{LO}\cdot m_{TX})^2})^2 \cdot e^{{-}2\alpha_{TX}^2/(\theta_{TX}\cdot \psi_1(\omega_{LO}\cdot m_{TX}/r_{TAP}))^2}.$$

Here $r_{TAP}$ is the radius of the transmitting aperture and $\omega _{LO}$ represents the waist radius of the LO beam. The $\psi _1(x)$ function, which accounts for the increased divergence angle due to the clipping, is a polynomial of 6th order with coefficients given in the Table 1 and is only valid in the range $0.2<x<2.0$. Considering the transmitting efficiency $\eta _{RX}$ for the RX beam from the receiver aperture to one segment of the photodiode and the relation

$$\theta_{TX}=\frac{\lambda}{\pi \cdot \omega_{LO}\cdot m_{TX}},$$
where $\lambda$ denotes the laser wavelength, the final expression of the power that is received by one segment of the photodiode is written as
$$\begin{aligned} P_{RX} \approx & \frac{2\pi^2 \cdot \eta_{RX} \cdot \omega_{LO}^2 \cdot m_{TX}^2 \cdot P_{TX} \cdot r_{RAP}^2}{L^2 \lambda^2}\cdot (1-e^{{-}r_{TAP}^2/(\omega_{LO}\cdot m_{TX})^2})^2 \\ & \quad \cdot e^{{-}2(\pi \alpha_{TX} \cdot \omega_{LO} \cdot m_{TX})^2/(\lambda \cdot \psi_1(\omega_{LO}\cdot m_{TX}/r_{TAP}))^2}. \end{aligned}$$

Tables Icon

Table 1. The coefficients of the polynomial functions $\psi _1(x)=\sum _{i=0}^{6} p_i\cdot x^i$ and $\psi _2(x)=\sum _{i=0}^{6} q_i\cdot x^i$. [3]

3.3 Heterodyne efficiency of interference

The heterodyne efficiency $\eta$ of the interference is another important factor in eq. (15) and raising the heterodyne efficiency means improving the overlap of RX and LO electric field. The expression of the heterodyne efficiency $\eta _{CPD}$ of the interference between the LO Gaussian beam and the RX beam, which has a constant intensity and phase within the aperture and zero magnitude outside, on the active area of a circular photodiode with radius $r_{PD}$ can be written as [3]

$$\eta_{CPD}=\frac{8\cdot \left| \int^{r_{PD}}_0 r\cdot e^{{-}r^2/\omega^2_{LO,PD}}\cdot J_0(2\pi r/\lambda \cdot \alpha_{RX,LO})dr\right|^2}{r_{PD}^2\cdot \omega^2_{LO,PD}\cdot (1-e^{{-}2r_{PD}^2/\omega^2_{LO,PD}})},$$
where $\omega _{LO,PD}$ represents the waist radius of the LO beam at the active area of the photodiode, $\alpha _{RX,LO}$ denotes the angle between directions of propagation of the RX and LO beams and $J_0$ is the Bessel function of first kind. When the direction of the LO beam coincides with that of the RX beam, i.e. $\alpha _{RX,LO}=0$, $\eta _{CPD}$ reaches the maximum and eq. (24) changes to
$$\eta_{CPD}=\frac{2\omega_{LO,PD}^2}{r_{PD}^2}\cdot \tanh\left(\frac{r_{PD}^2}{2\omega_{LO,PD}^2}\right).$$

Figure 3 shows that the heterodyne efficiency $\eta _{CPD}$ increases with the ratio of $\omega _{LO,PD}$ and $r_{PD}$. Moreover, $\eta _{CPD}$ is extremely close to 1 when $\omega _{LO,PD}\gg r_{PD}$.

 figure: Fig. 3.

Fig. 3. The graph of the heterodyne efficiency $\eta _{CPD}$ versus the ratio $\omega _{LO,PD}/r_{PD}$.

Download Full Size | PDF

If $\alpha _{RX,LO}\neq 0$, the heterodyne efficiency for a single segment can be derived approximately by the model [3]

$$\eta \approx \eta_{CPD}\cdot e^{{-}2\alpha_{RX,LO}^2/(\theta_{LO} \cdot \psi_2 (\omega_{LO,PD}/r_{PD})/m_{LO})^2},$$
where $\theta _{LO}$ is the divergence angle of the local Gaussian beam and the fitting function $\psi _2$, which is feasible for $x\in [0.2,2.0]$, is also a polynomial of 6th order with coefficients given in the second row of the Table 1.

According to the optical design, there are the following relations:

$$\omega_{LO,PD}=\omega_{LO}\cdot m_{LO}, r_{PD}=r_{RAP}\cdot m_{RX}, \alpha_{RX,LO}=\frac{\alpha_{RX}} {m_{RX}}, \theta_{LO}=\frac{\lambda}{\pi \omega_{LO}}, m_{TX}=\frac{m_{LO}}{m_{RX}}.$$

Thereupon, the full approximate expression of the heterodyne efficiency for a single segment of the photodiode is written as

$$\eta=\frac{2(\omega_{LO}\cdot m_{TX})^2}{r_{RAP}^2}\tanh\left(\frac{r_{RAP}^2}{2(\omega_{LO}\cdot m_{TX})^2}\right)e^{\frac{-2(\pi \omega_{LO}\cdot \alpha_{RX} \cdot m_{TX})^2}{(\lambda \cdot \psi_2(\omega_{LO}\cdot m_{TX}/r_{RAP}))^2}}$$

So far, the complete expression of the CNR for one channel of the photoreceiver has been derived by substituting eq. (28) and eq. (23) into eq. (17). Qualitatively, in addition to setting $P_{LO}$ for a single segment of the photodiode to be $0.198\ mW$, it’s vital for improving the CNR to increase the RX beam power and the heterodyne efficiency of interference between the RX and LO beams. In the following section, studying the CNR quantitatively will help determine the critical parameters of the optical bench design for future gravity missions.

4. Results and discussions

Before calculating results of the CNR, the values or ranges of the following parameters need to be ascertained:

$$\lambda;L;\omega_{LO};m_{TX};P_{TX};\eta_{RX};r_{RAP};r_{TAP};\alpha_{RX};\alpha_{TX}.$$

First of all, $1064\ nm$ wavelength laser with $100\ mW$ power and the mean distance $L=100\ km$ between two satellites are chosen in the preliminary plan for China’s future gravity missions. In the power budget of the optical design, $3\ mW$ laser is allocated for the frequency stabilization and around $1.6\ mW$ laser is assigned to be detected by eight segments of two quadrant photodiodes. A portion of the remaining laser is used as the TX beam without the clipping by the transmitting telescope and $P_{TX}$ can be described by

$$P_{TX}=\kappa \cdot (P_{laser}-P_{freq}-8\cdot P_{LO}),$$
where $\kappa$ represents the transmission efficiency from the fiber injector to the TX beam prior the clipping and is generally within the range of $[0.6,0.8]$ mainly due to the imperfection of real optical surfaces. In this paper, the mean value 0.7 is chosen for simplicity and then $P_{TX}$ is about $66.78\ mW$. Since $P_{LO}$ and $P_{TX}$ have been fixed, the split ratio of the beamsplitter BS1 in Fig. 1 satisfies $R:T\approx 1:40$ and $\eta _{RX}\approx (40/41)\times (1/8) \approx 0.122$ is the transmission efficiency of the RX beam from the receiver aperture to one segment of the photodiodes. It should be noted that because the absolute value of contributions of tiny slits between active segments to the CNR tends to be less than 1 dB-Hz and is same for different optical design parameters, it’s omitted in our calculations for convenience. Besides, because lenses of $12.7\ mm$ radius is used most frequently in laboratories, $r_{TAP}$ is set to be $12.7\ mm$. For the sake of simplicity, $\alpha _{RX}$ is set equal to $\alpha _{TX}$. Finally, there remains four variants: $\omega _{LO}, m_{TX}, r_{RAP}$ and $\alpha _{TX}$.

In the first step, we concentrate on the high case: $\alpha _{TX}=\alpha _{RX}=0$, which means that the directions of local and remote satellites coincide perfectly with the line-of-sight. Apart from that, four typical values {1.0, 1.5, 2.0, 2.5 mm} are chosen for the waist radius $\omega _{LO}$ of the LO beam. The two parameters $m_{TX}$ and $r_{RAP}$ are scanned in the ranges $[1,14]$ and $[1,14]mm$ respectively. Fig. 4 displays the distribution of the $C/N_0$ on the plane of $m_{TX}$ and $r_{RAP}$. It can be seen that the $C/N_0$ of almost all points in the diagram excesses 70 dB-Hz and the $C/N_0$ of partial regions is even above 115 dB-Hz when $\omega _{LO}$ is equal to $1.0\ mm$. As the waist radius of the LO beam increases, the maximum of the $C/N_0$ drops to about 105 dB-Hz. Besides, $m_{TX}$ that corresponds to the maximums of the $C/N_0$ decreases with $\omega _{LO}$ because of the clipping by the transmitting telescope with the fixed aperture radius of $12.7\ mm$.

 figure: Fig. 4.

Fig. 4. The distribution diagram of the $C/N_0$ on the plane of $m_{TX}$ and $r_{RAP}$. The four panels correspond to four discrete values of $\omega _{LO}$: 1.0, 1.5, 2.0 and 2.5 mm. Different colors(from blue to red) of points in the diagram denote the different amounts of the CNR.

Download Full Size | PDF

It’s time to discuss about the effect of misalignment angles of the local and remote satellites with respect to the line-of-sight on the $C/N_0$. Since the $1\ mm$ waist radius of the LO beam corresponds to the higher CNR and it’s convenient to use commercially available fiber injectors for future experiments, $\omega _{LO}$ is set to be $1.0\ mm$. Similarly, there are four choices for $\alpha _{RX}(=\alpha _{TX})$: $20,40,60,80\ \mu rad$ and the scanning ranges of $m_{TX}$ and $r_{RAP}$ are $[1,14]$ and $[1,14]mm$ as well. For convenience, we define

$$x=\frac{\omega_{LO}\cdot m_{TX}}{r_{TAP}};\ y=\frac{\omega_{LO}\cdot m_{TX}}{r_{RAP}}$$
and
$$F(\alpha_{TX})=e^{\frac{-2(\pi \omega_{LO}\cdot m_{TX}\cdot \alpha_{TX})^2}{(\lambda\cdot \psi_1(x))^2}}; G(\alpha_{RX})=e^{\frac{-2(\pi \omega_{LO}\cdot m_{TX}\cdot \alpha_{RX})^2}{(\lambda\cdot \psi_2(y))^2}}.$$

As mentioned earlier, the range of x and y should be $[0.2,2.0]$, where the fitting functions $\psi _1(x)$ and $\psi _2(y)$ are valid. In addition, the parameter space with $G(\alpha _{RX}) <0.13$ or $F(\alpha _{TX}) <0.13$ is also excluded. The scanning results of the $C/N_0$ are shown in Fig. 5. To make it look clearer, the ranges of coordinates for the X-axis of the two lower panels are changed to $[1,8]$. Obviously, the surviving regions are getting smaller due to the restriction of above-mentioned conditions and the maximum of the $C/N_0$ also decreases by about 20 dB-Hz as the angles $\alpha _{RX}$ and $\alpha _{TX}$ become larger. Based on this diagram, we can select $m_{TX}=4$ and $r_{RAP}=10\ mm$ for our optical bench design, meanwhile, the corresponding CNR for a single segment of the photoreceiver is close to 92 dB-Hz. Same as GRACE Follow-On mission, the radius of the photodiodes $r_{PD}$ is chosen to be $0.5\ mm$. Therefore, the transversal magnifications of the receiving and local telescopes are also determined: $m_{RX}=0.05,\ m_{LO}=0.2$.

 figure: Fig. 5.

Fig. 5. The distribution diagram of the $C/N_0$ on the plane of $m_{TX}$ and $r_{RAP}$. The four panels correspond to four discrete values of $\alpha _{TX}$: 20, 40, 60 and 80 $\mu rad$. Different colors(from blue to red) of points in the diagram denote the different amounts of the CNR. White areas with $x,y<0.2,\ x,y>2.0,\ F(\alpha _{TX})<0.13$ or $G(\alpha _{RX})<0.13$ have been excluded. The black point at (4,10) in the diagram denotes the parameter that is selected for the final design. Note that the ranges of coordinates for the X-axis of the two lower panels are changed to $[1,8]$ for clarity.

Download Full Size | PDF

The Table 2 shows the comparison of key parameters and the CNR of our optical design and the optical design of GRACE Follow-On mission [2,3]. It is obvious that the $C/N_0$ of our design is considerably larger even for the worst case mainly because of applications of a series of telescopes, shorter satellite separation and larger laser power. Finally, it comes to the conclusion that the optical bench design in this paper has more significant advantages than the optical design of GRACE Follow-On mission and it is a hopeful candidate for China’s future gravity missions.

Tables Icon

Table 2. The comparison of key parameters and the CNR of the optical bench design in this paper and the optical bench design in GRACE Follow-On mission.

5. Conclusion

In this work, we come up with a novel off-axis optical bench design which absorbs the virtues of the off-axis optical bench design of GRACE Follow-On mission and other on-axis optical designs. The detailed constructions of the design and the principles that the TX and RX beam are kept anti-parallel with DWS feedback loop and the TTL coupling noise is suppressed by introducing multiple sets of lens systems are described in detail. The expression of the carrier-to-noise ratio is derived and key parameters of the optical bench design are also ascertained. Finally, we confirm that the carrier-to-noise ratio for a single channel of the photoreceiver is larger than that of the optical bench design of GRACE Follow-On mission. The numerical simulations and experimental demonstrations of the optical bench design is still ongoing. The China’s future gravity missions can choose it as the optical bench with high CNR allowing for a robust phase tracking.

Funding

National Key Research and Development Program of China (2020YFC2200404, 2021YFC220144); CASC Youth Talent Support Program (201981); CAST Fund for Excellent Young Scholars (LXL2022).

Acknowledgments

The present work benefits immensely from many discussions with the AEI group and Dr.Vitali Müller in particular. Part of the work of Yichao Yang was done while he was visiting AEI through the LEGEND program of the MPG-CAS collaboration scheme in gravitational physics. Participation in the CNSA-ESA Bilateral Joint Working Group on geo-gravity Mission, MOST-EU-ESA Dialogue Forum on Space Sci Tech Cooperation is also constructive to our work.

Disclosures

The authors declare that we do not have any commercial or associative interest that represents a conflict of interest in connection with the work submitted.

Data availability

The data that support the findings of this study are available from the corresponding author upon reasonable request.

References

1. B. D. Tapley, S. Bettadpur, M. Watkins, and C. Reigber, “The gravity recovery and climate experiment: Mission overview and early results,” Geophys. Res. Lett. 31 (2004).

2. R. P. Kornfeld, B. W. Arnold, M. A. Gross, N. T. Dahya, W. M. Klipstein, P. F. Gath, and S. Bettadpur, “Grace-fo: The gravity recovery and climate experiment follow-on mission,” J. Spacecr. Rockets 56(3), 931–951 (2019). [CrossRef]  

3. V. Müller, “Design considerations for future geodesy missions and for space laser interferometry,” Ph.D. thesis, Gottfried Wilhelm Leibniz Universität Hannover, Hannover, Germany (2017).

4. O. Mandel, A. Sell, M. Chwalla, T. Schuldt, J. Krauser, D. Weise, and C. Braxmaier, “Architecture and performance analysis of an optical metrology terminal for satellite-to-satellite laser ranging,” Appl. Opt. 59(3), 653–661 (2020). [CrossRef]  

5. Y. Yang, K. Yamamoto, M. Dovale Álvarez, D. Wei, J. J. Esteban Delgado, V. Müller, J. Jia, and G. Heinzel, “On-axis optical bench for laser ranging instruments in future gravity missions,” Sensors 22(5), 2070 (2022). [CrossRef]  

6. P. Amaro-Seoane, H. Audley, S. Babak, et al., “Laser interferometer space antenna,” arXiv, arXiv:1702.00786 (2017). [CrossRef]  

7. Z. Luo, Z. Guo, G. Jin, Y. Wu, and W. Hu, “A brief analysis to taiji: Science and technology,” Results Phys. 16, 102918 (2020). [CrossRef]  

8. J. Luo, L.-S. Chen, H.-Z. Duan, et al., “Tianqin: a space-borne gravitational wave detector,” Classical Quantum Gravity 33(3), 035010 (2016). [CrossRef]  

9. R. B. Langley, “Gps receiver system noise,” GPS world 8, 40–45 (1997).

10. J. Kim, “Simulation study of a low-low satellite-to-satellite tracking mission,” Ph.D. thesis, The University of Texas at Austin (2000).

11. F. G. Cervantes, J. Livas, R. Silverberg, E. Buchanan, and R. Stebbins, “Characterization of photoreceivers for lisa,” Classical Quantum Gravity 28(9), 094010 (2011). [CrossRef]  

12. e.motion2 Team, “e.motion square: Earth system mass transport mission (square) - concept for a next generation gravity field mission,” (2014).https://dgk.badw.de/fileadmin/user_upload/Files/DGK/docs/b-318.pdf.

Data availability

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (5)

Fig. 1.
Fig. 1. The diagrammatic drawing of the optical bench design for future gravity missions.
Fig. 2.
Fig. 2. $C/N_0$ in dB-Hz as a function of the LO beam power for a single channel of the photoreceiver.
Fig. 3.
Fig. 3. The graph of the heterodyne efficiency $\eta _{CPD}$ versus the ratio $\omega _{LO,PD}/r_{PD}$.
Fig. 4.
Fig. 4. The distribution diagram of the $C/N_0$ on the plane of $m_{TX}$ and $r_{RAP}$. The four panels correspond to four discrete values of $\omega _{LO}$: 1.0, 1.5, 2.0 and 2.5 mm. Different colors(from blue to red) of points in the diagram denote the different amounts of the CNR.
Fig. 5.
Fig. 5. The distribution diagram of the $C/N_0$ on the plane of $m_{TX}$ and $r_{RAP}$. The four panels correspond to four discrete values of $\alpha _{TX}$: 20, 40, 60 and 80 $\mu rad$. Different colors(from blue to red) of points in the diagram denote the different amounts of the CNR. White areas with $x,y<0.2,\ x,y>2.0,\ F(\alpha _{TX})<0.13$ or $G(\alpha _{RX})<0.13$ have been excluded. The black point at (4,10) in the diagram denotes the parameter that is selected for the final design. Note that the ranges of coordinates for the X-axis of the two lower panels are changed to $[1,8]$ for clarity.

Tables (2)

Tables Icon

Table 1. The coefficients of the polynomial functions ψ 1 ( x ) = i = 0 6 p i x i and ψ 2 ( x ) = i = 0 6 q i x i . [3]

Tables Icon

Table 2. The comparison of key parameters and the CNR of the optical bench design in this paper and the optical bench design in GRACE Follow-On mission.

Equations (32)

Equations on this page are rendered with MathJax. Learn more.

| θ R X , o u t | = | θ R X , i n | / m R X ,
| θ T X , o u t | = | θ T X , i n | / m T X ,
| θ L O , o u t | = | θ L O , i n | / m L O .
θ T X , o u t = θ R X , i n
θ T X , i n = θ L O , i n .
m T X = | θ T X , i n | | θ T X , o u t | = | θ L O , i n | | θ R X , i n | = | θ R X , o u t | | θ R X , i n | | θ L O , i n | | θ L O , o u t | = m L O m R X .
P S D [ δ ϕ P M ] = 1   r a d 2 C / N 0 .
C / N 0 = c a r r i e r   R M S   p o w e r n o i s e   P S D = I c a r r i e r 2 P S D [ I n o i s e ] ,
I c a r r i e r 2 = 2 η P R X P L O η P D 2 ,
P S D [ I n o i s e ] = P S D [ I P R ] + P S D [ I S N ] + P S D [ I R I N ] .
P S D [ I P R ] = 25 p A 2 / H z ,   4 M H z < f < 20 M H z
P S D [ I S N ] = 2 e η P D P L O ,
P S D [ I R I N ] = P S D [ δ P L O ( f ) / P L O ] ( η P D P L O ) 2 .
C / N 0 = 2 η P R X P L O η P D 2 25 p A 2 / H z + 2 e η P D P L O + 1.0 × 10 15 / H z ( η P D P L O ) 2 .
C / N 0 = f ( P L O ) P R X η 1 p W ,
f ( P L O ) = 2 P L O 1 p W η P D 2 25 p A 2 / H z + 2 e η P D P L O + 1.0 × 10 15 / H z ( η P D P L O ) 2 .
C / N 0 = 63.5 + 10 l o g 10 ( P R X η 1 p W )   dB-Hz .
I R X 2 P T X π ω ( L ) 2 e 2 α T X 2 L 2 / ω ( L ) 2 ,
S = π r R A P 2 cos β R X π r R A P 2 ,
P R A P = 2 P T X r R A P 2 L 2 θ T X 2 e 2 α T X 2 / θ T X 2 ,
P R A P 2 P T X r R A P 2 L 2 θ T X 2 ( 1 e r T A P 2 / ( ω L O m T X ) 2 ) 2 e 2 α T X 2 / ( θ T X ψ 1 ( ω L O m T X / r T A P ) ) 2 .
θ T X = λ π ω L O m T X ,
P R X 2 π 2 η R X ω L O 2 m T X 2 P T X r R A P 2 L 2 λ 2 ( 1 e r T A P 2 / ( ω L O m T X ) 2 ) 2 e 2 ( π α T X ω L O m T X ) 2 / ( λ ψ 1 ( ω L O m T X / r T A P ) ) 2 .
η C P D = 8 | 0 r P D r e r 2 / ω L O , P D 2 J 0 ( 2 π r / λ α R X , L O ) d r | 2 r P D 2 ω L O , P D 2 ( 1 e 2 r P D 2 / ω L O , P D 2 ) ,
η C P D = 2 ω L O , P D 2 r P D 2 tanh ( r P D 2 2 ω L O , P D 2 ) .
η η C P D e 2 α R X , L O 2 / ( θ L O ψ 2 ( ω L O , P D / r P D ) / m L O ) 2 ,
ω L O , P D = ω L O m L O , r P D = r R A P m R X , α R X , L O = α R X m R X , θ L O = λ π ω L O , m T X = m L O m R X .
η = 2 ( ω L O m T X ) 2 r R A P 2 tanh ( r R A P 2 2 ( ω L O m T X ) 2 ) e 2 ( π ω L O α R X m T X ) 2 ( λ ψ 2 ( ω L O m T X / r R A P ) ) 2
λ ; L ; ω L O ; m T X ; P T X ; η R X ; r R A P ; r T A P ; α R X ; α T X .
P T X = κ ( P l a s e r P f r e q 8 P L O ) ,
x = ω L O m T X r T A P ;   y = ω L O m T X r R A P
F ( α T X ) = e 2 ( π ω L O m T X α T X ) 2 ( λ ψ 1 ( x ) ) 2 ; G ( α R X ) = e 2 ( π ω L O m T X α R X ) 2 ( λ ψ 2 ( y ) ) 2 .
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.