Probing phase structure of strongly coupled matter with holographic entanglement measures

We study the holographic entanglement measures such as the holographic mutual information, HMI, and the holographic entanglement of purification, EoP, in a holographic QCD model at finite temperature and zero chemical potential. This model can realize various types of phase transitions including crossover, first order and second order phase transitions. We use the HMI and EoP to probe the phase structure of this model and we find that at the critical temperature they can characterize the phase structure of the model. Moreover we obtain the critical exponent using the HMI and EoP.

The gauge/gravity duality, which relates a d-dimensional quantum field theory (QFT) to some classical gravitational theory in (d + 1)-dimensions, has been utilized to study various areas of physics such as condensed matter, quantum chromodynamics (QCD) and quantum information theory [1][2][3][4][5][6].This duality provides a powerful tool to investigate various phenomena in the strongly coupled QFT using their corresponding gravity duals.Studying some phenomena and quantities from a field theory point of view may be difficult while this is relatively simple to study them on the gravity side.For example the confinement-deconfinement phase transition of the QCD corresponds to the Hawking-Page phase transition on the gravity side which is the transition through black hole and thermal gas backgrounds [7,8].
In this paper we consider a holographic QCD model which is proposed in [9,10] intending to mimic the equation of state of the real QCD.As we know, QCD is strongly coupled at low energy and hence it is difficult to study it using perturbation techniques.Surprisingly, the gauge/gravity duality is extended to the more general cases such as the non-conformal field theory and hence one can study the QCD in the strongly coupled region by studying its dual gravitational theory [11][12][13][14][15][16][17][18].It is interesting to seek the holographic models which are dual to non-conformal field theories, such as QCD.One class of these models are top-down models which are directly constructed from string theory.One can find some examples of these models in [19][20][21][22][23][24][25].Another class of these models are bottom-up models such as the hard-wall model and the soft-wall model in which the gravitational theory are phenomenologically fixed to be in agreement with the lattice QCD data [26][27][28].Our considered holographic QCD model contains a dilaton field whose self-interaction potential is parameterized by some constants, and by choosing appropriate values for them the system possesses a crossover phase transition at a critical temperature and it is clearly seen that the thermodynamical properties obtained from this model are in complete agreement with the lattice QCD results [96].Moreover, by choosing other suitable values for the parameters in the dilaton potential, this model can also exhibit the first and second order phase transitions.Therefore this is a good model to study various phase structures of strongly coupled field systems.This model utilized for studying the phase structure of a system using the hydrodynamic transport coefficients [97], quasinormal modes [98], entanglement entropy [99] and complexity [100].
There are various physical quantities in the context of quantum information theory which we can study with the help of their corresponding gravity duals.One of them is the entanglement entropy (EE), a significant non-local quantity, which measures the quantum entanglement between subsystem A and its complement Ā for a given pure state.The EE suffers from short-distance divergence which satisfies an area-law and hence it is a scheme-dependent quantity in the UV limit [29,30].It is difficult to compute the EE in QFT.Fortunately, EE has a simple holographic dual in which the EE of a subregion A in the boundary field theory corresponds to the area of minimal surface extended in the bulk whose boundary coincides with the boundary of subregion A [31,32].Based on this prescription, many studies have been done in the literature to better understand the EE [33][34][35][36][37][38][39][40][41][42][43][44][45][46][47].However, for a mixed state, the EE is not a good measure of correlation and to do so one can study other quantities such as mutual information (MI) which measures the total (classical as well as quantum) correlation between two subsystems A and B which is defined as I(A, B) = S(A) + S(B) + S(A ∪ B) where S denotes the entanglement entropy of its associated subsystem [48][49][50][51][52][53][54][55][56][57].When the two subsystems A and B are complementary to each other (total system describes by a pure state), I(A, B) = 2S(A) = 2S(B).Since the divergent pieces in the EE cancel out, the MI is a finite scheme-independent quantity and it is always non-negative because of the subadditivity of the EE.
Another important quantity which measure the total correlation between two subsystems A and B in a mixed state is entanglement of purification (EoP) [58,59].In general, we can purify a mixed state ρ AB to a pure state |ψ AA BB by enlarging the Hilbert space.The EoP is defined by the minimum of the EE between two subsystems AA and BB among all possible purifications.From this definition it is clear that the EoP between two subsystems A and B reduces to the EE of subsystem A if one considers the subsystem B as complementary of subsystem A. There is a holographic prescription to calculate the EoP in which this quantity is dual to the minimal cross-section of entanglement wedge Ew of ρ AB [60,61].When ρ AB describes a pure state, Ew between subsystems A and B reduces to the holographic entanglement entropy (HEE) of subsystem A. There are other correlation measures such as reflected entropy, odd entropy and logarithmic negativity which are discussed in the literature and holographically all of them related to Ew [62][63][64][65][66][67][68].Various aspects of the EoP are discussed in the literature .
The remainder of this paper is organized as follows.In section II, we review the background and describe the thermodynamics of the considered model.In section III, after a short review on HEE, we introduce the MI and EoP and their holographic duals.In section IV we explain our numerical results and describe how the HMI and EoP characterize various phase structures of the system.Finally, we will conclude in section V with the discussion of our results.

