\section{\bf Introduction}\label{intro}
Let $M$ be a connected smooth semi-Riemannian manifold of dimension $n$ $(\geq 3)$ endowed with the semi-Riemannian metric $g$  with signature $(p, n-p)$ and the Levi-Civita connection $\nabla$. If $p = 1$ or $n-1$, then $M$ is Lorentzian and if $p = 0$ or $n$, then $M$ is Riemannian. Let  $R$, $S$ and $\kappa $ be respectively the Riemann-Christoffel curvature tensor of type (0,4), the Ricci tensor of type (0,2) and the scalar curvature of $M$. As a generalization of manifold of constant curvature, Cartan \cite{Cart26} introduced the notion of locally symmetric manifolds defined as $\nabla R =0$. Again, Cartan \cite{Cart46} generalized the concept of locally symmetric manifolds and introduced the notion of semisymmetric manifolds which were latter classified by Szab\'o (\cite{Szab82}, \cite{Szab84}, \cite{Szab85}). During the study of totally umbilical submanifolds of semisymmetric manifolds, Adam\'{o}w and Deszcz \cite{AD83} introduced the notion of pseudosymmetric manifolds which was latter on studied by many authors (cf. \cite{Desz92}, \cite{Desz93}, \cite{DGHS11}, \cite{DG98}, \cite{DHV08}, \cite{HV04}, \cite{HV07a}, \cite{HV07b}, \cite{SDHJK15} and also references therein). We note that the recent trend of modern mathematical research is the abstraction, generalization or extension, existence, characterization, classification and finally applications. We mention that during the last four decades generalized concept of pseudosymmetry by Deszcz has been characterized and classified with existence in several papers by Deszcz and his coauthors (see \cite{DD91a}, \cite{DD91b}, \cite{Desz93}, \cite{DGHZ15}, \cite{SDHJK15} and also references therein); and the pseudosymmetric spacetimes were studied by Deszcz and his co-authors (see, \cite{DDVV94}, \cite{DHV04}, \cite{DVV91}). It may be noted that G\"{o}del spacetime (\cite{DHJKS14}, \cite{Gode49}), Som-Raychaudhuri spacetime (\cite{SK16srs}, \cite{SR68}), Reissner-Nordstr\"{o}m spacetime \cite{Kowa06} and Robertson-Walker spacetimes (\cite{ADEHM14}, \cite{DDHKS00}, \cite{DK99}), are models of different pseudosymmetric type structures.\\
\indent It is known that the Robinson-Trautman metrics are the spacetimes which admit a geodesic, hypersurface orthogonal, shear-free and expanding null congruences. These metrics are the general relativistic analogues of Lienard-Wiechert solutions of Maxwell equations. Moreover, the Robinson-Trautman spacetimes generalize the assumption of spherical symmetry by having topologically equivalent two-spheres rather than strictly two spheres. From the physical point of view, the Robinson-Trautman metrics can be thought of as representing an isolated gravitationally radiating system. The Robinson-Trautman metric of Petrov type II is given by (\cite{RT62}, \cite{ste03}, \cite{zun75})
$$
ds^2 = -2(U^0 - 2 \gamma^0 r - \Psi^0_2 r^{-1})du^2 + 2 du dr - \frac{r^2}{2 P^2} d\zeta d\bar\zeta,
$$
where $U^0, \gamma^0, \Psi^0_2$ are constants and $P$ is a nowhere vanishing function of $\zeta$ and $\bar\zeta$. For the sake of simplicity, we write these constants $U^0, \gamma^0, \Psi^0_2$ as $a, b, q$ respectively, $P$ as $f$ and the coordinate $u$ as $t$. The above metric can be written as
\be\label{rtm}
ds^2 = -2(a - 2 b r - q r^{-1})dt^2 + 2 dt dr - \frac{2 r^2}{f^2} d\zeta d\bar\zeta.
\ee
It is easy to check that this metric can be written as a warped product metric (\cite{BO69}, \cite{Kruc57}) given by
$$ds^2 = d\bar s^2 + r^2 d\tilde s^2,$$
where $d\bar s^2 = -2(a - 2 b r - q r^{-1})dt^2 + 2 dt dr$ and $d\tilde s^2 = - \frac{2}{f^2} d\zeta d\bar\zeta$.\\
\indent The main object of the present paper is to investigate the pseudosymmetric type structures admitting by the Robinson-Trautman metric \eqref{rtm}. Lanczos potentials of the metric are studied by Ahsan and Bilal \cite{AZBM10}. In the study of differential geometry, pseudosymmetric structures play an important role due their applications in relativity and cosmology (see, \cite{DDVV94}, \cite{DHV04}, \cite{DHJKS14}, \cite{DVV91}, \cite{SK16srs}, \cite{SM17} and also references therein). It is noteworthy to mention that the present paper is unconventional as by considering a physically relevant metric, viz., Robinson-Trautman metric, we would like to investigate the curvature restricted geometric structures admitted by such metric. Then Robison-Trautman spacetime can be considered as a model of such pseudosymmetric structures.\\
\indent The paper is organized as follows. Section 2 deals with the different pseudosymmetric structures which are essential to investigate the curvature properties of the Robinson-Trautman metric. Section 3 is concerned with different components of several curvature tensors which are useful to investigate the pseudosymmetric type structures admitting by this metric. Section 4 is devoted to the conclusion as theorems and remarks. Section 5 deals with the investigation under which the energy momentum tensor is covariantly constant, Codazzi type and cyclic parallel. In the last section we have made a comparison (similarities and dissimilarities) between the curvature properties of Robinson-Trautman metric and Som-Raychaudhuri metric. It is shown that both the metric are 2-quasi-Einstein and pseudosymmetric Weyl conformal curvature tensor, whereas Robinson-Trautman metric is Deszcz pseudosymmetric but Som-Raychaudhuri metric is Ricci generalized pseudosymmetric.
\section{\bf Curvature restricted geometric structures}
 It is well known that a curvature restricted geometric structure is a geometric structure on $M$ obtained by imposing a restriction on its curvature tensors by means of covariant derivatives of first order or higher orders. We will now explain some useful notations and definitions of various curvature restricted geometric structures. For two symmetric $(0,2)$-tensors $A$ and $E$, their Kulkarni-Nomizu product $A\wedge E$ is defined as (see e.g. \cite{DGHS11}, \cite{Glog02}):
\begin{eqnarray*}
(A\wedge E)(X_1,X_2,X_3,X_4) &=& A(X_1,X_4)E(X_2,X_3) + A(X_2,X_3)E(X_1,X_4)\\
&-& A(X_1,X_3)E(X_2,X_4) - A(X_2,X_4)E(X_1,X_3),
\end{eqnarray*}
where $X_1,X_2,X_3,X_4 \in \chi(M)$, the Lie algebra of all smooth vector fields on $M$. Throughout the paper we will consider $X, Y, X_1, X_2, \cdots \in \chi(M)$.\\
\indent Again as an invariant under certain class of transformations on a semi-Riemannian manifold $M$, there arise various curvature tensors on $M$, such as conformal curvature tensor $C$, projective curvature tensor $P$, concircular curvature tensor $W$ and conharmonic curvature tensor $K$ (\cite{Ishi57}, \cite{YK89}) etc. In terms of Kulkarni-Nomizu product $C$, $W$, $K$ and the Gaussian curvature tensor $\mathfrak G$ are given by
$$C = R-\frac{1}{n-2}(g\wedge S) + \frac{r}{2(n-2)(n-1)}(g\wedge g),$$
$$W = R-\frac{r}{2n(n-1)}(g\wedge g),$$
$$K = R-\frac{1}{n-2}(g\wedge S) \ \mbox{ and } \ \mathfrak G = \frac{1}{2}(g\wedge g).$$
The projective curvature tensor $P$ is given by
$$
P(X_1, X_2, X_3, X_4) = R(X_1, X_2, X_3, X_4) - \frac{1}{n-1}[g(X_1, X_4)S(X_2, X_3)-g(X_2, X_4)S(X_1, X_3)].
$$
A $(0,4)$ tensor is called a generalized curvature tensor if it posses the symmetry like $R$, i.e., a $(0,4)$ tensor $D$ is said to be a generalized curvature tensor (\cite{DGHS11}, \cite{SDHJK15}, \cite{SK14}) if
$$D(X_1,X_2,X_3,X_4)+D(X_2,X_3,X_1,X_4)+D(X_3,X_1,X_2,X_4)=0,$$
$$D(X_1,X_2,X_3,X_4)+D(X_2,X_1,X_3,X_4)=0 \ \ \mbox{and}$$
$$D(X_1,X_2,X_3,X_4)=D(X_3,X_4,X_1,X_2).$$
We note that if $A$ and $E$ are both symmetric $(0,2)$-tensors, then $A\wedge E$ is a generalized curvature tensor, and linear combination of two generalized curvature tensors is also a generalized curvature tensor. Hence $C$, $W$ and $K$ are generalized curvature tensors but $P$ is not a generalized curvature tensor.\\
\indent Now for a $(0,4)$-tensor $D$ and a symmetric $(0,2)$-tensor $E$ we can define three endomorphisms $\mathcal E$, $\mathscr{D}(X,Y)$ and $X\wedge_E Y$ as (\cite{DDHKS00}, \cite{DGHS11})
$$
E(X, Y) = g(\mathcal E X, Y), \ \mathscr{D}(X,Y)X_1 = \mathcal D(X,Y)X_1 \ \mbox{and} \ (X\wedge_A Y)X_1 = A(Y,X_1)X-A(X,X_1)Y
$$
respectively, where $\mathcal D$ is the corresponding $(1,3)$-tensor of $D$, i.e.,
$$D(X_1,X_2,X_3,X_4) = g(\mathcal D(X_1,X_2)X_3, X_4).$$
\indent One can now easily operate $\mathscr{D}(X,Y)$ and $X\wedge_A Y$ on a $(0,k)$-tensor $H$, $k\geq 1$, and obtain two $(0,k+2)$-tensors $D\cdot H$ and $Q(A,H)$ respectively as follows (see \cite{DG02}, \cite{DGHS98}, \cite{DH03}, \cite{SDHJK15}, \cite{SK14}, \cite{SKgrt}, \cite{Tach74} and also references therein):
$$D\cdot H(X_1,X_2,\cdots,X_k,X,Y) = -H(\mathcal D(X,Y)X_1,X_2,\cdots,X_k) - \cdots - H(X_1,X_2,\cdots,\mathcal D(X,Y)X_k).$$
and
\beb
&&Q(A,H)(X_1,X_2, \ldots ,X_k,X,Y) = ((X \wedge_A Y)\cdot H)(X_1,X_2, \ldots ,X_k)\\
&&= A(X, X_1) H(Y,X_2,\cdots,X_k) + \cdots + A(X, X_k) H(X_1,X_2,\cdots,Y)\\
&& - A(Y, X_1) H(X,X_2,\cdots,X_k) - \cdots - A(Y, X_k) H(X_1,X_2,\cdots,X).
\eeb
In terms of local coordinates system, $D\cdot H$ and $Q(A,H)$ \cite{Desz03a} can be written as
\beb
(D\cdot H)_{i_1 i_2\cdots i_k j l} &=& -g^{pq}\left[D_{jli_1 q}H_{p i_2\cdots i_k} + \cdots + D_{jli_k q}H_{i_1 i_2\cdots p}\right],
\eeb
\beb
Q(A, H)_{i_1 i_2\cdots i_k j l} &=& A_{li_1}H_{j i_2\cdots i_k} + \cdots + A_{li_k}H_{i_1 i_2\cdots j}\\
																	& & -A_{ji_1}H_{l i_2\cdots i_k} - \cdots - A_{ji_k}H_{i_1 i_2\cdots l}.
