Short-Period and Long-Period Effects of Weak Gravitational Waves

We compute the long-term orbital variation of a test particle orbiting a central body acted upon by normal incident of plane gravitational wave. We use the tools of celestial mechanics to give the first order solution of canonical equations of long-period and short-period terms of the perturbed Hamiltonian of gravitational waves. We consider normal incident of plane gravitational wave and characteristic size of bound—two body system (earth’s satellite or planet) is much smaller than the wavelength of the wave and the wave’s frequency nw is much smaller than the particle’s orbital np. We construct the Hamiltonian of the gravitational waves in terms of the canonical variables (l,g,h,L,G,H) and we solve the canonical equations numerically using Runge-Kutta fourth order method using language MATHEMATICA V10. Taking Jupiter as practical example we found that there are long period perturbations on ω,Ω and i and not changing with revolution and the short period perturbations on a, e and changing with revolution during the interval of time (tt0 ) which is changing from 04π.

1. Introduction

The tremendous penetrating capacity of gravitational waves is a reason why it is so difficult to detect them. Another reason is the very low intensities by which they are emitted. A great number of theoreticians and experimentalists work on gravitational radiation not only because of its important applications but also because of its fundamental role in science. In the domain of general relativity, the particular solution of the field equations leading to the wave solution is usually preferred in the domain of the post-Newtonian (or weak field) approximation. The reason is that the wave contributes to its own source, besides the effect of the relatively strong field of source itself; greatly complicate the detection of the wave close to its source. Hence it is better to study the simple properties of the waves by looking far into the wave zone where the fields are weak. The effects of incident plane gravitational waves on the orbital motion of gravitational bound two-body system were investigated by several authors with different approaches and various features of both orbits and the waves  -  . In this work we treat the effect of GW on the orbital motion by semi-analytical treatment and we determine the perturbation on the orbital elements using new approach differ than the previous works. We construct the Hamiltonian of GW and solve the canonical equations in terms of the Delaunay elements $\left(l,g,h,L,G,H\right)$ , using the numerical method of Runge-Kutta and algebraic manipulations programs MATHEMATICA V10. The resulting theory includes short-period as long-period for practical application on planet Jupiter.

2. Hamiltonian of GW

According to the weak field approximation of general relativity, there exist in empty space gravitational waves generated by different sources far in the universe, and propagated with velocity c. The field equations for empty space in the linear weak field approximation

${\square }^{2}{h}_{ik}=0$

Has the plane wave as the simplest solution. ${\square }^{2}$ is d’Alembertian operator is given by

${\square }^{2}=\frac{{\partial }^{2}}{\partial {x}_{i}^{2}}-\frac{1}{{c}^{2}}\frac{{\partial }^{2}}{\partial {t}^{2}}\text{\hspace{0.17em}}\left(i=1,2,3\right)$

With the Z axis chosen in the direction of propagation ${h}_{ik}$ , ${\square }^{2}$ becomes

${\square }^{2}=\frac{{\partial }^{2}}{\partial {z}^{2}}-\frac{1}{{c}^{2}}\frac{{\partial }^{2}}{\partial {t}^{2}}$

and then the partial derivatives of ${h}_{ik}$ with respect to x and y vanish. Thus a plane gravitational wave is determined by two quantities ${h}_{12}$ and ${h}_{11}=-{h}_{22}$ in other words, gravitational waves are transverse waves whose polarization is determined by a symmetric tensor of the second rank in the XY plane, the sum of whose diagonal terms ${h}_{11}+\text{\hspace{0.17em}}{h}_{22}$ is zero. So the plane wave can be written a the sum of two components

${h}_{11}={h}_{+}\text{cos}\left({n}_{w}t+{\alpha }_{1}\right)$ (1)

${h}_{12}={h}_{×}\text{cos}\left({n}_{w}t+{\alpha }_{2}\right)$ (2)

where ${n}_{w}$ is the frequency of the wave, ${\alpha }_{1}$ and ${\alpha }_{2}$ are the phase difference, ${h}_{+}$ and ${h}_{×}$ are the amplitude of the wave in the two orthogonal directions in the transverse plane. From the equation of geodesic deviation we obtain the acceleration components of GW as in  .

$\frac{{\text{d}}^{2}{n}^{i}}{\text{d}{t}^{2}}=\frac{1}{2}\frac{{\partial }^{2}{h}_{ik}}{\partial {t}^{2}}{n}^{k}$ (3)

From which we can obtained the potential like function U of GW in terms of the orbital elements

$U=\frac{1}{2}{h}_{1}\left({y}^{2}-{x}^{2}\right)-{h}_{2}xy$ (4)