II. THE BACKGROUND AND THE THERMODYNAMICS OF THE SYSTEM
As we mentioned in the introduction, we are interested in studying the phase structure of the (3+1)-dimensional non-conformal field theory using the HMI and EoP.Therefore, we present a review on the five dimensional gravitational theory which is dual to the aforementioned field theory [9,10] and study the thermodynamics of this model.We consider black hole solutions which follow from the following 5-dimensional Einstein-dilaton action where G (5) N is the 5-dimensional Newton constant and we consider 16πG (5) N = 1 in the following discussion.g and R are the determinant of the metric and its corresponding Ricci scalar, respectively.φ denotes dilaton field and V (φ) is the dilaton potential which is considered in the following form [97,98] Each set of parameters {γ, b 2 , b 4 , b 6 } leads to a black hole solution.This family of black hole solutions has been introduced in [9,10] and used to mimic the equation of state of QCD.For small φ, one can expand V (φ) as follows where m 2 ≡ 2(b 2 −6γ 2 ).The first term is the negative cosmological constant (we fix the AdS radius to one).According to the AdS/CFT dictionary, the scalar field φ in the bulk is dual to a scalar operator in the dual boundary field theory O φ which is known as field-operator correspondence [2,6].The conformal dimension of the scalar operator, ∆, is related to the mass of the scalar field in the bulk through ∆(∆ − 4) = m 2 .Holographically, the presence of the dilaton field in the bulk is dual to a deformation of the boundary conformal field theory where Λ is an energy scale and ∆ is considered to be at range 2 ∆ < 4 corresponding to the relevant deformation of the CFT and satisfies the Breitenlohner-Freedman bound m 2 −4 [101,102].Different choices of the parameters {γ, b 2 , b 4 , b 6 } lead to different thermodynamical properties of this model.We consider three sets of parameters for dilaton potential, labeled by V QCD , V 2nd and V 1st which are given in Table I.The parameters for the V QCD are chosen to fit the lattice QCD (lQCD) data from [96] and from the QCD phase diagram we know that at zero chemical potential, the system exhibits a crossover phase transition at a critical temperature.The parameters for the V 1st and V 2nd are chosen so that the system exhibits the first and the second order phase transitions at a certain critical temperature, respectively.
The equations of motion following from the action (1) are By considering the following ansatz in the r = φ gauge [9, 10] the equations of motion (5) become simple and are given by where X = dX dφ .The strongly coupled field theory lives on the boundary at φ → 0 and the location of the black hole horizon is determined by the simple zero of the blackening function h(φ) In order to solve (7) we review the method proposed in [9,10] known as master equation formalism.If we consider G(φ) ≡ A (φ) the solution of the equations of motion can be expressed as where A H , B H , h H and h 1 are the integration constants which are obtained by applying appropriate boundary conditions at the black hole horizon (8) and the boundary of this background φ → 0 [9, 10] By manipulating the field equations ( 7) one can derive the following "nonlinear master equation" By solving this equation and obtaining G(φ), the metric coefficients A(φ), B(φ) and h(φ) can be obtained from (9).In order to solve (11) we need appropriate boundary conditions.From (11) one can find the following near-horizon expansion This expansion implies that Unfortunately, it is hard to solve (11) analytically and hence we use numerical methods to solve it.Our strategy is to specify a value for φ H and use the boundary conditions ( 12) and ( 13) for integrating (11).Each value of φ H leads to a unique black hole solution.We obtain numerically a family of black hole solutions for different values of φ H .