\eeb
\begin{defi}$($\cite{AD83}, \cite{Cart46}, \cite{Desz92}, \cite{SK14}, \cite{SKppsn}, \cite{SKppsnw}, \cite{SRK15}, \cite{Szab82}$)$
A semi-Riemannian manifold $M$ is said to be $H$-semisymmetric type if $D\cdot H = 0$ and it is said to be $H$-pseudosymmetric type if $\left(\sum\limits_{i=1}^k c_i D_i\right)\cdot H = 0$ for some scalars $c_i$'s, where $D$ and each $D_i$, $i=1,\ldots, k$, $(k\ge 2)$, are (0,4) curvature tensors. In particular, if $c_i$'s are all constants, then the manifold is called $H$-pseudosymmetric type manifold of constant type or otherwise non-constant type.
\end{defi}
\indent In particular, if $D = R$ and $T=R$ (resp., $S$, $C$, $W$ and $K$), then $M$ is called semisymmetric (resp., Ricci, conformally, concircularly and conharmonically semisymmetric). Again, if $i =2$, $D_1 = R$, $D_2 = \mathfrak G$ and $D= R$ (resp., $S$, $C$, $W$ and $K$), then $M$ is called Deszcz pseudosymmetric (resp., Ricci, conformally, concircularly and conharmonically pseudosymmetric). Especially, if $i =2$, $D_1 = C$, $D_2 = \mathfrak G$ and $D =C$, then $M$ is called a manifold of pseudosymmetric Weyl conformal curvature tensor. Again, if $i =2$, $D_1 = R$, $D_2 = \wedge_S$ and $D = R$, then $M$ is called Ricci generalized pseudosymmetric \cite{DD91a}.\\
\indent A semi-Riemannian manifold $M$ is said to be Einstein if its Ricci tensor is a scalar multiple of the metric tensor $g$. We note that in this case $S=\frac{\kappa}{n}g$ and $\kappa$ is scalar curvature. Again $M$ is called quasi-Einstein if at each point of $M$, rank of $(S - \alpha g)$ is less or equal to 1 for a scalar $\alpha$. Thus in this case $S = \alpha g + \beta \Pi \otimes \Pi$, where $\alpha$ and $\beta$ are scalars and $\Pi$ is an 1-form. In particular, if $\alpha = 0$, then a quasi-Einstein manifold is called Ricci simple. We note that G\"{o}del spacetime is Ricci simple \cite{DHJKS14} and perfect fluid spacetimes are quasi-Einstein (see, \cite{SB05}, \cite{SB06}). For quasi-Einstein spacetimes we refer the reader to see \cite{SYH09} and also references therein.
\begin{defi}
A semi-Riemannian manifold is said to be 2-quasi-Einstein if at each point of $M$, rank of $(S - \alpha g)$ is less or equal to 2 for a scalar $\alpha$.
\end{defi}
\indent It may be noted that Som-Raychaudhuri spacetime is 2-quasi-Einstein but not quasi-Einstein \cite{SK16srs}. For recent results on 2-quasi-Einstein manifolds, we refer the reader to see \cite{DGHZ-2016}, \cite{DGJZ-2016} and \cite{DGP-TV-2015} .
\begin{defi} \cite{Shai09}
A semi-Riemannian manifold $M$ is said to be pseudo quasi-Einstein if
\be\label{pqe}
S = \alpha g + \beta \Pi \otimes \Pi + \gamma E
\ee
holds on $U$ for some scalars $\alpha, \beta, \gamma$, an 1-form $\Pi$ and a trace free $(0,2)$-tensor $E$ such that the associated vector field $V$ corresponding to $\Pi$ satisfies $E(X,V) =0$.
\end{defi}
\indent In particular, if $E = \Phi \otimes \Phi$ (resp., $\Pi \otimes \Phi + \Phi \otimes \Pi$) then a pseudo quasi-Einstein manifold reduces to a generalized quasi-Einstein manifold by De and Ghosh \cite{DG04} (resp., by Chaki \cite{Chak01}).\\
\indent To get the decomposition of Riemann-Christoffel curvature tensor $R$, Deszcz (\cite{Desz03}, \cite{Desz03a}) introduced the notion of Roter type manifold, which was latter extended as generalized Roter type manifold \cite{SDHJK15}.
\begin{defi}
A semi-Riemannian manifold $M$ is said to be a Roter type manifold $($\cite{Desz03}, \cite{Desz03a}, \cite{DGHS11}, \cite{DGP-TV-2011}, \cite{DPSch-2013} and \cite{Glog-2007}$)$ $($resp., generalized Roter type manifold $($\cite{DGJPZ13}, \cite{DGJZ-2016}, \cite{DGP-TV-2015}, \cite{SDHJK15}, \cite{SKgrt}, \cite{SK16} and \cite{Saw-2015}$))$ if
$$R = N_1 (g\wedge g) + N_2 (g\wedge S) + N_3 (S\wedge S) \ \mbox{ and}$$
$$\left[\mbox{resp., } R = L_1 (g\wedge g) + L_2 (g\wedge S) + L_3 (S\wedge S) + L_4 (g\wedge S^2) + L_5 (S\wedge S^2) + L_6 (S^2\wedge S^2)\right]$$
holds for some $N_i, L_j \in C^{\infty}(M)$, $1\le i\le 3$ and $1\le j\le 6$.
\end{defi}
\begin{defi}(\cite{Bess87}, \cite{SKgrt})
A semi-Riemannian manifold $M$ is said to be $Ein(2)$ if $S^2$, $S$ and $g$ are linearly dependent, where $S^2$ is the second level Ricci tensor defined as $S^2(X,Y) = S(\mathcal S X, Y)$.
\end{defi}
\indent Between the class of Ricci symmetric manifolds and the class of manifolds of constant scalar curvature, Gray \cite{Gray78} presented two new classes of manifolds, namely, manifolds of Codazzi type Ricci tensor and manifolds of cyclic parallel Ricci tensor.
\begin{defi}
A semi-Riemannian manifold $M$ is said to be of Codazzi type Ricci tensor (resp. cyclic parallel Ricci tensor) (see, \cite{DHJKS14}, \cite{Gray78} and references therein) if
\[(\nabla_{X_1} S)(X_2, X_3) = (\nabla_{X_2} S)(X_1, X_3)\]
\[\left(\mbox{resp.} \ \ \ (\nabla_{X_1} S)(X_2, X_3) + (\nabla_{X_2} S)(X_3, X_1) + (\nabla_{X_3} S)(X_1, X_2) = 0 \ \right).\]
\end{defi}
Recently Mantica and Molinari showed \cite{MM12a} that the sufficient condition of Codazzi type for Derdzinski and Shen theorem \cite{DS83} can be replaced by a more general condition, and they introduced the notion of Riemann-compatible tensor along with  compatible tensors for other curvature tensors.
\begin{defi}\label{def2.8}
Let $D$ be a $(0,4)$-tensor and $E$ be a symmetric $(0, 2)$-tensor on $M$. Then $E$ is said to be $D$-compatible (\cite{DGJPZ13}, \cite{MM12b}, \cite{MM13}) if
\[
D(\mathcal E X_1, X,X_2,X_3) + D(\mathcal E X_2, X,X_3,X_1) + D(\mathcal E X_3, X,X_1,X_2) = 0
\]
holds. Again an 1-form $\Pi$ is said to be $D$-compatible if $\Pi\otimes \Pi$ is $D$-compatible.
\end{defi}
\indent Generalizing the concept of recurrent manifold (\cite{Ruse46}, \cite{Ruse49a}, \cite{Ruse49b}, \cite{Walk50}), recently Shaikh et al. \cite{SRK16} introduced the notion of super generalized recurrent manifold along with its characterization and existence by proper example.
\begin{defi}
A semi-Riemannian manifold $M$ is said to be super generalized recurrent manifold (\cite{SK14}, \cite{SRK16}, \cite{SKA16}) if
$$
\nabla R = \Pi \otimes R + \Phi \otimes (S\wedge S) + \Psi \otimes (g\wedge S) + \Theta \otimes (g\wedge g)
$$
holds on $\{x\in M: R \neq 0 \mbox{ and any one of } S\wedge S, g\wedge S \mbox{ is non-zero at $x$}\}$ for some 1-forms $\Pi$, $\Phi$, $\Psi$ and $\Theta$, called the associated 1-forms. Especially, if $\Phi = \Psi = \Theta = 0$ (resp., $\Psi = \Theta = 0$ and $\Phi = \Theta = 0$), then the manifold is called recurrent (\cite{Ruse46}, \cite{Ruse49a}, \cite{Ruse49b}, \cite{Walk50}) (resp., weakly generalized recurrent (\cite{SR11}, \cite{SAR13}) and hyper generalized recurrent (\cite{SP10}, \cite{SRK15})) manifold.
\end{defi}
\indent Again as a generalization of locally symmetric manifold and recurrent manifold, Tam$\acute{\mbox{a}}$ssy and Binh \cite{TB89} introduced the notion of weakly symmetric manifolds.
\begin{defi}
Let $D$ be a (0, 4)-tensor on a semi-Riemannian manifold $M$. Then $M$ is said to be weakly $D$-symmetric (\cite{TB89}, \cite{SK12}) if
\beb
&&\nabla_X  D(X_1, X_2, X_3, X_4) = \Pi(X) D(X_1, X_2, X_3, X_4) + \Phi(X_1) D(X_1, X_2, X_3, X_4)\\
&& \hspace{2cm} + \overline \Phi(X_2) D(X_1, X, X_3, X_4) + \Psi(X_3) D(X_1, X_2, X, X_4) + \overline \Psi(X_4) D(X_1, X_2, X_3, X)
\eeb
holds $\forall~ X, X_i \in \chi(M)$ $(i =1,2,3,4)$ and some 1-forms $\Pi, \Phi, \overline \Phi, \Psi$ and $\overline \Psi$  on $\{x\in M: R_x \neq 0\}$. In particular, if $\frac{1}{2}\Pi = \Phi = \overline \Phi = \Psi = \overline \Psi$, then the manifold is called Chaki $D$-pseudosymmetric manifold \cite{Chak87}.
\end{defi}
\indent It is noteworthy to mention that if $D$ is a generalized curvature tensor then the above defining condition of weak symmetry holds for $\frac{\Phi + \Psi}{2}$ in lieu of $\Phi, \overline \Phi, \Psi$ and $\overline \Psi$. It is also noted that the notion of Chaki pseudosymmetry is different from Deszcz pseudosymmetry. For details about the defining condition of weak symmetry and the interrelation between weak symmetry and Deszcz pseudosymmetry, we refer the reader to see \cite{SDHJK15} and also references therein.\\
\indent For a generalized curvature tensor $D$ and a symmetric (0, 2)-curvature tensor $Z$, one can define the curvature 2-form $\Omega_{(D)l}^m$ (\cite{Bess87}, \cite{LR89}) for $D$ and an 1-form $\Lambda_{(Z)l}$ \cite{SKP03} as
$$\Omega_{(D)l}^m= D_{jkl}^m dx^j \wedge dx^k \mbox{ and } \Lambda_{(Z)l} = Z_{lm} dx^m,$$
where $\wedge$ indicates the exterior product. Now $\Omega_{(D)l}^m$ (resp., $\Lambda_{(Z)l}$) are recurrent if
$$\mathcal D \Omega_{(R)l}^m  = \Pi \wedge \Omega_{(R)l}^m \mbox{ (resp., } \mathcal D \Lambda_{(S)l}  = \Pi \wedge \Lambda_{(S)l}),$$
where $\mathcal D$ is the exterior derivative and $\Pi$ is the associated $1$-form. Recently Mantica and Suh (\cite{MS12a}, \cite{MS13a}, \cite{MS14}) showed that $\Omega_{(D)l}^m$ are recurrent if and only if
\beb\label{man}
&&(\nabla_{X_1} D)(X_2,X_3,X,Y)+(\nabla_{X_2} D)(X_3,X_1,X,Y)+(\nabla_{X_3} D)(X_1,X_2,X,Y) =\\
&&\hspace{1in} \Pi(X_1) D(X_2,X_3,X,Y) + \Pi(X_2) D(X_3,X_1,X,Y)+ \Pi(X_3) D(X_1,X_2,X,Y)
\eeb
and $\Lambda_{(Z)l}$ are recurrent if and only if
$$(\nabla_{X_1} Z)(X_2,X) - (\nabla_{X_2} Z)(X_1,X) = \Pi(X_1) Z(X_2,X) - \Pi(X_2) Z(X_1,X)$$
for an 1-form $\Pi$.
\section{\bf Curvature properties of Robinson-Trautman metric}\label{com}
Setting $\zeta = x^3 + i x^4$, the metric \eqref{rtm} can be written as
\be\label{RTII}
ds^2 = -2(a - 2 b r - q r^{-1})dt^2 + 2 dt dr - \frac{r^2}{f^2} [(dx^3)^2+(dx^4)^2],
\ee
where $f$ is a function of the real variables $x^3$ and $x^4$.
Hence the metric tensor in coordinates $(t,r,x^3,x^4)$ is given by
$$g = \left(
\begin{array}{cccc}
 -2(a - 2 b r - q r^{-1}) & 1 & 0 & 0 \\
 1 & 0 & 0 & 0 \\
 0 & 0 & - \frac{r^2}{f^2} & 0 \\
 0 & 0 & 0 & - \frac{r^2}{f^2}
\end{array}
\right).$$
Again $g$ can be expressed as warped product $\bar g \times_r \tilde g$ with warping function $r$, where base metric $\bar g$ and fiber metric $\tilde g$ are given by
$$\bar g = \left(
\begin{array}{cc}
 -2(a - 2 b r - q r^{-1}) & 1 \\
 1 & 0
\end{array}
\right), \ \
\tilde g = \left(
\begin{array}{cc}
 - \frac{1}{f^2} & 0 \\
 0 & - \frac{1}{f^2}
\end{array}
\right).$$
From above it is clear that the metric is Lorentzian with signature $\{-,+,+,+\}$.\\
\indent The non-zero components of its Riemann-Christoffel curvature tensor $R$, Ricci tensor $S$ and scalar curvature $\kappa$ are given by
$$R_{1212}=-\frac{2 q}{r^3}, \ R_{1313}= R_{1414}=-\frac{2 \left(2 b r^2-q\right) \left(-a r+2 b r^2+q\right)}{f^2 r^2},$$
$$R_{1323}= R_{1424}=-\frac{2 b r^2-q}{f^2 r}, \ R_{3434}=\frac{r \left(-2 a r+4 b r^2+2 q+F r\right)}{f^4}.$$
$$S_{11}=-\frac{8 b \left(-a r+2 b r^2+q\right)}{r^2}, \ S_{12}=-\frac{4 b}{r}, \ S_{33}= S_{44}=-\frac{2 a-8 b r-F}{f^2}$$
$$\mbox{and } \kappa = -\frac{2 (-2 a+12 b r+F)}{r^2},$$
where $F = f_3^2+f_4^2-f (f_{33}+f_{44})$.
Also the non-zero components of $C$ and $P$ are given by
$$C_{1212}=\frac{2 a r-6 q-F r}{3 r^3}, \ \ C_{1313}= C_{1414}= -\frac{\left(-a r+2 b r^2+q\right) (2 a r-6 q-F r)}{3 f^2 r^2}$$
$$C_{1323}= C_{1424}=-\frac{2 a r-6 q-F r}{6 f^2 r}, \ \ C_{3434}=-\frac{r (2 a r-6 q-F r)}{3 f^4},$$
$$-P_{1212}= P_{1221}=\frac{2 \left(3 q-2 b r^2\right)}{3 r^3}, \ P_{1313}= P_{1414}=-\frac{2 \left(2 b r^2-3 q\right) \left(-a r+2 b r^2+q\right)}{3 f^2 r^2},$$
$$P_{1323}= P_{1424}= P_{2313}= P_{2414}=\frac{3 q-2 b r^2}{3 f^2 r},$$
$$P_{1331}= P_{1441}=-\frac{2 \left(-a r+2 b r^2+q\right) \left(-2 a r+2 b r^2+3 q+F r\right)}{3 f^2 r^2},$$
$$P_{1332}= P_{1442}= P_{2331}= P_{2441}=-\frac{-2 a r+2 b r^2+3 q+F r}{3 f^2 r},$$
$$P_{3434}= -P_{3443}=\frac{2 r \left(-2 a r+2 b r^2+3 q+F r\right)}{3 f^4}.$$
The non-zero components of $\nabla R$ and $\nabla S$ are given by
$$R_{1212,2}=\frac{6 q}{r^4}, \ R_{1213,3}= R_{1214,4}= R_{1313,2}= R_{1414,2}=\frac{2 \left(2 b r^2-3 q\right) \left(-a r+2 b r^2+q\right)}{f^2 r^3},$$
$$R_{1223,3}= R_{1224,4}= R_{1323,2}= R_{1424,2}=\frac{2 b r^2-3 q}{f^2 r^2}, \ R_{3434,3}=\frac{F_3 r^2}{f^4},$$
$$R_{2334,4}= -R_{2434,3}= -\frac{1}{2}R_{3434,2}=\frac{-2 a r+2 b r^2+3 q+F r}{f^4}, \ R_{3434,4}=\frac{F_4 r^2}{f^4},$$
$$S_{11,2}=\frac{8 b \left(-a r+2 b r^2+q\right)}{r^3}, \ S_{33,3}= S_{44,3}=\frac{F_3}{f^2},\ S_{33,4}= S_{44,4}=\frac{F_4}{f^2},$$
$$S_{12,2}=\frac{4 b}{r^2}, \ 2S_{23,3}= 2S_{24,4}= S_{33,2}= S_{44,2}=-\frac{2 (-2 a+4 b r+F)}{f^2 r}.$$
The non-zero components of $R\cdot R$, $Q(g,R)$ and $Q(S,R)$ are given by
$$-R\cdot R_{121323}= -R\cdot R_{121424}= R\cdot R_{122313}= R\cdot R_{122414}=\frac{\left(2 b r^2-3 q\right) \left(2 b r^2-q\right)}{f^2 r^4},$$
$$R\cdot R_{133414}= -R\cdot R_{143413}=\frac{2 \left(2 b r^2-q\right) \left(-a r+2 b r^2+q\right) \left(-2 a r+2 b r^2+3 q+F r\right)}{f^4 r^3},$$
$$R\cdot R_{133424}= -R\cdot R_{143423}= R\cdot R_{233414}= -R\cdot R_{243413}=\frac{\left(2 b r^2-q\right) \left(-2 a r+2 b r^2+3 q+F r\right)}{f^4 r^2}.$$
$$Q(g,R)_{121323}= Q(g,R)_{121424}= -Q(g,R)_{122313}= -Q(g,R)_{122414}=\frac{2 b r^2-3 q}{f^2 r},$$
$$-Q(g,R)_{133414}= Q(g,R)_{143413}=\frac{2 \left(-a r+2 b r^2+q\right) \left(-2 a r+2 b r^2+3 q+F r\right)}{f^4},$$
$$-Q(g,R)_{133424}= Q(g,R)_{143423}= -Q(g,R)_{233414}= Q(g,R)_{243413}=\frac{r \left(-2 a r+2 b r^2+3 q+F r\right)}{f^4},$$
$$-Q(S,R)_{121323}= -Q(S,R)_{121424}= Q(S,R)_{122313}= Q(S,R)_{122414}=\frac{4 a q+8 b^2 r^3-20 b q r-2q F}{f^2 r^3},$$
$$-Q(S,R)_{133414}= Q(S,R)_{143413}=\frac{2 \left(-a r+2 b r^2+q\right) \left(4 a b r^2+2 a q-16 b q r-2 b F r^2-q F\right)}{f^4 r^2},$$
$$Q(S,R)_{133424}= Q(S,R)_{143423}= Q(S,R)_{233414}= Q(S,R)_{243413}$$
$$=-\frac{-4 a b r^2-2 a q+16 b q r+2 b F r^2+q F}{f^4 r}.$$
\section{\bf Conclusion}\label{con}
 From above we can conclude that the metric \eqref{RTII} fulfills the following curvature restricted geometric structures.
