Meson effective mass in the isospin medium in hard-wall AdS/QCD model

We study a mass splitting of the light vector, axial-vector, and pseudoscalar mesons in the isospin medium in the framework of the hard-wall model. We write an effective mass definition for the interacting gauge fields and scalar field introduced in gauge field theory in the bulk of AdS space-time. Relying on holographic duality we obtain a formula for the effective mass of a boundary meson in terms of derivative operator over the extra bulk coordinate. The effective mass found in this way coincides with the one obtained from finding of poles of the two-point correlation function. In order to avoid introducing distinguished infrared boundaries in the quantization formula for the different mesons from the same isotriplet we introduce extra action terms at this boundary, which reduces distinguished values of this boundary to the same value. Profile function solutions and effective mass expressions were found for the in-medium ρ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho $$\end{document}, a1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a_1$$\end{document}, and π\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pi $$\end{document} mesons.

In holographic QCD the two phases of nuclear matter -the confining and deconfining ones -in the dual gravity theory side are described by the different metrics. In the bottom-up approach, when the quark matter is absent, the deconfining phase at the gravity side is described by the Schwarzschild AdS black hole (SAdS BH) metric [17]. The confined phase of the nuclear matter in the low temperature limit in the dual gravity theory is described by the thermal AdS space (tAdS) [40,41]. In [42] it was shown on the background geometry of the dual gravity for the confinement phase containing the quark matter fields as well. This metric is named the thermal charged AdS space (tcAdS) metric, and it can be obtained from the Reissner-Nordstrom black hole (RNAdS BH) metric by taking a zero of the mass of black hole and cutting of the fifth dimension at the infrared (IR) boundary. The deconfinement phase in dual gravity is described by the RNAdS BH metric. There occurs a Hawking-Page transition between these geometries as the confinement-deconfinement phase transition takes place in the dual field theory [43][44][45][46].
One of the simplifications in nuclear matter studies using holography is in considering a zero temperature limit. Another one is taking the quark number densities of the medium to be zero and leaving only a dependence on the isospin chemical potential. In the result of these simplifications the isospin medium is obtained for which the background geometry in the dual theory has no modification and the metric remains an ordinary AdS space [47,48]. Such a model is useful to separate the effects taking place due to isospin from the ones occurring under the influences of other quantities of the dense medium [1,34,38,39]. One of such effects is the splitting of the mass spectra of the meson states from the same isospin triplet in the medium due to their isospin interaction with the non-zero isospin of the medium. This effect was considered by a number of authors within the holographic QCD [30,[34][35][36][47][48][49][50][51]. In [47,48] in the hard-wall model framework the meson mass splitting effect was considered for the triplets of the light mesons ρ, π , and a 1 in the dense nuclear matter due to the isospin interaction in connection with the study of the pion condensation in the isospin medium, which had been carried out in [34]. In this paper we reconsider this effect because of the recent observation of a relation between the definition of the effective mass for a meson in the medium with non-zero isospin and a fixed infrared boundary of the hardwall model. Here we define an effective mass following the standard definition in field theory and show that this definition will require fixing the IR boundary differently for the mesons from the same meson isotriplet in the medium. In other words the splitting of the meson mass leads to the splitting of the IR boundary of the model, since the mass spectrum in this model is directly defined by the value of the IR boundary. It should be noticed that the IR boundary shift in the hard-wall model already is known from Refs. [37,52], where the authors deal with the bulk interaction of the fields, which contribute to the boundary meson mass.
The paper is organized as follows: In Sect. 2 we give a description for the dense and isospin mediums in the hardwall model. In Sect. 3 we reproduce the equation of motion for the ρ mesons in the isospin medium determine the effective mass of these mesons by means of holography. We avoid to distinguish the IR boundary introducing boundary terms. In Sect. 4 a similar analysis was made for the π and a 1 mesons. We discuss the relation between the mass shift and the boundary terms in Sect. 5.

Isospin medium in hard-wall model
Let us present at first a bulk foundation for the boundary isospin medium following earlier work [42,47,48,53]. The bulk field content for the isospin medium will be a simplified version of the one for the dense nuclear medium utilized in Ref. [47]. Since this medium is described by the SU (2) group, the bulk flavor gauge group should be chosen as SU (2) L × SU (2) R which then will be broken to SU (2) V . We introduce in the bulk of AdS space two gauge fields A (L) and A (R) transforming under the (1, 0) and (0, 1) representations of this flavor symmetry group, respectively. The field strength tensor for these fields is with where the flavor index a runs over a = 1, 2, 3 and T a are the generators in the spinor representation of the SU (2) L or SU (2) R group and are expressed in terms of Pauli matrices σ a by As a gauge condition on A z we shall choose the axial gauge, The action for the gravity and for the gauge fields will be written as where = −6/R 2 is the cosmological constant. The AdS/CFT correspondence establishes the following relations between the constants κ 2 , g 2 , the number of colors N c , and the radius R of space-time: 1/g 2 = N c /(4π 2 R) and 1/( M are divided into the background parts L M , R M and the fluctuations l M , r M of these parts: In the holographic description of the nuclear matter the bulk flavor group is U (2) L ×U (2) R and this matter in the boundary QCD is given by the bulk gauge fields L M , R M , while the fluctuations l M , r M in the bulk are needed for the description of the vector and axial-vector mesons in the boundary QCD. The homogeneous and isotropic matter at the boundary is described by the dual bulk fields L M and R M , which are not dependent on the space-time coordinates. Moreover, the only time components L 0 and R 0 of these fields are taken non-zero, because the only flavor diagonal element of them (a = 3) corresponds to the physical quantity of the boundary matter. From these components we compose the vector and axial-vector fields: The boundary value of V 3 0 maps to the isospin chemical potential of the medium u and d quarks 1 [53]. Following [47,48] here we shall deal with the case L a 0 = R a 0 , which means Lagrangian invariance under changing the left and right flavor groups SU (2) L ↔ SU (2) R . Obviously, A 3 0 = 0 for this case. This L ↔ R invariance in the dual boundary theory means that the medium particles, i.e. the nucleons, are required to be in the parity-even states. As is well known, the ground states of the nucleons are parity-even ones, and the first excited state of the nucleon can be a parity-even or a parity-odd state. In the parity-even state the nucleons have less energy than in the parity-odd one [15,16].
The field strength tensor for V 3 0 is a flavor diagonal matrix as well and hence it gets the following form for an Abelian field: In such a way the SU (2) symmetry of the V a M part of the gauge fields is broken down to two U (1) symmetries. Meanwhile, the fluctuations l M , r M remain non-abelian. Then the V 3 M part of the action (2) gets the form The holographic dual of the A (u),(d) V 3 0 field will be the u-and d-quarks of the boundary medium. Imposing the hard-wall cut-off on the bulk radial coordinate z makes these quarks confined ones in the boundary theory. It should be noted that the number densities of quarks thus defined (isoquarks) are zero (see [47]).
In the nuclear matter case the bulk geometry is a thermally charged AdS space (tcAdS) with radius R, which is the non-black brane solution of the Einstein-Maxwell system of equations obtained from (8) [42]: Here The constant q α is related to the number densities Q α of the medium u and d quarks, In the isospin medium case the number densities Q α are taken zero and consequently the metric (9) returns into the metric of ordinary AdS space: The radial coordinate z ranges in the limited area 0 < z ≤ z I R due to the hard-wall cut-off. It can be seen from a comparison of the geometries (9) and (11) in the dual boundary theory that there is an effect of the dense medium on the mesons which is provided by the f (z) function in the bulk metric tensor G M N . For instance, if the kinetic energy of the l M or r M field is E in the presence of medium fields A (α) 0 and is E in the absence of these medium fields, then E and E will be related by the relation E = G 00 /G 00 E = f (z)E. Since in the isospin medium case Q α = 0 and f (z) = 1 there is no modification of the AdS geometry by the medium fields and the effect of the metric (11) on the l M , and r M fields is the same as in free field case.
Solutions to the Maxwell equations obtained from the action (8) have the form In the isospin medium case the vector potential A (α) 0 is constant and is equal to the quark chemical potential μ α with a factor 2π 2 : Here we are going to deal with the medium in the confining phase of the medium u and d quarks, where the fundamental excitations are nucleons but not quarks. So, the solution (13) should be expressed in terms of the chemical potentials of the nucleons. Taking into account the quark content of nucleons as in [47] the isospin chemical potentials of nucleons may be defined as a sum of quark chemical potentials μ P = 2μ u +μ d and μ N = μ u + 2μ d . For isospin matter with two flavors the number densities of the nucleons are zero and V 3 0 equals √ 2π 2 (μ u − μ d ) in the case of the deconfinement phase and in the confinement phase one. Thus, V 3 0 in the dual boundary theory describes the homogeneous and constant isospin background field of a medium which is made of isonucleons.

Effective mass definition by use of holography
In this subsection applying the AdS/CFT correspondence we shall derive an effective mass formula for the vector meson minimally interacting with the external isospin field (14). It is well known that in AdS/QCD models the mesons in the dual boundary theory are described by the fluctuations of the gauge and scalar fields in the bulk. A bulk vector field fluctuation, which corresponds to the boundary vector meson, is composed of the left and right fluctuations l i μ and r i μ : where μ = 0, 1, 2, 3 is for the boundary coordinates and i = 1, 2, 3 is for the SU . As a condition for gauge fixing of the fifth components of the l M and r M fields the axial gauge is chosen, which takes values l i z = r i z = 0. The temporal component of the vector fluctuations will correspond to the fluctuations in the chemical potential (or the number density in the nuclear matter case) of medium particles. Setting this component zero (v i 0 = 0) we shall consider a medium with a constant isospin chemical potential. We suppose the bulk v μ field does not depend on the spatial coordinates x m . The action for the total vector [47,48] where is the field strength tensor. After taking the trace in (16) the action S V gets the form Here the m, n = 1, 2, 3 indices denote the boundary spatial coordinates. The terms in the second square bracket in (17) describe the interaction of the vector field fluctuations with the background field V 3 0 . In the dual boundary theory this interaction corresponds to the interaction of the vector meson with the isospin medium. Thus, the vector meson in the isospin medium in a holographic picture is described by the interacting gauge field theory in the bulk of AdS space.
The complex vector field ρ m is composed of the vector field components v i m as follows: As will be seen below the boundary values of the ρ a m bulk fields correspond to the neutral and charged vector mesons, respectively. For the free bulk gauge field ρ a m we can write the Fourier transformation: where the kinetic energies ω a of the ρ a fields ∂ 0 ∂ 0 ρ a n = ω 2 a ρ a n (a = 0, ±) (20) in the dual boundary theory correspond to the masses of the ρ mesons. Equation (20) is the one defining a mass in the free field theory. In interacting gauge field theory, where the interaction is included by minimal coupling, the mass of the interacting field is defined by means of covariant derivatives instead of ordinary ones in (20). Also, as was mentioned above the background gauge field V 3 0 does not change the background geometry (11). So, the interaction with the V 3 0 field still was not included into Eqs. (19)- (20) and the kinetic energies ω ± in these formulas are the ones for the free ρ a field. As is well known, the vacuum masses of all ρ mesons (charged and neutral ones) are equal. This degeneracy takes place for the masses ω a of the bulk ρ fields in the isospin background V 3 0 as well, since there is no metric change by this background: From the action (17) written in terms of the ρ a fields after taking into account (21) we will obtain the following equations of motion: Explicitly, Eq. (22) will be written as below: which have the same form as the ones obtained in [34] for the pions. These equations of motion were obtained in [48] introducing v i fields minimally interacting with the V 3 0 field. The Lagrangian for this system was constructed by means of the covariant derivative D M as follows:  (24) the action for this Lagrangian will be reduced to the following form: This action in terms of the ρ fields will split into three independent terms for each ρ a component: M ρ a n ) * (D (a)M ρ a n ).
Here the covariant derivative D This is a five-dimensional D'Alembert equation in curved space-time (11) for the vector field having a zero fifth component (ρ 5 = 0) and zero five-dimensional mass (m 5 = 0). Boundary terms which arise on obtaining (27) lead to the boundary conditions ∂ z ρ a n | z I R = 0 and ρ a n (z)| ε = 0, which are the same as the ones in the free field case. The explicit form of Eq. (27) is Eq. (22).
It is useful to rewrite Eq. (27) as a sum of its fourdimensional part and the z coordinate part where we divided it by √ −G. Now, let us write the D'Alembert equation in a four-dimensional curved spacetime for the massive vector field ρ a n minimally coupled with the constant external background gauge field V a μ : where the m * a is the mass of the ρ a n field. In our case the four-dimensional space-time is the UV boundary of the AdS space (11) and the massive vector field is the boundary value of the ρ field. In the AdS/CFT correspondence here the UV values of the ρ a fields correspond to the ρ mesons defined on this boundary and the equations of motion for these fields at the UV boundary correspond to the equations of motion for these mesons. Hence, we have to correlate Eq. (28), the equations in the UV limit, to the ones in (29). From this correspondence it is seen that the eigenvalues of the − operator correspond to the squared mass (m * a ) 2 of the ρ a vector fields in (29). So, we may accept the equality of the eigenvalues of this operator to the squared masses of boundary meson fields: Equation (30) is the determining formula of the effective mass for the vector field in the isospin medium by means of the derivative operator over the extra dimension.
In the AdS space geometry (11) the equations of motion (28) have got the form which is the same as (23). The operator D (a) μ D (a)μ has the following eigenvalues: In four-dimensional space-time, where Eq. (29) was written, the eigenvalues (32) determine the effective masses of the massive ρ a fields: according to the effective mass definition of vector field in the constant background field in four-dimensional gauge field theory. The two definitions of the effective mass (30) and (33) give same result and thus, the masses m * a are the masses of the ρ a mesons in the isospin medium.

Effective mass as an eigenvalue of profile function
Another way of determining the effective mass of mesons is using profile function solution. From the two-point correlation function study it is well known that the eigenvalues in the profile function are the mass spectrum of the vector meson in the vacuum case [54]. It is natural to extend this study to the constant isospin medium case. To this end let us write a solution to Eq. (23) for the ρ s a mode in the Kaluza-Klein Here m s a denotes the square root of the bracket in Eq. (23) for the s mode: and ω s a is the mass of the free s mode. This solution is the same as the one found in [54] for the free vector field V μ , except for the constant eigenvalues m s a in the arguments of the Bessel function, which were denoted by M n in the free vector field case. Of course, the distinction between the M n and m s a constants will not change the procedure of finding the poles by means of the two-point function Here ( p 2 ) was defined by means of the profile function V ( p, z) of the vector field as follows: In the free field case V ( p, z) has the form of ρ a s (ω, (34) and for this solution ( p 2 ) gets the form [54] ( Applying a Kneser-Sommerfeld expansion to Eq. (38) in [54] an explicit pole structure of the ( p 2 ) was found: where γ s 0 are the zeros of the Bessel function J 0 (z). For the solution (34) the ( p 2 ) in (39)  z ∂ z ρ a n operators for the s mode according to (33) and (31).
Thus, two ways of determination of the effective masses give the same result and the effective mass of the dual ρ mesons in the isospin medium are the following ones: Although the vacuum masses ω a of ρ mesons are equal, their in-medium masses m * a are different. The splitting in the mass spectrum of ρ mesons occur due to interaction with the isospin medium. After taking into account in (40) the equality of the vacuum mass (ω ± = ω 0 ), the splitting formula will look like The mass splitting formula (41) shows a decreasing of the mass for the positively charged meson and increasing it for the negatively charged one when the positive difference of chemical potentials μ P − μ N grows, and vise versa for the negative μ P − μ N . The isospin medium eliminates the isospin degeneration the of mass, and degeneration occurs when the chemical potentials of the medium protons and neutrons are equal. The splitting of the meson mass in the holographic approach was studied in [34,35,47,48], but the sign of the splitting obtained (41) is different from the one obtained in [34,47,48], and it agrees with the one obtained in [35,55] if we take the value of V 3 0 in (41) equal to the μ I introduced in [35]. It is obvious that the effective mass formula (42) does not obey the holographic definition (30). This distinction between Eqs. (41) and (42) is connected with the definition of the effective mass and will be clarified in the next subsection. The sign of the splitting has importance for meson condensation [55]. According to (41) the positively charged meson becomes massless at a positive value of the isospin medium, it condenses at the value V 3 0 = ω + , while a condensation of the negatively charged meson takes place at ω − = −V 3 0 .

Mass spectrum and boundary terms
Applying the Neumann boundary condition at the IR boundary ∂ z ρ a (z) | z I R = 0 on solution (34) yields the following formula [52] for the mass spectra m s a : which has same form as one in the vacuum and corresponds to mass spectra of excited states of the ρ mesons in the medium.
Here we have distinguished the values of the IR boundary for the differently charged meson, denoting them z 0 I R and z ± I R , respectively. Otherwise, the equality of the right hand sides of the first and second equalities in (43), would mean the equality of their left hand sides also, (ω s 0 ∓ √ 2π 2 (μ P −μ N ) = ω s 0 ), which is nonsense in the non-zero medium isospin case. However, when we introduce different z I R for different components, we get an interaction of fields defined in different space-times, which is not appropriate as well. In order to solve this dilemma we can add to the action new boundary terms, which will lead to the same Kaluza-Klein spectrum (43) for all ρ a components and therefore will ensure fixing z ± I R at the same value with z 0 I R . Obviously, new boundary terms should have a form not changing the equation of motion (22) and (23). We introduce such a kind of boundary terms at z = z I R in the following form: In terms of the v μ fields, S + and S − have the same expression with opposite sign: Then the total action for the ρ a fields will have the form where L (a) is the Lagrangian term producing the equation of motion (22) for the a isocomponent: When we derive the equations of motion (22) from the (47) the boundary terms arising on integration by parts have the same form (ρ a ) * ∂ z ρ a | z I R z U V for the different ρ a components. These boundary terms will be summed with the S ± boundary terms and the total boundary condition imposed on the ρ a solutions at the IR boundary ultimately will be written as below: All three boundary conditions give the same quantization formula -the same spectrum (43) but for the vacuum masses ω s 0 of the in-medium ρ mesons: and the same z I R will be chosen in all three quantization formulas (49). Thus, we conclude that only the vacuum mass part ω s 0 of the effective mass of the meson is subject to quantization. Now we can explain why there are two ways for the definition of the meson mass splitting (Eqs. (41) and (42)). The splitting formula (42) was obtained in the hard-wall model, which has one IR boundary for all isocomponents of the vector field. In this case all three equations in (23) were solved under the same boundary condition at IR and with the same value of z I R . This dictates that all coefficients in front of the last terms in Eq. (23) should be the same. This leads to the splitting formula (42). In fact, if we put back (42) into (23) we shall get the same expression for these coefficients, which equals ω 0 , i.e. we have the same three equations for ρ ± and ρ 0 . Obviously, in this case all m s a coefficients in the Bessel function solutions are the same, m s a = ω 0 , and the necessity does not arise to distinguish z I R in the mass spectrum formula (43). In (44) we obtained this spectrum formula introducing boundary terms, but the effective mass in our case obeys the defining relations (30) and (33).
The first state (s = 1) of the charged vector field tower (40) is the ground state of the ρ ± meson, which has a mass and for the neutral ρ meson we have the value of the mass in vacuum [52]: Finally, the normalized solution (34) for the ρ a meson attains the same form as the one in the vacuum case [52]:

Mass splitting formula for the a 1 and π mesons
In this section we shall define the effective mass for the axialvector and pion fields using the holographic correspondence as was done for the ρ mesons. Let us briefly present here the action and the equations of motion for the axial-vector and pion fields, all details of which can be found in [47,48]. The non-zero fluctuations of the axial-vector field, which is defined as were introduced into the model in order to describe the axial-vector mesons in the boundary QCD. In addition, a complex scalar field is introduced, which performs the chiral symmetry breaking SU (2) L × SU (2) R → SU (2) V of the model. The action for the field has the form Here D M = ∂ M − i L M + i R M and m 2 5 = −3. The field is written in the form = φ exp[i √ 2π a T a ]. The fluctuations π i describe the pions in the dual theory. In the isospin medium case, as was mentioned, the background geometry is given by (11), and thus the solution for the φ field has the simple form of the one in the vacuum case: where m q is the mass of the light quarks and σ is the value of the condensate. The action (53) contains an interaction of field with the axial-vector field a i μ and the total action of these fields will be the sum of S φ and the action S a for the fluctuations a i μ . The action S a , which describes the axialvector field interacting with the gauge field V 3 0 is obtained from the action (4) and has the form where , was defined as one for the vector field fluctuations v i . Since in the dual QCD theory the axial-vector a 1 mesons are described by the transverse field we divide the bulk axial-vector field a i μ into the transverseā i μ and the longitudinal χ parts: a i μ =ā i μ + ∂ μ χ i (∂ μā μ = 0). The longitudinal parts of the left and right fluctuations l i μ and r i μ are equal to + 1 2 ∂ μ χ i and − 1 2 ∂ μ χ i , respectively, and the vector field fluctuations v i μ turn out to be transverse. In an equivalent way an action for the transverse axial-vector fluctuationsā i μ can be obtained from the action (4) and has the form [47] The two actions (55) and (56) lead to the same equations of motion for theā field, as in the vector field case. It is recommended to incorporate the action for the longitudinal χ field with the action for the pseudoscalar field π derived from (53). In summary the action of χ and π fields has the following explicit form [47]: This action includes the interaction term of theā i field with the φ field, which is coming from the action (53) and obviously, this term will be present in the equation of motion for the a 1 field.

The a 1 mesons
Let us introduce the neutral a 0 1m and the charged a ± 1m components for the axial-vector fieldā m , which will correspond to the respective axial-vector a 1 meson in the holography: The Fourier transforms of these components are whereω 0 andω ± are the masses of free a 1 fields. The action for the a a 1 field will be written in the form M a a n 1 ) * × (D (a)M a a 1 n ) + 4g 2 φ 2 a a n 1 a a 1 n }, where the covariant derivative D This is a five-dimensional D'Alembert equation for the axialvector field interacting with the φ(z) field. Obviously, the fifth component and the five-dimensional mass M 5 of this field are taken to be zero. The boundary conditions, which will be imposed on a 1n are obtained from the boundary terms arising on deriving Eq. (61) and are the usual ones: ∂ z a a 1n | z I R = 0 and a a 1n (z)| ε = 0. It will be useful to rewrite Eq. (61) as the sum of its fourdimensional part and the fifth dimension part: On the QCD side there is a massive axial-vector meson, which is described by the a 1n field and has a shift in the mass due to the quark condensate background. We may take into account this mass shift by adding to the QCD Lagrangian an additional mass term proportional to the value of condensate, which will be denoted by 4g 2 φ 2 0 . The constant isospin background of the medium can be taken into account by introducing a constant gauge field V 3 0 . The interaction with this background will be taken into account by the minimal interaction. Let us write the D'Alembert equation for such an a 1n field in a four-dimensional curved space-time, which is, in fact, the UV boundary of the AdS space (11): Herem a is the mass of the a a 1n field not including the condensate shift and φ 0 is the value of the φ(z) field at some fixed value of z. Since the AdS/CFT correspondence of the fields takes place at the UV boundary we take as a fixed value of z the z U V , i.e. φ 0 = φ(z U V ). The condensate interaction term 4g 2 φ 2 0 G mn a a 1n in (63) could be included into the mass term by m 2 a = (m a ) 2 − 4g 2 φ 2 0 . Correlating Eqs. (62) to (63) we observe the equality of the mass term in (62) and the z-dependent part of Eq. (63), Equation (64) is the formula determining the effective mass of the axial-vector meson in the isospin medium by means of the derivative operator over the extra dimension defined in dual bulk theory. The equations of motion for the a a 1 field are obtained from the summary action of (63) and (64) and given by [47] The explicit form of Eq. (65) in the background geometry (11) after the replacement a a 1 (z) = za a 1 (z) will get the following form: Following earlier work [47,48] we consider at first the case of ignoring the condensate term 4g 2 φ 2 in the equations of motion (66). In this case we have no shift of the squared mass and a solution to Eq. (66) for the sth mode a a 1s (z) in the Kaluza-Klein decomposition a a 1n = s f as 1n (P 0 )a a 1s (P 0 , z) after applying UV boundary condition gets a form, which is the same as the one for the vector mesons: The mass spectram s a of the s states arē andω s ± =ω s 0 are the vacuum masses of the excited axialvector mesons. From (68) the mass splitting formula for the a ± 1 mesons will be written in the following form: It is obvious that the mass spectrum formula (43) is available for the axial-vector field case as well. The distinction in z I R can be compensated by introducing extra boundary terms S ± similar to the ones for the vector meson case: These boundary terms will shift back the value of the infrared boundary z ± I R to the same value z I R for all a a 1 mesons. In the case, when φ 2 = 0, Eq. (60) can be solved in the UV (z → 0) and in the IR (z → z I R ) limits. In order to find the mass spectrum in this case it is reasonable to apply the IR boundary condition to the asymptotic solution found in the IR limit. For the IR asymptotic solution we shall take the z → z I R limit from (60) and set z = z I R in the condensate term. Before performing this approximation let us compare numerically the last two terms in Eq. (60) when z → z I R . The approximate values are: 3 , g 2 = 4π 2 /N c ≈ 13.2. Then g 2 σ 2 (z I R ) 4 ≈ 0.06 (GeV) 2 and 1/(z I R ) 2 = 0.1 (GeV) 2 . Thus, the 1/z 2 term contributes twice more than the g 2 σ 2 z 4 term and so we may make an approximation in (60) by setting z = z I R only in the condensate term and keeping the term 1/z 2 variable. In this limit the condensate term in the equations (60) becomes constant and the IR asymptotic solution of these equations is expressed in terms of the Bessel function J 1 : Obviously, the UV boundary condition was applied to this solution. The mass spectrumm s a in (71) is expressed in terms of z I R : As is seen from (72) condensate term contributes to the mass of states with different s differently. The mass splitting formula for the spectra (72) will be written as usual: and it has the splitting sign as in the case of absence of the condensate. In the case of taking into account the condensate term the mass spectrum formula (41) is valid for the mass m s a shifted by the condensate: This time a shift in mass is not simply ∓V 3 0 as in the previous case and it has a contribution proportional to the condensate value φ(z I R ). This contribution depends on s, and therefore, z ± I R this time depends on s as well. This complicates the introduction of the IR boundary terms compensating the z ± I R distinguished with the z 0 I R , which would be unique for all s states. For this case we have to introduce the boundary term for each state s, which depends on the massω s ± of this state as well, and the total boundary term will be a sum of boundary terms of the s modes, We should compensate as well the shift in the mass spectrum of the neutral a 0 1 meson, which was caused by the condensate term. As is seen from Eq. (72) the condensate shift in them s 0 spectrum is different for the different s states. This means one should distinguish values of z I R , these being s dependent. In order to avoid the distinguished z I R we introduce an additional boundary term for the each Kaluza-Klein s mode and sum these terms: One should notice that in [52], where the axial-vector meson's mass spectrum including a condensate contribution was obtained in the vacuum case, an approximation was performed which shifts the boundary condition at z I R and, thus, the condensate contribution to the mass of meson was taken into account. To do a shift in such a way is equivalent to introducing the boundary term S 0 in our case. 2 4.2 π Mesons in the isospin medium Following [47], we introduce the charged π a and χ a fields by means of the π i and χ i fields, respectively, as below: In the AdS/CFT correspondence the UV boundary values of the π a fields (a = 0, ±) will be mapped to the pion fields of the dual QCD. The Fourier components for the π a and χ a fields are written as where ω a are the masses of the π a (and χ a ) fields in the absence isospin background V 3 0 , i.e. they are the vacuum masses and do not include the contribution of the condensate. As in the ρ a and a a 1 meson cases the vacuum masses ω a of the π a fields do not depend on μ P and μ N and their equality ω ± = ω 0 is right in the isospin medium case as well. The action (57) in terms of the π a and χ a fields splits into three independent parts: where The equations of motion for these scalar fields are obtained in [47] from the action (79) and are summarized as follows: and The boundary terms arising on obtaining these equations are and the boundary conditions imposed on the π a fields are the usual ones: ∂ z π a (z) | z I R = 0 and π a (z) | z U V = 0. Comparing the first equations with the second ones in (81) and in (82) one can establish the following relations between the ∂ z derivatives of the π and χ fields: Taking the derivative from Eqs. (81) and (82) and taking into account the relations (84), the equation of motion for the π a field will be separated from the one for the χ a field: The physical situation with the π fields is more complicated in comparison with the vector and axial-vector fields considered above. As is seen from the action (57) the π fields interact with three other fields: the background condensate field φ, the longitudinal axial-vector field χ , and the background isospin field V 3 0 . In the action (57) the interaction with the φ field gives a contribution to the mass term of the π field and this contribution is the same for all π a components. The interaction with the isospin field V 3 0 gives different contributions to the mass terms of the different π a components in the equations of motion and, thus, splits the mass of these fields, which can be inferred from Eq. (85). The interaction with the χ field lifts up the order in the equations of motion, which can be seen on taking the χ field to be zero in the action (57) or in Eqs. (81) and (82). In the case χ = 0 we have second order differential equations, instead of the third order one in Eq. (85). Another observation is that the φ(z) field modifies the AdS metric (11). It is easy to see this in the χ = 0 case of the equations for the π fields in (81) and (82). Choosing the following notation in the metric tensor: we get the modified background geometry: Then and the equations for the π fields in (81) and (82) get the form As is seen from (89) ignorance of the χ field in the action leads to the absence of the condensate term 4g 2 φ 2 in the equation of motion. This case is useful for a comparison with the mass splitting formula for pions obtained in [34,47,48]. Equation (89) is, in fact, the five-dimensional Klein-Gordon equation for the free π a field with zero five-dimensional mass term in the background geometry (87): Consequently, the equation of motion for the π a field interacting with the φ field in the background geometry (11) is equivalent to the one for the free π a field in the modified geometry (87). It is obvious that this equivalence does not depend on the isospin medium and is valid for the V 3 0 = 0 case as well. It should be noticed that the background geometry modification by the scalar field φ was observed in [56,57] also. Using (90) we can derive a formula determining the effective mass of the π fields. In order to make the correspondence similar to previous sections let us write the Klein-Gordon equation in a four-dimensional curved space-time with the new metric G for the massive scalar field π a interacting with the constant external gauge field V 3 0 : where m a is the mass of the π a field. If we are to consider Eq. (91) as the UV limit of Eq. (90), then the eigenvalues of the operator − 1 √ −G ∂ z √ −GG zz ∂ z in (90) will correspond to the squared mass (m a ) 2 of the 4D vector field π a in (91). So, we may admit an equality of the eigenvalues of this operator to the effective mass of the π field defined on this boundary: Equation (92) is the one determining the effective mass of the π field in the isospin medium by means of the eigenvalue of the derivative operator over the extra dimension. For the Kaluza-Klein mode π a s in the decomposition a solution to Eq. (83) is expressed in terms of the Bessel functions J 1 and Y 1 : Here the eigenvaluesm s a are the mass spectrum: and it has a splitting sign opposite to the one in [34,47,48]. For the solution (93) the spectrum (41) is available due to the Neumann boundary condition at the IR boundary and, obviously, in this case it can be introduced usual additional boundary terms S ± in order to avoid the z ± I R distinguished from the z 0 I R introduced in previous sections: It is interesting to consider the χ = 0 case as well. In this case Eq. (85) contains the condensate term 4g 2 φ 2 and can be solved in the UV and IR limits of z. As is well known, in the free case the equation of motion for the π a field was solved either in the chiral limit (m q = 0) [58,59] or in the zero condensate limit (φ = 0). In our isospin medium case Eq. (85) will be solved in the UV and IR limits of z, i.e. in the z → 0 and z → z I R limits. In order to find the UV limit solution we take the z → 0 limit from Eq. (85), which leads to the usual Bessel equation for the ∂ z π a (z): It is obvious all threeM s a get the same value in the V 3 0 → 0 limit.
Using (98) and the property Z −n = (−1) n Z n of the Bessel functions we can write the solutions π a s : The UV boundary condition on this solution leads to c 2 = 0. In the UV limit, when z → 0, the σ 2 z 4 term drops out and the solution to (85) will be the same as the one for Eq. (66). The general solution to Eq. (85) will be a linear combination of these asymptotic solutions. It can be shown that the constants c 2 and c 2 are zero in this case as well and the general solution has got the form π a s (z) = c 1 z J 1 (m s a z) + We shall apply the Neumann boundary condition at the IR boundary only on the IR asymptotic part solution (106), i.e. only on the second term in (106). This gives us a discrete mass spectrum in terms of the zeros α s 2 of the Bessel function J 2 : In order to verify our consideration ofM s a as an effective mass we should reduce Eqs. (85) to the form in Eq. (89) in the UV limit z → 0. Let us write Eq. (85) in terms of the induced metric (86): + (D a 0 D a 0 + 4g 2 φ 5/2 G 00 )π a = 4g 2 π a ∂ z (φ 5/2 G 00 ).
The right hand side of this equation in the UV limit gives zero due to the UV boundary condition π a | 0 = 0. Then the free constant obtained in the left hand side can be chosen zero and (107) will be written as +G 00 (D a 0 D a 0 + 4g 2 φ 5/2 G 00 )π a = 0.
This equation has the form of the equations in (89) and it can be considered as a Klein-Gordon equation in the background geometry with the metric G M N . Then the second term will be considered as an effective mass including the condensate contribution. This proves that the derivative operator over the extra dimension z determines the effective mass through the eigenvaluesM s a in this non-zero condensate term case as well.
In this case the additional boundary terms compensating for the distinguished z I R will be dependent on the condensate g 2 σ 2 z 4 I R and the mass of the state ω s a : Since the condensate term 4g 2 φ 2 shifts the masses of all pions, the boundary terms above include it also.

Discussion
Here we obtain a formula for the meson's effective mass in the isospin medium using holography. We observe that interaction with the isospin background leads to the necessity of distinguishing of the infrared boundaries for different mesons from the same isotriplet. To bring all infrared boundaries to the same value requires one to introduce additional boundary terms, which leads to a modification of the boundary condition in the presence of an isospin background field. Indeed, the modification of the boundary condition at the infrared boundary is known from Refs. [37,52]. In [52], the axial-vector field interacting with the condensate field gets a contribution in mass and this results in a modification in the boundary condition at z I R . The shift in the boundary condition which was made in this work is equivalent to adding a new boundary term proportional to the condensate field. Here we have an isospin field interaction with the ρ mesons, which contributes to the mass differently for the ρ + and ρ − mesons. This results in the introduction of different boundary terms for these mesons. In the case of a 1 and π mesons we have two fields (condensate field and isospin field), the interaction with which contributes to the mass of these mesons. In the result we have a double shift in the mass and this requires one to introduce boundary terms proportional to both these fields.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecomm ons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP 3 .