On the generation of magnetostatic solutions from gravitational two-soliton solutions of a stationary mass

In the paper, magnetostatic solutions of the Einstein–Maxwell field equations are generated from the gravitational two-soliton solutions of a stationary mass. Using the soliton technique of Belinskii and Zakharov (Sov Phys JETP 48:985, 1978, Sov Phys JETP 50:1, 1979), we construct diagonal two-soliton solutions of Einstein’s gravitational field equations for an axially symmetric stationary space-time and investigate some properties of the generated stationary gravitational metric. Magnetostatic solutions corresponding to the generated stationary gravitational solutions are then constructed using the transformation technique of Das and Chaudhuri (Pramana J Phys 40:277, 1993). The mass and the dipole moment of the source are evaluated. In our analysis we make use of a second transformation (Chaudhuri in Pramana J Phys 58:449, 2002), probably for the first time in the literature, to generate magnetostatic solutions from the stationary gravitational two-soliton solutions which give us simple and straightforward expressions for the mass and the magnetic dipole moment.


Introduction
In the paper, magnetostatic solutions of Einstein-Maxwell field equations are constructed from the stationary twosoliton solutions of Einstein's vacuum field equations in general relativity. For this purpose we first construct diagonal two-soliton solutions for a stationary axially symmetric space-time using the soliton technique of Belinskii and Zakharov [1,2]. The generated stationary gravitational solutions are then used to construct static magnetovac solutions of Einstein-Maxwell field equations. We have also investigated some properties of the generated stationary a e-mail: avinandac@gmail.com b e-mail: schaudhuri.phys@gmail.com gravitational two-soliton metric. The solutions are found to be asymptotically flat containing monopole, dipole and other higher mass multipoles. Under some restrictions on the constants, it is shown that the solutions reduce to the Schwarzschild [3], Kerr [4] and Kerr-NUT [5] metrics. When no such restrictions are imposed, the solutions may be interpreted as the nonlinear superposition of the Kerr metric with an external gravitational field. Our solution possesses an event horizon. The surface area of the horizon is evaluated.
Magnetostatic solutions are then constructed from the generated stationary gravitational two-soliton metric using the simple transformation of Das and Chaudhuri [6]. To relate the metric coefficients obtained by the soliton technique and those of the static magnetovac solutions, we have used a second transformation presented by Chaudhuri [7] by which magnetostatic solutions are obtained in a straightforward way from the gravitational two-soliton solutions.
The mass and the magnetic dipole moment of the source are evaluated. It may be mentioned here that the method followed in the paper to construct magnetostatic solutions deserves attention, although the transformation techniques are already known. The use of the second transformation makes our procedure of obtaining magnetostatic solutions from stationary gravitational two-soliton solutions very simple and probably was not used before in the literature for the purpose. The results are obtained in a compact form. It is expected that this is indeed an addition to the list of new solution generating techniques of Einstein/Einstein-Maxwell field equations.
In Sect. 2, the soliton technique of Belinskii and Zakharov [1,2] is briefly described. In Sect. 3, stationary axially symmetric vacuum solutions using the soliton technique are generated. Magnetostatic solutions corresponding to the generated stationary gravitational solutions are constructed in Sect. 4. Our conclusion follows in Sect. 5.