\begin{thm}
The Robinson-Trautman metric of Petrov type II given in \eqref{RTII} possesses the following curvature properties:
\begin{enumerate}[label=(\roman*)]
    \item It is Deszcz pseudosymmetric with $R\cdot R = \frac{q - 2 b r^2}{r^3} Q(g,R)$. Hence it is Ricci pseudosymmetric, conformally pseudosymmetric, projectively pseudosymmetric, concircularly pseudosymmetric and conharmonically pseudosymmetric of non-constant type with same associated scalar.

		\item  It satisfies
		$${r^3 \left(2 a r-6 q- F r\right)}\left[R\cdot R - Q(S, R)\right]= {2 \left(4 a q r+4 b^2 r^4-12 b q r^2-3 q^2-2 q F r\right)} Q(g,C).$$

		\item  It satisfies $C\cdot R = L Q(g, R)$, $L = -\frac{\left(2 a r-6 q - F r\right)}{6 r^3}$. Hence $C\cdot S = L Q(g, S)$, $C\cdot C = L Q(g, C)$ (i.e., pseudosymmetric Weyl conformal curvature tensor), $C\cdot W = L Q(g, W)$, $C\cdot P = L Q(g, P)$ and $C\cdot K = L Q(g, K)$.

		\item If $\left(4 q r (a-3 b r)+4 b^2 r^4-3 q^2-2 q F r\right)$ is nowhere zero, then it satisfies $C\cdot R - R\cdot C = L_1 Q(g,R) +L_2 Q(S,R)$ for
		$$L_1 = -\frac{(2 a r-6 q-F r) \left(2 \left(a q+4 b^2 r^3-6 b q r\right)-q F\right)}{3 r^2 \left(4 q r (a-3 b r)+4 b^2 r^4-3 q^2-2 q F r\right)}$$
		and
		$$L_2 = \frac{\left(q-2 b r^2\right) (-2 a r+6 q+F r)}{8 q r (3 b r-a)-8 b^2 r^4+6 q^2+4 q F r}$$

		\item  It satisfies $C\cdot R - R\cdot C = L_1 Q(g,C) +L_2 Q(S,C)$ for
		$$L_1 = \frac{2 (-2 a+12 b r+F)}{3 r^2} \mbox{ and } L_2 = 1$$

		\item If $-2 a r+6 q+F r$ is everywhere nonzero, then its conformal curvature 2-forms are recurrent for the 1-form
		$$\Pi = \left\{0, \frac{F-2 a}{-2 a r+6 q+F r}, \frac{F_3 r}{-2 a r+6 q+F r}, \frac{F_4 r}{-2 a r+6 q+F r}\right\}.$$

		\item If $-2 a+4 b r+F$ is everywhere nonzero, then this metric is a Roter type with $R = N_1 (S\wedge S) + N_2 (g\wedge S) + N_3 (g\wedge g),$ where
		$$N_1 = \frac{r \left(2 a r-6 q-F r\right)}{2 \left(-2 a+4 b r+F\right)^2},$$
		$$N_2 = \frac{4 a b r^2+6 a q+8 b^2 r^3-F \left(2 b r^2+3 q\right)-36 b q r}{r \left(-2 a+4 b r+F\right)^2}$$
		$$N_3 = -\frac{1}{r^3}\left[\-\frac{4 b r \left(6 a q+8 b^2 r^3-24 b q r-3 q F\right)}{\left(-2 a+4 b r+F\right)^2}+q\right].$$

		\item This is 2-quasi Einstein as $S-\alpha g$ is of rank 2 for $\alpha = \frac{2 a-8 b r-F}{r^2}$ or $\alpha = -\frac{4 b}{r}$.

		\item This is generalized quasi Einstein in the sense of Chaki for
		$$\alpha = -\frac{-2 a+8 b r+F}{r^2}, \ \beta = \frac{4 b \left(-2 a+4 b r+F\right)}{\Phi_1^2 r},$$
		$$\gamma = -2 a+4 b r+F, \ \ \Pi =\left\{\frac{q-a r}{\Phi_1 r^3},\frac{1}{\Phi_1 r^2},0,0\right\}, \ \ \Phi= \left\{\Phi_1,0,0,0\right\},$$
		where $\Phi_1$ being arbitrary scalar, $||\Pi|| = -\frac{4 b}{\eta (1)^2 r^3}$ and $||\Phi|| = 0$.

		\item This is generalized quasi Einstein in the sense of De and Ghosh for
		$$\alpha = -\frac{4 b}{r}, \ \beta = \gamma = \frac{-2 a+4 b r+F}{r^2},$$
		$$\Pi =\left\{0,0,\frac{r}{f},0\right\}, \ \ \Phi= \left\{0,0,0,\frac{r}{f}\right\},$$
		where $||\Pi|| = ||\Phi|| = -1$ and $\Pi$ is orthogonal to $\Phi$.

		\item This is pseudo quasi Einstein for
		$$\alpha = -\frac{F}{2 r^2}, \ \beta = -\frac{4 (a-6 b r)}{r^2}, \ \gamma = 1, \ \Pi =\left\{0,0,0,\frac{r}{f}\right\}.$$

		\item A (0,2) tensor of the form
$$\left(
\begin{array}{cccc}
 t(1,1) & t(1,2) & 0 & 0 \\
 t(1,2) & t(2,2) & 0 & 0 \\
 0 & 0 & t(3,3) & t(3,4) \\
 0 & 0 & t(3,4) & t(4,4)
\end{array}
\right),$$
where $t(i,j)$ being arbitrary scalar, is $R$-compatible as well as $C$-compatible, $W$-compatible and $K$-compatible, where $t(i,j)$ being arbitrary scalar. So the Ricci tensor is $R$-compatible as well as $C$-compatible. Again the vector of the form $\{t_1,0,0,0\}$ or $\{0,t_2,0,0\}$ or $\{0,0,t_3,0\}$ or $\{0,0,0,t_4\}$, $t_i$ being arbitrary scalar, is $R$-compatible as well as $C$-compatible.
		\end{enumerate}
