**International Journal of Astronomy and Astrophysics**

Vol.07 No.03(2017), Article ID:79254,9 pages

10.4236/ijaa.2017.73018

Path Integral Formulation for Wave Effect in Multilens System

Kazuhiro Yamamoto^{ }

Department of Physical Science, Hiroshima University, Higashi-Hiroshima, Japan

Copyright © 2017 by author and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

Received: July 30, 2017; Accepted: September 19, 2017; Published: September 22, 2017

ABSTRACT

A formula to investigate the wave effect in a multi-lens system is presented on the basis of a path integral formalism by generalizing the work by Nakamura and Deguchi (1999). The wave effect of a system with two lenses is investigated in an analytic way as a simple application to demonstrate usefulness of the formula and variety of wave effect in multi-lens system.

**Keywords:**

Cosmology, Gravitational Lensing, Wave Optics

1. Introduction

Recently several authors have investigated wave effect in gravitational lensing phenomenon, motivated by the possibility that the wave effect might be detected in future gravitational wave observatories [1] - [7] . The wave effect can be investigated by solving wave equation, in principle. However, it is difficult to solve the wave equation analytically for general configuration of lenses, excepting the simple case of a single Schwarzschild lens (e.g., [8] [9] , see also [10] and references therein, cf. [3] ). Indeed, investigation of the wave effect so far is restricted to special case of single lens model with the spherical symmetry, i.e., the Schwarzschild lens model and the singular isothermal sphere lens model (see [1] ).

Nakamura and Deguchi developed an elegant formalism for gravitational lens using the path integral approach [6] . The primary purpose of the present paper is to derive a useful formula to investigate the wave effect in general multi-lens system, by generalizing the formalism by Nakamura and Deguchi. We also apply it to a system with two Schwarzschild lenses to demonstrate usefulness of the formula. This paper is organized as follows: In Section 2, we present a generalized formula for multi-lens system. In Section 3, an application of the formula to a simple configuration with two lenses is considered. Section 4 is devoted to summary and conclusion. We use the convention $c=1$ .

2. Generalized Formulation

We start by reviewing the path integral formalism for gravitational lens phenomenon [6] . We consider the Newtonian spacetime with the metric

$\text{d}{s}^{2}=-\left(1+2U\left(r,\theta ,\phi \right)\right)\text{d}{t}^{2}+\left(1-2U\left(r,\theta ,\phi \right)\right)\left(\text{d}{r}^{2}+{r}^{2}\text{d}{\theta}^{2}+{r}^{2}si{n}^{2}\theta \text{d}{\phi}^{2}\right)\mathrm{,}$ (1)

where $U\left(r,\theta ,\phi \right)$ is the Newtonian potential. Propagation of massless field $\varphi $ is described by the wave equation

$\frac{\partial}{\partial {x}^{\mu}}\left(\sqrt{-g}{g}^{\mu \nu}\frac{\partial}{\partial {x}^{\nu}}\right)\varphi \left(r,t\right)=0,$ (2)

where g is the determinant of the metric and r represents the spatial coordinates $\left(r,\theta ,\phi \right)$ . We consider a monochromatic wave from a point source with the wave number k. We set

$\varphi \left(r,t\right)=\frac{A}{r}F\left(r\right){\text{e}}^{-ik\left(t-r\right)},$ (3)

where A is a constant, then Equation (2) yields

$\frac{{\partial}^{2}F}{\partial {r}^{2}}+2ik\frac{\partial F}{\partial r}+\frac{1}{{r}^{2}}\frac{1}{sin\theta}\frac{\partial}{\partial \theta}\left(sin\theta \frac{\partial F}{\partial \theta}\right)+\frac{1}{{r}^{2}}\frac{1}{si{n}^{2}\theta}\frac{{\partial}^{2}F}{\partial {\phi}^{2}}-4{k}^{2}U\left(r,\theta ,\phi \right)F=0.$ (4)

Assuming that the first term is negligible compared with the second term and $\theta \ll 1$ , Equation (4) reduces to

$i\frac{\partial F}{\partial r}=-\frac{1}{2}k{r}^{2}\left[\frac{1}{\theta}\frac{\partial}{\partial \theta}\left(\theta \frac{\partial F}{\partial \theta}\right)+\frac{1}{{\theta}^{2}}\frac{{\partial}^{2}F}{\partial {\phi}^{2}}\right]+2kU\left(r,\theta ,\phi \right)F.$ (5)

