Double Michelson/FabryPerot interferometer for laser and displacementnoisefree gravitationalwave detection
Abstract
In this paper we demonstrate that a double Michelson interferometer with FabryPerot cavities in its arms is able to perform laser and displacementnoisefree gravitationalwave (GW) detection if certain model assumptions are met. Assuming the input mirrors of a single Michelson/FabryPerot interferometer can be rigidly attached to beamsplitter on a central platform one can manipulate with interferometer’s response signals in a way to cancel laser noise and displacement noise of all test masses except the cental platform. A pair of symmetrically positioned Michelson/FabryPerot interferometers with common central platform can be made insusceptible to the later then, thus allowing complete laser and displacementnoisefree interferometry (DFI). It is demonstrated that the DFI response to GWs of the proposed interferometer is proportional to , where is the cavity halfbandwidth, that is the strongest DFI response allowed by general relativity.
pacs:
04.30.Nk, 04.80.Nn, 07.60.Ly, 95.55.YmI Introduction
It is widely known that the firstgeneration laser interferometric GW detectors suffer from the great amount of noises of various nature. The major limiting factors at low and middle frequencies are seismic noise and thermal noise which can be referred to the class of displacement noise of the test masses. At high frequencies photon shot noise is dominant. In the standard quantum limited (SQL) V.B. Braginsky and F.Ya. Khalili (1992) secondgeneration detectors, being under preparation, the cause of SQL is the fluctuating force of radiation pressure in the laser beam (backaction noise) pushing the interferometer mirrors in a random manner. Thus standard quantum limitation also aries due to displacement noise.
Each method of suppressing or eliminating of displacement noise proposed uptodate is only suited for control of only one kind of noise. For instance, active antiseismic isolation will definitely suppress seismic noise but is helpless against thermal noise or quantum radiation pressure. From the other hand, quantumnondemolition (QND) schemes of measurements C.M. Caves (1981); W.K. Unruh (1982); V.B. Braginsky and F.Ya. Khalili (1996) are able canceling backaction noise but are certainly not suited for dealing with seismic or thermal noise.
However, recently there has been proposed a revolutionary new method of displacement noise cancelation which simultaneously eliminates the information about all external fluctuating forces but leaves certain amount of information about the gravitational waves. The major idea is to construct such an interferometer that would respond differently to the motion of the test masses and the GWs. Then the proper linear combination of the interferometer responses will cancel the fluctuations of the test masses leaving a nonvanishing information about the GWs. One may find at least two different methods proposed uptodate.
The first one, described in a series of papers by S. Kawamura et al. S. Kawamura and Y. Chen (2004); Y. Chen and S. Kawamura (2006); Y. Chen et al. (2006), bears on the distributed nature of the GWs. This can be best explained from the viewpoint of a local observer (or the local Lorentz gauge). In such a reference frame interaction of the GW with a laser interferometer adds up to two effects M. Rakhmanov (2005); S.P. Tarabrin (2007). The first one is the motion of the test masses in the GW tidal forcefield. In this aspect GWs are indistinguishable from any nonGW forces since both are sensed by the light wave only at the moments of reflection from the test masses. We will refer to this as the localized nature of the forces acting on the test masses. If the linear scale of a GW detector is much smaller than the gravitational wavelength (the socalled longwave approximation) then the effect of the GW forcefield is of the order of , where is the absolute value of the GW amplitude. Relative motion of the test masses, separated by a distance , in any force field cannot be sensed by one of them faster than , thus resulting in the rise of terms of the order of describing time delays which take the light wave to travel between the masses. Second, GW directly couples to the light wave effectively changing the coordinate velocity of light, thus manifesting itself as an effective refraction index. Light wave traveling in such a (boundless) medium acquires the information about the GW in its phase gradually, therefore it is a distributed effect having the order in longwave approximation. Ultimately, from the viewpoint of local observer displacementnoisefree interferometry implies the cancelation of localized effects (GW and nonGW forces) leaving a nonvanishing information about the distributed effect (the direct coupling of the GW to light).
It was pointed in Refs. Y. Chen and S. Kawamura (2006); Y. Chen et al. (2006) that in order the GW detector to be a truly displacementnoisefree interferometer it should be also free from optical laser noise since the latter is indistinguishable from laser displacement noise. Their sum is usually called laser phase noise. Cancelation of laser phase noise in interferometric experiments is usually achieved by implementing the differential (balanced) schemes of measurements: in conventional interferometers (such as LIGO) it is the Michelson topology and in DFIs proposed in Ref. Y. Chen et al. (2006) it is the MachZehnder (MZ) topology.
Although DFI detectors that bear on distributed nature of GWs allow complete elimination of displacement noise, the “payment” for such a gain is the significant weakening of the GW response at low frequencies. This fact directly follows from the mechanism of noise cancelation: together with displacement noise we also cancel GW terms of the and orders, leaving only the order terms. In addition, as demonstrated in Ref. Y. Chen et al. (2006), only 3D (spacebased) configurations of DFIs allow the 2ndorder response, while the response of 2D (groundbased) detectors Y. Chen and S. Kawamura (2006) is even weaker — it is proportional to . For the groundbased detectors with m and m such a dramatic decrease of susceptibility to GWs () levels all the advantages of DFI in comparison with traditional detectors dominated by displacement noise.
The second type of DFI detectors proposed in Ref. S.P. Tarabrin and S.P. Vyatchanin (2008), however, is free from this drawback since the mechanism of noise cancelation does not utilize the notion of GW distributed nature. The basic element of the proposed toy model is a single FabryPerot (FP) cavity pumped through both its mirrors with orthogonally polarized light waves. Consider first the pump through one of the mirrors. Two responses will correspond to the input wave — the reflected wave and the transmitted one. The latter is proportional to the distance between the mirrors which may vary in time under the influence of displacement noise and the GW. The reflected wave, however, is somewhat different: besides the term proportional to the variable length of the cavity it also includes the component due to the prompt reflection of the incident wave from the input mirror. Thus, reflected and transmitted waves carry slightly different information about the GW and the motion of the test masses. Combining them properly one is able to exclude displacement noise of one of the cavity mirrors. In a similar way one can deal with the second pump through another mirror. Ultimately, the proper linear combination of all four responses allows elimination of displacement noise of both mirrors leaving a nonvanishing GW signal proportional to in longwave approximation. Unfortunately, since the mechanism of noise cancelation bears on the effect of prompt reflection which is nonresonant, the obtained DFI response is not amplified with the usual resonant factor associated with the accumulation of lowfrequency signal inside the cavity.
It may seem that the elimination of mirrors displacement noise at the order contradicts the forcelike behavior of the GWs at low frequencies. However, it has been demonstrated in Ref. S.P. Tarabrin and S.P. Vyatchanin (2008) that the loss of the resonant gain results in the sensitivity limitation by displacement noises of auxiliary optical elements such as lasers and detectors. In conventional nonDFI detectors these noises are negligible since they are finesse times smaller as compared to displacement noise of the mirrors. In a doublepumped FP cavity these residual fluctuations mean that the GWinduced displacements cannot be measured absolutely but only with respect to some reference test masses, in full agreement with the relativity principle. The fluctuations of these reference masses will ultimately limit the accuracy of coordinate measurements. Therefore, in a doublepumped FP cavity DFI response is not strictly speaking displacementnoisefree: we are able to cancel noise of the mirrors (the major “headache” in conventional detectors) but not the one of auxiliary optics (lasers, detectors, beamsplitters, etc). Another drawback of the doublepumped FP DFI is the inability to cancel laser noise since the singlecavityscheme is not balanced. Although one may propose several FPbased balanced optical setups, for instance, modifications of conventional Michelson/FabryPerot topology, in all such schemes DFI response will be still limited by residual displacement noise of some test mass(es). This is a fundamental restriction of the model in Ref. S.P. Tarabrin and S.P. Vyatchanin (2008).
In spite of the undoubtful theoretical “beauty” of DFI idea proposed by S. Kawamura et al., such a detector, if realized in a fullscale experiment, will surely face certain technical difficulties as pointed out in Refs. Y. Chen and S. Kawamura (2006); Y. Chen et al. (2006). For instance, it is hardly possible to have ideal 50/50 beamsplitters or lossless mirrors. Ultimately, one may expect fold reduction of displacement noise. In other words, although completely displacementnoisefree in theory, DFIs of the 1st type will suffer from various technical limitations in practice that will eventually leave a certain level of residual displacement noise. From this viewpoint the theoretical (fundamental) limitations of the 2nd type (incomplete) DFI schemes may not seem so dramatic as one could expect. Eventually, in both cases an experimentalist would have to deal with certain technical challenges and the only thing that matters is the total level of noise left, regardless its nature. And though suppression of residual displacement noise in DFIs of the 2nd type may turn a more complex problem than that in DFIs of the 1st type, the former still have a significant advantage — strong GW response at low frequencies.
Therefore, it is reasonable to consider such GW detectors that would have as simple topologies as possible simultaneously allowing as complete displacement noise cancelation as possible and having the strongest response to GWs that general relativity allows. The latter requirement can be derived straightforwardly. According to relativity principle absolute coordinate and velocity measurements are prohibited. This is another formulation of indistinguishability between the GW and nonGW forces at the 0th and 1st orders of . The difference rises only in the order at low frequencies, as described above. If and then such a displacementnoisefree GW response will be proportional to in spectral domain. In addition, this response can be further amplified with a FP cavity, for instance. Then one should expect the rise of the resonant multiplier, where is the cavity halfbandwidth. Ultimately, the strongest DFI GW response allowed by the first principles should be proportional to in spectral domain.
The optical setup satisfying the requirements of practical reasonableness and maximum completeness of displacement noise elimination, also restricted by condition of the strongest response derived above, does not, however, immediately follows from some basic principles. It is a matter of search at large, limited by practically reasonable configurations and assumptions. For instance, in this paper introducing several model assumptions we propose a pair of symmetrically positioned Michelson interferometers with FabryPerot cavities inserted into each arm of both interferometers, as a DFI GW detector with a reasonably simple optical setup and the strongest possible response.
First, consider a conventional LIGO topology (without power and signalrecycling mirrors). Let the end mirrors be partially transmittible. In this case an interferometer will produce three response signals: the reflected (lasernoisefree) one in conventional dark park and the transmitted ones in the arms. It is worth noting here that certain care is required when calculating the responses: each response signal should be evaluated in the proper reference frame of the detector that detects the corresponding signal S.P. Tarabrin and A.A. Seleznyov (2008). Otherwise, unmeasurable quantities may arise in the analysis. An experimentalist is able to measure the quadrature components of the interferometer responses and record them for further processing. From the set of transmitted signals quadratures it is possible, in principle, to construct a lasernoisefree linear combination. Therefore, at this stage we may obtain two signals (quadratures) free from laser phase noise: the reflected one and the combined transmitted one.
Due to the sophisticated frequency dependence of the FP cavities responses these two signals can be combined in turn to eliminate one of four differential mechanical degrees of freedoms associated with the test masses (beamsplitter, two input mirrors, two end mirrors and two end detectors). At this stage we introduce some restrictions into the optical scheme: (i) end detectors and end mirrors, and (ii) input mirrors and beamsplitter are assumed to be rigidly connected. The practical legitimacy of these assumptions seems questionable and is open for criticism, although, no basic principles forbid such a gedanken (thought) experiment. Under these restrictions we are left with only two differential degrees of freedom, one of which can be eliminated in a combination of two lasernoisefree signals. We choose to cancel displacement noise associated with the differential motion of the end platforms.
Ultimately, due to the symmetry of plane GW wavefront we are able to cancel the fluctuations of the central platform (with beamsplitter and input mirrors) if the similar interferometer is positioned symmetrically (see Fig. 3 in the text) and both interferometers have common central platform (this is another gedankenexperimentsupposion). In this case the singleinterferometer partial DFI responses will have GW term of the same sign but the fluctuations of the central platform will enter with different signs. Adding two singleinterferometer responses we cancel displacement noise of the latter. Then the obtained laser and displacementnoisefree GW response signal turns out to be proportional to amplified with the cavity resonant gain .
Ii Spacetime of accelerated observer in the gravitationalwave field
Let us first remind briefly the “tools” necessary for our further considerations.
As explained in Ref. S.P. Tarabrin and A.A. Seleznyov (2008), to obtain physically reasonable, i.e. measurable, quantities, calculations should be performed in the proper reference frames of the devices that produce the corresponding experimentally observed quantities. Since they are usually subjected to external fluctuative forces, they commit random motions and thus we have to deal with their proper noninertial reference frames. Corresponding tools for solving certain boundary electrodynamical problems in such reference frames have been developed in Ref. S.P. Tarabrin and A.A. Seleznyov (2008) and utilized in Ref. S.P. Tarabrin and S.P. Vyatchanin (2008). Therefore, in this paper we will not retell the content of these works in detail but will write several useful formulas in this section.
In particular, the spacetime of an observer having nongeodesic acceleration along the axis and falling in the weak, plane, ’+’polarized GW propagating along the axis takes the following form:
(1) 
Conditions of linearized theory and are assumed to be satisfied for all reasonable and . Without the loss of generality we may assume when considering onedimensional problems.
Consider a test mass which in a state of rest (no fluctuations and no GW) has the coordinate with respect to the observer (also in a state of rest). If the test mass is subjected to some external fluactuative force which moves it according to the motion law as seen from the laboratory frame (for instance, the one associated with the Earth surface), then its motion law with respect to the observer in spacetime (1) is:
(2a)  
(2b) 
where is the observer’s law of motion as seen from the same laboratory frame. It is assumed here that . In fact, Eq. (2b) is the coordinate transformation from laboratory frame to the observer’s frame. In spectral domain which will be widely used
It is also important to take into account the effects imposed by the GW and acceleration fields on the electromagnetic waves propagating in spacetime (1). It has been derived in Refs. S.P. Tarabrin (2007); S.P. Tarabrin and A.A. Seleznyov (2008) that the waves propagating in the positive and negative directions of the axis can be described by the following vector potentials:
(3) 
where
and
Both and describe the distributed effects: is responsible for the direct coupling between the GW and the electromagnetic wave and describe the redshift imposed on the electromagnetic wave by the noninertiality of the reference frame. Weak fields describe electromagnetic fluctuations (classical or quantum).
In this paper we will also deal with the motions along the axis. In this case all the formulas remain the same but the GW function should be taken with the opposite sign (this follows from the metric (1)).
Iii Responses of a FabryPerot cavity to the gravitational wave
Let us consider the operation of a single FabryPerot cavity as a GW detector (see Fig. 1).
Cavity is assembled of two movable mirrors and , both lossless and having the amplitude transmission coefficient , . We put distance between the mirrors in the absence of the gravitational wave and optical radiation to be equal to . The incident GW is assumed to be weak, plane, ’+’polarized and propagating along the axis. Then without the loss of generality we assume the cavity to be lying in the plane along one of the GW principal axes, coinciding with the axis.
Cavity is pumped by laser L whose center of mass commits a fluctuative motion along the axis as seen from the laboratory frame. Both mirrors and have associated displacement noise and . Finally, both the detectors and that measure the reflected wave and transmitted wave correspondingly fluctuate as and .
To evaluate the response signals of the cavity we should perform the calculations in the proper reference frames of detectors and , as pointed above: the reflected signal is measured with the first one and the transmitted signal is measured with the latter one.
iii.1 Evaluation of the reflected signal
Here we will derive the expression for the wave reflected from the cavity. Since it is detected by we will perform the calculation its proper reference frame. This section mostly repeats the similar considerations in Ref. S.P. Tarabrin and S.P. Vyatchanin (2008), therefore we will proceed without detailed comments.
The origin of the coordinate system is assumed to be set up at the center of mass of . Therefore, according to Eqs. (2a,2b) test masses of the system will have the following motion laws with respect to :
In last two equations we neglected the small distance (compared to the cavity length) between the optical bench where laser and detector are located and the input mirror.
Let the cavity be pumped by laser L through the input mirror with the input wave
(4) 
Strictly speaking, the argument of here should depend on itself like , but since is already the quantity of the 1st order of smallness we can neglect such dependence. The vacuumstate pump through mirror can be written as:
(5) 
Here is the “weak” field describing optical laser noise of the pump wave and is the “weak” field describing vacuum noise in the opposite input port. Remind, that both the laser and mirror are located at , where , thus input wave does not acquire distributed phase shift when it reaches mirror .
It is convenient to represent the optical field inside the cavity as a sum of two waves, and , running in the opposite directions:
(6) 
Here describe the phase shift accumulated by the light wave while circulating inside the cavity.
Output wave reflected from the cavity is:
(7) 
If detector is a quadratic amplitude detector then it measures the quantity proportional to (neglecting very small terms of the order of ). If detector is a balanced homodyne detector then it measures the quadratures of . In this case the reference oscillation can be produced by laser L. We will call the reflected signal below.
To obtain the reflected signal we substitute fields (4 – 7) into the set of boundary conditions (conditions of the electric field continuity along the surfaces of the mirrors) S.P. Tarabrin (2007); S.P. Tarabrin and S.P. Vyatchanin (2008):
and solve the system with the method of successive approximations (see Ref. S.P. Tarabrin and S.P. Vyatchanin (2008)). The required solution of the 1st order in spectral domain is:
(8) 
The following notations have been introduced:
having the following physical meaning: describes the resonant amplification of the input amplitude inside the cavity, describes the frequencydependent resonant amplification of the variation of the circulating light wave, and are the generalized coefficients of reflection (from a FP cavity) and transmission (through a FP cavity), is the mean amplitude of the optical wave inside the cavity running in the negative direction of the axis, is the mean amplitude of the wave reflected from the cavity and is the response to GW after a single round trip of light inside the cavity
It is also very useful to analyze the physical meaning of each summand in formula (8):