\end{thm}
\begin{rem}
From the value of the local components (presented in Section \ref{com}) of various tensors of the Robinson-Trautman metric \eqref{RTII}, we can easily conclude that the metric does not fulfill the following geometric structures:
\begin{enumerate}[label=(\roman*)]
\item not conformally symmetric and hence not locally symmetric, projectively symmetric, concircularly symmetric, conharmonically symmetric.
\item not conformally recurrent and hence not recurrent, projectively recurrent, concircularly recurrent, conharmonically recurrent.
\item not super generalized recurrent and hence not hyper generalized recurrent, weakly generalized recurrent.
\item not weakly symmetric for $R$, $C$, $P$, $W$ and $K$ and hence not Chaki pseudosymmetric for $R$, $C$, $P$, $W$ and $K$.
\item not weakly Ricci symmetric and hence not pseudo Ricci symmetric in the sense of Chaki and not Ricci symmetric.
\indent Scalar curvature is not constant and hence its Ricci tensor is not Codazzi type, not cyclic parallel.
\item $div R \ne 0$, $div C \ne 0$, $div P \ne 0$, $div W \ne 0$ and $div K \ne 0$.
\item not conformally semisymmetric and hence not semisymmetric, projectively semisymmetric, concircularly semisymmetric, conharmonically semisymmetric.
\item not Ricci generalized pseudosymmetric.
\item curvature 2-forms for $R$, $P$, $W$ and $K$ are not recurrent. The Ricci 1-forms are not recurrent.
\end{enumerate}
\end{rem}
\section{\bf Energy momentum tensor of Robinson-Trautman metric}\label{em-t}
From the Einstein's field equations, the energy momentum tensor, $T$ is given by
$$T= \frac{c^4}{8\pi G}\left[S-\left(\frac{\kappa}{2}-\Lambda\right)g\right],$$
where $c=$ speed of light in vacuum, $G=$ gravitational constant and $\Lambda=$ cosmological constant. Thus the components of the energy momentum tensor is calculated by considering the components of Ricci tensor and scalar curvature. Now the non-zero components of $T$ (upto symmetry) are given by
$$T_{11}=\frac{c^4 \left(-a r+2 b r^2+q\right) \left(-2 a+8 b r+F+\Lambda  r^2\right)}{4 \pi  G r^3}$$
$$T_{12}=\frac{c^4 \left(-2 a+8 b r+F+\Lambda  r^2\right)}{8 \pi  G r^2}, \ T_{33}= T_{44}=-\frac{c^4 r (4 b+\Lambda  r)}{8 \pi  f^2 G}.$$
Then the non-zero components of covariant derivative of $T$ are given by
$$T_{11,2}=-\frac{c^4 \left(-a r+2 b r^2+q\right) (-2 a+4 b r+F)}{2 \pi  G r^4}, \ T_{11,3}=\frac{c^4 F_3 \left(-a r+2 b r^2+q\right)}{4 \pi  G r^3},$$
$$T_{11,4}=\frac{c^4 F_4 \left(-a r+2 b r^2+q\right)}{4 \pi  G r^3}, \ T_{12,2}=-\frac{c^4 (-2 a+4 b r+F)}{4 \pi  G r^3},$$
$$T_{12,3}=\frac{c^4 F_3}{8 \pi  G r^2}, \ T_{12,4}=\frac{c^4 F_4}{8 \pi  G r^2},$$
$$T_{23,3}= T_{24,4}=-\frac{c^4 (-2 a+4 b r+F)}{8 \pi  f^2 G r}, \ T_{33,2}= T_{44,2}=\frac{b c^4}{2 \pi  f^2 G}.$$
\indent We note that from the above it is easy to check that $div T = 0$ which implies that the metric \eqref{RTII} satisfies energy condition. However it is obvious that $T$ is not covariantly constant as well as it is not Codazzi type and cyclic parallel. In this section we would like to investigate the conditions under which $T$ is parallel, $T$ is Codazzi type and $T$ is cyclic parallel.\\
\indent Now from the values of the components of $\nabla T$, we get
$$T_{11,2} + T_{12,1} + T_{21,1}=-\frac{c^4 \left(-a r+2 b r^2+q\right) (-2 a+4 b r+F)}{2 \pi  G r^4},$$
$$T_{11,3} + T_{13,1} + T_{31,1}=\frac{c^4 F_3 \left(-a r+2 b r^2+q\right)}{4 \pi  G r^3},$$
$$T_{11,4} + T_{14,1} + T_{41,1}=\frac{c^4 F_4 \left(-a r+2 b r^2+q\right)}{4 \pi  G r^3},$$
$$T_{12,2} + T_{22,1} + T_{21,2}=-\frac{c^4 (-2 a+4 b r+F)}{2 \pi  G r^3},$$
$$T_{12,3} + T_{23,1} + T_{31,2}=\frac{c^4 F_3}{8 \pi  G r^2}, \ T_{12,4} + T_{24,1} + T_{41,2}=\frac{c^4 F_4}{8 \pi  G r^2},$$
$$T_{23,3} + T_{33,2} + T_{32,3}= T_{24,4} + T_{44,2} + T_{42,4}=-\frac{c^4 (-2 a+2 b r+F)}{4 \pi  f^2 G r}$$
and
$$T_{12,1} - T_{11,2} = \frac{c^4 \left(-a r+2 b r^2+q\right) (-2 a+4 b r+F)}{2 \pi  G r^4},$$
$$T_{11,3} - T_{13,1} = \frac{c^4 F_3 \left(-a r+2 b r^2+q\right)}{4 \pi  G r^3}, \ \
T_{11,4} - T_{14,1} = \frac{c^4 F_4 \left(-a r+2 b r^2+q\right)}{4 \pi  G r^3},$$
$$T_{12,3} - T_{13,2} = T_{21,3} - T_{23,1} = \frac{c^4 F_3}{8 \pi  G r^2}, \ \ T_{12,4} - T_{14,2} = T_{21,4} - T_{24,1} = \frac{c^4 F_4}{8 \pi  G r^2},$$
$$T_{22,1} - T_{21,2} = \frac{c^4 (-2 a+4 b r+F)}{4 \pi  G r^3}, \ \ T_{33,2} - T_{32,3} =  T_{44,2} - T_{42,4} = \frac{c^4 (-2 a+8 b r+F)}{8 \pi  f^2 G r}.$$
\indent Hence from above we see that if the energy-momentum tensor $T$ is cyclic parallel or Codazzi type or covariantly constant, then $b = 0$. Therefore $T$ is cyclic parallel or Codazzi type or covariantly constant if and only if $b=0$ and $F=2a$.
Now we can state the following:
\begin{thm}
For the Robinson-Trautman metric, the following conditions are equivalent:\\
(i) the  energy-momentum tensor is covariantly constant\\
(ii) the  energy-momentum tensor is Codazzi type\\
(iii) the  energy-momentum tensor is cyclic parallel\\
(iv) $b=0$ and $F=2a$.
\end{thm}
\begin{exm}
We consider $a=b=0$ and $f(x^3, x^4) = e^{x^3+x^4}$ in \eqref{RTII}. Then $F=0$ and therefore the energy-momentum tensor $T$ is covariantly constant and hence cyclic parallel and Codazzi type.
\end{exm}
\indent Now since the Robinson-Trautman metric \eqref{RTII} is Ricci pseudosymmetric and $T$ is linear combination of $S$ and $g$, and $R\cdot g = Q(g,g) = 0$, then $R\cdot T =\frac{q - 2 b r^2}{r^3} Q(g, T)$. Thus we can state the following:
\begin{thm}
The energy-momentum tensor of Robinson-Trautman metric is pseudosymmetric type, i.e., $R\cdot T = \frac{q - 2 b r^2}{r^3} Q(g, T)$.
\end{thm}
\section{\bf Som-Raychaudhuri metric and Robinson-Trautman metric}\label{sr-rt}
Som-Raychaudhuri \cite{SR68} spacetime is a G\"odel type metric  (see \cite{DHJKS14}, \cite{RT83}, \cite{SK16srs} and references therein) and in cylindrical coordinates $(t, r, z, \phi)$, its line element is given by
\be\label{srm}
ds^2 = dt^2 - (r^2-a^2r^4) d\phi^2 - dr^2 - dz^2 + 2 a r^2 d\phi dt.
\ee
Recently, Shaikh and Kundu \cite{SK16srs} investigated the curvature restricted geometric structures admitting by the Som-Raychaudhuri spacetime and showed that such a spacetime is a 2-quasi-Einstein, generalized Roter type,  $Ein(3)$ manifold satisfying $R.R = Q(S,R)$, $C\cdot C = \frac{2a^2}{3} Q(g,C)$, and its Ricci tensor is cyclic parallel and Riemann compatible. In this section we study the comparison between the curvature properties of Som-Raychaudhuri metric and Robinson-Trautman metric.\\
\indent Physically, Som-Raychaudhuri \cite{SR68} spacetime is a stationary cylindrical symmetric solution of Einstein field equation corresponding to a charged dust distribution in rigid rotation whereas Robinson-Trautman metric admit a geodesic, hypersurface orthogonal, shear-free and expanding null congruences. Now from Section \ref{con} of this paper and Section 3 of \cite{SK16srs} we have the following comparison between Som-Raychaudhuri metric and Robinson-Trautman metric:\\
\textbf{A. Similarity:}\\
(i) Both are non-Einstein but 2-quasi-Einstein,\\
(ii) Both are pseudo quasi Einstein and generalized quasi-Einstein in the sense of Chaki as well as in the sense of De and Ghosh.\\
(iii) Both are of pseudosymmetric Weyl conformal curvature tensor.\\
(iv) Their Ricci tensors are Riemann compatible as well as conformal compatible, concircular compatible and conharmonic compatible.\\
\textbf{B. Dissimilarity:}\\
(i) Som-Raychaudhuri metric is Ricci generalized pseudosymmetric but Robinson-Trautman metric is Deszcz pseudosymmetric.\\
(ii) Som-Raychaudhuri metric is generalized Roter type but Robinson-Trautman metric is Roter type.\\
(iii) Ricci tensor of Som-Raychaudhuri metric is cyclic parallel whereas the scalar curvature of Robinson-Trautman metric is non-constant.\\
(iv) Robinson-Trautman metric is $Ein(2)$ but Som-Raychaudhuri metric is $Ein(3)$.


\section{Introduction} \label{sec:intro}
Fog networking is a novel paradigm in which computing, storage and communication functions are implemented at both cloud and edge nodes (ENs), such as base stations, of a wireless cellular system. As Fig. 1 shows,
content delivery can benefit from fog networking via edge caching (storing popular content at the ENs), as well as via cloud processing, which enables the delivery of content fetched from a central content library.

The information-theoretic analysis of edge caching in \cite{AliNiesen15,nad016fundamental,cao2016fundamental,hachem2016degrees} and of more general fog-assisted wireless networks, or Fog Radio Access Networks (F-RANs), in \cite{TandonSimeone16,sengupta2016cloud} makes the assumption that all files in the content library have the same timeliness constraint. Under this assumption, caching schemes in which all contents are allocated the same fraction of the ENs' caches were proven to be optimal or near-optimal in \cite{TandonSimeone16,sengupta2016cloud}. In practice, however, contents may have heterogeneous latency requirements; e.g. video chunks may be buffered to reduce the delay constraints, while information feeding an Augmented Reality (AR) application has stricter latency requirements. Reducing the delivery latency of a content type generally requires allocating a larger fraction of the ENs' cache capacity to it, which in turn increases the delivery latency of other contents.

As in \cite{AliNiesen15,nad016fundamental,cao2016fundamental,hachem2016degrees,TandonSimeone16,sengupta2016cloud}, in this paper, users are assumed to make simultaneous requests from a library of contents, which may be partially cached at the ENs during an offline caching phase. Unlike prior work, in which all request sets experience the same delivery coding latency, here we study the trade-offs among the latencies that are achievable across different request sets, when one allows arbitrary allocation of cache capacity at the ENs across files. Leveraging a fine-grained understanding of these trade-offs makes it possible to analyze individual content latency constraints, including the average latency for a content type under a probabilistic popularity model.

As in \cite{TandonSimeone16,sengupta2016cloud}, as well as in \cite{yi2016topological,cao2016fundamental}, delivery latencies are measured here in the high-Signal-to-Noise Radio (SNR) regime with respect to a reference interference-free system, yielding the performance metric of Normalized Delivery Time (NDT). In \cite{TandonSimeone16}, the minimum NDT was characterized for a system with two ENs and two users under the assumption that all users' requests should be guaranteed the same latency. Reference \cite{sengupta2016cloud} obtained upper and lower bounds that are within a multiplicative factor of 2 for any number of ENs and users.

\begin{figure}
\centering
\begin{tikzpicture}[font=\scriptsize]

\node[
  draw,
  align=center
  ] (cloud) at (-1.5,0) {$\{F_1,\dots,F_J\}$};
\node[above=0mm of cloud] {Cloud};

\begin{scope}[xshift=23mm,yshift=13mm]
\node[
  draw,
  cylinder, shape border rotate=90, minimum width=5.5mm, minimum height=4mm
  ] (muaa) at (-0.6,0) {};
\node at (-0.6,0) {$\mu_{1,1}$};

\node at (0,0) {$\cdots$};

\node[
  draw,
  cylinder, shape border rotate=90, minimum width=5.5mm, minimum height=4mm
  ] (muab) at (0.6,0) {};
\node at (0.6,0) {$\mu_{1,J}$};

\node[draw,fit=(muaa) (muab)] (ENa) {};
\node[above=0mm of ENa] {EN$_1$};

\node at (0,-0.5) {$\sum_{j=1}^J\mu_{1,j}\leq \mu J$};
\end{scope}

\begin{scope}[xshift=23mm,yshift=-13mm]
\node[
  draw,
  cylinder, shape border rotate=90, minimum width=5.5mm, minimum height=4mm
  ] (muba) at (-0.6,0) {};
\node at (-0.6,0) {$\mu_{2,1}$};

\node at (0,0) {$\cdots$};

\node[
  draw,
  cylinder, shape border rotate=90, minimum width=5.5mm, minimum height=4mm
  ] (mubb) at (0.6,0) {};
\node at (0.6,0) {$\mu_{2,J}$};

\node[draw,fit=(muba) (mubb)] (ENb) {};
\node[above=0mm of ENb] {EN$_2$};

\node at (0,-0.5) {$\sum_{j=1}^J\mu_{2,j}\leq \mu J$};
\end{scope}

\node[draw] (usera) at (5.2,1.4) {$F_{D_1}$};
\node[above=0mm of usera] {User 1};

\node[draw] (userb) at (5.2,-1.2) {$F_{D_2}$};
\node[above=0mm of userb] {User 2};

\draw[-latex] (cloud.east) -- node[above,sloped] {$C^f=r\log P$} (ENa.west);
\draw[-latex] (cloud.east) -- node[below,sloped] {$C^f=r\log P$} (ENb.west);
\draw[-latex] (ENa.east) -- (usera.west);
\draw[-latex] (ENa.east) -- (userb.west);
\draw[-latex] (ENb.east) -- (usera.west);
\draw[-latex] (ENb.east) -- (userb.west);

\node[draw, fill=white,
  ellipse,
  align=center,
  inner xsep=3mm,
  rotate=90] (channel) at (4.1,0) {Shared wireless\\ channel};

\end{tikzpicture}
\caption{Illustration of the F-RAN system under study for $M=2$ and $K=2$.}
\label{fig:fran}
\vspace{-12pt}
\end{figure}

Focusing on the special case with two ENs and two users as in \cite{TandonSimeone16}, the main contributions of this work are as follows: (\emph{i}) The performance metric of the \emph{NDT region} is introduced with the aim of analyzing the trade-off among the latencies achievable for individual users' requests under non-uniform cache partitions across files (Sec. II); (\emph{ii}) Novel achievable schemes are presented that yield an inner bound to the NDT region (Sec. IV); (\emph{iii}) Outer bounds on the NDT region are derived that conclusively characterize the NDT region (Sec. V); (\emph{iv}) Numerical results corroborate the analysis (Sec. VI).