Due to analogy of Equation (5) with the Schrödinger equation, using the path integral formulation of quantum mechanics, the solution can be written as follows,

$F\left(r,\theta ,\phi \right)={\displaystyle \int}\mathcal{D}\Theta \mathrm{exp}\left[ik{\displaystyle {\int}_{0}^{r}}\text{d}r\left(\frac{1}{2}{r}^{2}{\stackrel{\dot{}}{\Theta}}^{2}\left(r\right)-2U\left(r,\Theta \left(r\right)\right)\right)\right],$ (6)

where the dot denotes the differentiation with respect to $r$ and $\Theta $ is used to represent the variables $\left(\theta \mathrm{,}\phi \right)$ , which are related by $\Theta =\left(\theta cos\phi \mathrm{,}\theta sin\phi \right)$ . The expression (6) means the sum of all possible path $\Theta \left(r\right)$ , fixing the initial point (source’s position) and the final point (observer’s position).

Let us consider multi-lens system, including n lenses as shown in Figure 1, in which the source is located at the origin of coordinate and an observer is located at $\left(r\mathrm{,}\Theta \right)=\left({r}_{N}\mathrm{,}{\Theta}_{N}\right)$ . The radius between the source and the observer is discretized by N segments, and each discrete radius from the source is labeled by ${r}_{j}$ with j from 1 to N. Separation neighboring two discrete radii is $\epsilon $ . We assume that the n lenses are located at the radius $r={r}_{{l}_{m}}$ with m from 1 to n, and that the thin lens approximation is valid for the lenses. In this case the

Figure 1. Configuration of multi-lens system and coordinates for path integral formula.

Figure 2. Lensing system considered in section 3.

explicit expression of Equation (6) is written as Equation (24) in Appendix (see also [6] ). After some computation, the path integral expression reduces to (see Appendix for details)

#Math_24# (7)

with

$\psi \left({\Theta}_{{l}_{m}}\right)=2{\displaystyle {\int}_{{r}_{{l}_{m}}-\delta r}^{{r}_{{l}_{m}}+\delta r}}\text{d}rU\left(r,\Theta \right),$ (8)

where ${r}_{{l}_{m-1},{l}_{m}}={r}_{{l}_{m}}-{r}_{{l}_{m-1}}$ , ${r}_{{l}_{n+1}}={r}_{N}$ , and ${\Theta}_{{l}_{m}}$ is the variable on the mth lens plane. As expected, stationary condition of the phase of the integrand in Equation (7) reproduces the lens equation in the multi-lens system [10]

${\nabla}_{{\Theta}_{{l}_{{m}^{\prime}}}}{\displaystyle \underset{m=1}{\overset{n}{\sum}}}\left(\frac{{r}_{{l}_{m}}{r}_{{l}_{m+1}}}{2{r}_{{l}_{m},{l}_{m+1}}}{\left|{\Theta}_{{l}_{m}}-{\Theta}_{{l}_{m+1}}\right|}^{2}-\psi \left({\Theta}_{{l}_{m}}\right)\right)=0,$ (9)

for each ${m}^{\prime}=1,\cdots ,n$ , which hints at a way to the geometrical optics limit. The Gaussian approximation around a stationary solution yields the result in the geometrical optics limit (see also [6] ).

3. Application to a Simple Configuration

In this section we consider a simple case with two lenses. Figure 2 shows the configuration: The Schwarzschild lenses with mass ${M}_{1}$ and ${M}_{2}$ are located at the radius ${r}_{1}$ and ${r}_{2}$ , respectively. The model considered here is not general, because the source and the two lenses are arranged to be on a straight line. However, this simplification allows us to perform integration of Equation (7) analytically. We start by rewriting Equation (7)