B. The thermodynamics
In this subsection we would like to study the thermodynamics of the field theory which is dual with the black hole background described by metric (6).The Hawking temperature can be obtained from the metric ( 6) and the Beckenstein-Hawking formula for entropy density leads to At zero chemical potential the speed of sound and the specific heat are related to the temperature and the entropy density as One can obtain the free energy from the first law of thermodynamics, dF = −sdT .At zero chemical potential the free energy of a system with constant volume are given by We take T 0 to be zero.Since the metric of the zero temperature background (thermal gas) coincides with the black hole metric in the limit φ H → ∞, we expect that the free energy of the black hole background vanishes in this limit, F (T 0 = 0) = F (φ H → ∞) = 0 and hence the free energy is given by We plot the temperature T v.s.horizon φ H in FIG. 1. From the left diagram, related to V QCD , we observe that the temperature monotonically decreases to zero.In this case the black hole solutions are always thermodynamically stable 1 and a smooth crossover phase transition happens there.In the middle diagram where we consider V 2nd , there is a critical point at (φ 2nd c 4.27, T c 0.1568) which dT dφ H becomes zero at this point and it is a point of inflection.In subsubsection II B 2 we will discuss that the system possesses a second order phase transition at this point.In the Right diagram, corresponding to V 1st , it is shown that the temperature has a local minimum T min ≡ T m 0.2321 and a local maximum min and φ 1st max the black hole solutions are thermodynamically unstable.We expect a Hawking-Page phase transition happens at a temperature between T min and T max .We will discuss it later in subsubsection II B 3.
In the following, we will discuss the details of thermodynamics properties of the dual field theory for dilaton potentials V QCD , V 2nd and V 1st , respectively.

VQCD
In FIG. 2 we show the dependence of the entropy density s, the square of speed of sound c 2 s , specific heat C v and free energy F in terms of temperature T for V QCD .The dependence of c 2 s on T is in complete agreement with the 1 If the thermodynamic potential of a system is given by ψ(x 1 , ..., xn) which depends on some set of variables {x 1 , ..., xn}, the Hessian matrix H of the associated potential defined by . The system is called stable when the Hessian matrix is positive-definite.lattice QCD results [96].There is a critical temperature T c corresponding to the lowest dip of the c 2 s and it is clearly seen that s and ds dt ∝ C v are both continuous at temperature T Tc = 1.The free energy of the black hole is always negative and is less than the free energy of the thermal gas.Therefore, the thermodynamic system will always favor the black hole background and the Hawking-Page phase transition does not occur which means that the system has a crossover phase transition.From FIG. 2a it is seen that in the limit T → ∞, the value of c 2 s approaches to its conformal value c 2 s = 1 3 which is expected.This result is also valid for the potentials V 2nd and V 1st .

V 2nd
In FIG. 3 the dependence of s, c 2 s , C v and F on T for V 2nd are depicted.c 2 s v.s.T shows that at a critical temperature T c , c 2 s goes to zero but never becomes negative.F and dF dT = s are finite and continuous at the critical temperature T c but the value of d 2 F dT 2 ∝ C v diverges at T c and hence the phase transition is second order.The entropy density drops quickly as the temperature approaches to T c .Near the critical temperature, s and the C v of the system take typically the form where t ≡ |T −Tc| Tc and α is the specific heat critical exponent.In order to find the critical exponent, we focus on the region near the critical point and plot C v (T ) in this region in FIG. 4. By fitting a curve with the numerical result, the value of the critical exponent is found to be α = 0.67.This value is in complete agreement with the one reported in [9,98].To get this number, we also plot the linear log-log diagram for which the critical exponent is the slope of a line i.e. log(C v ) ∝ α log(t).To report how well our result is, we calculate relative error (RE) and root mean square (RMS) which are defined as where y fitted (i) is the value of fitted function y evaluated at ith data points x, y data is the corresponding value read from data and N is the number of data points.These numbers are reported in the caption of FIG. 4.