\section{System Model} \label{ssec:model}
\subsection{Model}
We consider an F-RAN architecture with $M$ edge nodes (ENs), which serve $K$ users over a shared wireless channel, see Fig. \ref{fig:fran}. As in prior works \cite{AliNiesen15,nad016fundamental,cao2016fundamental,hachem2016degrees,TandonSimeone16,sengupta2016cloud}, the system operates in two separate phases, namely an (offline) caching phase and an (online) delivery phase. In both phases, the content library $\mathcal{F}=\{F_1, \dots, F_J\}$ of $J\geq K$ files, where each file is of length $L$ bits, is fixed and static. The assumption of equally-sized files simplifies the treatment, as in prior work, and should be alleviated in future studies. In the caching phase, each EN $m$ can cache at most $\mu J L$ bits from the library, where $0\leq \mu \leq 1$ is referred to as the \emph{fractional cache capacity}.

The delivery phase consists of an arbitrary number of slots. In any slot, each user $k$ requests a file $F_{D_k}$ in $\mathcal{F}$ with index $D_k \in [1:J]$. We let $D=[D_1,\dots,D_K]$ denote the vector of requested files in a slot. We make the assumption that the requested files are distinct, as e.g. in~\cite{pedarsani2016online}. In future work, we plan to alleviate this limitation. The main goal of this work is understanding the trade-offs achievable among the delivery latencies that are achievable for different request vectors $D$. As we discuss in Sec. \ref{sec:eval}, this fine-grained understanding of the trade-offs among the delivery latencies for different requests can be used to study individual latency requirements for different files under a given popularity distribution.

The channel from the ENs to the users is defined by:
\begin{equation}
Y_k = \sum_{m=1}^M H_{m,k}X_m + Z_k, \label{eq:rx_signal}
\end{equation}
which is a standard quasi-static model, where $X_m\in\mathbb{C}^{n^e}$ is a codeword of length $n^e$ symbols transmitted by the EN $m$; $H_{m,k}\in\mathbb{C}$ is the channel coefficient from EN $m$ to user $k$; $Z_k$ is complex Gaussian additive noise with unitary power, i.i.d. over time and users and also independent of the channel coefficients; and $Y_k\in\mathbb{C}^{n^e}$ is the received signal of length $n^e$ symbols by user $k$. The channel coefficients are realizations of continuous random variables, and are i.i.d. over ENs and users. Using the notation $[1:K]=\{1,\dots,K\}$, let $\mathcal{H}=(H_{m,k})_{m \in [1:M], k \in [1:K]}$ denote the channel state information (CSI), which is assumed to be known throughout the network, {i.e.},\  at the cloud, the ENs and the users.

The cloud has orthogonal fronthaul links to each of the ENs. Using the parametrization in \cite{sengupta2016cloud}, the capacity is measured in bits per symbol, where a \emph{symbol} is a channel use of the wireless channel. Furthermore, as in \cite{sengupta2016cloud}, the fronthaul capacity is written as $C^f = r\log(P)$, where $r$ is defined as the \emph{fronthaul rate} and $P$ is the (high) SNR of the wireless edge links. The fronthaul rate describes the ratio between the fronthaul capacity and the high-SNR capacity of each EN-to-user when used with no interference from other links.

 In the \emph{caching phase}, EN $m$ stores an arbitrary function \begin{equation}
S_{m,j} = \pi_{m,j}^{c}(F_j)
\end{equation}
of each file $F_j, j \in [1:J]$,. We allow for an arbitrary partition of each EN's cache capacity. Denoting the entropy of the cached content for file $F_j$ at EN $m$ as $H(S_{m,j})=\mu_{m,j}L$, with $0 \leq \mu_{m,j} \leq 1$, we impose that the cache partition $\{\mu_{m,j}\}_{j \in [1:J]}$ satisfies the per-EN cache capacity constraint
\begin{equation} \label{eq:capacityconstraint}
\sum_{j=1}^J\mu_{m,j}\leq\mu J,
\end{equation}
for each $m$. We will refer to $\pi^c = (\pi^c_{m,j})_{m \in [1:M], j \in [1:J]}$ as the caching policy and to the matrix
\begin{equation} \label{eq:mumatrix}
\bm{\mu}=(\mu_{m,j})_{m\in[1:M], j \in [1:J]}
\end{equation}
as the \emph{cache partition matrix} for the given caching policy.

\begin{figure}
\centering
\begin{tikzpicture}[font=\scriptsize]

\draw[thick] (.7,0) -- node[at end, below,align=left,inner xsep=1pt] (ntime) {normalized\\ time} (8.5,0);
\node[left=0pt of ntime,inner xsep=0pt] {$\rightarrow$};

\draw (1,-.3) -- (1,.3);
\draw (3,-.3) -- (3,.3);
\draw (4.5,-.3) -- (4.5,.3);
\draw[latex-latex] (1,-.25) -- node[fill=white,inner ysep=0pt] {$\delta_D^f$} (3,-.25);
\path (1,.2) -- node[fill=white, inner ysep=0pt] {fronth.\ tx} (3,.2);
\draw[latex-latex] (3,-.25) -- node[fill=white, inner ysep=0pt] {$\delta_D^e$} (4.5,-.25);
\path (3,.2) -- node[fill=white, inner ysep=0pt] {edge tx} (4.5,.2);
\draw[decorate,decoration=brace] (1,.4) -- node[above=1pt] {Delivery of $D$} (4.5,.4);