$\begin{array}{c}F=\frac{k}{2\text{\pi}i}\frac{{r}_{1}{r}_{2}}{{r}_{1,2}}{\displaystyle \int}\text{\hspace{0.05em}}{\text{d}}^{2}{\Theta}_{1}\mathrm{exp}\left[ik\left(\frac{{r}_{1}r}{2{r}_{1,2}}{\left|{\Theta}_{1}-{\Theta}_{2}\right|}^{2}-\psi \left({\Theta}_{1}\right)\right)\right]\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\times \frac{k}{2\text{\pi}i}\frac{{r}_{2}{r}_{3}}{{r}_{2,3}}{\displaystyle \int}\text{\hspace{0.05em}}{\text{d}}^{2}{\Theta}_{2}\mathrm{exp}\left[ik\left(\frac{{r}_{2}{r}_{3}}{2{r}_{2,3}}{\left|{\Theta}_{2}-{\Theta}_{3}\right|}^{2}-\psi \left({\Theta}_{2}\right)\right)\right]\end{array}$ (10)

with

$\psi \left({\Theta}_{j}\right)=4G{M}_{j}ln\left(\left|{\Theta}_{j}\right|\right)\mathrm{,}$ (11)

for $j=1,2$ , where ${r}_{i,j}={r}_{j}-{r}_{i}$ and the position of an observer is $\left({r}_{3}\mathrm{,}{\Theta}_{3}\right)$ . From Eqution (10) we have