. This term states that the optical laser noise is indistinguishable from laser displacement noise so both always come together and their sum is usually called laser phase noise. In spectral domain reflected wave obviously contains laser phase noise multiplied by the generalized coefficient of reflection .

is the vacuum noise from the opposite input port which is transmitted through the cavity and comes with the corresponding coefficient of transmission in the reflected signal.

describes the part of reflected wave which flows out of the cavity and, obviously, represents the phase shift accumulated by optical wave circulating inside the cavity. This phase shift is proportional to the change of cavity length due to fluctuative motion of the mirrors and the change of optical path due to GW . Accumulation of the signal is described by the resonant factor .

is the phase shift due to the prompt reflection from the input mirror.

accounts for the motion of detector which receives the reflected signal. It is straightforward to verify that this summand will be canceled in the amplitude detection, evidently meaning that the amplitude detector is insusceptible to its own motion (since amplitude is measured, but not the phase).
It is worth noting that the same formula (8) one can derive without resorting to the noninertial frame by linearly combining the considerations in (i) the transversetraceless gauge for the GW response and (ii) the laboratory frame for response to the fluctuations of test masses. Such a simplification owns to the specific feature of a roundtrip scheme: the source and the receiver of optical radiation are located approximately at the same spacial position and thus their clocks could be synchronized almost perfectly.
iii.2 Evaluation of transmitted signal
Here we will derive the expression for the transmitted signal. Since the source L and the receiver are separated by a large distance, the effect of redshift due to the acceleration of detector arises which cannot be obtained solely in the laboratory frame. Therefore, we are forced to perform the calculations in the proper reference frame of detector .
We place the origin of the coordinates at the center of mass of . Then according to Eqs. (2a) and (2b) other test masses move as following:
Again we neglected the small constant distances separating laser L and input mirror , and mirror and detector .
From the viewpoint of detector laser L emits the input wave described by the following vector potential:
(9) 
This corresponds to the condition that in the immediate vicinity of the source optical wave acquires neither phase shift due to the motion of the source, nor phase shift due to both distributed effects.
The vacuum pump is:
(10) 
It is most convenient to write the wave traveling in the positive direction of the axis inside the cavity as following:
(11) 
Such a notation means that the phase is counted from . The motion of the input mirror will be taken into account in the boundary conditions.
The phase of the optical wave traveling in the negative direction will be counted from (remind, that both and vanish at ):
(12) 
Finally, the transmitted wave with the phase counted from () is:
(13) 
If detector is a quadratic amplitude detector then it measures the quantity proportional to (neglecting very small terms of the order of ). If detector is a balanced homodyne detector then it measures the quadratures of . In this case the reference oscillation should be produced by some local source. We call the transmitted signal below.
To obtain the transmitted signal we substitute fields (9 – 13) into the set of boundary conditions:
The required 1st order solution of this system in spectral domain is:
(14) 
where
Here is the forwardtrip GW response, is the mean amplitude of the wave traveling in the positive direction of the axis inside the cavity and is the mean amplitude of the transmitted wave.
Again it is useful to analyze the physical meaning of each summand in formula (14):