\draw (5,-.3) -- (5,.3);
\draw (6.2,-.3) -- (6.2,.3);
\draw (7.3,-.3) -- (7.3,.3);
\draw[latex-latex] (5,-.25) -- node[fill=white,inner ysep=0pt] {$\delta_{D'}^f$} (6.2,-.25);
\path (5,.2) -- node[fill=white, inner ysep=0pt] {fronth.\ tx} (6.2,.2);
\draw[latex-latex] (6.2,-.25) -- node[fill=white, inner ysep=0pt] {$\delta_{D'}^e$} (7.3,-.25);
\path (6.2,.2) -- node[fill=white, inner ysep=0pt] {edge tx} (7.3,.2);
\draw[decorate,decoration=brace] (5,.4) -- node[above=1pt] {Delivery of $D'$} (7.3,.4);

\end{tikzpicture}
\caption{Delivery consists of a fronthaul transmission and an edge transmission phase. The length of these phases depends on the files that are requested.}
\label{fig:timeline}
\end{figure}

Each slot of the \emph{delivery phase} consists of two subsequent subslots (Fig.~\ref{fig:timeline}). In the \emph{first subslot}, the cloud sends information on the requested files to the ENs on the fronthaul links, while in the \emph{second subslot} the ENs use the shared wireless channel to transmit to the users.
To elaborate, in the first subslot, the cloud sends a message $U_m$ to the EN $m$ on the fronthaul as a function of the demand vector, the files and the CSI
\begin{equation} \label{eq:pif}
U_m = \pi^f_m(D,\mathcal{F},\mathcal{H}).
\end{equation}
The first subslot has $n^f_D$ symbols, where we make explicit the dependence on the vector $D$, and the entropy of message $U_m$ must be bounded as $H(U_m) \leq C^{f}n^{f}_D$ in order to satisfy the fronthaul capacity constraints. We call $\pi^f=(\pi^f_1,\dots,\pi^f_M)$ the fronthaul policy.
In the second subslot, the ENs transmit a codeword $X_m$, of $n^e_D$ symbols, on the wireless channel as a function of the users' demand $D$, the cache content $S_m=\{S_{m,j}\}_{j \in [1:J]}$ of EN $m$, the fronthaul message $U_m$ to EN $m$ and the global CSI $\mathcal{H}$:
\begin{equation}
X_m = \pi^e_m(D, S_m, U_m, \mathcal{H}).
\end{equation}
We call $\pi^e=(\pi^e_1,\dots,\pi^e_M)$ the edge transmission policy. After receiving $Y_k$ in (\ref{eq:rx_signal}), user $k$ decodes the requested file as
\begin{equation}
\hat F_{D_k} = \pi^d_k(Y_k, D, \mathcal{H}),
\end{equation}
and we let $\pi^d=(\pi^d_1,\dots,\pi^d_K)$ denote the decoding policy. The error probability of a policy $\pi=(\pi^c, \pi^f, \pi^e, \pi^d)$ is defined as the worst-case error probability across requests and users
\begin{equation}
P_e = \max_{D} \max_{k\in[1:K]} P(\hat F_{k} \neq F_{D_k}).
\end{equation} A sequence of policies, parametrized by $L$ and $P$, is defined as \emph{feasible} if it satisfies the limit $\lim_{P\to\infty}\lim_{L\to\infty} P_e=0$.

\subsection{Problem Statement} \label{ssec:problem}
For any sequence of feasible policies $\pi$ parametrized by $L$ and $P$, we now define the high-SNR delivery time metric for each demand vector $D$. To this end, we introduce the normalized durations of the first and second subslots in the given transmission interval as
\begin{equation}  \label{eq:NDT_components1}
\delta^x_D = \lim_{P\to\infty} \lim_{L\to\infty} \frac{n^x_D}{L/\log P},
\end{equation}
where $x=f$ for the first subslot (fronthaul transmission) and $x=e$ for the second subslot (edge transmission). In (\ref{eq:NDT_components1}), the subslot durations are normalized by the high-SNR delivery time of a reference system in which each user is served on an interference-free dedicated channel by an EN, namely $L/\log P$. Note, in fact, that an interference-free channel has a high-SNR capacity of $\log P$ (see also \cite{sengupta2016cloud} for additional discussion). We refer to $\delta^f_D$ and $\delta^e_D$ as the \emph{fronthaul and edge NDTs}, respectively, for request $D$. The overall NDT for request $D$ is hence given by $\delta_D = \delta^e_D + \delta^f_D$.

We are interested in characterizing the region $\Delta^{*}(\mu,r)$ of all achievable NDT tuples $\delta = (\delta_D)_{D\in\mathcal{D}}$ under the per-EN capacity constraint~\eqref{eq:capacityconstraint}, which we refer to as \emph{NDT region}. We impose that the same NDT be achieved $\delta_D$ be achieved for all permutations of the vector $D$. This allows us to obtain a characterization that depends only on the subset of files that are requested. Henceforth, with a slight abuse of notation, $D$ represents a subset of $[1:J]$. As a result, the NDT region is contained in the positive orthant of $\mathbb{R}^{J\choose K}$.

To study the NDT region $\Delta^*(\mu, r)$ it is convenient to analyze also the region $\Delta^*(\bm{\mu},r)$ of all NDT tuples that are achievable with a given cache partition matrix $\bm{\mu}$ in (\ref{eq:mumatrix}).
By definition, we have
 \begin{equation} \label{eq:ndtasunion}
 \Delta^*(\mu,r) = \bigcup_{\bm{\mu}:\mathrm{ } \sum_{j=1}^J \mu_{m,j} \leq\mu, \mathrm{ }\forall m} \Delta^*(\bm{\mu}, r).
\end{equation}

As a first important observation is summarized in the following lemma.
\begin{lemma} \label{lem:convex}
	The NDT regions $\Delta^*(\mu,r)$ and $\Delta^*(\bm{\mu},r)$ are convex.
\end{lemma}
\emph{Sketch of proof}: The lemma follows from the fact that, given two NDT tuples $\delta$ and $\delta'$ that are achievable with cache partitions $\bm{\mu}$ and $\bm{\mu'}$, respectively, that satisfy~\eqref{eq:capacityconstraint}, the NDT tuple $\alpha \delta+(1-\alpha)\delta'$ (with entry-wise sum), for any $0 \leq \alpha \leq 1$, can also be achieved using cache sharing and file splitting for a cache partition $\alpha\bm{\mu}+(1-\alpha)\bm{\mu'}$ which satisfies~\eqref{eq:capacityconstraint}. A detailed proof can be found in Appendix~\ref{app:convex}. \qed

We finally note that the minimum NDT introduced in \cite{sengupta2016cloud}, corresponds to the minimum value $\delta$ in the NDT region $\Delta^*(\bm{\mu},r)$, with equal cache partition $\mu_{m,j}=\mu$, such that the equality $\delta=\delta_D$ holds for all request subsets $D$. In the rest of this paper, we focus on the special case $K=M=2$ and we write $\delta_D=\delta_{i,j}$ for any request subset $D=\{i, j\}$.

\begin{figure}

\begin{subfigure}{\linewidth}
\centering
\begin{tikzpicture}[font=\scriptsize]

  \node (g1) at (-1.5,.3) {$G_1:$};
  \node (g2) at (-1.5,-.3) {$G_2:$};
  \node[right=0mm of g1, minimum width=15mm,draw] (g1bits) {$\nu L$ bits};
  \node[right=0mm of g2, minimum width=15mm,draw] (g2bits) {$\nu L$ bits};

  \node[fit=(g1) (g2) (g1bits) (g2bits), draw] (cloud) {};
  \node[above=0mm of cloud] {Cloud};

  \begin{scope}[xshift=33mm,yshift=8mm]
  \node[
    ] (muaa) at (-0.6,0) {$G_1$};

  \node[draw,fit=(muaa)] (ENa) {};
  \node[above=0mm of ENa] {EN$_1$};

  \end{scope}

  \begin{scope}[xshift=33mm,yshift=-8mm]
  \node[
    ] (muba) at (-0.6,0) {$G_2$};

  \node[draw,fit=(muba)] (ENb) {};
  \node[above=0mm of ENb] {EN$_2$};

  \end{scope}

  \draw[-latex,very thick] (cloud.east) -- node[above=5mm] {$\delta^f=\nu/r$} (ENa.west);
  \draw[-latex,very thick] (cloud.east) -- node[below,sloped] {} (ENb.west);

\end{tikzpicture}
\caption{Hard-transfer fronthauling (HT).}
\label{fig:existingA}
\end{subfigure}

\medskip
\begin{subfigure}{\linewidth}
\centering
\begin{tikzpicture}[font=\scriptsize]

\begin{scope}[xshift=0mm,yshift=10mm]
  \node (g1) at (0,0) {$G_1:$};
  \node[below=1mm of g1] (g2) {$G_2:$};
  \node[right=0mm of g1, minimum width=15mm,draw] (g1bits) {$\nu L$ bits};
  \node[right=0mm of g2, minimum width=15mm,draw] (g2bits) {$\nu L$ bits};
  \node[draw,fit=(g1) (g2) (g1bits) (g2bits)] (ENa) {};
\node[above=0mm of ENa] {EN$_1$};
\end{scope}

\begin{scope}[xshift=0mm,yshift=-10mm]
  \node (g1) at (0,0) {$G_1:$};
  \node[below=1mm of g1] (g2) {$G_2:$};
  \node[right=0mm of g1, minimum width=15mm,draw] (g1bits) {$\nu L$ bits};
  \node[right=0mm of g2, minimum width=15mm,draw] (g2bits) {$\nu L$ bits};
  \node[draw,fit=(g1) (g2) (g1bits) (g2bits)] (ENb) {};
\node[above=0mm of ENb] {EN$_2$};
\end{scope}

\node[ ] (useraa) at (4.2,0.7) {$G_1$};
\node[draw,fit=(useraa)] (usera) {};
\node[above=0mm of usera] {User 1};

\node[ ] (userbb) at (4.2,-1.25) {$G_2$};
\node[draw,fit=(userbb)] (userb) {};
\node[above=0mm of userb] {User 2};

\draw[-latex,very thick] (ENa.east) -- node[above=3mm] {$\delta^e=\nu$} (usera.west);
\draw[-latex,very thick] (ENa.east) -- (userb.west);
\draw[-latex,very thick] (ENb.east) -- (usera.west);
\draw[-latex,very thick] (ENb.east) -- (userb.west);

\end{tikzpicture}
\caption{Zero-forcing beamforming (ZF).}
\label{fig:existingB}
\end{subfigure}

\medskip
\begin{subfigure}{\linewidth}
\centering
\begin{tikzpicture}[font=\scriptsize]

  \node (g1) at (-1.5,.3) {$G_1:$};
  \node (g2) at (-1.5,-.3) {$G_2:$};
  \node[right=0mm of g1, minimum width=15mm,draw] (g1bits) {$\nu L$ bits};
  \node[right=0mm of g2, minimum width=15mm,draw] (g2bits) {$\nu L$ bits};

  \node[fit=(g1) (g2) (g1bits) (g2bits), draw] (cloud) {};
  \node[above=0mm of cloud] {Cloud};

  \begin{scope}[xshift=33mm,yshift=8mm]
  \node[minimum size=4mm
    ] (muaa) at (-0.6,0) {};

  \node[draw,fit=(muaa)] (ENa) {};
  \node[above=0mm of ENa] {EN$_1$};

  \end{scope}

  \begin{scope}[xshift=33mm,yshift=-8mm]
  \node[
     minimum size=4mm
    ] (muba) at (-0.6,0) {};

  \node[draw,fit=(muba)] (ENb) {};
  \node[above=0mm of ENb] {EN$_2$};

  \end{scope}

\node[ ] (useraa) at (5.2,.8) {$G_1$};
\node[draw,fit=(useraa)] (usera) {};
\node[above=0mm of usera] {User 1};

\node[ ] (userbb) at (5.2,-0.8) {$G_2$};
\node[draw,fit=(userbb)] (userb) {};
\node[above=0mm of userb] {User 2};

  \draw[-latex,very thick] (cloud.east) -- node[above=4mm] {$\delta^f=\nu/r$} (ENa.west);
  \draw[-latex,very thick] (cloud.east) -- node[below,sloped] {} (ENb.west);

\draw[-latex,very thick] (ENa.east) -- node[above=3mm] {$\delta^e=\nu$} (usera.west);
\draw[-latex,very thick] (ENa.east) -- (userb.west);
\draw[-latex,very thick] (ENb.east) -- (usera.west);
\draw[-latex,very thick] (ENb.east) -- (userb.west);

\end{tikzpicture}
\caption{Soft-transfer fronthauling with zero-forcing beamforming (ST+ZF).}
\label{fig:existingC}
\end{subfigure}

\medskip
\begin{subfigure}{\linewidth}
\centering
\begin{tikzpicture}[font=\scriptsize]

\begin{scope}[xshift=0mm,yshift=10mm]
  \node (g1) at (0,0) {$G_{1,1}:$};
  \node[below=1mm of g1] (g2) {$G_{1,2}:$};
  \node[right=0mm of g1, minimum width=15mm,draw] (g1bits) {$\nu L$ bits};
  \node[right=0mm of g2, minimum width=15mm,draw] (g2bits) {$\nu L$ bits};
  \node[draw,fit=(g1) (g2) (g1bits) (g2bits)] (ENa) {};
\node[above=0mm of ENa] {EN$_1$};
\end{scope}

\begin{scope}[xshift=0mm,yshift=-10mm]
  \node (g1) at (0,0) {$G_{2,1}:$};
  \node[below=1mm of g1] (g2) {$G_{2,2}:$};
  \node[right=0mm of g1, minimum width=15mm,draw] (g1bits) {$\nu L$ bits};
  \node[right=0mm of g2, minimum width=15mm,draw] (g2bits) {$\nu L$ bits};
  \node[draw,fit=(g1) (g2) (g1bits) (g2bits)] (ENb) {};
\node[above=0mm of ENb] {EN$_2$};
\end{scope}

\node[ ] (useraa) at (5.2,0.7) {$G_{1,1}$, $G_{2,1}$};
\node[draw,fit=(useraa)] (usera) {};
\node[above=0mm of usera] {User 1};

\node[ ] (userbb) at (5.2,-1.25) {$G_{1,2}$, $G_{2,2}$};
\node[draw,fit=(userbb)] (userb) {};
\node[above=0mm of userb] {User 2};

\draw[-latex,very thick] (ENa.east) -- node[above=3mm] {$\delta^e=3\nu$} (usera.west);
\draw[-latex,very thick] (ENa.east) -- (userb.west);
\draw[-latex,very thick] (ENb.east) -- (usera.west);
\draw[-latex,very thick] (ENb.east) -- (userb.west);

\end{tikzpicture}
\caption{X-channel interference alignment (X-IA).}
\label{fig:existingD}
\end{subfigure}
\caption{Illustration of constituent delivery strategies.}
\label{fig:existing}
\end{figure}

\begin{figure*}

\hfill
\begin{subfigure}{.45\linewidth}
\centering
\begin{tikzpicture}[font=\scriptsize,yscale=0.7]

\draw (0,2) rectangle (8,3);
\node[anchor=east] at (0,2.6) {$F_i$};
\draw (1.8,2) -- (1.8,3);
\draw (3.2,2) -- (3.2,3);
\draw[latex-latex] (0,1.8) -- node[below] {$\mu_i$} (1.8,1.8);
\draw[latex-latex] (0,2.5) -- node[fill=white,align=center,inner sep=0pt] {Cache\\ in EN$_1$} (1.8,2.5);
\draw[latex-latex] (1.8,2.5) -- node[fill=white,align=center,inner sep=0pt] {HT to\\ EN$_1$} (3.2,2.5);
\draw (6.2,2) -- (6.2,3);
\draw (4.8,2) -- (4.8,3);
\draw[latex-latex] (6.2,1.8) -- node[below,pos=0.5] {$\mu_i$} (8,1.8);
\draw[latex-latex] (6.2,2.5) -- node[fill=white,align=center,inner sep=0pt] {Cache\\ in EN$_2$} (8,2.5);
\draw[latex-latex] (4.8,2.5) -- node[fill=white,align=center,inner sep=0pt] {HT to\\ EN$_2$} (6.2,2.5);

\draw (0,0) rectangle (8,1);
\node[anchor=east] at (0,0.5) {$F_j$};
\draw (3.2,0) -- (3.2,1);
\draw[latex-latex] (0,1.2) -- node[above] {$\mu_j$} (3.2,1.2);
\draw[latex-latex] (0,0.5) -- node[fill=white] {Cache in EN$_1$} (3.2,0.5);
\draw (4.8,0) -- (4.8,1);
\draw[latex-latex] (4.8,1.2) -- node[above,pos=0.5] {$\mu_j$} (8,1.2);
\draw[latex-latex] (4.8,0.5) -- node[fill=white] {Cache in EN$_2$} (8,0.5);

\draw[rounded corners,dashed] (0.2,-0.2) -- (0.2,4) -- (7.8,4) -- (7.8,-0.2) -- (5,-0.2) -- (5,3.4) -- (3,3.4) -- (3,-.2) -- (0.2,-0.2);
\node at (4,3.7) {X-IA};

\draw[rounded corners,dashed] (3.4,-.6) rectangle (4.6,3.2);
\node at (4,-.4) {ST+ZF};

\end{tikzpicture}
\caption{Achievable strategy for $r\leq 1$, $\mu_i<1/2$ and $\mu_j<1/2$.}
\label{fig:ach1}
\end{subfigure}
\hfill
\begin{subfigure}{.45\linewidth}
\centering
\begin{tikzpicture}[font=\scriptsize,yscale=0.7]

\draw (0,2) rectangle (8,3);
\node[anchor=east] at (0,2.6) {$F_i$};
\draw (2.4,2) -- (2.4,3);
\draw[latex-latex] (0,1.8) -- node[below] {$\mu_i$} (2.4,1.8);
\draw[latex-latex] (0,2.5) -- node[fill=white] {Cache in EN$_1$} (2.4,2.5);
\draw (4,2) -- (4,3);
\draw[latex-latex] (2.4,2.5) -- node[fill=white,align=center,inner sep=0pt] {HT to\\ EN$_1$} (4,2.5);
\draw (5.6,2) -- (5.6,3);
\draw[latex-latex] (4,2.5) -- node[fill=white,align=center,inner sep=0pt] {HT to\\ EN$_2$} (5.6,2.5);
\draw[latex-latex] (5.6,1.8) -- node[below,pos=0.5] {$\mu_i$} (8,1.8);
\draw[latex-latex] (5.6,2.5) -- node[fill=white] {Cache in EN$_2$} (8,2.5);

\draw (0,0) rectangle (8,1);
\node[anchor=east] at (0,0.5) {$F_j$};
\draw (1.8,0) -- (1.8,1);
\draw[latex-latex] (0,1.2) -- node[above] {$\mu_j$} (6.2,1.2);
\draw[latex-latex] (0,0.5) -- node[fill=white,align=center,inner sep=0pt] {Cache\\ in EN$_1$} (1.8,0.5);
\draw[latex-latex] (1.8,0.5) -- node[fill=white,align=center] {Cache in EN$_1$ and in EN$_2$} (6.2,0.5);
\draw (6.2,0) -- (6.2,1);
\draw[latex-latex] (1.8,-0.2) -- node[below] {$\mu_j$} (8,-0.2);
\draw[latex-latex] (6.2,0.5) -- node[fill=white,align=center,inner sep=0pt] {Cache\\ in EN$_2$} (8,0.5);

\draw[rounded corners,dashed] (0.2,-0.8) -- (0.2,3.8) -- (7.8,3.8) -- (7.8,-0.8) -- (4.2,-0.8) -- (4.2,3.2) -- (3.8,3.2) -- (3.8,-.8) -- (0.2,-0.8);
\node at (4,3.5) {X-IA};

\end{tikzpicture}
\caption{Achievable strategy for $r\leq 1$ and $\mu_i\leq 1/2\leq\mu_j$.}
\label{fig:ach2}
\end{subfigure}
\hfill{}

\vspace{5mm}

\hfill
\begin{subfigure}[b]{.45\linewidth}
\centering
\begin{tikzpicture}[font=\scriptsize,yscale=0.7]

\draw (0,2) rectangle (8,3);
\node[anchor=east] at (0,2.6) {$F_i$};
\draw (2.4,2) -- (2.4,3);
\draw[latex-latex] (0,3.2) -- node[above,pos=.6] {$\mu_i$} (5.6,3.2);
\draw[latex-latex] (0,2.5) -- node[fill=white] {Cache in EN$_1$} (2.4,2.5);
\draw[latex-latex] (2.4,2.5) -- node[fill=white,align=center,inner sep=0pt] {Cache in\\ EN$_1$ and in EN$_2$} (5.6,2.5);
\draw (5.6,2) -- (5.6,3);
\draw[latex-latex] (2.4,1.8) -- node[below,pos=0.4] {$\mu_i$} (8,1.8);
\draw[latex-latex] (5.6,2.5) -- node[fill=white] {Cache in EN$_2$} (8,2.5);

\draw (0,0) rectangle (8,1);
\node[anchor=east] at (0,0.5) {$F_j$};
\draw (1.8,0) -- (1.8,1);
\draw[latex-latex] (0,1.2) -- node[above,pos=.6] {$\mu_j$} (6.2,1.2);
\draw[latex-latex] (0,0.5) -- node[fill=white,align=center,inner sep=0pt] {Cache\\ in EN$_1$} (1.8,0.5);
\draw[latex-latex] (1.8,0.5) -- node[fill=white,align=center,inner sep=0pt] {Cache in\\ EN$_1$ and in EN$_2$} (6.2,0.5);
\draw (6.2,0) -- (6.2,1);
\draw[latex-latex] (1.8,-0.2) -- node[below] {$\mu_j$} (8,-0.2);
\draw[latex-latex] (6.2,0.5) -- node[fill=white,align=center,inner sep=0pt] {Cache\\ in EN$_2$} (8,0.5);

\draw[rounded corners,dashed] (0.2,-0.4) -- (0.2,4.6) -- (7.8,4.6) -- (7.8,-0.4) -- (5.8,-0.4) -- (5.8,4) -- (2.2,4) -- (2.2,-.4) -- (0.2,-0.4);
\node at (4,4.4) {X-IA};

\draw[rounded corners,dashed] (2.6,-.8) rectangle (5.4,3.8);
\node at (4,-.6) {ZF};

\end{tikzpicture}
\caption{Achievable strategy for $r\leq 1$ and $\mu_i> 1/2$.}
\label{fig:ach3}
\end{subfigure}
\hfill
\begin{subfigure}[b]{.45\linewidth}
\centering
\begin{tikzpicture}[font=\scriptsize,yscale=0.7]

\draw (0,2) rectangle (8,3);
\node[anchor=east] at (0,2.6) {$F_i$};
\draw (4.3,2) -- (4.3,3);
\draw[latex-latex] (0,1.8) -- node[below,pos=.5] {$\mu_i$} (4.3,1.8);
\draw[latex-latex] (0,2.5) -- node[fill=white,align=center,inner sep=0pt] {Cache in\\ EN$_1$ and in EN$_2$} (4.3,2.5);

\draw (0,0) rectangle (8,1);
\node[anchor=east] at (0,0.5) {$F_j$};
\draw (5.6,0) -- (5.6,1);
\draw[latex-latex] (0,1.2) -- node[above,pos=.5] {$\mu_j$} (5.6,1.2);
\draw[latex-latex] (0,0.5) -- node[fill=white,align=center,inner sep=0pt] {Cache in\\ EN$_1$ and in EN$_2$} (5.6,0.5);

\draw[rounded corners,dashed] (0.2,-.2) rectangle (4.1,3.7);
\node at (2.05,3.4) {ZF};

\draw[rounded corners,dashed] (4.5,-.2) rectangle (7.8,3.7);
\node at (6.15,3.4) {ST+ZF};

\end{tikzpicture}
\caption{Achievable strategy for $r>1$.}
\label{fig:ach4}
\end{subfigure}
\hfill{}

\caption{Achievable strategies}
\end{figure*}

\section{Preliminaries} \label{ssec:existing}
Here we review delivery strategies, see Fig.~\ref{fig:existing}, for the fronthaul and edge channels from \cite{sengupta2016cloud}, which will be used as ingredients in the next section to propose a more general caching and delivery policy. (1) \emph{Hard-transfer fronthauling} (HT, Fig.~\ref{fig:existingA}): As shows, via the HT fronthaul delivery strategy, the cloud delivers a fraction $\nu$ of one of the requested files, say $G_1$, to EN 1 and a fraction of the other file $G_2$ to EN 2 on the respective fronthaul links. (2) \emph{Zero-forcing beamforming} (ZF, Fig.~\ref{fig:existingB}): If both ENs have both requested messages $G_1$ and $G_2$, or a fraction $\nu$ thereof, available in the respective caches, the edge delivery strategy of cooperative ZF beamforming can be carried out on this fraction to deliver $G_i$ to user $i$, yielding parallel interference-free channels to both users. (3) \emph{Soft-transfer fronthauling with zero-forcing beamforming} (ST+ZF, Fig.~\ref{fig:existingC}): With the fronthaul-edge delivery strategy, the cloud implements ZF beamforming and transmits the resulting baseband signals to the ENs in quantized form. The ENs simply forward the quantized signals over the shared wireless channel \cite{simeone2009downlink}. (4) \emph{X-channel interference alignment} (X-IA, Fig.~\ref{fig:existingD}): If the ENs cache different fractions $\nu$ of each requested file, the resulting channel model for the delivery for this fraction is an X-channel, for which interference alignment (IA) edge delivery strategies were presented in~\cite{maddah2008communication}.

\begin{lemma} \label{lem:existing}\cite{TandonSimeone16}
The following fronthaul and edge NDTs
are achievable using the delivery strategies summarized above.

\noindent{\em HT:} Let $G_1$ and $G_2$ be messages of $\nu L$ bits that are available in the cloud. HT requires the fronthaul NDT $\delta^f = \frac{\nu}{r}$ to transmit $G_1$ to EN$_1$ and $G_2$ to EN$_2$.

\noindent{\em ZF:} Let both ENs have messages $G_1$ and $G_2$ of $\nu L$ bits available. ZF requires the edge NDT $\delta^e = \nu$ to transmit $G_1$ to user $1$ and $G_2$ to user $2$.

\noindent{\em ST+ZF:} Let $G_1$ and $G_2$ be messages of $\nu L$ bits that are available in the cloud. ST+ZF requires the fronthaul and edge NDTs $\delta^f = \frac{\nu}{r}\quad\text{and}\quad\delta^e = \nu$ to transmit $G_1$ to user $1$ and $G_2$ to user $2$.

\noindent{\em X-IA:} Let $G_{i,1}$ and $G_{i,2}$ be messages of $\nu L$ bits that are available at EN$_i$, $i=1,2$. X-IA requires the edge NDT $\delta^e = 3\nu$ to transmit $G_{1,1}$ and $G_{2,1}$ to user $1$ and  $G_{1,2}$ and $G_{2,2}$ to user $2$.
\end{lemma}

\section{Achievable NDT Region} \label{sec:inner}
In this section, we present achievable strategies that yield an inner bound on the NDT region $\Delta^{*}(\mu,r)$. To this end, we consider policies with cache partitions $\bm{\mu}$ such that the two ENs cache the same number of bits for each file, i.e., $\mu_{1,j}=\mu_{2,j}=\mu_j$. As a result, each file $F_j$ is generally allocated a different cache fraction $\mu_j$ at the ENs. We will show in the next section that this restriction comes with no loss of optimality.
\begin{theorem} \label{th:inner}
An inner bound on the NDT region is given by the inclusion \begin{equation} \Delta^*(\mu,r) \supseteq\Delta^{(\mathsf{in})}(\mu,r) = \bigcup_{\substack{\bm{\mu}: \mathrm{ } \mu_{1,i}=\mu_{2,i}=\mu_i,\\ \sum_{j=1}^J\mu_{j}\leq\mu}} \Delta^{(\mathsf{in})}(\bm{\mu}, r),\end{equation}where the region
\begin{equation}
\Delta^{(\mathsf{in})}(\bm{\mu},r)= \left\{\delta_D \middle| \delta_{i,j}\geq\delta_{i,j}^{(\mathsf{in})}(\bm{\mu},r), \forall \{i,j\}\in\mathcal{D}\right\}
\end{equation}
is included in $\Delta^*(\bm{\mu},r)$, and we have
\begin{align} \label{eq:deltainner}
\delta_{i,j}^{(\mathsf{in})}(\bm{\mu},r) =
\begin{cases}
\delta_{i,j}^{(\mathsf{in},1)}(\bm{\mu},r),\quad &\text{if } r\leq 1, \mu_i<\frac{1}{2} \text{ and } \mu_j < \frac{1}{2}, \\
\delta_{i,j}^{(\mathsf{in},2)}(\bm{\mu},r),\quad &\text{if } r\leq 1 \text{ and }  \\
&\ \hspace{-6.6mm}(\mu_i\leq\frac{1}{2}\leq\mu_j, \text{ or } \mu_j\leq\frac{1}{2}\leq\mu_i), \\
\delta_{i,j}^{(\mathsf{in},3)}(\bm{\mu},r),\quad &\text{if } r\leq 1, \mu_i>\frac{1}{2} \text{ and } \mu_j>\frac{1}{2}, \\
\delta_{i,j}^{(\mathsf{in},4)}(\bm{\mu},r),\quad &\text{if } r>1,
\end{cases}
\end{align}
with the definitions
\begin{align}
\delta_{i,j}^{(\mathsf{in},1)}(\bm{\mu},r) =&\ 1 + \frac{1}{r} - \left(\frac{1}{r} - 1\right)\max\{\mu_i, \mu_j\} \notag \\
& - \frac{1}{r}\min\{\mu_i, \mu_j\}, \\
\delta_{i,j}^{(\mathsf{in},2)}(\bm{\mu},r) =&\ \frac{3}{2} + \frac{1}{r}\left(\frac{1}{2}-\min\{\mu_i,\mu_j\}\right), \\
\delta_{i,j}^{(\mathsf{in},3)}(\bm{\mu},r) =&\ 2 - \min\{\mu_i,\mu_j\}, \\
\delta_{i,j}^{(\mathsf{in},4)}(\bm{\mu},r) =&\ 1 + \frac{1}{r} - \frac{1}{r}\min\{\mu_i,\mu_j\}.
\end{align}
\end{theorem}
In the remainder of this section, we present the achievable strategies that yield the inner bound in the previous theorem at an intuitive level. The detailed proof of Theorem~\ref{th:inner} is given in Appendix~\ref{app:inner}.
To this end, we will present two different caching policies for the cases $r\leq 1$ and $r>1$. Note that the caching policy cannot depend on the demand $D=\{i,j\}$, unlike the delivery policy. In the following, we set $\mu_i\leq \mu_j$ without loss of generality.

\noindent \emph{Caching policy for $r\leq 1$}: As seen in Figures~\ref{fig:ach1}--\ref{fig:ach3}, each file $F_i$ is cached so that the bits indexed by $1,\dots,\mu_i L$ are stored in EN$_1$ and bits $(1-\mu_i)L,\dots, L$ are stored in EN$_2$, {i.e.},\  we minimize the overlap in the cached content in EN$_1$ and EN$_2$ by storing the first part of the file in EN$_1$ and the last part of the file in EN$_2$. If $\mu_i\geq 1/2$ some overlap will occur and some bits will be stored in both ENs.

\noindent \emph{Caching policy for $r>1$}: As seen in Figure~\ref{fig:ach4}, each file $F_i$ is cached so that bits $1,\dots,\mu_i L$ in both EN$_1$ and EN$_2$, {i.e.},\  we cache only the first part of the file and we maximize the overlap between the content that is cached in the ENs.

\noindent \emph{Delivery strategy for $\mu_j<1/2$ and $r\leq 1$}: This case is illustrated in Figure~\ref{fig:ach1} and achieves $\delta_{i,j}^{(\mathsf{in},1)}(\bm{\mu},r)$. The delivery proceeds in three phases: a) we use HT to deliver bits $\mu_iL,\dots\mu_jL$ and $(1-\mu_j)L,\dots,(1-\mu_i)L$ of file $F_i$ to EN$_1$ and EN$_2$, respectively; b) we use X-IA to transmit bits $1,\dots,\mu_jL$ and $(1-\mu_j)L,\dots,L$ of the files from the ENs to the users; c) we use ST+ZF to transmit bits $\mu_jL,\dots,(1-\mu_j)L$ directly from the cloud. Note that the strategy in~\cite{sengupta2016cloud} does not require step a). In fact, interestingly, the optimal policy in~\cite{sengupta2016cloud} did not make any use of HT fronthauling. The optimality results presented in the next section demonstrate that, instead, when $\mu_i\neq\mu_j$, the joint use of both HT and ST are instrumental in achieving the optimal NDT performance.