V1st
In FIG. 5, the dependence of s, c 2 s , C v and F on T for V 1st are shown.From these diagrams one can see that in the range of temperature T m < T < 1.272T m there are three branches of solutions where two of them correspond to the stable black hole solution (solid blue curves) and one of them corresponds to the unstable black hole solution (dashed cyan curve).There are no unstable solutions at T < T m and 1.272T m < T .The dashed curves show that the system indicates the Gregory-Laflamme instability [103,104] which can be seen from the behavior of C v .In [105][106][107] it is shown that in the absence of conserved charges related to gauge symmetries, the Gregory-Laflamme instability is equivalent to the negative values of C v which is depicted clearly in panel 5c.From the panels 5a and 5c we observe that the negative values of C v corresponds to the negative c 2 s .F (T ) shows that the phase transition occurs at critical temperature T c 1.05T m in which the two stable branches of solutions cross each other at this point.Although F is continuous across the phase transition, dF dT = s is not continuous which means the system possesses a first order phase transition at T c .Since s/T 3 counts the effective number of DOF, this number increases suddenly when we approach to T c from the left [99].

III. THE HOLOGRAPHIC ENTANGLEMENT MEASURES
In this section we would like to study the HMI and EoP in the black hole background described by metric (6) and investigate these entanglement measures behavior near the phase transition temperature.These quantities are known in quantum information theory and measure total, both classical and quantum, correlation between subsystems A and B for the total system described by mixed state ρ AB .

A. The holographic mutual information
When the system describes by a pure state, the EE is a unique measure which determines the quantum entanglement between subsystem A and its complementary Ā.If we consider a pure quantum system described by the density matrix ρ = |ψ ψ| and divide the total system into two subsystems A and its complement Ā, the entanglement between these subsystems is measured by the EE which is defined as where ρ A = T r Ā(ρ) is the reduced density matrix for the subsystem A. It is shown that the EE contains short-distance divergence which satisfies an area law and hence the EE is a scheme-dependent quantity in the UV limit [29,30].It B A FIG. 6: Left: A simplified sketch of a strip region A with width l and length L. ΓA is the RT-surface of the region A and φ * is the turning point of this surface.Right: A simplified sketch of two strip regions A and B with equal size l which are separated by the distance l .When l is small enough, the minimal surface of A ∪ B are given by Γ l ∪ Γ l +2l and when l is large enough, the minimal surface of A ∪ B are given by 2Γ l .The minimal surfaces, the dashed curves, are denoted by Γ.
is difficult to calculate the EE using QFT techniques.However, in the framework of the gauge/gravity duality, there is a simple prescription to compute entanglement entropy in terms of a geometrical quantity in the bulk [31,32].According to this prescription, the holographic entanglement entropy (HEE) is given by where Γ A is a codimension-2 minimal hypersurface, called Ryu-Takayanagi surface (RT-surface), in the bulk whose boundary coincides with the boundary of region A. In [99] the behavior of the HEE is studied in three dilaton potentials, by choosing suitable parameters to the scalar self interaction potential, and they showed that the HEE can be characterize the crossover/ phase transition.
If there are two disjoint subsystems on the boundary entangling region, one of the most important quantity to study is the mutual information which measures the total correlation between the two subsystems, including both classical and quantum correlations, in a mixed state [109].The mutual information between two disjoint subsystems A and B is defined as a linear combination of entanglement entropy where S A , S B and S A∪B denote the entanglement entropy of the region A, B and A ∪ B, respectively.From ( 23) one can conclude that the mutual information is a finite quantity since the divergent pieces in the entanglement entropy cancel out and it is always positive because of the subadditivity of the entanglement entropy, S A + S B ≥ S A∪B .We consider the two symmetric disjoint subsystems both rectangular strips of size l which are separated by the distance l on the boundary, see FIG. 6.Using holographic prescription we can compute easily the HEE of the individual subsystems A and B. In order to compute S A∪B , we have two possible configurations.When the separation distance is large enough (disconnected configuration), the two subsystems A and B are completely disentangled and we have S A∪B = S A + S B = 2S(l), and hence the mutual information vanishes I(A, B) = 0. On the other hand, when the two subsystems A and B are close enough to each other (connected configuration), S A∪B = S(l ) + S(l + 2l) and we get I(A, B) > 0. One can assume that the transition of the mutual information from positive values to zero occurs at the distance which we call x d .To summarize, I(A, B) is given by B. The holographic entanglement of purification When the system is described by a mixed state, another important quantity to study is the entanglement of purification which measures the total (quantum and classical) correlation between two disjoint subsystems.In order to define the EoP, we consider a mixed bipartite system with density matrix ρ AB .we can always purify this mixed state into a pure state |ψ AA BB by adding auxiliary degrees of freedom to the Hilbert space as H A ⊗H B → H A ⊗H B ⊗ FIG.7: The gray region shows the entanglement wedge dual to ρAB and Σ min AB , the green curve, is the entanglement wedge cross section between subregions A and B.Here we only show the connected configuration.For disconnected configuration the entanglement wedge cross section vanishes Ew = 0.The RT-surfaces, the dashed curves, are denoted by Γ.
H A ⊗ H B such that the total density matrix in enlarged Hilbert space is given by ρ AA BB = |ψ AA BB ψ AA BB |.This pure state is called a purification of ρ AB if we have ρ AB = T r A B (|ψ AA BB ψ AA BB |).The EoP is then defined by minimizing the entanglement entropy S AA over all purifications of ρ AB [58] where S AA is the entanglement entropy corresponding to the density matrix In general, it is a difficult task to compute the EoP in the context of the QFT.Holographically, it has been conjectured that the EoP is dual to the entanglement wedge cross-section E w of ρ AB which is defined by [60,61] where Σ min AB is the minimal surface in the entanglement wedge E w (ρ AB ) that ends on the minimal surface Γ A∪B , the dashed blue line in FIG. 7. As a result, we have [60,61]

