Gravitational waves observation : brief comments

In preceding papers we have shown the fundamental aspects of the General Relativity (GR), of the emission and detection of gravitational waves (GW). With the same objective we analyze the two recent observations of the GW done by the Laser Interferometer Gravitational-Wave Observatory (LIGO) in the United States of America. These GW observations known as the GW150914 and GW151226 events are emitted by binary-star systems of black-holes (BBH). We present the basic principles of the laser interferometric technique that today is considered as the only one able to detect with certainty the GW. Using a simple relativistic approach we explain approximately the observed GW in the spiral stage.


Introduction
It is important to begin this paper remembering that solving the Einstein's equations for a static distribution of mass [1 -4] we can obtain the unperturbed metric tensor g µν.When this mass distribution is slightly perturbed we write g µν = g

(o)
µν + h µν where h µν represents a small correction of the initial gravitational field.In absence of matter g µν is given by Minkowsky tensor g (o) µν = (1,1,1,-1).All relativistic accepted gravitation [1 -6] theories predict the existence of gravitational waves (GW) and that they are expected to have extremely small amplitudes with h µν ∼ 10 −21 .This can be seen in recent articles [5,6], in which the fundamental equations of the Einstein gravitation theory are constructed, experimental tests [5 -7] are shown, and there are deductions of the basic equations predicting the emission of gravitational waves (GW) [8].These equations have also been used to estimate the intensities, wave amplitudes (or strain) h µν of the GW generated by some transmitter systems [9,10] like, for instance,binary stars, neutron star pulsations, the precession of deformed neutron stars with oscillating quadrupoles, rotating bars and cataclysmic processes that give rise to supernovae.Many different detection techniques have been proposed to detect GW [3,11,12]: laser interferometric techniques, resonant solids, fluctuations in the distance between Earth and Moon, crust oscillations of the Earth, normal modes of vibrations in solids in the form of rectangles, forks and rings, spinning rods and spinning tubes with fluids inside.In Section 2 we present the basic principles of the Laser Interferometry 13 technique.Interferometers are the most sensible detecting systems: theoretically they are able to detect GW with amplitudes h ≈ 3 X 10 −22 , about 300 times smaller than estimates for any other detectors.In addition, another advantage of interferometers is that they can detect waves with any frequency.Only in 1990 it was arrived at the conclusion that GW could be detected with certainty only by interferometers [9, 10 -12].Detection of GW by another systems would be fortuitous.This occurred only after 20 years of intense researches with different techniques in USA, Japan, Italy, Australia, Holland,... Interferometric detectors were first suggested in the early 1960s and the 1970s [8 -12].In Section 3 we analyzed the two recent GW observations performed by the two interferometric detectors (LIGO) installed in the United States of America.These GW have been created in the merge of a binary black hole (BBH) system.As it will be seen in Section 2 the entire merge evolution of the BBH can be divided into three stages: "inspiral", "merger"(or "plunge") and "ringdown".To calculate exactly the GW emitted in the complete merge process it is necessary to solve the full equations of general relativity theory (GRT).This can be done only in numerical relativity simulations.Some simplified algebraic models can be used, for instance, in the "inspiral"stage.In Section 4, using a simple relativistic approach, we explain approximately the amplitude ("strain") of the GW signals detected by the LIGO in the "inspiral"stage.Of course, these estimations are non rigorous, only instructive.

