\section{Introduction}

Over the last years, studies on the black hole (BH) entropy have revealed  close connections between the  thermodynamic properties and the event horizon of a BH \cite{Hawking:1971vc,Bardeen:1973gs,Bekenstein:1973ur}. Among all different kinds of BHs provided with different geometries and thermodynamic features, Schwarzschild BH represents the simplest case, where part of the radiation is absorbed by the BH mass \cite{Hawking:1975vcx,Wald:1999vt}. Other interesting BH solutions include Reissner-Norstr\"om BH, whose thermodynamics is similar to that of regular BH \cite{Myung:2007av, Ghaffarnejad:2013wpa,Good:2020qsy}, Ho\v{r}ava-Lifshitz BH \cite{Cai:2009qs,JahaniPoshteh:2021clv}, characterized by rich thermodynamic properties  and so forth.

Moreover, asymptotic BH solutions to the de Sitter space can be obtained from Einstein's equations with a positive cosmological constant $(\Lambda>0)$ \cite{Gibbons:1977mu,Liu:2016urf}. In such a case, it has been shown that the surface gravity of the BH horizon would determine the temperature of particles emitted from the BH \cite{Pappas:2017kam}. The same, however, happens even for the cosmological event horizon, so that a thermal equilibrium may occur only if the two surfaces coincide \cite{Davies:1988dma,Chatterjee:2020gfc,Donnay:2019zif}. Further, BHs that are asymptotic to the anti-de Sitter (AdS) space can be found as solutions to Einstein's field equations with $\Lambda<0$ \cite{Anabalon:2012tu}. As in the case of the asymptotically flat space, the entropy and temperature of AdS BH are equal to 1/4 of the event horizon area, whereas, differently from the flat space case, such objects admit, at a given temperature, a stable equilibrium with radiation, and a positive specific heat \cite{Hawking:1982dh}. Throughout recent years, the physics of asymptotically AdS BHs has gained a renewed interest due to the AdS/CFT duality \cite{Maldacena:1997re,Lunin:2001jy,Hubeny:2014bla}. In this context, particular attention was given  to the study of thermal field theories living on the AdS boundary and, from the bulk perspective, to the several phase transitions that these types of BH exhibit. 

Furthermore, in treating the cosmological constant as the thermodynamic pressure, it has been shown that the thermodynamics of a charged/rotating AdS BH exactly coincides with that of Van der Waals's fluid \cite{Rajagopal:2014ewa,Roy:2021ucl,Capozziello:2004ej}. Subsequently, an asymptotically AdS solution to Einstein’s equations was obtained by matching the BH thermodynamic parameters with those of a particular class of polytropic gas \cite{Setare:2015xaa}. Then, an additional AdS BH solution was found in the thermodynamic framework of modified Chaplygin gas \cite{Debnath:2019mzs,Kamenshchik:2001cp}.

Motivated by those findings, in the present study, we focus on the class of logotropic models, whose  thermodynamic features permit to heal astrophysical issues related to dark matter distribution in galaxies, and unify the cosmological dark sector \cite{Chavanis:2015paa,Chavanis:2015eka,Chavanis:2016pcp,Benaoum:2021pqf}. As a prominent result, these models can be generalized to the well-known Anton-Schmidt fluid \cite{ANTON1997449,ASfluid}. These scenarios have been recently proposed in the cosmological context as a unified dark energy model \cite{Capozziello:2017buj,Capozziello:2018mds,Boshkayev:2019qcx,Odintsov:2018obx}\footnote{For alternative approaches to dark energy, see also \cite{DAgostino:2019wko,DAgostino:2021vvv,Capozziello:2019cav,Capozziello:2022wgl,Capozziello:2022rac,Bajardi:2022tzn,DAgostino:2022tdk}.}. In this respect, the Anton-Schmidt fluid has been also studied in the Tolman-Oppenheimer-Volkov formalism \cite{Tolman:1939jz,Oppenheimer:1939ne} to obtain analytical solutions for a static and spherically symmetric BH \cite{Capozziello:2020szy}. In particular, from the relation between the Anton-Schmidt free parameters and the BH mass, one can find spacetime solutions describing Schwarzschild-de Sitter BH and naked singularities. Thus, it appears natural to investigate the thermodynamic consequences to check whether the inclusion of a logotropic and/or Anton-Schmidt's equation of state (EoS) may lead to reasonable results in the BH description.

In this paper, we search for a BH solution to Einstein's field equations, whose corresponding thermodynamics coincides with that of logotropic models. Starting from Anton-Schmidt's EoS, we propose a general treatment to obtain an asymptotic Schwarzschild-AdS metric, which reproduces the thermodynamic properties of the involved fluids. To determine the most suitable metric functions, we present a general method involving any possible density term contribution. Moreover, in order to have a physical BH, we invoke the existence of an event horizon and investigate under which circumstances the  energy conditions may hold. We also explore the possibility of violating the strong energy condition, in order to find additional physical properties. We then discuss the physical consequence of this recipe in view of the free constants emerging from the integration procedure. With the aim of distinguishing among different thermodynamic BHs, we consider the optical properties of our solutions and show how to construct an effective refractive index, following the standard procedure adopted for static and spherically symmetric spacetimes. In particular, we show that the net dependence of the refractive index on the underlying spacetime can lead to different outcomes.  The refractive index increases significantly under the choice of particular constant values, whereas the asymptotic regime is investigated in terms of density, showing the limit to Schwarzchild-AdS. Hence, we explore possible departures with respect to the case without the logotropic fluid, corresponding to a pure Schwarzschild-AdS case.

The paper is organized as follows. After this introduction, in Sec.~\ref{sec:A-S BH} we introduce Anton-Schmidt's EoS and its limit to logotropic models. There we postulate the metric ansatz for a static, spherically symmetric metric that is consistent with an asymptotic AdS spacetime. We thus analyze the thermodynamic properties of Anton-Schmidt's BH in terms of its mass, temperature and entropy. In Sec.~\ref{sec:energy conditions}, we constrain BH solutions requiring the presence of an event horizon and checking the validity of the energy conditions. In particular, we discuss how the violation of the strong energy condition may lead to a metric solution containing a factor that can be associated with a refractive index. Finally, in Sec.~\ref{sec:conclusions}, we summarize our findings and draw the conclusions of our work. In this study, we use Planck units $c=\hbar=G=1$.


