Entanglement entropy and complexity in the holographic model of superfluid

We numerically study the holographic entangle-ment entropy and complexity conjectured with the volume in the holographic superfluid with full backreaction which can realize first and second order phase transitions. Our results show that both the entanglement entropy and complexity exhibit the behaviors characterizing the type of the transition. For the first order phase transition, there is a fast drop of the entanglement entropy and a fast jump of the complexity at the critical temperature, while both of them are continuous but non-differentiable at the second order phase transition point. These suggest that both of holographic entanglement entropy and complexity may be used as a good probe to the type of superfluid phase transition. Moreover, at a fixed temperature in the superfluid phase, we observe that the increasing superfluid velocity increases the entanglement entropy but decreases the complexity. Interestingly, we find that, for the condensation operators O+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {O}_{+}$$\end{document} and O-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {O}_{-}$$\end{document}, the dependence of holographic entanglement entropy on the backreaction is inconsistent and so is the dependence of holographic complexity on the superfluid velocity in the normal phase, which indicates that the entanglement entropy and complexity may reflect deep physics about the difference between the two operators in the superfluid dual system.


Introduction
The anti-de Sitter/conformal field theories (AdS/CFT) correspondence [1][2][3][4] is a very convenient and efficient tool to model strongly correlated systems in terms of a gravity dual.Since its discovery, this duality has made great development and been extensively used in various areas in physics.The a e-mail: laichuyu@gzhu.edu.cn(corresponding author) b e-mail: panqiyuan@126.comquantum information has important application in the study of quantum gravity and quantum field theory.Two quantities of the boundary field theory, which play important roles in quantum information, are entanglement entropy and complexity.Both of them have gravitational descriptions by certain geometric quantities in the bulk spacetime in holographic framework.
The entanglement entropy is a most studied idea in the last decades.It is directly related to the degrees of freedom in the quantum entangled state [5][6][7][8].In the light of the AdS/CFT correspondence, the holographic entanglement entropy describes the entanglement properties of the boundary field theory and might be eventually used to understand the quantum structure of spacetime.Ryu and Takayanagi [9,10] had provided an effective way to calculate the holographic entanglement entropy.The Ryu-Takayanagi proposal states that the entanglement entropy of the CFT in a closed region A with its complement is proportional to the area of a minimal surface γ A extended into the bulk from the boundary of the region A, and a simple and elegant formula is given as where G N is the Newton constant of the general gravity.
The complexity is recently introduced as a new measure of information of spacetime that can be useful to understand the black hole interior where the entanglement entropy can not cover [11].In quantum field theory, the complexity is defined as the minimal number of elementary operations (gates) that are required to obtain a desired target state from a reference state.There are two quintessential holographic proposals of the complexity.The first one is called "complexity=volume" (CV) conjecture [12,13] which relates the complexity of the boundary state to the volume of a codimension-one maximal hypersurface anchored to the given boundary time.The second one is the "complexity=action" (CA) conjecture [14,15] which relates the complexity of the boundary state to the onshell bulk gravitational action within the Wheeler-DeWitt (WDW) patch.Inspired by the Ryu-Takayanagi proposal of the entanglement entropy, Alishahiha constructed the holographic dual to the complexity for a subsystem on the boundary within the framework of the CV proposal, which is referred to as the holographic subregion complexity [16].Specifically, the complexity for a subregion A equals the volume enclosed by the Ryu-Takayanagi surface γ A where L is the AdS radius.In this work we will focus on the time-independent holographic subregion complexity, i.e., the minimal surface γ A totally outside of the black hole horizon.Among many sectors where the AdS/CFT correspondence has been applied successfully, the holographic realization of high-temperature superconductors has gained plenty of interest for their potential applications to the condensed matter physics.The simplest initial holographic superconductor model is constructed by applying a scalar field and a Maxwell field coupled in an AdS black hole background [17,18].The main idea is the spontaneous U (1) symmetry breaking induced by the condensation of the scalar hair close to the horizon of the black hole in low temperature.Since then, kinds of more complete holographic superconductor models have been explored, and many of their interesting features which resemble that of the real superconductors have been revealed, for reviews, see Refs.[19,20] and references therein.The entanglement entropy was firstly introduced as a tool to describe the phase transitions in the swave holographic superconductor in Ref. [21].By now a lot of works have been carried out for investigating the properties of phase transitions from the entanglement entropy side in various holographic superconductor models and in different backgrounds [22][23][24][25][26][27][28][29][30][31][32][33][34][35].And the entanglement entropy turns out to be a good probe to the critical point and the order of the phase transition in the holographic superconductor systems.In recent years, as there is a deep connection between the holographic entanglement entropy and holographic complexity, the complexity has also been discussed by using Alishahiha's subregion CV conjecture in different types of holographic superconductors, such as s-wave model [36,37], p-wave model [38], d-wave model [39] and so on [40][41][42][43].These studies show that the holographic complexity may be used as another independent probe to the physics of the phase transition as the holographic entanglement entropy does, though it was found that the complexity behaves in different ways from the entanglement entropy.Specially for the holographic Stückelberg superconductor [41] which allows both first and second order phase transitions to occur, the complexity behaves differently crossing the second order and first order superconducting phase transition point, which indicates that the complexity may also be used to determine the type of transition.
In this paper, we will generalize the investigation on the entanglement entropy and subregion complexity to the holographic model of superfluid with full backreaction, which has not been discussed as far as we know.Supercurrent in superconducting materials is a well studied phenomenon in condensed matter systems.The supercurrent solution can be holographically constructed by performing a deformation of the superconducting black hole, i.e., turning on the spatial component of the U(1) gauge field in the bulk theory [44,45].And it was found that the superfluid phase transition can switch from second order to first order when the velocity of the superfluid component increases relative to the normal component.In Ref. [46], considering various scalar field mass in AdS 5 , it was observed that the Cave of Winds phase structure exists for some special mass in the superfluid model.Along this line, further efforts in the s-wave holographic superfluid have been made in the related research and discussion [47][48][49][50], and more general and complex models have also been explored, such as p-wave superfluid [51][52][53][54][55] and s + p-wave superfluid [56].Besides, considering the holographic superfluid models away from the probe limit, the gravity backreaction also brings some interesting features to the system.In the p-wave superfluid system, it was argued that the order of the phase transition depends on the backreaction [57].The authors of Ref. [58] further studied the phase diagram modified by the backreaction in the supercurrent solution in the AdS black hole background.Their results showed that the backreaction has different impacts on the condensation behaviors for different operators, i.e., with the increase of the strength of the backreaction, the critical values of the superfluid velocity to accommodate the first order phase transition decreases for the operator O + but increases for the operator O − .Recently, the study has been enlarged to the backreacted superfluid model with the competition between the p-wave solution and the p + i p solution [59].Starting with the simple s-wave backreated holographic superfluid model, we would like to carry out a further analysis of this model from the aspects of the entanglement entropy and the complexity.It is meaningful to examine how the two quantities respond to the different types of phase transitions and characterize the phase structure under the influence of the backreaction, which may help us to enhance and supplement the description of the superfluid phase transition.We expect that the entanglement entropy and the complexity can suggest more deeply physical insights so as to better understand the underlying mechanism in the superfluid dual system.
The framework of this paper is as follows.In Sect.2, we briefly review the construction of a fully backreacted holographic superfluid in the four-dimensional AdS black hole spacetime and show the condensations and phase transitions.In Sect.3, behaviors of the holographic entanglement entropy and the holographic subregion complexity are investigated in detail to analyze the phase transitions in the superfluid model.The summary and discussions are included in Sect. 4.

