Boson star with parity-odd symmetry in wormhole spacetime

In this paper, we revisit the model of bosonic matter in the form of a free complex scalar field with a nontrivial wormhole spacetime topology supported by a free phantom field. We obtain a new type of boson star with wormhole solutions, in which the complex scalar field possess full parity-odd symmetry with respect to the two asymptotically flat spacetime regions. When the size of the throat is small, the behavior of boson stars with wormhole approaches that of boson stars. When the size of the throat is intermediate, the typical spiraling dependence of the mass and the particle number on the frequency of the boson stars is replaced by a loop structure. However, as the size becomes relatively large, the loop structure will also disappear. In particular, The complex scalar field could form two boson stars with opposite phase differences with respect to the two spacetime regions in the limit of vanishing throat size. We analyze the properties of this new type of boson stars with wormhole and further show that the wormhole spacetime geometry.


I. INTRODUCTION
Wormholes are fascinating hypothetical structures in space that have captured the imagination of scientists and science-fiction writers alike.The concept of a wormhole was first introduced in 1916 by the German physicist Ludwig Flamm.However, it wasn't until 1935 that A. Einstein and N. Rosen developed the theory of wormholes as part of their work on general relativity, also known as an Einstein-Rosen bridge [1], which is a theoretical tunnel or shortcut that connects two separate points in spacetime.These tunnels can be used to travel vast distances through space or even time, bypassing the need to traverse the vast distances of the universe at sub-light speeds.However, such a Einstein-Rosen wormhole would not be traversable for matter or for light [2,3].Later, M. Morris and K. Thorne [4] found that creating or finding a humanly traversable wormhole would require the existence of exotic matter with negative energy density, which violates null energy conditions (NEC) [5].Several types of such traversable wormholes have actually been discovered in [6][7][8][9], and a phantom scalar field is selected as the potential candidate for such exotic matter, which has a reversed sign in its kinetic term.
Previously, most of the research on boson stars focused on topologically trivial spacetime.In [37,38], the investigation focused on mixed systems that included boson stars with wormholes in the core of the star.The bosonic matter that comprised these systems consisted of an ordinary complex boson field with self-interaction, which could exhibit symmetric or asymmetrical distribution across the two asymptotically flat regions.As a result, the wormhole spacetime could be either symmetric or asymmetric.Furthermore, the rotation of the scalar field induces a symmetric and asymmetric rotating wormhole spacetime as well [39,40].In these works, both the symmetric and asymmetric wormhole solutions the complex scalar field exhibit the same sign difference with respect to the two asymptotically flat spacetime regions connected by the wormhole.Thus, it is an interesting question whether there exist complex scalar fields with opposite sign differences at the two ends of a wormhole.In this paper, we reinvestigate the model of bosonic matter in the form of a free complex scalar field with a nontrivial wormhole spacetime topology supported by a free phantom field.A new kind of boson star with wormhole solutions is obtained, in which the complex scalar field could possess full parity-odd symmetry with respect to the two asymptotically flat spacetime regions.In particular, the scalar field that form two boson stars with respect to the two spacetime regions has opposite phase differences in the limit of vanishing throat size.It is noteworthy that this new kind of solution still possesses symmetric wormhole spacetime structure.
The paper is organized as follows.In Sec.II , we review the model of a complex scalar field and a phantom scalar field field coupled with gravity.We present the boundary conditions and quantities of interest in Sec.III, discuss the numerical results in Sec.IV.The last section is devoted to conclusion.