\section{Logotropic black holes}
\label{sec:A-S BH}

Let us start by considering Einstein's field equations with the cosmological constant in the form
\begin{equation}
G_{\mu\nu}+\Lambda g_{\mu\nu}=8\pi T_{\mu\nu}\,,
\label{eq:EFE}
\end{equation}
where $G_{\mu\nu}\equiv R_{\mu\nu}-\frac{1}{2}Rg_{\mu\nu}$ is the Einstein tensor, $g_{\mu\nu}$ is the metric tensor, and $T_{\mu\nu}$ is the stress-energy tensor of the source fluid.
According to recent studies 
\cite{Kubiznak:2012wp,Gunasekaran:2012dq}, in the extended phase space one can interpret $\Lambda$ as a thermodynamic pressure, namely\footnote{Alternatively, it is possible to work out the same recipe adopting the conjugate variable of pressure, namely the volume \cite{Kubiznak:2012wp}.}
\begin{equation}
P=-\dfrac{\Lambda}{8\pi}=\frac{3}{8\pi l^2}\,,
\label{eq:thermo P}
\end{equation}
where $l$ is the AdS curvature constant.
Our aim is to construct an asymptotic AdS BH whose thermodynamics matches that of Anton-Schmidt's fluid with pressure given by \cite{ASfluid,Capozziello:2017buj}
\begin{equation}\label{ASfluid}
    P=A\left(\frac{\rho}{\rho_\ast}\right)^{-n}\ln\left(\frac{\rho}{\rho_\ast}\right),
\end{equation}
where the density $\rho$ is normalized to a reference density $\rho_\ast$, while $A>0$ and $n\neq -1$ are constants. 
The strategy is  starting from Eq.~\eqref{ASfluid}, which generalizes the logotropic models with $n=0$. In so doing, we recover the logotropic thermodynamics as a limiting case of the Anton-Schmidt fluid\footnote{For the sake of completeness, generalized versions of logotropic models can be as well investigated \cite{Capozziello:2018mds,Benaoum:2021pqf,Chavanis:2022vzi}. However, the latter have been recently criticized in \cite{Boshkayev:2021uvk}.}. 

We shall thus consider the static and spherically symmetric line element
\begin{equation}
ds^2=f(r,\rho)dt^2-\frac{dr^2}{f(r,\rho)}-r^2d\Omega^2\,,
\label{eq:metric}
\end{equation}
where 
\begin{equation}\label{Cgeneralization}
    f(r,\rho)=-\frac{2M}{r}+\frac{r^2}{l^2}-h(r,\rho)\,.
\end{equation}
Here, $M$ is the BH mass, and $h(r,\rho)$ is an unknown auxiliary function to determine. 
Clearly, in the case of $h(r,\rho)\rightarrow0$, the postulated metric represents an asymptotic AdS spacetime. In this scheme, the term $h(r,\rho)$ represents a correction when it is considered different from zero. 
Using Eq.~\eqref{eq:thermo P}, we can rewrite Eq.~\eqref{Cgeneralization} as
\begin{equation}\label{metricAnsatzgeneral}
    f(r,\rho)=-\frac{2M}{r}+\frac{8}{3}\pi r^2 P-h(r,\rho)\,.
\end{equation}
In what follows, we seek a class of metrics such that the BH thermodynamics predicted by the latter coincide exactly with Anton-Schmidt's EoS, containing the logotropic models in the limit $n=0$. 
The cosmological pressure $P$, affecting the BH thermodynamics, can be therefore associated with the pressure of the fluid $P(\rho)$. In this way, the function $h(r,\rho)$ accounts for the metric correction that occurs by considering an EoS different from the case of a pure cosmological constant\footnote{Although degenerating with a pure cosmological constant, the case of dark fluid \cite{Luongo:2018lgy,DAgostino:2022fcx,Belfiglio:2022egm} has not been explored here.}, namely $P=-\rho$.


\subsection{Black hole thermodynamics}\label{BHth}

To find the most suitable form of $h(r,\rho)$, we start from the standard BH entropy in terms of the horizon radius $r_h$ and area $A$ as \cite{Bekenstein:1973ur,Bekenstein:1974ax,Bekenstein:1975tw}
\begin{equation}
S=\frac{A}{4}=\pi r_h^2\,.
\label{eq:entropy}
\end{equation}
We can thus relate the BH thermodynamic properties to the parameters of the Anton-Schmidt and logotropic fluids. In particular, the BH mass could be obtained from the definition of horizon radius, namely $f(r_h,\rho)=0$:
\begin{equation}\label{mass}
    M=\frac{4}{3}\pi r_h^3 P-\frac{r_h}{2}h(r_h,\rho)\,.
\end{equation}
As the conjugate variable of the pressure, the volume can be calculated as
\begin{equation}\label{volume}
    V=\frac{\partial E}{\partial P}=\frac{\partial M}{\partial P}=\frac{4\pi r_h^3}{3}-\frac{r_h}{2}\frac{\partial h(r_h,\rho)}{\partial \rho}\left(\frac{\partial P}{\partial \rho}\right)^{-1}.
\end{equation}
According to \cite{Hawking:1975vcx}, the temperature is obtained by studying the surface gravity of the BH as
\begin{equation}\label{temperature}
    T=\frac{1}{4\pi}\frac{\partial f(r,\rho)}{\partial r}\bigg|_{r_h}=2r_h P-\frac{1}{4\pi r_h}\frac{\partial (rh(r,\rho))}{\partial r}\bigg|_{r_h}\,.
\end{equation}
From the first law of thermodynamics, $dE=TdS-PdV$, with $dE=dM=\rho dV$, and assuming the following integrability condition
\begin{equation}
\frac{\partial^2 S}{\partial T \partial V}=\frac{\partial^2 S}{\partial V \partial T}\,,
\end{equation}
one finds
\begin{equation}\label{int1stlaw}
    S=\left(\dfrac{\rho+P}{T}\right)V\,.