Holographic model of superfluid with backreactions
We work with the simplest background geometry of a black hole and a U (1) gauge field coupling with a charged scalar field.The lagrangian density reads where R is the Ricci scalar, is the strength of the electromagnetic field.ψ denotes a scalar field with charge q and mass m, A represents the gauge field.The AdS radius L and the charge q are set to unity hereafter.And γ here is a parameter measuring the strength of the backreaction.
We are interested in studying the fully backreacted holographic superfluid model, thus we consider the ansatz for the line element of the four-dimensional AdS black hole with the form and the gauge field and the charged scalar field to be, respectively, The function C(r ) in the metric ansatz is introduced by considering a nonzero x component of the Maxwell field in order to take the possibility of DC supercurrent into account.The function B(r ) gives the position of the event horizon of the black hole, i.e.B(r h ) = 0.The Hawking temperature of the background is expressed as Then, the equations of motion derived from the action yield There is a useful scaling symmetry of the above equations where λ is a real positive number.For this superfluid phase transition system, there are two classes of the solutions for Eq. ( 7), i.e., ψ = 0 for the normal phase and ψ = 0 for the superfluid phase.In order to employ the shooting method to solve the field equations numerically, the adequate boundary conditions on the black hole horizon r = r h and the asymptotic AdS boundary r → ∞ are imposed.Around the horizon we can make the Taylor's expansion Note that a regular event horizon at r = r h needs B(r h ) = 0.And the requirement of regularity of the matter fields on the horizon leads to φ(r h ) = 0. Near the asymptotic AdS boundary (r → ∞), the various fields asymptotically behave like where ± = (3 ± √ 9 + 4m 2 )/2 is the conformal dimension of the operator O ± .According to AdS/CFT correspondence, μ and ρ represent the chemical potential and charge density in the dual field theory respectively.ξ and σ respectively describe the current density and the dual current source.The coefficients ψ − and ψ + can be regarded as the vacuum expectation value of the operator O in the dual field theory while the other is the source.In our analysis, we set the mass-squared of the scalar field m 2 = −2, then both ψ − and ψ + are normalizable.Thus we have two choices of boundary conditions, i.e., ψ − = 0, ψ + = O + and ψ + = 0, ψ − = O − , in which ψ + = 0 (ψ − = 0) corresponds to a vanishing source so as to describe a spontaneous symmetry breaking.For the rest of our discussion we would like to consider the both cases in order to explore the properties of the holographic entanglement entropy and complexity for the two different condensation operators O + and O − .
Using the scaling symmetry (8), we can rescale the relevant quantities as which allow us to employ the dimensionless quantities the computation of the condensates.Moreover, for the mathematical simplicity, we make a coordinate transformation from r -coordinate to zcoordinate by defining z = r h /r .The region r h < r < ∞ now corresponds to 1 > z > 0.
The numerical results of the condensates of the operators O + and O − with the changes of the temperature and the superfluid velocity given by the ratio σ/μ are plotted in Fig. 1.The scalar field begins to condensate as the temperature reduces to a critical value T c , and this implies the system undergoes phase transition to a superfluid phase.In particular, from the pictures, it can be seen that there is a threshold value (presented in orange) of the superfluid velocity σ/μ, below which the condensate O ± grows continuously and monotonically from zero near the critical point, signaling a typically second order phase transition.While we enhance σ/μ across the threshold, the condensate becomes multivalued, signaling a first order transition.That is to say, the type of superfluid phase transition is affected by the superfluid velocity, and when the superfluid velocity σ/μ increases beyond a translating value, the second order phase transition will switch to the first order one.
For the second order phase transition, the critical temperature T c can be directly read off from Fig. 1, at which the condensation appears.However, for the first order phase transition, the condensate is multi-valued near the transition point.In order to determine the critical temperature, we calculate the free energy of the system which is related to the product of the on-shell Euclidean bulk action and the tem-perature T as = −T S E [60,61].And in our model, where h is the determinant of the induced metric on the AdS boundary, and K is the trace of the extrinsic curvature K μν .
Using the equations of motion (7) and the boundary behaviors (10), we obtain the expression for the free energy of the boundary theory as with V 2 = dxdy which will be set to be unit and B 3 is a coefficient in (10).
As an example, fixing the backreaction parameter γ = 0.2, we plot the results of the free energy for operator O + with σ/μ = 0.15, 0.25 and for operator O − with σ/μ = 0.35, 0.45 in Fig. 2.Here is scaled by the chemical potential μ to be dimensionless.The red lines are for the superfluid phases while the blue lines represent the normal phases.The left two panels demonstrate the free energy for the two branches when the system undergoes the second order phase transition.We observe that below the critical temperature T c which is marked by the dashed vertical lines, the superfluid phase has lower free energy than the normal phase and becomes physically favorable.For the right two panels, the free energy develops a swallow tail, which is a typical characteristic for first order phase transitions.It is obvious that the temperature T = 0.0382021 and T = 0.165526 respectively representing the appearance of the condensates O + and O − , are not the physical critical point.The first order transitions happen at the interaction points T c = 0.0385111 for O + and T c = 0.166974 for O − , and below which the superfluid phase is thermodynamically favorable.Thus, from the free energy, we can obtain T c for the first order phase transitions.The results of T c for the various superfluid velocity we chosen in Fig. 1 in the case of γ = 0.2 are listed in Table 1.Obviously, the critical temperature T c decreases as the superfluid velocity σ/μ increases, which indicates that larger superfluid velocity makes the scalar hair more difficult to be developed, and is correct for other values of the backreaction parameter.
Furthermore, we also concern about the dependence of the condensates on the backreaction of the matter fields measured by the parameter γ .In Fig. 1, three values of the backreaction parameter have been considered, i.e., γ = 0.1 (left two panels), γ = 0.2 (middle two panels), and γ = 0.3 (right two  and O − , the backreaction has distinct effects on the translating value of σ/μ to accommodate the first order phase transition.Concretely, when the strength of backreaction increases, the translating superfluid velocity σ/μ decreases for O + , which shows that the sufficiently strong backreaction can trigger the first order transition of the system even for a small superfluid velocity.However, for O − , the translating superfluid velocity σ/μ increases with the increasing strength of backreaction, which is just opposite to the case of O + .Our findings is consistent with the results reported in Ref. [58].It is of interest to study the holographic entanglement entropy and holographic complexity in this backreacted model of superfluid phase transition.We expect the two quantities would reflect different properties in the two types of phase transition, and provide more detailed physics about the two condensation operators O + and O − .

Holographic entanglement entropy and subregion complexity
In this section, we will investigate the holographic entanglement entropy and holographic complexity in a strip geometry of the superfluid dual model.We consider a strip subregion A which is defined by where l is the size of region A and R is a regulator which can be set to be infinity.The induced metric of the hypersurface γ A whose boundary is the same as the strip reads Then, using the proposal given by Eq. ( 1), we have the holographic entanglement entropy connecting with the area of γ A as We now apply the Hamiltonian approach to impose the following minimality condition on the surface area where z * is defined as the turning point of the smooth extremal surface such that dz dx | z * = 0, then C * denotes C(z * ).Integrating the condition (17), we obtain which satisfies x(z * ) = 0 and x( → 0) = l/2 with a UV cutoff .Subsequently, we can deduce the resulting holographic entanglement entropy as with R = dy.The 1 term is the UV divergent term and known as the area law, which will not change since the new solution after the operator condensation still asymptotically approaches to AdS space near the AdS boundary.Subtracting this divergence from S A , we have the finite term s which is physically important.On the other hand, according to Ref. [16], the complexity for the subregion A is holographically related to the volume surrounded by the aforementioned minimal surface (Ryu-Takayanagi surface) γ A .So it can be evaluated as Note that the complexity C A also includes a universal term c and a divergent term which is associated with a function of z * as F(z * )/ 2 .Though it is hard to find a general form of the divergence analytically, the value of F(z * ) in various situations can be found numerically [37] so that we can confirm the singular part in C A and pick up the finite term c.Due to the presence of the scaling symmetry (8), the strip width l, the entanglement entropy s, and the subregion complexity c can be transformed into Therefore, it is useful to work with the following dimensionless quantities to study the physics Fixing the width l μ 2 − σ 2 = 2 since the other choices will not qualitatively modify our results, we present the holographic entanglement entropy s and subregion complexity c as a function of temperature T for different superfluid velocity σ/μ with the backreaction γ = 0.1, 0.2 and 0.3 in Fig. 3 (for operator O + ) and Fig. 4 (for operator O − ), where the dashed curves represent the normal state and the solid curves represent the superfluid state.We can clearly observe that the entanglement entropy and complexity exhibit behaviors characterizing two types of phase transition.For the second order phase transition with σ/μ below the translating value to accommodate the first order phase transition, both the entropy s and complexity c are continuous but have a discontinuous first derivative at the critical point where the superfluid phase (solid curves) intersect with the the normal phase (dashed curves).And under the same conditions, the discontinued slopes of the both quantities occurs at the same critical temperature, at which the condensation appears as shown in Fig. 1.The entropy s monotonously decreases as the temperature decreases, and its value in the superfluid phase is always smaller than in the normal phase.This is comprehensible because the formation of the Cooper pairs in the superfluid phase suppresses the number of the effective degrees of freedom, and the entanglement entropy is just a measure of the degrees of freedom in the field theory.Contrary to s, the complexity c monotonously increases with reducing T , so the superfluid phase has a larger c than the normal phase, which can be interpreted as the quantum state of the system becoming more complicated in the superfluid phase since the complexity describes the difficulty of turning a certain reference state into a target state.
For the first order phase transition with σ/μ above the translating point, it is interesting to note that both two quantities have a multi-valued region around the critical temperature T c , and the similar multi-valued behaviors appear in the same temperature region in the free energy diagrams and the condensation diagrams.According to the free energy of the system, below the critical temperature, the lower entanglement entropy and larger complexity are physically important in the multi-valued region.Thus, the real process of entropy s has a discontinuous drop from the dashed curve to the lowest part of the solid curve at the critical point and where the complexity c has a discontinuous jump from the dashed curve to the highest part of the solid curve.Similar phenomena was also found in the first order transition for the holographic Stückelberg superconductor [41].This is an important distinguishing feature of the first order transition from the second order one in our model.Except the suddenly drop or jump, the entropy s is a monotonic decreasing function but the complexity c is a monotonic increasing function with respect to T , which agree with the case of the second order transition.
In addition, we find that the critical temperature T c decreases with the increasing value of the superfluid velocity σ/μ, that is to say, the larger superfluid velocity makes the scalar hair harder to be developed.This observation is consistent with the previous findings in the condensates in Fig. 1 as well.
Therefore, both of the holographic entanglement entropy and complexity can capture the superfluid phase transition in the backreated dual model and also be a good probe to the type of the phase transition.We argue that the performances of the entropy and complexity crossing the two types of phase transition can be perceived as qualitative characteristic in general holographic model.
To see clearer of the effect of the superfluid velocity σ/μ on the entanglement entropy and complexity, the two quantities as a function of σ/μ in the superfluid phase with chosen values of the backreaction parameter γ are respectively displayed in Fig. 5 for the operator O + at T √ μ 2 −σ 2 = 0.028 and in Fig. 6 for O − at T √ μ 2 −σ 2 = 0.145.We observe that, for a fixed backreaction γ , the larger superfluid velocity results in a larger entanglement entropy s, but makes for a smaller complexity c.
All the above features are shared by the case of operators O + and O − .Considering the disagreement between the two condensation operators about the effect of the backreaction on the translating superfluid velocity σ/μ shown in Fig. 1, we expect that the holographic entanglement entropy and complexity for different operators would perform different properties.Comparing the results in Figs. 3 and 4, we are surprised to find that in the normal phase, the larger superfluid velocity results in a smaller complexity c in low temperature region before the condensate O + but a larger c in relatively high temperature region before the condensate O − , which can not be seen in the entropy s since it monotonously increases as σ/μ increase for both the cases.However, from Figs. 5 and  6, we obtain that with the fixed temperature and superfluid velocity, the increasing value of the backreaction parameter γ increases the entropy s for operator O + but decreases s for O − , while the complexity c for both the operators always rises with the strengthening backreaction.These observations indicate that the holographic entanglement entropy and complexity contain different information about the superfluid system, meanwhile both the two quantities may capture the variance between the condensation operators O + and O − and have the potential to help us better understand the physical mechanism behind the two condensates in the dual model.

Conclusions
In this paper, we have numerically investigated the holographic entanglement entropy and complexity associated to a strip-shaped subregion for the superfluid phase transition in the four-dimensional AdS black hole with full backreaction.This holographic superfluid model can realize first and second order phase transitions which depends on the superfluid velocity.Moreover, for the different condensates O + and O − , the strength of the backreaction has different effects on the translating value of superfluid velocity to accommodate the first order transition of the system.
Our results show that regardless of either the operators O + or O − , both the holographic entanglement entropy and complexity present particular features distinguishing the two types of phase transitions, though the two quantities behave differently.For the second order phase transition with σ/μ below the translating value, both the entropy s and complexity c are continuous but not differentiable at the critical temperature T c .With the decreasing temperature, the entropy s monotonously decreases while the complexity c increases so that the superfluid phase always has a smaller s but larger c than the normal phase.Some physical interpretations can be given to the behaviors of the both quantities.As temperature goes down, the Cooper pairs forming in the condensed phase makes the reduction of the number of degrees of freedom which causes the drop of the entanglement entropy.And a larger complexity at lower temperatures is related to the quantum state of the system becoming more complicated across the phase transition.For the first order phase transition with σ/μ above the translating value, both s and c have a multi-valued area around the phase transition point, and the real processes of them are discontinuous at the critical temperature T c , i.e., from the normal state to the superfluid state, the entropy s experiences a suddenly drop and the complexity experiences a suddenly jump.The discontinuous change of the entropy may be related to some kind of non-trivial reorganization of the degrees of freedom of the system, which probably includes extra degrees of freedom coming from other entanglement apart from the entangled pairs.And the introduction of this extra degrees of freedom may be due to some In addition, we find that under the same conditions, the holographic entanglement entropy and complexity give the same critical temperature of the second order phase transition, at which the condensation appears as shown in Fig. 1.For the first order phase transitions, near the critical temperature the superfluid phase of both the holographic entanglement entropy and complexity have a multi-valued region, and the similar multi-valued behaviors appear in the same temperature region in the free energy diagrams and the con-densation diagrams.And we can see that the larger superfluid velocity makes the critical temperature lower, which agrees well with the observation in the condensations.It indicates that the spatial component of the gauge field will hinder the phase transition.Therefore, the transition temperature of the system is quantitatively in agreement with the insights from holographic entanglement entropy and holographic complexity.
We also present the entanglement entropy and complexity in the superfluid phase as a function of the superfluid velocity σ/μ with different values of the backreaction parameter, i.e., γ = 0.1, 0.2 and 0.3 at It is shown that when the value of the superfluid velocity σ/μ increases, the entanglement entropy always rises and the complexity drops, which is independent of the backreaction.
On the other hand, by comparing our results of the holographic entanglement entropy and complexity for the operators O + and O − , it is interesting to find that there exist some inconsistent performances of both the quantities for the two operators.One is that in the normal phase, the increasing superfluid velocity decreases the complexity in low temperature region before the condensate of operator O + , but increases c in relatively high temperature region before the condensate of operator O − .The other is that the increasing value of the backreaction parameter increases the entropy s for the operator O + but decreases s for O − .It implies that both the holographic entanglement entropy and complexity have the potential to offer richer physics about the difference between the operators O + and O − in the dual system, and the information reflected by the complexity is different from the information captured by the entanglement entropy.The physical origins behind them and non-trivial insights from the field theory side deserve further study.
The Cave of Winds is an interesting phenomenon appearing in some holographic superfluid system, such as the fivedimensional gravity dual model and the model in Lifshitz spacetime.Thus, we hope to examine how the holographic entanglement entropy and complexity describe the Cave of Winds phase structure in the near future.And, it would also be interesting to extend our numerical analysis to the holographic superfluid model with the competition and coexistence of different orders, so as to further understand the deep physics involved in the rich phase diagram in holographic superfluid.
It is noticed that here we computed the holographic complexity via Alishahiha's subregion CV conjecture.As is known that there are also other interesting holographic proposals of the complexity, which have been discussed in various holographic models.In Ref. [62], the authors investigated the complexity of formation, which defined in Ref. [63] within the framework of the CV proposal, in the first holographic superconductor model proposed by Hartnoll et al. [17,18], and analyzed the full time evolution of the complexity from the normal phase to superconducting phase.It was found that the complexity of formation can also be a good probe to the superconducting phase transition.In Ref. [64], by considering the CA proposal, the holographic complexity was studied in a holographic QCD model proposed by Gubser et al. [65], which can realize three types of phase transition, crossover or first and second order.Their results indicated that the growth rate of the holographic complexity exhibits the behavior characterizing the type of phase transition.Investigating the complexity by other holographic proposals in the holographic superfluid will be another interesting direction.We leave this issue as future works.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material.If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/.Funded by SCOAP 3 .SCOAP 3 supports the goals of the International Year of Basic Sciences for Sustainable Development.

Fig. 1 3 Fig. 2
Fig. 1 The condensates O + and O − in terms of temperature in the backreacted holographic superfluid model.The left two panels are for the backreaction γ = 0.1, the middle ones are for γ = 0.2, and the right ones are for γ = 0.3

Fig. 3 Fig. 4
Fig.3The holographic entanglement entropy and subregion complexity of the operator O + versus temperature for different values of the superfluid velocity σ μ with the strength of the backreactions γ = 0.1 (left two panels), 0.2 (middle two panels) and 0.3 (right two panels)

Fig. 5 Fig. 6
Fig.5 The holographic entanglement entropy and subregion complexity of the operator O + as a function of the superfluid velocity σ μ with the strength of the backreactions γ = 0.1 (left two panels), 0.2 (middle two panels) and 0.3 (right two panels), as T √ μ 2 −σ 2 = 0.028 in the superfluid phase

Table 1
The critical temperature T c for operators O + and O − with a number of superfluid velocity σ/μ and the fixed backreaction parameter γ = 0.2 It is interesting to find that, for the condensates O +