II. THE MODEL
We consider the Einstein-Hilbert action including the Lagrangian for the complex scalar field and the phantom scalar field field where R is the Ricci scalar.The term L p and L m are the Lagrangians defined by with Here Ψ and Φ represent the complex scalar field and the phantom field, respectively.By varying the action (1) with respect to the metric, we can obtain the Einstein equations with stress-energy tensor and the matter field equations by varying with respect to the scalar field and phantom field and We consider the general static spherically symmetric solution with a wormhole, and adopt the Ansätzes as follows, see e.g.[38]: where A and B are functions of radial coordinate r, h = r 2 + r 2 0 with the throat parameter r 0 .and r ranges from positive infinity to negative infinity.The two limits r → ±∞ correspond to two distinct asymptotically flat spacetime.Furthermore, we assume stationary complex scalar field and phantom field in the form Here, ψ is only a real function of the radial coordinate r, and the constant ω is referred to as the synchronization frequency.Moreover, the phantom field Φ is also a real function and is independent of the time coordinate t.
Substituting the above Ansätzes ( 8) and ( 9) into the Einstein equations ( 4) leads to the following field equations Variation of the action with respect to the complex scalar field and to the phantom field leads to the equations These five equations are divided into three groups: three of these equations ( 10), ( 11) and ( 13) are solved together, yielding a coupled system of three ODEs on the three unknown functions A, B and ψ.After obtaining the solution for these three functions, ones can solve Eq. ( 14) to determine the value of φ.The remaining one equation ( 12) is treated as the constraint and used to check the numerical accuracy of the method.Especially, The above derivative in Eq. ( 14) happens to be zero, which means it can be equal to a constant.Integrating the last equation, ones can leads to Here the constant of integration D denotes the scalar charge of the phantom field.By taking Eq. ( 15) into Eq.( 12) the scalar charge D can be expressed into the following form The system of coupled non-linear equations ( 10), ( 11) and ( 13) presents a formidable challenge in seeking a general analytical solution, necessitating the use of numerical methods to obtain the solutions to differential equations.However, prior to numerically solving the aforementioned equations, it is worthwhile to investigate two special solution cases.Firstly, when the complex scalar field vanishes, one can derive the solution for an Ellis wormhole.Secondly, as the throat size r 0 approaches zero, two mini-boson stars emerge about the r 0 = 0 symmetry and correspond to two distinct asymptotically flat spacetimes.These two special cases provide valuable insight into the behavior of the system and can inform subsequent numerical analyses.

III. BOUNDARY CONDITIONS AND QUANTITIES OF INTEREST
To solve the coupled non-linear equations (10), (11) and (13) in asymptotically flat spacetime, we give appropriate boundary conditions for the functions as follows at infinity (r → ∞).At the origin (r = 0), the symmetry of functions can be used to distinguish between two types of solutions.For even parity solutions, we require the following boundary conditions The above boundary conditions were employed in [38] to investigate the scenario of symmetric solutions.Considering that the complex scalar with parity-odd symmetry, we impose The ADM mass M is the key quantities we are interested in, which is encoded in the asymptotic expansion of metric components The action of the complex scalar field is invariant under the U (1) transformation ψ → e iα ψ with a constant α.
According to Noether's theorem, there is a conserved current corresponding to the complex scalar field: Integrating the timelike component of the above conserved currents on a spacelike hypersurface S, ones could obtain Noether charge: The global charge Q is equivalent to the number of particles in the complex scalar field.

IV. NUMERICAL RESULTS
In this work, all the numbers are dimensionless as follows For convenience, we set µ 0 = 1 and κ = 0.2.We transform the radial coordinates by the following equation The use of Eq. ( 24) map the infinite region (−∞, ∞) to the finite region (−1, 1).Next, what we need is to discretize the differential equations on a grid.This allows the ordinary differential equations to be approximated by algebraic equations.The grid with 500 points covers the integration region −1 ≤ x ≤ 1.In our work, relative errors are less than 10 −6 .First, as an example, we plot in Fig. 1 two typical profiles of our numerical results for the scalar field ψ as a function of x.In the left panel, the throat size r 0 = 1 is fixed and the curves with several values of the frequency ω are shown.We can see that the scalar field ψ has an antisymmetric structure about the origin x = 0. ψ max (the maximum value of the field function ψ) increases as the frequency ω decreases.Moreover, the extreme points converge towards the origin x = 0 as the frequency decreases.In the right panel, the frequency ω = 0.95 is fixed and curves for different values of throat size r 0 are displayed.In particular, we represent the two boson star solutions with black dashed lines, and these two solutions form an anti-symmetric structure at the two ends of the wormhole.Ones can see that at positions where x is large, the curves for different values of r 0 do not differ much and basically overlap.In the region where x is small, the value of the scalar field increases as r 0 decreases.When the limit r 0 → 0, the solutions that are very close to two boson stars at the two ends of the wormhole except near the origin.