IV. NUMERICAL RESULT
In order to study the HMI and EoP in the three cases of the dilaton potentials, we use ( 24) and ( 26) and compute the area of the minimal surfaces and the entanglement wedge cross section for the metric (6).It is hard to do analytical calculation in the background ( 6) and hence we use numerical methods to obtain the HMI and EoP.We compute the HMI and EoP for two different values of l and l i.e. (l = 0.05, l = 0.01) and (l = 0.05, l = 0.005).In FIGs.8 and 9, we plot the HMI and EoP in terms of temperature for the dilaton potentials V QCD , V 2nd and V 1st .From these figures one can see that in the high temperature limit, the values of the HMI and EoP approaches to a constant values in the three cases of potentials which is equivalent to their conformal values.This behavior has been seen for the growth rate of holographic complexity density in [100].From the left panels in FIGs.8 and 9, V QCD case, we observe that the HMI and EoP have smooth behavior.In this case the HMI and EoP and their derivatives with respect to the temperature are continuous which show that there is a crossover phase transition at the critical temperature.The dependence of the HMI and EoP on temperature in the V 2nd case, middle panels, show that although at the critical temperature T c the HMI and EoP are both finite, their derivatives with respect to the temperature show a power law divergence in the vicinity of T c .This indicates that the system possesses a second order phase transition at T c .We will study the critical exponent using the HMI and EoP later on.In the right panels of the FIGs.8 and 9, we plot the HMI and EoP for the V 1st potential.Similar to the thermodynamic quantities in the range of temperature T m < T < 1.272T m which we mentioned in subsection II B, there are three branches of solutions which two of them correspond to the thermodynamically stable black hole solutions (solid blue curves) and one of them corresponds to the unstable black hole solution (dashed cyan curve).At the critical temperature T c 1.05T m the HMI and EoP and their derivative with respect to the temperature are not continuous at T c which means the system possesses a first order phase transition at T c .In short, we conclude that the behavior of the HMI and EoP suggests that one can use these quantities to characterize the type of phase transition.The interesting point is that although the shape of the HMI and EoP functions alter by changing l and l , their behavior does not change in the vicinity of T Tc = 1 and hence we can use them for probing the phase structure of the strongly coupled matter.This is in complete agreement with the result obtained for the second order phase transition in [83].

