text:
\section{Introduction}

	The notion of \emph{generalized surfaces} introduced by Young~\cite{young1937generalized-cur,young1942generalized-sur,young1942generalized-surII}, and known today as Young's measures, rests on the fundamental idea to consider functions as graphs.
	%

	Young realized that the weak convergence of the graphs of a sequence of functions carries
	substantially more information than the weak convergence of the functions themselves. In fact, this was the cornerstone leading to the
	following fundamental principle of Young measures:
	Let $(u_k)_{k\in\mathbb{N}} \subset \mathrm{L}^1(\Omega;\mathbb{R}^N)$ be a uniformly bounded sequence. Then (up to passing to a subsequence) there exists a family of probability measures $\{\nu_x\} \subset \mathrm{Prob}(\mathbb{R}^N)$, parameterized by $x$ in $\Omega$, such that
	\begin{equation}\label{eq:Y1}
	g(u_{k}) \rightharpoonup u_g \quad \text{for all $\mathrm{L}^1(\Omega)$}\quad \text{for some $g \in \mathrm{C}_c(\mathbb{R}^N)$}\,,
	\end{equation}
	where
	\begin{equation}\label{eq:Y2}
	u_g = \,\thetant-_{\Omega} g(z) \;\mathrm{d} \nu_x(z)\,.
	\end{equation}

	In the field of applications, the seminal work of M\"uller \cite{muller1999variational-mod} (among many others) adds up to a fair amount of applications in optimal design where the framework of
	Young measures plays a fundamental role in their development.
	However, as pointed out by Tartar~\cite{tartar1995beyond-young-me} and Pedregal~\cite{pedregal2006multi-scale-you}, Young measures have their own drawbacks and
	limitations. The \emph{first} important limitation is their incapability to keep track of concentration of mass (mass carried by the sequence may escape
	to infinity, in $\mathbb{R}^N$, while leaving the limit in~\eqref{eq:Y1} unchanged). To solve this issue, DiPerna \& Majda adapted Young's ideas using a compactification of $\mathbb{R}^N$. They
	extended classical Young measures to what today is known as generalized Young measures \cite{diperna1987oscillations-an} (see also \cite{alibert1997non-uniform-int}), which are capable of representing more general limits than~\eqref{eq:Y1}. More precisely, extending~\eqref{eq:Y1} to the representation of limits $\mu_f$ of the form
	\begin{equation}\label{eq:Y3}
	f(u_k) \overset{*}\rightharpoonup \mu_f \quad \text{in $\mathcal M(\Omega)$ \quad among integrands satisfying $\tilde f \in \mathrm{C}(\beta \mathbb{R}^N)$},
	\end{equation}
	with $\tilde f(z) \coloneqq f(z)/(1 + |z|)$, and where $\beta \mathbb{R}^N$ is the Stone--\v Cech compactification of $\mathbb{R}^N$. Here, we omit the representation
	formula for generalized Young measures as it is substantially more involved than~\eqref{eq:Y2}; we shall postpone this to the Appendix where for the convenience of the reader we give a
	brief sketch of the construction and properties of (the different notions of) Young measures. The \emph{second} main drawback of Young's
	construction is the failure to record patterns such as the direction or speed where oscillation (and/or concentration) occurs. This is easily illustrated by the following
	one dimensional example. Fix a positive real $\alpha$ and consider the purely oscillatory sequence
	\[
	u_k(x) = \sin(k^\alpha x), \qquad k = 1,2,\dots
	\]
	Clearly, the choice of $\alpha$ significantly changes the length-scale period at which oscillations occur as $k$ tends to $\infty$. However, regardless of the
	choice of $\alpha$, a change of variables argument shows the associated Young measure to this sequence is given by the family $\{\nu_x\}_{x \in \mathbb{R}}$ where $\nu_x = \nu_0$ for
	almost every $x  \in \Omega$, and where $\nu_0$ is the probability measure satisfying
	\[
	\dpr{\nu_0,g} = \frac{1}{2\pi}\int_0^{2\pi} g(\sin y) \;\mathrm{d} y \quad \text{for all $g \in \mathrm{C}_c(\mathbb{R})$}\,.
	\]
	Hence, there is a need to extend the notion of Young measure to one that incorporates the dependence of the parameter $\alpha$. To record this information, Pedregal
	considered the joint Young measure $\upsilon$ associated to the sequences of pairs
	\[
	w_k(x) \coloneqq (\area{\chi_k},u_k) : x \mapsto (\area{kx},u_k(x))\,,
	\]
	where $\area{x}$ denotes the equivalence class of the vector $x \in \mathbb{R}^d$ in the $d$-dimensional flat torus $Z \coloneqq \mathbb{R}^d/\mathbb{Z}^d$. {After performing a slicing argument, $\upsilon_x$ decomposes into} $\upsilon_x = \pi_\# \nu_x \otimes \nu_{x,\xi}$, where $\pi : Z \times \mathbb{R}^N \to Z$ is the canonical projection on the torus. {Pedregal} introduced the resulting family $\{\nu_{x,\xi}\}$ of probability measures on $\mathbb{R}^N$ as the associated \emph{two-scale Young measure}, in turn, designed to represent weak-limits of the form
	\begin{equation}
	g(w_k) \rightharpoonup U_g \quad \text{in $\mathrm{L}^1(\Omega)$}\quad \text{for integrands $g \in \mathrm{C}_c(Z \times \mathbb{R}^N)$}\,.
	\end{equation}
	Similarly to the case of generalized Young measures, a sketch of Pedregal's construction is further discussed in the Appendix.

	We are now in a position to give a rough description of the content and goals of this work.

	\subsection{Main results}

	The \emph{first} goal of this paper is to introduce \emph{generalized multi-scale} Young measures (\enquote{multi-scale*} throughout the text for short) in the following sense. We combine the ideas of
	DiPerna \& Majda with the approach from Pedregal to construct a {new type of Young measures}, capable of dealing with oscillation-concentration effects while also quantifying the length-scales where these phenomena occur.
	For simplicity we shall restrict only to two-scale* Young measures. Effectively, we introduce a measure-theoretic tool to represent weak-$*$ limits of the form
	\begin{equation}\label{eq:Y5}
	f(\chi_k,u_k) \overset{*}\rightharpoonup \mu_f \quad \text{in $\mathcal M(\Omega)$ \quad amongst integrands $f$ satisfying $\tilde f \in \mathrm{C}(Z \times \beta \mathbb{R}^N)$}\,,
	\end{equation}
	where $\tilde f(\xi,z) \coloneqq f(\xi,z)/(1 + |z|)$. Next, we give the rigorous definition
	and state some of the main properties of two-scale* Young measures.

	\begin{definition}[two-scale* Young measures]
		A four-tuple $\bm{\nu} = (\nu,\lambda,\rho,\nu^\infty)$ is called a two-scale* Young measure on $\Omega$ with values in $E$ (a finite dimensional euclidean space) provided that
		\begin{enumerate}[(i)]
			\item $\nu$ is a weak-$*$ measurable from $\cl \Omega \times Z$ into the set of probability measures over $E$ such that the map
			$(x,\xi) \mapsto \dpr{\nu_{x, \xi}, |\,\sbullet\,|}$ belongs to $\mathrm{L}^1(\cl\Omega \times Z)$,
			\item $\lambda$ is a bounded positive measure on $\cl \Omega$,
			\item $\rho$ is a weak-$*$ $\lambda$-measurable map from $\cl \Omega$ to the
			set of probability measures on the $d$-dimensional torus,
			\item $\nu^\infty$ is a weak-$*$ ($\lambda \otimes \rho_x)$-measurable from
			$\cl \Omega \times Z$ into the set of probability measures over $\mathbb{S}_E$ (the sphere of radius one in $E$).\footnote{The semi-product $\lambda \otimes \nu_x$, between a positive measure $\lambda \in \mathcal M(\Omega)$ and a weak-$*$ $\lambda$-measurable map $\nu : \Omega \to \mathrm{Prob}(K) : x \mapsto \nu_x$, is the measure of  $\Omega \times K$ defined as
				\[
				(\mu \otimes \nu_x )(U) \coloneqq \int_\Omega \int_K \chi_U(x,z) \;\mathrm{d} \nu_x(z) \;\mathrm{d} \mu(x), \qquad \text{for all $U \in \B(\Omega \times K)$}.
				\]
			}
		\end{enumerate}
		We denote the set of two-scale* Young measures by $\Y^2(\Omega;E)$.
	\end{definition}

	\subsubsection{Representation via two-scale* Young measures}
	Here and in what follows \enquote{$\epsilon \searrow 0$} will denote a sequence of positive reals converging to zero, which heuristically shall represent a microscopic length-scale.

	\begin{theorem}[representation]\label{thm:representation} Let $(\mu_\epsilon)_\epsilon\subset \mathcal M(\Omega;E)$ be a sequence of vector-valued measures with uniformly bounded total variation, i.e.,
		\[
		\sup_\epsilon |\mu_\epsilon|(\Omega) < \infty\,.
		\]
		The following representation result holds up to taking a subsequence of $(\mu_\epsilon)_\epsilon$.

		There exists a two-scale* Young measure $\bm{\nu} = (\nu,\lambda,\rho,\nu^\infty) \in \Y^2(\Omega;E)$ satisfying
		the following fundamental property. Let $f : \Omega \times \mathbb{T}^d \times E \to \mathbb{R}$ be a continuous integrand for which the recession function
		\[
		f^\infty(x,\xi,z) \coloneqq \lim_{\substack{(x',\xi',z') \to (x,\xi,z)\\t \to \infty}} \frac{f(x',\xi',tz')}{t} \qquad (x,\xi,z) \in \cl \Omega \times \mathbb{T}^d \times E,
		\]
		exists. Then, there exists a Radon measure $\mu_f$ on $\cl \Omega$ such that
		\[
		f(\,\,\sbullet\,\,, \,{\,\sbullet\,/\epsilon}\,,\, u_\epsilon\,) \, \mathscr L^d \restrict \Omega \; \overset{*}\rightharpoonup \; \mu_f \quad \text{as measures on $\cl \Omega$}\,,
		\]
		and $\mu_f$ is characterized by the values
		\begin{align*}
		\mu_f(\omega) & = \int_{\omega} \Big( \int_{\mathbb{T}^d} \Big( \int_E f(x,\xi , z) \;\mathrm{d} \nu_x(z) \Big) \;\mathrm{d} \xi \Big)  \;\mathrm{d} x \\
		& 	\qquad + \; \int_{\cl\omega} \Big(\int_{\mathbb{T}^d} \Big(\int_{\Bbb_E} f^\infty(x,\xi,z) \;\mathrm{d} \nu_x^\infty(z) \Big) \;\mathrm{d} \rho_x(\xi) \Big) \;\mathrm{d} \lambda(x)\,,
		\end{align*}
		where $\omega$ ranges among all Borel subsets of $\cl \Omega$. In this case we say that $(\mu_\epsilon)_\epsilon$ generates the two-scale*  Young measure $(\nu,\lambda,\rho,\nu^\infty)$, in symbols
		\[
		\mu_\epsilon \overset{\Ybf^2}{\to} \bm{\nu}\,.
		\]
	\end{theorem}

	The \emph{second} objective, is to endow two-scale* Young measures with a measure-theoretic toolbox tailored for applications in the calculus of variations. Our hope is to
	lay a transparent framework which casts the geometrical meaning of the blow-up methods (introduced by Fonseca \& M\"uller~\cite{fonseca1992quasi-convex-in}) into the context of two-scale analysis~\cite{bensoussan2011asymptotic-anal,nguetseng1989a-general-conve,allaire1992homogenization-}.
	Based on a localization principle, this comprehends the representation of integral functionals arising from {$\Gamma$}-convergence (see~\cite{de-giorgi1975sulla-convergen,de-giorgi1984g-operators-and}) in the context of homogenization of PDE-constrained structures~\cite{braides1996homogenization-,muller1987homogenization-,braides2000mathcal-a-quasi,matias2015homogenization-}.

	Formally, this toolbox consists of establishing the following properties:

	\subsubsection{Fundamental properties of two-scale* Young measures}
	\begin{enumerate}
		\item Compactness. In a natural way, two-scale* Young measures are elements of due dual of the class of integrands $\E(\Omega;E)$ (see~Section~\ref{integrands}). In Proposition~\ref{cor:Compactnessofyoungmeasures}
		we show the sequential \emph{ weak-$*$ compactness} (with respect to the weak-$*$ topology of $\E(\Omega;E)^*$) of {uniformly bounded} subsets $\mathcal{Y}$ of $\Y^2(\Omega;\E)$, that is, for sets such that
		\[
		\sup \; \setBB{ \int_{\cl \Omega} \int_{Z}\dpr{\nu_{\xi,x}, |\,\sbullet\,|} \;\mathrm{d} \xi \;\mathrm{d} x \; + \; \lambda(\cl \Omega) }{(\nu,\lambda,\rho,\nu^\infty) \in \mathcal{Y}} < \infty\,.
		\]
		This result is a fundamental step towards the proof of Theorem~\ref{thm:representation}.
		\item Localization. The relevance of Propositions~\ref{lem:locatregular} and~\ref{prop:loc_singular} is briefly explained as follows.  If a sequence $(\mu_\epsilon)_\epsilon$ generates a Young measure $(\nu,\lambda,\rho,\nu^\infty)$,
		then at $(\mathscr L^d + \lambda)$-almost every $x_0 \in \Omega$ we may find a \emph{blow-up sequence} of the original sequence (at $x_0$) that generates a (global) \emph{tangent two-scale* Young measure}
		\[
		\bm{D\nu(x_0)} = (\nu_{x_0},D\lambda,\Gamma_\#^{\xi_0}\rho_{x_0},\nu^\infty_{x_0}) \quad \text{for some $\xi_0 \in Z$, and $D\lambda \in \Tan_1(\lambda,x_0)$}\,.
		\]
		Hence, extending the concept of {tangent} Young measure  introduced in~\cite{rindler2011lower-semiconti} by Rindler; see also \cite{rindler2012lower-semiconti}.
	\end{enumerate}

	\subsubsection{Barycenter measures and second-scale convergence}

	Given a sequence $(u_j)_{k \in \mathbb{N}}$ that generates a two-scale* Young measure $\bm \nu$, the representation Theorem~\ref{thm:representation} yields  that
	\[
	u_j \overset{*}\rightharpoonup [\bm{\nu}]  \quad \text{as measures on $\cl \Omega$}\,,
	\]
	where $[\bm \nu]$ is the barycenter of $\bm{\nu}$ (see Definition~\ref{def:bary}). In this sense our notion of barycenter of a two-scale* Young measure coincides with
	the one for generalized Young measures. We then define a weak-$*$ $(\mathscr L^d \restrict \cl \Omega + \lambda^s)$-measurable map $\bbpr{\bm{\nu}} : \cl \Omega \to \mathcal M(Z;E)$, called the \emph{second-scale barycenter} of $\bm{\nu}$, which has the property that
	\[
	\text{$\mu_\epsilon$ two-scale converges to $(\mathscr L^d + \lambda^s) \otimes \bbpr{\bm{\nu}}_x$}\,,
	\]
	where the convergence above shall be understood as an extension of Nguetseng's concept of two-scale convergence~\cite{nguetseng1989a-general-conve,allaire1992homogenization-} (see Definitions~\ref{def:tsc} and~\ref{def:bary2}).
	In fact, via the compactness of two-scale* Young measures, we give a fairly short proof of the compactness of uniformly bounded sequences with respect to two-scale convergence;
	see Corollary~\ref{cor:tsc}.

	%

	\subsubsection{Structure of PDE-constrained two-scale* Young measures}

	Consider a homogeneous linear partial differential operator of order $k$ on $\mathbb{R}^d$ (with constant coefficients)  of the form
	\[
	\mathcal A = \sum_{\substack{|\alpha| = k\\\alpha \in \mathbb{N}^d}}  A_\alpha\,\partial^\alpha, \qquad A_\alpha \in \textrm{Lin}(E;F).
	\]
	A vector-valued measure $\mu \in \mathcal M(\Omega;E)$ is called $\mathcal A$-free provided that
	\[
	\mathcal A\mu = 0 \quad \text{in the sense of distributions on $\Omega$}.
	\]
	We say that two-scale* Young measure $\bm{\nu} \in \Y^2(\Omega;E)$ is $\mathcal A$-free if it is generated by a sequence of (asymptotically) $\mathcal A$-free measures. The set of such Young measures is denoted by $\Y^2_\mathcal A(\Omega;E)$. We establish the following rigidity properties of $\mathcal A$-free two-scale* Young measures:

	\begin{enumerate}
		\item The second-scale inherits the $\mathcal A$-free constraint. In Proposition~\ref{prop:secondbarycenterisafree} we show that the PDE-constraint is inherited also to the second-scale barycenter of any $\mathcal A$-free two-scale* Young measure (the
		constraint holds trivially for the barycenter), that is,
		\[
		\mathcal A\bbpr{\nu}_{x} = 0 \quad \text{in the sense of distributions on $Z$},
		\]
		for $(\mathscr L^d + \lambda^s)$-almost evert $x \in \cl \Omega$. The corresponding version of this result in terms of two-scale convergence is contained in Corollary~\ref{cor:Afree}.
		\item Based on the recent {developments~\cite{de-philippis2016on-the-structur} concerning} the structure of PDE-constrained measures, it further holds (see Corollary~\ref{cor:GuidoFilip})
		at $\lambda^s$-almost every $x \in \Omega$ that
		\[
		\frac{\;\mathrm{d} \bbpr{\bm{\nu}}_{x}}{\;\mathrm{d} |\bbpr{\bm{\nu}}_{x}^s|}(\xi) \in \Lambda_{\mathcal A} \quad \text{for $|\bbpr{\bm{\nu}}_{x}^s|$-almost every $\xi \in Z$},
		\]
		Here $\Lambda_\mathcal A$ is the so-called \emph{wave-cone} associated to $\mathcal A$ defined as
		\[
		\Lambda_{\mathcal A} \coloneqq \bigcup_{|\eta| =1} \ker \mathbb{A}^k(\eta) \subset E,
		\]
		and which consists of all Fourier amplitudes (vectors $z \in E$) where $\mathcal A$ is not elliptic with respect to one-directional oscillations.\footnote{Following standard notation, the $k$-homogeneous map
			\[
			\eta \mapsto \mathbb{A}^k(\eta) \coloneqq \sum_{|\alpha|=k} \eta^\alpha A_\alpha \in \textrm{Lin}(E,F), \qquad \eta \in \mathbb{R}^d,
			\]
			is the \emph{principal symbol associated to $\mathcal A$}. Using the Fourier transform it is immediate to verify that a vector $v \in E$ belongs to $\ker \mathbb{A}^k(\eta)$ --- for some $\eta \in \mathbb{R}^d \setminus \{0\}$ --- if and only if the one-directional function
			\[
			x \mapsto v \exp(2\pi\textrm{i}x \cdot \eta) \quad \text{is $\mathcal A$-free on $\mathbb{R}^d$}.
			\]}
		\item The support of the non-biting part of $\bm \nu$ (in the sense of Chacon; see Lemma~\ref{chacon}) is constrained by the differential constraint of $\mathcal A$. More specifically, in Lemma~\ref{lem:support} we show that
		\[
		\supp \, (\nu^\infty_{x, \xi}) \subset \spn \{\Lambda_{\mathcal A}\} \cap \partial \Bbb_E \quad \text{for $(\lambda^s \otimes \rho_x)$-almost every $(x,\xi) \in {\Omega} \times Z$}\,.
		\]
	\end{enumerate}

	\subsubsection{Applications to homogenization}

	We conclude our exposition solving a particular case of homogenization for PDE-constrained measures. We start by considering a family $\{I^\epsilon\}_{\epsilon >0}$ of functionals of the form
	\begin{equation}\label{eq:hom1}
	I^\epsilon(\mu) \coloneqq \int_\Omega f (x,x/\epsilon,\mu^{ac}(x)) \;\mathrm{d} x
	+ \int_\Omega f^\infty(x,x/\epsilon,\mu_s(x)) \;\mathrm{d} |\mu^s|(x),
	\end{equation}
	defined for measures $\mu = \mu^{ac} \mathscr L^d \restrict \Omega + \mu_s |\mu^s|\in \mathcal M(\Omega;\mathbb{R}^N)$.

	The integrand $f : \Omega \times \mathbb{T}^d \times \mathbb{R}^N \to \mathbb{R}$ is assumed to be a continuous integrand with linear-growth at infinity. Further we will require that $f(x,\xi,\,\sbullet\,)$ is convex for every $x,\xi \in \Omega \times \mathbb{T}^d$. The candidate measures $\mu \in \mathcal M(\Omega;\mathbb{R}^N)$ are assumed to satisfy the PDE-constraint
	\begin{equation}\label{eq:A}
	\mathcal A\mu  = 0 \quad \text{in the sense of distributions on $\Omega$}.
	\end{equation}

	Our goal is to show that as $\epsilon \searrow 0$, the rapidly oscillating variable ${x}/{\epsilon}$ averages out and the functionals $I^\epsilon$ \emph{converge} (in the context of $\Gamma$-convergence, detailed below) to an \enquote{homogenized} integral
	\begin{equation}\label{eq:hom2}
	I^\text{hom}(\mu) = \int_\Omega f_{*\mathcal{A}}(x,\mu^{ac}(x)) \;\mathrm{d} x
	+ \int_\Omega (f_{*\mathcal{A}})^\infty
	(x,\mu_s(x)) \;\mathrm{d} |\mu^s|(x),
	\end{equation}
	where the integrand $f_{*\mathcal{A}}$ is characterized by means of the cell minimization problem

	%

	%
	\begin{equation}\label{eq:defhom2}
	\begin{split}
	f_{*\mathcal A}(x,z) \coloneqq \inf\setBB{\int_{Q}& f(x,y,z + w(y)) \;\mathrm{d} y}{ \\
		& w \in \mathrm{C}^\infty(Q;E) \cap \ker \mathcal A, \int_Q w(y) \;\mathrm{d} y = 0}.
	\end{split}
	\end{equation}

	%
	%

	Let us briefly recall the notions of $\Gamma$-convergence and homogeneous envelope which will be required to give sense to our problem.
	Let $\{\epsilon_j\}_j$ be a sequence of positive numbers converging to zero. The \emph{$\Gamma$-limit inferior} of the \emph{sequence} of functionals $\{I^{\epsilon_j}\}_j$ with respect to the weak-$*$ convergence of measures is defined as
	\[
	\Gamma-\liminf  I^{\epsilon_j}{(\mu)}\coloneqq \inf\setB{\liminf_{j \to \infty}I^{\epsilon_j}(\mu_j)}{\text{$\mu_j \overset{*}\rightharpoonup \mu$ in $\mathcal M(\Omega;\mathbb{R}^N)$}}.
	\]
	We say that {a functional }$I$ is the \emph{$\Gamma$-limit inferior} of the \emph{family} of functionals $\{I^\epsilon\}_{\epsilon > 0}$ if
	\[
	I = \Gamma-\liminf_{j \to \infty} I^{\epsilon_j} \quad \text{for every sequence $\epsilon_j \searrow 0$}.
	\]
	In this case, we write
	\[
	I = \Gamma-\liminf_{\epsilon \searrow 0} I^\epsilon.
	\]

	%

	%
	%

	The main homogenization result is contained in the following Theorem.

	\begin{theorem}[$\Gamma$-$\liminf$ of $I^\epsilon$]\label{thm:hom} Let $f : \Omega \times \mathbb{T}^d \times E \to \mathbb{R}$ be a continuous integrand with linear-growth at infinity. Further assume that $f(x,\xi,\,\sbullet\,)$ is convex for all $x \in \Omega$ and all $\xi \in \mathbb{T}^d$. Then, the $\Gamma$-$\liminf$ of the family of functionals
		\begin{align*}
		\mu \mapsto I^\epsilon(\mu) =
		\int_\Omega & f(x,x/{\epsilon},\mu^{ac}(x)) \;\mathrm{d} x \\
		&  + \int_\Omega f^\infty(x,{x}/{\epsilon},\mu_s(x)) \;\mathrm{d} |\mu^s|(x)\,,
		\qquad \mu \in \mathcal M(\Omega;E) \cap \ker \mathcal A,
		\end{align*}
		with respect to the weak-$*$ convergence in $\mathcal M(\Omega;E)$, is given by the homogenized functional
		\[
		I^\textnormal{hom}(\mu) = \int_\Omega f_{*\mathcal{A}}(x,\mu^{ac}(x)) \;\mathrm{d} x + \int_\Omega (f_{*\mathcal{A}})^\infty(x,\mu_s(x)) \;\mathrm{d} |\mu^s|(x),
		\]
		defined for measures $\mu \in \mathcal M(\Omega;E) \cap \ker \mathcal A$.
	\end{theorem}

	%

	%
	%

	%

	%
	%
	%

	%
	%

	%

	%

	%
	%
	%
	%
	%
	%
	%

	%

	%

	%

	%

	%

	%

	%

	\section{Preliminaries and Notation}\label{sec:preliminaries}

	Here and in what follows $\Omega \subset \mathbb{R}^d$ is an open and bounded set with $\mathscr L^d(\partial \Omega) = 0$ (where $\mathscr L^d$ denotes the $d$-dimensional Lebesgue measure).
	To avoid cumbersome definitions we shall simply write $Z$ to denote the $d$-dimensional torus $\mathbb{T}^d$, and $Q$ to denote the closed $d$-dimensional unit cube $[0,1]^d$. To distinguish the $d$-dimensional Lebesgue measure between two locally $d$-dimensional euclidean spaces $E$ and $F$ we will often write
	$\mathscr L^d_E$ and $\mathscr L^d_F$ respectively. We denote the indicator function of a set $A$ by $\chi_A$. If $E,F$ are two Banach spaces, we denote by $\textrm{Lin}(E;F)$ the space of linear maps from $E$ to $F$.

	\subsection{Geometric measure theory}
	Let $X$ be a locally convex space. We denote by $\mathrm{C}_c(X)$ the space of compactly supported and continuous functions on $X$, and by $\mathrm{C}_0(X)$ we denote its completion with respect to the $\|\,\sbullet\,\|_\infty$ norm.
	The space $\mathrm{C}_c(X)$ is not a complete normed space in the usual sense, however, it is a complete metric space as the inductive union of Banach spaces $\mathrm{C}_0(K_m)$ where $K_m \subset X$ are compact and $K_m \nearrow X$.
	By the Riesz representation theorem, the space $\mathcal M(X)$ of bounded signed Radon measures on $X$ is the dual of $\mathrm{C}_0(X)$; a local argument of the same theorem states that the space $\mathcal M_\mathrm{loc}(X)$ of signed Radon measures on $X$ is the dual of $\mathrm{C}_c(X)$.

	We notate by $\mathcal M_\mathrm{loc}^+(X)$ the subset of non-negative measures.
	Since $\mathrm{C}_0(X)$ is a Banach space, the Banach--Alaoglu theorem and its characterizations hold and in particular bounded sets of $\mathcal M(X)$ are weak-$\ast$ metrizable.
	On the other hand, the local compactness of $\mathrm{C}_c(X)$ permits the existence of a complete and separable metric on $\mathcal M_\mathrm{loc}(X)$ with the property that convergence with respect to that metric is equivalent to the weak-$\ast$ convergence in $\mathcal M_\mathrm{loc}(X)$  (see Remark 14.15 in \cite{mattila1995geometry-of-set}). In a similar manner, for a finite dimensional euclidean space $E$, $\mathcal M(X;E)$ and $\mathcal M_\mathrm{loc}(X;E)$ will denote the spaces of {$E$-valued} bounded Radon measures and {$E$-valued} Radon measures respectively.

	The space $\mathcal M(X)$ is a normed space endowed with the \emph{total variation} norm
	\[
	|\mu|(X;E) \coloneqq \sup\setBB{\int_X \varphi \;\mathrm{d} \mu}{\phi \in \mathrm{C}_0(X;E), \|\phi\|_\infty \le 1}.
	\]
	The set of all positive Radon measures on $X$ with total variation equal to one is denoted by
	\[
	\mathrm{Prob}(X) \coloneqq \setB{\nu \in \mathcal M^+(X)}{\nu(X) = 1};
	\]
	the set of \emph{probability measures} on $X$

	The \emph{push-forward} of a measure $\mu \in \mathcal M(\Omega;E)$, with respect to a Borel map $T : \Omega \to \Omega'$, is formally defined through the \emph{change of variables} formula \enquote{$T_\# \mu = \mu \circ T^{-1}$} as follows. For every $\phi \in \mathrm{C}(\Omega';E)$, we define the measure $T_\#\mu \in \mathcal M(\Omega;E)$ via
	\[
	\int_{\Omega'} \phi \;\mathrm{d}(T_\#\mu) \coloneqq \int_\Omega \phi \circ T \;\mathrm{d} \mu.
	\]

	We define the linear action of a measure $\mu \in \mathcal M(X)$ on a function {$\phi \in \mathrm{C}_0(X)$}

	by the paring $\langle \phi,\mu \rangle =  \int_X \phi \;\mathrm{d} \mu$.
	If $X = Z$ is the $d$-dimensional torus (which is a compact manifold), then any map {$g\in \mathrm{C}_c(Z) =  \mathrm{C}_0(Z)$}
	can be represented by a periodic and continuous function on the $d$-dimensional semi-closed unit cube; in this case we write
	\[
	\dprb{ g,\mu} = \int_{Z}g \;\mathrm{d} \mu \coloneqq \int_{[0,1[^d} g \;\mathrm{d} \mu \quad \text{for $\mu \in \mathcal M(Z)$}.
	\]

	For a positive measure $\lambda \in \mathcal M(X)$ we write $\mathrm{L}_\lambda(X;F)$ to denote the set (space provided that $F$ is a space) of $\lambda$-measurable functions with values on $F \subset E$. For $p \in [1,\infty]$ we write
	\[
	\mathrm{L}^p_\lambda(X;F) \coloneqq \setBB{f \in \mathrm{L}(X;F)}{\int_X |f|_E^p \;\mathrm{d} \lambda < \infty},
	\]
	the space of $\lambda$-measurable functions with values on $F \subset E$ that are $p$-integrable. We will also use the short-hand notations $\mathrm{L}^p(X) \coloneqq \mathrm{L}^p(X;\mathbb{R})$ and
	\[
	\mathrm{L}_{\lambda,\mathrm{loc}}^p(X;F) = \setBB{f \in \mathrm{L}(X;F)}{\text{$\int_K |f|_E^p \;\mathrm{d} \lambda < \infty$ for all $K \Subset X$}}.
	\]
	If $F$ is an euclidean space, Riesz' representation theorem tells us that every vector-valued measure $\mu \in \mathcal M_\mathrm{loc}(X;F)$ can be written as
	\[
	\mu \; = \; f|\mu| \quad \text{for some $f \in \mathrm{L}^\infty_{|\mu|,\mathrm{loc}}(\Omega;\mathbb{S}_F)$};
	\]
	this decomposition is often referred as the \emph{polar decomposition} of $\mu$. The set of points $x \in X$ where
	\[
	 \lim_{r \downarrow 0} \int_{B_r(x_0)} |f(x) - f(x_0)| \;\mathrm{d} |\mu|(x) = 0,
	\]
	is called the set of $|\mu|$-Lebesgue points; this set has full $|\mu|$-measure. Another resourceful representation of a measure is given by its \emph{Radon--Nykod\'ym--Lebesgue decomposition}
	\[
	\mu \; = \; \frac{\;\mathrm{d} \mu}{\;\mathrm{d} \mathscr L^d} \, \mathscr L^d \; + \; \frac{\;\mathrm{d} \mu}{\;\mathrm{d} |\mu^s|} |\mu^s| \; = \; \mu^{ac} \mathscr L^d \; + \; \mu_s |\mu^s|
	\]
	where as usual $\mu^{ac} \in \mathrm{L}^1_\mathrm{loc}(X;F)$, $|\mu^s| \perp \mathscr L^d$, and $\mu_s \in \mathrm{L}_{|\mu^s|,\mathrm{loc}}(X;\mathbb{S}_F)$.

	Let $E,F$ be open or closed subsets of an euclidean finite dimensional space and let $\mu$ be a non-negative Radon measure on $E$.
	A map $\nu : E \mapsto \mathcal M(F) : x \mapsto \nu_x$ is said to be \emph{weak-$*$ $\mu$-measurable} if the map $x \mapsto \nu_x(B)$ is $\mu$-measurable for all Borel sets $B \in \mathfrak{B}(F)$.
	A simple method to check the $\mu$-measurability of such a measure valued map $x \mapsto \nu_x$ is to test the $\mu$-measurability of the map $x \mapsto \int_F g(x,z) \;\mathrm{d} \nu_x(z)$ for every $\mathfrak{B}(E) \times \mathfrak{B}(F)$-measurable function $g : E\times F \to \mathbb{R}$.
	The set of all weak-$*$ $\mu$-measurable maps $x \mapsto \nu_x$ endowed with the norm $|\nu|_\infty(E) \coloneqq \esssup_{(E,\mu)} |\nu_x|(F)$ conforms a Banach space which will shall denote by $\mathrm{L}^\infty_{\mu,\star}(E;\mathcal M(F))$. In the particular case that $\mu = \mathscr L^d$ we shall simply write $\mathrm{L}^\infty_{\star}(E;\mathcal M(F))$.

	Given a weak-$*$ $\mu$-measurable map $\nu \in \mathrm{L}^\infty_{\mu,\star}(E;\mathcal M(F))$ we can define the {generalized product} of $\mu$ and $\nu_x$ which is the measure taking the values
	\[
	(\mu \otimes \nu_x )(U) \coloneqq \int_E \int_F \chi_U(x,z) \;\mathrm{d} \nu_x(z) \;\mathrm{d} \mu(x), \qquad \text{for all $U \in \mathfrak{B}(E \times F)$}.
	\]
	The main step in the construction of Young measures (and multi-scale Young measures) relies on a well-known disintegration result for measures $\lambda \in \mathcal M^+(E \times F)$.
	Merely, it establishes a condition under which $\lambda = \mu \otimes \nu_x$, where $\mu$ is the push-forward of $\lambda$ under the projection onto $E$ and $\nu_x$ is a weak-$*$ $\mu$-measurable map of probability measures.
	This is recorded in the next theorem (for a proof see Theorem 2.28 in \cite{ambrosio2000functions-of-bo}).

	\begin{theorem}[disintegration]\label{thm:desintegration}
		Let $\lambda \in \mathcal M^+(E \times F)$ and let $\pi : E \times F \to E$ be the projection on the first factor.
		Assume that the push-forward measure $\mu \coloneqq \pi_\#\lambda \in \mathcal M(E)$ is a finite Radon measure.
		Then there exists a weak-$*$ $\mu$-measurable map $\nu \in \mathrm{L}^\infty_{\mu,\star}(E;\mathcal M(F))$, uniquely defined up to equivalence classes, such that $\lambda = \mu \otimes \nu_x$.
		Moreover,
		\[
		\nu_x\in \mathrm{Prob}(F) \quad \text{for $\mu$-almost every~$x$ in $E$}.
		\]
	\end{theorem}

	We close this section by introducing the notion of \emph{probability tangent measure} as introduced in \cite[Sec. 2.7]{ambrosio2000functions-of-bo}. Let $\{r_j\}_{j \in \mathbb{N}}$ be an infinitesimal sequence of positive real numbers ($r_j \searrow 0$). A \emph{local blow-up sequence} of a measure $\mu \in \mathcal M(\Omega;\mathbb{R}^N)$ at a point $x_0 \in \Omega$ is a sequence of (normalized) measures of the form
	\[
	\tau_j = \frac{1}{|\mu|(\cl{Q_{r_j}(x_0)})^{-1}}\,{T^{x_0,r_j}_\# \mu} \in \mathcal M(\cl Q),
	\]
	which are well-defined provided that $|\mu|(Q_{r_j}(x_0)) > 0$. A weak-$*$ limit $\tau \in \mathcal M(\cl{Q})$ of a local blow-up sequence is said to be a {local tangent measure}. We write
	\[
	\tau \in \Tan_1(\mu,x_0) \subset \mathrm{Prob}(\cl{Q}),
	\]
	to denote the set of all probability tangent measures (\emph{tangent measures} for short). \footnote{The term \enquote{tangent measure} refers to a more general object than the one referred to here. For the purposes of this paper, in particular the development of the localization principles in Section~\ref{sec:5}, it is technically more convenient to work with local tangent measures. However, our conclusions are compatible with the more general notion of \emph{tangent measure} introduced by Preiss \cite{preiss1987geometry-of-mea}.} At $|\mu|$-almost every $x_0 \in \Omega$, all tangent measures of $\mu$ at $x_0$ are constant multiples of a positive measure. More precisely, if $\mu =  f|\mu|$ is the polar decomposition of $\mu$, then at every $|\mu|$-Lebesgue point $x_0 \in \Omega$ it holds
	\[
	\Tan_1(\mu,x_0) = f(x_0) \Tan_1(|\mu|,x_0).
	\]

	In particular, every tangent measure $\tau \in \Tan_1(\mu,x_0)$ can be written as
	\[
	\tau = \frac{\;\mathrm{d} \mu}{\;\mathrm{d} |\mu|}(x_0)|\tau|.
	\]

	%

	%
	%
	%

	%
	%

	%

	%

	\subsection{Integrands}\label{integrands}

	Let $n$ be a non-negative integer. Throughout this section and the rest of the paper we will consider continuous integrands $f : \Omega \times Z^{n-1} \times E \to \mathbb{R}$ (recall that $\Omega \subset \mathbb{R}^d$ is an open and bounded set, $Z = \mathbb{T}^d$ is the $d$-dimensional torus, and $E$ is a euclidean space of finite dimension endowed with a norm $|\,\sbullet\,|=|\,\sbullet\,|_E$). Elements of $Z^{n-1}$ shall be denoted by $\bm{\xi} = (\xi_1,\dots,\xi_{n-1})$.

	Consider the transformation
	\[
	(Tf)(x,\bm{\xi},\hat z) \coloneqq (1 - |\hat z|)\cdot f\left(x,\bm{\xi},\frac{\hat z}{1 - |\hat z|}\right), \quad \text{$(x,\bm{\xi},\hat z) \in \Omega \times Z^{n-1} \times \Bbb_E$},
	\]
	defined for $f \in \mathrm{C}(\Omega \times Z \times \mathbb{R}^N)$ where $\Bbb_E$ is the open unit ball in $E$. The vector space
	\begin{align*}
	\E^n(\Omega;E) \coloneqq \setB{&f \in \mathrm{C}(\Omega \times Z^{n-1} \times E)}{\text{$(Tf)$ extends} \\& \qquad  \text{to a continuous function in $\mathrm{C}\big(\cl{\Omega \times Z^{n-1} \times \Bbb_E}\big)$}}
	\end{align*}
	endowed with the norm $\| f \|_{\E^n(\Omega;E)} \coloneqq \| Tf \|_{\infty}$, $\E^n(\Omega ;E)$ is a Banach space. Moreover, $T$ is a \emph{compactification} in the sense that
	\[
	T : \E^n(\Omega ;E) \to \mathrm{C}\big(\cl{\Omega \times Z^{n-1} \times \Bbb^E}\big)
	\]
	is an isometry of spaces with inverse
	\[
	(Sg)(x,\bm{\xi},z) \coloneqq (1 + |z|)\cdot g\left(x,\bm{\xi},\frac{z}{1 + |z|}\right), \quad \text{$(x,\bm{\xi}, z) \in \Omega \times Z^{n-1} \times E$}.
	\]
	The map $T : \E^n(\Omega;E) \to \mathrm{C}\big(\cl{\Omega \times Z^{n-1} \times \Bbb^E}\big)$ induces, through duality, an isomorphism $T^*$ (with inverse $S^*$) of the dual spaces; $T^*$ and $S^*$ are also isometries.

	Hence, a subset $X \subset (\E^n(\Omega;E))^*$ is sequentially weak-$\ast$ closed if and only if its image under $S^*$,
	$S^*X \subset \mathcal M\big(\cl{\Omega \times Z^{n-1} \times \Bbb^E}\big)$, is sequentially weak-$\ast$ closed in the sense of measures.

	Every $f \in \E^n(\Omega;E)$ has linear-growth at infinity, meaning there exists a positive constant $M$ (in this case given by $M = \|Tf\|_\infty$) for which $|f(x,\bm{\xi},z)| \le M(1 + |z|)$ for all $(x,\bm{\xi}, z) \in \Omega \times Z^{n-1} \times E$. This allows one to define a regularization at infinity: for $f\in \E^n(\Omega;E)$, we define the \emph{strong recession function} of $f$ as the limit
	\[
	f^\infty(x,\bm{\xi},z) \coloneqq \lim_{\substack{x' \to x\\ \bm{\xi}' \to \bm{\xi}\\t \to \infty}}\frac{f(x',\bm{\xi}',tz)}{t},
	\qquad \text{for $(x,\bm{\xi}, z ) \in \cl\Omega \times Z^{n-1}\times E$.}
	\]

	The continuity of $Tf$ ensures the limit is well-defined,
	and, in fact, $Tf(x,\bm{\xi},\hat z) = f^\infty(x,\bm{\xi},\hat z)$ for all $(x,\bm{\xi},\hat z) \in \cl \Omega \times Z^{n-1} \times \partial\Bbb_E$. Observe that $f^\infty$ is positively one-homogeneous in the $z$-variable and hence it can be recovered from the extended values of $Tf$.

	To complement our notation, we also define

	\[
	\E^0(E) \coloneqq  \setB{h \in \mathrm{C}(E)}{\text{$h^\infty$ exists}}
	\]
	{where $h^\infty$ is the strong recession function above without $x$ or $\xi$ dependance.}

	Observe that if $h \in \E^0(E)$, then the integrand $(\phi \otimes g \otimes h)(x,\bm{\xi},z) \coloneqq \phi(x)g(\bm{\xi})h(z)$

	belongs to the space $\E^n(\Omega;E)$ whenever $\phi \in \mathrm{C}(\cl \Omega)$ and $g \in \mathrm{C}(Z^{n-1})$.

	A function $g : Z^{n-1} \to \mathbb{R}$ is called \emph{upper (lower) semicontinuous} if its $Q$-periodic extension to $(\mathbb{R}^d)^{n-1}$ is  upper (lower) semicontinuous.
	Under this convention, the function $\xi \mapsto \xi$ with $\xi \in (0,1]$ is \textit{not} lower semicontinuous on $Z$ while the function $\xi \mapsto \xi$ for $\xi \in [0,1)$ \textit{is}.
	\begin{figure}[h]
		\captionsetup{justification=raggedright}
		\centering
		\begin{tikzpicture}
		\draw[->] (0,-0.2)--(0,1.2);
		\draw(-1,0)--(2,0);

		\draw(-1,0)--(0,1);
		\draw(0,0) -- (1,1);
		\draw(1,0) -- (2,1);
		\draw (0,0) circle (1.5pt);
		\draw (-1,0) circle (1.5pt);
		\draw (1,0) circle (1.5pt);
		\filldraw(0,1) circle (1.5pt);
		\filldraw(1,1) circle (1.5pt);
		\filldraw(2,1) circle (1.5pt);
		\end{tikzpicture}
		\hspace{2cm}
		\begin{tikzpicture}
		\draw[->] (0,-0.2)--(0,1.2);
		\draw(-1,0)--(2,0);

		\draw(-1,0)--(0,1);
		\draw(0,0) -- (1,1);
		\draw(1,0) -- (2,1);
		\draw (0,1) circle (1.5pt);
		\draw (1,1) circle (1.5pt);
		\draw (2,1) circle (1.5pt);
		\filldraw(-1,0) circle (1.5pt);
		\filldraw(0,0) circle (1.5pt);
		\filldraw(1,0) circle (1.5pt);
		\end{tikzpicture}

		\caption{Lower semicontinuity of periodic functions, the functions $\xi \mapsto \xi$ for $\xi \in (0,1]$ and $\xi \in [0,1)$. } \label{fig:Lscofperiodic}
	\end{figure}
	In the following we say that a {function} $f:\Omega \times Z^{n-1} \times E \to \mathbb{R}$ is upper (lower) semicontinuous if $f(\,\sbullet\,,\xi,\,\sbullet\,)$ is upper (lower) semicontinuous on $\Omega \times E$ for all $\bm{\xi} \in Z^{n-1}$ and $f(x,\,\sbullet\,,z)$ is upper (lower) semicontinuous on $Z^{n-1}$.

	In general, the strong recession function of a Borel (or even continuous) integrand  might fail to exist. Instead, one can always define the \emph{upper} and \emph{lower recession} functions by setting
	\[
	f^\#(x,\bm{\xi},z) \coloneqq \limsup_{\substack{x' \to x\\ \bm{\xi}' \to \bm{\xi}\\t \to \infty}}\frac{f(x',\bm{\xi}',tz)}{t},
	\qquad \text{for $(x,\bm{\xi}, z ) \in \cl\Omega \times Z^{n-1}\times E$}
	\]
	and
	\[
	f_\#(x,\bm{\xi},z) \coloneqq \liminf_{\substack{x' \to x\\ \bm{\xi}' \to \bm{\xi}\\t \to \infty}}\frac{f(x',\bm{\xi}',tz)}{t},
	\qquad \text{for $(x,\bm{\xi}, z ) \in \cl\Omega \times Z^{n-1}\times E$}.
	\]
	In the same way as $f^\infty$, $f^\#$ and $f_\#$ are positively $1$-homogeneous on the $z$-variable. Observe that $f^\#$ (respectively $f_\#$) is nothing else than the upper (lower) semicontinuous regularization of $Tf$ on $\cl{\Omega \times Z^{n-1} \times \Bbb_E}$. This observation is the key argument behind the following approximation result

	(its proof follows the same arguments given in the proof \cite[Lemma~2.3]{alibert1997non-uniform-int} in the context of classical generalized Young measures).

	\begin{proposition}\label{prop:aprox}

		Let $f : \Omega \times Z^{n-1} \times E \to \mathbb{R}$ be an upper semicontinuous integrand and assume there exists $M>0$ such that $f(x,\xi,z) \geq -M (1+|z|)$ for all $(x,\xi,z) \in \Omega\times Z\times E$.
		Then, there exists a non-increasing sequence of functions $(f_m)_{m \in \mathbb{N}} \subset \E^n(\Omega;E)$ such that
		\begin{align*}
		\inf_{m > 0} f_m = \lim_{m \to \infty} f_m & = f \qquad  \textrm{pointwise on $\;\Omega \times Z^{n-1} \times E$, and} \\
		\inf_{m > 0} f_m^\infty = \lim_{m \to \infty} f_m^\infty & = f^\#  \qquad \textrm{pointwise on $\;\cl \Omega \times Z^{n-1} \times E$}.
		\end{align*}
	\end{proposition}

	\begin{remark} The analogous statement holds for lower semicontinuous integrands $f$ and their lower recession function $f_\#$ {for a non-decreasing sequence of approximating functions.}
	\end{remark}

	\section{Generalized multi-scale Young measures}\label{sec:3}

	Our construction extends the exposition in \cite{kristensen2010characterizatio} which itself goes back to the seminal works of DiPerna \& Majda \cite{diperna1987oscillations-an} and Alibert \& Bouchitt\'e \cite{alibert1997non-uniform-int}.

	Heuristically, $n$ will represent the number of hierarchical scales of the vector
	\[
	\underbrace{(x,\xi_1,\dots,\xi_{n-1})}_{n}
	\]
	where $x \in \Omega$ is the \emph{macroscopic} variable and the $\xi_i$'s conform a hierarchical family of periodic \emph{microscopic} scales.
	The indexing corresponding to the microscopic scales reflects a disassociation between the $i$-th scale and the finer $(i + 1)$-th scale. Mathematically, this is reflected by the homogenization with respect to scales
	\[
	\bigg(x,\underbrace{\frac{x}{\epsilon_1}}_{\xi_1},\dots,\underbrace{\frac{x}{\epsilon_{n-1}}}_{\xi_{n-1}}\bigg)
	\]
	where we assume that $\epsilon_{i +1} \ll \epsilon_i$ for each $i = \{1,\dots,n-2\}$.
	We are now ready to introduce the notion of multi-scale Young measure.

	We are now ready to give the precise definition of multi-scale* Young measure and state their main properties.

	\begin{definition}[$n$-scale* Young measure]
		A four-tuple $ \bm{\nu} =(\nu, \lambda, \rho, \nu^\infty)$
		is called a \textit{generalized $n$-scale* Young measure} on $\cl \Omega$ with values on $E$ provided that
		\begin{enumerate}
			\item[\textnormal{(i)}] $\lambda$ is a positive measure on $\cl \Omega$,
			\item[\textnormal{(ii)}] $\rho$ is a weak-$*$ $\lambda$-measurable map from $\cl \Omega$ into the set $\mathrm{Prob}(Z^{n-1})$ of probability measures over the product of $(n-1)$ copies of the $d$-dimensional torus, and
			\item[\textnormal{(iii)}] the map $(x,{\xi}) \mapsto \dpr{|\,\sbullet\,|,\nu_{x,{\xi}}}$ belongs to  $\mathrm{L}^1(\Omega \times Z^{n-1})$.
		\end{enumerate}
		Additionally, $\nu$ and $\nu^\infty$ are weak- and weak-$*$ measurable maps on $\cl \Omega \times Z$ respectively:
		\begin{enumerate}
			\item[\textnormal{(iv)}] $\nu$ is a weak measurable map from $\cl \Omega \times Z^{n-1}$ into the set $\mathrm{Prob}(E)$ of probability measures over $E$, and
			\item[\textnormal{(v)}] $\nu^\infty$ is a weak-$*$ $(\lambda \otimes \rho_x)$-measurable map from $\cl \Omega \times Z^{n-1}$ into the set $\mathrm{Prob}(\partial \Bbb_E)$ of probability measures supported on the unit sphere in $E$.
		\end{enumerate}

		The set of all $n$-scale Young measures is denoted by $\Y^n(\Omega;E)$.

	\end{definition}

	\subsection{Construction}

	We shall restrict our analysis to two-scale* Young measures  $\Y^2(\Omega;E)$. The results and ideas behind the proofs extend analogously to $n$-scale* Young measures.

	Two-scale* Young measures conform a set of dual objects to the space of integrands $\E^2(\Omega;E)$ in the following way:
	For $f \in \E^2(\Omega;E)$ and $\bm{\nu} \in \Y^2(\Omega,E)$, we define a bilinear product by setting
	\begin{equation}\label{eq:Yproduct}
	\begin{split}
	\ddprb{f, \bm{\nu}} & \coloneqq \int_{\Omega} \bigg(\int_{Z} \dprb{f(x,\xi,\,\sbullet\,),\nu_{x,\xi}} \;\mathrm{d} \xi \bigg) \;\mathrm{d} x \\
	& \; \qquad
	+ \int_{\cl\Omega}\bigg(\int_{Z} \dprb{f^\infty(x,\xi,\,\sbullet\,),\nu_{x,\xi}^\infty}\;\mathrm{d} \rho^\nu_x(\xi)\bigg) \;\mathrm{d} \lambda^\nu (x).
	\end{split}
	\end{equation}
	By the definition of two-scale* Young measure and the  the linear-growth of the elements of $\E^2(\Omega;E)$ we can estimate its norm by\begin{align*}
	|\ddprb{f,\bm{\nu}}| & \le \|Tf\|_\infty\bigg(\int_{\Omega \times Z} \dprb{1 + |\,\sbullet\,|,\nu_{x,\xi}} \;\mathrm{d} \xi \;\mathrm{d} x + \int_{\cl \Omega} \;\mathrm{d} \lambda^\nu(x)\bigg) \\
	& =  \|Tf\|_\infty\big(\mathscr L^d(\Omega) + \|g^\nu\|_{\mathrm{L}^1(\Omega\times Z)} + \lambda^\nu(\cl \Omega)\big),

	\end{align*}
	where we used the short-hand notation $g^\nu(x,\xi) \coloneqq \dprb{|\,\sbullet\,|,\nu_{x,\xi}}$.

	It follows that $\ddprb{\,\sbullet\,,\bm{\nu}} \in \E^2(\Omega;E)^*$, and, in this sense, we shall identify $\Y^2(\Omega;E)$ with a subset of $\E^2(\Omega;E)^*$. For $\bm{\nu}_j,\bm{\nu} \in \Y^{{2}}(\Omega;E)$ we write $\bm{\nu}_j \overset{*}\rightharpoonup \bm{\nu}$ in $\E^2(\Omega;E)$ if
	\[
	\ddprb{f,\bm{\nu}_j} \to \ddprb{f,\bm{\nu}} \quad \text{for all $f \in \E^2(\Omega;E)$}.
	\]
	In this case we say that $\bm{\nu}_j$ weak-$\ast$ converges to $\bm{\nu}$ as two-scale* Young measures.

	The following weak-$*$ semicontinuity results hold.

	\begin{lemma}[semicontinuity] Let $f:\Omega \times Z \times E \to \mathbb{R}$ be a lower semicontinuous integrand with linear-growth at infinity. Then, the functional
		\begin{align}
		\bm{\nu} \mapsto \ddprb{f,\bm{\nu}}_\# \coloneqq \int_{\Omega} \bigg(& \int_Z \dprb{f(x,{\mathbf \xi},\,\sbullet\,),\nu_{x,{\mathbf \xi}}} \;\mathrm{d} \xi \bigg) \;\mathrm{d} x \\
		&+\int_{\cl \Omega} \bigg(\int_Z \dprb{f_\#(x,\xi,\,\sbullet\,),\nu^\infty_{x,\xi}} \;\mathrm{d} \rho_x(\xi)\bigg) \;\mathrm{d} \lambda(x)
		\end{align}
		is weak-$*$ lower semicontinuous in $\Y^2(\Omega;E)$.

		Similarly, if $f:\Omega \times Z \times E \to \mathbb{R}$ is a Borel integrand that is upper semicontinuous and has linear growth at infinity, then the functional
		\begin{align}
		\bm{\nu} \mapsto \ddprb{f,\bm{\nu}}^\# \coloneqq \int_{\Omega} \bigg(& \int_Z \dprb{f(x,{\mathbf \xi},\,\sbullet\,),\nu_{x,{\mathbf \xi}}} \;\mathrm{d} \xi \bigg) \;\mathrm{d} x  \\
		& + \int_{\cl \Omega} \bigg(\int_Z \dprb{f^\#(x,\xi,\,\sbullet\,),\nu^\infty_{x,\xi}} \;\mathrm{d} \rho_x(\xi)\bigg) \;\mathrm{d} \lambda(x)
		\end{align}
		is weak-$*$ upper semicontinuous in $\Y^2(\Omega;E)$.
	\end{lemma}
	\begin{proof} Since the two statements are equivalent modulo taking $(-f)$ in place of $f$, we shall only argue the case of upper-semicontinuity. Let $f : \Omega \times Z \times E \to \mathbb{R}$ be an upper-semicontinuous integrand and let $\bm{\nu}_j \overset{*}\rightharpoonup \bm{\nu}$ in $\Y^2(\Omega;E)$. Further let $(f_m)_{m\in\mathbb{N}} \subset \E^2(\Omega;E)$ be the monotone approximating sequence provided by Proposition~\ref{prop:aprox}. Fix $m \in \mathbb{N}$, then by continuity of the pairing and the properties of the approximating sequence we get
		\begin{align*}
		\limsup_{j \to \infty} \, \ddprb{f,\bm{\nu}_j}^\# & \le \lim_{j \to \infty} \,  \ddprb{f_m, \bm{\nu}_j}  \\
		& =  \ddprb{f_m,\bm{\nu}}.
		\end{align*}
		The conclusion then follows by letting $m \to \infty$ in the inequality above and arguing with the monotone convergence theorem.
	\end{proof}
	An immediate consequence of this result is that the weak-$*$ continuity of the map $\bm{\nu} \mapsto \ddprb{f,\bm{\nu}}$ can be extended from $\E^2(\Omega;E)$ to the  larger class
	\[
	\mathbf{R}^2(\Omega;E) \coloneqq \setB{f \in \mathrm{C}(\Omega \times Z \times E)}{f^\infty \equiv f^\# \equiv f_\#}\,
	\]
	of continuous integrands possessing a \emph{strong recession function}. In particular $z$-convex integrands with linear growth at infinity belong to this class:

	\begin{proposition} Let $f : \Omega \times Z \times E \to \mathbb{R}$ be a continuous integrand with linear growth at infinity. If $f(x,\xi,\,\sbullet\,)$ is convex for every $(x,\xi) \in \Omega\times Z$, then $f \in \mathbf{R}^2(\Omega;E)$.
	\end{proposition}
	\begin{proof}
		The proof follows from Remark~2.4 in \cite{alibert1997non-uniform-int}.
	\end{proof}

	If $\bm{\nu} \in \Y^2(\Omega;E)$, duality yields
	\[
	\dprb{g,S^*\bm{\nu}}
	= \ddprb{Sg,\bm{\nu}}\quad \text{for all $g \in \mathrm{C}\big(\cl{\Omega \times Z \times \Bbb_E}\big)$}.
	\]
	In other words, the following  diagram
	\begin{equation}\label{diagram}
	\begin{tikzcd}
	\Y^2(\Omega;E)   \arrow[r, hook] \arrow[d,"S^*"]
	& \E^2(\Omega;E)^*
	\arrow[d,"S^*"] \\
	S^*(\Y^2(\Omega;E)) \arrow[r,hook]
	& \mathcal M(\cl{\Omega\times Z \times \Bbb_E})

	\end{tikzcd}
	\end{equation}
	is commutative.
	Likewise, every $\mu = S^*\bm{\nu} \in S^*(\Y^2(\Omega;E))$ verifies the identity
	\[
	\dprb{Tf,\mu} = \ddprb{f,\bm{\nu}} \quad \text{for all $f \in \E^2(\Omega;E)$}.
	\]
	Using that $S^*$ is an isometry of Banach spaces, it can be deduced the isomorphism lowers to a weak-$*$ isomorphism
	$S^*:(\E^2(\Omega;E)^*,\text{w-$*$}) \to (\mathcal M(\cl{\Omega \times Z \times \Bbb_E}),\text{w-$*$})$.

	%
	%

	It is then straightforward  to check that
	\[
	S^*\bm{\nu}_j \overset{*}\rightharpoonup \mu \quad \text{in $\mathcal M(\cl{\Omega \times Z \times \Bbb_E}) \quad \Leftrightarrow \quad \bm{\nu}_j \overset{*}\rightharpoonup T^*\mu$ \quad in $\E^2(\Omega;E)$}.
	\]
	\begin{remark}[topological isomorphism]\label{rem:iso} Topological properties are preserved under isomorphisms, in particular

		\begin{gather*}
		S^*[\Y^2(\Omega;E)] \; \text{is sequentially weak-$\ast$ closed in $\mathcal M(\cl{\Omega \times Z \times \Bbb_E})$} \\
		\Leftrightarrow\\
		\Y^2(\Omega;E) \; \text{is sequentially weak-$\ast$ closed in $\E^2(\Omega;E)^*$}\,.\\

		\end{gather*}
	\end{remark}

	The following characterization will play a fundamental role in proving the weak-$*$ compactness of uniformly bounded sets of two-scale* Young measures. We follow the presentation given in Lemma 2 in \cite{kristensen2010characterizatio} for $\Y^1(\Omega;E)$.
	\begin{lemma}\label{lem:representationformula}
		The set $S^*[\Y^2(\Omega;E)]$ consists precisely of all positive measures $\mu \in \mathcal{M}^+(\cl{\Omega \times Z \times \Bbb_E})$ satisfying following property: for all $\phi \in \mathrm{C}(\Omega)$ and all $g \in \mathrm{C}(Z)$ it holds that

		\begin{align}\label{eq:representation}
		\int_{\cl{\Omega \times Z\times \Bbb_E}} \, (\phi \otimes g)(x,\xi)\,(1-|\hat z|) \;\mathrm{d}\mu (x,\xi,\hat z) = \int_\Omega \phi(x) \;\mathrm{d} x \int_{Z}g(\xi) \;\mathrm{d} \xi\,.
		\end{align}

	\end{lemma}
	\begin{proof}
		\proofstep{Necessity.} Fix $\mu = S^*\bm{\nu}$.
		Applying $\mu$ on the function $Tf = \phi \otimes g \otimes (1 -|\,\sbullet\,|)$ and using that $f = \phi \otimes g \otimes \chi_{E}$ (and hence $f^\infty \equiv 0$) we get
		\[
		\dprb{ Tf, \mu } = \dprb{ Tf, S^*\bm{\nu} } = \ddprb{ f, \bm{\nu} } =\int_\Omega \int_{Z} \phi(x) g(\xi) \;\mathrm{d} \xi \;\mathrm{d} x,
		\]
		which is (\ref{eq:representation}).
		The positivity of $\mu$ follows from  the positivity of $\nu$ and $\nu^\infty$ and the definition of $S$.

		\proofstep{Sufficiency.} For the reverse statement, fix a measure $\mu \in \mathcal{M}^+(\cl{\Omega \times Z \times \Bbb_E})$ satisfying~(\ref{eq:representation}). We want to find $\bm{\nu} \in \Y^2(\Omega;E)$ satisfying  $\mu = S^*\bm{\nu}$.

		By disintegration (see Theorem~\ref{thm:desintegration}), we may decompose $\mu$ as a double semi-product $\mu = \hat \lambda \otimes \hat \rho_x \otimes \hat \nu_{x,\xi}$
		where $\hat\lambda \in \mathcal M^+(\cl\Omega)$, $\hat\rho \in \mathrm{L}^\infty_{\hat\lambda,\star}(\cl\Omega;\mathrm{Prob}(Z))$ and $\hat \nu \in \mathrm{L}^\infty_{\hat\lambda \otimes \hat\rho_x,\star}(\cl\Omega\times Z;\mathrm{Prob}( \cl{\Bbb_E}))$.
		Let us further set
		\[
		u(x,\xi) \coloneqq \dpr{1 - |\,\sbullet\,|,\hat\nu_{x,\xi}} \in \mathrm{L}^{1}_{\hat\lambda \otimes \hat\rho_x }(\cl \Omega \times Z)\,.
		\]
		Thus, we may re-write~\eqref{eq:representation} as the equivalence $u \, (\hat \lambda \otimes \hat \rho_x) \equiv \mathscr L^d_\Omega \otimes \mathscr L^d_Z$.

		On the one hand, this gives $u \, [(\hat \lambda^s \otimes \hat\rho_x) + (\mathscr L^d_\Omega \otimes \hat \rho_x^s)] \equiv 0$, which, in turn, is equivalent to
		\begin{equation}\label{eq:tec1}
		\supp\,(\hat \nu_{x,\xi}) \subset \partial \Bbb_E \quad \text{$[(\hat \lambda^s \otimes \hat\rho_x) + (\mathscr L^d_\Omega \otimes \hat \rho_x^s)]$-almost everywhere}.
		\end{equation}
		On the other hand, the same equivalence yields
		\begin{equation}\label{eq:tec2}
		\lambda^{ac}(x) \cdot \rho_x^{ac}(\xi) \cdot  u(x,\xi) =  1 \quad \text{$\mathscr L^d_\Omega \otimes \mathscr L^d_Z$-almost everywhere}.
		\end{equation}
		In particular $u$ is $\mathscr L^d_\Omega \otimes \mathscr L^d_Z$-measurable.

		The goal now is to exhibit a $\bm{\nu} = (\nu,\lambda,\rho,\nu^\infty) \in \Y^2(\Omega;E)$ for which $S^*\bm{\nu}= \mu$.

		\textbf{Construction of $\bm{\nu}$.} Let us being by defining  the weak-$*$ $(\mathscr L^d_\Omega \otimes \mathscr L^d_Z)$-measurable map
		\[
		(x,\xi) \mapsto \nu_{x,\xi}\in \mathrm{Prob}(E)\,,
		\]
		where each $\nu_{x,\xi}$ is the probability measure satisfying
		\[
		\dpr{\nu_{x,\xi},h} \coloneqq \frac{1}{u(x,\xi)}  \int_{\Bbb_E}Th(\hat z) \;\mathrm{d}\hat \nu_{x,\xi}(\hat z) \qquad \text{for all $h \in \mathrm{C}(\beta E)$}\,.
		\]
		That $\nu$ is a weak-$*$ measurable map follows from the measurability of $u$ and the properties of $\hat \nu \restrict \partial \Bbb_E$ in terms of the measure $(\hat \lambda \otimes \hat \rho_x)$. To check that each $\nu_{x, \xi}$ is indeed a probability measure ($\mathscr L^d_\Omega \otimes \mathscr L^d_Z$-almost everywhere) follows by testing with $\chi_E$ and using the definition of $u$ (recall that $T\chi_E \equiv 1 - |\,\sbullet\, |$ as functions on $\Bbb_E$). Moreover, $\dpr{1 + |\,\sbullet\,|,\nu_{x,\xi}} = u(x,\xi)^{-1} \hat\nu_{x, \xi}(\Bbb_E)$ and hence by~\eqref{eq:tec2} we infer the map $(x,\xi) \mapsto \dpr{1 + |\,\sbullet\,|,\nu_{x,\xi}}$ is integrable on $\Omega \times Z$.

		We define the remaining  $\lambda,\rho$ and $\nu^\infty$ as follows. First, we set
		\[
		\lambda \coloneqq \bigg(\int_Z \hat\nu_{x,\xi}(\partial \Bbb_E) \;\mathrm{d} \rho_x(\xi)\bigg) \, \hat \lambda \in \mathcal M^+(\cl \Omega)\,.
		\]
		Once this positive measure has been defined, we define a map $\rho$ from $\cl \Omega$ into the set $\mathrm{Prob}(Z)$ of probability measures over the $d$-dimensional torus by setting
		\[
		x \mapsto \rho_x \coloneqq \bigg(\int_Z \hat\nu_{x,\xi}(\partial \Bbb_E)\;\mathrm{d}\hat\rho_x(\xi)\bigg)^{-1} m_x\,\hat\rho_x\,.
		\]
		where we have used the short-hand notation $m_x(\xi) \coloneqq \hat\nu_{x,\xi}(\partial \Bbb_E)$. Since by definition $m_x$ is a $\hat \rho_x$-measurable map, we infer that $\rho$ is
		weak-$*$ $\lambda$-measurable.

		Lastly, we define a map $\nu^\infty$ from $\cl \Omega \times Z$ into $\mathrm{Prob}(\partial \Bbb_E)$ by setting
		\[
		(x,\xi) \mapsto \nu_{x,\xi}^\infty \in \mathrm{Prob}(\partial \Bbb_E)\,,
		\]
		where each $\nu_{x,\xi}^\infty$ is given (in terms of duality) by
		\[
		\dpr{h,\nu_{x,\xi}} \coloneqq \,\thetant-_{\partial \Bbb_E}h(\hat z) \;\mathrm{d} \hat \nu_{x,\xi}(\hat z) \qquad \text{for all $h \in \mathrm{C}(\partial \Bbb_E)$}\,.
		\]
		That $\nu^\infty$ is a weak-$*$ $(\lambda \otimes \rho_x)$-measurable map is then a consequence of the weak-$*$ measurability of $\hat \nu$, and the way $\lambda$ and $\rho$ are defined.
		Altogether these properties imply  $\bm{\nu} \coloneqq (\nu,\lambda,\rho,\nu^\infty) \in \Y^2(\Omega;E)$.

		\textbf{Pre-image property $(S^* \bm{\nu} = \mu$).} Let $f \in \E^2(\Omega;E)$ be a fixed but arbitrary integrand, later we shall exploit this choice through duality. By construction of $\bm{\nu}$ we get that
		\begin{align*}
		\ddprb{f,\bm{\nu}} & = \int_{\Omega}\int_Z \dpr{f(x,\xi,\,\sbullet\,),\nu_{x,\xi}} \;\mathrm{d} \xi \;\mathrm{d} x + \int_{\cl \Omega}\int_Z\dpr{f^\infty(x,\xi,\,\sbullet\,),\nu^\infty_{x,\xi}} \;\mathrm{d}\rho_x(\xi) \;\mathrm{d} \lambda(x)\\
		& = \int_{\Omega}\bigg(\int_Z\bigg(\int_{\Bbb_E}Tf(x,\xi,z)\;\mathrm{d}\hat\nu_{x,\xi}(z)\bigg)\,\hat\rho_x^{ac}(\xi)\;\mathrm{d}\xi\bigg)\,\hat\lambda^{ac}(x) \;\mathrm{d} x \\
		& \quad + \int_{\Omega}\bigg(\int_Z\bigg(\,\thetant-_{\partial \Bbb_E} Tf(x,\xi,z)\;\mathrm{d}\hat\nu_{x,\xi}(z)\bigg)\,m_x(\xi)\,\hat\rho^{ac}_x(\xi) \;\mathrm{d}\xi\bigg)\,\hat\lambda^{ac}(x) \;\mathrm{d} x \\
		& \qquad + \int_{\Omega}\bigg(\int_Z\bigg(\int_{\partial \Bbb_E} Tf(x,\xi,z)\;\mathrm{d}\hat\nu_{x,\xi}(z)\bigg)\;\mathrm{d}\hat\rho_x^s(\xi)\bigg)\,\hat\lambda^{ac}(x) \;\mathrm{d} x \\
		& \quad\qquad + \int_{\cl\Omega}\bigg(\int_Z\bigg(\int_{\partial \Bbb_E} Tf(x,\xi,z)\;\mathrm{d}\hat\nu_{x,\xi}(z)\bigg)\,\;\mathrm{d}\hat\rho_x(\xi)\bigg)\,\hat\lambda^{s}(x)\,,
		\end{align*}
		where in passing to the last equality we used that
		$f^\infty(x,\xi,\,\sbullet\,) \equiv Tf(x,\xi,\,\sbullet\,)$ as functions over $\partial \Bbb_E$. Furthermore, since $m_x(\xi) = \hat\nu_{x,\xi}(\partial \Bbb_E)$, the first two lines of the last equality above add up to
		\[
		\int_{\cl{\Omega \times Z \times \Bbb_E}} Tf \;\mathrm{d}(\hat\lambda^{ac}\mathscr L^d_\Omega\otimes\hat\rho_x^{ac}\mathscr L^d_Z \otimes\hat\nu_{x,\xi})\,.
		\]
		On the other hand, using~\eqref{eq:tec1} we may re-write the last two lines in the expression of $\ddprb{f,\bm{\nu}}$ as
		\[
		\int_{\cl{\Omega \times Z \times \Bbb_E}} Tf  \;\mathrm{d}(\hat\lambda^{ac}\mathscr L^d_\Omega\otimes\hat\rho_x^s\otimes\hat\nu_{x,\xi})
		+ \int_{\cl{\Omega \times Z \times \Bbb_E}} Tf  \;\mathrm{d}(\hat\lambda^s\otimes\hat\rho_x\otimes\hat\nu_{x,\xi})\,.
		\]
		Regrouping these three summands together we deduce
		\[
		\ddprb{f,\bm{\nu}}= \int_{\cl{\Omega \times Z \times \Bbb_E}} Tf \;\mathrm{d}\mu = \dpr{Tf,\mu} \qquad \text{for all $f \in \E^2(\Omega;E)$}\,.
		\]
		Equivalently,  by a duality argument,
		\[
		\dprb{\Phi,\mu} = \ddprb{S\Phi,\bm{\nu}} = \dprb{\Phi,S^*\bm{\nu}} \quad \text{for all $\Phi \in \mathrm{C}(\cl{\Omega \times Z \times \Bbb_E})$}\,.
		\]
		Thence $\mu \equiv S^*\bm{\nu}$ as measures in $\cl{\Omega \times Z \times \Bbb_E}$. This proves the sufficiency.
	\end{proof}

	%
	%
	%
	%

	%
	%

	A direct consequence of this characterization and Remark~\ref{rem:iso} is the following fundamental property of Young measures. Here and in what follows, we shall write
	\[
	{[\,\sbullet\,]_{E}} \coloneqq \chi_{\cl \Omega} \otimes \chi_Z \otimes (1 + |\,\sbullet\,|_E).
	\]

	\begin{proposition}[compactness of two-scale Young measures]\label{cor:Compactnessofyoungmeasures} The set of two-scale* Young measures $\Y^2(\Omega;E)$ is sequentially weak-$\ast$ closed in $\E^2(\Omega;E)^*$.
		Moreover, each subset $\mathcal Y \subset \Y^2(\Omega;E)$ satisfying

		\[
		\sup \; \setBB{\int_{\cl \Omega}\int_Z \dpr{1 + |\,\sbullet\,|,\nu_{x, \xi}} \;\mathrm{d} \xi \;\mathrm{d} x + \lambda(\cl \Omega)}{\bm{\nu} = (\nu,\lambda,\rho,\nu^\infty) \in \mathcal{Y}} < \infty
		\]

		is pre-compact with respect to the relative weak-$\ast$ topology on $\Y^2(\Omega;E)\subset \E^2(\Omega;E)^*$.
	\end{proposition}
	\begin{proof}
		To verify that $\Y^2(\Omega;E)$ is sequentially weak-$*$ closed it suffices (by Remark~\ref{rem:iso} and Lemma~\ref{lem:representationformula}) to observe that~\eqref{eq:representation} is a closed property with respect to the sequential weak-$*$ convergence of measures in $\mathcal M^+(\cl \Omega \times Z \times \Bbb_E)$. Let $\mathcal{Y}$ as in the assumptions and observe that $T(1 + |\,\sbullet\,|) =  \chi_{\Bbb_E}$ so that
		\[
		\sup_{\mu \in S^*[\mathcal{Y}]} \, |\mu|(\cl \Omega\times Z \times \Bbb_E) = \sup_{\bm{\nu} \in \mathcal{Y}} \, \ddprb{{[\,\sbullet\,]_{E}},\bm{\nu}} < \infty.
		\]
		The Banach--Alaoglu theorem tells us that $S^*[\mathcal{Y}]$ is pre-compact with respect to the weak-$*$ topology of measures. Hence, again by Remark~\ref{rem:iso}, $\mathcal{Y}$ is
		weak-$*$ pre-compact with respect to the relative weak-$*$ topology of $\E^2(\Omega;E)^*$.
	\end{proof}

	We close this section with an important separability property. The proof of this result follows from a straightforward adaptation of the arguments given in the proof of Lemma~3 in \cite{kristensen2010characterizatio}.
	\begin{lemma} \label{lem:density}
		There exists a countable family of non-negative integrands $\{\phi_m \otimes g_m\otimes h_m\}_{m\in \mathbb{N}} \subset \E^2(\Omega;E)$
		that separates $\Y^2(\Omega;E)$. That is, if $\bm{\nu}, \bm{\sigma}\in \Y^2(\Omega;E)$, then
		\[
		\ddprb{ f_m, \bm{\nu}} = \ddprb{ f_m, \bm{\sigma}} \;\; \text{for all $m \in \mathbb{N}$ \quad $\Rightarrow$ \quad $\bm{\nu} \equiv \bm{\sigma}$ \, in \, $\Y^2(\Omega;E)$}.
		\]
		Moreover the family can be chosen so that each $h_m$ is uniformly Lipschitz on $E$ and each $g_m$ is uniformly continuously differentiable.
	\end{lemma}

	\subsection{Generating sequences}

	Vector-valued Radon measures can be naturally identified with a Young measure via the (compact) embedding
	\[
	\mathcal M(\Omega;E) \hookrightarrow \Y^2(\Omega;E) : \mu \mapsto \bm{\nu}_\mu \coloneqq (\delta_{\mu^{ac}},|\mu^s|,\mathscr L^d_Z,\delta_{\mu_s})\,.
	\]

	Given an integrand $f \in \mathbf{R}^2(\Omega;E)$ and a positive real $\epsilon$, we may define a functional
	\[
	I^\epsilon_f(\mu) \coloneqq \int_{\Omega} f(x,x/\epsilon,\mu^{ac}(x)) \;\mathrm{d} x
	+  \int_{\Omega } f^\infty(x,x/\epsilon,\mu_s(x)) \;\mathrm{d} |\mu^s|(x)\,,
	\]
	defined on vector-valued Radon measures $\mu \in \mathcal M(\Omega;\mathbb{R}^N)$.
	\begin{definition}[generating sequence]
		Let $\epsilon \searrow 0$ be an infinitesimal sequence of positive real numbers. We say that a sequence $(\mu_\epsilon)_\epsilon \subset \mathcal M(\Omega;\mathbb{R}^N)$ \emph{generates} the two-scale* Young measure $\bm{\nu} \in \Y^2(\Omega;E)$ if and only if
		\begin{align*}
		I^\epsilon_f(\mu_\epsilon)  \to  \ddprb{f,\nu} \quad \text{for all $f \in \E^2(\Omega;E)$\,.}
		\end{align*}
		In the case the domain of convergence is understood we simply write $\mu_\epsilon \overset{\Ybf^2}{\to} \bm{\nu}$.
	\end{definition}

	%

	%

	%

	\subsection{Proof of Theorem~\ref{thm:representation}}
	We are now ready to give the proof of Theorem~\ref{thm:representation} which asserts that every uniformly bounded sequence of measures generates (up to a subsequence) a two-scale* Young measure. It is worthwhile to mention the proof is not an immediate consequence of the compactness of Young measures (Corollary~\ref{cor:Compactnessofyoungmeasures}) since the microscopic variable \enquote{$x/\epsilon$} does not appear in the bi-linear pairing $\ddprb{\,\sbullet\,,\,\sbullet\,}$. Instead, the argument relies on the careful inspection of the limiting two-scale* Young measures, Proposition~\ref{lem:representationformula}, and the topological equivalence of Remark~\ref{rem:iso}.

	\begin{proof} In the context of the notation introduced above, we may re-formulate the statement of Theorem~\ref{thm:representation}  as follows: let $(\mu_\epsilon)\subset \mathcal M(\Omega;E)$ a sequence of measures with uniformly bounded variation. Then, there exists  a two-scale* Young-measure $\bm{\nu} \in \Y^2(\Omega;E)$ satsifying (up to a subsequence)
		\begin{equation}
		\label{eq:goal1}
		I_f^\epsilon(\mu_\epsilon) \to \ddprb{f,\bm{\nu}} \quad \text{for all $f \in \E^2(\Omega;E)$}.
		\end{equation}

		\proofstep{Step~1.}  Let $\delta > 0$. For an integrand $\Phi \in \mathrm{C}(\cl{\Omega \times Z \times \Bbb_E})$ define a continuous linear functional by setting
		\[
		L_\delta(\Phi)  \coloneqq I^\delta_{(S\Phi)}(\mu)\,.
		\]
		By the definition of the map $f \mapsto I^\delta_f$ with $\epsilon = \delta$ we obtain the estimate $| \langle L_\epsilon, \Phi \rangle | \leq \|\Phi\|_{\infty}( \mathscr L^d(\Omega)+ \|\mu_\epsilon\|)$ for every $\epsilon > 0$ of the infinitesimal sequence. The Banach--Alaoglu theorem and the Riesz representation theorem then yield the existence of a subsequence $(\varepsilon_k)_{k \in \mathbb{N}}$
		and a measure $\mu_L\in \mathcal M\big(\cl{\Omega \times Z \times \Bbb_E}\big)$ satisfying
		\[
		L_{\varepsilon_k} \overset{*}\rightharpoonup \dpr{\mu_L,\,\sbullet\,} \quad \text{as functionals in $\mathrm{C}(\cl{\Omega \times Z \times \Bbb_E})^*$}.
		\]

		Testing this convergence with an integrand of the form $\Phi =  \phi \otimes g \otimes (1 - |\,\sbullet\,|)$
		and using that $S\Phi = \phi \otimes g \otimes \chi_{E}$ we further deduce
		\[
		L_\epsilon(\Phi) = \int_{\Omega} \phi(x) g(x/\epsilon_k) \;\mathrm{d} x  \to  \dpr{L,\Phi}\,.
		\]
		This gives
		\[
		\int_{\cl{\Omega \times Z \times \Bbb_E}} \phi(x)g(\xi)(1 - |\hat z|) \;\mathrm{d} L(x,\xi,\hat z) = \int_{\Omega}\phi(x)\;\mathrm{d} x \int_Z  g(\xi) \;\mathrm{d} \xi
		\]
		for all $\phi \in \mathrm{C}(\cl \Omega)$ and $g \in \mathrm{C}(Z)$.
		We apply Lemma~\ref{lem:representationformula} to the measure $L$, to find a two-scale* Young measure  $\bm{\nu} \in \Y^2(\Omega;E)$ with $S^*\bm{\nu} = L$.

		\proofstep{Step~2.} Using the commutative diagram~\eqref{diagram} and the identity {$S \circ T = \id_{\E^2(\Omega;E)}$} we conclude that
		(recall that  $Tf(x,\xi,\,\sbullet\,) \equiv f^\infty(x,\xi,\,\sbullet\,)$ as functions over $\partial \Bbb_E$)
		\[
		I^\epsilon_f(\mu_\epsilon) = L_\epsilon(Tf) \to \dpr{L,Tf} = \dpr{S^*\bm{\nu},Tf} = \ddprb{\bm{\nu},f} \qquad \text{for all $f \in \E^2(\Omega;E)$}\,.
		\]
		This proves~\eqref{eq:goal1}. \end{proof}

	%

	%

	\subsection{Barycenter measures and two-scale convergence}

	We now turn to the concept of two-scale convergence. Following~\cite{nguetseng1989a-general-conve} we extend this notion to the two-scale convergence of measures as follows.

	\begin{definition}[two-scale convergence]\label{def:tsc} Let $\epsilon \searrow 0$ be a sequence of infinitesimal real numbers. Let also $\kappa \in \mathcal M^+(\cl \Omega)$ be a positive measure, and $\theta$ be a {weak-$*$} $\kappa$-measurable map from $\cl \Omega$ into $\mathcal M(Z;E)$.  We say that the sequence of measures $(\mu_\epsilon)_\epsilon \subset \mathcal M(\Omega;E)$ two-scale converges (as $\epsilon \downarrow 0$) to the generalized product measure $\mu = \kappa \otimes \theta_x$ if and only if

		\[
		\int_\Omega \Psi\bigg(x,\frac{x}{\epsilon}\bigg)\;\mathrm{d} \mu_\epsilon(x)  \to \int_{\cl \Omega}\bigg(\int_{Z} \Psi(x,\xi) \;\mathrm{d} \theta_x(\xi) \bigg) \;\mathrm{d} \kappa(x) \quad \text{for all $\Psi \in \mathrm{C}(\cl \Omega \times Z)$}.
		\]

	\end{definition}

	This {limit} concept is linked to the notion of \emph{barycenter} and \emph{second-scale barycenter} of a two-scale* Young measure (defined below) which will be significant for our techniques.

	\begin{definition}[barycenter]\label{def:bary}
		Let $\bm{\nu}  = (\nu,\lambda,\rho,\nu^\infty) \in \Y^2(\Omega;E)$. We define the \emph{barycenter} of $\bm{\nu}$ to be the $E$-valued  measure in $\mathcal M (\cl\Omega; E)$ defined as
		\begin{align*}
		[\bm{\nu}] \coloneqq \bigg(\int_{Z} \dprb{ \id_{E},\nu_{x,\xi} } \;\mathrm{d} \xi \bigg) \mathscr L^d_\Omega  +  \bigg( \int_{Z} \dprb{ \id_{\partial \Bbb_E},\nu^\infty_{x,\xi} } \;\mathrm{d} \rho_x(\xi)\bigg) \lambda\,.
		\end{align*}

	\end{definition}
	Let $\bm{\nu} = (\nu,\lambda,\rho,\nu^\infty) \in \Y^2(\Omega;E)$ be fixed. By the Radon--Nikodym decomposition theorem, applied to the measure $\lambda$, there exist a partition of $\cl \Omega = \mathrm{Reg}_{\bm{\nu}}(\Omega) \cup \mathrm{Sing}_{\bm{\nu}}(\Omega)$ by subsets satisfying
	%

	\[
	\mathscr L^d(\mathrm{Sing}_{\bm{\nu}}(\Omega)) = \lambda^s(\mathrm{Reg}_{\bm{\nu}}(\Omega)) = 0\,.
	\]
	Moreover,
	\[
	\mathrm{Sing}_{\bm{\nu}}(\Omega) = \setBB{x \in \cl \Omega}{\frac{\;\mathrm{d} \mathscr L^d}{\;\mathrm{d} |\lambda^s|}(x) = 0}.
	\]
	We are now in position to give the notion of second-scale barycenter.

	\begin{definition}[second-scale barycenter]\label{def:bary2}
		Let $\bm{\nu} = (\nu,\lambda,\rho,\nu^\infty) \in \Y^2(\Omega;E)$. The \emph{second-scale barycenter} of $\bm{\nu}$ is the weak-$\ast$ $(\mathscr L^d_\Omega + \lambda^s)$-measurable map $\bbpr{\bm{\nu}} : x \mapsto \bbpr{\bm{\nu}}_x \in \mathcal M(Z;E)$ where $\bbpr{\bm{\nu}}_x$ is the measure defined by the values on Borel subsets $V \subset Z$ as
		\[
		\bbpr{\bm{\nu}}_{x}(V) =
		\begin{cases} \displaystyle \int_V \dprb{\id_{E},\nu_{x,\xi}} \;\mathrm{d} \xi \; + \;  \lambda^{ac}(x) \int_V \dprb{\id_{\partial \Bbb_E},\nu_{x,\xi}^\infty} \;\mathrm{d} \rho_{x}(\xi) & \text{if
			$x \in \mathrm{Reg}_{\bm{\nu}}(\Omega)$}\\
		\\
		\displaystyle \int_V \dprb{\id_{\partial \Bbb_E},\nu_{x,\xi}^\infty} \;\mathrm{d} \rho_{x}(\xi) & \text{if $x \in \mathrm{Sing}_{\bm{\nu}}(\Omega)$}
		\end{cases}\,.
		\]
	\end{definition}
	The barycenter $[\bm{\nu}]$ can be recovered by integration from the second-scale barycenter. Indeed,
	$[\bm{\nu}] = \bbpr{\bm{\nu}}(Z) \, (\mathscr L^d_\Omega + \lambda^s)$ and

	in particular
	\begin{equation}\label{eq:density2}
	\bbpr{\nu}_x(Z) = \begin{cases}\displaystyle
	\frac{\;\mathrm{d} [\nu]}{\;\mathrm{d} \mathscr L^d}(x) & \text{$\mathscr L^d$-almost everywhere in $\Omega$}\\
	\\
	\displaystyle \frac{\;\mathrm{d} [\nu]}{\lambda^s}(x) & \text{$\lambda^s$-almost everywhere in $\cl \Omega$}\\
	\end{cases}\,.
	\end{equation}
	Moreover, if $\mu_\epsilon \overset{\Ybf}{\to}^2 \bm{\nu}$, then
	\begin{equation}\label{eq:1-2}
	\mu_\epsilon \overset{*}\rightharpoonup [\bm{\nu}] \quad \text{as measures on $\cl \Omega$}
	\end{equation}
	and
	\begin{equation}\label{eq:1-3}
	\text{$\mu_\epsilon$ two-scale converges to $(\mathscr L^d_\Omega + \lambda^s) \otimes \bbpr{\bm{\nu}}_x$}.
	\end{equation}

	\begin{corollary}[compactness of two-scale convergence]\label{cor:tsc} Let $\epsilon \searrow 0$ be an infinitesimal sequence of real numbers and let $(\mu_\epsilon)_\epsilon \subset \mathcal M(\Omega;E)$ be a sequence of measures with uniformly bounded total variation.

		Then there exists a subsequence $(\epsilon_k)_{k \in \mathbb{N}}$, a positive measure $\kappa \in \mathcal M(\cl \Omega)$, and a weak-$*$ $\kappa$-measurable map $\theta$ from $\cl \Omega$ into $\mathcal M(Z;E)$ such that
		\[
		\mu_{\epsilon_k} \; \text{two-scale converges to $\kappa \otimes \theta_x$}\,.
		\]
		Moreover $(\kappa \otimes |\theta_x|)(\cl \Omega \times Z) \le \liminf_{\epsilon \downarrow 0} |\mu_\epsilon|(\Omega)$.
	\end{corollary}
	\begin{proof}Apply Theorem~\ref{thm:representation} and~\eqref{eq:1-3} to the sequence $(\mu_\epsilon)_\epsilon$. The first conclusion follows by setting  $\kappa = \mathscr L^d_\Omega + \lambda^s$ and $\theta = \bbpr{\bm{\nu}}$. The lower semicontinuity of the norms follows from the fact that
		\[
		\mu_\epsilon \overset{*}\rightharpoonup \bpr{\bm{\nu}} = \theta(Z)\, \kappa\quad \text{as measures in $\cl \Omega$}\,.
		\]
	\end{proof}

	\begin{remark}[classical two-scale convergence] Notice that if a sequence of functions $(u_\epsilon)_\epsilon$ is equi-integrable (or if its is uniformly bounded in $\mathrm{L}^p$ for some $1 < p \le \infty$) and the sequence $u_\epsilon \mathscr L^d_\Omega$ generates a two-scale* Young measure $\bm{\nu} = (\nu,\lambda,\rho,\nu^\infty) \in \Y^2(\Omega;E)$, then $\lambda \equiv 0$ and hence
		\[
		\mu_\epsilon \; \text{two-scale converges to $u\, (\mathscr L^d_\Omega \otimes \mathscr L^d_Z)$}
		\]
		where
		\[
		u(x,\xi) \coloneqq \dpr{\id_E, \nu_{x,\xi}}\,.
		\]
		Thus the space of all such two-scale limits can be identified with $\mathrm{L}^1(\Omega \times Z;\mathbb{R}^N)$. In particular, our definition of two-scale convergence extends Nguetseng's original definition of two-scale convergence in $\mathrm{L}^p$-spaces \cite{nguetseng1989a-general-conve} (see also \cite{allaire1992homogenization-}).

	\end{remark}

	\subsubsection{Weighted barycenter measures}
	It will be often resourceful to interpret the pairing $\ddprb{\,\sbullet\,,\,\sbullet\,}$ as a measure over $\cl \Omega$. Let us introduce some additional notation by extending the definitions of barycenter and two-scale barycenter of a two-scale* Young measure $\bm{\nu} \in \Y^2(\Omega;E)$.

	\begin{definition}Let $f \in \mathbf{R}^2(\Omega;E)$. The $f$-\emph{barycenter} of $\bm{\nu} = (\nu,\lambda,\rho,\nu^\infty)$ is the vector-valued measure
		\begin{align*}
		[f,\bm{\nu}] \coloneqq  \bigg(\int_{Z} \dprb{f(x,\xi,\,\sbullet\,),\nu_{x,\xi} } \;\mathrm{d} \xi \bigg) \mathscr L^d_\Omega
		+ \bigg( \int_{Z} \dprb{f^\infty(x,\xi,\,\sbullet\,),\nu^\infty_{x,\xi} } \;\mathrm{d} \rho_x(\xi)\bigg) \lambda\,.
		\end{align*}
	\end{definition}
	Using this notation we get $\bpr{f,\bm{\nu}}(\cl \Omega) = \ddprb{f, \bm{\nu}}$, and $[\bm{\nu}]=\big(\bpr{(\id_{\mathbb{R}^N})_j,\bm{\nu}}\big)_{j=1,\dots,N}$.
	We also define a weighted \emph{second-scale $f$-barycenter}.

	\begin{definition}[second-scale $f$-barycenter] Let $f \in \mathbf{R}^2(\Omega;E)$. We define the $f$-\emph{barycenter} of $\bm{\nu} = (\nu,\lambda,\rho,\nu^\infty)$ as the weak-$*$ $(\mathscr L^d + \lambda^s)$-measurable map $\bbpr{f,\bm{\nu}} : x \mapsto \bbpr{ f, \bm{\nu}}_{x} \in \mathcal M(Z;{E})$ where is the measure defined by its values on Borel subsets $V \subset Z$ as

		\begin{equation}\label{def:weighted}
		\bbpr{ f, \bm{\nu}}_{x}(V) =
		\begin{cases}
		\displaystyle \int_V \dprb{f(x, \xi,\,\sbullet\,),\nu_{x,\xi}}\, \;\mathrm{d} \xi \\
		\displaystyle\qquad + \; \lambda^{ac}(x) \int_Z \dprb{f^\infty(x,\xi,\,\sbullet\,),\nu_{x,\xi}^\infty}\;\mathrm{d} \rho_{x}(\xi) \quad & \text{if $x \in \mathrm{Reg}_{\bm{\nu}}(\Omega)$}\\
		\\
		\displaystyle \int_V \dprb{f^\infty(x, \xi,\,\sbullet\,),\nu_{x,\xi}^\infty} \;\mathrm{d} \rho_{x}(\xi) & \text{if $x \in \mathrm{Sing}_{\bm{\nu}}(\Omega)$}
		\end{cases}.
		\end{equation}
	\end{definition}
	As we have seen before for the barycenter measures, a similar integral property holds for the weighted barycenters, namely $\bpr{f,\bm{\nu}} = \bbpr{f,\bm{\nu}}(Z) \, (\mathscr L^d_\Omega + \lambda^s)$ which in particular entails the identities
	\begin{equation}\label{eq:density3}
	\bbpr{f,\bm{\nu}}_x(Z) = \begin{cases}\displaystyle
	\frac{\;\mathrm{d} [f,\bm{\nu}]}{\;\mathrm{d} \mathscr L^d}(x) & \text{$\mathscr L^d$-almost everywhere in $\Omega$}\\
	\\
	\displaystyle \frac{\;\mathrm{d} [f,\bm{\nu}]}{\lambda^s}(x) & \text{$\lambda^s$-almost everywhere in $\cl \Omega$}\\
	\end{cases}\,.
	\end{equation}

	%

	%
	%

	%

	%

	%

	%

	%
	%

	%

	%

	%

	\subsection{Heuristics and some generic examples}

	Let us begin by recalling the criterion for $\mathrm{L}^1$-weak compactness due to Dunford \& Pettis: a sequence $(w_\epsilon)_\epsilon \subset \mathrm{L}^1(\Omega)$ is $\mathrm{L}^1$-weak relatively compact if and only if it is equi-integrable, that is, whenever
	\[
	\lim_{R \to \infty} \bigg(\limsup_{\epsilon \downarrow 0} \int_{\Omega \cap \{|w_\epsilon| \ge R\}} |w_\epsilon| \bigg)= 0.
	\]
	At those regions where a sequence fails to be equi-integrable, but its weak-$*$ limit remains absolutely continuous ---which we call \enquote{continuous concentration}, we speak of a \emph{biting limit} of the sequence:
	\begin{lemma}[Chacon]\label{chacon} Let $(w_\epsilon) \subset \mathrm{L}^1(\Omega)$ be a uniformly bounded sequence. Then there exist $w \in \mathrm{L}^1(\Omega)$, a subsequence $(w_{\epsilon_j})_{j \in \mathbb{N}}$, and a non-increasing sequence of measurable subsets $(K_m)_{m \in \mathbb{N}}$ of $\Omega$ such that
		\begin{enumerate}
			\item $\mathscr L^d(K_m) \to 0$ as $m \to \infty$,
			\item $w_{\epsilon_j} \rightharpoonup w$ in $\mathrm{L}^1(\Omega \setminus K_m)$  for all $m \in \mathbb{N}$.
		\end{enumerate}
		Provided that the these hold, we say that $w$ is the biting limit of $(w_{\epsilon_j})_{j \in \mathbb{N}}$, and the set
		\[
		B_u \coloneqq \setBB{x \in \Omega}{x \notin \bigcap_{m \in \mathbb{N}} K_m}
		\]
		is called the set of biting points of $(w_{\epsilon_j})_{j \in \mathbb{N}}$.
	\end{lemma}

	Given a sequence $(\mu_\epsilon)_\epsilon$ that generates a two-scale* Young measure $\bm{\nu} = (\nu,\lambda,\rho,\nu^\infty)$ it is possible to understand, qualitatively speaking, how each element of the four-tuple can be understood in terms of the generating sequence. In the following arguments we will assume with a slight abuse of notation that $(\mu_\epsilon = u_\epsilon)_\epsilon \subset \mathrm{L}^1(\Omega;E)$.
	A simple consequence of the representation of two-scale* Young measures is that $\lambda$ and $\rho_x$ are supported (with respect to the $x$-variable) at non-biting points of the sequence that generates $\bm\nu$. Moreover, this two measures carry the mass of the sequence $(|\mu|_\epsilon)_\epsilon$ for different length-scales:

	\begin{enumerate}[(i)]
		\item \emph{The measure $\lambda$ quantifies the limit mass carried by $(|\mathscr L^d_\Omega u_\epsilon|)_\epsilon$ in the set where it fails to be equi-integrable.}

		By Theorem 2.9 in~\cite{alibert1997non-uniform-int}, the biting limit of the sequence is the measure $\dprb{|\,\sbullet\,|,\tilde\nu{\;\mathrm{d} x}} \mathscr L^d_\Omega$ where $\upsilon = (\tilde\nu,\tilde \lambda,\tilde \nu^\infty) \in \Y^1(\Omega;E)$ is the Young measure generated by $(u_\epsilon)_\epsilon$. Hence, we conclude from the representation of two-scale*
		Young measures that
		\[
		|u_\epsilon \mathscr L^d_\Omega| - \dprb{|\,\sbullet\,|,\tilde\nu{\;\mathrm{d} x}} \mathscr L^d_\Omega \overset{*}\rightharpoonup \tilde \lambda \equiv \lambda \quad \text{as measures on $\cl\Omega$}.
		\]

		\item \emph{The probability measure $\rho_{x_0}$ quantifies, at a given non-biting point $x_0$ of $(|u_\epsilon \mathscr L^d_\Omega|)_\epsilon$, the homogenized  mass carried by the sequence about $x_0$ in the following sense: If $A \subset Z$, then
			\[
			\int_U \chi_A\bigg(\frac{x}{\epsilon}\bigg)|u_\epsilon| \;\mathrm{d} x \to \int_U \rho_x(A) \;\mathrm{d} \lambda(x) \quad \text{for all Borel subsets $U \subset \cl \Omega$}.
			\]}

		The argument is a direct consequence of the representation of two-scale* Young measures and point (i).
		%

	\end{enumerate}

	To give a better understanding of how the different components of the limiting Young measure are connected to the features of a limiting sequence, we present the following examples, that also emphasize the possibility of concentration for $\rho$ and $\lambda$.
	\begin{figure}[h]
		\captionsetup{justification=raggedright}
		\centering
		\begin{tikzpicture}
		\draw (1.5,1)--(1.5,0.95);
		\draw (-1.5,1)--(-1.5,0.95);
		\draw (-1.5,1)--(1.5,1);
		\node at (1.5,0.95)[anchor = north] {\scriptsize{$1$}};
		\node at (-1.5,0.95)[anchor = north] {\scriptsize{$-1$}};
		\draw (0,1)--(0,0.95);
		\node at (0,0.95) [anchor=north]{\scriptsize{$0$}};

		\draw[->] (0,1)--(0,3);

		\draw (0.5,1)--(0.5,0.95);
		\node at (0.5,0.95)[anchor=north]{\scriptsize{$\epsilon^\alpha$}};

		\draw (0,2.5)--(0.5,2.5)--(0.5,1);
		\node at (0,2.5)[anchor=east]{\scriptsize{$\epsilon^{-\alpha}$}};
		\fill[ gray, opacity =0.3] (0,1) rectangle (0.5,2.5);

		\node at (0,0.4){~};
		\end{tikzpicture}
		\hspace{0.8cm}
		\begin{tikzpicture}

		\draw (0,1)--(3,1);
		\draw (3,1)--(3,0.95);
		\node at (3,1)[anchor = north] {\scriptsize{$1$}};
		\draw (1.5,1)--(1.5,0.95);
		\node at (1.5,1)[anchor = north] {\scriptsize{$\frac 12$}};

		\draw (0,1)--(0,0.95);
		\node at (0,0.95) [anchor=north]{\scriptsize{$0$}};

		\draw[->] (0,1)--(0,3);

		\node at (0,2)[anchor=east]{\scriptsize{$\epsilon^{-1}$}};

		\node at (0.7,2)[anchor=south]{\scriptsize{$\epsilon^{2}$}};

		\draw (0,1)rectangle(0.1,2);
		\fill[ gray, opacity =0.3] (0,1)rectangle(0.1,2);
		\draw (0.3,1)rectangle(0.4,2);
		\fill[ gray, opacity =0.3] (0.3,1)rectangle(0.4,2);
		\draw (0.6,1)rectangle(0.7,2);
		\fill[ gray, opacity =0.3] (0.6,1)rectangle(0.7,2);
		\draw (0.9,1)rectangle(1,2);
		\fill[ gray, opacity =0.3] (0.9,1)rectangle(1,2);
		\draw (1.2,1)rectangle(1.3,2);
		\fill[ gray, opacity =0.3] (1.2,1)rectangle(1.3,2);

		\node at (0.75,1)[anchor=north]{\scriptsize{$\epsilon^{-1}$ many}};

		\node at (0,0.4){~};
		\end{tikzpicture}
		\hspace{0.8cm}
		\begin{tikzpicture}
		\draw (1.5,1)--(1.5,0.95);
		\draw (-1.5,1)--(-1.5,0.95);

		\draw (-1.5,1)--(1.5,1);
		\node at (1.5,1)[anchor = north] {\scriptsize{$1$}};
		\node at (-1.5,1)[anchor = north] {\scriptsize{$-1$}};
		\draw (0,1)--(0,0.95);
		\node at (0,0.95) [anchor=north]{\scriptsize{$0$}};

		\draw[->] (0,1)--(0,3);
		\draw (0.5,1)--(0.5,0.95);
		\node at (0.5,0.95)[anchor=north]{\scriptsize{$\frac 12 \epsilon$}};

		\draw (0.5,1)--(0.5,2.75)--(0.65,2.75)--(0.65,1);
		\node at (0,2.75)[anchor=east]{\scriptsize{$\epsilon^{-2}$}};
		\node at (0.595,2.75)[anchor=south]{\scriptsize{$\epsilon^{2}$}};
		\fill[ gray, opacity =0.3] (0.5,1) rectangle (0.65,2.75);

		\node at (0,0.4){~};

		\end{tikzpicture}

		\caption{Sketches for the functions in Examples~\ref{expl:Singularpart}(a),~\ref{expl:Singularpart}(b) and~\ref{expl:xi0isthere} respectively.} \label{fig:Expl}
	\end{figure}

	For the sake of simplicity we will discuss examples in dimensions $d = 1$ and $E = \mathbb{R}$; their correspondent versions to higher dimensions are easily constructed by adding invariant directional measures or by mimicking similar constructions along transversal directions. We cover generic examples for each of the qualitative \emph{first-scale/second-scale} scenarios:
	\begin{table}[h!]
		\centering
		\begin{tabular}{| c | c | c |}
			\hline
			macro-scale ($x$) & micro-scale ($\xi$) & correlation   \\ [0.5ex]
			\hline
			abs. cont. & abs. cont. & $\mu_\epsilon$ equi-integrable   \\
			\cline{2-3}
			& concentration/dissipation & scale of phenomena $\gg \epsilon$  \\
			\hline
			concentration/dissipation & abs. cont. & scale of phenomena $\le \epsilon$   \\
			\cline{2-3}
			& concentration/dissipation & scale of phenomena $\gg \epsilon$   \\   [1ex]
			\hline
		\end{tabular}
		\caption{Qualitative properties of the elements of a two-scale* Young measure.}
		\label{table:1}
	\end{table}

	\begin{example}[singularity on the macro-scale]\label{expl:Singularpart}
		Fix $\alpha > 0$ and let $\Omega = (-1,1)$. We consider the family of functions

		\[
		u_{\varepsilon}(x) =  \varepsilon^{-\alpha} \chi_{(0,\varepsilon^\alpha)}(x), \quad \epsilon > 0.
		\]

		Let us assume, up to taking a sequence of $\epsilon$'s that $u_\epsilon \overset{\Ybf^2}{\to} \bm{\nu} = (\nu,\lambda,\rho,\nu^\infty)$ (in this case every subsequence generates the same Young measure). The following observations are easy to check:
		\begin{enumerate}
			\item \emph{Pure concentration in the $x$-variable.} On the one hand $|u_\varepsilon| \mathscr L^1 \overset{*}\rightharpoonup \delta_0$. On the other hand, using $f(x,y,z)= |z|$ as a test function the Young measure representation yields
			\[
			|u_{\varepsilon_k}|  \mathscr L_\Omega^d \overset{*}\rightharpoonup \bigg(\int_{Z}\langle| \,\sbullet\, | ,\nu_{x,\xi} \rangle \;\mathrm{d}\xi\bigg) \mathscr L_\Omega^d + \lambda\,,
			\]
			whence we deduce $\nu_{x,\xi} = \delta_0$ for $\mathscr L^1 \otimes \mathscr L^1$-almost every $(x,\xi)$ in $[-1,1] \times Z$ and $\lambda \equiv \delta_0$.  Thus, it suffices to characterize $\bm\nu$ at $x = 0$.
			\item Testing with an integrand of the form $f(x,\xi,z) = \phi(x)g(\xi)|z|$ we see through a change of variables that
			\begin{align}
			\int_{-1}^{1} \phi(x) g(x/\varepsilon)h(u_\varepsilon(x))  \;\mathrm{d} x
			& = \int_{0}^{\epsilon^\alpha} \varepsilon^{-\alpha}\phi(x)g(x/\varepsilon)h(1)\;\mathrm{d} x	 \nonumber\\
			& = \int_0^{1} \phi(\varepsilon^\alpha y )g(\varepsilon^{\alpha-1} y)h(1) \;\mathrm{d} y  \label{eq:explsingulartorus}  \rightarrow c(\alpha)\,.

			\end{align}
			\begin{enumerate}
				\item \emph{Pure singularity on the micro-scale:} if $\alpha > 1$,  the Young measure representation and the limit above give
				\[
				\int_Z g(\xi)\dpr{h,\nu_{0,\xi}^\infty} \;\mathrm{d} \rho_0(\xi)  = c(\alpha) =  g(0)h(1) \quad \Rightarrow.
				\]
				In conclusion $\rho_0 = \delta_0$, $\nu_{0,0}^\infty = \delta_1$, and
				\[
				\bm{\nu} = (\delta_0,\delta_0,\delta_0,\delta_1)\,.
				\]
				\item \emph{Absolute continuity in the second-scale:} if, in turn $\alpha \le 1$, we get
				\[
				\int_Z g(\xi)\dpr{h,\nu_{0,\xi}^\infty} \;\mathrm{d} \rho_0(\xi) = c(\alpha) = \phi(0)h(1)\int_Z g(\xi) \;\mathrm{d} \xi.
				\]
				Therefore $\rho_0 = \mathscr L^1_Z$ and $\nu_{0,0}^\infty = \delta_1$, which altogether yields
				\[
				\bm{\nu} =  (\delta_0,\delta_0,\mathscr L^1_Z,\delta_1)\,.
				\]

			\end{enumerate}
		\end{enumerate}
	\end{example}
	\begin{remark}
		Observe that when oscillations of the sequence $(u_\varepsilon)_\epsilon$ happens at a coarser length-scale than $\{\epsilon\}$, then the the two-scale Young measure does not provide more information than the classical Young measure. On the other hand, if oscillations of $(u_\varepsilon)_\epsilon$  takes part at a finer scale than $\{\varepsilon\}$, then the two-scale* Young measure cannot be recovered from the classical Young measure.
	\end{remark}

	\begin{example}[non-biting limit] Fakir's construction (also known as \emph{Fakir's carpet}) provides a good way to produce continuous concentrations. The idea is to create \emph{many} small concentrations which \emph{diffuse} before each of them can actually gain mass.

		Let $\Omega = (0,1)$. In order to showcase the sensitivity of Young measures we shall consider a \enquote{bi-directional half-carpet} which is generated by the functions
		\[
		u_{\frac 1k}(x) \coloneqq k \sum_{i=0}^{k/2} (-1)^i\chi_{[\frac ik, \frac ik + \frac 1{k^2}]}(x), \quad  k = 2,4,6,\dots
		\]
		Similarly to the example above, let us already assume that $u_{\frac 1k} \overset{\Ybf^2}{\to} \bm\nu$. The following observations can be deduced directly from the representation of two-scale* Young measures:

		\begin{enumerate}
			\item \emph{First-scale analysis:} $|u_\varepsilon|\, \mathscr L^1_\Omega\overset{*}\rightharpoonup \mathscr L^1\restrict [0,1/2]$. There exist no singular points in the first variable which is encoded by the equality of measures $\lambda^s \equiv 0$. However, every point in the interval $[0,1/2]$ is a non-biting point of the sequence $(u_{\epsilon_k})_k$ (clearly the sequence fails to  be equi-integrable at any subset of this interval) and hence
			\[
			\lambda \equiv \mathscr L^1\restrict [0,1/2] \quad \text{as measures on $[0,1]$}\,.
			\]
			Moreover,
			\[
			x \mapsto \int_Z \dpr{\nu_{x,\xi},|\,\sbullet\,|} \;\mathrm{d} \xi \equiv 0 \quad \Rightarrow \quad \nu_{x,\xi} = \delta_0 \quad \text{$\mathscr L^1_\Omega \times \mathscr L^1_Z$-almost everywhere}.
			\]

			\item \emph{Second-scale analysis:} let $f = \phi \otimes g \otimes h \in \mathbf{R}^2(\Omega;\mathbb{R})$. The limit of the energies $I^{\epsilon_k}_f(u_{\epsilon_k})$ as $k \to \infty$ can be computed by Riemann-integral partial sums as

			\begin{align*}
			\int_{0}^{1} \phi(x) g(kx)h(u_{\frac 1k}(x))  \;\mathrm{d} x
			& = k \sum_{i=0}^{k/2} \int_{\frac i k}^{\frac i k + \frac 1{k^2}}\phi(x) g(kx)h((-1)^i) \;\mathrm{d} x\\
			& =  \sum_{i=0}^{k/2} \int_{i}^{i + \frac 1{k}}\phi(y/k) g(y)h((-1)^i) \;\mathrm{d} y \\
			& \sim {g(k^{-1})} \sum_{i=0}^{k/2}\phi(i/k) h((-1)^i) \;\mathrm{d} y \\
			& \rightarrow  g(0)\bigg(\frac{h(1)}{2}\int_0^\frac{1}{2} \phi(x) \;\mathrm{d} x + \frac{h(-1)}{2}\int_0^\frac{1}{2} \phi(x) \;\mathrm{d} x\bigg)\,.
			\end{align*}
			The representation of two-scale* Young measures and a density argument then give
			\[
			\bm{\nu} = (\delta_0,\mathscr L^1 \restrict [0,1/2],\delta_0,\frac 12 \delta_{-1} + \frac 12 \delta_1)\,.
			\]

			\begin{remark}
				As it can be seen from the representation of Young measures, at biting-points $x \in [0,1]$ of the sequence $\{u_\epsilon\}$, the correspondent probability measures $\rho_x$ must be the uniform measure $\mathscr L^1\restrict(0,1]$.
			\end{remark}

		\end{enumerate}

	\end{example}

	The precise representation of a Young measure does not only depend on  the generating sequence $u_\epsilon$ but is also strongly influenced by the speed of oscillation $\epsilon$.
	This is an interesting and important feature of the compactness, that will also play a role in the localization principles below, and is emphasized in the following example.
	\begin{example}[non-uniqueness] \label{expl:xi0isthere}
		Let again $d=N=1$ and $\Omega = (0,1)$. For a fixed $\epsilon > 0$ and arbitrary $c,d \in (0,1]$, we consider the function
		\[
		u_\epsilon (x) \coloneqq \frac 1{\epsilon^{2}} \chi_{(a + b\epsilon, a + b \epsilon +\epsilon^2)}(x).
		\]
		\emph{First-scale:} If $\epsilon \searrow 0$, then for every subsequence $\epsilon_1\searrow 0$ and
		\begin{equation}\label{eq:non}
		u_{\epsilon_1} \overset{\Ybf}{\to} \upsilon = (\delta_0,\delta_a,\delta_1).
		\end{equation}

		\noindent\emph{Second-scale:} We may assume that $\area{\frac a{\epsilon_1}} \to \xi_1 \in Z$. Here, we recall that $\area{x} \in Z$ stands for the equivalence class of $x \in \mathbb{R}$ in the one-dimensional torus $Z$. Testing with an integrand $f = \phi \otimes g \otimes h$ gives

		\begin{align*}
		\int_\Omega f(x,x/{\epsilon_1},u_{\epsilon_1}) \;\mathrm{d} x
		& =  \int_{a + b {\epsilon_1}}^{a + b  \epsilon_1 +\epsilon_1^2} \phi(x)g(x/\epsilon_1)h(\epsilon_1^{-2}) \;\mathrm{d} x \\
		& = \epsilon h(\epsilon_1^{-2})\int_{\frac{a}{\epsilon_1} + b}^{\frac{a}{\epsilon_1} +b + \epsilon_1} \phi(\epsilon_1 y)g(y) \;\mathrm{d} y\\
		& \sim  \phi(a)g\big(\area{{a}/{\epsilon_1}} + b + \epsilon_1\big) \cdot \epsilon_1^2 h(\epsilon_1^{-2})\\
		& \rightarrow \phi(a)g(\xi_1 + b)h^\infty(1).
		\end{align*}
		From this we infer that $u_{\epsilon_1} \overset{\Ybf^2}{\to} \bm{\nu}_1$, where
		\[
		\bm{\nu}_1 =  (\delta_0,\delta_a,\rho^1 = \delta_{\xi_1 + b},\delta_1).
		\]
		%

		Notice that $\xi_1$ does not depend on $\{u_{\epsilon_1}\}$, but solely on subsequence $\{\epsilon_1\}$.

		Hence, the choice of a different subsequences generates a range of two-scale* Young measures  (compare this with the uniqueness of \eqref{eq:non}). Indeed, if $\{\epsilon_2\}$ is another subsequence satisfying
		\[
		\xi_2 \coloneqq \lim_{\epsilon_2 \downarrow 0} \area{\frac a {\epsilon_2}} \neq \xi_1,
		\]
		then $u_{\epsilon_2} \overset{\Ybf^2}{\to} \bm{\nu}_2 = (\delta_0,\delta_a,\rho^2 = \delta_{\xi_2 + b},\delta_1)$.
		However, passing to a subsequence does not entirely \emph{forget}  in the sense that it is possible to relate $\bm\nu_1$ and $\bm \nu_2$ by a translation in the torus:
		\[
		\Gamma^{\xi_1}_\# \rho^1 \equiv \Gamma_\#^{\xi_2} \rho^2 \equiv \delta_b\,,
		\]
		where $\Gamma^\eta : \xi \mapsto \xi - \eta$ is a translation map in  $\mathbb{R}^d$ (which also determines a translation in the $d$-dimensional torus).
		In fact this translation in the second-scale also occurs at biting points of the sequence. However, this is not reflected in the two-scale* Young measure since $\mathscr L^d_Z$ is an invariant measure under translation, that is,
		\[
		\Gamma^\xi_\# \mathscr L^d_Z \equiv \mathscr L^d_Z\quad \text{for all $\xi \in Z$}.
		\]

	\end{example}%

	%

	%

	%

	%

	%

	%

	\section{Localization principles}\label{sec:5}

	In this section we treat the (measure theoretic) \emph{differentiation} of Young measures which confirms the observation that the convergence
	\[
	\mu_\epsilon \overset{\Ybf^2}{\to} \bm\nu \in \Y^2(\Omega;E)
	\]
	is in fact local (with respect to the macroscopic variable $x$).
	We show that at a point $x_0 \in \Omega$, the information carried by $\bm\nu$ can be recovered by simply looking at the homogeneous Young measures $\bm\sigma \in \Y^2(Q;E)$ generated by blow-ups of the generating sequence $(\mu_\epsilon)_\epsilon$ at $x_0$. With the mathematical thrust of introducing Young measures as a serving tool, we establish localization principles at both continuity and singular points.
	Next, we recall some facts about the push-forward of blow-up maps on measures.

	Throughout this section we will indistinctly use the zero-extension map $\mathcal M(\Omega;E) \hookrightarrow \mathcal M_\mathrm{loc}(\mathbb{R}^d;E)$ to identify measures defined on $\Omega \subset \mathbb{R}^d$ with measures defined on the whole space $\mathbb{R}^d$.
	Fix $\mu \in \mathcal{M}(\Omega;\mathbb{R}^N)$ and consider the map $T^{(x_0,r)}(x) := (x - x_0)/r$,
	which blows-up $B_r(x_0)$, the open ball around $x_0 \in \Omega$ with radius $r > 0$, into the open unit ball $B_1 \subset \mathbb{R}^d$. Following the definition of push-forward, it is easy to check that
	\[
	T^{(x_0,r)}_\# \mu(B) := \mu(x_0+rB ) \quad \text{whenever $B \subset \frac 1r(\Omega-x_0)$ is a Borel set.}
	\]
	A simple calculation shows that the Radon--Nykodym decomposition of $T^{(x_0,r)}_\#\mu$ can be re-written in terms of $\mu$ as
	\begin{equation}\label{eq:blow_r}
	T^{(x_0,r)}_\#\mu = a \, \mathscr L^d + s \, T_\#^{(x_0,r)} |\mu^s|,
	\end{equation}
	where the densities $a \in \mathrm{L}^1_\mathrm{loc}(\mathbb{R}^d;\mathbb{R}^N)$ and $s \in \mathrm{L}^\infty_{T_\#^{(x_0,r)} |\mu^s|}(\mathbb{R}^d;\mathbb{S}^{N-1})$ are defined by the rules
	\begin{align*}
	a(y) & \coloneqq r^d\frac{\;\mathrm{d} \mu}{\;\mathrm{d} \mathscr L^d}(x_0 + ry), \\
	s(y)& \coloneqq \frac {\;\mathrm{d}\mu}{\;\mathrm{d} |\mu^s|} (x_0+ry).
	\end{align*}

	Retaking the notation of the last section, we write $\Gamma^{\xi_0}:\xi \mapsto \xi - \xi_0$ to denote the translation in $\mathbb{R}^d$ by a vector $\xi_0 \in \mathbb{R}^d$. To avoid a more intricate
	notation, we shall also write $\Gamma^{\xi_0}$ (with $\xi_0 \in [0,1)^d$) to denote the same translation when restricted to $Z$, that is, $\Gamma^{\xi_0} : Z \to Z : {\xi} \to \area{\xi - \xi_0}$.
	In this way, the push-forward action $\Gamma_\#^{\xi_0} : \mathcal M(Z) \to \mathcal M(Z)$ defines an automorphism of spaces.
	In particular, if $g\in L^1_\rho(Z)$, then
	\begin{align}\label{eq:translationinvariance}
	\int_{Z} g(\xi - \xi_0)   \;\mathrm{d} \rho(\xi)   =  \int_{Z}  g(\xi)   \;\mathrm{d} (\Gamma^{\xi_0}_{\#}\rho)(\xi).
	\end{align}
	Notice that if $\rho$ is a uniform measure in $Z$, then $\rho$ is translation invariant:
	\[
	\Gamma^{\xi_0}_{\#}\rho \equiv \rho\quad \text{for all $\xi_0 \in Z$)}.
	\]

	To avoid any possible confusion we shall write $y \in Q$ (or $\mathbb{R}^d$) to denote the blow-up variable; we keep the notation $x \in \Omega$ for the macroscopic scale.

	\subsection{Localization at regular points}

	\begin{lemma}\label{lem:locatregular}
		Let $\bm{\nu} = (\nu,\lambda,\rho,\nu^\infty) \in \Y^2(\Omega;E)$ be a two-scale* Young measure generated by a sequence  $(\mu_\epsilon)_\epsilon \subset \mathcal M(\Omega;E)$.
		Then, at $\mathscr L^{d}$-almost every $x_0\in \Omega$ there exists  (up to a translation in $Z$) a regular tangent two-scale* Young measure ${\bm{D\nu}}$ $= (D\nu,D\lambda,D\rho,D\nu^\infty) \in \Y^2(Q;E)$ of $\bm\nu$ at $x_0$. That is, there exists vector $\xi_0 = \xi(x_0)\in Z$ such that
		\begin{align}
		D\lambda & = \lambda^{ac}(x_{0}) \, \mathscr L^d_Q \in \Tan_1(\lambda,x_0),\\
		D\rho_y & = \Gamma^{\xi_0}_{\#}\rho_{x_0} \quad \text{for $D \lambda$-almost every $y\in \mathbb{R}^d$}.
		\end{align}
		Moreover,  $\{D\nu_{y,\xi}\},\{D\nu^\infty_{y,\xi}\}$ are homogeneous Young measures in the sense that
		\begin{align}
		D\nu_{y,\xi} & =  \nu_{x_0,\xi + \xi_0}\quad \text{for $(\mathscr L^d \otimes \mathscr L^d_Z)$-almost every $(y,\xi) \in \mathbb{R}^d \times Z$},\\
		D\nu^\infty_{y,\xi} & = \nu^\infty_{x_0,\xi + \xi_0}\quad \text{ for $(D\lambda \otimes D\rho_y)$-almost every $(y,\xi) \in \mathbb{R}^d \times Z$}.
		\end{align}

	\end{lemma}

	%

	\begin{proof}

		Let $\{f_m \coloneqq \phi_m \otimes g_m \otimes h_m \,|\,m\in \mathbb{N}\} \subset \E^2(Q;E)$ be the restriction to $\cl Q \times Z \times E$ of the dense subset provided by Lemma~\ref{lem:density}; without loss of generality assume that $g_1 \equiv \chi_Z \in \mathrm{C}(Z)$. Let also $\{\xi_k\}_{k \in \mathbb{N}}$ be a countable dense subset of $Z$.

		\emph{Step~1. Selection of regular points.}
		%
		%

		First, let us define a set of full $\mathscr L^d$-measure where we aim to show the assertions of the lemma. To do this we first list three Lebesgue-type properties which are satisfied for $\mathscr L^d$-almost every $\tilde x \in \Omega$: \\
		\textit{Lebesgue property 1.} Lebesgue points of the measure $\lambda$, that is,
		\begin{equation}\label{eq:x_01}
		\lim_{r \downarrow 0} \frac{\lambda(Q_r(\tilde x))}{r^d} = \lambda^{ac}(x_0), \quad \lim_{r \downarrow 0} \frac{\lambda^s(Q_r(\tilde x))}{r^d} = 0.
		\end{equation}
		\textit{Lebesgue property 2.} Lebesgue points of the map
		\begin{align}\label{eq:x_03}
		x \mapsto \int_Z  \langle |\,\sbullet\,|_E, \nu_{x,\xi} \rangle \;\mathrm{d} \xi.
		\end{align}
		\textit{Lebesgue property 3.} Lebesgue points of the family of weighted barycenter measures  $\{\bpr{f_{k,m}, \bm\nu}\}_{k,m \in \mathbb{N}}$ where $f_{k,m} \coloneqq \phi_m \otimes g_m \circ \Gamma^{\xi_k}\otimes h_m$. Recall from~\eqref{eq:density3} that {being a Lebesque point for all elements of the family} is equivalent to
		\begin{align}\label{eq:x_02}
		\frac{\;\mathrm{d} [f_{k,m},\bm\nu]}{\;\mathrm{d} \mathscr L^d}(\tilde x) = \bbpr{f_{k,m},\bm\nu}_{\tilde x}(Z)
		\quad \text{for all $k,m \in \mathbb{N}$}.

		\end{align}

		We shall show the conclusions of the lemma hold for all
		\[
		x_0 \in R \coloneqq \set{\tilde x \in \Omega}{\text{$\tilde x$ satisfies \eqref{eq:x_01}-\eqref{eq:x_02}}},
		\]
		which is a set of full $\mathscr L^d$-measure in $\Omega$.

		\emph{Step~2. Blow-up sequence.} As before, we write $\langle z \rangle$ to denote the equivalence class of a vector $z \in \mathbb{R}^d$ in the $d$-dimensional torus $Z$. Since $Z$ is a compact manifold we may assume (up to passing to a  subsequence $(\mu_{\epsilon_i})_{i \in \mathbb{N}}$) that
		\begin{equation}\label{eq:limitpoint}
		\Big\langle \frac{x_0}{\epsilon_i} \Big \rangle \to  \xi_0 \in Z \quad \text{as $i \to \infty$}.
		\end{equation}
		Let
		$r_j \downarrow 0$ (with $r_1 = 1$) be an infinitesimal sequence of radii and consider, for fixed $j \in \mathbb{N}$,
		the blow-up sequence
		\[
		\gamma_{\delta_{i,j}} \coloneqq \frac 1{r_{j}^{d}}T_\#^{(x_0,r_{j})} \mu_{\epsilon_{i}} \in \mathcal M(Q;E), \quad i \in \mathbb{N},
		\]
		where
		\[
		\delta_{i,j} \coloneqq \frac{\epsilon_{i}}{r_{j}},\quad i \in \mathbb{N},
		\]
		is the readjusted blow-up length-scale sequence (this conforms again an infinitesimal sequence).

		Since $(\mu_{\epsilon_i})_{i \in \mathbb{N}}$ is uniformly bounded in $\mathcal M(\Omega;\mathbb{R}^N)$, we also have
		\[
		\sup_{i \in \mathbb{N}}\gamma_{\delta_{i,j}}(Q)
		< \infty
		\quad \text{for each $j \in \mathbb{N}$.}
		\]

		For $j = 1$, we use the compactness result in Theorem~\ref{thm:representation} to find a subsequence $\{1(i)\} \subset \{i\}$

		and $\bm{\sigma}^{(1)} \in \Y^2(Q;E)$
		such that
		\[
		\gamma_{\delta_{1(i),1}} \overset{\Ybf^2}{\to} \bm{\sigma^{(1)}} \quad \text{on $\cl Q$ \quad as \quad $i \rightarrow \infty$.}
		\]
		Recursively, for each $2 \le j \in \mathbb{N}$, we may find a sequence $\{j(i)\}\subset \{(j-1)(i)\}$ and a two-scale* Young measure $\bm{\sigma^{({j})}}\in \Y^2(Q;E)$ such that
		\[
		\gamma_{\delta_{j(i),j}} \overset{\Ybf^2}{\to} \bm{\sigma^{({j})}} \qquad \text{ as \quad $i \to \infty$.}
		\]

		\emph{Step 3. Characterization of $\bm{\sigma^{({j})}}$.}
		In this step we fix $j \in \mathbb{N}$.
		Let $f = \phi \otimes g \otimes h \in \E^2(Q;\mathbb{R}^N)$ with $g \in \mathrm{C}^1(Z)$ and $h$ uniformly Lipschitz.
		A change of variables and the decomposition in~\eqref{eq:blow_r} yields
		\begin{equation}\label{eq:loctest1}
		\begin{split}
		\ddprb{f, \bm{\sigma^{(j)}}} & = \lim_{i \to \infty}\bigg(
		\int_Q \phi(y)\,g(\,y/\delta_{j(i),j}\,)\,h(\,\mu_{\epsilon_{j(i)}}^{ac}(x_0 + r_jy)\,) \;\mathrm{d} y \\
		&+ \int_{\cl Q}  \phi(y)g\bigg(\frac{y}{\delta_{j(i),j}}\bigg)h^\infty\bigg(\frac{\;\mathrm{d} \mu_{\epsilon_{j(i)}}}{\;\mathrm{d} |\mu^s_{\epsilon_{j(i)}}|}(x_0 + r_jy)\,\bigg) \;\mathrm{d} \big( r_j^{-d}
		T_\#^{(x_0,r_j)}|\mu^s_{\epsilon_{j(i)}}|)(y)\bigg) \\
		& =  \frac 1{r_j^{d}}\lim_{i \to \infty} \bigg(\int_{Q_{r_j}(x_0)} \big(\varphi\circ T^{(x_0,r_j)}\big)(x)  \, g\bigg(\frac{x - x_0}{\epsilon_{j(i)}}\bigg)h(\,{\;\mathrm{d} \mu^{ac}_{\epsilon_{j(i)}}}(x)\,) \;\mathrm{d} x \\
		&+  \int_{\cl{Q_{r_j}(x_0)}}  \big(\varphi\circ T^{(x_0,r_j)}\big)(x)  \,g\bigg(\frac{x - x_0}{\epsilon_{j(i)}}\bigg)h^\infty\bigg(\frac{\;\mathrm{d} \mu_{\epsilon_{j(i)}}}{\;\mathrm{d} |\mu^s_{\epsilon_{j(i)}}|}(x)\bigg) \;\mathrm{d}|\mu^s_{\epsilon_{j(i))}} |(x)\bigg).

		\end{split}
		\end{equation}

		Setting $C \coloneqq \limsup_{\epsilon \downarrow 0}|\mu_{\epsilon}|(\overline{Q}) < \infty$ (from the original sequence) we may estimate the limiting behavior of the difference
		\begin{align*
		\bigg{|} \int_{\cl{Q_{r_j}(x_0)}} \phi \bigg(\frac{x-x_0} {r_j}\bigg)\bigg[g\bigg( \frac{x- x_0}{\varepsilon_{j(i)}}\bigg) - g\bigg(\frac{x}{\varepsilon_{j(i)}} - \xi_0\bigg)\bigg] \;\mathrm{d} |\mu_j|(x) \bigg|

		\end{align*}
		by
		\begin{align} \label{eq:xi0appears}
		C \|\phi\|_\infty \cdot \|Dg\|_\inft
		\cdot \bigg|\dprBB{\frac{x_0}{\epsilon_{{j(i)}}}} - \xi_0\bigg| = \mathrm{\textup{O}}(i) \to 0 \quad \text{as $i \to \infty$},
		\end{align}
		where to see that the last term vanishes as $i \to \infty$ we have used~\eqref{eq:limitpoint}.

		Using this  we may re-write~\eqref{eq:loctest1} as
		\begin{equation}\label{eq:loc2}
		\begin{split}
		\ddprb{f , \bm{\sigma^{(j)}}} & = \frac 1{r_j^{d}} \lim_{i \to \infty} \bigg(\int_{Q_{r_j}(x_0)} (\varphi \circ T^{(x_0,r_j)})(x) \, (g \circ \Gamma^{\xi_0}))\bigg(\frac{x}{\epsilon_{j(i)}}\bigg) \, h(\,{\;\mathrm{d} \mu^{ac}_{\epsilon_{j(i)}}}(x)\,) \;\mathrm{d} x \\
		& + \int_{\cl{Q_{r_j}(x_0)}}  (\varphi \circ T^{(x_0,r_j)})(x) (g \circ \Gamma^{\xi_0}))\bigg(\frac{x}{\epsilon_j}\bigg)h^\infty\bigg(\frac{\;\mathrm{d} \mu_{\epsilon_{j(i)}}}{\;\mathrm{d} |\mu^s_{\epsilon_{j(i)}}|}(x)\bigg) \;\mathrm{d}|\mu^s_{\epsilon_{j(i)}}| (x)\bigg) \\
		& = r^{-d}_j\ddprb{\varphi \circ T^{(x_0,r_j)} \otimes  g \circ \Gamma^{\xi_0} \otimes h, \bm{\nu}},
		\end{split}
		\end{equation}
		where in passing to the last equality we have used that $\mu_{\epsilon} \overset{\Ybf^2}{\to} \bm{\nu}$ in $\Y^2_\mathrm{loc}(\mathbb{R}^d;E)$.
		Applying  this to $|\,\sbullet\,|_{E}$ yields, together with (\ref{eq:x_01}) and (\ref{eq:x_03}),
		\begin{align}\label{eq:tris}
		\sup_{j \in \mathbb{N}} \ddprb{\phi \otimes g \otimes |\,\sbullet\,|_{E} , \bm{\sigma^{(j)}}}

		< \infty \quad \text{for all $\phi \in \mathrm{C}_c(\mathbb{R}^d)$ and $g \in \mathrm{C}^1(Z)$}\,.
		\end{align}
		We are then in position to apply the following global version of Corollary \ref{cor:Compactnessofyoungmeasures} (whose proof relies on a localization argument): there exists  a subsequence of $(r_j)_{j \in \mathbb{N}}$ (not relabeled) satisfying
		\[
		\text{$\bm{\sigma^{(j)}} \overset{*}\rightharpoonup \bm{D\nu(x_0)}$ \; in \; $\E^2(Q;E)^*$,  \; for some $\bm{D\nu(x_0)}\in  \Y^2(Q;E)$}\,.
		\]

		\emph{Step 4: Characterization of $\bm{\sigma}$.} Fix $m \in \mathbb{N}$ and let $j \in \mathbb{N}$ be an arbitrary positive integer.
		Let us write $\tilde f_m \coloneqq \phi_m \otimes g_m \circ \Gamma^{\xi_0} \otimes h_m$. For an arbitrary positive real number $\eta > 0$, we may use the uniform continuity of $g_m$ (recall that $g_m \in \mathrm{C}^1(Z)$) and the density of the set of points $\{\xi_k\}_{k \in \mathbb{N}}$ in $Z$ to find a sufficiently large $k = k(\eta)$ with the following property (here we use the positivity of $\phi_m$ and $h_m$):
		\begin{gather*}
		|\xi_0 - \xi_{k}| = \mathrm{\textup{O}}(\eta),
		\end{gather*}
		and
		\begin{align*}
		|\tilde f_m - f_{k,m}|(x,\xi,z) & \le \|g_m \circ \Gamma^{\xi_0} - g_m \circ \Gamma^{\xi_k} \|_\infty \cdot \phi_m(x)\cdot h_m(z) \\
		& \le \mathrm{\textup{O}}(\eta) \, [{\varphi \otimes \chi_Z \otimes |\,\sbullet\,|_E(x,\xi,z)}],
		\end{align*}
		where $\mathrm{\textup{O}}(\eta) \to 0$ as $\eta$ tends to zero. In particular,
		\begin{align*}
		\limsup_{j \in \mathbb{N}} \; & \Big|{r_j}^{-d} \bpr{\tilde f_m,\bm\nu}(\cl{Q_{r_j}(x_0)}) -  {r_j}^{-d} \bpr{f_{k(\eta),m},\bm\nu}(\cl{Q_{r_j}(x_0)})\Big| \\

		& \le \mathrm{\textup{O}}(\eta) \cdot \|\varphi_m\|_\infty \cdot \textrm{Lip}(h_m) \cdot \limsup_{j \in \mathbb{N}} r_j^{-d}\bpr{\phi_m \otimes \chi_Z \otimes |\,\sbullet\,|_E,\bm{\nu}}(\cl{Q_{r_j}(x_0)}) \\
		&  \lesssim_{(m)} \mathrm{\textup{O}}(\eta) \cdot \liminf_{j \to \infty} \bigg( 1 + \frac{\lambda(\cl Q_{r_j}(x_0))}{r_j^d} \bigg).
		\end{align*}

		Testing~\eqref{eq:loc2} with $f_m$ (not to be confused with $\tilde f_m$), it follows from the estimate above and~\eqref{eq:x_01} that
		\begin{equation}\label{eq:x_04}
		\begin{split}
		\ddprb{f_m, \bm{\sigma^{(j)}}}
		& = {r_j}^{-d} \ddprb {\phi_m \circ T^{(x_0,r_j)} \otimes g_m \circ \Gamma^{\xi_0} \otimes h_m, \bm{\nu}}  \\
		& = {r_j}^{-d} \bpr{\tilde f_m,\bm\nu}(\cl Q_{r_j}(x_0)) \\

		& = {r_j}^{-d} \bpr{f_{k(\eta),m},\bm\nu}(\cl Q_{r_j}(x_0)) + \mathrm{\textup{O}}(\eta).

		\end{split}
		\end{equation}

		%

		Letting $j \to \infty$ at both sides of~\eqref{eq:x_04},~\eqref{eq:x_02} and the weak-$*$ convergence $\bm{\sigma^{(j)}} \overset{*}\rightharpoondown \bm{D \nu(x_0)}$ imply

		\begin{align*}
		\ddprb{f_m,\bm{D \nu(x_0)}}  = \lim_{j \to \infty}  \ddprb{f_m, \bm{\sigma^{(j)}}}  = \bbpr{f_{k(\gamma),m},\nu}_{x_0}(Z) + \mathrm{\textup{O}}(\eta).

		\end{align*}
		Following analogous arguments to the ones in~\eqref{eq:x_04} and using $|\xi_0 - \xi_{k(\gamma)}| = \mathrm{\textup{O}}(\gamma)$, we may let $\eta \searrow 0$ to deduce\begin{equation}\label{eq:x_05}
		\begin{split}
		\ddprb{f_m,\bm{D \nu(x_0)}} & = \bbpr{\tilde f_m,\nu}_{x_0}(Z) \\
		& = \int_Q\bigg(\int_Z \dpr{f_m,\nu_{x_0,\xi_0 + \xi}} \;\mathrm{d}(\Gamma^{\xi_0}_\#\mathscr L^d_Z)(\xi) \\
		& \qquad + \frac{\;\mathrm{d} \lambda^\nu}{\;\mathrm{d} \mathscr L^d}(x_0)\int_Z\dpr{f_m^\infty,\nu_{x_0,\xi_0 + \xi}^\infty} \;\mathrm{d}(\Gamma^{\xi_0}_\#\rho_{x_0})(\xi)\bigg)\;\mathrm{d} y.
		\end{split}
		\end{equation}
		The sought assertion follows from the arbitrariness of $m$ on Step~4 (see~\eqref{eq:x_05}) and Lemma~\ref{lem:density} which translates into the equivalence
		\[
		\sigma \equiv \Big(\nu_{x_0,\xi_0 + \xi},\lambda^{ac}(x_0) \, \mathscr L^d,\Gamma^{\xi_0}_\#\rho_{x_0},\nu_{x_0,\xi_0 + \xi}^\infty\Big)_{y \in \cl Q, \xi \in Z}
		\]
		as two-scale* Young measures in the set $\Y^2(Q;\mathbb{R}^N)$. Notice we have used that the $d$-dimensional Lebesgue is uniformly distributed and hence $\Gamma_\#^{\xi_0} \mathscr L_Z^d \equiv \mathscr L_Z^d$. This proves the desired result.

		%

		%
		%

		%

	\end{proof}

	\subsection{Localization at singular points}

	\begin{proposition}\label{prop:loc_singular}Let $\bm{\nu} = (\nu,\lambda,\rho,\nu^\infty) \in \Y^2(\Omega;E)$ be two-scale* Young measure which is generated by a sequence of measures $(\mu_{\epsilon_i})_{i \in \mathbb{N}} \subset \mathcal M(\Omega;\mathbb{R}^N)$.
		Then, there exists a set $S \subset \Omega$ with full $\lambda^s$-measure that satisfies the following property:
		at every $x_0 \in S$ there exists (up to a translation in $Z$) a local tangent two-scale* Young measure $\bm{D\nu}= (D\nu,D \lambda,D \rho,D\nu^\infty)\in \Y^2(Q;\mathbb{R}^N)$ of $\bm\nu$ at $x_0$. That is, there exists $\xi_0 = \xi(x_0) \in Z$ such that
		\begin{align}
		D \lambda & \in \Tan_1(\lambda^s,x_0)\,, \\
		D\rho_y & = \Gamma^{\xi_0}_{\#}\rho_{x_0} \quad \text{for $D \lambda$-almost every $y\in Q$}\,.

		\end{align}
		Moreover,  $D\nu$ is concentrated on zero and $D\nu^\infty$ is homogeneous in the sense that
		\begin{align}
		D\nu_{y,\xi} & =  \delta_0 \qquad \quad \text{for $(\mathscr L^d_Q \otimes \mathscr L^d_Z)$-almost every $(y,\xi) \in Q \times Z$},\\
		D\nu^\infty_{y,\xi} & = \nu^\infty_{x_0,\xi + \xi_0} \quad  \text{for $\big(D \lambda \otimes D\rho_y\big)$-almost every $(y,\xi) \in Q \times Z$}.
		\end{align}

	\end{proposition}

	\begin{proof}
		The structure of the proof is very similar to the one of the localization principle for regular points. The main difference lies in the scaling of the blow up and, as a consequence of that, different terms  vanish.

		Let $\epsilon \searrow 0$ and $(\mu_{\epsilon}) \subset \mathcal{M} (\overline {\Omega};E)$ be such that $\mu_j \overset{\Ybf^2}{\to} \bm{\nu}$. We select a countable dense family $\{\xi_k\}_{k \in \mathbb{N}}$ in $Z$ and we write $\{f_m = \phi_m \otimes g_m \otimes h_m\,|\, m\in \mathbb{N} \} \subset \E^2(Q,E)$ to denote the restriction to $\cl Q \times Z \times E$ of the dense subset introduced in Lemma \ref{lem:density}. Additionally, this time we will assume without loss of generality that $h_1 \equiv \chi_{E}$. Recall we may assume the $h_m$'s to be Lipschitz continuous.

		\emph{Step~1. Selection of the singular set $S$.} We shall consider points $\tilde x \in \Omega$ satisfying the density estimate

		\begin{gather} \label{eq:x01}
		\lim_{r\downarrow 0}\displaystyle \frac { r^d + \int_{Q_r(\tilde x)} \int_Z \dpr{|\,\sbullet\,| \;\mathrm{d} \xi,\nu_{x,\xi}} \;\mathrm{d} \xi + \int_{Q_r(\tilde x)} \lambda^{ac}(x) \;\mathrm{d} x}{\lambda^s(Q_r(\tilde x))} =0\,,

		\end{gather}
		and which are
		\begin{equation}
		\text{$\lambda^s$-Lebesgue points of $\bpr{f_{k,m},\bm\nu} \in \mathcal M(\cl\Omega)$ \quad for all $k,m \in \mathbb{N}$}\,.
		\end{equation} Here, these barycenter measures are parametrized by the $x$-homogenous integrands
		\[
		f_{k,m} \coloneqq \chi_\Omega \otimes g_m\circ \Gamma^{\xi_k} \otimes h_m \quad k,m \in \mathbb{N}.
		\]
		In particular, by~\eqref{eq:density3} and the Lipschitz continuity of the $h_m$'s, at such points $\tilde x \in \Omega$ it holds
		\begin{align}\label{eq:x02}
		\frac{\;\mathrm{d} [f_{k,m},\bm\nu]}{\;\mathrm{d} \lambda^s}(\tilde x) = \bbpr{f_{k,m},\bm\nu}_{\tilde x}(Z) = \bbpr{f_{k,m}^\infty,\bm\nu}_{\tilde x}(Z)

		\quad \text{for all $k,m \in \mathbb{N}$}.

		\end{align}
		For the rest of the proof we fix a point $x_0 \in S \coloneqq \setb{\tilde x \in \Omega}{\text{$\tilde x$ satisfies \eqref{eq:x01}-\eqref{eq:x02}}}$, which is a set of full $\lambda^s$-measure in $\Omega$.

		\proofstep{Step~2. Blow-up sequence.} Since $Z$ is a compact manifold, we may again restrict to a subsequence of the generating sequence $(\mu_{\epsilon_i})_{i \in \mathbb{N}}$  and find $\xi_0 = \xi_0(x_0) \in Z$ such that
		\[
		\Big\langle \frac{x_0}{\epsilon_i} \Big \rangle \to \xi_0 \in Z \quad \text{as $i \to \infty$}.
		\]

		%

		%

		For a positive radius $r>0$,  we set $c_r  \coloneqq\lambda^s(Q_{r}(x_0))^{-1}$. By the compactness properties of measures we may find a weak* convergent blow-up sequence in the sense there exists a sequence of infinitesimal radii $(r_j)_{j \in \mathbb{N}}$ such that
		\begin{equation}\label{eq:Lambda}
		\Lambda_j \coloneqq c_{r_j} T^{(x_0,r_j)}_\# \lambda^s \overset{*}\rightharpoonup D\lambda \in \Tan_1(\lambda^s,x_0) \quad \text{and} \quad \gamma(\cl Q) = 1.
		\end{equation}

		For the sake of simplicity let us write $c_j \coloneqq c_{r_j}$. We consider, for fixed $j \in \mathbb{N}$, the $i$-indexed sequence
		\[
		\gamma_{i,j} \coloneqq c_{j_0}T_\#^{(x_0,r_{j_0})} \mu_{\epsilon_i}, \quad \text{where \;
			$\delta_{i,j} \coloneqq \frac{\epsilon_i}{r_{j}}$} \searrow 0\quad \text{as $i \to \infty$}\,.
		\]
		After an iterative procedure as the one for regular points, for each $j \in \mathbb{N}$ we may find a sequence $j(i)$ with the following properties: $\{(j+1)(i)\}_{i \in \mathbb{N}} \subset \{j(i)\}_{i \in \mathbb{N}}$ for all $j \in \mathbb{N}$, and
		\[
		\gamma_{j(i),j} \overset{\Ybf^2}{\to} \bm{\sigma^{({j})}} \; \text{for some $\bm{\sigma^{({j})}} \in \Y^2(Q;E)$}, \qquad j \in \mathbb{N}.
		\]
		Moreover, up to passing to a subsequence of $\{j\}_{j \in \mathbb{N}}$, we may assume there exists $\bm{D\nu(x_0)} \in \Y^2(Q;E)$ for which
		\begin{equation}\label{eq:limiting}
		\bm{\sigma^{({j})}} \overset{*}\rightharpoonup \bm{D\nu(x_0)} \quad \text{in $\E^2(Q;\mathbb{R}^N)^*$}.
		\end{equation}

		\emph{Step 3: Characterization of $\bm{\sigma^{({j})}}$.}
		For fixed $\phi \otimes g \otimes h \in \mathrm{C}(\overline {Q}) \times \mathrm{C}^1(Z) \times \E^0(\mathbb{R}^N)$ we deduce, by the same  change of variables used in (\ref{eq:loctest1}), that
		\begin{equation}
		\begin{split}
		\ddprb{f, & \bm{\sigma^{(j)}} }  = \lim_{i \to \infty}\bigg(
		\int_Q \phi(y)\,g({y}/{\,\delta_{j(i)}}\,)\,h\big(\,c_jr_j^d\mu^{ac}_{\epsilon_{j(i)}}(x_0 + r_jy)\,\big) \;\mathrm{d} y \\
		&+ \int_{\cl Q}  \phi(y)g\bigg(\frac{y}{\delta_{j(i)}}\bigg)h^\infty\bigg(\frac{\;\mathrm{d} \mu_{\epsilon_{j(i)}}}{\;\mathrm{d} |\mu^s_{\epsilon_{j(i)}}|}(x_0 + r_jy)\bigg) \;\mathrm{d} \big(c_{j}
		T_\#^{(x_0,r_j)}|\mu^s_{\epsilon_{j(i)}}|)(y)\bigg) \\
		& =  \lim_{i \to \infty} \frac 1{r_j^{d}}\bigg(\int_{Q_{r_j}(x_0)} \big(\varphi\circ T^{(x_0,r_j)}\big)(x)  \, g\bigg(\frac{x - x_0}{\epsilon_{j(i)}}\bigg)h\big(c_jr_j^d{\mu^{ac}_{\epsilon_{j(i)}}}(x)\big) \;\mathrm{d} x \\
		&+  \int_{\cl{Q_{r_j}(x_0)}}  \big(\varphi\circ T^{(x_0,r_j)}\big)(x)  \,g\bigg(\frac{x - x_0}{\epsilon_{j(i)}}\bigg)h^\infty\bigg(\frac{\;\mathrm{d} \mu_{\epsilon_{j(i)}}}{\;\mathrm{d} |\mu^s_{\epsilon_{j(i)}}|}(x)\bigg) \;\mathrm{d}(c_jr_j^d|\mu^s_{\epsilon_{j(i)}} |)(x)\bigg).

		\end{split}
		\end{equation}
		Making use of~\eqref{eq:x01} and the estimate~\eqref{eq:xi0appears} as it was used in~\eqref{eq:loc2}, and the fact that $\mu_{\epsilon_{j(i)}} \overset{\Ybf^2}{\to} \bm{\nu}$, we can re-write the limit on the right hand side in terms of a weighted barycenter measure as (here we use that $h^\infty(c \,\sbullet\,) = ch^\infty$ for all $c \ge 0$)
		\begin{equation}\label{eq:sloc2}
		\begin{split}
		\ddprb{f , & \bm{\sigma^{(j)}}}  =  \frac 1{r_j^{d}}\lim_{i \to \infty} \bigg(\int_{Q_{r_j}(x_0)} (\varphi \circ T^{(x_0,r_j)})(x) \, (g \circ \Gamma^{\xi_0}))\bigg(\frac{x}{\epsilon_{j(i)}}\bigg) \, h\big(c_jr_j^d{\mu^{ac}_{\epsilon_{j(i)}}}(x)\big) \;\mathrm{d} x \\
		& + \int_{\cl{Q_{r_j}(x_0)}}  (\varphi \circ T^{(x_0,r_j)})(x) (g \circ \Gamma^{\xi_0}))\bigg(\frac{x}{\epsilon_{j(i)}}\bigg)h^\infty\bigg(c_jr_j^d\frac{\;\mathrm{d} \mu_{\epsilon_{j(i)}}}{\;\mathrm{d} |\mu^s_{\epsilon_{j(i)}}|}(x)\bigg) \;\mathrm{d}|\mu^s_{\epsilon_{j(i)}}| (x)\bigg) \\
		& = r_j^{-d}[\varphi \circ T^{(x_0,r_j)} \otimes  g \circ \Gamma^{\xi_0} \otimes h(c_jr_j^d \,\sbullet\,),\bm{\nu}](\cl{Q_{r_j}(x_0)})
		\end{split}
		\end{equation}

		\proofstep{Step~4. Characterization of $\bm{D\nu(x_0)}$.} The last step consists on characterizing the limiting Young measure  in~\eqref{eq:limiting} as $j$ tends to infinity. Here, we expect the biting part of $\bm{\sigma}$ to vanish since we are performing a blow-up at a purely singular point $x_0 \in S$. Arguing as in Step~4 of the proof of the localization at regular points, we cast the $\mathrm{\textup{O}}(\eta)$-approximation of the dense family $\{\xi_k\}$ behind~\eqref{eq:x_04}; this time with the integrands $\tilde f_{k,m} \coloneqq \phi_m \circ T^{(x_0,r_j)} \otimes g_m \circ \Gamma^{\xi_k} \otimes h_m(c_j r_j^d \,\sbullet\,)$. For positive $\eta$ we find $k(\eta) \in \mathbb{N}$ sufficiently large so that that (using the density of the family $\{\xi_k\}_{k \in \mathbb{N}}$ in $Z$, and the uniform continuity of $g_m$)
		\[
		|\xi_0 - \xi_k(\eta)| = \mathrm{\textup{O}}(\eta) \quad \text{and} \quad |\tilde f_m - \tilde f_{k,m}| \le \mathrm{\textup{O}}(\eta) \,  [{\varphi \otimes \chi_Z \otimes (1 + c_jr_j^d|\,\sbullet\,|_E) (x,\xi,z)}],
		\]
		where $\tilde f_m  \coloneqq \phi_m \circ T^{(x_0,r_j)} \otimes g_m \circ \Gamma^{\xi_0} \otimes h_m(c_j r_j^d \,\sbullet\,)$.
		Using these estimates in~\eqref{eq:sloc2} we deduce by~\eqref{eq:x01} that
		\begin{equation}\label{eq:?}
		\begin{split}
		\ddprb{f_m, \bm{\sigma^{(j)}}}
		&

		= {r_j}^{-d} \bpr{\tilde f_{k(\eta),m},\bm\nu}(\cl{Q_{r_j}(x_0)}) + \mathrm{\textup{O}}(\eta).
		\end{split}
		\end{equation}

		For the next step let  $\eta$ be an arbitrary positive real. We claim  that taking then the limit as $j\to \infty$ at both sides of the equality yields
		\begin{equation}\label{eq:yacasi}
		\begin{split}
		\lim_{j \to \infty} \ddprb{f_m, \bm{\sigma^{(j)}}} & =  h_m (0)\int_Q \phi_m(y) \bigg(\int_Z g_m \circ \Gamma^{k(\eta)}(\xi) \;\mathrm{d} \xi \bigg) \;\mathrm{d} y  + \\ & \qquad  \bbpr{\tilde f_{k(\eta),m},\bm{\nu}}_{x_0}(Z)\cdot \int_{\cl{Q}}\phi_m(y) \;\mathrm{d} (D\lambda)(y)  +  \mathrm{\textup{O}}(\eta)\,.
		\end{split}
		\end{equation}
		First, we deal with the absolutely continuous part of the barycenter measure in~\eqref{eq:?} as follows. For all $\mathscr L^d$-almost every $x \in \cl{Q_{r_j}(x_0)}$ we estimate
		(recall that $h_1 \equiv \chi_{\mathbb{R}^N}$)

		\[
		\left| \frac{\;\mathrm{d} \bpr{\tilde f_{k(\eta),m} - \phi_m\circ T^{(x_0,r_j)} \otimes g_m \circ \Gamma^{\xi_{k(\eta)}} \otimes h_m(0)h_1,\bm\nu}}{\;\mathrm{d} \mathscr L^d}(x)\right| \le  c_jr_j^dC(m)\cdot \bbpr{|\,\sbullet\,|_E,\bm{\nu}}_{x}(Z),
		\]
		where $C(m) \coloneqq \|\phi_m\|_\infty\cdot\|g_m\|_\infty \cdot\textrm{Lip}(h_m)$. This estimate leads to the bound
		\begin{align*}
		\,\thetant-_{{Q_{r_j}(x_0)}}& \bigg|\frac{\;\mathrm{d} \bpr{\tilde f_{k(\eta),m} - \phi_m\circ T^{(x_0,r_j)} \otimes g_m \circ \Gamma^{\xi_{k(\eta)}} \otimes h_m(0)h_1,\bm\nu}}{\;\mathrm{d} \mathscr L^d}(x)\bigg| \;\mathrm{d} x\\& \lesssim_{(m)}
		\frac {\int_{Q_{r_j}(x_0)} \int_Z \dpr{|\,\sbullet\,|_E \;\mathrm{d} \xi,\nu_{x,\xi}} \;\mathrm{d} \xi + \int_{Q_{r_j}(x_0)} \lambda^{ac}(x) \;\mathrm{d} x}{\lambda^s(Q_{r_j}(\tilde x))} = \mathrm{\textup{O}}(r_j).
		\end{align*}
		Here, to reach the last equality we have used that $x_0 \in S$ satisfies~\eqref{eq:x01}. Thus, when passing to limit $j \to \infty$, we may substitute the absolutely continuous part of $r_j^{-d}\bpr{\tilde f_{k(\eta)}, \bm{\nu}}$ in~\eqref{eq:?} by the integrable function
		\[
		x \mapsto \frac{1}{r^d_j}\bigg( \int_Z\dprb{\tilde f_{k(\eta),m}(x,\xi,\,\sbullet\,),\delta_0} \;\mathrm{d} \xi \bigg)\,.
		\]
		We now deal with the singular  part and the passing to the limit. By~\eqref{eq:limiting} and the density identity from~\eqref{eq:x02} (applied to the point $\tilde x = x_0 \in S$), we may let $j$ tend to infinity at both sides of~\eqref{eq:yacasi} to deduce (recall that $h_m^\infty(c_jr_j^d \,\sbullet\,) = c_jr_j^dh_m^\infty$)
		\begin{align*}
		\ddprb{f_m,\bm{D\nu(x_0)}} &  =  \lim_{j \to \infty} \bigg(\,\thetant-_{Q_{r_j(x_0)}}\int_Z \dpr{\tilde f_{k(\eta),m}(x,\xi,\,\sbullet\,),\delta_0} \;\mathrm{d} \xi \;\mathrm{d} x \\
		& \qquad  +  \frac{1}{r_j^d}\int_{\cl{Q_{r_j}(x_0)}}\frac{\;\mathrm{d} \bpr{\tilde f_{k(\eta),m}^\infty,\bm\nu}}{\;\mathrm{d} \lambda^s}(x) \;\mathrm{d} \lambda^s(x) \bigg) + \mathrm{\textup{O}}(\eta) \\
		&  =  h_m (0)\int_Q \phi_m(y) \bigg(\int_Z g_m \circ \Gamma^{\xi_{k(\eta)}}(\xi) \;\mathrm{d} \xi \bigg) \;\mathrm{d} y  + \\ & \qquad  \lim_{j \to \infty} \,\thetant-_{\cl{Q_{r_j}(x_0)}} \phi_m(x)\frac{\bpr{\tilde f^\infty_{k(\eta),m},\bm{\nu}}}{\;\mathrm{d} \lambda^s}(x) \;\mathrm{d}\lambda^s(x) + \mathrm{\textup{O}}(\eta)\\
		& =  h_m (0)\int_Q \phi_m(y) \bigg(\int_Z g_m \circ \Gamma^{\xi_{k(\eta)}}(\xi) \;\mathrm{d} \xi \bigg) \;\mathrm{d} y  + \\ & \qquad  \bbpr{\tilde f_{k(\eta),m},\bm{\nu}}_{x_0}(Z)\cdot \int_{\cl{Q}}\phi_m(y) \;\mathrm{d} (D\lambda)(y)  + \mathrm{\textup{O}}(\eta)\,,
		\end{align*}
		where we have used $\Lambda_j = c_j T^{(x_0,r_j)}_\#\lambda^s \overset{*}\rightharpoonup D\lambda$ on $\cl Q$ in passing to the last equality. This proves the claim.

		To conclude we observe that $\eta$ has so far been chosen arbitrarily. Therefore, by the dominated convergence theorem, the identity~\eqref{eq:yacasi} implies
		\begin{equation}
		\begin{split}
		\ddprb{f_m,\bm{D\nu(x_0)}} & = h_m (0)\int_Q \phi_m(y) \bigg(\int_Z g_m \circ \Gamma^{\xi_0}(\xi) \;\mathrm{d} \xi \bigg) \;\mathrm{d} y \\ & \qquad +  \bbpr{\tilde f_{m},\bm{\nu}}_{x_0}(Z)\cdot\int_{\cl{Q}}\phi_m(y) \;\mathrm{d} (D\lambda)(y) \\
		& = \int_Q\bigg( \int_Z\dprb{f_m(x,\,\sbullet\,,\,\sbullet\,),\delta_0}  \;\mathrm{d} \xi\bigg)  \;\mathrm{d} y \\
		& \qquad  +\int_{\cl Q}\bigg(\int_Z\dpr{f_m^\infty(x,\,\sbullet\,,\,\sbullet\,),\nu_{x_0,\xi_0 + \xi}^\infty} \;\mathrm{d}(\Gamma^{\xi_0}_\#\rho_{x_0})(\xi) \bigg) \;\mathrm{d} (D\lambda)(y)\,.
		\end{split}
		\end{equation}
		Here, we have used that
		$\Gamma^{\xi_0}_\# \mathscr L^d_Z = \mathscr L^d_Z$ to reach the second equality. Since $\{f_m\}_{m \in \mathbb{N}}$ separates $\E^2(Q;\mathbb{R}^N)$,
		Lemma~\ref{lem:density} gives
		\[
		\bm{D\nu(x_0)} = (\delta_0,D\lambda,D\rho,D\nu^\infty), \quad D\lambda \in \Tan_1(\lambda^s,x_0),
		\]
		where $D\rho(x) = \Gamma^{\xi_0}_\#\rho_{x_0}$ and $D\nu^\infty(y,\xi) = \nu_{x_0,\xi_0 + \xi}^\infty$ for all $y \in \cl Q$. This finishes the proof.
	\end{proof}

	%

	%

	%

	%

	\section{PDE-constrained Young measures}\label{sec:5a}

	Let $k \in \mathbb{N}$ and let $E,F$ be finite-dimensional real vector spaces. For $\alpha \in \mathbb{N}^d$ we define its modulus $|\alpha| \coloneqq \alpha_1 + \dots + \alpha_d$. We shall consider general homogeneous differential operators of order $k$ on $\mathbb{R}^d$ from $E$ to $F$, that is, operators of the form
	\[
	\mathcal A = \sum_{\substack{\alpha \in \mathbb{N}^d\\|\alpha| = k}}  A_\alpha\,\partial^\alpha, \qquad A_\alpha \in \textrm{Lin}(E;F).
	\]
	A vector-valued measure $\mu \in \mathcal M(\Omega;E)$ is called $\mathcal A$-free provided that
	\[
	\mathcal A\mu = 0 \quad \text{in the sense of distributions on $\Omega$}.
	\]
	We say that $\bm\nu \in \Y^2(\Omega;E)$ is an $\mathcal A$-free two scale* Young measure if it is generated by a sequence of (asymptotically) $\mathcal A$-free measures:

	\begin{definition}[$\mathcal A$-free Young measures] Let $1 < p < \frac{d}{d-1}$. A two-scale* Young measure $\bm\nu \in \Y(\Omega;E)$ is called $\mathcal A$-free if there exist a sequences $\epsilon \searrow 0$ and $(\mu_\epsilon)_\epsilon \subset \mathcal M(\Omega;E)$ such that
		\[
		\mathcal A\mu_\epsilon \to 0 \;\; \text{strongly in $\mathrm{W}^{-k,p}(\Omega)$ \quad and \quad $\mu_\epsilon \overset{\Ybf^2}{\to} \bm\nu$}.
		\]

	\end{definition}

	%

	%

	%

	%

	%

	%

	%

	%
	%

	%
	%

	%

	%

	%

	%

	%

	%

	%

	%

	%

	%

	%

	%

	%

	%

	%

	\subsection{Rigidity properties of $\mathcal A$-free two-scale* Young measures}

	Clearly, the barycenter $[{\bm{\nu}}]$ of an $\mathcal A$-free two-scale* Young measure is $\mathcal A$-free. This same property is inherited to second-scale barycenters $\bbpr{\bm{\nu}}$ as it is portrayed in the next proposition. To deal with a possible abuse of notation about the domain of partial differential operators, we shall write $\mathcal A_\xi$ to denote the action of $\mathcal A$ on $\mathcal{D}'(Z;E)$, i.e.,
	\[
	(\mathcal A_\xi\eta)[g] = \dprb{\eta,\mathcal A^*g} \quad \text{for all $\eta \in \mathcal{D}'(Z;E)$ and $g \in \mathrm{C}^\infty_\mathrm{per}(Q;F)$}\,.
	\]

	\begin{proposition}\label{prop:secondbarycenterisafree} Let $\bm{\nu} = (\nu,\lambda,\rho,\nu^\infty) \in \Y^2(\Omega;E)$ be an $\mathcal A$-free two-scale* Young measure. Then at $(\mathscr L^d + \lambda^s)$-almost every $x_0 \in \Omega$ it holds that

		\[
		\mathcal A_\xi\bbpr{\bm{\nu}}_{x_0} = 0 \quad \text{in the sense of distributions on $Z$}\,.
		\]
	\end{proposition}
	\begin{proof}
		Let $\phi \in \mathrm{C}^\infty_c(\Omega)$ and let $g \in \mathrm{C}^\infty(Z)$ be arbitrary functions.
		Let $(\mu_\epsilon)_\epsilon$ be a sequence of $\mathcal A$-free measures that generates the Young measure $\bm{\nu}$ (on $\Omega$) and let $T^\varepsilon(x) = x/\varepsilon$ be the re-scaling by the factor $\epsilon$.
		As a consequence of the product rule there exist constants $c_{\alpha ,\beta}$ such that
		\[
		\mathcal A^*(\phi \cdot \epsilon^k (g\circ T^\epsilon)) = (-1)^k\sum_{\substack{|\alpha|\le k\\\beta \le \alpha}} \epsilon^{k - |\beta|}c_{\alpha ,\beta}A^T_\alpha  (\partial^{\alpha - \beta} \phi)(\partial^\beta g) \circ T^\epsilon,
		\]
		where we write $\beta \le \alpha$ if and only if $\beta_i \le \alpha_i$ for every $i = 1,\dots,d$.
		Observe that $c_{\alpha,\alpha} = 1$ for every $\alpha$, and in particular it follows that
		\[
		\| \mathcal A^*(\phi \cdot  \epsilon^k (g\circ T^\epsilon)) - \phi \cdot (\mathcal A^* g)\circ T^\epsilon\|_\infty \to 0, \quad \text{as $\epsilon \searrow 0$}.
		\]
		Hence, since $\mathcal A \mu_\epsilon \rightarrow 0$ in strongly in $\mathrm{W}^{-k,p}(\Omega) \overset{c}{\embed} \mathrm{C}_0(\Omega)^*$, we deduce from the convergence above that
		\begin{equation}\label{eq:discarde}
		\begin{split}
		0 & = \lim_{\epsilon \searrow 0} \int_\Omega \mathcal A^*[\phi \cdot \epsilon^k (g\circ T^\epsilon)](x) \;\mathrm{d} \mu_\epsilon(x) \\
		& =  \lim_{\epsilon \searrow 0} \int_\Omega \phi(x)(\mathcal A^*g)(x/\epsilon) \;\mathrm{d} \mu_\epsilon(x) \ \\
		& = \ddprb{\phi \otimes \mathcal A_\xi^*g \otimes \id_{E},\bm{\nu}}.
		\end{split}
		\end{equation}
		Using the tensor structure of the integrand and the property~\eqref{def:weighted} of the second-scale barycenters we obtain
		\begin{align*}
		0 & = \int_{\cl\Omega}\int_{Z}  \;\mathrm{d} \bbpr{\phi \otimes \mathcal A_\xi^* g \otimes \id_{\mathbb{R}^N},\bm{\nu}}_x(y) \;\mathrm{d} (\mathscr L^d + \lambda^s)(x) \\
		& = \int_{\Omega} \phi(x) \bigg(\int_Z (\mathcal A_\xi^*g)(\xi) \;\mathrm{d} \bbpr{\bm\nu}_x(\xi)\bigg) \;\mathrm{d} (\mathscr L^d + \lambda^s)(x).
		\end{align*}

		Since the choice of $\phi \in \mathrm{C}^\infty_c(\Omega)$ was arbitrary, this identity must hold locally. Namely, for $(\mathscr L^d + \lambda^s)$-a.e.~$x_0 \in \Omega$ it holds
		\[
		\dprb{g,\mathcal A_\xi \bbpr{\bm{\nu}}_{x_0}} =
		\int_{Z} (\mathcal A_\xi^*g)(\xi) \;\mathrm{d} \bbpr{\bm{\nu}}_{x_0}(\xi) =0 \quad \text{for all $g \in \mathrm{C}^\infty(Z)$}.
		\]
		By the distributional definition of derivative this is equivalent to
		\[
		\mathcal A_\xi \bbpr{\bm{\nu}}_{x_0} = 0 \quad \text{in the sense of distributions on $Z$}.
		\]
		This proves the assertion.
	\end{proof}

	\begin{corollary}[differential rigidity of the second-scale]\label{cor:Afree} Let $\epsilon \searrow 0$ and let $(\mu_\epsilon)_\epsilon \subset \mathcal M(\Omega;E)$ be a sequence of asymptotically $\mathcal A$-free measures that two-scale converges to a limit $\lambda \otimes \theta_x$. That is, such that
		\[
		\text{$\mathcal A\mu_\epsilon \to 0$ \;\; strongly in $\mathrm{W}^{-k,p}(\Omega)$}\,.
		\]
		Then, at $\lambda$-almost every $x_0 \in \Omega$, it holds
		\[
		\mathcal A_\xi \theta_{x_0} = 0 \quad \text{in the sense of distributions on $Z$.}
		\]
	\end{corollary}
	\begin{proof}
		This follows directly from~\eqref{eq:1-3} and the previous proposition.
	\end{proof}

	\begin{corollary}[structure of $\mathcal A$-free two-scale* Young measures]\label{cor:GuidoFilip}
		Let $\bm{\nu}\in \Y^2_\mathcal A(\Omega;E)$ be an $\mathcal A$-free two-scale* Young measure. Then, at $\lambda^s$-almost every $x_0 \in \Omega$, the following differential inclusion holds:
		\[
		\frac{\;\mathrm{d} \bbpr{\bm{\nu}}_{x_0}}{\;\mathrm{d} |\bbpr{\bm{\nu}}_{x_0}^s|}(\xi) \in \Lambda_{\mathcal A} \quad \text{for $|\bbpr{\bm{\nu}}_{x_0}^s|$-almost every $\xi \in Z$}.
		\]
	\end{corollary}
	\begin{proof}
		By Proposition \ref{prop:secondbarycenterisafree} we have that for $(\mathscr L^d + \lambda^s)$-almost every $x_0 \in \Omega$, the second-scale barycenter  $\bbpr{\bm{\nu}}_{x_0}$ is an $\mathcal A_\xi$-free measure on $Z$. Since locally, $\mathcal A_\xi$ and $\mathcal A$ coincide as operators in $\mathbb{R}^d$, it follows from \cite[Theorem~1.1]{de-philippis2016on-the-structur} that
		\[
		\frac{\;\mathrm{d} \bbpr{\bm{\nu}}_{x_0}}{\;\mathrm{d} |\bbpr{\bm{\nu}}_{x_0}^s| }(\xi)    \in \Lambda_{\mathcal A} \text{ for } |\bbpr{\bm{\nu}}_{x_0}^s| \text{-a.e. $\xi \in Z$.}
		\]
		This finishes the proof.

	\end{proof}

	The next lemma asserts the support of the purely singular part of an $\mathcal A$-free measure cannot be arbitrary. In fact, it must be contained in the smallest vectorial space containing the wave cone $\Lambda_{\mathcal A}$.

	\begin{lemma}\label{lem:support} Let $\bm{\nu} = (\nu,\lambda,\rho,\nu^\infty) \in \Y^2(\Omega;E)$ be an $\mathcal A$-free two-scale* Young measure.

		Then the support of the purely singular part of $\bm\nu$ is contained in $\spn \{\Lambda_\mathcal A\}$, that is,
		\[
		\supp \, (\nu^\infty_{x, \xi}) \subset \spn \{\Lambda_{\mathcal A}\} \cap \partial \Bbb_E \quad \text{for $(\lambda^s \otimes \rho_x)$-almost every $(x,\xi) \in {\Omega} \times Z$}\,.
		\]
	\end{lemma}
	\begin{proof}
		Recall that if $(\mu_\epsilon)_\epsilon$ generates the two-scale* Young measure $\bm \nu$, then the same sequence generates the generalized Young measure $\bm\upsilon = (\upsilon,\lambda,\upsilon^\infty)$ where $\upsilon^\infty$ is the weak-$*$ $\lambda$-measurable map $x \mapsto \upsilon_x^\infty \in \mathrm{Prob}(\partial \Bbb_E)$ and each probability measure $\upsilon_x$ is defined by duality as
		\begin{equation}\label{eq:com3}
		\dprb{h^\infty,\upsilon_x^\infty} = \int_Z \bigg(\int_{\partial \Bbb_E} h^\infty(\hat z) \;\mathrm{d} \nu^\infty_{x, \xi}(\hat z)\bigg) \;\mathrm{d} \rho_x(\xi) \quad \text{for all $h^\infty \in \E^0(E)$}\,.
		\end{equation}
		Moreover, if $(\mu_\epsilon)_\epsilon$ is originally an asymptotically $\mathcal A$-free sequence, then by definition $\bm{\upsilon}$ is an $\mathcal A$-free generalized Young measure (see~\cite{arroyo-rabasa2017lower-semiconti} for the corresponding definition). The localization principle~\cite[Proposition~2.25]{arroyo-rabasa2017lower-semiconti} (for generalized Young measures at singular points) yields the existence of a set $S \subset \Omega$ with full $\lambda^s$-measure and satisfying the following property: at each $x_0 \in S$ there exists an $\mathcal A$-free sequence $(u_\epsilon)_\epsilon \subset \mathrm{L}^1(Q;E)$ ---depending on the point $x_0$--- such that
		\begin{equation}\label{eq:com1}
		u_\epsilon \overset{\Ybf}{\to}  \bm{\upsilon(x_0)} = (\delta_0,\lambda_{x_0},\upsilon_{x_0}) \in \Y^1(Q;E), \quad \tilde \lambda_{x_0}(\partial Q) = 0,
		\end{equation}
		and $\tilde \lambda_{x_0} \in \Tan_1(\lambda,x_0)$ is a probability measure (hence a non-zero positive measure). Moreover, by Theorem~\cite[Theorem~1.1]{de-philippis2016on-the-structur}, we may further assume
		${\dprb{\id_{\partial \Bbb_E},\upsilon_{x_0}}} = {|\dprb{\id_{\partial \Bbb_E},\upsilon_{x_0}}|} \cdot {\;\mathrm{d} [\upsilon]}/{\;\mathrm{d} |[\upsilon]|^s}(x_0) \in \Lambda_\mathcal A$.
		In particular
		\begin{equation}\label{eq:com2}
		[\bm{\upsilon(x_0)}] \in \mathcal M(Q;\spn\{\Lambda_\mathcal A\}) \text{\quad for every $x_0 \in S$}\,.
		\end{equation}

		By properties~\eqref{eq:com1}-\eqref{eq:com2} we may apply~\cite[Lemma~3.2]{arroyo-rabasa2017lower-semiconti} to each generalized young measure $\bm{\upsilon(x_0)}$. Using~\eqref{eq:com2} once more we deduce
		\[
		\text{$\supp\,(\upsilon_{x_0}^\infty) \subset \spn\{\Lambda_\mathcal A\} \cap \partial \Bbb_E$ for every $x_0 \in S$}\,.
		\]
		Hence, by~\eqref{eq:com3}, $\supp\,(\nu^\infty_{x,\xi}) \subset  \spn\{\Lambda_\mathcal A\} \cap \partial \Bbb_E$ for $(\lambda^s \otimes \rho_x)$-almost every $(x,\xi) \in Q \times Z$. This finishes the proof.
	\end{proof}

	%

	%

	%
	%
	%
	%

	%

	%

	%

	%
	%
	%

	%
	%

	%

	%

	In the introduction we have defined the $\mathcal A$-free homogeneous envelope for integrands $\mathrm{C}(Z \times E)$ which are convex in their second argument. This definition is nothing else than a simplified representation of the (general) definition of $\mathcal A$-free homogeneous envelope defined for arbitrary integrands:

	\begin{definition}
		Let $h : Z \times E \to \mathbb{R}$ be a continuous integrand. The $\mathcal A$-free homogeneous envelope of $h$ is the integrand
		\begin{equation}\label{eq:hom}
		\begin{split}
		h_{\A-\mathrm{hom}}(z) & \coloneqq \inf\setBB{\,\thetant-_{Q_R} h(y,z + w(y)) \;\mathrm{d} y }{ R \in \mathbb{N}, \\
			& \qquad  \quad \qquad  w \in \mathrm{C}_\mathrm{per}^\infty(Q_R;E) \cap \ker \mathcal A, \int_{Q_R} w \;\mathrm{d} y = 0}\,, \quad z \in E\,.
		\end{split}
		\end{equation}
	\end{definition}
	The following relation about the commutability of the recession operation and the homogenization of an integrand holds.

	\begin{proposition}[recession regularization vs. homogenization] \label{prop:exinftyhom}

		Let $h \in \mathrm{C}(Z \times \mathbb{R}^N)$ be an integrand with linear-growth at infinity. Then
		\[
		(h^\#)_{{\A-\mathrm{hom}}} \ge (h_{{\A-\mathrm{hom}}})^\#.
		\]
		\begin{proof}

			Fix a vector $z \in E$. Let also $R \in \mathbb{N}$ and $w \in \mathrm{C}_\mathrm{per}^\infty(Q_R;E)$ as in~\eqref{eq:hom}. Recall that $h(\xi,t \,\sbullet\,)/t \le M(1 + |\,\sbullet\,|) \in \mathrm{L}^1_\mathrm{loc}(\mathbb{R}^d)$, hence we may use Fatou's lemma and the definition of $h^\#$ to obtain
			\begin{equation}\label{eq:inequality}
			\begin{split}
			\,\thetant-_{Q_R} h^\#(\xi,A + w(y))\;\mathrm{d} y  & =   \,\thetant-_{Q_R} \limsup_{t \to \infty} \frac{h(y,tz + tw(y))}{t} \;\mathrm{d} y \\
			& \geq \limsup_{t \to \infty} \,\thetant-_{Q_R} \frac{h(y,tz + tw(y))}{t} \;\mathrm{d} y \\
			& \geq  \limsup_{t \to \infty} \frac{h_{{\A-\mathrm{hom}}}(tz)}{t} = (h_{{\A-\mathrm{hom}}})^\#(z).
			\end{split}
			\end{equation}
			In passing to the first inequality we have used the linearity of both $\mathcal A$ and the mean value operation to ensure that $tw \in \mathrm{C}^\infty(Q_R;E) \cap \ker \mathcal A$, $\int_{Q_R} tw = 0$.
			Taking the infimum over such $w$'s first, and subsequently over all $R \in \mathbb{N}$ in~\eqref{eq:inequality} gives
			\[
			(h^\#)_{{\A-\mathrm{hom}}}(z) \ge (h_{{\A-\mathrm{hom}}})^\#(z),
			\]
			as desired.
		\end{proof}
	\end{proposition}

	\section{Convex homogenization}\label{sec:6}

	%

	\subsection{Jensen-type inequalities}
	This section is devoted to the study of (Jensen) integral inequalities satisfied by the {$\mathcal A$}-homogeneous envelope of convex integrands with respect to arbitrary $\mathcal A$-free two-scale* Young measures. The plan is to establish Jensen type inequalities at the first- and second-scale; naturally involving the first and second barycenters. Bridging these two inequalities into an \emph{homogenized}  Jensen inequality is, in turn, the key argument towards the proof of Theorem~\ref{thm:hom}. We close this section with an open problem and a discussion about non-convex integrands.

	\subsubsection{First-scale Jensen's inequality}

	\begin{proposition}\label{lem:J1}
		Let $\bm{\nu}= (\nu,\lambda,\rho,{\nu}^\infty) \in \Y^2(\Omega;E)$ be an $\mathcal A$-free two-scale* Young measure and let $h \in \mathrm{C}(Z \times E) \cap \mathbf{R}^2(\Omega;E)$ be an integrand that is convex in its second argument, that is, $h(\xi,\,\sbullet\,)$ is convex for all $\xi \in Z$. Then,
		\begin{enumerate}[(i)]
			\item at every regular point $x \in \mathrm{Reg}_{\bm{\nu}}(\Omega)$ it holds
			\begin{equation}\label{eq:J1r}
			\begin{split}
			h_{*\mathcal A}\bigg(\frac{\;\mathrm{d} \bpr{\bm{\nu}}}{\;\mathrm{d} \mathscr L^d}(x)\bigg) \leq \int_Z & h\bigg(\xi,\frac{\;\mathrm{d} \bbpr{\bm{\nu}}_x}{\;\mathrm{d} \mathscr L^d}(\xi)\bigg) \;\mathrm{d} \xi  \\
			& + \int_Z h^\infty\bigg(\xi,\frac{\;\mathrm{d} \bbpr{\bm{\nu}}^s}{\;\mathrm{d} \rho_x}(\xi)\bigg) \;\mathrm{d} \rho_x (\xi)\,,
			\end{split}
			\end{equation}
			\item and, at every singular point $x \in \mathrm{Sing}_{\bm{\nu}}(\cl\Omega)$,
			\begin{equation}\label{eq:J1s}
			(h^\infty)_{*\mathcal A}\bigg(\frac{\;\mathrm{d} \bpr{\bm{\nu}}_x}{\;\mathrm{d} \lambda^s}(x)\bigg) \le \int_Z h^\infty\bigg(\xi,\frac{\;\mathrm{d} \bbpr{\bm{\nu}}_x}{\;\mathrm{d} \rho_x}(\xi)\bigg) \;\mathrm{d} \rho_x(\xi)\,.
			\end{equation}
		\end{enumerate}
	\end{proposition}
	\begin{proof}
		Let  $\phi$ be a non-negative mollifier on $Z$ with $\int_Z \phi = 1$. Set $\phi_\delta \coloneqq \delta^{-d}\phi(x/\,\sbullet\,)$ so that $\phi_\delta \mathscr L^d_Z$ is a probability measure {on $Z$}.
		We define, for fixed $\delta > 0$ and $x \in \Omega$, the mollified second-order barycenter
		\[
		v_\delta(\xi) \coloneqq (\phi_\delta \ast \bbpr{\bm{\nu}}_x)(\xi)\,.
		\]
		Recall from  \cite[Remark 2.2]{arroyo-rabasa2017lower-semiconti} that convergence of mollified measures is strengthened to area-convergence, that is, $v_\delta$ area-converges to $\bbpr{\nu}_x$ on $Z$ (as $\delta \downarrow 0$). Hence, by \cite[Theorem~5]{kristensen2010characterizatio} we obtain
		\begin{equation}\label{eq:join1}
		\int_Z h(\xi,v_\delta(\xi)) \;\mathrm{d} \xi \to \int_Z h\bigg(\xi,\frac{\;\mathrm{d} \bbpr{\bm{\nu}}_x}{\;\mathrm{d} \mathscr L^d}(\xi)\bigg)
		+ \int_Z h^\infty\bigg(\xi,\frac{\;\mathrm{d} \bbpr{\bm{\nu}}^s_x}{\;\mathrm{d} |\bbpr{\bm{\nu}}_x^s|}(\xi)\bigg) \;\mathrm{d} |\bbpr{\bm{\nu}}_x^s|(\xi)\,.
		\end{equation}
		On the other hand, by the properties of mollifiers and the differential rigidity of the second-scale barycenter proved in Corollary~\ref{cor:Afree}, it holds that every test function $\tilde v_\delta \coloneqq v_\delta - \bbpr{\bm{\nu}}_x(Z)$ is $\mathcal A_\xi$-free and has zero mean value. This property of the $\tilde v_\delta$'s enables us to use the definition of the $\mathcal A$-homogeneous envelope (of a convex integrand) which yields
		\begin{equation}\label{eq:join2}
		\int_Z h(\xi,v_\delta(\xi)) \;\mathrm{d} \xi = \int_Z h(\xi,\tilde v_\delta(\xi) + \bbpr{\bm{\nu}}_x(Z)) \;\mathrm{d} \xi \ge h_{*\mathcal A}(\bbpr{\bm{\nu}}_x(Z))\,.
		\end{equation}
		From \eqref{eq:join1} and~\eqref{eq:join2} we conclude
		\begin{equation*}
		\int_Z h\bigg(\xi,\frac{\;\mathrm{d} \bbpr{\bm{\nu}}_x}{\;\mathrm{d} \mathscr L^d}(\xi)\bigg)
		+ \int_Z h^\infty\bigg(\xi,\frac{\;\mathrm{d} \bbpr{\bm{\nu}}^s_x}{\;\mathrm{d} |\bbpr{\bm{\nu}}_x^s|}(\xi)\bigg) \;\mathrm{d} |\bbpr{\bm{\nu}}_x^s|(\xi) \ge h_{*\mathcal A}(\bbpr{\bm{\nu}}_x(Z))\,.
		\end{equation*}
		Thus, taking into account that $|\bbpr{\bm{\nu}}_x^s| = |\dpr{\id_{\partial \Bbb_E},\nu^\infty_{x,\xi}}| \rho_x^s$ and $h^\infty(\xi,\dpr{\id_{\partial \Bbb_E},\nu^\infty_{x,\xi}}) = h^\infty(\xi,0) = 0$ for $\rho_x^*$, where $\rho_x^*$ is the singular part of $\rho_x^s$ with respect to $|\bbpr{\bm{\nu}}_x^s|$, we get (using the 1-homogeneity of $h^\infty$) the refined estimate
		\begin{equation}\label{eq:join3}
		\int_Z h\bigg(\xi,\frac{\;\mathrm{d} \bbpr{\bm{\nu}}_x}{\;\mathrm{d} \mathscr L^d}(\xi)\bigg)
		+ \int_Z h^\infty\bigg(\xi,\frac{\;\mathrm{d} \bbpr{\bm{\nu}}^s_x}{\;\mathrm{d} \rho_\nu^s}(\xi)\bigg) \;\mathrm{d} \rho_x^s(\xi) \ge h_{*\mathcal A}(\bbpr{\bm{\nu}}_x(Z))\,.
		\end{equation}
		In particular, if $x \in \mathrm{Reg}_\nu(\cl\Omega)$, we apply the inequality~\eqref{eq:join3} to deduce
		\begin{align*}
		\int_Z h\bigg(\xi,\frac{\;\mathrm{d} \bbpr{\bm{\nu}}_x}{\;\mathrm{d} \mathscr L^d}(\xi)\bigg)
		+ \int_Z h^\infty\bigg(\xi,\frac{\;\mathrm{d} \bbpr{\bm{\nu}}^s_x}{\;\mathrm{d} \rho_x^s}(\xi)\bigg) \;\mathrm{d} \rho_x^s(\xi) & \stackrel{\eqref{eq:density2}}{\ge} h_{*\mathcal A}\bigg(\frac{\;\mathrm{d} \bpr{\bm{\nu}}}{\;\mathrm{d} \mathscr L^d}(x)\bigg).
		\end{align*}
		This proves~\eqref{eq:J1r}. If on the other hand $x \in \mathrm{Sing}_\nu(\cl \Omega)$, we use~\eqref{eq:join3} (with $h^\infty$ in place of $h$) to get
		\begin{align*}
		\int_Z h^\infty\bigg(\xi,&\frac{\;\mathrm{d} \rho_x}{\;\mathrm{d} \mathscr L^d}(\xi)\dpr{\id_{E},\nu_{x,\xi}^\infty} \bigg) \;\mathrm{d} \xi +
		\int_Z h^\infty\bigg(\xi,\dpr{\id_{\partial \Bbb_E},\nu_{x,\xi}^\infty} \bigg)\;\mathrm{d} \rho_x^s (\xi) \\[5pt]
		& \qquad = \int_Z h^\infty\bigg(\xi,\frac{\;\mathrm{d} \bbpr{\bm{\nu}}_x}{\;\mathrm{d} \rho_x}(\xi)\bigg) \;\mathrm{d} \rho_x (\xi) \stackrel{\eqref{eq:density2}}{\ge} (h^\infty)_{*\mathcal A}\bigg(x,\frac{\;\mathrm{d} \bpr{\bm{\nu}}}{\;\mathrm{d} \lambda^s}(x)\bigg).
		\end{align*}

		This proves~\eqref{eq:J1s} and the proof is completed.
	\end{proof}

	\subsubsection{Second-scale Jensen inequalities}

	\begin{proposition}[at regular points]\label{lem:J2r} Let $\bm{\nu}= (\nu,\lambda,\rho,{\nu}^\infty) \in \Y^2(\Omega;E)$ be an $\mathcal A$-free two-scale* Young measure and let $h \in \mathrm{C}(Z \times E) \cap \mathbf{R}^2(\Omega;\E)$ be an integrand that is convex in its second argument, that is, $h(\xi,\,\sbullet\,)$ is convex for all $\xi \in Z$. Then, for all regular points $x \in \mathrm{Reg}_{\bm{\nu}}(\cl\Omega)$ the following inequalities hold:

		\begin{enumerate}[(i)]
			\item at $\mathscr L^d_Z$-almost every $\xi$ in $Z$,
			\begin{equation}\label{eq:J2r}
			h\bigg(\xi,\frac{\;\mathrm{d} \bbpr{\bm{\nu}}_x}{\;\mathrm{d} \mathscr L^d}(\xi)\bigg) \le \frac{\;\mathrm{d} \bbpr{h,\bm{\nu}}_x}{\;\mathrm{d} \mathscr L^d}(\xi),
			\end{equation}
			\item and, at $\rho_x^s$-almost every $\xi$ in $Z$,
			\begin{equation}\label{eq:J2r'}
			h^\infty\bigg(\xi,\frac{\;\mathrm{d} \bbpr{\bm{\nu}}_x}{\;\mathrm{d} \rho_x^s}(\xi)\bigg) \le \frac{\;\mathrm{d} \bbpr{h^\infty,\bm{\nu}}_x}{\;\mathrm{d} \rho_x^s}(\xi).
			\end{equation}
		\end{enumerate}

	\end{proposition}
	\begin{proof}Let $x \in \mathrm{Reg}_\nu(\Omega)$ and $\xi \in Z$ so that ${\;\mathrm{d} \bbpr{\bm{\nu}}_{x}}/{\;\mathrm{d} \mathscr L^d}(\xi)$ exists (note that this property is satisfied at $\mathscr L^d_Z$-almost every $\xi$ in $Z$). Observe  that since $h$ is $\xi$-uniformly Lipschitz in its second argument, ${\;\mathrm{d} \bbpr{h,\bm{\nu}}_{x}}/{\;\mathrm{d} \mathscr L^d}(\xi)$ also exists at such $\xi$'s. We have
		\[
		\frac{\;\mathrm{d} \bbpr{\bm{\nu}}_{x}}{\;\mathrm{d} \mathscr L^d}(\xi) =\dpr{\id_{E},\nu_{x,\xi}} + \frac{\;\mathrm{d} \lambda}{\;\mathrm{d} \mathscr L^d}(x)\cdot \frac{\;\mathrm{d}\rho_x}{\mathscr L^d}(\xi)\cdot \dpr{\id_{\partial \Bbb_E},\nu_{x,\xi}^\infty}.
		\]

		By the classical Jensen's inequality and the positive $1$-homogeneous character of $h^\infty$ we further obtain
		\begin{align*}
		\frac{\;\mathrm{d} \bbpr{h,\bm{\nu}}_x}{\;\mathrm{d} \mathscr L^d}(\xi) & = \dpr{h(\xi,\,\sbullet\,),\nu_{x,\xi}}  + \frac{\;\mathrm{d} \lambda}{\;\mathrm{d} \mathscr L^d}(x)\cdot \frac{\;\mathrm{d}\rho_x}{\mathscr L^d}(\xi)\cdot \dpr{h^\infty(\xi,\,\sbullet\,),\nu_{x,\xi}^\infty}  \nonumber \\
		& \ge h(\xi,\dpr{\id_{E},\nu_{x,\xi}}) + h^\infty\bigg(\xi,\frac{\;\mathrm{d} \lambda}{\;\mathrm{d} \mathscr L^d}(x)\cdot \frac{\;\mathrm{d}\rho_x}{\mathscr L^d}(\xi)\cdot \dpr{\id_{\partial\Bbb_E},\nu_{x,\xi}^\infty}\bigg).
		\end{align*}
		Now, let us recall the following sub-additive property satisfied by convex functions and their recession functions: every convex function $g : E \to \mathbb{R}$
		$g(z_1 + z_2) \le g(z_1) + g^\infty(z_2)$ for all $z_1,z_2 \in E$. It follows directly from this observation and the inequality above that
		\[
		\frac{\;\mathrm{d} \bbpr{h,\bm{\nu}}_x}{\;\mathrm{d} \mathscr L^d}(\xi) \ge  h\bigg(x,\xi,\frac{\;\mathrm{d} \bbpr{\bm{\nu}}_x}{\;\mathrm{d} \mathscr L^d}(\xi)\bigg)\,.
		\]

		This proves the first inequality in~\eqref{eq:J2r}.

		For the proof of the second inequality we let $\xi \in Z$ be such that ${\;\mathrm{d} \bbpr{\bm{\nu}}_x}/{\;\mathrm{d} \rho_x^s}(\xi)$ exists, hence also ${\;\mathrm{d} \bbpr{h,\bm{\nu}}_x}/{\;\mathrm{d} \rho_x^s}(\xi)$ exists ($h^\infty$ is $\xi$-uniformly Lipschitz continuous in its second argument). In fact, these two measure-theoretic derivatives are given by the parings $\dpr{\id_{\partial\Bbb_E}, \nu^\infty_{x,\xi}}$ and $\dpr{h(\xi,\,\sbullet\,), \nu^\infty_{x,\xi}}$ respectively. Once more, a simple application of Jensen's classical inequality gives
		\begin{align*}
		\frac{\;\mathrm{d} \bbpr{h,\bm{\nu}}_x}{\;\mathrm{d} \rho_x^s}(\xi) & = \dpr{h^\infty(\xi,\,\sbullet\,),\nu^\infty_{x,\xi}}  \\
		& \geq h^\infty(\xi,\dpr{\id_{\partial\Bbb_E},\nu_{x,\xi}^\infty})
		\end{align*}
		Since ${\;\mathrm{d} \bbpr{\bm{\nu}}_x}/{\;\mathrm{d} \rho_x^s}(\xi)$ and $\dpr{h(\xi,\,\sbullet\,), \nu^\infty_{x,\xi}}$ exist $\rho_x^s$-almost everywhere in $Z$, this proves~\eqref{eq:J2r'}.\end{proof}

	%

	%

	\begin{proposition}[at singular points] Let $\bm{\nu}= (\nu,\lambda,\rho,{\nu}^\infty) \in \Y^2(\Omega;E)$ be an $\mathcal A$-free two-scale* Young measure and let $g \in \mathrm{C}(Z \times \mathbb{R}^N)$ be continuous integrand. Further, assume that $g(\xi,\,\sbullet\,)$ is convex and positively 1-homogeneous for all $\xi \in Z$. Then, at every $x \in \mathrm{Sing}_{\bm\nu}(\Omega)$ it holds that
		\begin{equation}\label{eq:J2s}
		g\bigg(\xi,\frac{\;\mathrm{d} \bbpr{{\bm\nu}}_x}{\;\mathrm{d} \rho_x}(\xi)\bigg) \le \frac{\;\mathrm{d} \bbpr{g,{\bm\nu}}_x}{\;\mathrm{d} \rho_x}(\xi) \quad \text{for $\rho_x$-almost every $\xi \in Z$}\,.
		\end{equation}
	\end{proposition}
	\begin{proof}
		The proof can be reproduced by following the exact same ideas in the proof of the second part of the proof of Proposition~\ref{lem:J2r}.
	\end{proof}

	We are now in place to prove the \emph{homogeneous} version of Jensen's inequality that involves the $\mathcal A$-$\hom$ envelope (of convex integrands).

	\begin{theorem}[homogenized Jensen's inequality]\label{thm:J}
		Let $\bm{\nu}= (\nu,\lambda,\rho,{\nu}^\infty) \in \Y^2(\Omega;E)$ be an $\mathcal A$-free two-scale* Young measure and let $h \in \mathrm{C}(Z \times E) \cap \mathbf{R}^2(E)$ be an integrand that is convex in its second argument, that is, $h(\xi,\,\sbullet\,)$ is convex for all $\xi \in Z$. Then,
		\begin{enumerate}
			\item at every regular point $x \in \mathrm{Reg}_{\bm\nu}(\Omega)$ it holds that
			\begin{equation}\label{eq:J1rr}
			\begin{split}
			h_{*\mathcal{A}}\bigg(\frac{\;\mathrm{d} \bpr{{\bm\nu}}}{\;\mathrm{d} \mathscr L^d}(x)\bigg) \leq \frac{\;\mathrm{d} \bbpr{h,{\bm\nu}}_x(Z)}{\;\mathrm{d} \mathscr L^d}(x)\,,
			\end{split}
			\end{equation}
			\item and, at every singular point  $x \in \mathrm{Sing}_{\bm\nu}(\Omega)$,
			\begin{equation}\label{eq:J1ss}
			(h^\infty)_{*\mathcal{A}}\bigg(\frac{\;\mathrm{d} \bpr{{\bm\nu}}}{\;\mathrm{d} \lambda^s}(x)\bigg) \le  \frac{\;\mathrm{d} \bbpr{h,{\bm\nu}}_x(Z)}{\;\mathrm{d} \lambda^s}(x)\,.
			\end{equation}
		\end{enumerate}
	\end{theorem}
	\begin{proof}
		The proof is a direct consequence of Propositions~\ref{lem:J2r},~\ref{lem:J2r}, and~\ref{lem:J1}.
	\end{proof}

	\subsubsection{Comments on the non-convex case}

	In general, even at singular points $x \in\mathrm{Sing}_{\bm{\nu}}\cap \cl \Omega$, one \emph{cannot} expect the second-scale differential inclusion
	\begin{equation}\label{eq:fail}
	\frac{\;\mathrm{d} \bbpr{\bm{\nu}}_{x_0}}{\;\mathrm{d} \mathscr L^d_Z}(\xi) \in \Lambda_{\mathcal A} \quad \text{for $\mathscr L^d_Z$-almost every $\xi \in Z$}
	\end{equation}
	to hold.
	Instead, we believe the following weaker statement holds under mild non-degeneracy assumptions on $\mathcal A$ (for instance if $\mathcal A$ satisfies Murat's constant rank condition).
	\begin{conjecture}\label{con:sera} Let $\bm{\nu} \in \Y^2(\Omega;\mathbb{R}^N)$ be an $\mathcal A$-free two-scale* Young measure and let $h \in \mathrm{C}(Z \times \mathbb{R}^N)$ be a positively 1-homogeneous for all $\xi \in Z$.
		Then, for all $x_0 \in \mathrm{Sing}_{\bm{\nu}}$ it holds that
		\begin{equation}\label{eq:conjecture}
		\bbpr{h,\bm\nu}_{x_0}(Z)
		\ge h_{{\A-\mathrm{hom}}}(\bbpr{\bm\nu}_{x_0}(Z))
		\end{equation}
	\end{conjecture}
	This is a \emph{powerful} inclusion, and, in fact, it is the key inequality towards the characterization of the homogenization of \emph{non-convex}  integrals. The proof of this conjecture is relatively simple provided that $\mathcal A$ is a first-order operator (see~\cite{matias2015homogenization-}).

	For operators of general order, the veracity of this conjecture seems to be linked to the \emph{structural}  properties of tangent $\mathcal A$-free measures.

	\subsubsection{The conjecture in $\mathrm{BD}$}
	A natural candidate to study the elastic perfectly plastic behavior of materials is the space of
	functions whose linearized strains are measures. Formally, these deformations belong to the space of \emph{functions of  bounded deformation} which is defined as
	\[
	\mathrm{BD}(\Omega) \coloneqq \setb{u \in \mathrm{L}^1(\Omega;\mathbb{R}^d)}{Eu \coloneqq Du + Du^T \in \mathcal M(\Omega;\mathbf{M}_{d \times d}^\mathrm{sym})}.
	\]
	In the $\mathcal A$-free context, a measure $\mu \in \mathcal M(\Omega;\mathbf{M}_{d \times d}^\mathrm{sym})$ is locally a symmetric gradient if and only if it satisfies (see~\cite[Example~3.10(e)]{fonseca1999mathcal-a-quasi}) the second order PDE-constraint
	\[
	\curl \curl \mu \coloneqq \bigg(\sum_{i=1}^d \partial_{ik} \mu_{i}^j+\partial_{ij} \mu_{i}^k-\partial_{jk} \mu_{i}^i-\partial_{ii} \mu_{j}^k\bigg)_{jk} \qquad j,k=1,\ldots,d\,.
	\]
	In a forthcoming~\cite{arroyo-rabasa2018rigidity-of-tan} paper we shall give a positive answer to Conjecture~\ref{con:sera} in the case $\bm{\nu}$ is $\mathbf{BDY}^2(\Omega)$ two-scale* Young measure, that is, when
	\[
	Eu_j \overset{\Ybf^2}{\to} \bm{\nu} \qquad (u_j)_{j \in \mathbb{N}} \subset \mathrm{BD}(\Omega).
	\]

	\subsection{Proof of the Theorem~\ref{thm:hom}} By the definition of $\Gamma$-limit it is enough to verify the conclusion of the theorem for infinitesimal sequences. Let $\epsilon_j\searrow 0$ and $\mu_n \overset{*}\rightharpoonup \mu$ in $\mathcal M(\cl \Omega;E) \cap \ker \mathcal A$.
	Observe that up to extracting a subsequence we may assume without loss of generality that
	\[
	\liminf_{j \to \infty} I^{\epsilon_j}      (\mu_j) = \lim_{j \to \infty} I^{\epsilon_j}(\mu_j)
	\]
	and
	\[
	\text{and} \quad \mu_n \overset{\Ybf^2}{\to} \bm{\nu} \quad \text{for some $\bm{\nu}= (\nu,\lambda,\rho,\nu^\infty)\in \Y^2(\Omega;E)$}\,.
	\]

	The Young-measure representation and the Radon--Nykodym differentiation give
	\begin{align*}
	\liminf_{j \to \infty} I^{\epsilon_j}(\mu_j) & = \ddprb{f,\bm{\nu}} \\
	&  = \int_{\cl \Omega} \bbpr{f,\bm{\nu}}_x(Z) \;\mathrm{d}(\mathscr L^d + \lambda^s)(x) \\
	& =  \int_{\Omega} \frac{\bbpr{f,\bm{\nu}}_x(Z)}{\;\mathrm{d} \mathscr L^d}(x) \;\mathrm{d} x + \int_{\cl \Omega} \frac{\bbpr{f,\bm{\nu}}_x(Z)}{\;\mathrm{d} \lambda^s}(x) \;\mathrm{d} \lambda^s(x)\,.

	\end{align*}
	An $x$-point-wise application of the Jensen inequalities in Theorem~\ref{thm:J} with the family of integrands $\{h(\xi,z)\}_x \coloneqq \{f(x,\xi,z)\}_x$ on the right-hand side above yields (notice that it suffices to perform this at all $x \in \cl \Omega$ where the conclusion of Theorem~\ref{thm:J} holds)
	\begin{align*}
	\liminf_{j \to \infty} I^{\epsilon_j}(\mu_j) & \ge \int_{\Omega} f_{*\mathcal{A}}\bigg(x,\frac{\;\mathrm{d} \mu}{\;\mathrm{d} \mathscr L^d}(x)\bigg) \;\mathrm{d} x + \int_{\cl \Omega} (f^\infty)_{*\mathcal{A}}\bigg(x,\frac{\;\mathrm{d} \bpr{\nu}_x}{\;\mathrm{d} (\lambda^\nu)^s}(x)\bigg) \;\mathrm{d} \lambda^s(x) \\
	& \ge \int_{\Omega}f_{*\mathcal{A}}\bigg(x,\frac{\;\mathrm{d} \mu}{\;\mathrm{d} \mathscr L^d}(x)\bigg) \;\mathrm{d} x + \int_{\cl \Omega} (f^\infty)_{*\mathcal{A}}\bigg(x,\frac{\;\mathrm{d} \mu}{\;\mathrm{d} |\mu^s|}(x)\bigg) \;\mathrm{d} |\mu^s|(x)\,.
	\end{align*}
	Here, to reach the second inequality we have used that $|\dpr{\id_{\mathbb{R}^N},\nu_x}|(\lambda^\nu)^s \equiv |\mu^s|$ and that $\frac{\;\mathrm{d} |\mu^s|}{\;\mathrm{d} \lambda^s}(x) = 0$ for $\lambda^*$-a.e. $x \in \cl \Omega$ where $\lambda^*$ is the singular part of $\lambda^s$ with respect to $|\mu^s|$.

	Lastly, it follows from Proposition~\ref{prop:exinftyhom} that $(f^\infty)_{*\mathcal{A}} \ge (f_{*\mathcal{A}})^\infty$ whence we conclude
	\[
	\liminf_{j \to \infty} I^{\epsilon_j}(\mu_j) \ge\int_{\Omega} f_{*\mathcal{A}}\bigg(x,\frac{\;\mathrm{d} \mu}{\;\mathrm{d} \mathscr L^d}(x)\bigg) \;\mathrm{d} x + \int_{\cl \Omega} (f_{*\mathcal{A}})^\infty\bigg(x,\frac{\;\mathrm{d} \mu}{\;\mathrm{d} |\mu^s|}(x)\bigg) \;\mathrm{d} |\mu^s|(x)\,.
	\]
	This finishes the proof.\qed

redpajama_set_name: RedPajamaArXiv