A. Global charges of the solutions
In this subsection, we will initially examine the global charges of the resulting wormhole solutions.In Fig. 2, we exhibit the mass M and Noether charges Q versus the frequency ω with the different values of throat size r 0 , the black curve corresponds to the boson star solution.It is interesting that the curve does not form a spiral.Instead, When r 0 starts to increase from 0 value, the curve will first form a loop structure.Then, as r 0 continues to increase, the multi-valued curve will become a single-valued curve.In the first branch, when the frequency is fixed, we find that as r 0 increases, both the mass M and charge Q also increase, which have similar behavior.
In Fig. 3, we exhibit the scalar charge D versus the frequency ω with the different values of throat size r 0 .We can see that as r 0 increases, the scalar charge also increases, which means that more phantom matter is needed.For curves with larger r 0 values, the charge decreases monotonically with decreasing frequency ω.

B. Geometric properties of wormhole
To explain that the metric (8) does indeed describe a wormhole, one can make use of geometrical embedding diagrams.It is helpful to consider a two-dimensional space at a fixed time and θ.The resulting two-dimensional spatial hypersurface of the wormhole spacetime can then be embedded in a three-dimensional Euclidean space, where the embedding diagram can be used to visualize the wormhole geometry.This technique allows us to better understand the topology and properties of the wormhole solution.
To characterize the geometry of the throat, we begin by examining the embeddings of planes with θ = π/2.Using cylindrical coordinates (ρ, ϕ, z), The metric on this plane can be expressed by the following formula Comparing the two equations above, we then obtain the expression for ρ and z, ρ(r) = B(r)e −A(r) h(r), z(r) = ± B(r)e −A(r) − dρ dr Here ρ corresponds to the circumferential radius, which corresponds to the radius of a circle located in the equatorial plane and having a constant coordinate r.The function ρ(r) has extreme points, where the first derivative is zero.
When the second derivative of the extreme point is greater than zero, we call this point an throat, which corresponds to a minimal surface.When the second derivative of the extreme point is less than zero, we call this point a equator, which corresponds to a maximal surface.
To illustrated geometric properties of throats, In Fig. 4, we show two dimensional view of the isometric embedding of the equatorial plane of symmetric solutions with throat parameter r 0 = 1 and several values of the complex scalar frequency ω.When the frequency is high, the wormhole has only one throat with a small radius.As the frequency increases, the radius of the throat also increases.At near ω = 0.26, two new throats begin to appear, and the wormhole then has three throats.As the frequency further increases, the radii of the three throats will first increase and then decrease.In particular, the middle throat will become less and less visible, although it always exists, but tends to disappear.In order to get further insight into the geometry of the wormholes , we present two three-dimensional plot in Fig. 5, including one with three throats and two equators with ω = 0.15, and the other with only one throat with ω = 0.4.

V. CONCLUSIONS
In this paper, we reconsidered the model of a complex boson field and a phantom field minimally coupled to Einstein gravity and investigated the domain of existence and physical characteristics of wormhole solutions embedded in a complex bosonic matter field with parity-odd symmetry.We have presented the relationship between the conserved charges and frequency, which depends on the throat size.Furthermore, we have compared our results to those of the single boson star solution.Surprisingly, the curve does not exhibit a spiral shape.Instead, as the throat size r 0 increases from zero, the curve initially forms a loop structure.As r 0 continues to increase, the multi-valued curve gradually becomes a single-valued curve.Unlike the even-symmetric case studied previously, the solutions with a complex odd-symmetric bosonic matter field cannot degenerate into boson stars.
There are some interesting extensions of our work which we plan to investigate in future projects.First, in order to achieve a wormhole solution in general relativity, the presence of exotic matter is required.However, under modified gravity models, it is possible to construct wormhole solutions without the need for exotic matter.Therefore, studying the existence of such solutions under modified gravity is interesting.Furthermore, we are planning to investigate the stability of these solutions.Finally, considering that boson stars have both ground state and excited state solutions, we can explore the excited state solutions in the future.

FIG. 1 :
FIG. 1: The radial distribution of the complex scalar field.Left: The distribution with several values of the frequency ω for r0 = 1.Right: The distribution with several values of the throat size r0 for ω = 0.95.

5 FIG. 2 :
FIG.2:The ADM mass M and Noether charge Q as a function of the frequency ω with several values of the throat size r0.

5 FIG. 3 :FIG. 4 :
FIG.3:The scalar charge D of the phantom field as a function of the frequency ω with several values of the throat size r0.

FIG. 5 :
FIG. 5: Geometric properties of throats: Isometric embeddings of the equatorial plane of symmetric solutions with throat parameter r0 = 1.(a) Three dimensional plot with three throats for ω = 0.15.(b) Three dimensional plot with a single throat for ω = 0.4.