A. Critical exponent
When the system enjoys the second order phase transition at a critical temperature T c , we saw that the HMI and EoP are continuous at T c while their slope with respect to temperature are not continuous.How the behavior of these quantities are near the critical temperature is an important question and we will discuss it in the following.We focus on the vicinity of the critical temperature in the HMI and EoP diagrams for V 2nd and this suggests that the slope of these plots for the near of T c can be fitted with a function of the form t −θ where t ≡ |T −Tc| Tc and hence the number θ, called the critical exponent, describes the variation of the HMI and EoP with respect to temperature.In FIG. 10, we plot the slope of the HMI and EoP for l = 0.05 and l = 0.01 near the critical temperature where we have defined the slope of a quantity y with respect to temperature T as where i represents the ith point in the corresponding data points.By fitting a curve with the numerical result, the value of the critical exponent is obtained 0.668 and 0652 from the HMI (left panel) and the EoP (right panel), respectively.These values are in agreement with one obtained from the behavior of the specific heat in [9,98].We also plot the linear log-log diagrams for the HMI and EoP where the critical exponent is the slope of a line i.e. log( dy dT ) ∝ θ log(t), and y denotes the HMI or the EoP.We calculate the RE and RMS which are reported in the caption of FIG. 10.

V. CONCLUSION
In this paper, we consider a non-conformal field theory at finite temperature which has holographic dual.This holographic model is used to mimic the equation of state of QCD by introducing a nontrivial dilaton field whose corresponding potential break the conformal symmetry.The dilaton potential are given by four parameters and the values of these parameters can be chosen to reproduce the lattice QCD results in which the system exhibits a crossover phase transition.Moreover, different choices of these parameters lead to different thermodynamical properties of this model.We choose three set of parameters for dilaton potential, labeled by V QCD , V 2nd and V 1st in which the system exhibits respectively the crossover, first and the second order phase transitions at a certain critical temperature.In this model, we study some thermodynamical quantities i.e. entropy density s, speed of sound c 2 s , specific heat C v and free energy F whose behaviors confirm the phase structure of the system.We calculate the HMI and EoP for a symmetric configuration including two disjoint strip with equal width l which are separated by the distance l .Our result show that the behavior of the HMI and EoP can determine the type of the phase transition and hence we can use these entanglement measures to characterize the phase structures of strongly coupled matters.Moreover, at the second order phase transition, we focus on the near critical point and obtain the critical exponent using the EoP and HMI which is in agreement with the specific heat critical exponent.
Several problems call which we leave for further investigations.There are other quantum information quantities such as reflected entropy, odd entanglement entropy and logarithmic negativety and one can study them to probe various phase structures of the strongly coupled matters.It would also be interesting to study the quantum information quantities in the holographic QCD model at finite temperature including chemical potential which can generate more complete phase diagram of QCD [110].
background and the thermodynamics of the system A. The review of the background B. The thermodynamics 1.V QCD 2. V 2nd 3. V 1stIII.The holographic entanglement measuresA.The holographic mutual information B. The holographic entanglement of purification

FIG. 2 :
FIG. 2: The square of speed of sound c 2 s (a), entropy density s (b), specific heat Cv (c) and free energy F (d) v.s.temperature T for VQCD potential.

FIG. 3 : 5 FIG. 4 :
FIG. 3: The square of speed of sound c 2 s (a), entropy density s (b), specific heat Cv (c) and free energy F (d) v.s.temperature T for V 2nd potential.

FIG. 5 :
FIG. 5: The square of speed of sound c 2 s (a), entropy density s (b), specific heat Cv (c) and free energy F (d) v.s.T for V1st potential.The blue solid curves correspond to the stable black hole solution and the dashed cyan curves correspond to the unstable solution.

FIG. 8 :
FIG.8: Top row: The HMI v.s.temperature T for the dilaton potentials VQCD (left), V 2nd (middle) and V1st (right) for fixed l = 0.05 and l = 0.01.Bottom row: The HMI v.s.temperature T for the dilaton potentials VQCD (left), V 2nd (middle) and V1st (right) for fixed l = 0.05 and l = 0.005.In the right panels the solid blue curves correspond to the thermodynamically stable black hole solutions and the dashed cyan curves corresponds to the unstable solutions.

FIG. 9 :- 19 - 5 FIG. 10 :
FIG.9: Top row: The EoP v.s.temperature T for the dilaton potentials VQCD (left), V 2nd (middle) and V1st (right) for fixed l = 0.05 and l = 0.01.Bottom row: The EoP v.s.temperature T for the dilaton potentials VQCD (left), V 2nd (middle) and V1st (right) for fixed l = 0.05 and l = 0.005.In the right panels the solid blue curves correspond to the thermodynamically stable black hole solutions and the dashed cyan curves corresponds to the unstable solutions.

TABLE I :
The three scalar potential VQCD, V 2nd and V1st