_{1}

^{*}

It’s proposed to use a global seismic antenna (GSA) as a gravitational telescope, arbitrary “quiet” seismic stations are its elements, and aperture of GSA must be of the order 10,000 km. The relative displacements of various points of the Earth are detected by GSA, these displacements are described as quasi-harmonic elliptical signals generated by gravitational waves, their amplitude ≈ 2.5 * 0
^{?15} m. It is found that these waves cause deformation (strain) of the order h ≈ 10
^{?21}. Pulsars are a natural source of periodic waves. The fact of confident registration of gravitational wave is confirmed by detection of quasi-harmonic signals in the frequency band near 6.023 Hz for 90 hours (confidence probability of detection is close to 1). It is found that a small part of the rotation energy of associated pulsar (ε ≈ 10
^{?5}) is expended on the radiation corresponding to the gravitational wave.

Gravitational waves (GW-waves) are the inevitable consequence of many theories of gravity [

Gravitational telescopes (gravitational antennas) are created for the direct detection of gravitational waves. There are two types of gravitational wave detectors [

The low-frequency mechanical vibrations of a massive body caused by gravitational wave are measured in the first type. Detectors of this type use a massive metal bar, cooled to a low temperature. Weber bar is a famous example of such detector.

From the currently valid detectors on this principle operate spherical antenna MiniGRAIL, as well as antenna ALLEGRO, AURIGA, EXPLORER and NAUTILUS.

Detector of another type uses laser interferometry to measure gravitational wave induced motion between separated “free” masses.

This principle is applied in experiments LIGO, GEO600, TAMA-300 and VIRGO.

There are also projects for the detection of gravitational waves using seismographs [

It’s proposed to use a global seismic antenna (GSA) as a gravitational telescope, and arbitrary “quiet” seismic stations are its elements, its aperture must be of the order 10,000 km.

Demonstration of the GSA possibilities was carried by the results of registration of gravitational radiation at a frequency f ≈ 6 Hz.

General relativity (GR) predicts the existence of gravitational radiation (GW-radiation) as a perturbation of the gravitational field.

The corresponding gravitational wave (GW-wave) moves at the speed of light and is described by two independent components, arranged at an angle 45˚ to each other (

The simplest type: periodic compression and stretching of the body in two antiphase directions (

GW-wave leads to a relative deformation of the body, therefore, the absolute value of the deformation depends on the size of the deformable body.

For example, the Earth must be deformed into an ellipsoid in the field of gravitational radiation, stretched (compressed) perpendicular to the direction of the incoming wave, and the degree of stretching (compression) varies with the frequency of the incoming of gravitational radiation (

As a result, seismic waves are excited by the action of the tidal force F in the body of the Earth’s [^{−21} the value of earth surface displacements:

Such displacements are registered using seismometers and better use the network of seismometers located on different continents. Such network is a global seismic antenna (GSA).In this paper we consider the variant of GSA, based on 19 seismic stations of the International Monitoring System (IMS) of the Comprehensive Nuclear Test Ban Treaty (CTBT), placed on different continents (

To detect displacements of the Earth’s surface, threshold sensitivity of GSA must be:

Registration of displacements is carried against the background seismic noise and the main restriction on sensitivity of GSA is the natural seismic noise.

Natural seismic fields are generated by natural processes in the Earth’s interior, and have different physical nature. The field of seismic noise has two main components: diffuse and coherent [

The diffuse component is generated by a large number of simultaneously acting randomly distributed in space and unrelated sources. This is the result of the spontaneous seismic emission. Corresponding diffuse field accurately described as a homogeneous Gaussian field with a small radius of spatial correlation (not more than 10 - 20 km at frequencies greater than 1 Hz). Their spectral density has a “smooth” form without defined peaks at any frequencies [

Coherent component is generated by strong sources of noise, localized in space.

Fields of storm microseisms have clearly coherent nature, center frequency f ≈ 0.2 Hz. Local seismic sources (a distance of less than 100 km) cause the appearance of coherent components at frequencies greater than 1 Hz and they are specific to each station.

Placements of IMS stations are selected so as to minimize the coherent component, and diffuse component would be decisive.

It follows that the seismic noises at the stations practically are not correlated with each other, if the stations are spaced apart at distances of several thousand kilometers.

They do not have common sources of seismic noise.

Selected stations (

Power spectrum of noise at these stations are generally close to or lower than the corresponding curve for a “quiet” conditions, as shown in

The following characteristics of seismic noise at frequencies above 1 Hz will be considered further [

-normal distribution of seismic noise;

-noises on seismic stations are independent, if stations are spaced apart from one another at distances of several hundred kilometers or more;

-power spectral density

-power spectral density decreases with increasing frequency

-stationarity interval of seismicnoise ≈ 6 - 8 hours in the frequency band less than 0.1 Hz.

15 (fifteen) seismic arrays (SA) are the basis of the seismic network (_{sensor} times, N_{sensor}―number of SA sensors.

Frequencies more than 1 Hz allow to achieve the highest sensitivity of GSA and they are preferred for the detection of assumed periodic sources because many periods of pulsars are less 2 sec (T_{rot} < 2 s, rotation frequency f_{rot} > 0.5 Hz and f_{pul} > 1 Hz).

The response of test mass under the action of gravitational waves on the surface of the Earth [

where

In the particular case it is possible to convert the Equation (1), for example, if GW-wave falls on seismometer (detector) with an angular frequency

If the source is periodic, then the matrix A(t) and the force F(t) are periodic.

Therefore, solutions of the system (Equations (2)) are also periodic functions:

and they can be expanded in a Fourier series:

First approximations of these solutions are harmonic components (n = 1):

For an arbitrary direction of propagation of the GW-wave the useful signals are represented as:

orthonormal basis

_{E} = 6,371,000 m;

Signal (Equation (3))―the ellipse in the plane perpendicular to the direction p,

Canonical equations of ellipses, describing motion of the mass element in the coordinate system of the detector, are given in many works, in particular [

Ellipticity of solutions of Equations (2) is a property of stable solutions of linear differential equations with periodic coefficients.

The particular case

If GW-wave falls in the direction z, then motion near (x_{0}, y_{0}) is described by equations:

a)

b)

Point of circle

Motion of point

General Relativity is the geometrical theory, the principle equivalence of inertial and gravitational mass is observed. The bars of equal length are stretched identically and it does not depend on the material from which they are made.

Equation (1) describes the primary effect of GW-wave on matter.

Registered seismic process can be represented as:

where

Because the third component of the signal is zero, it is sufficient to consider the projection of seismic processes on the plane

where

Registered data are represented as the vector of observations:

where

N_{st}―number of stations seismic network; N―size of the sample data;

∆t―sampling interval time,

Detection of gravitational signal is based on the selection of one from two alternative hypotheses:

H_{0}:

H_{1}:

Optimal detection of signals (Equation (3)) is based on the evaluation of log-likelihood ratio:

where

Since seismic noises are independent for different stations, then:

where

Such representation greatly simplifies the evaluation of log-likelihood ratio, and allows using the processing methods of seismic array data.

The formation of total antenna beam (beam forming) is one of the basic methods of seismic data processing:

where τ_{i}―delays of arrival of the wave at different stations.

Result of the beam forming is elliptical signal, if elliptical signals are at stations of GSA.

The optimal functional of detection of a priori unknown signal is represented as follows:

where

After the simplification of this expression can be shown that the optimal functional of detection evaluates the energy of gravitational radiation coming from an arbitrary point on the celestial sphere

Expression (Equation (5)) is an estimation of energy of the elliptically polarized seismic process (EPSP), (

Estimation of energy is carried out on the time fragment_{0} and duration T.

Next, we consider estimation

To estimate the energy 32 directions (beams) p were chosen:

and you must find the p, for which the maximum value of

Sources of harmonic signals, naturally, are associated with pulsars. The radio emission signals (RES) of pulsar PSR 1919 + 21 are shown in

The visible frequency (

The center frequency f_{cent} = 0.77 Hz significantly differs from the visible frequency, shift (

Instability of RES forms is observed for many pulsars [

Center frequency can differ a few percent from the visible frequency, specified in the catalogs of pulsars.

Pulse repetition periods of the observed pulsars lie in the range of ≈1.6 ms to ≈4.3 s. The emission of pulsars is generally strongly polarized, the degree of polarization of the radio emission is close to 100%, almost circular polarization radio emission is observed for some pulsars [

Gravitational wave radiation (GWR) of the pulsar is determined by the quadrupole moment of the source. Given a priori ignorance of the quadrupole moment, we can assume that the complexity of gravitational radiation is similar to the complexity of pulsar radio emission. Consequently, the spectrum of GWR is the convolution of the spectra of the modulating function and aoriginal signal and spectrum of GWR expands and “floats” near the center frequency.

Most of the energy of the signal (

The detection needs to produce in the frequency bands (

Detection of useful harmonic signals is carried out in bands having a width

Evaluation of seismic noise dispersion in the frequency band

where

Using this equation for total antenna beam in each array, the evaluation of seismic noise dispersion in the frequency band

If the signal

where τ_{i}―corresponding of signal arrival delays, and evaluation of the amplitude is carried out according to the formula:

Dispersion of the noise is decreased by V times due to the processing of seismic data for a long time:

where N_{st}―number of stations of GSA;

N_{uncorr}―number of uncorrelated time values;

T_{an}―duration of analysis fragment;

Δf―band of frequency analysis.

For example, dispersion of evaluation of amplitude:

The amplitudes of the velocity and displacement are related:

Corresponding standard deviation of the noise amplitude (displacement):

This is an upper rough estimation of GSA sensitivity threshold near 6 Hz, and it characterizes the detection capabilities of GSA. This estimate is obtained theoretically using known characteristics and properties of the seismic noises.

More accurate estimate can be obtained if one considers that the signal is not only the harmonic, but also has an elliptical polarization.

Optimum processing functional (Equation (4)) tuned to detect a harmonic elliptically polarized signals for different arbitrary sources.

The energy of noise can be represented as the sum of the energies of the polarized and unpolarized components, elliptical polarization of noise represents a strict condition, and it leads to a significant reduction of a level of the elliptical polarized component:

Processing of GSA data was performed on a time interval 00^{00} 27.02.2009 - 18^{00} 2.03.2009.

By averaging the values of the energies at a frequency f by space, we estimate the energy, received by the antenna at a given frequency

Estimating the energies on time fragments (shift = 1 hour), you can get a diagram of spectral-time analysis, ^{−27} (m/s)^{2} to 35 ´ 10^{−27} (m/s)^{2}.

Analyzing _{0} = 6.023 Hz.

Estimating the significance level (the error of the 1^{st} kind), it is possible to confirm the hypothesis that the fact of energy ejection at the frequency f_{0} = 6.023 Hz is not random.

It is also possible to evaluate the average energies (

Spectrum

Using the values of

The amplitude of the ejection of

where

The error of the 1^{st} kind:

where ξ = N (0,1)―normally distributed random variable.

The value α is close to zero,the confidence probability of detection

The minimum values

Achieved threshold of GSA sensitivity at frequencies close to 6 Hz:

This sensitivity of GSA allows detecting deformation of the Earth (

The average value of energy of the detected signal at the frequency f = 6.023 Hz for 90 hours is equal to ^{−13} m/s and corresponding displacement:

Thus, it is experimentally confirmed (Equation (11)), that relative displacements (≈10^{−15} m) of different points of the Earth are detected using the optimal processing data of global seismic antenna.

The semimajor axis is more than twice the semiminor axis for the corresponding detected signals [

and therefore the eccentricity

These motions are an objective reality, because a false alarm is practically zero and confidence probability P_{conf} ≈ 1.

Therefore search of possible sources of the corresponding signals is performed further in this paper.

It should be noted that the detection of GW-signals (Equation (4)) is based on the phasing of the antenna at the source beyond Earth assuming that the propagation speed of the wave equal to the speed of light, and that signals were detected during a prolonged period of time (90 hours).

Also, the local maxima are seen at the frequencies f = 6.12 Hz and f = 6.26 Hz with a smaller SNR, than at the frequency f = 6.023 Hz (

Errors of 1^{st} kind

Local deformation caused by signal (Equation (3)) can be estimated:

Estimation of D_{loc}, h_{0} can be obtained (if f = f_{0}), using Equation (5):

and local deformation (strain) for f_{0} = 6.023 Hz can be estimated:

The relative change in the distance between stations is

The phase shift between the stations is unknown, in the worst case, ΔL_{i,j} is less than double amplitude of the detected signal.

Deformation space can be estimated by the change in the distance between any stations on the network, for example between stations ZALV and TXAR (

where

From a comparison of the two estimates of deformation D_{gl}, D_{loc}, it follows, that the velocity of gravitational waves and the speed of light have an equal order.

Equality of velocity of detected waves and the speed of light is not used in the evaluation D_{gl} and therefore it correct to use to estimate the velocity of propagation of the detected waves and you can use for this the analogue of the formula from the article [

The energy of a full rotation of the pulsar:

where J―moment of inertia about the axis of rotation;

f_{r}―rotation frequency of pulsar;

f_{pul}―frequency of GW-radiation.

The maximum of energy, emitted per unit of frequency interval, in the form of gravitational radiation is equal [

The spectral energy density of the gravitational-wave radiation (upper bound) at a distance r_{pul} [

Such estimates (Equation (15)) are obtained, assuming isotropy of radiation of pulsar (the energy is evenly distributed over the surface of a sphere of radius r_{pul}).

A more accurate estimate will be by considering directional pattern of gravitational radiation from the pulsar.

This radiation is determined by the quadrupole moment of a source that is not known a priori. The quadrupole moment exists, if the body is not symmetrical, and therefore, in general, the GW- radiation of pulsar is not isotropic and has a significant selectivity (concentration) of the radiation in space [

factor of pulsar radiation can be obtained a more accurate estimate

For a signal from the registered pulsar it must be satisfied [

This is a necessary condition for registration GW-radiation from real sources, this is a simple test for the reliability of detection.

Otherwise, if

The total energy that gravitational wave carry past a unit surface area of detector is [

Using Parseval’s theorem, the energy flow per unit area [

where

The integrand expression in (18) is the amount of energy transferred by GW-waves per unit frequency interval and per unit area (the analogue of the equation (37.31) in [

If the source is periodic (for example, if_{0}―frequency of gravitational radiation), in this case the spectral energy density of the gravitational-wave radiation is represented as:

Using Equations (5), (13), the spectral energy density of the gravitational-wave radiation can be represented as:

The value of energy of the detected signal equals _{0} = 6.023 Hz and according to the Equation (21):

Relation of the radiation source with the pulsar (pulsar association) carried out by the frequency of the radiation source and position recorded on the celestial sphere [

List of pulsars (rotation frequency f_{rot} ≈ 3 Hz, f_{pul} ≈ f_{0}) is shown in

Expression

Considering heterogeneity of the gravitational radiation, it will be spent an even smaller part of the rotational energy on this radiation, i.e.,

Similarly, you can perform the detection of other sources of GW-waves.

For example, candidates for the signals can be seen near 6.12 Hz and 6.26 Hz, SNR ≈ 6.5 and 7, respectively.

Pulsar | f_{rot} (Hz) | r_{pul}, kpc |
---|---|---|

B0450+55 | 2.93487997706 | 1.18 |

J0511−6508 | 3.10499386985 | 2.17 |

J0945−4833 | 3.015812519280 | 2.71 |

J1046+0304 | 3.06493262635 | 2.25 |

J1604−7203 | 2.92909024976 | 2.56 |

J1739−3951 | 2.92592314444 | 1.13 |

B1800−27 | 2.990292674144 | 3.62 |

J1820−0509 | 2.964537059313 | 2.81 |

J1852−2610 | 2.97320725542 | 2.26 |

B2020+28 | 2.912037613413 | 2.10 |

B2048−72 | 2.92966262614 | 1.05 |

B2123−67 | 3.0696382358 | 2.75 |

There are two types of gravitational signals: pulse and continuous quasi-harmonic.

LIGO detectors are aimed at the detection of signals of the first type [

Complementarity between LIGO and GSA is evident.

If the sensitivity of LIGO is increased to h ≈ 10^{−22} and more advanced methods of processing strain information are used then it is possible to detect gravitational waves from pulsars (frequency band of 35-350 Hz).

There are more than 75 pulsars of this frequency band that are located on distance less 6 kpc.

Detecting continuous signals (second type) by two methods will significantly increase the confidence of detecting gravitational waves.

By using the optimal methods of processing GSA data, 7 signals (see Appendix A) are detected with a confidence probability P_{conf} > 0.97. These signals lie on the planes which are perpendicular to the direction of the movement on certain points of the celestial sphere and have an elliptical polarization close to 90% (eccentricity ε > 0.866). The detected GSA signals are elliptically polarized, related to transverse waves and do not have longitudinal components.

Physical consistency of General Relativity or any other theory of gravity is checked for existence of a hypothetical scalar component of gravitational radiation [

These results support the basic provisions of General Relativity about gravitational waves:

-existence of two “canonical” polarizations h_{+} and h_{×};

-transverse nature of gravitational waves.

The required sensitivity of gravitational telescope has been achieved for the detection of gravitational waves (h ≈ 10^{−21}).

Relative displacements (≈10^{−15} m) of the different points of the Earth were detected; corresponding signals lie on the planes which are perpendicular to the direction of the radiation source. These signals have a high degree of elliptical polarization.

Detection of signals for a long time (90 hours) confirms the fact of registration of the GW-wave in the frequency band near 6.023 Hz.

Confidence probability of detection of GW-wave is close to 1.

Characteristics of detected signals (amplitude, center frequency) vary in time, i.e. GW-signal is quasi-har- monic signal with smoothly varying parameters. The source of these GW-signals is periodic (generally quasi- periodic), which is typical for pulsars.

The pulsar J0945-4833 (f_{rot} ≈ 3.01 Hz) is the most probable source of detected GW-signals, it’s expended about ε ≈ 10^{−5} of the rotational pulsar energy on gravitational radiation.

Described universal approach can be used for the detection of pulsars with f_{rot} > 0.5 Hz.

It is possible to increase the sensitivity of the gravitational telescope and to detect GW-waves with h < 10^{−21} by increasing the number of stations and accumulation intervals.

Alexander V. Lukanenkov, (2016) Gravitational Telescope. Journal of High Energy Physics, Gravitation and Cosmology,02,209-225. doi: 10.4236/jhepgc.2016.22019

The association of pulsars

Detected signals characterized by the following parameters:

The maximum value of

Pulsars, related with the detected signals, are found in a pulsar catalogue (for example, ATNF pulsar catalogue) by the rules:

1.

2. Right ascension (α_{pul}) and declination (δ_{pul}) of pulsar must be within the limits:

Then, pulsar is selected with rules:

Closeness in frequency

7 (seven) signals were detected from pulsars with confidence probability P_{conf} > 0.97 according to these rules, the frequencies of signals :1.457 Hz, 1.574 Hz, 1.662 Hz, 1.896 Hz, 2.033 Hz, 3.75 Hz, 6.023 Hz [

For example, signal has the characteristics φ_{0} = 110˚, λ_{0} = −60˚ if f_{0} = 6.023 Hz [