\end{equation}
Therefore, plugging Eqs.~\eqref{ASfluid}, \eqref{eq:entropy}, \eqref{volume} and \eqref{temperature} into Eq.~\eqref{int1stlaw} and considering the solution for a generic $r\geq r_h$, we obtain
\begin{equation}\label{diffeq}
    8\pi r^2(P-2\rho) P_\rho+6(\rho+P) h_\rho-3(rh)' P_\rho=0\,,
\end{equation}
where the subscript $\rho$ and the prime  denote the partial derivatives with respect to the density and radial coordinate, respectively. 
Starting from the theoretical setup presented in \cite{Rajagopal:2014ewa}, we seek a solution of Eq.~\eqref{diffeq} for a  generic $r$ by implementing a general method that makes use of  combinations of linearly independent functions of the density. In particular, a similar approach has been employed in  \cite{Setare:2015xaa, Debnath:2019mzs}, but imposing \emph{a priori} the functional expressions for $R_i(\rho)$. In our treatment, we relax this hypothesis  as illustrated in more detail in appendix \ref{appendixA}. We thus write
\begin{equation}\label{addLinComb}
    h(r,\rho)=\sum_{i}X_i(r)R_i(\rho)\,,
\end{equation}
where the coefficients $X_i$ depend on the radial coordinate, $r$. In this way,  Eq.~\eqref{diffeq} takes the formal expression
\begin{equation}\label{eqs}
    \sum_j\xi_j(r)F_j(\rho)=0\,,
\end{equation}
where $j$ labels the linearly independent density functions, $F_j(\rho)$, accounting for the information on the fluid EoS, and  the relative weights, $\xi(r)$. 
In this way, we obtain a system of equations, whose solution provides the coefficients of $h(r,\rho)$, i.e. $X_i(r)$. 
We report all the calculations in the case of Anton-Schmidt's fluid in appendices \ref{appendixA} and \ref{appendixB}. The final result leads to
\begin{equation}
\label{addtermfinal}
    h(r,\rho)=\frac{c_1}{r}+\frac{8}{3}\pi r^2P(\rho)+c_2r^{\frac{2}{B_2}-1}(a(\rho))^{1/B_2}\,.
   
\end{equation}   
where $c_1$, $c_2$ and $B_2$ are free parameters, whereas  
\begin{equation}
    a(\rho)\equiv \exp\left(\int \frac{ P_\rho }{P+\rho} \,d\rho\right).
    \label{eq:a}
\end{equation}
Clearly, the above solution leads to a $f(r,\rho)$ that, once plugged into Eq.~\eqref{Cgeneralization}, implies a solution to the vacuum field equations plus a constant term that mimics the constant arising from the AdS assumption made above. 

In the following, we analyze the solutions resulting from Eq.~\eqref{addtermfinal} and the relative energy conditions.


\section{Metric solutions and energy conditions}
\label{sec:energy conditions}

Inserting $h(r,\rho)$ from Eq.~\eqref{addtermfinal} into Eq.~\eqref{metricAnsatzgeneral}, we find the metric solution
\begin{equation}\label{generalmetric}
    f(r,\rho)=-\frac{2M+c_1}{r}-c_2r^{\frac{2}{B_2}-1}a(\rho)^{1/B_2}\,.
\end{equation}
Here, $c_1$ has the only effect to modify the value of the BH mass without affecting the underlying physics\footnote{Even if we do not consider changing the mass sign, repulsive gravity effects are possible in arbitrary spacetimes, see e.g. \cite{Luongo:2014qoa,Luongo:2015zaa,Giambo:2020jjo}.}, so that we can safely set  $c_1=0$. Also, to lighten the notation, we redefine $B_2=\beta$.

In order to obtain an event horizon, we impose the condition $f(r_h,\rho)=0$, namely
\begin{equation}\label{radiusAS}
    r_h=\left(\frac{2M}{-c_2}\right)^{\frac{\beta}{2}}\frac{1}{\sqrt{a(\rho)}}\,,
\end{equation}
requiring the condition $c_2<0$. Therefore, for simplicity, hereafter we set $c_2=-1$. 

A further condition is
$f(r,\rho)\rightarrow 0$ as approaching the horizon from outside, i.e. for $r\rightarrow r_h$. 
We thus require $f(r,\rho)>0$ in the range $r>r_h$, namely  $f(r_h+\epsilon,\rho)>0$ for very small $\epsilon$: 
\begin{equation}
    0<-\frac{2M}{r_h+\epsilon}+(r_h+\epsilon)^{\frac{2}{\beta}-1}a(\rho)^{1/\beta}
    \sim f(r_h,\rho)+\epsilon\left[\frac{2M}{r_h^2}+r_h^{\frac{2}{\beta}-2}\left(\frac{2}{\beta}-1\right)a(\rho)^{1/\beta}\right]. 
\end{equation}
Since $f(r_h,\rho)=0$ by definition, the last inequality becomes
\begin{equation}\label{condpositivehorizon}
  \frac{2M}{r_h^2}+r_h^{\frac{2}{\beta}-2}\left(\frac{2}{\beta}-1\right)(a(\rho))^{1/\beta} >0\,.
\end{equation}
Now we shall insert the expression for the horizon radius \eqref{radiusAS} into \eqref{condpositivehorizon}, obtaining
\begin{equation}
\begin{split}
    &(2M)^{1-\beta}a(\rho)+(2M)^{1-\beta}\left(\frac{2}{\beta}-1\right)a(\rho)>0\,,
\end{split}
\end{equation}
which implies $\beta>0$ in order to have a BH.

It is worth mentioning that the horizon existence is mathematically plausible, although, physically,  the corresponding compact object can also be different from a BH\footnote{In other words, any BH is provided with a horizon, but not all objects exhibiting horizons are BHs.}. To better clarify this point, one may consider the Hartle-Thorne metric \cite{Papapetrou:1966zz}, whose well-known horizon  describes massive and compact stars \cite{Stephani:2003tm}. However, the size of these objects is usually larger than the size of the event horizon, so the latter has not a precise meaning. Moreover, the Hartle-Thorne metric possesses its internal counterpart, which is meaningless for the event horizon, albeit mathematically one can always calculate it in analogy to the Sun and Earth. For the aforementioned reasons, we intend to clarify that our BH solution is that of a thermodynamic BH, i.e. a BH-like counterpart exhibiting the same thermodynamics of logotropic models.
Nevertheless, this does not limit the validity of our solution. In principle, we can  apply our outcomes to some precise cases, e.g. compact objects, quite different from genuine BH configurations. In fact, our recipe may be extended to other cases that may require, in astrophysics, an event horizon to determine.

