It was produced processing of the LIGO 14.09.2015 registration data. It was established that the chirp signals (signals of merger of black holes) are absent in the data. This is proved by using a coherent filtering and also two-stage causal filtering. Soliton-like signals such as wavelet “Mexican hat” are found by filtering based on the Butterworth filter. These signals are different in polarity, and their spectra are quite different. Doubts about the authenticity of the detection of chirp signals, which announced on February 11, 2016, are justified by the results the conducted analysis.
In 1916, Albert Einstein predicted the existence of gravitational waves [
The discovery of the binary pulsar system PSR B1913 + 16 by Hulse and Taylor [
Experiments to detect gravitational waves began with Weber [
It was announced the observation of gravitational waves from merging black holes thanks to the two signals registered according to observatories LIGO Hanford (H1) and the LIGO Livingston (L1) in [
This paper presents the results of processing of registration data posted on the website LIGO ( https://losc.ligo.org/events/GW150914) [
The time of arrival of gravitational waves:
−Tdetect = 1,126,259,462.39 = September 14, 2015, 09:50:45.39 UTC, respective signals detected in the frequency band [35 ¸ 350 Hz], the events have a combined signal-to-noise ratio SNRannounced ≈ 24 [
These data were presented in February-May, 2016 on the website [
Because of this it is now the time of arrival of gravitational waves is specified without hundredths of a second on the website [
Tdetect = 1,126,259,462 = September 14 2015, 09:50:45 UTC.
It was announced in [
It is strange that removed hundredths of a second, no time of signal arrivals at the L1 and H1, but announced about the delay between the arrival of the GW- wave at different observatories to within a few thousandths of a second.
Time series (fragments) strains h L ( t ) , h H ( t ) are contained in the files:
L-L1_LOSC_4_V1-1126259446-32.txt, H-H1_LOSC_4_V1-1126259446-32.txt.
The corresponding random processes have a bandwidth [10 ¸ 2048 Hz].
Fragments h L ( t ) , h H ( t ) have a duration of 32 seconds, strain time series centered at
GPS 1126259462 = September 14 2015, 09:50:45 UTC.
Sampling frequency f d = 4096 Hz .
Start Time of the fragment Tstart = 1,126,259,446 = September 14 2015, 09:50:29 UTC, it is 16 seconds before the time of signal arrival Tdetect = 16.39 sec.
Further seconds and its share are relative T 0 = 09 : 50 : 45 or T start 0 = 09 : 50 : 29 .
On the page “Data release for event GW150914” ( https://losc.ligo.org/events/GW150914/) are given:
-results of processing carried out by the members of LIGO Collaboration (
-signals s H ( t ) , s L ( t ) (
Experimental templates, shown in
S ( t ) = A ( t ) cos ( φ 0 + 2 π ( f 0 t + β t 2 2 ) ) , (1)
where A ( t ) ―the signal amplitude;
f0―the carrier frequency;
φ0―the initial phase;
β―the speed of change of frequency;
f 0 = ( F max + F min ) / 2 ;
β = ( F max − F min ) / T s ;
Fmax, Fmin―minimum and maximum value of the signal frequency;
Ts―signal duration.
These signals are called linear frequency modulated (LFM) signals.
1) Linear-frequency modulated signal (LFM-signal), LFM-signal = Chirp signal.
Chirp signal stays the chirp signal after passing through the band filter.
Filtering is defined by the convolution:
Y ( t ) = X ( t ) ⊗ h f ( t ) , (2)
where h f ( t ) ―impulse response of a filter; X ( t ) ―input process; Y ( t ) ―the filtered process; ⊗ is the sign of the convolution.
2) The integral of the square of the signal remains practically unchanged if the signal in the filter band. This is a consequence of Parseval’s theorem.
F Y ( f ) = F X ( f ) ⋅ H ( f ) . (3)
| H ( f ) | ≈ 1 ⇒ ∫ t − T / 2 t + T / 2 X ( τ ) 2 d τ ≈ ∫ F max − Δ f F max + Δ f | F X ( f ) | 2 ⋅ 2 d f ≈ ∫ F max − Δ f F max + Δ f | F Y ( f ) | 2 ⋅ 2 d f ≈ ∫ t − T / 2 t + T / 2 Y ( τ ) 2 d τ , (4)
where H ( f ) ―the transfer function of the filter, H ( f ) = F ( h f ) ―a spectrum filter’s impulse response, | H ( f ) | ―amplitude-frequency characteristic(AFC) of the filter.
3) The spectrum of the signal after filtering remains practically unchanged if the signal is filtered in band, | H ( f ) | ≈ 1 .
4) If the registered signals S H ( t ) and S L ( t ) and they coincide with the delay and inversion (
s H ( t ) = − s L ( t + τ ) , x L ( t ) = s L ( t ) + n L ( t ) , x H ( t ) = s H ( t ) + n H ( t ) ,
s L f ( t ) = s L ( t ) ⊗ h f ( t ) , s H ( t ) = s H ( t ) ⊗ h f ( t ) ,
H ( f ) = F h f ( f ) , F s L f ( f ) = F s L ( f ) ⋅ H ( f ) , F s H f ( f ) = F s H ( f ) ⋅ H ( f ) ,
| H ( f ) | ≈ 1 ⇒ F s L f ( f ) ≈ − F s H f ( f ) ⇒ | F s L f ( f ) | ≈ | F s H f ( f ) | (5)
Amplitude spectra ( | F s H f ( f ) | ≈ | F s L f ( f ) | ) are almost equal after filtering in the filter band ( | H ( f ) | ≈ 1 ).
Detection of gravitational chirp signal is based on the selection of one from two alternative hypotheses:
H0: x H ( t ) = n H ( t ) , x L ( t ) = n L ( t ) , ―GW-wave chirp signal is absent;
H1: x H ( t ) = n H ( t ) + s H ( t ) , x L ( t ) = n L ( t ) + s L ( t ) , ―GW-wave chirp signal is present.
From the mathematical theory of statistics follows [
The matched filter has an impulse response equal to the inverted templates (
T L ( t ) = s L ( T template − t ) , T H ( t ) = s H ( T template − t ) , (6)
Ttemplate―template duration (duration of impulse response).
The matched filter of input records in the band [35 ¸ 350 Hz] with the signals, represented in
h 1 f ( t ) = h L f ( t ) ⊗ T L ( t ) , h 2 f ( t ) = h H f ( t ) ⊗ T H ( t ) , (7)
where h L f ( t ) = h L ( t ) ⊗ h f ( t ) , h L f ( t ) = h H ( t ) ⊗ h f ( t ) ―filtered processes after bandpass filtering in the band [35 ¸ 350 Hz].
The procedure for testing the hypothesis of the presence of chirp signals in the records:
1) Produced filtering h H ( t ) and h L ( t ) and find h H f ( t ) and h L f ( t ) .
2) The convolution is produced in accordance with (Equation (7)) and find h 1 f ( t ) and h 2 f ( t ) .
3) Modules h 1 f ( t ) and h 2 f ( t ) are compared with thresholds (P1, P2).
| h 1 f ( t ) | > P 1 ⇒ true H 1 else H 0 | h 2 f ( t ) | > P 2 ⇒ true H 1 else H 0
The signal/noise ratio is determined by the formula:
SNR = A max / σ n . (8)
The results of the matched filter in accordance with the Equations (7) are shown in
Such a small value of the SNR ( SNR real ≪ SNR announced = 24 ) leads to the inevitable conclusion about doubtfulness of the detected signals, which are stated in [
For example, a mixtures of signal and noise were formed to demonstrate the matched filtering efficiency:
x L ( t ) = h L ( t ) ⊗ h f ( t ) + s L ( t − τ s ) σ L σ s k SNR , x H ( t ) = h H ( t ) ⊗ h f ( t ) + s H ( t − τ s ) σ H σ S k SNR , (9)
where h f ( t ) ―impulse response of a band filter [35 ¸ 350 Hz];
h L ( t ) , h H ( t ) ―time series of strains;
σ s ―rms signal s ( t ) ;
σ L ( σ H ) ―rms noise n L ( t ) ( n H ( t ) ) ;
σ = Disp ( h ( t ) ) , kSNR―factor depending on SNR;
τ s ―time of signal arrival.
Mixtures of signal and noise x L ( t ) , x H ( t ) are calculated for the following parameters: time of model signal arrival τ s ≈ 15 sec , k SNR = 1.0 .
Matched filtering results x L ( t ) , x H ( t ) with model signals s H ( t ) and s L ( t ) shown in
The signal/noise ratio at the output SNRout is determined by the SNRin and gain at the expense of processing (BSNR):
SNR out = B SNR ⋅ SNR in . (10)
For formed processes x L ( t ) , x H ( t ) signal/noise ratio at the input in the simulation:
SNR inmod ≈ 2.323 .
It is clearly distinguished model signals to the 15-th second after matched filtering on
Input signal/noise ratio for gravitational signals in
SNR in L ≈ 5.5 , SNR in H ≈ 5.25
At such values SNRin by the Equation (10)
SNR out > 8.9 * 5.25 / 2.323 = 20 , (11)
that more than doubled when the model signals.
If the real signals are in
A small rise seen near 16-th second, corresponding SNRreal ≈ 3, as a result of matched filtering (the upper part of
Thus, the hypothesis H1 is not confirmed.
This means: the chirp signals are absent in the original records (
If the chirp signals are absent in registration data (
For this purpose, the optimum detector is used, which calculates the likelihood ratio:
L ( h H / s ) = ln P ( h H / s ) P ( h H / n ) , L ( h L / s ) = ln P ( h L / s ) P ( h L / n ) . (12)
From (Equation (12)) it follows that detection of the unknown signal is based on the computation of the functional in the first approximation [
U ( t ) = ∫ f l f h G t + T ( f ) G t − T ( f ) d f , (13)
where G t + T ( f ) , G t − T ( f ) ―energy process spectra before and after the current time t, f l , f h ―lower and upper frequency range of the signal.
Looking at
- Confident detection signals according to the observatory L1, H1;
- signal arrived first at L1 and ≤ 10 ms later at H1;
- SNR ≥ 7;
- arrival time at L1 earlier 16.4 sec.
This is done using a two-stage causal filtering:
1st step. Butterworth filter, 35 - 350 Hz band, the filter order = 8.
2nd step. Butterworth filter, 60 - 450 Hz band, the filter order = 4 (to notch components at frequencies near 32 Hz and 60 Hz).
According to [
However, these wave phases are not observed in the time interval 16.39 sec - 16.45 sec in
This apparent contradiction, and thereby violated №1 invariance property. Similarly, we can be sure of violating the invariance property №2.
In more detail the detected signals are represented in
Signals (
They have the property of solitons:
the higher the frequency, the greater the amplitude,
max A L > max A H , f L > f H
Spectra of these signals are shown in
It follows from
max S H ( f ) ≈ 165 Hz , max S L ( f ) ≈ 2 0 5Hz ;
Δ F = F max L − F max H ≈ 40 Hz = 205 − 165 Hz .
Spectrum of Signal L1 is a higher frequency than spectrum of signal H1 and invariants (№3, 4) are not performed.
This is a clear contradiction to the declared data about detected signals [
It means: s H ( t ) ≠ − s L ( t + τ ) !
Different types of filtration are used in detection of signal with a priori unknown form: bandpass, whitening(rectify), etc. [
Physically realizable filter whitening is based on the use of optimal filtering Wiener-Kolmogorov theory [
In [
Whitening (rectify) is proposed to conduct in the frequency domain in [
Y ( t ) = F − 1 ( F X ( f ) F h ( f ) ) , (14)
where F X ( f ) ―Fourier transform of X ( t ) ;
F h ( f ) ―transfer function of the whitening filter;
F − 1 ―inverse Fourier transform.
The transfer function of the whitening filter is in [
R whitening ( w ) = 1 G n ( w ) , Im whitening ( w ) = 0. (15)
This function has a zero phase.
The filter (Equation (14), Equation (15)) is physical be unrealizable, because any physically realizable filter performs a phase shift [
For example, the Butterworth filter is physically realizable (causal) filter.
If the real (true) signals shown in
However, the detection of the signal is observed on the 0.395 sec on the results of optimal detector (
Also according to [
The difference between the times of the signal arrivals indicates a contradiction with the assumption of trueness of image of waveform in
The difference, which is equal to 60 milliseconds, is the result of violation of the principle of causality in the processing of registration data [
The reason is that the transfer function (Equation (15)) defines a physically unrealizable (not causal) filter, its impulse response contains non-zero values for t < 0, and the use of such a filter leads to false detection of signals, shifts the time of their arrival and distortion of their forms (
Chirp signals are not visible after causal filtering (
1) Chirp signals are absent in data LIGO 14.09.2015.
This is established by means of a matched filter and two-stage filtration, therefore signals of the merger of black holes are absent in these data.
2) Signals of type “soliton” are distinguished by means of bandpass filtering (two steps).
Two pulse signals of different polarity were found in these records.
The corresponding wavelet is: “Mexican hat” or “Sombrero”.
3) Optimal detector detects these signals, SNR ≥ 7.
4) Signal arrived first at L1and ≤10 ms later at H1.
5) The signal spectra have a maximum of:
max S H ( f ) ≈ 16 5Hz , max S L ( f ) ≈ 205 Hz .
The spectra differ significantly as a result of a violation of invariance of spectra of signals at a filtration (properties №3, 4).
6) It is found that the detected signals at L1 and H1:
s H ( t ) ≠ − s L ( t + τ )
7) Physical unrealizability (not causal) filter is used to whitening the registration data in [
8) It is necessary to be based on the principle of causality when processing signals in physics experiments.
P.S. This article describes a repeatable physical experiment.
Initial data are known: parameters of used Butterworth filters, time intervals of registration and also file names.
Data of registration can be downloaded from website:
https://losc.ligo.org/events/GW150914.
This experiment can be repeated to check the main conclusions.
Lukanenkov, A.V. (2017) Experimental Confirmation of the Doubts about Authenticity of Event GW150914. Journal of Applied Mathematics and Physics, 5, 538-550. https://doi.org/10.4236/jamp.2017.52046