\noindent \emph{Delivery strategy for $\mu_i\leq 1/2\leq\mu_j$ and $r\leq 1$}: This case is illustrated in Figure~\ref{fig:ach2} and achieves $\delta_{i,j}^{(\mathsf{in},2)}(\bm{\mu},r)$. The delivery proceeds in two phases: a) we use HT to deliver bits $\mu_iL,\dots L/2$ and $L/2,\dots,(1-\mu_i)L$ of file $F_i$ to EN$_1$ and EN$_2$, respectively; and b) we use X-IA to deliver both complete files from the ENs to the users. We remark that this scenario is not relevant for the special case from~\cite{sengupta2016cloud}. We emphasize the important role of HT for deriving an achievable strategy, which is used here but not in~\cite{sengupta2016cloud},

\noindent \emph{Delivery strategy for $\mu_i>1/2$ and $r\leq 1$}: This case is illustrated in Figure~\ref{fig:ach3} and achieves $\delta_{i,j}^{(\mathsf{in},3)}(\bm{\mu},r)$. The delivery proceeds in two phases: a)  we use X-IA to deliver bits $1,\dots,(1-\mu_i)L$ and $\mu_iL,\dots,L$ of the files from the ENs to the users; and b) we use ZF to transmit bits $(1-\mu_i)L,\dots,\mu_iL$ from the ENs to the users. Note that this strategy does not make use of the fronthaul during the delivery.

\noindent \emph{Delivery strategy for $r > 1$}: The first $\min\{\mu_i, \mu_j\}L$ bits of both files, which are stored at both ENs, are delivered using ZF. The remaining $(1-\min\{\mu_i,\mu_j\})L$ bits are delivered using ST+ZF. The strategy is illustrated in Fig.~\ref{fig:ach4} and achieves $\delta_{i,j}^{(\mathsf{in},4)}(\bm{\mu},r)$.

