# Full text of "Search for a stochastic background of 100-MHz gravitational waves with laser interferometers"

## See other formats

Search for a stochastic background of 100-MHz gravitational waves with laser interferometers Tomotada Akutsu/'Q Seiji Kawamura, 2 Atsushi Nishizawa, 3 Koji Arai, 2 Kazuhiro Yamamoto, 2 Daisuke Tatsumi, 2 Shigeo Nagano, Erina Nishida, 5 Takeshi Chiba, 6 Ryuichi Takahashi, 7 Naoshi Sugiyama, 7 Mitsuhiro Fukushima, 2 Toshitaka Yamazaki, 2 and Masa-Katsu Fujimoto 2 O 1 Department of Astronomy, University of Tokyo, Bunkyo, Tokyo 113-0033, Japan <n ; < qq 1 Kyoto University, Kyoto 606-8501, Japan o oo O 2 National Astronomical Observatory of Japan, Mitaka, Tokyo 181-8588, Japan 3 Graduate School of Human and Environmental Studies, ^National Institute of Information and Communications Technology, Koganei, Tokyo 184-8795, Japan cr 5 Graduate School of Humanities and Sciences, Ochanomizu University, Bunkyo, Tokyo 112-8610, Japan 6 Department of Physics, College of Humanities and Sciences, > q ', Nihon University, Tokyo 156-8550, Japan 7 Graduate School of Science, Nagoya University, Nagoya 464-8602, Japan (Dated: August 8, 2008) Abstract ^ ■ This letter reports the results of a search for a stochastic background of gravitational waves (GW) at 100 MHz by laser interferometry. We have developed a GW detector, which is a pair of 75-cm baseline synchronous recycling (resonant recycling) interferometers. Each interferometer has a strain sensitivity of ~ 10~ 16 Hz -1 / 2 at 100 MHz. By cross-correlating the outputs of the two interferometers within 1000 seconds, we found /ifog^gw < 6 x 10 25 to be an upper limit on the energy density spectrum of the GW background in a 2-kHz bandwidth around 100 MHz, where a flat spectrum is assumed. PACS numbers: Valid PACS appear here 1 Recently, Cruise and Ingley reported on a detector for gravitational waves (GW) at 100 MHz l|. Their GW detector is a pair of waveguide loop cavities, each of which has a strain sensitivity of ~ 10~ 14 Hz -1 / 2 at the frequency. Except for this, no experiments were attempted to directly detect GWs at very high frequencies (above 100 kHz), while many theories predict a stochastic gravitational-wave background (GWB) in a broad range of frequencies, 10~ 18 — 10 10 Hz. At very high frequencies, a relatively large GWB is predicted by some models of the early universe and compact astronomical objects (references are summarized in our previous paper [2]). Although the amount of the cosmic GWB is indirectly limited by not only the helium-4 abundance due to big-bang nucleosynthesis [3|, but also measurements of the cosmic microwave background |4| , direct search experiments for a GWB at very high frequencies should be significant. We have developed a more sensitive detector for 100-MHz GWs using laser interferome- ters. The detector is a pair of synchronous recycling interferometers, where the synchronous recycling (or resonant recycling) technique was proposed by Drever in the 1980s |5|. In our previous papers 0, E| , we showed that this interferometer is suitable to detect a GWB at very high frequency with high signal-to-noise ratio (SNR), and that the SNR can be im- proved by cross-correlating the outputs of the two interferometers. In this letter, we report the first results of the search for a stochastic GWB at 100 MHz with the GW detector. Synchronous recycling interferometer. — The interferometer has a resonant response to GWs at a specific frequency p], IsJ]. GW signals are enhanced in a recycling cavity (see Fig. [1]), which is formed by a recycling mirror (RM), a transfer mirror (TM), and two end mirrors (EMI and EM2). The size of the recycling cavity determines the resonant frequency, where the signal enhancement is proportional to the laser power kept in the cavity. At the entrance of the interferometer, a laser beam is divided into two orthogonal directions by a beamsplitter (BS). Thus two beams are incident on the RM, which is a beamsplitter but with relatively high reflectivity. When the laser frequency is stabilized to the recycling cavity, the two beams passing through the RM are resonant in the cavity by circulating many times along a common path in opposite directions (clockwise and counterclockwise). At the same time, the two circulating beams will experience differential phase shifts due to quadrupole components of GWs. The phase difference is maximized for the GWs at the same frequency as the free- spectral range z/fsr, the inverse of the round-trip period of the circulating beams. The phase 2 Recycling cavity^ £^ End mjrror 2 Recycling mirror Laser ; Beamsplitter ~£t Photodetector FIG. 1: Schematic view of a synchronous recycling interferometer. GW signals are enhanced in the recycling cavity, and detected with the photodetector. difference is enhanced as the laser power builds up depending on the finesse of the recycling cavity. The beams that left the cavity are recombined at the BS so that the differential components (GW signals) are detected with the photodetector (PD). Experimental setup. — We have developed two synchronous recycling interferometers, hereafter called IFO-1 and IFO-2. For each interferometer (see Fig. [2]), we use a Nd:YAG continuous-wave laser with a wavelength of 1064 nm and a laser power of 0.5 W. The laser beam passes through an electro-optic phase modulator (EOl) and a Faraday isolator (FI), and then enters the interferometer. The recycling cavity is designed to have a baseline length (distance from the RM to the EMI or EM2 [15j) of L ~ 75 cm so that the GW response is maximized at z/fsr = c/(4L) ~ 100 MHz, where c is the speed of light 16]. Because this experiment is the first step in the direct detection of a GWB at 100 MHz, both interferometers are constructed in the air, and each recycling cavity is designed to have a finesse of ~ 100; each RM has relatively low reflectivity (nominal 98.5%). For calibration, we use the E02 to simulate GWs by modulating the phases of the circulating beams in the cavity. The size of the Sagnac interferometer, which is formed by the BS, the RM, and two steering mirrors, is relatively small (12.5-cm square optical path), and thus its GW response is insignificant compared to that of the recycling cavity. The laser frequency is stabilized to the recycling cavity by the Pound-Drever-Hall tech- nique [91]. This technique requires phase- modulation sidebands spaced by a radio frequency (RF) from the laser-source (carrier) frequency z/ in the optical frequency domain. The RF 3 PS EM2 LO: 14.696840 MHz Freq. stabilization Qdaq FIG. 2: (Color online.) Schematic view of the experimental setup of one of the interferometers. EO: electro-optic phase modulator; FI: Faraday isolator; DAQ: data acquisition system. The laser frequency is stabilized to the recycling cavity by the Pound-Drever-Hall technique. GW signal sidebands are once converted to intermediate-frequency (IF) signals at the PD2. Then the IF signals are mixed with a local oscillator (LO), and converted to audio- frequency (AF) signals. The AF signals are recorded with the DAQ. The E02 is used to simulate GW signal sidebands for calibration. sidebands at v ± / RF are induced at the EOl, where the laser light is phase- modulated & t /rf = 85.4 MHz. The PD1 detects the light reflected from the cavity and produces a photocurrent, which contains RF signals modulated by the relative deviation between the laser and the cavity. We correct the relative deviation using signals demodulated from the RF signals. The target GW signals are converted to electrical signals at intermediate frequencies (IF) ~ 15 MHz with the PD2, since it is difficult to make a low-noise photodetector that can respond to signals at very high frequencies (~ 100 MHz). The PD2 produces IF signals at /if = few — /rf in response to the beat between the RF sidebands (also used for the laser stabilization) and signal sidebands (representation of the GW signals in the optical frequency domain) at z/ ± /gw> where /gw is the GW frequency. A small fraction of the RF sidebands leaks to the PD2, since the splitting ratio of the BS is not exactly balanced in the realistic case. Otherwise, non-differential components including the RF sidebands are in 4 100 Frequency [MHz] 101 FIG. 3: (Color online.) Strain sensitivities of the interferometers estimated from the IF signals. The solid red and dashed green lines represent the sensitivities of IFO-1 and IFO-2, respectively. principle completely reflected into the PD1. The PD2 is designed to have a band-path filter centered at /if ~ 15 MHz with about 1-MHz bandwidth. Thus the conversion coefficient from GWs to IF signals contains the filter response of the PD2 as well as the frequency response of the recycling cavity. We find the strain sensitivity of each interferometer is about 10 -16 Hz -1 ' 2 around 100 MHz (Fig. [3]) [17]. We estimate the sensitivity from IF signals. For calibration of the interfer- ometer outputs, we estimate the conversion coefficient from the IF signals Vjp to the GW amplitudes h by modulating the phases of the circulating beams with the E02 driven by calibration signals V ca ,\ [l8|. The coefficient is the product of the response of the recycling cavity and the band-path filter of the PD2, and is estimated by A(f GW )C(f GW ). (1) V lF (fw) M/if) where A(rad/V) is the measured modulation efficiency of the E02, and C (strain/rad) is the calculated conversion coefficient from the phase modulation to the simulated GWs. The term C is a function of the distance of the E02 from the RM, and it is ~ 25 cm in our experiment. The IF signals vary too quickly to be sampled with an inexpensive data acquisition (DAQ) system. We convert the IF signals to recordable audio-frequency (AF) signals at /af = few - (/ c - A//2) with a local oscillator (LO) at f LO = f c - Af/2 - f RF , where we choose Af = 6.32 kHz as a signal bandwidth to be recorded, and f c = 100.1 MHz as a center frequency of the bandwidth. They yield /lo = 14.696840 MHz. For example, GWs at 100.1 MHz corresponds to AF signals at 3.16 kHz. Cross-correlation analysis. — Using the outputs of the two interferometers, we have per- formed a cross- correlation analysis to reduce uncorrelated noises between them and improve the SNR, the ratio of the GW signals to the interferometer noises. The analysis method is similar to the method used in LIGO 10 11] . We assume that a GWB is isotropic, unpolarized, stationary, and Gaussian, and it is so small that the interferometer outputs are dominated by their noises rather than GW signals. The GWB is often character- ized by a normalized energy density spectrum per unit logarithmic frequency interval fi gw (/) = p^ 1 dp(f)/d\n f, where p(f) is the cumulative energy density of GWB included below / Hz, and p c = 3HqC 2 / (8tcG) is the critical energy density of the universe; here G is the Newton constant, and Hq = h\oo x lOOkm/s/Mpc is the Hubble constant. In this letter, we also use the form hl 00 Q gw (f), which is independent of the value of hi 00 . We define a cross-correlation statistic: '12 1 I™ x* x (f)x 2 (f)Q(f)df, (2) J — oo where x\ and x 2 are Fourier components of the signal outputs from IFO-1 and IFO-2, respectively; T is the observation time period; Q is the optimal filter that optimizes the SNR of an expectation value (ensemble average) of Z\ 2 estimated from available data (the exact definition of Q will be given later in Eq.([5])). The expectation value of Z\ 2 and its variance are respectively written as ^ s (zi2) = H SZ *^^mqw> (3) 4 = (Z 2 l2 ) - (Z 12 ) 2 ~ — / dfPi(|/|)P 2 (i/|)|Q(/)i 2 , (4) J —00 where Pi and P 2 are the one-sided power spectral densities (PSD) of the noises in IFO-1 and IFO-2, respectively [h]]; and 712 is called the reduced overlap reduction function. As the usual overlap reduction function in the low-frequency limit Q] , 712 represents the reduction of the signal correlation caused by the distance between the two interferometer sites and the alignment of their arms. In our experiment, 712 ~ 0.93 is nearly constant around 100 MHz, because the two recycling cavities are co-aligned and almost co-located (the distance is ~ 10cm) j^ . 6 The SNR of the estimation is defined as Hz / y&z- To optimize the SNR, the filter Q is chosen as n(f\ - y n gw(l/l)7i2(/) (r .s Qn \f\ 3 Pi(\f\m\f\y [) where K is a normalization constant [6|. Using this optimal filter, we find that the SNR is written as snr - 3H <> s/t \ r .„ -rumuw) i 1/2 SNR -T^ VT [J_JwmlmW)\ ■ (6) Thus the SNR in principle increases proportional to VT. The observation time period used for the cross-correlation analysis is 1070.5 seconds. The data record is divided into N = 439 segments. For each segment, the cross-correlation and its uncertainty are calculated based on Eqs.© and (j4"l); we will refer to the calculated ones as Zu and a|, respectively. The ensemble average \iz defined in Eq.Q is estimated by a weighting average: N A(n) z HZ = ^' V 2 (n) ' ^ ' ' n=\ °Z where the superscript "(n)" indicates that the quantity is calculated from the n-th segment n, = 1, 2, ... , N); a is the uncertainty of fiz and is written as N 1 T27; n=l a Z -1 (8) As the integration domain in Eq.(T2]), we use a range from 2.08 kHz to 4.19 kHz for AF signals, which corresponds to a 2-kHz bandwidth around 100.1 MHz for GW signals. Because the optimal filter in Eq.(jSJ) contains fi gw (/) itself, we need to assume its spectrum in advance. We assume that the spectrum will be flat in such a narrow bandwidth. The variation of fiz itself with respect to the observation time period is shown in FigJH The two curves represent ^±1.65(7^, and the area enclosed by the curves is a two-sided 90% confidence interval of \xz- After the 1070.5-second observation, we obtain (iz = 4.9 x 10~ 12 with = 3.7 x 10~ 12 . The 90% confidence interval [fiz — 1.65<7 M ,/iz + 1.65(5"^] includes Hz = 0; in other words, there is a possibility of Q gw = at 100 MHz. Instead, we consider an upper limit on the amount of the stochastic GWB at 100 MHz. We define the upper limit as a one-sided 90% confidence level; in terms of Hz, the upper limit corresponds to fiz + 1.28a M . Then we find h^Q^ < 6 x 10 25 as an upper limit on the stochastic GWB at around 100.1 MHz from the direct search experiment. Note that 1 1 1 1 i i , , , , i i i i i i i i i i i i i i : _ , , , , , , , . i . . , , i , , , , . ... t .... ' 200 400 600 800 1000 1200 Observation time [sec] FIG. 4: (Color online.) Variation of fiz with respect to the observation time period is shown as filled red circles. The area enclosed by black curves is its two-sided 90% confidence interval. this is also an upper limit on the correlated noises between the two interferometers at this frequency. Conclusions. — We searched for a stochastic GWB at 100 MHz by laser interferometry. The GW detector is a pair of synchronous recycling interferometers. Each interferometer has a strain sensitivity of ~ 10 -16 Hz -1 ' 2 to GWs at 100 MHz. Using the two interferometers, we directly searched for a stochastic GWB centered at 100.1 MHz with 2-kHz bandwidth in 1070.5 seconds. We performed a cross-correlation analysis to improve the SNR of the search. We found h^Q^ < 6 x 10 25 to be an upper limit on the energy density of a stochastic GWB at 100 MHz. We plan to improve the GW detector by increasing the finesse of each recycling cavity up to about 4.5 x 10 4 . For this purpose, each cavity will be constructed in a vacuum with high-reflectivity mirrors in future. Then each interferometer will have a strain sensitivity of about 4.7 x 10" 21 Hz~ 1 / 2 . For about a one-year observation, we should obtain a tighter upper limit as hl 00 fl gw ~ 2.8 x 10 14 around 100 MHz by a cross-correlation analysis with these two interferometers. This research is supported by Grant-in- Aid for Scientific Research (A) 17204018 from the Ministry of Education, Culture, Sports, Science and Technology. 8 * Electronic address: takutsu@g ravity.mtk.nao.ac.jp| [1] A. M. Cruise and R. M. J. Ingley, Class. Quantum. Grav. 23, 6185 (2006). [2] A. Nishizawa et al., Phys. Rev. D 77, 022002 (2008). [3] M. Maggiore, Phys. Rep. 331, 283 (2000). [4] T. L. Smith, E. Pierpaoli, and M. Kamionkowski, Phys. Rev. Lett. 97, 021301 (2006). [5] R. W. P. Drever, Gravitational radiation, edited by N. Deruelle and T. Piran (North-Holland, Amsterdam, 1983), pp.321-338. [6] A. Nishizawa et al., arXiv:0801.4149, "Optimal Location of Two Laser-interferometric Detec- tors for Gravitational Wave Backgrounds at 100 MHz". [7] J.-Y. Vinet, B. Meers, C. N. Man, and A. Brillet, Phys. Rev. D 38, 433 (1988). [8] B. J. Meers, Phys. Rev. D 38, 2317 (1988). [9] R. W. P. Drever, J. L. Hall, F. V. Kowalski, J. Hough, G. M. Ford, and H. Ward, Appl. Phys. B 31, 97 (1983). [10] B. Allen and J. D. Romano, Phys. Rev. D 59, 102001 (1999). [11] B. Abbott et al., Astrophys. J. 659, 918 (2007). [12] N. Christensen, Phys. Rev. D 46, 5250 (1992). [13] E. E. Flanagan, Phys. Rev. D 48, 2389 (1993). [14] K. X. Sun, M. M. Fejer, E. Gustafson, and R. L. Byer, Phys. Rev. Lett. 76, 3053 (1996). [15] The cavity's optical path is crossed so that it encloses nearly zero area so as to be insensitive to the Sagnac effect [141 ]. [16] One can construct a narrowband audio-frequency GW detector by adding delay lines or Fabry- Perot cavities in the recycling cavity Q] . [17] There was a concern on environmental electromagnetic noises around 100 MHz. We use a loop antenna to measure the noise features, and choose a quiet band for the GW search. [18] The calibration signal generator induces noises as shown in Fig. [3] (the common peaks for both IFOs at about 100.1 MHz in this case), even when the generator is in a stand-by mode. We power off the generator during the 1000-second data taking described below. [19] The square root of each PSD is the strain sensitivity. [20] The overlap reduction function, j(f), defined in our previous paper (| is related to 712 in the 9 following form: -f(f) = 712 (/) sin I ^ 1J L - l ^FSR 7T / 2 ^FSR . Both functions reduce to the usual overlap reduction function in the low-frequency limit: / — > 0. 10