\subsection{Matching the energy conditions}


In order to guarantee the energy conditions to hold and to investigate their effects in the framework of our solution, we shall now  compute the stress-energy tensor induced by Eq.~\eqref{eq:metric}. In the tetrad formalism, the stress-energy tensor can be written as
\begin{equation}
    T^{\mu\nu}=\varrho e_0^\mu e_0^\nu+\sum_{i=1}^3p_ie_i^\mu e_i^\nu\,,
\end{equation}
where $\varrho$, $p_1$, $p_2$ and $p_3$ are \cite{Rajagopal:2014ewa,Setare:2015xaa,Debnath:2019mzs} 
\begin{align}
    \varrho&=-p_1=\frac{1-f-rf'}{8\pi r^2}+P\,, \label{varrhogeneral} \\
    p_2&=p_3=\frac{rf''+2f'}{16\pi r}-P\,. \label{transpressuregeneral}
\end{align}
Here, $p_1$ is the longitudinal pressure, whereas $p_2$ and $p_3$ are the transversal pressures that will be denoted by $p_{tr}$. The term $P$ in Eqs.~\eqref{varrhogeneral} and \eqref{transpressuregeneral} is provided by the cosmological constant  in  Einstein's field equations \eqref{eq:EFE}. 
It is worth  remarking that $\rho$ is the density of the Anton-Schmidt fluid, while $\varrho$ is the density arising in the spacetime whose thermodynamics is emulated by such a fluid.
After simple manipulations, Eqs.~\eqref{varrhogeneral} and \eqref{transpressuregeneral} read
\begin{align}
\varrho&=P-\frac{a(\rho)^{1/\beta}}{4\pi \beta}r^{\frac{2}{\beta}-3}+\frac{1}{8\pi r^2}\,, \\
p_{tr}&=\left(\frac{1}{\beta}-\frac{1}{2}\right)\frac{a(\rho)^{1/\beta}}{4\pi \beta}r^{\frac{2}{\beta}-3}-P\,,
\end{align}

We can now restrict the possible solutions of Eq.~\eqref{generalmetric}, investigating when the energy conditions are satisfied. Specifically, we have:
\begin{itemize}
    \item the \textit{null energy condition} (NEC): $\varrho+p_i\ge0$ for $i=1,2,3$. In our case $p_2=p_3$ and $p_1=-\varrho$. Thus, the NEC becomes $\varrho+p_{tr}\ge 0$;
    \item the \textit{weak energy condition} (WEC): $\varrho\ge0$. 
    \item the \textit{strong energy condition} (SEC): $\varrho+\sum_{i}p_i\ge0$, which, in our case, becomes $p_{tr}\ge0$;
    \item the \textit{dominant energy condition} (DEC): $\rho\ge|p_i|$ that becomes $\rho\ge|p_{tr}|$.
\end{itemize}We can prove that there is only one metric function $f(r,\rho)$ satisfying all of them. In particular, the WEC reads
    \begin{equation}\label{WEC}
    P-\frac{a(\rho)^{1/\beta}}{4\pi \beta}r^{\frac{2}{\beta}-3}+\frac{1}{8\pi r^2}\ge0\,, 
\end{equation}
while the SEC becomes
\begin{equation}\label{SEC}
    \left(\frac{1}{\beta}-\frac{1}{2}\right)\frac{a(\rho)^{1/\beta}}{4\pi \beta}r^{\frac{2}{\beta}-3}-P\ge0\,.
\end{equation}
Moreover, the NEC condition becomes 
\begin{equation}\label{NEC}
    \varrho+p_{tr}=\frac{1}{8\pi r^2}+\left(\frac{1}{\beta}-\frac{3}{2}\right)\frac{a(\rho)^{1/\beta}}{4\pi \beta}r^{\frac{2}{\beta}-3}\ge0\,.
\end{equation}
The inequality for the DEC depends on the sign of the longitudinal and transverse pressures. In our case, if the SEC is satisfied, the DEC becomes
\begin{equation}
    2P+\frac{1}{8\pi r^2}-\left(\frac{1}{\beta}+\frac{1}{2}\right)\frac{a(\rho)^{1/\beta}}{4\pi \beta}r^{\frac{2}{\beta}-3}\ge0\,.
\end{equation}


Let us check when the WEC and the SEC are satisfied considering four different ranges of values for $\beta$.
\begin{enumerate}[i.]
\item $0<\beta<2/3$. In this case, the term $\frac{a(\rho)^{1/\beta}}{4\pi \beta}r^{\frac{2}{\beta}-3}$ in the WEC dominates when $r$ is high, since $\frac{2}{\beta}-3>2$. However, this term is also negative, therefore the WEC is not satisfied for high radii.
\item $2/3<\beta\le2$. One can easily prove that, in this case, the WEC is satisfied when
\begin{equation}\label{WECinequality}
    P(\rho)r^2-\frac{a(\rho)^{1/\beta}}{4\pi \beta}r^{\frac{2}{\beta}-1}+\frac{1}{8\pi}\ge0\,.
\end{equation}
Since $0\le\frac{2}{\beta}-1<2$, the WEC is satisfied both in the limits $r\rightarrow0$ and $r\rightarrow\infty$. However, it may happen that the WEC is not satisfied in a certain finite range of $r$. Nevertheless, ignoring this possibility, the SEC is certainly not satisfied for high radii, since the dominant term becomes $-P$.
\item $\beta>2$. For small radii, we have a negative term in the WEC dominating over the positive ones. Thus, the WEC is not satisfied for small radii.
\item $\beta=2/3$. The WEC and the SEC read, respectively, 
 \begin{align}
     P-\frac{3}{8\pi}a(\rho)^{3/2}+\frac{1}{8\pi r^2}\ge0\,, \\
     \frac{3}{8\pi}a(\rho)^{3/2}-P\ge0\,.
 \end{align}
 Since they must be satisfied for all the radii, we conclude that  
 \begin{equation}\label{energyconditiontot}
     P(\rho)=\frac{3}{8\pi}a(\rho)^{3/2}.
 \end{equation}
  Thus, $\varrho=-p_1=\frac{1}{8\pi r^2}$ and $p_{tr}=0$, which implies that also the NEC and the DEC are satisfied.
 Moreover, we note that the asymptotic AdS spacetime is exactly recovered in this case. Indeed, inserting Eq.~\eqref{energyconditiontot} into Eq.~\eqref{generalmetric}, we obtain
\begin{equation}
    f(r,\rho)=-\frac{2M}{r}+\frac{8}{3}\pi P(\rho)r^2\,.
\end{equation}
\end{enumerate}
 
 Now, we shall study the restrictions on the parameters of the Anton-Schmidt fluid emulating the situation $\beta=2/3$ just described. From Eq.~\eqref{energyconditiontot}, since $a(\rho)$ is positive, it is clear that the pressure must be positive. By applying the logarithm on both sides of Eq.~\eqref{energyconditiontot}, we obtain
\begin{equation}\label{u7}
    \frac{3}{2}\int\frac{ P_\rho}{P+\rho}d\rho=\ln\left(\frac{8\pi}{3}P\right),
 \end{equation}
where we have used the definition in  Eq.~\eqref{eq:a}. 
Then, it is possible to solve Eq.~\eqref{u7} by simply computing the derivative with respect to $\rho$. Thus, we immediately find
\begin{itemize}
    \item[--] a physical solution, namely $P=const$,
    \item[--] a unphysical solution, namely $P=2\rho$.
\end{itemize}
Even though the first case appears appealing, it just represents the widely-studied trivial cosmological constant case. It appears clear that, in order to have it from Eq.~\eqref{ASfluid}, one needs $\rho=const$, having both the pressure and density to be perfectly constant. On the other hand, the second case  would imply a sound speed faster than light. 
Consequently, in the case of the Anton-Schmidt and logotropic models, an interesting scenario would arise from the violation of at least one energy condition. 

In view of the aforementioned considerations, in what follows we analyze the case that corresponds to relaxing the SEC. 


\subsection{Relaxing the strong energy condition}


As previously stated, it appears interesting to relax the SEC. Previously, we computed  possible values of $h(r,\rho)$ in Eq.~\eqref{addtermfinal}, in order for the WEC and the SEC to hold. In so doing, we proved that two solutions for $P$ satisfy both of them, implying automatically that the NEC and the DEC hold as well. 
Since the only plausible outcome is the one with constant pressure, it is of utmost importance to investigate a wider range of possible solutions for Eq.~\eqref{generalmetric}.

We limit ourselves to those  ranges where  WEC is satisfied\footnote{As SEC is not satisfied in this ranges, $p_{tr}<0$ and consequently, the DEC coincides with the NEC, Eq.~\eqref{NEC}. So, the NEC is necessary and sufficient to prove the DEC.}, having $2/3<\beta\le2$.
Then, the NEC (as the DEC) becomes
\begin{equation}\label{DECineq}
    \left(\frac{1}{\beta}-\frac{3}{2}\right)\frac{a(\rho)^{1/\beta}}{4\pi \beta}r^{\frac{2}{\beta}-3}+\frac{1}{8\pi r^2}\ge0\,.
\end{equation}
The above condition is never satisfied when $2/3<\beta<2$, since the first term is negative in this range and it dominates over the second. The only possibilities are then $\beta=2/3$ or $\beta=2$. 
In the first case, we have exactly the same situation studied before satisfying the SEC, i.e. Eq.~\eqref{energyconditiontot}. 
Instead, in the case $\beta=2$, we find that the WEC, the NEC and the SEC are all satisfied when  $\frac{a(\rho)^{1/\beta}}{4\pi \beta}\le\frac{1}{8\pi}$, implying
\begin{equation}\label{energycondpartial}
   \exp\left(\frac{1}{2}\int\frac{P_\rho }{\rho+P}\, d\rho\right)\le1\,.
\end{equation}
Applying the logarithm to both sides, we obtain
\begin{equation}
    \int\frac{P_\rho}{\rho+P}\,d\rho\le0 \,.
\end{equation}
which is satisfied as long as the integral is upper-bounded. 
In the case of the Anton-Schmidt and pure logotropic fluids, we respectively have
\begin{equation}\label{derivativesign}
   \frac{P_\rho}{P+\rho}=\frac{\frac{A}{\rho_\ast}\left(\frac{\rho}{\rho_\ast}\right)^{-n-1}\left[1-n\ln\left(\frac{\rho}{\rho_\ast}\right)\right]}{\rho+A\left(\frac{\rho}{\rho_\ast}\right)^{-n}\ln\left(\frac{\rho}{\rho_\ast}\right)} \,,
\end{equation}
and
\begin{equation}\label{derivativesignBIS}
   \frac{P_\rho}{P+\rho}=\frac{A}{\rho\left[\rho+A\ln\left(\frac{\rho}{\rho_\ast}\right)\right]} \,.
\end{equation}



The condition $P=\frac{|\Lambda|}{8\pi}$ ensures the positivity of $P(\rho)$ and, thus, $\rho >\rho_\ast$. The sign of the right side of Eq.~\eqref{derivativesign} depends on the value of $n$:
\begin{itemize}
    \item[--] If $n>0$, the numerator is positive for $\rho<\rho_\ast e^{1/n}$ and negative for $\rho>\rho_\ast e^{1/n}$. As a consequence, the function $\ln a(\rho)$ reaches its maximum at $\rho=\rho_\ast^{1/n}$, being upper-bounded, so that and condition \eqref{energycondpartial} is  satisfied.
    
    \item[--] If $-1<n\le0$, Eq.~\eqref{derivativesign} and Eq.~\eqref{derivativesignBIS} are both always positive, and a maximum for $\ln a(\rho)$ cannot be found. However, it is easy to show that, when $\rho\gg\rho_\ast$, $\ln a(\rho)$ behaves as $\alpha-\frac{1}{\rho^{1+n}}$, where $\alpha$ is an integration constant. Since $n>-1$, then $\ln a(\rho)$ has an asymptote identified with $\alpha$, and thus it is upper-bounded also in this case.
    
    \item[--] If $n<-1$, one can easily prove that Eq.~\eqref{derivativesign} is always positive. However, in this case, the second term of the denominator dominates over the first and the whole expression behaves as $\rho^{-1}$ for large $\rho$. This means that $\ln a(\rho)$ increases always as $\ln(\rho)$,  so it is not upper-bounded and the condition \eqref{energycondpartial} is never satisfied.
\end{itemize}



\section{Optical properties of logotropic black holes}

In order to study optical properties of the above BH solution, we first focus on the validity of Eq.~\eqref{energycondpartial}, holding for  $\rho>\rho_\ast$, which requires $n>-1$ and thus $\sqrt{a(\rho)}\le1$. Since $n>-1$ ensures that $a(\rho)$ is upper-bounded, we can always find an integration constant for $a(\rho)$ such that the WEC is satisfied. 
Hence, we select this integration constant such that $a(\rho)=1$ at its maximum value. This occurs when $\rho\to\infty$, if $n<0$, and when $\rho=\rho_\ast e^{1/n}$, if $n>0$ (i.e. when $P_\rho=0$). As the first case is unphysical, since it would imply an infinite value of the pressure, we limit our prescription to $n>0$, selecting the function $a(\rho)$ such that $a(\rho_\ast e^{1/n})=1$.
In this way, we ensure that $\sqrt{a(\rho)}\le1$ is satisfied for each $\rho\ne\rho_\ast$ and that $\sqrt{a(\rho)}=1$ when $P$ reaches its maximum.

Therefore, by virtue of Eq.~\eqref{generalmetric}, one easily finds
\begin{equation}
    f(r,\rho)=\sqrt{a(\rho)}-\frac{2M}{r}\,,
    \label{schwsolution}
\end{equation}
implying that the metric \eqref{eq:metric} reads
\begin{equation}\label{pseudoschwmetric}
\begin{split}
    ds^2=\ &\sqrt{a(\rho)}\left(1-\frac{2M}{\sqrt{a(\rho)}r}\right)dt^2\\
    &-\frac{dr^2}{\sqrt{a(\rho)}\left(1-\frac{2M}{\sqrt{a(\rho)}r}\right)}-r^2d\Omega^2\,.
\end{split}
\end{equation}
Adopting the following coordinate transformations
\begin{subequations}
\begin{align}
t&\to (\sqrt{a(\rho)})^{-{1\over2}}t\,,\label{timerescaling}\\
r&\to (\sqrt{a(\rho)})^{{1\over2}}r\,,\label{rrescaling}
\end{align}
\end{subequations}
and preserving the angles, Eq.~\eqref{pseudoschwmetric} describes a Schwarzschild BH with a mass rescaled by
\begin{equation}
    M\rightarrow M(\sqrt{a(\rho)})^{-3/2}\,,
\end{equation}

The stress-energy tensor arising from this metric is 
\begin{equation}\label{pseudoschw stressenergy}
    T^{\mu\nu} =\left(P+\frac{1-\sqrt{a(\rho)}}{8\pi r^2}\right)(e_0^\mu e_0^\nu-e_1^\mu e_1^\nu)-P(e_2^\mu e_2^\nu+e_3^\mu e_3^\nu)\,,
\end{equation}
where the standard Schwarzschild case is recovered as $\sqrt{a(\rho)}\rightarrow1$ and the stress-energy tensor reduces to 
\begin{equation}
    T^{\mu\nu}=\frac{|\Lambda|}{8\pi}\left(e_0^\mu e_0^\nu-\sum_{i=1}^3e_i^\mu e_i^\nu\right)=\frac{|\Lambda|}{8\pi}g^{\mu\nu}\,.
\end{equation}
The latter represents a source for Einstein's field equations that cancels with the term $|\Lambda|g^{\mu\nu}$, thus leading to a vacuum solution and allowing to recover completely the Schwarzschild solution.
In this respect, it appears evident that an Anton-Schmidt fluid is thermodynamically equivalent to a Schwarzschild BH when its density is $\rho=\rho_\ast e^{1/n}$. As the fluid pressure takes its maximum value $P=\frac{A}{n e}$, $P$ being associated to the cosmological pressure, we can fix $A$ as $A=\frac{n e|\Lambda|}{8\pi}$, for $n>0$.



\subsection{The effective refractive index}

We here investigate the properties of our solution arising from the condition $\sqrt{a(\rho)}\ne1$, which implies modifications of the Schwarzschild metric. In particular, we relate this effect to optical properties of the spacetime in presence of a  medium made by logotropic fluids.

To do so, we study the refractive index that could lead to different results due its double interpretation, namely the \textit{optical refractive index} \cite{Perlick:2010zh,Gibbons:2008rj}, $n_o$, and the \textit{Fermat refractive index} \cite{Yi_2011}, $n_F$.
Specifically, we check whether the effects of a logotropic fluid medium are significant to change the optical configuration around a BH. We also show that our procedure is general and can be adapted to other thermodynamic models. 

\subsubsection{Optical refractive index}

The effective optical refractive index, $n_o$, emerges by modifying the metric as \cite{Gibbons:2008rj}
\begin{equation}\label{opticalmetric}
   ds^2= f^2(\hat{r})\left[dt^2-n_o^2\left(d\hat{r}^2+\hat{r}^2d\Omega^2 \right)\right]\,,
\end{equation}
and we can  find plausible transformations to match Eq.~\eqref{pseudoschwmetric} with Eq.~\eqref{opticalmetric}. For the sake of simplicity, we only consider the time rescaling \eqref{timerescaling} neglecting the rescaling of the radius \eqref{rrescaling}. One can easily prove that the refraction indices are the same up to a mass rescaling $M\to M(\sqrt{a(\rho)})^{-1/2}$. Equating the components of the two metrics, we thus obtain
\begin{align}
    f(\hat{r})&=\left(1-\frac{2M}{\sqrt{a(\rho)}r}\right)^{1/2}\,, \label{condtime} 
    \\
    f(\hat{r})n_o(\hat{r})d\hat{r}&=\frac{dr}{\sqrt{\sqrt{a(\rho)}\left(1-\frac{2M}{\sqrt{a(\rho)}r}\right)}}\,, \label{condlong} 
    \\
    f(\hat{r})n_o(\hat{r})\hat{r}&=r\,. \label{condangular}
\end{align}
Dividing Eq.~\eqref{condlong} by Eq.~\eqref{condangular}, we get the following differential equation:
\begin{equation}
    \frac{d\hat{r}}{\hat{r}}=\frac{dr}{r{\sqrt{\sqrt{a}\left(1-\frac{2M}{\sqrt{a}r}\right)}}}\,,
\end{equation}
which, once integrated, gives
\begin{equation}\label{opt solution}
    K\hat{r}^{\sqrt[4]{a}}=r\frac{\sqrt{a}}{M}-1+\sqrt{r\frac{\sqrt{a}}{M}}\sqrt{r\frac{\sqrt{a}}{M}-2}\,,
\end{equation}
where $K$ is an integration constant. Then, inverting Eq.~\eqref{opt solution} yields
\begin{equation}
    r(\hat{r})=\frac{M}{\sqrt{a}}\frac{\left(1+K\hat{r}^{\sqrt[4]{a}}\right)^2}{2K\hat{r}^{\sqrt[4]{a}}}\,.
\end{equation}
so that, from Eq.~\eqref{condtime}, we have
\begin{equation}
    f(\hat{r})=\frac{K\hat{r}^{\sqrt[4]{a}}-1}{K\hat{r}^{\sqrt[4]{a}}+1}\,.
\end{equation}
Hence, the optical refractive index could be computed from Eq.~\eqref{condangular} as
\begin{equation}
    n_o(\hat{r})=\frac{M}{\sqrt{a}}\frac{\left(1+K\hat{r}^{\sqrt[4]{a}}\right)^3}{2K\hat{r}^{\sqrt[4]{a}+1}\left(K\hat{r}^{\sqrt[4]{a}}-1\right)}\,.
\end{equation}
The Schwarzschild case is recovered when $\sqrt{a}=1$. To be consistent with  the notation of \cite{Perlick:2010zh}, in the Schwarzschild case, we set $K=\frac{2}{M}$, thus leading to
\begin{equation}\label{opticalindex}
    n_o(\hat{r})=\frac{1}{\sqrt{a}\hat{r}^{1-\sqrt[4]{a}}}\frac{\left(1+\frac{M}{2\hat{r}^{\sqrt[4]{a}}}\right)^3}{1-\frac{M}{2\hat{r}^{\sqrt[4]{a}}}}.
\end{equation}

Two cases of interest occur for small and large radii. In particular, for small radii, the index of refraction diverges at the horizon radius \cite{Perlick:2010zh,Yi_2011}, i.e. $r_h=\frac{2M}{\sqrt{a}}$, as
\begin{equation}
    n_o(r\to r_h)\sim \frac{2^{\frac{1}{2}+a^{-1/4}}}{M^{a^{-1/4}}}\frac{r_h}{\sqrt{r-r_h}}\,.
\end{equation}

Instead, the case of large radii is not straightforward. In particular,  for $\sqrt{a}\ne1$, the optical refractive index seems to tend to zero, rather than to the unity:
\begin{equation}
    n_o(\hat{r}\to \infty)\sim\frac{1}{\sqrt{a}\hat{r}^{1-\sqrt[4]{a}}}.
\end{equation}
\begin{figure}
    \centering
    \includegraphics[scale=0.67]{OpticalIndex.pdf}
    \caption{Effective optical refractive index as a function of the external radial coordinate $r$ (distance from the BH center) in a logarithmic scale. The curves correspond to different values of $\sqrt{a(\rho)}$: $\sqrt{a}=1$ (solid), corresponding to the Schwarzschild case); $\sqrt{a}=0.95$ (dashed); $\sqrt{a}=0.9$ (dot-dashed); $\sqrt{a}=0.85$ (dotted). Units of $2M=1$ were adopted.}
    \label{fig: optical index}
\end{figure}
Such a feature is shown in Fig.~\ref{fig: optical index}, where it appears evident that $n_o$ seems to break down at a given $r$, as it does not tend to unity like in the Schwarzschild case. The physical explanation for this behavior is related to the definition of $n_0$. Indeed, at very large radii, the density $\rho$ cannot act as a medium, since the corresponding BH would act as a fully transparent point-like object whose index of refraction cannot depart from $n_o=1$. This can be seen if one does not fix $\sqrt{a(\rho)}$ to a given value, assuming that $\rho$ distributes providing a limiting case for the  Minkowski spacetime. To see that, it is straightforward to notice that, for varying $\rho$, when $\hat r\rightarrow\infty$, $a(\rho)\rightarrow0$ in the denominator of Eq.~\eqref{pseudoschwmetric}, since there is no evidence of logotropic fluids at very large distances. To guarantee asymptotic flatness of Eq.~\eqref{pseudoschwmetric}, one requires $\sqrt{a(\rho)}r\rightarrow\infty $ as both $r\rightarrow\infty$ and $P\rightarrow0$. This can be generalized to any fluid whose density does not explicitly depends on the radial coordinate. Indeed, if one considers $\rho=\rho(\hat r)$, it would be possible to explore the above case without the need to have
$\sqrt{a(\rho)}r\rightarrow\infty $ when both $r\rightarrow\infty$ and $P\rightarrow0$.

The above considerations about the optical refractive index imply that its validity still holds, but it highly depends on the  functional evolution of $\rho$ in terms of $\hat r$. To overcome this issue, one can search for an optical index that, on the contrary, does not take into account how the EoS evolves as a function of the radial coordinate. This subject is investigated in the following, where we deal with effective refractive index resulting from the application of  the Fermat principle.


\subsubsection{Refractive index from the Fermat principle}

As stated above, a simpler approach to investigating the refractive index involves the use of the Fermat principle. This approach turns out to be quite different than $n_0$, because it does not depend on the radial distance $\hat r$.
Its use spans within several gravitational lensing contexts  \cite{Liu:2015wma,Walters:2010gk,Yi_2011}, and it is constructed by considering light rays, i.e. $ds^2=0$, with constant angles. Thus, from Eq.~\eqref{pseudoschwmetric}, we have
\begin{equation}
    dt=\frac{dl}{\sqrt{1-\frac{2M}{\sqrt{a}r}}}\,,
\end{equation}
where $dl^2$ is the spatial part of Eq.~\eqref{pseudoschwmetric}. For a path $\Gamma$, the Fermat principle  reads
\begin{equation}
    \delta\int\frac{d\Gamma}{\sqrt{1-\frac{2M}{\sqrt{a(\rho)}r}}}=0\,.
\end{equation}
Hence, the corresponding refractive index is obtained by comparing the latter with the Fermat principle in the flat case $\delta\int n_F \,d\Gamma=0$:
\begin{equation}\label{Fermat refindex}
    n_F(r)=\left(1-\frac{2M}{\sqrt{a(\rho)}r}\right)^{-{1\over2}}.
\end{equation}
\begin{figure}
    \centering
    \includegraphics[scale=0.67]{refIndex.pdf}
    \caption{Fermat refractive index as a function of the distance from the BH center $r$ in a logarithmic scale. The different curves correspond to the values of $\sqrt{a(\rho)}$ as in Fig.~\ref{fig: optical index}.}
    \label{refindex}
\end{figure}
\noindent The behavior of the latter, for different values of $\sqrt{a}$, is shown in Fig.~\ref{refindex}. We notice that, for large radii, the Fermat refractive index tends to  $n_F=1$, regardless of the value $\sqrt{a(\rho)}$, consistently  with the fact that, far from the BH, the influence of the fluid cannot be perceived. For small radii, $n_F$ diverges at the event horizon $r_h$ as
\begin{equation}
    n_F(r\to r_h)\sim\sqrt{\frac{r_h}{r-r_h}}.
\end{equation}
\subsubsection{Horizon radius shift}
From both the optical and Fermat refractive indices, it appears evident that the BH event horizon is shifted when $\sqrt{a(\rho)}\ne1$ from $2M$ to $\frac{2M}{\sqrt{a(\rho)}}$. The deviation of the horizon radius from the Schwarzschild case can be expressed as a function of the Anton-Schmidt density as
\begin{equation}\label{horizonradius}
    \Delta r_h(\rho)=\frac{r_h(\rho)-2M}{2M}=\frac{1}{\sqrt{a(\rho)}}-1\,.
\end{equation}
\begin{figure}
    \centering
    \includegraphics[scale=0.67]{MassForDensitynorotatelabel.jpg}
    \caption{Deviation of the horizon radius as a function of the density of Anton-Schmidt's fluid. The values on the $x$ and $y$ axes are expressed in powers of $10^{4}$ and $10^{-6}$, respectively. The parameter $n$ is fixed to $0.1$, while the different curves correspond to different values of $\rho_\ast$: $\rho_\ast=0.5 A$ (solid), $\rho_\ast=A$ (dotted) and $\rho_\ast=2A$ (dashed). The Schwarzschild solution is recovered at the minima of the curves, i.e. when $\rho/\rho_\star=e^{1/n}\simeq 22\cdot10^3.$ The top-right subplot enhances the behavior of $r_h-1$ around such minimum.}
    \label{massscaling}
\end{figure}
A plot of the latter quantity is provided in Fig.~\ref{massscaling}. It is interesting to notice how different values of the Anton-Schmidt parameters (namely, different ratios $\rho_\ast/A$) affect the behavior of the horizon radius in the optical framework. Namely, the greater $\rho_\ast/A$, the smaller deviations of the horizon radius from the Schwarzschild one occur as $\rho$ deviates from $\rho_\ast e^{1/n}$.




\section{Outlook and perspectives}
\label{sec:conclusions}

In this paper, we investigated a class of asymptotic AdS BH metrics whose thermodynamics matches that of logotropic models. For this purpose, we considered Anton-Schmidt's fluid, which contains the pure logotropic case as a limiting case. We thus developed a general method to obtain a metric solution and, by requiring the presence of a BH horizon, we found plausible metric functions describing the physical scenario under study. We then studied the most suitable values of the free parameters associated with our class of solutions. 

In particular, taking into account the  energy conditions for the source fluid, we showed that a BH solution of the  Einstein field equations is obtained for an asymptotic AdS spacetime, leading to two specific cases: i) a pure constant pressure, corresponding to a constant density of the Anton-Schmidt fluid; ii) an exotic outcome, which should be physically discarded. In this respect, we checked whether, relaxing the SEC, more general classes of solutions may be found. In particular, our analysis demonstrates that an Anton-Schmidt fluid with positive pressure is capable of emulating the thermodynamics of a spherically-symmetric compact object as soon as the pressure of the fluid is maximized. 