where ${h}_{1}=\frac{1}{2}\frac{{\partial }^{2}{h}_{11}}{\partial {t}^{2}}$ and ${h}_{2}=\frac{1}{2}\frac{{\partial }^{2}{h}_{12}}{\partial {t}^{2}}$ , and we have

$x=r\left(\mathrm{cos}\Omega \mathrm{cos}\left(f+\omega \right)-\mathrm{cos}i\mathrm{sin}\Omega \mathrm{sin}\left(f+\omega \right)\right)$ (5)

$y=r\left(\mathrm{sin}\Omega \mathrm{cos}\left(f+\omega \right)+\mathrm{cos}i\mathrm{cos}\Omega \mathrm{sin}\left(f+\omega \right)\right)$ (6)

where $i,\Omega ,\omega$ and $f$ are the inclination, longitude of node, argument of perigee and the true anomaly of an orbit respectively. Substituting Equations ((1), (2), (5) and (6)) into Equation (4) we obtain

$\begin{array}{c}U=ϵ{r}^{2}\left\{\frac{{s}^{2}}{8}\left[\mathrm{cos}\left(2\Omega +{n}_{w}t+{\alpha }_{1}\right)-\mathrm{cos}\left(2\Omega -{n}_{w}t-{\alpha }_{1}\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-{A}_{g}\left[\mathrm{sin}\left(2\Omega +{n}_{w}t+{\alpha }_{2}\right)+\mathrm{sin}\left(2\Omega -{n}_{w}t-{\alpha }_{2}\right)\right]\right]\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\left(\frac{{c}^{2}}{16}+\frac{c}{8}-\frac{1}{16}\right)\left[\mathrm{cos}\left(2f+2\omega +2\Omega +{n}_{w}t+{\alpha }_{1}\right)+\mathrm{cos}\left(2f+2\omega +2\Omega -{n}_{w}t-{\alpha }_{1}\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-{A}_{g}\left[\mathrm{sin}\left(2f+2\omega +2\Omega +{n}_{w}t+{\alpha }_{2}\right)+\mathrm{sin}\left(2f+2\omega +2\Omega -{n}_{w}t-{\alpha }_{2}\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\mathrm{sin}\left(2f+2\omega -2\Omega +{n}_{w}t+{\alpha }_{2}\right)+\mathrm{sin}\left(2f+2\omega -2\Omega -{n}_{w}t-{\alpha }_{2}\right)\right]\right]\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\left(\frac{{c}^{2}}{16}-\frac{c}{8}+\frac{1}{16}\right)\left[\mathrm{cos}\left(2f+2\omega -2\Omega +{n}_{w}t+{\alpha }_{1}\right)+\mathrm{cos}\left(2f+2\omega -2\Omega -{n}_{w}t-{\alpha }_{1}\right)\right]\right\}\end{array}$ (7)

where

$c=\mathrm{cos}i,\text{\hspace{0.17em}}s=\mathrm{sin}i,\text{\hspace{0.17em}}ϵ=\frac{1}{2}{n}_{w}^{2}{h}_{×}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{and}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{A}_{g}=-\frac{{h}_{×}}{{h}_{+}}$ (8)

The Hamiltonian will be in the form

$H={H}_{0}+ϵ{H}_{1}$

${H}_{0}$ is the unperturbed Hamiltonian given by

${H}_{0}=-\frac{{\mu }^{2}}{2{L}^{2}}$ (9)

$\mu =GM$ the gravitational constant and $L=\sqrt{\mu a}$ . The perturbed Hamiltonian ${H}_{1}$ due to the gravitational wave is given by

$\begin{array}{c}{H}_{1}=\frac{{L}^{4}}{{\mu }^{2}}\frac{{r}^{2}}{{a}^{2}}\left\{\frac{{s}^{2}}{8}\left[\mathrm{cos}\left(2\Omega +{n}_{w}t+{\alpha }_{1}\right)-\mathrm{cos}\left(2\Omega -{n}_{w}t-{\alpha }_{1}\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-{A}_{g}\left[\mathrm{sin}\left(2\Omega +{n}_{w}t+{\alpha }_{2}\right)+\mathrm{sin}\left(2\Omega -{n}_{w}t-{\alpha }_{2}\right)\right]\right]\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\left(\frac{{c}^{2}}{16}+\frac{c}{8}-\frac{1}{16}\right)\left[\mathrm{cos}\left(2f+2\omega +2\Omega +{n}_{w}t+{\alpha }_{1}\right)+\mathrm{cos}\left(2f+2\omega +2\Omega -{n}_{w}t-{\alpha }_{1}\right)\end{array}$

$\begin{array}{c}\text{\hspace{0.17em}}\text{\hspace{0.17em}}-{A}_{g}\left[\mathrm{sin}\left(2f+2\omega +2\Omega +{n}_{w}t+{\alpha }_{2}\right)+\mathrm{sin}\left(2f+2\omega +2\Omega -{n}_{w}t-{\alpha }_{2}\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\mathrm{sin}\left(2f+2\omega -2\Omega +{n}_{w}t+{\alpha }_{2}\right)+\mathrm{sin}\left(2f+2\omega -2\Omega -{n}_{w}t-{\alpha }_{2}\right)\right]\right]\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\left(\frac{{c}^{2}}{16}-\frac{c}{8}+\frac{1}{16}\right)\left[\mathrm{cos}\left(2f+2\omega -2\Omega +{n}_{w}t+{\alpha }_{1}\right)+\mathrm{cos}\left(2f+2\omega -2\Omega -{n}_{w}t-{\alpha }_{1}\right)\right]\right\}\end{array}$ (10)

3. Equations of Motion

The Delaunay canonical variables $\left(l,g,h,L,G,H\right)$ are defined in terms of the Keplerian elements $\left(a,e,i,\Omega ,\omega ,M\right)$ by

$\begin{array}{l}l=M=n\left(t–\tau \right);\text{\hspace{0.17em}}L=\sqrt{\mu a}\\ g=\omega ;\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}G=\sqrt{\mu a}\left(1-{e}^{2}\right)\\ h=\Omega ;\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}H=G\mathrm{cos}i\end{array}$ (11)

The considered systems of differential equations in terms of canonical elements under the effect of GW are

$\frac{\text{d}{L}_{i}}{\text{d}t}=-\frac{\partial H}{\partial {l}_{i}};\text{\hspace{0.17em}}\frac{\text{d}{l}_{i}}{\text{d}t}=\frac{\partial H}{\partial {L}_{i}}\text{\hspace{0.17em}}\left(i=1,2,3,4\right)$ (12)

Then

${l}^{·}=\frac{\partial H}{\partial L};\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{g}^{·}=\frac{\partial {H}_{1}}{\partial G};\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{h}^{·}=\frac{\partial {H}_{1}}{\partial H}$ (13)

${L}^{·}=-\frac{\partial {H}_{1}}{\partial M};\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{G}^{·}=-\frac{\partial {H}_{1}}{\partial g};\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{H}^{·}=-\frac{\partial {H}_{1}}{\partial h}$ (14)

Using Equations ((11), thus (13) and (14)) yield

${a}^{·}=2\sqrt{\frac{a}{\mu }}{L}^{·}$ (15)

${e}^{·}=\frac{\left(1-{e}^{2}\right)}{n{a}^{2}e}{L}^{·}-\frac{\sqrt{\left(1-{e}^{2}\right)}}{n{a}^{2}e}{G}^{·}$ (16)

${i}^{·}=\frac{\mathrm{cot}i}{n{a}^{2}\sqrt{\left(1-{e}^{2}\right)}}{G}^{·}-\frac{\mathrm{csc}i}{n{a}^{2}\sqrt{\left(1-{e}^{2}\right)}}{H}^{·}$ (17)

${\Omega }^{·}=\frac{\partial {H}_{1}}{\partial H}=\frac{\partial {H}_{1}}{\partial i}\frac{\partial i}{\partial H}$ (18)

${\omega }^{·}=\frac{\partial {H}_{1}}{\partial G}=\frac{\partial {H}_{1}}{\partial e}\frac{\partial e}{\partial G}+\frac{\partial {H}_{1}}{\partial i}\frac{\partial i}{\partial G}$ (19)

${M}^{·}=\frac{\partial {H}_{0}}{\partial L}+\frac{\partial {H}_{1}}{\partial L}=\frac{\partial {H}_{0}}{\partial L}+\frac{\partial {H}_{1}}{\partial a}\frac{\partial a}{\partial L}+\frac{\partial {H}_{1}}{\partial e}\frac{\partial e}{\partial L}$ (20)

We solve this equations numerically using Runge-Kutta fourth order method, the mathematical program written by language of MATHEMATIC V10.

4. Solution and Results

We expand the trigonometric functions of the true anomaly $\frac{{r}^{2}}{{a}^{2}}$ , $\frac{{r}^{2}}{{a}^{2}}\mathrm{cos}2f$ and $\frac{{r}^{2}}{{a}^{2}}\mathrm{sin}2f$ of Hamiltonian (10) in series of the trigonometric functions of

mean anomaly $M$ up to order three in the eccentricity  . Changing the independent variable from the time $t$ to the mean anomaly $M$ , substituting the Hamiltonian into the equations of motion and solving the equations numerically to obtain the first order perturbation of GW on Jupiter. Table 1 represents the orbital elements of Jupiter in the epoch 2000, considering Bursts of GW with

frequency $56×{10}^{-6}$ Hz and amplitude $15×{10}^{-22}$ . ${\alpha }_{1}=0$ and ${\alpha }_{2}=\frac{\text{π}}{2}$ . The

variations of Jupiter’s elements due to the considering effect of gravitational waves are presented in Table 2, during the interval between zero and 4π as in Figure 1.

Table 1. The orbital elements of Jupiter Reference date 12.00 UT 1 Jan. 2000.

Table 2. The perturbed orbital elements of Jupiter due the incident of GW.

The variation of the perturbed semi-major X-axis represents the variation of time from 0 to 4π in degree and Y-axis represents δa in radian.
The Variation of the perturbed inclination X-axis represents the variation of time from 0 to 4π in degree and Y-axis represents δi in radian
The variation of the perturbed eccentricity X-axis represents the variation of time from 0 to 4π in degree and Y-axis represents δe in radian
The variation of the perturbed Ω X-axis represents the variation of time from 0 to 4π in degree and Y-axis represents δΩ in radian
The variation of the perturbed ω X-axis represents the variation of time from 0 to 4π in degree and Y-axis represents δω in radian

Figure 1. The short and long period perturbations of the orbital elements (semi-major, eccentricity, inclination, Ω, ω) of Jupiter during the interval of time $\left(t-{t}_{0}\right)$ , from 0 to 4π.

5. Conclusion

The resulting solution for our model from Equations (15) to (20) includes short-period as well as long-period terms for practical application it appears more expedient to calculate the perturbations on a revolution by revolution. The averaging process over the fast variable or short-period comes after the differentiated with respect to slow variable or long-period. The short-period perturbation on semi-major a, the eccentricity e, and the mean anomaly $M$ increase with each revolution during the time interval $\left(t-{t}_{0}\right)$ from zero to 4π. The long-period perturbation on the inclination i, Ω, and ω takes long time to change with revolution in the same interval as seen in Figure 1. This paper provides a good approach for long-period orbital evolution studies for planets orbiting under the effect of the gravitational waves.

Conflicts of Interest

The authors declare no conflicts of interest.

Cite this paper

Youssef, M. (2017) Short-Period and Long-Period Effects of Weak Gravitational Waves. International Journal of Astronomy and Astrophysics, 7, 230-237. doi: 10.4236/ijaa.2017.73019.

  Bertotti, B. (1973) Is the Solar System Gravitationally Closed? The Astrophysical Journal Letters, 14, 51.  Rudenko, V.N. (1975) Test Bodies under the Effect of Gravitational Radiation. Soviet Astronomy Letters, 19, 270.  Mashoon, B. (1978) On Tidal Resonance. Astrophysical Journal, 223, 285. https://doi.org/10.1086/156262  Mashoon, B. (1987) Wave Propagation in a Gravitational Field. Physics Letters A, 122, 299-304.  Turner, M.S. (1979) Influence of a Weak Gravitational Wave on a Bound System of Two Point-Masses. Astrophysical Journal, 233, 685. https://doi.org/10.1086/157429  Nelson, L.A. and Chau, W.Y. (1982) Orbital Perturbations of a Gravitationally Bound Two-Body System with the Passage of Gravitational Waves. Astrophysical Journal, 254, 735. https://doi.org/10.1086/159785  Ivanshchenko, A.V. (1987) The Variation of the Keplerian Elements of a Planetary Orbit under the Action of a Gravitational Wave. Soviet Astronomy Letters, 31, 76.  Thorne, K.S. and Braginsky, V.B. (1976) Gravitational Wave Burst from the Nuclei of Distant Galaxies and Quasars. The Astrophysical Journal Letters, 204, L1. https://doi.org/10.1086/182042  Lorenzo, I. (2014) Orbital Effect of a Monochromatic Plane GW with Ultra-Low Frequency Incident on a Gravitationally Bound Two-Body System. ScienceOpen Research, V2.  Straumann, N. (1984) General Relativity and Relativistic Astrophysics. Springer-Verlag, Berlin. https://doi.org/10.1007/978-3-642-84439-3  Smart, W.M. (1961) Celestial Mechanics. John Wiley, New York. 