describes laser noise transmitted through the cavity. The summand with means that the initial phase shift due to laser displacement is additionally redshifted with the GW and the motion of observer (detector ), because GW and acceleration fields change the rate of laser clock with respect to detector clock. It should be also noted that even the amplitude detector will be susceptible to ( in time domain), since phase modulation is transformed into amplitude modulation in a FP cavity.

is the vacuum noise reflected from the cavity into the transmitted port.

describes the total variation of the phase accumulated inside the cavity.

accounts for the displacement of the receiver. If detector is the amplitude one, this term becomes unmeasurable.
Iv Michelson/FabryPerot topology
The major disadvantage of a single cavitybased GW detector is the significant level of laser noise which dominates over other noises in practice. To cancel laser noise one should implement a balanced optical setup, for instance, a Michelson interferometer tuned to darkport regime.
iv.1 Response signals of a Michelson/FabryPerot interferometer
Let us consider a Michelsontype interferometer with FP cavities in its arms (see Fig. 2).
Laser L which randomly moves along the axis as pumps the interferometer with the input wave . Upon arrival to 50/50beamsplitter BS which may fluctuate along the  and axes as and respectively, the input wave is splitted into two waves and which pump horizontal and vertical arms respectively. FabryPerot cavity in the horizontal arm is assembled of two mirrors and which fluctuate along the axis as and . The similar cavity in the vertical arm is assembled of mirrors and which fluctuate along the axis as and . Both the cavities may produce reflected and transmitted waves. Reflected waves and return towards the beamsplitter and interfere. Assume the interferometer is tuned to dark port regime. This means that the reflected waves interfere destructively and the mean optical power returns towards the laser. However, the weak timedependent (signal) part penetrates into the dark port and falls on detector which fluctuates along the axis as . Let the latter one operate as balanced homodyne detector with the reference oscillation produced by laser L. Dark port of also produces the vacuum pump . Transmitted waves and are measured with the corresponding detectors: in horizontal arm it is and in vertical arm it is moving randomly as and correspondingly. Both detectors may operate either as amplitude or homodyne ones. In the latter case reference oscillations should be produced by some local sources (see below). Both detector ports also produce vacuum pumps and .
Let us now write the equations for input and output waves in an interferometer. We will not write the expressions for the waves explicitly; one can obtain them straightforwardly making obvious changes in the formulas from the previous section. At beamsplitter the relation between the input waves is as following:
The relation between the reflected waves is:
At this stage we do not define the reference frame, therefore, fields and coordinates of the test masses should be specified explicitly for this or that frame.
These equations can be solved straightforwardly. However, we do not need to do this since we already know the solution for a single cavity (8). First, we need to write explicitly the expressions for weak reflected fields and . To do this we use the first pair of BS boundary conditions to obtain in spectral domain:
(15a)  
(15b)  
(15c) 
Now let us specify the reference frame. Since both reflected waves will ultimately end up at detector , it is necessary to work in its proper frame. Another way is to use the laboratory frame which implementation is justified with the roundtrip situation: laser L and detector can be approximately considered as located at the same spatial position. In any case, substitution Eqs. (15a – 15c) into formula (8) we obtain:
Here we assumed that the ’+’polarized GW is perfectly aligned along the interferometer arms. We also neglected the term proportional to since the signal wave does not include “strong” mean component. For simplicity let us assume that both the cavities have equal detunings and bandwidths. This results in and .
The boundary condition for reflected waves dictates that:
Substituting here the obtained expressions for and we obtain:
(16) 
One can note that the obtained signal is very similar to the one of a single cavity. Namely, the following degrees of freedom are equivalent from this viewpoint: , , . The latter relation means that the beamsplitter effectively cuts all laser phase noise introducing, however, its own displacement noise.
In an experiment homodyne detector measures the quadrature components of . However, keeping this in mind, we will deal with the field amplitude itself, since calculations with quadratures result in very cumbersome formulas, while not changing the physical meaning of the ultimate results.
Now let us derive the transmitted signals. Since they are detected by two different devices, each of the signals should be calculated in the proper reference frame of the corresponding detector. Keeping in mind that and should be explicitly specified for each of these reference frame, we substitute Eqs. (15a – 15c) into formula (14) and obtain:
Two different regimes of detection of transmitted signals are possible.

Resonant regime. Let both cavities be tunes to resonance. In this case amplitude detector measures