Laser interferometric technique
The laser interferometric detector (LID) is basically a laser large Michelson [11 -13] interferometer with perpendicular arms where there are three mirrors M, M1 and M2 as shown in Figure 1.M1 and M2 mirrors are attached to suspended blocks that can swing freely like pendulums.With the aid of a highly potent monochromatic laser the interferometer measures the relative displacement of the mirrors which would be generated by GW.
Since 1990 several interferometric detectors were being built and designed: GEO600 (cooperation between Germany and England), LIGO (USA), VIRGO (collaboration between Italy and France), TAMA 300 (Japan), AIGO (Australia), LISA (spatial design, NASA, ESA) and LISC (LISA-in- ternational space project).As there is a huge amount of technical information involved in such major projects, we suggest to readers to search at Google using the words "Interferometric Detectors".The LIGO, VIRGO and GEO600 are very similar in concept.All projects have L-shaped facilities with multi-kilometer-long arms (4 km for LIGO, 3 km for VIRGO and 600 m for GEO600) with evacuated tubes that contain laser beams monitoring the positions of precision mirrors using interferometry.There is an intense collaboration between these international groups.According to Einstein's theory, the relative distance of the mirrors along the two arms changes very slightly when a GW passes by.The interferometers are set up in such a way that change in the lengths of the arms as small as ∼ 10 −18 m (a thousandth the diameter of the atomic nucleus) can be detected.
Let us see the principle of operation of these detectors.Thus (see Fig. 1) consider the GW incident on the z-axis direction (perpendicular to the plane of the figure) and the polarizing axes (+) along the x-axis (passing through the laser, M and M2) and y (passing by M and M1).Following our article [6c], let us see how this GW changes the distance between two particles located along the x axis (y = z = 0) at a point x = -dx/2, and the other at x = dx/2 .Supposing that the GW have a single polarization (+) with an amplitude or strain h the distance d (x,t) (calculated for delayed times t -xc) between the particles is given by, Cattani e Bassalo e4202-3 where k= ωc = 2π/λ and ω and λ, respectively, are the angular frequency and wavelength of the GW.Thus, if the points are separated by a finite distance with coordinates x = ± o /2 from (1) we obtain the distance (t), where o is the distance before the GW reach the particles.Thus, according to (2) the distance change ∆ (t) = (t)o generated by G is given by showing that the amplitude ∆ of the displacement is given by Under analogous conditions for the two points along the y axis, ∆ is also given by (4).Note that for o → 0 the deformation amplitude is given by ε = ∆ / o → h/2, which is the shearing amplitude [5] along the x axis or y axis, as it is expected.The distance changes ∆ of these arms have opposite phases: while one arm shortens the other lengthens and vice versa [6c].
According to (4) the maximum ∆ is obtained when o = λ/2, that is, (∆ ) max = λh/2π.This implies that a minimum amplitude h min of one GW which would be able to create ∆ is given by h min ≈ ∆ /(λ/2).Distance changes ∆ (t) create fluctuations in interference fringes.The laser light is divided into two parts by the mirror M which propagate along the interferometer arms.The light beams are reflected by M1 and M2, return, coherently recombine again in M, and are detected by the photosensitive device.In (4) the distance o would be the undisturbed "optical path length".From (4) we see that the maximum ∆ value is obtained when o = λ/2, or (∆ ) max = o h/π.For GW with f = 1 kHz, (∆ ) max is obtained when o = 150 km.Multiple passes are made in each arm to increase the optical path thereby obtaining an "effective optical length "L.A simplified scheme is shown in Fig. 1 with the beam back and forth between the mirrors by successive reflections.This process is known as "delay line"(optical delay line).Assuming h = 10 −21 and f = 1 kHz, λ = c/f =300 km, we found that (∆ ) max = h λ /2π = 5 X 10 −15 cm, i.e. 100 times less than the proton radius .
Due to the extremely small intensities of GW [5,6], the interferometer must be optically perfect and extremely well isolated from the rest of the world (seismic isolation, cosmic rays, electromagnetic fields, etc.).Lasers must be very powerful, extremely monochromatic and stable, the mirrors must have high reflectivity, the light path must be done over extreme high vacuum tubes, etc.
The measurement accuracy is mainly limited by fluctuations in the interference fringes, the number of detected photons as they mimic the effect of changes in the optical path.Suppose at a point P (in the mirror M) the intensity of the light from one arm is A and the intensity from the other arm is exp (i2π∆L/λ em ), where ∆L = L 1 -L 2 is difference of the effective paths of the two light beams along the arms 1 and 2, respectively, and λ em is the wavelength of the laser light.Thus, the light intensity at the detector is given by The distribution of the number N of photons in the interference fringes is then given by [13] Due to the incident GW the distances originally undisturbed L 1 and L 2 are changed by dL 1 and dL 2 , respectively, resulting in a total variation d(∆L) = dL 1 -dL 2 .The effect of d(∆L) on N is given by As the statistical error of N is N 1/2 , the accuracy in the d(∆L) measurement can be estimated replacing dN by N 1/2 into (7), (8) where we see that d(∆L) is minimum when sin (2π∆L/ λ em ) = 1.With this condition (6) becomes written as N = N max /2 = N o , where N o is the mean intensity.Under these conditions, (8) gives The coefficient N o depends on the laser power P, the time interval ∆τ of the duration of the measurement, and the detection efficiency ε of the photons according to equation [7,9], where ω em is the angular frequency of the laser light.Substituting ( 10) into ( 9) we have Using ( 4) we have seen that the maximum amount of displacement due to the GW is achieved when the path length was equal to λ/2.Under these conditions the photons remain a time λ/2c within the interferometer arm.This means that the measurement time ∆τ must be equal to or greater than the travel time λ/2c.So, putting ∆τ ≈ λ/2c into (11) we get According to (4) we conclude that the minimum strain h min of GW that generates a change d(∆L) in an effective length must obey the relationship h min = d(∆L)/L = d(∆L)/(λ/2) from which we obtain, using (12), Using ( 13) and P = 100 W, λ em = 500 nm, f = 1kHz and ε = 0.3, we verify that h min ≈ 310 −22 .
As in (13) we supposed that f = 1 kHz and the optical path L = λ/2 = c /2f = 150 km the light ray must make a number n = 150/L b of passages where L b is the length of the interferometer arm.In the arm the travel time of the photon ∆τ is equal to ∆τ = λ /2c = 5ms.In VIRGO, for instance, as L b = 3 km we verify that n = 50.In these conditions we see that the optical path variation d(∆L) = h min L = 4.5 10 −15 cm and the fluctuation dN in the number of photons is dN ≈ N 1/2 o .According to (9 ), it is given by dN = N Note that the interferometers can in principle detect GW with any frequency.Eq.( 13) is a simple estimation of the minimum strain h min taking into account only a few factors such: λ em , laser power P, ε detection efficiency of the light signal, etc ...However, there are many other factors that limit sensitivity of interferometers and that depend on frequency f of the GW.Very accurate estimates and simulations have been made to estimate the h min by the observatories GEO 600, LIGO, VIRGO, TAMA 300, AIGO, LISA and LISC taking into account numerous effects.These are seismic, thermal and thermo-elastic in the pendulums that hold the mirrors, radiation pressure, acoustic, electromagnetic, cosmic rays, distortion of the surfaces of the mirrors due to laser coating of mirrors, "shot noise", etc.We believe instructive, from a didactical point of view, to display in one figure how reliable effects contribute to obtain the sensitivity curve sensitivity curve h min (f)/sqrt(Hz) as function of f[Hz].This is the case of the VIRGO interferometer [10,11] that was projected [14,15] in 1990 and today is installed at Cascina (Pisa-Italy); each arm of the interferometer has 3 km.
As it will be shown in Section 3, the current Advanced LIGO and VIRGO Interferometers are able to detect signals with h min ≈ 3 X 10 −22 which is theoretically 300 times smaller than the estimates for any other GW detectors [6 -12].