Furthermore, we analyzed the optical properties of our solutions. To do so, we followed the standard recipe to obtain an effective refractive index in spherical coordinates. We thus adopted two definitions, namely the optical and the Fermat refractive indexes. In the first case, we showed the limits of the model predictions, which mostly require knowing how the density evolves as a function of the radial distance, and provides an unphysical behavior for large radii. Moreover, to overcome the aforementioned issues, we considered the Fermat refractive index. 
In so doing, we discussed the functional dependence on $r$, with particular regard to its asymptotic values. Likely divergences at small radii were discussed, whereas a direct comparison with the Schwarzschild case was prompted, showing where the deviations from the latter case are much more evident. Quite clearly, the corresponding optical effects imply that the refractive indexes increase due to the presence of the thermodynamic medium constituted by Anton-Schmidt's fluids. 

Future works will focus on alternative fluids characterized by thermodynamic effects that are mainly different from logotropic models. In particular, it would be interesting to take into account real fluids and show, for instance, how transition phases could affect these scenarios. Another plausible extension would be to promote our treatment to more general compact objects \cite{Cardoso:2019rvt} and/or on accretion disk contexts \cite{DAgostino:2022ckg,Kurmanov:2021uqv,Boshkayev:2021chc}. A crucial point deserving future efforts will be investigating the relation between the thermodynamic properties of the fluid and the BH, in order to enable a richer use of thermodynamics and make the extended phase approach more robust and fully justified from a physical perspective. 



\begin{acknowledgements}
S.C. and R.D. acknowledge the support of Istituto Nazionale di Fisica Nucleare (INFN), \textit{iniziative specifiche} MOONLIGHT2 and QGSKY. A.L. and O.L. would like to thank Roberto Giambò and Stefano Mancini for helpful discussions. O.L. acknowledges the Ministry of Education and Science of the Republic of Kazakhstan, Grant: IRN AP08052311. 
\end{acknowledgements}


