Mixed-State Entanglement and Transport in Einstein-Maxwell-Axion-Horndeski Theory

We present a comprehensive study exploring the relationship between transport properties and measures of quantum entanglement in the Einstein-Maxwell-Axion-Horndeski theory. By using holographic duality, we study the entanglement measures, holographic entanglement entropy (HEE) and entanglement wedge cross-section (EWCS), and transport coefficients, for this model and analyze their dependence on free parameters which we classify into action parameter, observable parameters and axion factor. We find contrasting behaviors between HEE and EWCS with respect to observable parameters (charge and temperature), and the axion factor, indicating that they capture different types of quantum correlations. We also find that HEE exhibits positive correlation with both charge and thermal excitations, whereas EWCS exhibits a negative correlation with charge-related conductivities and thermal fluctuations. Furthermore, we find that the Horndenski coupling term, as the modification to standard gravity theory, does not change the qualitative behaviors of the conductivities and the entanglement measures.

However, calculating quantities related to quantum entanglement poses challenges due to the exponential increase in degrees of freedom and the large Hilbert space dimension in quantum systems [1,[9][10][11].
Recently, it is found in condensed matter theory that entanglement can serve as a distinct and novel tool to detect both thermal and quantum phase transitions [1,2,41].Quantum phase transitions occur at absolute zero temperature due to quantum fluctuations, whereas thermal phase transitions result from thermal fluctuations.The transport properties of a system, such as electrical and thermal conductivity, often provide insight into its underlying physics.Notably, many quantum and thermal phase transitions demonstrate unique signatures during these transport measurement processes, enabling us to identify the nature of the transitions.In the holographic duality theory, entanglement has also been generally recognized as a novel diagnostic tool for identifying both thermal and quantum phase transitions.Specifically, it has been established that during holographic quantum and thermal phase transitions, the entanglement measures of the system proves crucial in determining the nature of the transition [42][43][44][45].
While the gravity models without modification terms have been extensively studied, the exploration of modified gravity models, which could exhibit more interesting phenomena in entanglement, is still lacking.Modified gravity theories aim to solve issues from cosmological observations or gravity renormalization by incorporating correction factors to general relativity [46][47][48][49][50][51][52][53].One such theory is higher derivative gravity, which includes high-order derivatives of the Riemann curvature tensor [54][55][56][57].However, the complexity of solving higher-order partial differential equations involved in determining the entanglement measures hinders the exploration of modified models [58].
To address this challenge and investigate the relationship between transport properties and entanglement properties in modified gravity models, we adopt the Einstein-Maxwell-Axion-Horndeski model.The Horndeski model, which is the most general scalar-tensor theory involving a non-minimal coupling between a scalar field and the Einstein tensor, has gained attention as a cosmological model for dark energy and inflation [42,[59][60][61][62][63][64][65][66][67][68][69][70][71][72][73].This model is particularly interesting for several reasons.First, the Horndeski model is characterized by a non-minimal coupling between a scalar field and the Einstein tensor, which can significantly affect both the transport properties and the entanglement entropy prescription.This makes it a good candidate for exploring how modifications to gravity can impact these properties.Second, the inclusion of axions breaks the translation symmetry of the system, leading to finite DC conductivities.This feature makes the system more realistic and allows us to examine the transport properties in greater detail.Third, the analytical solution of this system enables us to calculate the DC conductivities analytically, and to work out entanglement-related quantities more easily.By starting with the Horndeski model with axion fields as our first step towards a more general study of transport properties and entanglement in the presence of modifications, we can gain valuable insights into how these phenomena are affected by non-minimal couplings and broken symmetries.
In this Einstein-Maxwell-Axion-Horndeski model, we adopt the entanglement entropy prescription of Dong [58] due to the presence of modification, i.e., the Horndenski coupling.This prescription accounts for the effects of the Horndenski coupling and simplifies the entanglement calculations by involving only ordinary second-order equations of motion.
Moreover, we extend the prescription to the EWCS by adjusting the definition of the minimum cross-section based on the modified entanglement entropy prescription.This allows us to explore the entanglement properties in the context of the Horndeski model with axions.
In the subsequent sections, we present the Horndeski spacetime model and introduce holographic information-related quantities, along with the calculation of the equation of motion and holographic thermodynamic conductivities (Section II).We study HEE and

EWCS and discuss their relation to transport properties (Sections III A, III B and III C).
Finally, we conclude with a summary in Section IV.

GRAPHIC QUANTUM INFORMATION
In this section, we start by introducing Horndenski gravity theories with axions and explaining their theoretical framework.Then, we delve into the concepts of HEE and EWCS, emphasizing their significance in quantum information theory.Finally, we discuss the com-putation of these Holographic Information-related Quantities, including the methods and techniques used to calculate them in gravity theories.
A. Einstein-Maxwell-Axion-Horndeski gravity model The four-dimensional Einstein-Maxwell-Horndeski gravity model with two free axions ϕ i can be described by the following Lagrangian density, where κ is a constant that we set to 1 for concreteness, and Λ is the cosmological constant, which acts as a variable in this model.In addition to the usual gravitational terms, this model includes a scalar field χ, with the derivative of the scaler filed χ coupled to the Einstein tensor G ab .ζ is the constant measuring the kinematic energy of the scalar field, and γ, is the Horndeski coupling.The constants ζ and γ are related by, where g = 1 l is the inverse of the AdS radius.It may seem that to turn off the coupling between the scalar field and the gravity, one simply needs to set γ and ζ to zero, returning the action to that of the AdS-RN-Axion theory.However, even in this limit, the modification does not vanish entirely.According to (4), a finite term will still remain due to the χ term with a 1/γ dependence.In other words, this modification will always persist and cannot be switched off.This is the characteristic of the Horndeski model.
The charged Horndeski planar black hole solution in asymptotic AdS spacetime is given by [64,70,74], where In these equations, β is a constant related to the scalar field χ, and q is charge of the black hole.The Hawking temperature and the cosmological constant are determined by, Here, r h is the radius of the horizon of black hole.In this article, we examine the adjustable parameters (ζ, γ, β, q, T, k) and group them into three categories: 1. Action parameters: These are the parameters ζ, γ, and β, which determine the action or theory of the system.
2. Observable parameters: These include the parameters T and q, which are related to observable quantities, the temperature and the charge of the dual system.
3. Axion factor: The parameter k, representing the axion factor, plays a crucial role in the transport properties.
To study the holographic dual of a strongly coupled field theory in Minkowski spacetime, it is necessary to redefine the coordinates as (t → 1 g 2 t, x → 1 g x, y → 1 g y) to ensure that the new metric of the asymptotic AdS spacetime takes the standard form at the boundary: Additionally, to simplify the calculations of the EWCS, a change of variable from r to z (r → r h z ) is performed.In this new coordinate system, the black hole horizon is located at z = 1, while the asymptotic boundary of the spacetime is at The thermodynamics of AdS black holes in Horndeski theory has been extensively studied in previous works [69,73].These studies used the Wald formula (14) to analyze the thermodynamic properties of black holes [75].The Wald entropy formula ( 14) is commonly used to analyze the thermodynamic properties of black holes [75].In this Horndenski theory, the entropy density s is given by: In addition to thermodynamics, various transport coefficients, known as conductivities, play an important role in understanding the response of a system to external electric and thermal perturbation.The direct current (DC) thermoelectric conductivities are commonly investigated, which involve constructing a radially conserved current linking the physical properties at the boundary to information present on the black hole horizon [64,[76][77][78][79][80].
This can be achieved through perturbating the system, allowing the derivation of radially conserved electric and heat currents on the horizon [42,64].Based on these calculations, the DC conductivity matrix reads, where σ DC and κDC are the DC electric and thermal conductivities, respectively, while α DC and ᾱDC are the DC thermoelectric conductivities.With these analytical expressions for the DC thermoelectric conductivities, we can easily determine how these quantities depend on the action Parameters (ζ, γ, β), observable parameters (q, T ), and axion Factor k.The relationships between the DC conductivities and parameters are summarized in the first three rows of TABLE I.
Remind that, according to [79], in the AdS-RN-Axion model,where the modification terms vanishes, the DC conductivities are given by, Therefore, it is deduced that, the relationship between the DC conductivities and parameters are similar to the case of the current model.
In light of the preceding overview of the EMAH model, we will now proceed to introduce the holographic quantum information measures.

Cross-section
The von Neumann entropy is a fundamental concept in quantum mechanics that quantifies the entanglement between a subsystem, denoted as A, and its complement in a pure state |ψ⟩ [81].It serves as the quantum counterpart of classical information entropy and is commonly referred to as the entanglement entropy (EE).The EE of subsystem A can be computed using the following formula, where ρ A represents the reduced density matrix of subsystem A, obtained by tracing out subsystem B from the total density matrix ρ A = Tr B (|ψ⟩⟨ψ|).
The holographic duality of entanglement entropy, known as holographic entanglement entropy (HEE), was initially proposed by Ryu and Takayanagi in 2006 based on the holographic principle, marking a significant breakthrough in the field of quantum information [7].In holography, the entanglement entropy is related to the geometry of a certain region in the gravitational theory.Specifically, the entanglement entropy S A of a certain subregion A in the quantum system is proportional to the area of a minimal surface that spans the boundary of region A in the gravitational theory, where γ A is the minimal surface and G N is the Newton's gravitational constant.
The EE measures the amount of entanglement between two subsystems in a pure state.
It quantifies the degree to which the state cannot be factored into separate states for each subsystem.For a pure state, the entanglement entropy is non-zero, while for a fully factorizable state (no entanglement), the entropy is zero.However, when dealing with mixed states, which are ensembles of different pure states, the entanglement entropy is not a good diagnostic since it can be non-zero even in the absence of entanglement.To overcome this limitation, one can use the mutual information (MI).The MI quantifies the amount of shared information between two subsystems, regardless of whether it arises from entanglement or from classical correlations.It is defined as the difference between the entanglement entropies of the combined region A ∪ B and the individual regions A and B: The MI provides a more complete picture of the correlations in a quantum system, as it captures both classical and quantum correlations.In the case of a pure state, the MI is equal to the entanglement entropy, but for mixed states, it can also capture classical correlations that do not involve entanglement.Hence, MI is not a perfect measure of entanglement in mixed states, as it can be dominated by the entanglement entropy.
An alternative technique is the entanglement of purification (EoP), which involves a double minimization procedure to purify the extra degrees of freedom of mixed states [24].The holographic EoP has been proposed as proportional to the area of the minimal cross-section of the entanglement wedge (EWCS, see Fig. 1) [31].Later, the EWCS, has been proposed to be dual to other mixed state entanglement measures such as logarithmic negativity [30,33].Building upon the earlier discussion on holographic quantum information measures, we will now delve into their computational methods in the Axion-Horndeski model.In this scenario, the entanglement entropy prescription undergoes substantial modifications, resulting in a more intricate calculation procedure.

C. Computation of the Holographic Information-related Quantities
Now, we will discuss the computation methods for holographic information-related quantities.Our discussion will start by explaining how to obtain the minimum surfaces for HEE and EWCS.Once we have established this foundation, we will then elaborate on the methods for determining the minimal cross-section within the entanglement wedge.
The entanglement entropy of Horndeski can be evaluated by substituting the Lagrangian density (1) into the entanglement entropy prescription proposed by Dong in [58], The Greek letters are indices of tensors in the d-dimensional spacetime.Several key terms are defined as Here, n µ and m µ are unit vectors perpendicular and transversal to the surface, respectively.The Latin letters (a), (b) and (i), (j) represent the internal indices in the tangent space.The induced metric on the tangent space is denoted by G ab , while ε ab and K aij represent the Levi-Civita tensor and the extrinsic curvature, respectively.The coefficient q ζ corresponds to an anomaly coefficient associated with each term in the expansion [58].
However, it is important to note that in the case of Horndeski Lagrangian density (1), only a first-order term of the Riemann tensor is present.Calculations using this Lagrangian density will yield results that are equivalent to the Wald entropy formula, which is given by, By extracting each curvature-related term from the Lagrangian density, one obtains the following expression for the entanglement entropy [82], The h ij represents the induced metric on the surface.We can treat (15) as an action and obtain the equation of the minimum surface through the variational principle However, the explicit expression of equation of motion of ( 16) is quite complicated.To simplify our computations, we follow the computational methods in [29] and consider an infinite strip along the y axis as the subregion (see Fig. 2).By parameterizing the surface with θ as θ = arctan(z/x), we can derive the equations of motion more easily.To obtain the minimum surface, we need to solve the equations of z(θ) and x(θ) within the range of θ ∈ [0, π].However, since there is a mirror symmetry about θ = π/2, it suffices to work within θ ∈ [0, π/2].We treat the integrand in equation ( 15) as a Lagrangian density of the minimum surface.Resultantly, the minimum surface can be obtained by solving the equations of motion from this Lagrangian density, Here, z ′ and x ′ are derivatives with respect to θ.By taking the variational derivative of (17) with respect to z(θ) and x(θ) directly, we can derive the equations of motion for z and x. Since these equations are highly non-linear ordinary differential equations, solving them can be challenging.One approach is to use the Chebyshev collocation and Newton iteration algorithm.This iterative method is commonly used to solve boundary value problems and involves approximating the solution by a polynomial that satisfies the differential equations at a set of collocation points.The algorithm then updates the polynomial coefficients using Newton's method until a satisfactory solution is obtained.Although this method can be computationally intensive, it has proven to be effective in solving similar types of problems.
After obtaining the minimum surface, the HEE can be calculated through numerical integration using the Lagrangian density (17) of the minimum surface.However, the HEE diverges at the AdS boundary (z = 0), so we need to identify and subtract the diverging term to obtain a finite value.The diverging term is given by, which needs to be subtracted so as to render a finite value of HEE.
Next, we elaborate on how to compute EWCS.To determine the EWCS, we first identify the minimum cross section within the entanglement wedge.This involves finding two minimum surfaces and then locating the minimum cross-section that connects them.In order to parametrize the minimum cross-section, we use the variable z, which has been verified as a good parametrization in our numerics.The Lagrangian density (17) under this parametrization is given by, The equation of motion for x(z) is obtained by variating (19) with respect to x(z).In our previous work [41], we proposed an algorithm to identify the minimal cross-section based on the geometric fact that the minimal cross-section must be perpendicular to the minimum surfaces.This observation significantly speeds up the computation process.However, when considering the modified Lagrangian where the scalar field also takes part in the dual information-related quantities, we need to examine whether the perpendicular condition still holds.Essentially, the perpendicular condition comes from the term and is valid for any form of Lagrangian once this term appears as a factor.Therefore, the algorithm in [41] still applies to the current situation.Finally, we can work out the area of this minimum cross-section by integrating (19) over z.
Besides the computational methods discussed earlier, a primary focus of this paper is on exploring the thermodynamic and transport properties of the Horndenski gravity model with axions and their connections to holographic information-related quantities.In the subsequent section, we provide detailed explanations of these holographic information-related quantities and thermodynamic and transport properties.Second, when considering the observable parameters, black hole charge q and the Hawking temperature T , it has been observed that the HEE monotonically increases with both of these parameters.This observation is depicted in Fig. 4.This behavior can be understood from both the gravity side and the dual CFT side.From the gravity side, the increase in HEE with black hole charge and temperature can be attributed to the growth of the black hole horizon.As the charge q increases, the horizon area expands, leading to an increase in the Bekenstein-Hawking entropy.The HEE, which is proportional to the area of the minimal surface in the bulk theory holographically dual to the boundary CFT, follows a similar trend and thus increases with the charge.Similarly, an increase in temperature leads to a larger black hole horizon and subsequent increase in the HEE.This result is consistent with our expectation that larger black holes have more entanglement entropy.From the dual CFT side, the monotonically increasing HEE with both charge and temperature can be understood by considering the behavior of the dual conformal field theory.The HEE measures the quantum entanglement between degrees of freedom in the CFT.As the charge and temperature increase, the system becomes more excited, leading to an enhanced entanglement between its constituents.This increased entanglement is reflected in the larger value of the HEE.
The HEE versus the axion factor k is depicted in Fig. 5.As observed, HEE exhibits a monotonically increasing trend with increasing k.In the context of holographic duality, k is commonly associated with the transport properties of the system.By breaking the translational symmetry, k allows for the emergence of finite DC conductivities.This implies that as k is increased, the system experiences a smaller DC resistivity, accompanied by an increase in HEE.From the perspective of gravity, it can be argued that the increase in HEE with k is due to the modification in the gravitational background caused by the axion field.Similar results have been observed in [41] and the mechanisms for this monotonic behavior can be attributed to the relationship between the horizon radius and the parameter k.On the dual field theory side, the enhanced EE can be attributed to the extra degrees of freedom aside from the confined charged degrees of freedoms, such as the enhanced lattice vibrations.
Generally, we would expect HEE to exhibit the same pattern as the thermodynamics entropy ( 6) when the width of the subregion is large.All these behaviors have been summarized in TABLE I.

B. Entanglement Wedge Cross-section
We will now investigate how varying different parameters affects the EWCS, while keeping the configuration (a, b, c) fixed.It is important to note that the results we obtained are independent of the specific choices of (a, b, c).Therefore, we will only present a few characteristic plots.In the case of EWCS, g is a free parameter, allowing us to consider γ and ζ as two independent variables.Next, we examine the behavior of the EWCS in relation to the action parameters ζ, γ, β, as well as the observable parameter q, T , and axion factor k.
Let us start with ζ, which is related to the magnitude of the kinetic energy of the scalar field.We observe that as ζ increase, the EWCS decreases (see the left plot of Fig. 6).
From the perspective of gravity, this can be understood as a consequence of the scalar field possessing more kinetic energy.With a larger ζ, the gravitational effects induced by the scalar field become stronger.This increased gravitational impact leads to modifications in the spacetime geometry, resulting in a decrease of the EWCS and reducing the entanglement across the entanglement wedge.Moving on to γ, which represents the strength of the coupling between gravity and the scalar field, we find that the EWCS increases as γ increases (see the middle plot of Fig. 6).
From the gravity side, this can be understood as a consequence of stronger gravitational effects due to the scalar field.As γ increases, the coupling between the scalar field and gravity becomes stronger, resulting in enhanced entanglement across the entanglement wedge.
Consequently, the EWCS exhibits an increase.
Lastly, we consider the effect of β, which is embedded within the definition of the scalar field χ.As β increases, we observe a decrease in the EWCS (see the right plot of Fig. 6).
Remind that, HEE exhibits a monotonically increasing behavior.This suggests that HEE exhibit a opposite behavior with β when compared with EWCS.This means that, EWCS as a mixed state entanglement measure, captures distinct features of the entanglement from the EE.Now, we study how the observable parameters T and q affects the EWCS.The EWCS is found to decrease with increasing temperature T and charge q, as depicted in Fig. 7.
This phenomenon is in contrast to the behavior of HEE which typically increases with temperature and charge.
From a gravitational viewpoint, we can gain an intuitive understanding of the EWCS based on the reasoning behind why HEE increases with q and T .As temperature and charge increases, the black hole horizon expands, causing HEE to resemble thermal entropy more closely.In thermal systems, different subregions should be disentangled, as the total density matrix can be approximately represented by the direct product of each subregion's density matrix.This concept is easily understood from a gravity perspective: when both minimum surfaces approach the horizon, the cross-section becomes smaller due to the narrowing entanglement wedge.
On the dual CFT side, a possible understanding for the decreasing EWCS with increasing temperature and charge can be derived from the behavior of entanglement in strongly coupled systems.At higher temperatures, the thermal fluctuations in the system become more significant.These fluctuations can disrupt or break the entanglement between particles, leading to a decrease in overall entanglement.Increased temperatures generally introduce more randomness and decoherence into the system, which can result in the loss of entanglement.Additionally, the presence of a higher charge may lead to stronger interactions between particles.These interactions can also contribute to decoherence and disrupt the entanglement present in the system.The increased charge density may cause stronger repulsion, which can break or reduce the entanglement between particles.One example where charge density and interactions can affect entanglement is in the context of the Hubbard model, which describes interacting particles on a lattice.In this model, increasing charge density can lead to stronger on-site interactions between particles, which can influence the entanglement properties of the system.
The EWCS versus the axion factor k also depicts a similar trend as the HEE.It exhibits a monotonically increasing behavior with increasing k, as can be seen in Fig. 8.This observation suggests that mixed state entanglement measures, such as the EWCS, also show an enhanced entanglement when the system possesses a larger axion factor k. From the gravity side, the increase in EWCS can be understood as a consequence of the modification in the gravitational background caused by the axion field.The interplay between the axion field and the gravitational dynamics affects the entanglement characteristics of the system, leading to a higher value of EWCS as k increases.This phenomenon has been studied in [41], where a similar monotonic behavior of the entanglement measures with k was observed.
The specific mechanisms underlying this relationship between the axion factor and the entanglement measures may be attributed to the interplay between the horizon radius and the parameter k, as well as other relevant properties of the gravitational background.On the dual field side, the enhanced mixed state entanglement can be attributed to the effects of localization and confinement resulting from the breaking of translational symmetry by the axion field.The breaking of translational symmetry induces a more localized and confined behavior for charged degree of freedom of the quantum states, meanwhile also induces a enhanced extra degrees of freedom such as lattice vibration, resulting in an increased mixed state entanglement.
In the study of modified gravity models, it has been observed that the Horndeski model with axion and the AdS-RN-Axion model (without modification) exhibit similar behaviors in terms of their DC conductivities.This means that the inclusion of the modification term does not significantly alter the DC conductivities in these models.Furthermore, it has also been found that the HEE and EWCS are quantitatively the same for both models2 .These findings suggest that when axions are minimally coupled to gravity, the modification term does not have a significant impact on the quantum information, such as the DC conductivities and HEE/EWCS, and their relationships.The relationship between transport properties and entanglement related quantities is a topic of great interest in the study of condensed matter physics.Transport properties, such as DC electrical conductivities, play a crucial role in characterizing the behavior of materials, distinguishing between metals and insulators, and probing phase transitions.
Recent research has revealed intriguing connections between transport properties and entanglement properties, particularly in the context of condensed matter systems and holographic frameworks.In certain scenarios, the entanglement properties exhibit extremal behaviors across phase transitions and in strongly correlated systems.These findings suggest that there may be underlying connections between the transport properties and entanglement properties.
In the current Horndeski model with axion, we also systematically analyze the relationship between transport properties and entanglement in this model.To illustrate this more clearly, we present all the results in Table I.In Table I, the rows and columns represent physical quantities and the free parameters of the system, respectively.In the columns, for better comparison, we highlight the action parameter part in blue, observable parameters in orange, and the axion factor in green.In the rows, the first three rows represent the trans-port coefficients of the system, while the last three rows represent the quantities related to entanglement and the thermodynamic entropy density.
Let us now examine the observable parameters, namely q and T .Increasing the value of q leads to an increase in charge-related conductivities (σ DC and α DC ), while the thermal conductivities decrease.This can be understood in the context of condensed matter systems, such as insulators, where the primary mode of heat transfer is through lattice vibrations (phonons) rather than free electron movement.In such systems, an increase in charge density does not enhance thermal conductivity.In fact, high charge densities can even decrease thermal conductivity by increasing electron-electron scattering, which hinders the flow of heat.Additionally, increasing q is always associated with higher EE.This is because higher charge densities result in more charged degrees of freedom, leading to increased entanglement as well as the thermal entropy, as explained in previous sections.However, the mixed state entanglement EWCS, decreases with increasing q, which can result from the localization induced by strong correlation between electrons.Meanwhile, the thermal conductivities exhibit the opposite trend.This is understandable since enhanced heat transfer tends to disrupt entanglement.
For the temperature T , the heat-related DC conductivities (κ DC and α DC ) increase with increasing T , while σ DC decreases.This behavior can be explained by the fact that higher temperatures facilitate more efficient heat transfer, resulting in larger values of κDC and α DC .The decreasing behavior of σ DC indicates that the system is in a typical metallic phase.
Turning our attention to the entanglement-related quantities, we find that increasing the temperature leads to an increase in the entanglement entropy (HEE) due to the contribution of thermal entropy, while the EWCS decreases due to the disruption of thermal effects.
Based on the above observations, it appears that the HEE is aligned with the thermoelectric conductivity α DC .Which means that HEE detects both the heat and charged degrees of freedoms.Meanwhile, EWCS is negatively related to the charged degrees of freedom in the sense that, on the q direction, σ DC and α DC both increases with increasing charge density, while EWCS decreases due to the strong correlation of the charged degrees of freedom.Also, EWCS is aligned with σ DC on the temperature direction, but the mechanisms are different at that, σ DC decreases as a metallic charge transport behavior, but the thermal effects tends to disrupt the mixed state entanglement, i.e., EWCS.
Next, let us focus on the axion factor k. All three DC conductivities decrease as the value of k increases.This decrease is largely attributed to the process of momentum relaxation, which not only affect momentum but also impacts the flow of charge and energy.Keep in mind that this pattern has been observed in several other models, thereby suggesting that it might be a common occurrence within such axion structures, particularly within those that are minimally coupled.Conversely, as k increases, the entanglement-related quantities such as the HEE, and the EWCS all increase.This seems counterintuitive, as a reduction in the flow of charge -especially when charged degrees of freedom are predominant -tends to significantly decrease entanglement.We present this as evidence that in this model, despite decreasing charge mobility and associated momentum and energy, the charged degrees of freedom does not dominate the system's entanglements.It suggests that the augmentation of other degrees of freedom, such as phonon, takes place, resulting in higher entanglement.

IV. DISCUSSION
This work presents a comprehensive study of holographic entanglement entropy (HEE) and entanglement wedge cross-section (EWCS) in the Horndeski gravity model with axions.Leveraging holographic duality, we investigate the interplay between these measures of quantum entanglement and the transport properties described by DC thermoelectric conductivities.Our analysis reveals several key insights: • The HEE increases with the action parameters ζ, β, black hole charge q, and temperature T .This aligns with expectations from black hole thermodynamics, where entropy increases with charge and temperature.The EWCS exhibits the opposite trend, decreasing with q and T .
• The contrasting behavior of HEE and EWCS suggests that they capture different types of entanglement.HEE appears more closely related to the degrees of freedom associated with charge and thermal excitations.In contrast, EWCS seems more influenced by disruptive effects like thermal fluctuations that can destroy quantum correlations.
• All DC conductivities decrease with the axion factor k, indicating momentum relaxation effects that hinder charge and energy transport.However, both HEE and EWCS increase with k, highlighting the complex relationship between transport and entanglement.
• The overall trends indicate that EWCS shows a negative relation with the chargerelated conductivity σ DC and α DC , while HEE parallels the thermoelectric conductivity α DC .This implicates EWCS is negatively related to charged degrees of freedom, while HEE reflects both charge and thermal excitations.
These findings elucidate the essential properties of the Horndeski holographic model and demonstrate the value of entanglement measures in decoding transport phenomena.The contrasting behavior of HEE and EWCS with respect to various parameters emphasizes the need for multiple diagnostic tools to fully characterize quantum information.Moreover, our findings reveal that the dependence of these quantities on the system parameters remains unchanged when comparing scenarios with and without modification terms.This suggests that the observed behavior is a general characteristic of axions minimally coupled to gravity.
While this work focuses on a specific model, the methodology of correlating entanglement and transport properties can offer broader physical insights.Extending the analysis to other holographic models could reveal universal patterns tied to phenomena like momentum relaxation, metal-insulator transitions, and electron correlations.An interesting extension will be to the modification that the axions are not minimally coupled to the Horndenski gravity [84].We expect that this non-minimal coupling can significantly impact the relationship between the entanglement and transport properties.Moreover, modifications like higher derivative terms are also an important topic to further study.However, this could bring in intricate numerical challenges such as the calcultion of fourth-order differential equations.
We are working on these directions and will report our findings in future publications.

FIG. 1 :
FIG.1: Visualization of the minimal cross-section (green surface on the left plot and red string on the right plot) that cuts through the entanglement wedge.

FIG. 2 :
FIG. 2: A visualization of the minimal surface (yellow) corresponding to an infinite strip (red) along the y axis.

FIG. 3 :
FIG. 3: The influence of action Parameters ζ and β on holographic entanglement entropy (HEE).Left: Variation of HEE with different values of ζ.Right: Variation of HEE with different values of β.

FIG. 4 :
FIG.4:The influence of observable parameters q and T on HEE.Left: Variation of HEE with different values of q.Right: Variation of HEE with different values of T .

1 , 1 , 15 ,FIG. 6 :
FIG. 6: The influence of action Parameters ζ, γ and β on the EWCS.Left: Variation of EWCS with different values of ζ.Right: Variation of EWCS with different values of γ.Right: Variation of EWCS with different values of β.

TABLE I :
The dependence of conductivities, EWCS, HEE, and thermodynamic entropy on variations of different parameters.An upward arrow indicates an increase with the parameter, while a downward arrow indicates a decrease with the parameter.C. The Relationship Between the Transport Properties and the Entanglement-