\section{Introduction}

The synthesis and characterization of molecule-based magnetic systems has been an intense research area for decades,
prompted by the need for information storage devices and advances in quantum technologies (see for example Refs.~\cite{gaita2019molecular,troiani2011molecular,stamp2009spin,mcadams2017molecular}).
The motivations for targeting such complexes stem from their physical-chemical properties
ranging from long coherence time~\cite{atzori2016quantum, atzori2016room, bader2014room, graham2014influence, atzori2018structural}
to manipulation possibilities~\cite{bayliss2020optically,carretta2021perspective,li2021manipulation,nelson2020cnot,thiele2014electrically,
hussain2018coherent}.
In the presence of spin-switchable units, the local spin state of archetypical
 iron(II) or cobalt(II) ($3d^6$ and $3d^7$, respectively) ions can be controlled using
 external fields. Important structural modifications are then observed, with changes in bond distances of up to 0.2 \AA~(\textit{i.e.} 10 \%)\cite{hauser2004ligand}
 and deep changes in  charge distribution (up to 0.5 electrons)~\cite{kepenekian2009primary}. The latter are responsible for
hysteretic behavior in materials,\cite{kenepekian2009magnetic}
and have long been underestimated whereas they are the main characteristics of valence tautomers\cite{tezgerevska2014valence}.
Let us stress that the transition between tautomers is not necessarily
accompanied by a modification of the spin multiplicity, in contrast with spin-transition systems.
More recently, transition metal ions combined with organic ligands have been considered as possible
targets in the development of molecular-based quantum units of information, \textit{e.g.} qubits or qudits (see for example Refs.~\cite{gaita2019molecular,troiani2011molecular,stamp2009spin,mcadams2017molecular}). Theoretical studies have also revealed the potential interest of radical ligands for the manipulation of quantum information via the entanglement of local spin degrees of freedom~\cite{roseiro2022combining,roseiro2022Qubit}.
A prerequisite is the binding of radical ligands to paramagnetic metal centers
without losing their open-shell character.
In this respect, oxoverdazyl-based ligands have proven to fulfill such requirement,
giving rise to a wealth of magnetic coupling schemes\cite{oms2010beyond,brook2010strong,barclay2003verdazyl}.
Furthermore, the flexibility and well-established redox activity of organic-based compounds make
such materials particularly interesting.
Not only can inter-unit interactions be modulated with speculated
spin-crossover behaviour\cite{rota2010toward},
but the crystal field generated by several open-shell ligands may give rise
to unusual and puzzling scenarios.