The soliton technique
The soliton technique discovered by Belinskii and Zakharav (BZ) [1,2] is a very useful method of generating exact solutions of Einstein's gravitational field equations. Several authors have presented a number of exact solutions of Einstein's field equations using this technique. Since the method of BZ is well known in the literature, here we describe the technique in brief. For a detailed study of the soliton technique one may consult references [1,2,[8][9][10][11][12].
A general stationary axially symmetric line element in (r, z) coordinates may be represented by where a, b take values 1, 2 and t, φ = x 1 , x 2 , respectively. The metric coefficients g ab and ν are functions of r, z only. The Einstein field equations for the metric (1) are [2] U r + V z = 0, where U = rg r g −1 and V = rg z g −1 The subscripts r, z denote partial differentiations. The metric coefficient ν can be determined from Eqs. (3) and (4) once g is known and the integrability condition of ν is Eq. (2). For physically acceptable solutions we have the supplementary condition [9] det g o = −r 2 , where g 0 is a (2 × 2) matrix associated with g ab . Soliton solutions to Eq. (2) are obtained by solving the following two differential equations [1,2]: where λ is a complex spectral parameter and the eigenfunction is a two-dimensional matrix which satisfies the relation Once is known, the generated metric coefficients are obtained from the following relations [2,9,12]: where m k oc and ω k are arbitrary constants and μ k 's are the pole trajectories. The indices k and l can take values from 1 to N . Here N is the number of solitons i.e. the number of poles that appears in the scattering matrix. For N = 1, it is called one-soliton solutions, N = 2, is the two-soliton solutions and so on.
The new metric g ab given by Eq. (10) is a solution to Eq. (2) but, in general, it does not satisfy Eq. (6). A physically acceptable g = g ph ab is defined by the relation [9] For general N -soliton solutions the function ν is given by [9,12] where C N are arbitrary constants and ν 0 is the ν-function of the static metric. The expression for ν 0 is given in Sect. 3.

The generated stationary gravitational solutions
In this section, two-soliton solutions of Einstein's vacuum field equations for an axially symmetric stationary metric are constructed using a Laplace's solution as seed. The solutions are found to be asymptotically flat and contain monopole, dipole and other higher mass multipoles. We consider a general static axially symmetric metric where ψ = ψ(r, z) is a Laplace solution and ν 0 , the νfunction of the static metric, is given by [9] ν 0 (ψ) = −ψ + where the subscripts r and z denote partial differentiations. We take Laplace's solution ψ as seed in the general static metric and assume g 0 as satisfying Eq. (6).
Since the seed solution is diagonal, we take the eigenfunction as a 2 × 2 diagonal matrix and write in the form such that F k = ψ when λ = 0.
The function F k = F k (r, z, λ) satisfies the following relations [9]: For two-soliton solutions there are two poles μ 1 and μ 2 .
Since the solutions are required along the pole trajectories λ = μ k , F k, should have two values, F 1 and F 2 .
In prolate spheroidal coordinates (x, y) defined by we take the Laplace solution ψ in the form where K , z 1 and α 0 are constants. Assuming that the arbitrary constant ω k appearing in Eq. (14) is real having values ω 1 and ω 2 for two-soliton solutions, the poles μ 1 and μ 2 in (x, y) coordinates are found to be With ψ given by Eq. (24), F 1 and F 2 , in prolate spheroidal coordinates, are obtained: Using Eqs. (10)- (17) and after cumbersome mathematical calculations, physically acceptable metric coefficients g ph ab are found to be The stationary metric (1) is thus completely given by Eqs. (28)-(37).
The asymptotic expansion of g ph 11 is obtained: The static limit of the solutions (28)-(30) is obtained with m 2 01 = 0, m 1 02 = 0. In the static limit one obtains This may be interpreted as the Schwarzschild metric embedded in an external gravitational field. When the parameter α 0 , representing the strength of the external field, vanishes, i.e. α 0 = 0, one recovers the Schwarzschild solution [3] g ph If we redefine the constants by and substitute α 0 = 0 in Eq. (28), we obtain Equation (43) is the Kerr-NUT solution [5] with NUT parameter h. Its reduction to the Kerr metric [4] is only a trivial task and can be obtained by substituting h = 0, where p and q are defined by The generated solutions presented in Eqs. (28)-(33) may thus be interpreted as the nonlinear superposition of a Kerr metric with an external gravitational field. The parameter α 0 is related to the strength of the external gravitational field.
The event horizon of the constructed stationary metric is at x = x h = 1. The infinite red shift surface (i.r.s.) depends on the values of y i.e. x i.r.s. = x(y) and in general, x i.r.s. > x h . At the poles (y = ±1), these two surfaces coincide. The event horizon is thus always surrounded by the infinite red shift surface similar to the Kerr metric.
The surface area of the event horizon is found to be In the absence of any external field, if one substitutes c 3 = 0 and K = m, one obtains the area of the event horizon corresponding to the Schwarzschild metric.

Static magnetic dipole solution
The static magnetic dipole solution corresponding to the generated two-soliton stationary metric (28)-(33) may be obtained using the technique of Das and Chaudhuri [6]. In the following we describe in brief the technique of Das and Chaudhuri to generate magnetostatic solutions from stationary gravitational fields.
The superscripts m and s refer to the magnetovac and stationary cases, respectively. Comparing the magnetostatic field equations (48)-(51) of the metric (47) with the stationary gravitational field equations (54)-(57) corresponding to the metric (53), one readily finds that static magnetovac solutions may be generated from the stationary axisymmetric vacuum solutions and vice versa via the transformations To obtain physically acceptable solutions, the imaginary quantity i in Eq. (58) is eliminated by the parameter change technique of Das and Banerji [15,16]. The generation of magnetostatic solutions from stationary gravitational fields is thus complete. From Eq. (58), it appears that once ω is known, A 3 can be evaluated. Since the metric (1) does not contain any term involving ω, it is very difficult to evaluate ω directly from the soliton solutions of Belinskii and Zakharov. However, one can find ω from the solutions of BZ as follows.
We consider a stationary metric where f, ω and γ s are the metric coefficients. Chaudhuri has shown that the metric coefficients f and ω for the line element (59) are related to the metric coefficients g ab of the metric (1) by [7] f = −g 11 , (60) ω = −(2g 11 ) −1 (g 12 + g 21 ). (61) For a diagonal metric, ( 0 ) 12 = ( 0 ) 21 , and we have g 12 = g 21 . One thus obtains and The metric coefficient g 22 is found to be From Eqs. (60), (63), (64) and under the coordinate transformation (23), it is easy to show that det g 0 = −r 2 . This satisfies the supplementary condition (6). From Eq. (60) and considering the line elements (53) and (59) we obtain the expression for e u : The other metric coefficient ω is obtained from Eq. (62) and is found to be where D is given by Eq. (33). In order to get real magnetostatic solutions we make the substitutions c 1→ic 1 and c 2 → ic 2 in Eqs. (65) and (66), and using the transformation (58) we obtain the Schwarzschild, Kerr and Kerr-NUT metrics. The constructed solutions are well behaved at spatial infinity containing monopole, dipole and other higher mass multipole moments. The area of the event horizon is also evaluated.
Magnetostatic solutions corresponding to the stationary gravitational solutions (28)-(36) are constructed in Eqs. (67)-(69). The asymptotic behaviour of the solutions is also discussed. In the absence of a magnetic field the solutions do not correspond to the Schwarzschild metric in a pure vacuum limit; rather one obtains a mass monopole term M = 2m. Bonnor also presented a magnetostatic solution [17] by his own method which also does not reduce to the Schwarzschild form in the absence of the magnetic parameter. The elegance of our technique is that a stationary gravitational two-soliton solution directly goes over to a static magnetovac solution and vice versa by two simple transformations, and the solutions are obtained in a compact form.
For a particular value of the magnetic parameter α (i.e. α = 1), the mass of the source vanishes and one obtains the dipole moment of a massless source. This corresponds to a finite current loop at large distances. The magnetic dipole moment directly depends on the strength of the current loop. Thus the study on magnetic dipole moment bears some relevance for the analysis of real astrophysical objects which are assumed to be usually charged due to a significant charge separation during their evolution.