\section{Characterization of the NDT Region} \label{sec:outer}
In this section, we present an outer bound on the NDT region $\Delta^*(\mu,r)$ and we prove that the inner bound from the previous section is in fact tight, hence characterizing the NDT region. The first result of this section provides an outer bound on the achievable NDT tuple region for a fixed, and generic, cache partition $\bm{\mu}$.
\begin{theorem} \label{th:outer} For any cache partition $\bm{\mu}$, we have the outer bound $\Delta^*(\bm{\mu},r)\subseteq\Delta^{(\mathsf{out})}(\bm{\mu}, r)$, where
\begin{equation}
\Delta^{(\mathsf{out})}(\bm{\mu}, r) = \left\{ \delta\ \middle|\ \delta_{i,j} \geq \delta_{i,j}^{(\mathsf{out})}(\bm{\mu},r), \forall \{i,j\}\in\mathcal{D}  \right\},
\end{equation}
with
\begin{align} \label{eq:deltaouter}
\delta_{i,j}^{(\mathsf{out})}(\bm{\mu},r) =
\begin{cases}
\max_{\ell=1,\dots,3}\left\{\delta_{i,j}^{(\mathsf{out},\ell)}(\bm{\mu},r)\right\},\quad &\text{if } r\leq 1, \\
\delta_{i,j}^{(\mathsf{out},4)}(\bm{\mu},r),\quad &\text{if } r>1,
\end{cases}
\end{align}
and the definitions
\begin{multline}
\delta_{i,j}^{(\mathsf{out},1)}(\bm{\mu},r) = 1 + \frac{1}{r} - \min\left\{\mu_{1,i}, \mu_{2,i}, \mu_{1,j}, \mu_{2,j}\right\} \\
-\frac{1}{2}\left(\frac{1}{r}-1\right)\left(  \mu_{1,i} +  \mu_{2,i} + \mu_{1,j} + \mu_{2,j}  \right), \label{eq:thconverseA}
\end{multline}
\begin{multline}
\delta_{i,j}^{(\mathsf{out},2)}(\bm{\mu},r) = \frac{3}{2} + \frac{1}{2r} - \min\left\{\mu_{1,i}, \mu_{2,i}, \mu_{1,j}, \mu_{2,j}\right\} \\
- \frac{1}{2}\left(\frac{1}{r}-1\right) \min\{\mu_{1,i} + \mu_{2,i}, \mu_{1,j} + \mu_{2,j}\}, \label{eq:thconverseB}
\end{multline}
\begin{equation}
\delta_{i,j}^{(\mathsf{out},3)}(\bm{\mu},r) = 2 - \min\left\{\mu_{1,i}, \mu_{2,i}, \mu_{1,j}, \mu_{2,j}\right\}, \label{eq:thconverseC}
\end{equation}
\begin{equation}
\delta_{i,j}^{(\mathsf{out},4)}(\bm{\mu},r) = 1 + \frac{1}{r} - \frac{1}{r}\min\left\{\mu_{1,i}, \mu_{2,i}, \mu_{1,j}, \mu_{2,j}\right\}. \label{eq:thconverseD}
\end{equation}

\end{theorem}
The proof of Theorem~\ref{th:outer} is given in Appendix~\ref{app:outer}.
Using the outer bound in the previous theorem, we show that the inner bound from the previous section is tight. This result implies that, in order to exhaust the NDT region, it is sufficient to consider cache partitions in which $\mu_{1,j}=\mu_{2,j}$ for all files $j \in [1:J].$
\begin{theorem} \label{th:region} The NDT region is given as $\Delta^*(\mu,r) = \Delta^{(\mathsf{in})}(\mu,r)$.
\end{theorem}
The proof of Theorem~\ref{th:region} is given in Appendix~\ref{app:region}.

\section{Numerical Example} \label{sec:eval}
Consider a set-up in which the set of popular files is partitioned into two disjoint classes as $\mathcal{F}=\mathcal{F}_{1}\cup\mathcal{F}_{2}$, where class $\mathcal{F}_{i}$ has $J_i$ files. We first illustrate the NDT region derived above and we then discuss how this can be used to obtain optimal trade-offs among the individual latencies of different files under a popularity distribution.

We first illustrate a slice of the NDT region in which we impose that the same NDT $\delta_{(i),(j)}$ be achieved for all subsets $D$ for which one file is in the class $\mathcal{F}_i$ and the other in $\mathcal{F}_j$. Recall that we assume that the requested files are distinct (even if requested from the same class). The considered slice of the NDT region is three-dimensional with axes given by $\delta_{(1),(1)}$, $\delta_{(2),(2)}$ and $\delta_{(1),(2)}$. To further reduce the dimensionality, we let $\delta_{(2),(2)}$ be arbitrary, so as to focus only on the plane $(\delta_{(1),(2)}, \delta_{(1),(1)})$. In order to evaluate the boundary of this slice of the NDT region, it can be argued that it is sufficient to consider cache partitions such as all files within the same class, which we denote as $\mu_{(1)}$ and $\mu_{(2)}$ for the files of class $1$ and $2$, respectively. With this choice, the cache capacity constraint~\eqref{eq:capacityconstraint} reduces to $J_1\mu_{(1)} + J_2\mu_{(2)} \leq \mu(J_1+J_2)$.

The slice of the NDT region at hand is illustrated in Figure~\ref{fig:ndtregion} for $r=1/5$, $J_1=J_2$ and various values of $\mu$. The figure also indicates the values of the cache allocations $(\mu_{(1)},\mu_{(2)})$ that are required to obtain various points on the boundary of the region as well as the delivery strategy that should be used at various segments of the boundary. As it can be seen, the slice of the NDT region is a polyhedron and each linear portion of the boundary corresponds to a different delivery strategy as indicated in the figure (see Sec. III for a correspondence between strategies and NDT tuples $\delta^{(\mathrm{in},\ell)}$).  For instance, it is seen that, for $\mu=3/8$, one has to use a different strategy depending on the operating point: as $\delta_{(1),(2)}$ increases, one needs to switch between the strategies that achieve the NDTs $\delta_{i,j}^{(\mathsf{in},1)}$ and $\delta_{i,j}^{(\mathsf{in},2)}$.

Finally, we consider individual latency constraints for different files under a given popularity profile. To this end, we let $a$ and $1-a$ denote the probabilities that a file is requested from class $\mathcal{F}_{1}$ and from class $\mathcal{F}_{2}$, respectively. We then have that the probability $p_{11}$ that two files from class $\mathcal{F}_{1}$ are selected is $p_{11}=a^2$; the probability $p_{12}$ that one file from each class is requested is $p_{12}=2a(1-a)$; and the probability $p_{22}$ that two files from class $\mathcal{F}_{2}$ are requested is $p_{22}=(1-a)^2$. The average latency for a file from a given class is:
\begin{align}
\bar\delta_{(1)} = \mathbb{E}[\delta_{(1)}] &= \frac{p_{11}}{p_{11}+p_{12}}\delta_{(1),(1)} +  \frac{p_{12}}{p_{11}+p_{12}}\delta_{(1),(2)}, \\
\bar\delta_{(2)} = \mathbb{E}[\delta_{(2)}] &= \frac{p_{22}}{p_{22}+p_{12}}\delta_{(1),(1)} +  \frac{p_{12}}{p_{22}+p_{12}}\delta_{(1),(2)}.
\end{align}
Note that $\bar\delta_{(i)}$ is the average latency for files of class $\mathcal{F}_{(i)}$ when averaged over the second requested file.

In Figure~\ref{fig:average}, we illustrate the optimal trade-off between the average NDTs $\bar\delta_{(1)}$ and $\bar\delta_{(2)}$ that arises from adjusting the cache allocations among the two classes. In the figure we have set $\mu=3/8$, $r=1/5$, $J_1=J_2$ and considered various values of $a$. The figure confirms that obtaining lower average delivery latencies for some files entails a larger average delivery latencies for other files due to the limited cache capacities.

\begin{figure}
\centering
\begin{tikzpicture}
\tikzstyle{point}=[draw,fill,circle,minimum size=3pt,inner sep=0pt]
\begin{axis}[
	xlabel=$\delta_{(1),(2)}$,
	ylabel=$\delta_{(1),(1)}$,
	xmin=0.3,
	xmax=5,
	ymin=0.3,
	ymax=5,
	font=\scriptsize,
]
\addplot[
  line width=.3mm,color=red, loosely dashed,
  ]
table[
  header=false,x index=0,y index=1,
  ]
{matlab/two_classes2ext.csv}
node [pos=0.12,rotate=90,anchor=north] {$\mu=1/4$};
\node[point,red] (pa) at (axis cs: 4,1.5) {};
\node[red,align=center] at (axis cs: 4.3,1.7) {$\left(\frac{1}{2}, 0\right)$};
\node[point,red] (pb) at (axis cs: 3.75,3.75) {};
\node[red,align=center] at (axis cs: 4.1,3.75) {$\left(\frac{1}{4}, \frac{1}{4}\right)$};
\addplot[
  line width=.3mm,color=green!50!black, loosely dashed,
  ]
table[
  header=false,x index=0,y index=1,
  ]
{matlab/two_classes3ext1.csv}
node [pos=0.15,rotate=90,anchor=north] {$\mu=3/8$};
\addplot[
  line width=.3mm,color=green!50!black, dotted,
  ]
table[
  header=false,x index=0,y index=1,
  ]
{matlab/two_classes3ext2.csv};
\node[point,green!50!black] (pc) at (axis cs: 4,1.25) {};
\node[green!50!black,align=center] at (axis cs: 3.7,1.3) {$\left(\frac{3}{4}, 0\right)$};
\node[point,green!50!black] (pd) at (axis cs: 2.75,1.5) {};
\node[green!50!black,align=center] at (axis cs: 3.05,1.7) {$\left(\frac{1}{2}, \frac{1}{4}\right)$};
\node[point,green!50!black] (pe) at (axis cs: 2.625,2.625) {};
\node[green!50!black,align=center] at (axis cs: 3, 2.625) {$\left(\frac{3}{8}, \frac{3}{8}\right)$};
\addplot[
  line width=.3mm,color=black, dotted,
  ]
table[
  header=false,x index=0,y index=1,
  ]
{matlab/two_classes4ext.csv}
node [pos=0.075,rotate=90,anchor=north] {$\mu=1/2$};
\node[point,black] (pf) at (axis cs: 4,1) {};
\node[black,align=center] at (axis cs: 4,0.8) {$\left(1, 0\right)$};
\node[point,black] (pg) at (axis cs: 1.5,1.5) {};
\node[black,align=center] at (axis cs: 1.85,1.7) {$\left(\frac{1}{2}, \frac{1}{2}\right)$};
\addplot[
  line width=.3mm,color=blue, solid,
  ]
table[
  header=false,x index=0,y index=1,
  ]
{matlab/two_classes6ext.csv}
node [pos=0.07,rotate=90,anchor=south] {$\mu=3/4$};
\node[point,blue] (ph) at (axis cs: 1.5,1) {};
\node[blue,align=center] at (axis cs: 1.5,0.8) {$\left(1,\frac{1}{2}\right)$};
\node[point,blue] (pi) at (axis cs: 1.25,1.25) {};
\node[blue,align=center] at (axis cs: 0.9,1.25) {$\left(\frac{3}{4}, \frac{3}{4}\right)$};
\end{axis}
\end{tikzpicture}
\vspace{-4mm}
\caption{Slice of the NDT region as a function of the fractional cache capacity $\mu$ ($r=1/5$, $J_1=J_2$). The labels indicate the cache allocations $(\mu_{(1)},\mu_{(2)})$ that are required at specific points. The line styles indicate the strategy to be used, with the dashed line corresponding to $\delta_{i,j}^{(\mathsf{in},1)}$, the dotted lines to $\delta_{i,j}^{(\mathsf{in},2)}$, and the solid line to $\delta_{i,j}^{(\mathsf{in},3)}$.} \label{fig:ndtregion}
\end{figure}

\begin{figure}
\centering
\begin{tikzpicture}
\tikzstyle{point}=[draw,fill,circle,minimum size=3pt,inner sep=0pt]
\begin{axis}[
	height=60mm,
	width=85mm,
	xlabel=$\bar\delta_{(1)}$,
	ylabel=$\bar\delta_{(2)}$,
	xmin=1,
	xmax=4.5,
	ymin=1,
	ymax=4,
	font=\scriptsize,
]

\addplot[
  line width=.3mm,color=green, dashed,
  ]
table[
  header=false,x index=0,y index=1,
  ]
{matlab/two_classes_average_1.csv};
\addlegendentry{$a=0.1$};

\addplot[
  line width=.3mm,color=blue, solid,
  ]
table[
  header=false,x index=0,y index=1,
  ]
{matlab/two_classes_average_5.csv};
\addlegendentry{$a=0.5$};

\addplot[
  line width=.3mm,color=red, dotted,
  ]
table[
  header=false,x index=0,y index=1,
  ]
{matlab/two_classes_average_9.csv};
\addlegendentry{$a=0.9$};

\end{axis}
\end{tikzpicture}
\vspace{-4mm}
\caption{Optimal trade-off between the average delivery latencies for the files of two classes for different popularity profiles defined by $a$ ($\mu=3/8$,$r=1/5$, $J_1=J_2$).}
\label{fig:average}
\end{figure}

\section{Conclusions}
This work characterized the set of delivery latencies supported by an F-RAN with two ENs and two users in the high SNR regime, when allowing for any cache partition across the files in a set of popular contents. Various aspects call for further investigation, including the explicit minimization of the average delivery latency as a function of the content popularity profile, the extension of the main results to any number of ENs and users (see~\cite{sengupta2016cloud} for the case of uniform file popularity) and the derivation of an extended NDT region in which the same contents may be requested by multiple users.

\bibliographystyle{IEEEtran}