Recently, it has been suggested that the spectroscopy of a cobalt(II) ion coordinated to
open-shell radicals may not be readily interpreted.
At the crossroads of exchange coupled and spin-crossover systems, this
compound has questioned the traditional pictures emerging from a metal ion, either high-spin or low-spin, in the electrostatic field of neighbouring ligands~\cite{fleming2020valence,roseiro2022combining}.
Wavefunction-based calculations supported that the spin states are characterized
by combinations of various local spin
states on the metal center, a feature of entanglement
named \textit{spinmerism}.
 \begin{figure
    \centering
    \includegraphics[width=8cm]{ROSEIRO_FIG1.pdf}
    \caption{ {[Fe(dipyvd)$_2$]$^{2+}$ compound \textbf{(1)} from Ref\cite{brook2018anelectron}.}}
    \label{fig:molecule_intro}
\end{figure}
This phenomenon can be seen as an extension of mesomerism to spin degrees of freedom with entanglement
in between two local sub-parts of a molecule.
Since the ligand field includes Coulomb and direct exchange contributions in a complex built on spin-coupled partners, one may wonder whether different local spin states may coexist on the metal ion.
The introduction of radical ligands may indeed disrupt the assumption of a given spin state on the metal centre
in the description of ground and excited states.

Inspired by these observations, the pseudo-octahedral iron(II)-oxoverdazyl
compound [Fe(dipyvd)$_2$]$^{2+}$, with dipyvd = 1-
isopropyl-3,5-dipyridil-6-oxoverdazyl (\textbf{1}) was considered as another prototype
to be examined (see Figure~\ref{fig:molecule_intro})\cite{brook2018anelectron}.
Magnetic and spectroscopic properties were reported in the literature
and complemented by density functional theory-based calculations.
Iron(II) being a spin-crossover metal ion, any intermediate crystal field
should be appropriate to observe the coexistence of a $S_{Fe} = 0$ (low spin)
and $S_{Fe} = 2$ (high spin). From the Fe-N bond distances values (\textit{ca.} 1.9 \AA),
a low-spin iron center is expected and was indeed confirmed\cite{brook2018anelectron}.
However, the $d^6$ Tanabe-Sugano diagram exhibits
a crossing between excited triplet and quintet states (see Figure~\ref{fig:d6_TS_diagram}).
Therefore,
\textit{excited sate spinmerism} might be anticipated, with total spin states
resulting from combinations of local $S_{Fe} = 1$ and $S_{Fe} = 2$.
For all these reasons, \textit{ab initio} calculations were carried out to inspect the
spectroscopy of compound \textbf{1}  (shown in Figure~\ref{fig:molecule_intro}). The eigen-states were constructed using
localized molecular orbitals (LMOs) to allow for projections onto the local spin states
 of the metal and ligands, $S_{Fe}$ and  $S_{L}$, respectively.

The numerical results presented in this work highlight the presence of a strong quantum entanglement between the local spin states, \textit{i.e.} a \textit{spinmerism} effect,
of the metal ion and the radical ligand environment.
This evidence reshuffles the standard views in molecular magnetism and simultaneously opens up new perspectives for technological development at the interface of physics and chemistry. The variability of local spin states could provide a pathway to encode quantum information on synthetic molecular systems, paving the way for the generation of molecular spin-qubits.
 \begin{figure
    \centering
    \includegraphics[width=\columnwidth]{fig2_tanabesuganod6.png}
    \caption{ {Low part of the Tanabe-Sugano diagram of an iron(II) d$^6$ ion in an octahedral field. The critical ligand field
    value marking spin-crossover is
    indicated by a vertical dotted line. The crossing between the excited $S_{Fe} = 1$ and
    $S_{Fe} = 2$  states is highlighted.}}
    \label{fig:d6_TS_diagram}
\end{figure}

\section{Computational details}

The presence of multiple open-shells on both metal ion and ligands strongly invites the use of a wavefunction-based method,
such as the complete active space self-consistent field (CASSCF) method. The complete active space should include six electrons from the iron center and two from the oxoverdazyl ligands in seven molecular orbitals (CAS[8,7]).
However, the active space was reduced to CAS[6,6] by inspecting the occupation numbers of the active molecular orbitals (MOs).
The CAS[6,6]SCF sextuplet MOs were then transformed into localized molecular orbitals (LMOs)\cite{DOLOcode} which are shown in Figure~\ref{fig:MOs_cas6_6}. Such transformation allows for a reading of the wavefunctions following a Lewis-like description.
In the low-energy spectrum of \textbf{1}, we checked that the mostly metal d-type LMOs remained doubly occupied (low-spin iron(II)) and
the active space was even reduced to CAS[2,2].
All CASSCF calculations were performed on the available Xray structure without any geometry optimization and used the MolCAS 8.0 package~\cite{MOLCAS}.
The iron and first coordination sphere atoms were described with 4s3p2d and 3s2p1d basis sets, respectively.
Smaller basis sets 3s2p were used for all other atoms whilst hydrogen atoms were depicted with a 1s basis set.
The dynamical correlation and polarization effects were included following the Difference Dedicated Configuration Interaction (DDCI) method~\cite{DDCIref1,DDCIref2} as implemented in the CASDI code\cite{CASDIcode}. Given a set of MOs, the structures of the spin states can be directly compared from the CI amplitudes.
Depending on the classes of excitations involved in the CI expansion, different levels (CAS + S, CAS + DDC2 and CAS + DDCI) can be reached beyond the CAS pictures (CAS[2,2] or CAS[6,6]).
This variational method that follows a step-by-step construction of the wavefunction has produced a wealth of interpretations and rationalizations
with a systematic relaxation of the wavefunction (so-called "fully decontracted method").
The resulting CI
eigenfunctions were projected onto the local spin states of the oxoverdazyl ligands pair ($S_{L} = 0$ or $1$)
and the iron ion ($S_{Fe} = 0$, $1$ or $2$). This procedure allows for a decomposition into the different entangled metal-ligand contributions.
To implement these projections and conduct our local spin analysis, the open access package \textit{QuantNBody}~\cite{codeQuantNBody} was used.
This numerical python toolbox has been recently developed by one of us (SY) to facilitate the manipulation
of quantum many-body operators and wavefunctions.
Based on this tool, matrix representations (in the many-electron basis) of the local metal and ligands spin operators $\hat{S}^2_{Fe}$ and $\hat{S}^2_{L}$ were built and diagonalized to access the local spin subspaces. This approach makes it possible to design spin projectors (in the many-electron basis) to target specific local spins contributions for the metal and the ligands in the multi-reference wavefunction.

\section{Results and Discussion}

First, CASSCF calculations were conducted on \textbf{1}. Inspections based on a
CAS[8,7]  highlight the presence of a MO
mostly localized on a 3d iron atomic orbital with occupation number 1.99.
Thus, the active space was reduced down to [6,6]. The CAS[6,6]SCF MOs were localized either on the metal center
or on each individual dipyvd ligand (see Figure~\ref{fig:MOs_cas6_6}).
 \begin{figure
    \centering
    \includegraphics[width=\columnwidth]{fig3_MOs_CAS.jpg}
    \caption{ {CAS[6,6] LMOs of \textbf{1} generated from a CAS[6,6]SCF calculation for the sextuplet spin state.}}
    \label{fig:MOs_cas6_6}
\end{figure}
Then, the electronic structure of the environment was inspected following
the procedure recently developed on a [Co(dipyvd)$_2$]$^{2+}$ compound (cobalt(II)).~\cite{fleming2020valence,roseiro2022combining}
Based on a fictitious [Zn(dipyvd)$_2$]$^{2+}$
(closed-shell divalent metal ion) analogue of the cobalt(II) compound, it was shown that
the triplet-singlet energy difference is of the order of 2 cm$^{-1}$.~\cite{roseiro2022combining}
Considering the geometry changes moving to compound \textbf{1}, a similar inspection  was conducted by
substituting iron(II) by zinc(II).
A 1.9 cm$^{-1}$ triplet-singlet
energy difference was computed at the CAS[2,2] + DDCI level.
The absence of any quantitative exchange coupling between the radical
ligands  suggests that any projection of the total spin states of \textbf{1} may
simultaneously involve the local $S_L = 0$ and $S_L=1$ states.
Evidently, the nature of the metal center is
likely to modify the exchange coupling constant value,
a particular issue we wanted to inspect.

Due to the system size and the number of open-shells, a CAS[6,6] + DDCI level of calculation is out of reach.
Thus, all our conclusions are based on CAS[6,6] + DDC2
excitation energies and
corresponding wavefunctions analysis, as summarized in Table~\ref{table:DDCI2_results}.
\begin{table}[h!]
\centering
\begin{tabular}{|   m{1.5cm} |   m{1.8cm}  |   m{1.7cm}  |  m{1.7cm} |}
 \hline
 Energy ($cm^{-1}$)& $2S_{total}+1$ & $S_{Fe}$ & $S_{L}$ \\
 \hline\hline
  5326 & 7 & 100\% \textit{Q} & 100\% \textit{T} \\
  \hline
  3323 & 5 &     { 62\%\textit{T} \par 21\% \textit{Q} }    &   {  11\% \textit{S}  \par  72\% \textit{T} } \\
  \hline
  2906 & 5 &     { 17\% \textit{T}  \par  72\% \textit{Q}} & {  39\% \textit{S}  \par   50\% \textit{T}} \\
  \hline
  120 & 1 &     { 84\% \textit{S}} & {  84\% \textit{S}  } \\
  \hline
  0 (reference) & 3 & 87\% \textit{S} & 87\% \textit{T} \\
 \hline
\end{tabular}
\caption{\textbf{Energy spectrum of \textbf{1}. Quintets and  heptuplet energies are
calculated at the CAS[6,6] + DDC2 level.
The ground state energy is used as a reference energy. Spin multiplicities $2S_{total}+1$, local spin proportions ($S_{Fe}$ and $S_L$) are given. $S$, $T$ and $Q$ correspond to singlet, triplet and quintet, respectively.
}}
\label{table:DDCI2_results}
\end{table}
First, the ground state is  triplet  ($2S_{total} +1 = 3$),
dominated by  the local spin states  $S_{Fe} = 0$ and $S_L = 1$, respectively, followed by a singlet  lying $120$ cm$^{-1}$ above. This picture is consistent
with a $d^6$ metal ion in a high-field environment. Let us mention that the wavefunction includes 13\% of charge transfers contributions ($S_{Fe} = 1/2$ and $S_L=1/2$).
Since the CAS[6,6]SCF occupation numbers of the mostly d-type LMOs are larger than 1.98,
the singlet-triplet energy difference was further
confirmed from CAS[2,2] + DDCI calculations (CAS[2,2]SCF triplet MOs) and turned out to be $120$ cm$^{-1}$. This value is in reasonable agreement
with the reported one which may vary depending on the extraction from broken-symmetry
calculations\cite{brook2018anelectron}.
Thus, the low-energy part of the spectrum of \textbf{1} can be viewed as
two organic radicals coupled through a closed-shell iron(II) (see Figure~\ref{fig:spectrum_fevdz}).
Evidently, a Heisenberg Hamiltonian $\hat{H} = -2J \hat{s}_{1}  \hat{s}_{2}$
can be derived from the singlet-triplet energy difference. The
strong ferromagnetic exchange coupling constant $J = 60$ cm$^{-1}$  accounts for the low-temperature magnetic properties experimentally observed.
 \begin{figure
    \centering
    \includegraphics[width=\columnwidth]{fig4_spectrum_fevdz.png}
    \caption{ {Energy spectrum of \textbf{1} reflecting a Heisenberg behaviour in the low-energy part, whereas \textit{excited sate spinmerism} is observed higher in energy. The excited quintet states $Q_1$ and $Q_2$
    exhibit a strong mixing between $S_{Fe} = 1$ and $S_{Fe} = 2$ local spin states.}}
    \label{fig:spectrum_fevdz}
\end{figure}

Based on this preliminary observation suggesting a $S_L = 1$ crystal field, let us then concentrate on the third ($2 S_{total} + 1 = 5$)
and fourth ($2 S_{total} + 1 = 7$) excited states
given in Table~\ref{table:DDCI2_results}. The iron(II) spin state is successively dominated by
a $S_{Fe} = 1$ (62\%) and $S_{Fe} = 2$ (100\%).
Such picture obtained from the
spin projections of the wavefunctions is consistent with a high-field $d^6$ Tanabe-Sugano diagram~\cite{tanabe1954onthe}.
However, the quintet excited state $Q_2$ in Figure~\ref{fig:spectrum_fevdz}
 exhibits a non-negligible  contribution  from
the local  $S_{Fe} = 2$  spin state (21\%). Such superposition of
triplet and quintet metal spin states in the absence of spin-orbit coupling is
a manifestation of the recently reported \textit{spinmerism} effect\cite{roseiro2022combining}.
Our analysis supports the appearance of a similar phenomenon
in \textbf{1}, named as
\textit{excited state spinmerism}.
This manifestation results from the open-shell
character of the environment and the energy crossing between the $S_{Fe} = 1$ and $S_{Fe} = 2$ excited states in
the high-field regime of the iron(II)  $d^6$ Tanabe-Sugano diagram.
The traditional resonance in the Tanabe-Sugano diagram is lifted from the intimate open-shell nature of the
crystal field.
As expected, the spectroscopy exhibits a second close-in-energy quintet state resulting from this mixing (see $Q_1$ in Figure~\ref{fig:spectrum_fevdz}).
The latter is found 2906 cm$^{-1}$ above the ground state, with a dominant $S_{Fe} = 2$ character (72\%), and
a significant mixing between the environment $S_L = 0$ (39\%) and   $S_L = 1$ (50\%) spin state.
The entanglement between the spin states of the metal ion and its coordination sphere
is manifested in the excited state of the high-field regime of the Tanabe Sugano diagram. In agreement with a recent inspection based on a model Hamiltonian~\cite{roseiro2022Qubit}, the emergence of the two quintet  states $Q_1$ and $Q_2$ in Figure~\ref{fig:spectrum_fevdz}
from avoided crossing stresses the importance of not only the
magnitude of the crystal field (Coulomb contributions), but also the
open-shell character of such field (decisive exchange contributions).

\section{Conclusion}
The spin states structures of a coordination compound [Fe(dipvdz)$_2$]$^{2+}$ built on a spin-crossover
ion (iron(II)) and two  radical ligands (oxoverdazyl) were analyzed from DDCI wavefunctions calculations.
The procedure is based on the generation  of localized molecular
orbitals and spin projections
onto the local spin states. A ground state triplet $S_{total} = 1$ characterized by a $S_{Fe} = 0$ local spin state
is found, a reflection of a high-field regime.
A strong ferromagnetic interaction $J = 60$ cm$^{-1}$   is calculated
featuring the coupling of organic radical spin holders through a low-spin metal ion.
Even though the low-energy part of the spectrum can be rationalized by a Heisenberg spin Hamiltonian,
the magnetic picture delivered by a d$^6$ Tanabe-Sugano diagram is deeply reshuffled in
the next nearest excited states. The flexibility afforded by the open-shell character of the environment gives rise
to marked superpositions of local spin states $S_{Fe} = 1$
and $S_{Fe} = 2$  in the $S_{total} = 2$ excited states.
The traditional quasi-degeneracy in the high-field regime is lifted by 417 cm$^{-1}$
with significant contributions from both $S_{L} = 1$ and $S_{L} = 0$ on the ligand pair.
This observation, which we named \textit{excited state spinmerism}, extends a phenomenon that was reported in a
cobalt(II) analogue.
The prerequisite for the manifestation of such entanglement is fulfilled
from the presence of an iron(II) ion and radical organic ligands.
This particular class of coordination chemistry compounds combining versatile local spin states
not only enlarges the traditional pictures in molecular magnetism but
might become original targets for spin-qubit generation.

\section*{Acknowledgements}
This work was supported by the Interdisciplinary Thematic Institute SysChem via the IdEx Unistra (ANR-10-IDEX-0002) within the program Investissement d’Avenir.
P. R. acknowledges the Ecole Doctorale de
Sciences Chimiques de Strasbourg, EDSC222, and the french minister for financial support.
D.J.R.B acknowledges the support of the National Science Foundation Grant CHE-1900491.
The authors would like to thank Pr. C. Train, Dr. G. Novitchi and Dr. S. Stoian for useful discussions.

\section*{AUTHOR DECLARATIONS}

\subsection*{CONFLICTS OF INTEREST}
The authors have no conflicts to disclose.

\subsection*{ DATA AVAILABILITY }
The data that support the findings of this study are available from the corresponding author upon reasonable request.


\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}