Observations of Gravitational Waves by LIGO
As we have mentioned above, interferometric detectors were first suggested in the early 1960s and the 1970s.Long-baseline broadband laser interferometers began to be proposed in 1989 -2000.Combinations of these detectors made joint observations from 2002 through 2011, setting upper limits on a variety of gravitational-wave while LIGO and VIRGO became the first significantly more sensitive network of advanced detectors to begin observations [16,17].
In Figure 2   A passing gravitational wave (Figure 3) effectively alters the arm lengths such that the measured Figura 3: Schematic incidence of the plane GW on the Ligo detector.Note the tensorial nature of the incident GW given by h αβ (x,t) [4,8].(Publicly available image at https://www.ligo.caltech.edu/)difference is ∆L(t) = δL x -δL y = h(t)L, where h(t) is the gravitational-wave strain.
Figure 4 shows for the (Advanced) LIGO the strain sensitivity h min / √ Hz as a function of the gravitational frequency f in Hz.LIGO detector is able to detect signals with h min ∼ 7 X 10 −23 .The current strain sensitivity curves h min / √ Hz for VIRGO and LIGO are very similar [18].

Binary Black Hole (BBH)
For many years, proving the existence of BBHs was quite difficult [19] because of nature of the BHs themselves, and the limited means of detection available.However, in one event that a pair of BHs were to merge, an immense amount of energy should be given off as GW, with distinct waveforms that can be calculated using the GR.Therefore, during the late 20th and early 21st centuries, BBHs became of great interest scientifically as a potential source of such waves, and a means by which gravitational waves could be proven to exist.BBH mergers would be one of the strongest known sources of gravitational waves in the Universe, and thus offer a good chance of directly detecting such waves [19].These sources have been finally observed in the events known as GW150914 [16] and GW151226 [20].
The entire merge evolution of the BBH can be divided into three stages: "inspiral", "merger"(or "plunge") and "ringdown".The "inspiral"is the first stage of the BBH life which resembles a gradually shrinking orbit and takes a longer time, as the emitted GWs are very weak when the BHs are distant from each other.This phase contributes only with a small fraction of the GW energy emitted in the entire process.As the BH orbit shrinks, the speed increases , and the GW emission increases.When the BHs are close the GW, they cause the orbit to shrink rapidly.In the final fraction of a second the BH can reach extremely high velocity.This is followed by a plunging orbit and the BH will "merge"once they are close enough.At this time the GW amplitude reaches its peak (it is only clearly observed in the GW151226 event [20]).Once merged, the single hole settles down to a stable form, via a stage called "ringdown", where any distortion in the shape is dissipated as more gravitational waves.To calculate exactly the GW emitted in the complete merge process it is necessary to solve the full equations of general relativity.This can be done only in numerical relativity simulations [16,17,19,20].Some simplified algebraic models can be used, for instance, in the "inspiral"stage.In Section 4 using the general relativity theory we estimate approximately the amplitude ("strain") of the GW signals detected by the LIGO in the "inspiral"stage.

Event GW150914 observation
The existence of stellar-mass BBHs (and gravitational waves themselves) was finally confirmed on September 14, 2015 at 09:50:45 UTC when the two LIGO detectors H1 and L1 have detected simultaneously a wave transient gravitational-wave signal.The signal sweeps upwards in frequency from 35 to 250 Hz with a peak gravitational-wave strain h = 1.0 X 10 −21 .This transient signal measured in H1 and L1 is seen in Figure 5 [12].The signals h(t) can be divided into three phases : "inspiral", "merger"(or "plunge") and "ringdown", that will be clearly shown in Figure 7.The merging BBH with a pair of BH with masses ∼30 M Θ was ∼ 1.3 billion light years away and in its final moments of spiraling inward and merging released around 3 M Θ as gravitational energy, peaking at a rate of 3.6 10 49 W -more than the combined power of all light energy radiated by all stars on the observable universe put together [19,20].
In Figure 6 we compare the strain h(t) measured by H1 and L1 .We verify that they are very similar, in good agreement within the experimental errors.
Observing the evolution of the gravitational-wave signal h(t) one could deduce [16,8] that the most plausible explanation for this is that a binary system formed by black holes (BH) with masses m 1 and m 2 ( Figs. 7) initially moving in a circular orbit begin to move in spiral around the center of mass and finally merging.With this model the predicted strain h(t) of the "gravitational luminosity"was calculated by the LIGO collaborators and is shown in Fig. 7 compared with the observed strain that has a peak |h| ∼ 10 −21 .Figure 7 clearly shows that the signals can be divided into three phases : "inspiral", "merger"and "ringdown".In the initial inspiral phase the BH move freely, in the merger phase they begin to "plunge"and, finally, the "ringdown"is the post-fusion phase.

Event GW151226 observation
On December 26, 2015 at 03:38:53 UTC the twin LIGO detectors H1 and L1 have detected simultaneously [20] a second wave transient gravitational-wave signal GW151226 shown in Figure 8.The signal persisted in the LIGO frequency band for approximately 1 s, from t = -1.0 up t = 0, increasing in frequency and amplitude over about 55 cycles from 35 to 450 Hz, and reached a peak gravitational strain Figura 6: Comparison between the gravitational-wave signals h(t) measured by H1 and L1. 16(Publicly available image at https://www.ligo.caltech.edu/)|h| ∼ 3.4 X 10 −22 .In this interval, that is, from t = 1.0 s up to t ≈ 0 s the BBH is practically in the "inspiral"stage.At t ≈ 0 s the BHs "merge"and the GW amplitude reaches its peak as can be seen in the third row of the Fig. 8 [Signal-to-noise (SNR) time series].For t > 0, we have the "ringdown"stage.
Comparing Figs.5 -6 and 8 we see that they are very similar.That is, observing the evolution of the gravitational-wave signal h(t) one could also deduce that the most plausible explanation is that a BBH with masses m 1 ≈ 14.2 M Θ and m 1 ≈ 7.5 M Θ (see Fig. 7) initially moving in a circular orbit begin to move in spiral around the center of mass and finally merging.With this model the predicted strain h(t) of the "gravitational luminosity"was calculated by the LIGO collaborators and is shown in Fig. 9.According to these estimations for three different times in the "inspiral"stage (seen in the first row of Fig. 9) we have h ∼ 0.17 10 −21 , h ∼ 0.18 10 −21 and finally, close to the merging point t ∼ 0 s, h ∼ 0.34 10 −21 .
From Figs.4-7 and Figs.8,9 we observe a good agreement within the experimental errors between theory and measurements.Comparing orbital with observed GW frequencies one can see [16,20] that the emitting objects of the binary systems can only be very compact BH.These results confirm that GW have finally been detected and that they are generated by merging BBH system.In next Section using the GTR we estimate approximately the GW strain h(t) in the "inspiral"stage of a BBH system.In order to explain the gravitational-wave strain h(t), we assumed that it was due to a binary system of two BH with masses m 1 and m 2 initially moving in a circular orbit that begin to move in spiral around the center of mass and finally merging [16].(Publicly available image at https://www.ligo.caltech.edu/)

Strain estimation for a binary-star system
We remark that our calculations are not rigorous; they have only a didactical objective.As we have said before, to describe exactly the detected signals of the GW in the complete merge BBH process of the BBH ("inspiral", "plunge"and "ringdown") it is necessary to use a completely relativistic model that can be treated only with numerical simulati- To explain the gravitational-wave strain h(t), we assumed that is was due a binary system (two black holes) with masses m 1 and m 2 initially moving in a circular orbit that begin to move in spiral around the center of mass and finally merging. 20(Publicly available image at https://www.ligo.caltech.edu/) ons [16. 17, 19, 20].Our approach describes only approximately the "inspiral"phase of the observed BBH.
Let us begin calculating the GW strain h(t) emitted by a generic binary-star system.We take R as the distance between the detector and the GW source that has a dimension D. According to preceding papers [4,8], the emitted GW in the radiation zone (R » D) is represented by the second-rank tensor where R is the radius and Q αβ (t) is the tensor quadrupole moment of the source given by For a binary system (Figure 10) composed by stars with masses m 1 and m 2 separated by a distance r one can show [9] that where µ = m 1 m 2 /(m 1 + m 2 ) and ω is the orbital angular frequency.Under these conditions, and using ( 14) and ( 16), one sees that the GW would be given by Figura 10: Binary-star system with masses m 1 and m 2 in the plane (x,y).
showing that the GW frequency is ω g = 2ω.Note that the GW is a tensor function, and not a scalar function as an electromagnetic wave.From (17) we verify that the strain h(t) could be written as h(t) ∼ h cos(2ωt) = h cos(ω g t) , where the GW angular frequency is ω g = 2ω.

Strain estimation
To obtain the effective strain h of the GW that arrive at the detector it is necessary to take into account the contributions of all h αβ tensor components (14) to the gravitational luminosity = L GW = dE/dt = total gravitational energy/s.The luminosity GW is the total energy lost by the source irradiated in all directions arriving on a spherical surface with radius R .Thus, according to reference 4,8 the energy flux Φ n = (d 2 E/dt 2 dA) (energy/m 2 s) that passes by the area dA = dA n = R 2 dA n = dΩ n is given by where ...
Using (18) we can calculate the Luminosity L GW : where the brackets < f(t) is the time average of f(t).
As the energy of the GW in the radiation zone is transported by a plane wave with amplitude h and orbital frequency ω the average flux value < Φ n > is given by: 4 Thus, with (19) and (20) let us define the following protocol to obtain h: first calculate L GW with (19) and second calculate h from L GW given by (20), from which we obtain h 2 = (8πG/ω 2 c 3 )(L GW /4πR 2 ). (21) For a binary-star we can show [6b] that From ( 21) and ( 22) we have h 2 = (2G/ω 2 c 3 ) (L GW /R 2 ),so that It is important to note that the relation ω 2 r 3 = 2GM can be used only when relativistic effects are very small, that is, when v 1 /c << 1 and v 2 /c << 1.

Strain h(ω) for the Event GW151226
In this event we have [20] R ≈ 440 Mpc and the BH masses, m 1 ≈ 14.2 M Θ and m 1 ≈ 7.5 M Θ .Thus, with these parameters (24) becomes h(ω) ∼ 0.7310 −23 ω 2/3 . (26) During signal detection the GW frequencies f g increased from 35 to 450 Hz.According to Fig. 9 in the "inspiral"stage we see that for 3 different times, h ∼ 0.17 10 −21 , h ∼ 0.18 x 10 −21 and close to the merging time t ∼ 0 s, h ∼ 0.34 x 10 −21 .As the GW frequencies f g are two times larger than the orbital frequencies, for initial orbital frequencies

Figura 1 :
Figura 1: Scheme of a typical large Michelson interferometer to be used in the detection of gravitational waves.(Fig.10.5 of Ref. 11) it is shown a simplified diagram of the Advanced LIGO detector (see details in reference 12), a modified Michelson interferometer: each arm is formed by two mirrors connected to test masses, separated by L x and L y = L = 4 km (see reference

Figura 2 :
Figura 2: Simplified diagram of the Advanced LIGO detector not in scale with the arms along the axis x and y with lengths L x = L y = L = 4 km .Insert (a): locations of the LIGOs at Hanford (H1) and Livingstone (L1) distant ∼3,000 km with arms in different orientations.(Publicly available image at https://www.ligo.caltech.edu/)

Figura 4 :
Figura 4: The strain sensitivity h min / √ Hz for LIGO as a function of the gravitational frequency f in Hz. (Publicly available image at https://www.ligo.caltech.edu/)

Figura 5 :
Figura 5: The gravitational-wave event detected simultaneously by the two LIGO interferometers H1 and L1 on September 14, 2015 at 09:50:45 UTC.Detailed description of this figure is given in reference 20.It is fundamental to note that the signals detected in H1 and L1 are, within the experimental errors, the same.Note that the frequencies (Hz) shown in Fig.5 are the frequencies ω g of the detected GW. (Publicly availble image at https://www.ligo.caltech.edu/)

Figura 7 :
Figura 7: In order to explain the gravitational-wave strain h(t), we assumed that it was due to a binary system of two BH with masses m 1 and m 2 initially moving in a circular orbit that begin to move in spiral around the center of mass and finally merging[16].(Publicly available image at https://www.ligo.caltech.edu/)

Figura 8 :
Figura 8: GW151226 observed by the LIGO Hanford (left column) and Livingston (right column) detectors, where times are relative to December 26, 2015 at 03:38:53.648UTC.It is fundamental to note that the signals detected in H1 and L1 are, within the experimental errors, the same [20].(Publicly available image at https://www.ligo.caltech.edu/)

Figura 9 :
Figura 9: Estimated gravitational-wave strain from GW151226 projected onto the LIGO Livingston detector with times relative to December 26, 2015 at 03:38:53.648UTC.To explain the gravitational-wave strain h(t), we assumed that is was due a binary system (two black holes) with masses m 1 and m 2 initially moving in a circular orbit that begin to move in spiral around the center of mass and finally merging.20 (Publicly available image at https://www.ligo.caltech.edu/)