$\begin{array}{c}F=\frac{k}{i}\frac{{r}_{1}{r}_{2}}{{r}_{1,2}}{\displaystyle {\int}_{0}^{\infty}}\text{\hspace{0.05em}}\text{d}{\theta}_{1}{\theta}_{1}^{1-4ikG{M}_{1}}\mathrm{exp}\left[ik\frac{{r}_{1}{r}_{2}}{2{r}_{1,2}}\left({\theta}_{1}^{2}+{\theta}_{2}^{2}\right)\right]{J}_{0}\left(\frac{k{r}_{1}{r}_{2}}{{r}_{1,2}}{\theta}_{1}{\theta}_{2}\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\times \frac{k}{i}\frac{{r}_{2}{r}_{3}}{{r}_{2,3}}{\displaystyle {\int}_{0}^{\infty}}\text{\hspace{0.05em}}\text{d}{\theta}_{2}{\theta}_{2}^{1-4ikG{M}_{2}}\mathrm{exp}\left[ik\frac{{r}_{2}{r}_{3}}{2{r}_{2,3}}\left({\theta}_{2}^{2}+{\theta}_{3}^{2}\right)\right]{J}_{0}\left(\frac{k{r}_{2}{r}_{3}}{{r}_{2,3}}{\theta}_{2}{\theta}_{3}\right),\end{array}$ (12)

where $\left|{\Theta}_{j}\right|={\theta}_{j}$ and ${J}_{0}\left(y\right)$ is the Bessel function of the first kind. Integration with respect to ${\theta}_{1}$ can be performed (see [10] )

#Math_44# (13)

where $\alpha $ is a real number which represents a constant phase and ${}_{1}{F}_{1}\left(a\mathrm{,}b\mathrm{;}y\right)$ is the Kummer’s function. With the use of the definition of the Bessel function

${J}_{0}\left(z\right)={\displaystyle \underset{L=0}{\overset{\infty}{\sum}}}\frac{{\left(-1\right)}^{L}}{{\left(L!\right)}^{2}}{\left(\frac{z}{2}\right)}^{2L},$ (14)

we have (see [11] )

$\begin{array}{c}F={\text{e}}^{i{\alpha}^{\prime}}{\text{e}}^{\text{\pi}kG\left({M}_{1}+{M}_{2}\right)}\Gamma \left(1-2ikG{M}_{1}\right)z{\displaystyle \underset{L=0}{\overset{\infty}{\sum}}}\frac{{\left(-i\right)}^{L}}{{\left(L!\right)}^{2}}\Gamma \left(1+L-2ikG{M}_{2}\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\times {\left(xz\right)}^{L}{}_{2}F{}_{1}\left(1-2ikG{M}_{1},1+L-2ikG{M}_{2},1;1-z\right),\end{array}$ (15)

where we defined

$z=\frac{{r}_{3}\left({r}_{2}-{r}_{1}\right)}{{r}_{2}\left({r}_{3}-{r}_{1}\right)},$ (16)

$x=\frac{k{r}_{2}{r}_{3}{\theta}_{3}^{2}}{2\left({r}_{3}-{r}_{2}\right)},$ (17)

${\alpha}^{\prime}$ is a real constant and ${}_{2}{F}_{1}\left(a\mathrm{,}b\mathrm{,}c\mathrm{;}y\right)$ is the Hypergeometric function. In the limit ${\theta}_{3}=0\left(x=0\right)$ , Eqution (15) reduces to

$\begin{array}{c}F={\text{e}}^{i{\alpha}^{\prime}}{\text{e}}^{\text{\pi}kG\left({M}_{1}+{M}_{2}\right)}\Gamma \left(1-2ikG{M}_{1}\right)\Gamma \left(1-2ikG{M}_{2}\right)z\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\times {}_{2}F{}_{1}\left(1-2ikG{M}_{1},1-2ikG{M}_{2},1;1-z\right),\end{array}$ (18)

and we have

${\left|F\right|}^{2}=\frac{4\text{\pi}kG{M}_{1}}{1-{\text{e}}^{-4\text{\pi}kG{M}_{1}}}\frac{4\text{\pi}kG{M}_{2}}{1-{\text{e}}^{-4\text{\pi}kG{M}_{2}}}{z}^{2}{\left|{}_{2}{F}_{1}\left(1-2ikG{M}_{1},1-2ikG{M}_{2},1;1-z\right)\right|}^{2}.$ (19)

We consider the coincidence limit that the distance between the two lenses becomes zero, i.e., ${r}_{1}={r}_{2}$ . In this case, from previous investigation (e.g., [10] ), $F$ should be

$\mathcal{F}\left(x\right)\equiv {\text{e}}^{\text{\pi}kG\left({M}_{1}+{M}_{2}\right)}\Gamma \left(1-2ikG{M}_{1}\right){}_{1}F{}_{1}\left(1-2ikG\left({M}_{1}+{M}_{2}\right),1;-ix\right),$ (20)

excepting a constant phase factor. We compare our result with Eqution (20). First, let us consider the limit ${r}_{1}={r}_{2}$ , i.e., $z=0$ of Eqution (15). Using the mathematical formula

#Math_61# (21)

in the case ${\theta}_{3}=0\left(x=0\right)$ , we can easily have

$\underset{z\to 0}{lim}{\left|F\right|}^{2}=\frac{4\text{\pi}kG\left({M}_{1}+{M}_{2}\right)}{1-{\text{e}}^{-4\text{\pi}kG\left({M}_{1}+{M}_{2}\right)}}={\left|\mathcal{F}\left(0\right)\right|}^{2}\mathrm{.}$ (22)

This is the expected result. Note that ${\left|F\right|}^{2}$ is regarded as the magnification factor. Figure 3 plots $R=\left|F/\mathcal{F}\left(0\right)\right|$ , as a function of z, where we fixed $kG{M}_{1}=kG{M}_{2}=1$ . Thus the maximum magnification depends on lens configuration significantly. Figure 4 plots $R=\left|F/\mathcal{F}\left(x\right)\right|$ as a function of x when fixing $z=0.5$ and $kG{M}_{1}=kG{M}_{2}=1$ , which indicates that interference pattern also depends on lens configuration. In these figures we have fixed

Figure 3. $R=\left|F/\mathcal{F}\left(0\right)\right|$ as a function of z with fixing $kG{M}_{1}=kG{M}_{2}=1$ .

Figure 4. $R=\left|F/\mathcal{F}\left(x\right)\right|$ as a function of x with fixing $z=0.5$ and $kG{M}_{1}=kG{M}_{2}=1$ .

$kG{M}_{1}=kG{M}_{2}=1$ , but result depends significantly on the parameters $kG{M}_{1}$ , $kG{M}_{2}$ and $z$ , as shown in Figure 4, which suggests variety of wave effect depending on lens configuration.

4. Summary and Conclusion

In the present paper, we have presented a general formula to investigate wave effect in multi-lens system, which has been derived with the use of path integral approach. The formula is expressed in terms of integration with respect to variables of lens planes. It is difficult to perform the integration in general cases, but a system with two Schwarzschild lenses is an example for which the integration can be performed in an analytic way. The model considered in the present paper is simplified and limited, however, it suggests variety of wave effect in multi-lens phenomenon. It is required to develop a numerical method to perform integration in general lens configuration in future work.

Acknowledgements

This work is supported in part by Grant-in-Aid for Scientific research of Japanese Ministry of Education, Culture, Sports, Science and Technology, No. 15740155.

Cite this paper

Yamamoto, K. (2017) Path Integral Formulation for Wave Effect in Multilens System. International Journal of Astronomy and Astrophysics, 7, 221-229. https://doi.org/10.4236/ijaa.2017.73018

References

- 1. Takahashi, R. and Nakamura, T. (2003) Wave Effects in the Gravitational Lensing of Gravitational Waves from Chirping Binaries. The Astrophysical Journal, 595, 1039-1051. https://doi.org/10.1086/377430
- 2. Seto, N. (2004) Strong Gravitational Lensing and Localization of Merging Massive Black Hole Binaries with LISA. Physical Review D, 69, Article ID: 022002 https://doi.org/10.1103/PhysRevD.69.022002
- 3. Yamamoto, K. and Tsunoda, K. (2003) Wave Effect in Gravitational Lensing by a Cosmic String. Physical Review D, 68, Article ID: 041302. https://doi.org/10.1103/PhysRevD.68.041302
- 4. De Paolis, F., Ingrosso, G., Nucita, A.A. and Qadir, A. (2002) A Note on Gravitational Wave Lensing. Astronomy and Astrophysics, 394, 749-752.https://doi.org/10.1051/0004-6361:20021258
- 5. Baraldo, C., Hosoya, A. and Nakamura, T.T. (1999) Gravitationally Induced Interference of Gravitational Waves by a Rotating Massive Object. Physical Review D, 59, Article ID: 083001. https://doi.org/10.1103/PhysRevD.59.083001
- 6. Nakamura, T.T. and Deguchi, S. (1999) Wave Optics in Gravitational Lensing. Progress of Theoretical Physics Supplement, 133, 137-153. https://doi.org/10.1143/PTPS.133.137
- 7. Nakamura, T.T. (1998) Gravitational Lensing of Gravitational Waves from Inspiraling Binaries by a Point Mass Lens. Physical Review Letters, 80, 1138-1141. https://doi.org/10.1103/PhysRevLett.80.1138
- 8. Deguchi, S. and Watson, W.D. (1986) Wave Effects in Gravitational Lensing of Electromagnetic Radiation. Physical Review D (Particles and Fields), 34, 1708-1718. https://doi.org/10.1103/PhysRevD.34.1708
- 9. Deguchi, S. and Watson, W.D. (1986) Diffraction in Gravitational Lensing for Compact Objects of Low Mass. Astrophysical Journal, Part 1, 307, 30-37. https://doi.org/10.1086/164389
- 10. Schneider, P., Ehlers, J. and Falco, E.E. (1992) Gravitational Lenses. Springer-Verlag, Berlin.
- 11. Magnus, W., Oberhettinger, F. and Soni, R.P. (1966) Formulas and Theorems for the Special Functions of Mathematical Physics. Springer-Verlag, Berlin. https://doi.org/10.1007/978-3-662-11761-3

A. Derivation of Equation (7)

In this Appendix, we review derivation of Equation (7) from Equation (6). We consider the configuration depicted as Figure 1. The source is located at the origin of coordinate and the position of an observer is specified by $\left(r,\Theta \right)=\left({r}_{N},{\Theta}_{N}\right)$ . The space between the source and the observer is discritized by $N-1$ planes. The radius of each plane is labeled by ${r}_{j}$ for $j=1,\cdots ,N-1$ . Note that ${r}_{N}$ specifies the plane of the observer. ${\Theta}_{j}$ is (angle) variable on the jth plane. We consider the system with n lenses, and assume that mth lens is located at the radius ${r}_{{l}_{m}}$ . Assuming the validity of the thin lens approximation, we introduce the two dimensional potential

$\psi \left({\Theta}_{{l}_{m}}\right)=2{\displaystyle {\int}_{{r}_{l}-\delta r}^{{r}_{l}+\delta r}}\text{d}rU\left(r,\Theta \right),$ (23)

for $m=1,\cdots ,n$ , respectively. In this case, the path integral formula (6) can be written as

$F=\left[{\displaystyle \underset{j=1}{\overset{N-1}{\prod}}}{\displaystyle \int}\frac{{\text{d}}^{2}{\Theta}_{j}}{{A}_{j}}\right]\mathrm{exp}\left[ik\left(\u03f5{\displaystyle \underset{j=1}{\overset{N-1}{\sum}}}\frac{{r}_{j}{r}_{j+1}}{2}{\left|\frac{{\Theta}_{j+1}-{\Theta}_{j}}{\u03f5}\right|}^{2}-{\displaystyle \underset{m=1}{\overset{n}{\sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\psi \left({\Theta}_{{l}_{m}}\right)\right)\right],$ (24)

where the normalization is chosen

${A}_{j}=\frac{2\text{\pi}i\u03f5}{k{r}_{j}{r}_{j+1}},$ (25)

so that $F=1$ in the limit $\psi =0$ . Then, Equation (24) is rephrased as

$\begin{array}{c}F=\left[{\displaystyle \underset{j=1}{\overset{{l}_{1}-1}{\prod}}}{\displaystyle \int}\frac{{\text{d}}^{2}{\Theta}_{j}}{{A}_{j}}\right]\mathrm{exp}\left[ik\left(\u03f5{\displaystyle \underset{j=1}{\overset{{l}_{1}-1}{\sum}}}\frac{{r}_{j}{r}_{j+1}}{2}{\left|\frac{{\Theta}_{j+1}-{\Theta}_{j}}{\u03f5}\right|}^{2}\right)\right]\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\times \left[{\displaystyle \underset{j={l}_{1}}{\overset{{l}_{2}-1}{\prod}}}{\displaystyle \int}\frac{{\text{d}}^{2}{\Theta}_{j}}{{A}_{j}}\right]\mathrm{exp}\left[ik\left(\u03f5{\displaystyle \underset{j={l}_{1}}{\overset{{l}_{2}-1}{\sum}}}\frac{{r}_{j}{r}_{j+1}}{2}{\left|\frac{{\Theta}_{j+1}-{\Theta}_{j}}{\u03f5}\right|}^{2}-\psi \left({\Theta}_{{l}_{1}}\right)\right)\right]\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\times \cdots \\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\times \left[{\displaystyle \underset{j={l}_{n}}{\overset{N-1}{\prod}}}{\displaystyle \int}\frac{{\text{d}}^{2}{\Theta}_{j}}{{A}_{j}}\right]\mathrm{exp}\left[ik\left(\u03f5{\displaystyle \underset{j={l}_{n}}{\overset{N-1}{\sum}}}\frac{{r}_{j}{r}_{j+1}}{2}{\left|\frac{{\Theta}_{j+1}-{\Theta}_{j}}{\u03f5}\right|}^{2}-\psi \left({\Theta}_{{l}_{n}}\right)\right)\right],\end{array}$ (26)

where $\u03f5$ is the separation between two neighboring planes. With the use of the following equality, which can be proven by the mathematical induction,

$\underset{j={l}_{m}}{\overset{{l}_{m+1}-1}{\sum}}}{r}_{j}{r}_{j+1}{\left|{\Theta}_{j+1}-{\Theta}_{j}\right|}^{2}=\u03f5\frac{{r}_{{l}_{m}}{r}_{{l}_{m+1}}}{{r}_{{l}_{m+1}}-{r}_{{l}_{m}}}{\left|{\Theta}_{{l}_{m+1}}-{\Theta}_{{l}_{m}}\right|}^{2}+{\displaystyle \underset{j={l}_{m}+1}{\overset{{l}_{m+1}-1}{\sum}}}{r}_{j}^{2}\frac{{r}_{j+1}-{r}_{{l}_{m}}}{{r}_{j}-{r}_{{l}_{m}}}{\left|{\Theta}_{j}-{u}_{{l}_{m},j}\right|}^{2$

with

${u}_{{l}_{m},j}=\frac{{r}_{{l}_{m}}{\Theta}_{{l}_{m}}+\left(j-{l}_{m}\right){r}_{j+1}{\Theta}_{j+1}}{j\left({r}_{j+1}-{r}_{{l}_{m}}\right)},$ (27)

we have

#Math_98# (28)

which is equivalent to Eqution (7).

Submit or recommend next manuscript to SCIRP and we will provide best service for you:

Accepting pre-submission inquiries through Email, Facebook, LinkedIn, Twitter, etc.

A wide selection of journals (inclusive of 9 subjects, more than 200 journals)

Providing 24-hour high-quality service

User-friendly online submission system

Fair and swift peer-review system

Efficient typesetting and proofreading procedure

Display of the result of downloads and visits, as well as the number of cited articles

Maximum dissemination of your research work

Submit your manuscript at: http://papersubmission.scirp.org/

Or contact ijaa@scirp.org