In this article, we investigate the nonparaxial propagation properties of the chirped Airy Gaussian vortex(CAiGV)beams in uniaxial crystals orthogonal to the optical axis analytically and numerically. We discuss how...In this article, we investigate the nonparaxial propagation properties of the chirped Airy Gaussian vortex(CAiGV)beams in uniaxial crystals orthogonal to the optical axis analytically and numerically. We discuss how the linear chirp parameters, the quadratic chirp parameters, and the Gaussian factors influence the nonparaxial propagation dynamics of the CAiGV beams. The intensity, the energy flow, the beam center, and the angular momentum of the CAiGV beams are deeply investigated. It is shown that the Gaussian factors have a great effect on the intensity and the centroid positions of the CAiGV beams. With the Gaussian factors increasing, the intensity of CAiGV beams decreases rapidly. The main lobes of the transverse intensity distribution of the CAiGV beams are similar to triangles.展开更多
The concept of partially coherent nonparaxial modified Bessel Gauss (MBG) beams is proposed. Based on the generalized Rayleigh-Sommerfeld diffraction integral, the analytical propagation equations of nonparaxial MBG...The concept of partially coherent nonparaxial modified Bessel Gauss (MBG) beams is proposed. Based on the generalized Rayleigh-Sommerfeld diffraction integral, the analytical propagation equations of nonparaxial MBG beams in free space are derived and analysed, and some special cases are discussed. In particular, under the paraxial approximation our results reduce to the corresponding paraxial ones. Numerical calculation examples are given to illustrate the dependence of intensity and spectral degree of coherence on the beam order m, ζ and f parameters, and to compare the difference between the paraxial and nonparaxial results.展开更多
Analytical expressions for the three components of the nonparaxial propagation of a Hermite-Laguerre-Gaussian (HLG) beam in uniaxial crystal orthogonal to the optical axis are derived. The intensity distribution of ...Analytical expressions for the three components of the nonparaxial propagation of a Hermite-Laguerre-Gaussian (HLG) beam in uniaxial crystal orthogonal to the optical axis are derived. The intensity distribution of an HLG beam and its three components propagating in a uniaxial crystal orthogonal to the optical axis are demonstrated by numerical examples. Although the y and z components of an HLG beam in the incident plane are both equal to zero, they emerge upon propagation inside the uniaxial crystal. Moreover, the beam profile of the x component is relatively stable and the beam profiles of the y and z components have the same evolution law. If the ratio of the extraordinary refractive index to the ordinary refractive index is larger than unity, the beam profile of the HLG beam is elongated in the x direction and generally rotates clockwise. Otherwise, the beam profile of the HLG beam is elongated in the y direction and generally rotates anticlockwise. This research is beneficial to the optical trapping and nonlinear optics involved in the rotation of a beam profile.展开更多
This paper derives the closed-form expressions for nonparaxial phase flipped Gaussian (PFG) beams propagating in free space, through a knife edge and an aperture, which enable us to study nonparaxial propagation pro...This paper derives the closed-form expressions for nonparaxial phase flipped Gaussian (PFG) beams propagating in free space, through a knife edge and an aperture, which enable us to study nonparaxial propagation properties of PFG beams and to compare nonparaxial results with paraxial ones. It is found that the f parameter, offsetting distance of the knife edge and truncation parameter affect the nonparaxial beam propagation properties. Only under certain conditions the paraxial approximation is applicable. The results are illustrated by numerical examples.展开更多
Based on the vectorial Rayleigh-Sommerfeld integral formulae, this paper derives the analytical nonparaxial propagation equation of a super-Lorentz-Gauss (SLG) SLGo1 mode beam in free space. The far field expression...Based on the vectorial Rayleigh-Sommerfeld integral formulae, this paper derives the analytical nonparaxial propagation equation of a super-Lorentz-Gauss (SLG) SLGo1 mode beam in free space. The far field expression and the scalar paraxial result are treated with special cases of the general formulae. According to the obtained analytical representation, the nonparaxial propagation properties of the SLG01 mode beam are illustrated and analysed with numerical examples. This research provides an approach to investigate the propagation of the SLG01 mode beam within the framework of the nonparaxial regime.展开更多
Based on the vectorial Ftayleigh-Sommerfeld integral formula and the complex Gaussian expansion of the hardedge aperture function, an analytical propagation expression for a nonparaxial vectorial off-axis Lorentz beam...Based on the vectorial Ftayleigh-Sommerfeld integral formula and the complex Gaussian expansion of the hardedge aperture function, an analytical propagation expression for a nonparaxial vectorial off-axis Lorentz beam passing through a rectangular aperture is derived. The unapertured case, the far field expression and the scalar paraxial result are also presented as special cases of the general formulae, respectively. Some numerical examples are also given to show the propagation characteristics of a nonparaxial vectorial off-axis Lorentz beam through a rectangular aperture. It is indicated that the f parameter, the off-axis displacement and the truncation parameter all play an important role in determining nonparaxial propagation behaviour.展开更多
The nonparaxial propagation of radially polarized chirped Airy beams(RPCAiBs) in uniaxial crystal orthogonal to the optical axis is analytically studied. The effects of the first and the second order chirp factors, an...The nonparaxial propagation of radially polarized chirped Airy beams(RPCAiBs) in uniaxial crystal orthogonal to the optical axis is analytically studied. The effects of the first and the second order chirp factors, and the ratio of the extraordinary refractive index to the ordinary refractive index on the nonparaxial evolution of RPCAi Bs in different observation planes are investigated in detail. The results show that when one parameter changes, different components behave differently, and even the same component has various behavior in different directions. The initial slope of RPCAi Bs in the x-direction varies more with the first order chirp factor than that in the y-direction. Meanwhile, with the second order chirp factor becoming larger, the depth of the focus of the y-component in the x-direction decreases while that in the y-direction has no difference. In addition, the different ratios of the extraordinary refractive index to the ordinary refractive index change the difference of the depth of the focus between the x-and the y-components.展开更多
Starting from the Rayleigh-Sommerfeld diffraction integral, this paper studies the spectral behavior in Young's experiment illuminated by nonparaxial partially coherent light and compares with the paraxial case, wher...Starting from the Rayleigh-Sommerfeld diffraction integral, this paper studies the spectral behavior in Young's experiment illuminated by nonparaxial partially coherent light and compares with the paraxial case, where the influence of nonparaxiality of partially coherent light on the spectral shifts and spectral switches is stressed. It is shown that there is a spectral shift in the nonparaxial case relative to the paraxial one and the critical position changes, at which the spectral switch occurs. The ratio of the waist width to the central wavelength ω0/λ0 and relative spatial correlation length △ affect the spectral difference. The smaller ω0/λ0 is, the larger the difference between the nonparaxial and paraxial results appears. The effect of relative spatial correlation length △ is relatively small.展开更多
A new type of beam called the NPEP-BG(a Bessel-Gaussian(BG)beam with a spiral phase term of a new powerexponent-phase(NPEP))is theoretically studied in this paper.The results show that the number of singularities of t...A new type of beam called the NPEP-BG(a Bessel-Gaussian(BG)beam with a spiral phase term of a new powerexponent-phase(NPEP))is theoretically studied in this paper.The results show that the number of singularities of the phase and side lobes of the intensity of this beam were equal to the topological charges(TCs)and the beam has the characteristics of self-healing during propagation.The NPEP-BG beam combined the partial characteristics of the new power-exponentphase vortex(NPEPV)and the Bessel-Gaussian beam.At the focus of the beam,the characteristic like a perfect vortex beam was present,in which the focal radius is stable and independent of the topological charge.There are multi-focal spots around a ring in the focal plane,meaning that the NPEP-BG beam has a potential for application in multi-particle manipulation.展开更多
Based on the integral representation of the Bessel functions and the generating function of the Tricomi function, an analytical expression of the Wigner distribution function (WDF) for a coherent or partially cohere...Based on the integral representation of the Bessel functions and the generating function of the Tricomi function, an analytical expression of the Wigner distribution function (WDF) for a coherent or partially coherent Bessel Gaussian beam is presented. The reduced two-dimensional WDFs are also demonstrated graphically, which reveals the dependence of the reduced WDFs on the beam parameters.展开更多
The direct acceleration of electrons by using two linearly polarized crossed Bessel-Gaussian (BG) beams with equal frequency and amplitude in vacuum is proposed and studied. It is shown that two linearly polarized B...The direct acceleration of electrons by using two linearly polarized crossed Bessel-Gaussian (BG) beams with equal frequency and amplitude in vacuum is proposed and studied. It is shown that two linearly polarized BG beams of the same order (0 or 1) with a π-rad phase difference have a resultant non-zero longitudinal electric field on the z-axis and can be used, in principle, to accelerate electrons.展开更多
The formalism of generalized diffraction integral for paraxial misaligned optical systems is used to investigate the propagation of the Modified Bessel-Gaussian (MBG) beam through a misaligned thin lens. The propertie...The formalism of generalized diffraction integral for paraxial misaligned optical systems is used to investigate the propagation of the Modified Bessel-Gaussian (MBG) beam through a misaligned thin lens. The properties of the propagation of MBG beam traveling through this misaligned ABCD optical system are discussed. A special case of misaligned circular thin lens is illustrated analytically and numerically. The shape of the MBG beam at the exit of the misaligned optical system is unchanged;however the center of the beam is shifted from the propagation axis in correlated manner with the design parameters of the optical system.展开更多
Starting from the vectorial Rayleigh-Sommerfeld integrals, the free-space propagation expressions for vectorial Hermit-Lagucrre-Gaussian (HLG) beams beyond the paraxial approximation are derived. The far-field expre...Starting from the vectorial Rayleigh-Sommerfeld integrals, the free-space propagation expressions for vectorial Hermit-Lagucrre-Gaussian (HLG) beams beyond the paraxial approximation are derived. The far-field expressions and the scalar paraxial results are given as special cases of our general expressions. The intensity distributions of vectorial nonparaxial HLG beams are studied and illustrated with numerical examples.展开更多
To reduce the walk-off angle of the extraordinary third-harmonic ultraviolet wave at 355 nm generated by type ⅡKTiOPO4 and type I β-BaB2O4 optical crystals,and the Gaussian output beam of a Q-switched Nd:YAG laser,...To reduce the walk-off angle of the extraordinary third-harmonic ultraviolet wave at 355 nm generated by type ⅡKTiOPO4 and type I β-BaB2O4 optical crystals,and the Gaussian output beam of a Q-switched Nd:YAG laser,a simple theoretical model was developed based on a rotatable BK7 plate of variable thickness.By rotating the plate up to 35°along the beam direction,we reduced the walk-off angle up to ~ 13%.The same phenomenon is predicted by the model,confirming the performance of the model.It is found that,due to the walk-off effect,the intensity profile of the thirdharmonic generation beam is slightly degraded.To compensate for the observed phenomena and further reduce the walkoff,we used a combination of a convex lens and an axicon to transform the beam profile of the interacting fundamental and second-harmonic generation waves to the zero-order Bessel-Gaussian form.As a result,the walk-off is decreased to~ 48.81 mrad,providing ~ 30%relative reduction.By using the same BK7 plate rotated up to 35° along the third-harmonic beam direction,the walk-off angle is further reduced to 38.9 mrad.Moreover,it is observed that the beam profile of the emerged Bessel-Gaussian third-harmonic generation beam remains unchanged with no degradation.展开更多
In this study,we report on the fractional Talbot ffect of nonpar axial self-accelerating beams in a multilevel electromagnetically induced transparency(EIT)atomic configuration,which,to the best of our knowledge,is th...In this study,we report on the fractional Talbot ffect of nonpar axial self-accelerating beams in a multilevel electromagnetically induced transparency(EIT)atomic configuration,which,to the best of our knowledge,is the first study on this subject.The Tallbot ffect originates from superposed eigenmodes of the Helmholtz equation and forms in the EIT window in the presence of both linear and cubic susceptibilities.The Talbot ffect can be realized by appropriately selecting the cofficients of the beam components.Our results indicate that the larger the radial difference between beam components,the stronger the interference between them,the smaller the Tallbot angle is.The results of this study can be useful when studying optical imaging.optical measurements,and optical computing.展开更多
The nonparaxial property of the chirped pulsed beam is analyzed both quantitatively and qualitatively. Through the qualitative investigation of the paraxial approximation condition, we show there are chirp-induced cha...The nonparaxial property of the chirped pulsed beam is analyzed both quantitatively and qualitatively. Through the qualitative investigation of the paraxial approximation condition, we show there are chirp-induced changes in the nonparaxial propagation of the chirped pulsed beam. A quantitative nonparaxial correction was developed by use of the perturbational technic and the Fourier transform for a few-cycle chirped pulsed beam with relative small chirp parameter. It was shown that the nonparaxial corrections were enhanced near the leading or trailing edge of pulse depending on weather the chirp parameter is positive or negative. An example for pulsed Gaussian beam driven by a chirped Gaussian pulse is shown in the numerical result to confirm our analysis.展开更多
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11775083 and 11374108)
文摘In this article, we investigate the nonparaxial propagation properties of the chirped Airy Gaussian vortex(CAiGV)beams in uniaxial crystals orthogonal to the optical axis analytically and numerically. We discuss how the linear chirp parameters, the quadratic chirp parameters, and the Gaussian factors influence the nonparaxial propagation dynamics of the CAiGV beams. The intensity, the energy flow, the beam center, and the angular momentum of the CAiGV beams are deeply investigated. It is shown that the Gaussian factors have a great effect on the intensity and the centroid positions of the CAiGV beams. With the Gaussian factors increasing, the intensity of CAiGV beams decreases rapidly. The main lobes of the transverse intensity distribution of the CAiGV beams are similar to triangles.
基金Project supported by the National High Technology Development Program of China (Grant No 823070) and the National Natural Science Foundation of China (Grant No 10574097).
文摘The concept of partially coherent nonparaxial modified Bessel Gauss (MBG) beams is proposed. Based on the generalized Rayleigh-Sommerfeld diffraction integral, the analytical propagation equations of nonparaxial MBG beams in free space are derived and analysed, and some special cases are discussed. In particular, under the paraxial approximation our results reduce to the corresponding paraxial ones. Numerical calculation examples are given to illustrate the dependence of intensity and spectral degree of coherence on the beam order m, ζ and f parameters, and to compare the difference between the paraxial and nonparaxial results.
基金supported by the National Natural Science Foundation of China (Grant Nos. 10974179 and 61178016)the Natural Science Foundation of Zhejiang Province,China (Grant No. Y1090073)the Key Project of the Education Commission of Zhejiang Province,China (Grant No. Z201120128)
文摘Analytical expressions for the three components of the nonparaxial propagation of a Hermite-Laguerre-Gaussian (HLG) beam in uniaxial crystal orthogonal to the optical axis are derived. The intensity distribution of an HLG beam and its three components propagating in a uniaxial crystal orthogonal to the optical axis are demonstrated by numerical examples. Although the y and z components of an HLG beam in the incident plane are both equal to zero, they emerge upon propagation inside the uniaxial crystal. Moreover, the beam profile of the x component is relatively stable and the beam profiles of the y and z components have the same evolution law. If the ratio of the extraordinary refractive index to the ordinary refractive index is larger than unity, the beam profile of the HLG beam is elongated in the x direction and generally rotates clockwise. Otherwise, the beam profile of the HLG beam is elongated in the y direction and generally rotates anticlockwise. This research is beneficial to the optical trapping and nonlinear optics involved in the rotation of a beam profile.
基金Project supported by the National Natural Science Foundation of China(Grant No10574097)
文摘This paper derives the closed-form expressions for nonparaxial phase flipped Gaussian (PFG) beams propagating in free space, through a knife edge and an aperture, which enable us to study nonparaxial propagation properties of PFG beams and to compare nonparaxial results with paraxial ones. It is found that the f parameter, offsetting distance of the knife edge and truncation parameter affect the nonparaxial beam propagation properties. Only under certain conditions the paraxial approximation is applicable. The results are illustrated by numerical examples.
基金Project supported by the National Natural Science Foundation of China (Grant No.10974179)the Natural Science Foundation of Zhejiang Province of China (Grant No.Y1090073)
文摘Based on the vectorial Rayleigh-Sommerfeld integral formulae, this paper derives the analytical nonparaxial propagation equation of a super-Lorentz-Gauss (SLG) SLGo1 mode beam in free space. The far field expression and the scalar paraxial result are treated with special cases of the general formulae. According to the obtained analytical representation, the nonparaxial propagation properties of the SLG01 mode beam are illustrated and analysed with numerical examples. This research provides an approach to investigate the propagation of the SLG01 mode beam within the framework of the nonparaxial regime.
基金supported by the Scientific Research Fund of Zhejiang Provincial Education Department,China(Grant No 20060677)
文摘Based on the vectorial Ftayleigh-Sommerfeld integral formula and the complex Gaussian expansion of the hardedge aperture function, an analytical propagation expression for a nonparaxial vectorial off-axis Lorentz beam passing through a rectangular aperture is derived. The unapertured case, the far field expression and the scalar paraxial result are also presented as special cases of the general formulae, respectively. Some numerical examples are also given to show the propagation characteristics of a nonparaxial vectorial off-axis Lorentz beam through a rectangular aperture. It is indicated that the f parameter, the off-axis displacement and the truncation parameter all play an important role in determining nonparaxial propagation behaviour.
基金Projected supported by the National Natural Science Foundation of China(Grant Nos.11775083 and 11374108)the Science and Technology Program of Guangzhou Province,China(Grant No.2019050001)the Special Funds for the Cultivation of Guangdong College Students’Scientific and Technological Innovation,China(Grant No.pdjh2020a0149)。
文摘The nonparaxial propagation of radially polarized chirped Airy beams(RPCAiBs) in uniaxial crystal orthogonal to the optical axis is analytically studied. The effects of the first and the second order chirp factors, and the ratio of the extraordinary refractive index to the ordinary refractive index on the nonparaxial evolution of RPCAi Bs in different observation planes are investigated in detail. The results show that when one parameter changes, different components behave differently, and even the same component has various behavior in different directions. The initial slope of RPCAi Bs in the x-direction varies more with the first order chirp factor than that in the y-direction. Meanwhile, with the second order chirp factor becoming larger, the depth of the focus of the y-component in the x-direction decreases while that in the y-direction has no difference. In addition, the different ratios of the extraordinary refractive index to the ordinary refractive index change the difference of the depth of the focus between the x-and the y-components.
基金suported by the National Natural Science Foundation of China (Grant No 10574097)the Outstanding Young Researcher Foundation of Sichuan Province,China (Grant No 03ZQ026-061)+1 种基金the Applied and Basic Research Foundation of Sichuan Province,China (Grant No 05J Y029-102)the Natural Science Foundation of Yibin University,China (Grant No 2008B04)
文摘Starting from the Rayleigh-Sommerfeld diffraction integral, this paper studies the spectral behavior in Young's experiment illuminated by nonparaxial partially coherent light and compares with the paraxial case, where the influence of nonparaxiality of partially coherent light on the spectral shifts and spectral switches is stressed. It is shown that there is a spectral shift in the nonparaxial case relative to the paraxial one and the critical position changes, at which the spectral switch occurs. The ratio of the waist width to the central wavelength ω0/λ0 and relative spatial correlation length △ affect the spectral difference. The smaller ω0/λ0 is, the larger the difference between the nonparaxial and paraxial results appears. The effect of relative spatial correlation length △ is relatively small.
基金Project supported by the Program of Shanghai Academic/Technology Research Leader(Grant No.20SR014501)the Strategic Priority Research Program of the Chinese Academy of Sciences(Grant No.XDB1603)the Youth Innovation Promotion Association of the Chinese Academy of Sciences(Grant No.2019247)。
文摘A new type of beam called the NPEP-BG(a Bessel-Gaussian(BG)beam with a spiral phase term of a new powerexponent-phase(NPEP))is theoretically studied in this paper.The results show that the number of singularities of the phase and side lobes of the intensity of this beam were equal to the topological charges(TCs)and the beam has the characteristics of self-healing during propagation.The NPEP-BG beam combined the partial characteristics of the new power-exponentphase vortex(NPEPV)and the Bessel-Gaussian beam.At the focus of the beam,the characteristic like a perfect vortex beam was present,in which the focal radius is stable and independent of the topological charge.There are multi-focal spots around a ring in the focal plane,meaning that the NPEP-BG beam has a potential for application in multi-particle manipulation.
文摘Based on the integral representation of the Bessel functions and the generating function of the Tricomi function, an analytical expression of the Wigner distribution function (WDF) for a coherent or partially coherent Bessel Gaussian beam is presented. The reduced two-dimensional WDFs are also demonstrated graphically, which reveals the dependence of the reduced WDFs on the beam parameters.
基金Project supported by the National Natural Science Foundation of China (Grant No 10574097).
文摘The direct acceleration of electrons by using two linearly polarized crossed Bessel-Gaussian (BG) beams with equal frequency and amplitude in vacuum is proposed and studied. It is shown that two linearly polarized BG beams of the same order (0 or 1) with a π-rad phase difference have a resultant non-zero longitudinal electric field on the z-axis and can be used, in principle, to accelerate electrons.
文摘The formalism of generalized diffraction integral for paraxial misaligned optical systems is used to investigate the propagation of the Modified Bessel-Gaussian (MBG) beam through a misaligned thin lens. The properties of the propagation of MBG beam traveling through this misaligned ABCD optical system are discussed. A special case of misaligned circular thin lens is illustrated analytically and numerically. The shape of the MBG beam at the exit of the misaligned optical system is unchanged;however the center of the beam is shifted from the propagation axis in correlated manner with the design parameters of the optical system.
基金Project supported by the National Natural Science Foundation of China (Grant No 10574097).
文摘Starting from the vectorial Rayleigh-Sommerfeld integrals, the free-space propagation expressions for vectorial Hermit-Lagucrre-Gaussian (HLG) beams beyond the paraxial approximation are derived. The far-field expressions and the scalar paraxial results are given as special cases of our general expressions. The intensity distributions of vectorial nonparaxial HLG beams are studied and illustrated with numerical examples.
文摘To reduce the walk-off angle of the extraordinary third-harmonic ultraviolet wave at 355 nm generated by type ⅡKTiOPO4 and type I β-BaB2O4 optical crystals,and the Gaussian output beam of a Q-switched Nd:YAG laser,a simple theoretical model was developed based on a rotatable BK7 plate of variable thickness.By rotating the plate up to 35°along the beam direction,we reduced the walk-off angle up to ~ 13%.The same phenomenon is predicted by the model,confirming the performance of the model.It is found that,due to the walk-off effect,the intensity profile of the thirdharmonic generation beam is slightly degraded.To compensate for the observed phenomena and further reduce the walkoff,we used a combination of a convex lens and an axicon to transform the beam profile of the interacting fundamental and second-harmonic generation waves to the zero-order Bessel-Gaussian form.As a result,the walk-off is decreased to~ 48.81 mrad,providing ~ 30%relative reduction.By using the same BK7 plate rotated up to 35° along the third-harmonic beam direction,the walk-off angle is further reduced to 38.9 mrad.Moreover,it is observed that the beam profile of the emerged Bessel-Gaussian third-harmonic generation beam remains unchanged with no degradation.
基金the National Natural Science Foundation of China(Grant Nos.61805068,11747046,and 61875053)China Postdoctoral Science Foundation(CPSF)(Grant No.2017M620300)+1 种基金the Natural Science Fund of Shaanxi Province(Grant No.2018JQ6002)the Science and Technology Department of Henan Province(Grant No.202102210111).
文摘In this study,we report on the fractional Talbot ffect of nonpar axial self-accelerating beams in a multilevel electromagnetically induced transparency(EIT)atomic configuration,which,to the best of our knowledge,is the first study on this subject.The Tallbot ffect originates from superposed eigenmodes of the Helmholtz equation and forms in the EIT window in the presence of both linear and cubic susceptibilities.The Talbot ffect can be realized by appropriately selecting the cofficients of the beam components.Our results indicate that the larger the radial difference between beam components,the stronger the interference between them,the smaller the Tallbot angle is.The results of this study can be useful when studying optical imaging.optical measurements,and optical computing.
基金This work was supported by the National Natural Sci-ence Foundation of China under(No.60278013). the team Project of Natural Science Foundation of Guang-dong Province(No.20003061), the Foundation of National High Technology Development Program
文摘The nonparaxial property of the chirped pulsed beam is analyzed both quantitatively and qualitatively. Through the qualitative investigation of the paraxial approximation condition, we show there are chirp-induced changes in the nonparaxial propagation of the chirped pulsed beam. A quantitative nonparaxial correction was developed by use of the perturbational technic and the Fourier transform for a few-cycle chirped pulsed beam with relative small chirp parameter. It was shown that the nonparaxial corrections were enhanced near the leading or trailing edge of pulse depending on weather the chirp parameter is positive or negative. An example for pulsed Gaussian beam driven by a chirped Gaussian pulse is shown in the numerical result to confirm our analysis.