\section{Introduction and motivation}

This paper is motivated by problems arising in multi-period valuation in insurance, but the applicability of the results extends beyond the insurance setting.
Aggregate insurance liability cashflows typically correspond to payments several years into the future. However, insurance regulation requires that capital requirements are calculated for a one-year horizon, considering cashflows during the year and end-of-year liability values.
This implies that liability values must be computed recursively, backwards in time, starting from the year of the most distant liability payments.
Solving such backward recursions leads to major computational challenges. Consequently, simple and explicit (standard) formulas with questionable conceptual soundness have become the preferred alternative.
With this paper we say that there are alternatives to commonly encountered formulas that give rise to explicit formulas that are easy to use in practice and that retain both economic interpretability and the conceptual soundness of the original principles for insurance valuation. Mathematically, the explicit formulas follow from combining the original principles for insurance valuation with widely applicable large portfolio asymptotics.

We begin by presenting the motivating insurance problem.
Let $(C^{n_0},\mathbb{F}^{n_0})$ be a stochastic processes in discrete time with $(C^{n_0}_t)_{t=1}^{T}$ denoting the discounted payments at times $t$ due to insurance claims, and $\mathbb{F}^{n_0}=(\mathcal{F}^{n_0}_t)_{t=0}^T$, with $\mathcal{F}^{n_0}_0=\{\Omega,\emptyset\}$, denoting a filtration to which $C^{n_0}$ is adapted. $C^{n_0}$ describes discounted claim payments over time in a so-called runoff situation, when no new contracts are written. The number $n_0$ is intended as a measure of volume or exposure, for instance the number of contracts that may generate future claims payments.
We want to assign a value to this liability cashflow. In order to do this we consider a sequence $(C^{n},\mathbb{F}^n)_{n\geq 1}$ of such stochastic processes, with $\Vert C^{n} \Vert \to \infty$ almost surely as $n\to\infty$, where $\Vert \cdot \Vert$ denotes the Euclidean norm on $\mathbb{R}^T$.
An example could be $C^{n}$ given by
\begin{align*}
C^{n}_t = \sum_{k=1}^{M_n} 1_{\{D_k=t\}}Z_k, \quad t=1,\dots,T,
\end{align*}
where $M_n$ denotes the total number of claims payments and $M_n\to\infty$ almost surely as $n\to\infty$,
$D_k$ denotes the time and $Z_k$ denotes the size of the $k$th claims payment. For this example, the filtration could be the one generated by the discounted cashflow $C^{n}$ or a larger filtration also including information about the number of payments $\sum_{k=1}^{M_n} 1_{\{D_k=t\}}$ at each point in time, or more.

We assume nonrandom sequences with terms $a_n\in (0,\infty)$ and $b_n\in \mathbb{R}^T$ such that there is convergence in distribution
\begin{align*}
X^n=a_n^{-1}(C^{n}-b_n)\stackrel{d}{\to} X \quad \text{as } n\to\infty.
\end{align*}
The value $V^n_0(C^{n})$ of the insurance liabilities in a multi-period model is the result of applying a suitable functional to the pair $(C^{n},\mathbb{F}^n)$. Natural valuation functionals satisfy the property $V^n_0(C^{n})=a_nV^n_0(X^n)+\sum_{s=1}^{T}b_{n,s}$. Convergence $X^n\stackrel{d}{\to} X$ suggests $V^n_0(X^n)\to V_0(X)$ and therefore the approximation
\begin{align}\label{eq:value_approximation}
	V^n_0(C^{n})\approx a_nV_0(X)+\sum_{s=1}^{T}b_{n,s}.
\end{align}
However, the flow of information over time is an essential ingredient in valuation and it is not true that convergence $X^n\stackrel{d}{\to} X$ implies the convergence $V^n_0(X^n)\to V_0(X)$ (even if $\mathbb{F}^n$ is taken to be the filtration generated by $C^{n}$). Consequently, one of the main objectives of the present paper is to determine sharp sufficient conditions for the convergence of $V^n_0(X^n)$ to $V_0(X)$ in terms of an appropriate mode of convergence of $(X^n,\mathbb{F}^n)$ to $(X,\mathbb{F})$, similar to so-called extended weak convergence introduced by Aldous in \cite{Aldous-81} and conditional weak convergence studied by Sweeting in \cite{Sweeting-89}.
The importance of the approximation \eqref{eq:value_approximation} is because $V^n_0(C^{n})$ is typically very difficult to compute numerically whereas $V_0(X)$ is easier to compute numerically and, more importantly, in the case of a Gaussian limit model, is given by an explicit expression in terms of conditional variances of components of $X$ (note that Gaussian vectors have the rare feature that conditional variances of one component given a subset of components are nonrandom).

The paper is organized as follows. Section \ref{sec:preliminaries} introduces notation and basic properties of conditional distributions and risk measures. Section \ref{sec:compv} presents the main contents of the paper and, following a general presentation of the mathematical setup, gives economic motivation of key quantities in Section \ref{sec:economicmotivation} and presents the main results in Section \ref{sec:mainresults}. All proofs together with auxiliary results are found in Section \ref{sec:proofs}.

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

$\mathbb{N}=\{1,2,\dots\}$, $\mathbb{R}$ denotes the real numbers and $\mathbb{R}_+=[0,\infty)$. Whenever relevant, for $T\in \mathbb{N}$, a vector $x\in \mathbb{R}^T$ is assumed to be column vector and its transpose $x^{{\mathrm{T}}}$ a row vector. For $x=(x_1,\dots,x_T)^{{\mathrm{T}}}\in \mathbb{R}^T$, ${\Vert x \Vert}^2=x^{{\mathrm{T}}}x$ and we let $x_{\leq t}=(x_1,\dots,x_t)^{{\mathrm{T}}}$ and $x_{> t}=(x_{t+1},\dots,x_T)^{{\mathrm{T}}}$.
For $d\in\mathbb{N}$ and a Borel set $A\in\mathcal{B}(\mathbb{R}^d)$, $\mathcal{P}(A)$ denotes the set of probability measures on $A$. For $\mu\in\mathcal{P}(A)$, $\operatorname{supp}(\mu)$ denotes its support.
We consider a probability space $(\Omega,\mathcal{F},\operatorname{P})$.
For a $\sigma$-algebra $\mathcal{F}_t\subset \mathcal{F}$, $L^0(\mathcal{F}_t,\operatorname{P})$ denotes the vector space of all real-valued $\mathcal{F}_t$-measurable random variables, and for $p\in (0,\infty)$, $L^p(\mathcal{F}_t,\operatorname{P})$ denotes the subset $\{Y\in L^0(\mathcal{F}_t,\operatorname{P}):\operatorname{E}[|Y|^p]<\infty\}$.

For an $\mathbb{R}^d$-valued random variable $Y$ we let $\mathcal{L}(Y)$ denote its distribution, i.e.~the induced probability measure $\operatorname{P}(Y\in \cdot)$ on $\mathbb{R}^d$. Given $\sigma$-algebras $\mathcal{F}_t\subset \mathcal{F}$ and a random variable $Y$, a regular conditional distribution $Q_{\mathcal{F}_t,Y}$ is a version of $\operatorname{P}(Y\in \cdot\mid \mathcal{F}_t)$ which forms a probability kernel from $(\Omega,\mathcal{F}_t)$ to $(\mathbb{R}^d,\mathcal{B}(\mathbb{R}^d))$.
If $\mathcal{F}_t=\sigma(Z)$ for an $\mathbb{R}^{d'}$-valued random variable $Z$, then $Q_{\mathcal{F}_t,Y}$ is an $\mathcal{F}_t$-measurable random measure $\omega\mapsto\kappa(Z(\omega),\cdot)$ on $\mathbb{R}^d$, where $\kappa$ is a probability kernel from $(\mathbb{R}^{d'},\mathcal{B}(\mathbb{R}^{d'}))$ to $(\mathbb{R}^d,\mathcal{B}(\mathbb{R}^d))$, see Theorem 6.3 in Kallenberg \cite{Kallenberg-02}. The notation $\mathcal{L}(Y\mid Z=z)$ means $\kappa(z,\cdot)$. In particular, $z\mapsto\mathcal{L}(Y\mid Z=z)$ is well defined on $\operatorname{supp}(\mathcal{L}(Z))$.

We define the conditional $p$-quantile of $Y\in L^0(\mathcal{F},\operatorname{P})$ given $\mathcal{F}_t$ as the random variable $F_{Y\mid \mathcal{F}_t}^{-1}(p)$ given by
\begin{align*}
\omega\mapsto \min\{m\in \mathbb{R}:Q_{\mathcal{F}_t,Y}(\omega,(-\infty,m])\geq p\}.
\end{align*}
If $Q_{\mathcal{F}_t,Y}(\omega,\cdot)=\kappa(Z(\omega),\cdot)$, then $F_{Y\mid Z=z}^{-1}(p)$ means $\min\{m\in \mathbb{R}:\kappa(z,(-\infty,m])\geq p\}$ and is well defined on $\operatorname{supp}(\mathcal{L}(Z))$.
Similar to ordinary quantiles, the conditional quantiles have the property $F_{aY+b\,\mid \mathcal{F}_t}^{-1}(p)=aF_{Y\mid \mathcal{F}_t}^{-1}(p)+b$ if $a\in\mathbb{R}_+$ and $b\in L^0(\mathcal{F}_t,\operatorname{P})$ (in particular if $b\in\mathbb{R}$).
Monetary conditional risk measures appear naturally for multi-period valuations, see e.g.~Chapter 11 in F\"ollmer and Schied \cite{foell2016}.
The risk measures Value-at-Risk and Average Value-at-Risk, conditional on $\mathcal{F}_t$ and for $u\in (0,1)$, are given by
\begin{align*}
\operatorname{V@R}_{u}(Y\mid \mathcal{F}_t)&=F_{-Y\mid \mathcal{F}_t}^{-1}(1-u), \\
\operatorname{AV@R}_{u}(Y\mid \mathcal{F}_t)&=\frac{1}{u}\int_0^{u}\operatorname{V@R}_{v}(Y\mid\mathcal{F}_t)dv.
\end{align*}
In \cite{foell2016}, $\operatorname{V@R}_{u}(Y\mid \mathcal{F}_t)$ is defined as $\operatorname{essinf}\{m_t\in L^0(\mathcal{F}_t,\operatorname{P}):\operatorname{P}(Y+m_t<0\mid\mathcal{F}_t)\leq u\}$ but the two definitions of $\operatorname{V@R}_{u}(Y\mid \mathcal{F}_t)$ are equal almost surely. For $p\geq 1$, $\operatorname{V@R}_{u}(Y\mid \mathcal{F}_t)\in L^p(\mathcal{F}_t,\operatorname{P})$ if $Y\in L^p(\mathcal{F},\operatorname{P})$, and similarly for $\operatorname{AV@R}_{u}(Y\mid \mathcal{F}_t)$.
Let $\mathcal{P}([0,1])'$ denote the subset of $\mathcal{P}([0,1])$ consisting of measures $\mu\in\mathcal{P}([0,1])$ that either have a bounded density with respect to Lebesgue measure or satisfy $\operatorname{supp}(\mu)\subset [a,b]$ for some $0<a<b<1$.
Both $\operatorname{V@R}_{u}(Y\mid \mathcal{F}_t)$ and $\operatorname{AV@R}_{u}(Y\mid \mathcal{F}_t)$ can be expressed as
\begin{align}\label{eq:quantriskmeas}
\rho(Y\mid \mathcal{F}_t)=\int_0^1 F_{-Y\mid \mathcal{F}_t}^{-1}(p)\mu(dp), \quad \mu\in \mathcal{P}([0,1])'.
\end{align}
$\operatorname{V@R}_{u}(Y\mid \mathcal{F}_t)$ corresponds to $\mu(dp)=\delta_{1-u}(dp)$ (a unit point mass at $1-u$) and $\operatorname{AV@R}_{u}(Y\mid \mathcal{F}_t)$ corresponds to $\mu(dp)=u^{-1}I\{p\in [1-u,1]\}(dp)$ (a bounded density).
Conditional risk measures of the form \eqref{eq:quantriskmeas} satisfy $\rho(aY+b\mid \mathcal{F}_t)=a\rho(Y\mid \mathcal{F}_t)-b$ if $a\in\mathbb{R}_+$ and $b$ is $\mathcal{F}_t$-measurable, called positive homogeneity and conditional cash additivity. Moreover, they are monotone: $\widetilde{Y}\geq Y$ implies $\rho(\widetilde{Y}\mid \mathcal{F}_t)\leq \rho(Y\mid \mathcal{F}_t)$. Throughout the paper (in)equalities between random variables should be interpreted in the almost sure sense.

\section{Convergence of multi-period valuations}\label{sec:compv}

Fix $T\in\mathbb{N}$ and let $X,X^1,X^2,\dots$ be random vectors in $\mathbb{R}^T$ and let $Y,Y^1,Y^2,\dots$ be random vectors in $(\mathbb{R}^d)^T$, $d\in\mathbb{N}$.
Suppose that $\mathcal{L}(X^n,Y^n)\stackrel{w}{\to} \mathcal{L}(X,Y)$ as $n\to\infty$, where $\mathcal{L}(X,Y)$ is Gaussian.
Set $\mathcal{F}_t=\sigma((X,Y)_{\leq t})$, $\mathcal{F}^n_t=\sigma((X^n,Y^n)_{\leq t})$. $(X,Y)$ and $(X^n,Y^n)$ are adapted discrete-time stochastic processes with respect to the filtrations $(\mathcal{F}_t)_{t=0}^T$ and $(\mathcal{F}^n_t)_{t=0}^T$, respectively, where $\mathcal{F}_0=\mathcal{F}^n_0=\{\Omega,\emptyset\}$.
For each $n$, $X^n=(X^n_t)_{t=1}^T$ corresponds to a discounted incremental cashflow in a multi-period model with $T$ periods and time points $0,1,\dots,T$.
For each $n$, $Y^n=(Y^n_t)_{t=1}^T$ corresponds to a stochastic process that provides additional information, additional to the information provided by the process $X^n$. Taking $Y^n$ to be nonrandom means that no such additional information is considered and that $(\mathcal{F}^n_t)_{t=0}^T$ is the natural filtration generated by $X^n$.
The discounted value at time $t$ of the cashflow occurring after time $t$ is denoted by $V^n_t(X^n)$. Since no cashflows occur after time $T$, $V^n_T(X^n)=0$. The process of values of the cashflow $X^n$ is $(V^n_t(X^n))_{t=0}^T$. Both $X^n_t$ and $V^n_t(X^n)$ are $\mathcal{F}^n_t$-measurable.
The processes $(V_t(X))_{t=0}^T$ and $(V^n_t(X^n))_{t=0}^T$ are defined backward recursively by
\begin{align}
& V_T(X)=0, \quad V_t(X)=\phi_t(X_{t+1}+V_{t+1}(X)), \quad t<T, \label{eq:VtVal}\\
& V^n_T(X^n)=0, \quad V^n_t(X^n)=\phi^n_t(X^n_{t+1}+V^{n}_{t+1}(X^n)), \label{eq:VntVal} \quad t<T,
\end{align}
where $\phi_t : L^1(\mathcal{F}_{T})\to L^1(\mathcal{F}_{t})$ and $\phi^n_t : L^1(\mathcal{F}^n_{T})\to L^1(\mathcal{F}^n_{t})$ are mappings that are made precise below. In particular, the mappings $\phi_t$ and $\phi^n_t$ will satisfy
\begin{align}\label{eq:phitab}
\phi_t(aY+\widetilde{Y})=a\phi_t(Y)+\widetilde{Y}, \quad
 \phi^n_t(aY+\widetilde{Y})=a\phi^n_t(Y)+\widetilde{Y} \quad \text{for } a\in\mathbb{R}_+, \widetilde{Y}\in L^1(\mathcal{F}_t).
\end{align}
These properties are called positive homogeneity (the effect of multiplication by a nonnegative scalar) and conditional cash additivity.
As a consequence, for $a\in\mathbb{R}_+, b\in \mathbb{R}^T$,
\begin{align}\label{eq:Vtab}
V_t(aX+b)=aV_t(X)+\sum_{s=t+1}^Tb_s, \quad
V^n_t(aX^n+b)=aV^n_t(X^n)+\sum_{s=t+1}^Tb_s.
\end{align}
Theorems \ref{thm:convergentvaluations} and \ref{thm:CondExpVar} presented below essentially say that if $\mathcal{L}(X^n,Y^n)\stackrel{w}{\to} \mathcal{L}(X,Y)$ as $n\to\infty$ and also in terms of conditional distributions (see e.g.~\eqref{eq:contconv}) together with a uniform integrability assumption (see e.g.~\eqref{eq:condxnui}), then $\lim_{n\to\infty}V^n_0(X^n)=V_0(X)$, where the values $V^n_0(X^n)$ and $V_0(X)$ are defined with respect to the filtrations generated by $(X^n,Y^n)$ and $(X,Y)$, respectively. The essential point is that $V^n_0(X^n)$ is typically impossible to compute analytically and only with significant difficulties numerically, whereas $V_0(X)$ has an explicit expression that is straightforward to compute in terms of the mean vector and covariance matrix of the Gaussian weak limit (see e.g.~\eqref{eq:V0expression}).
Therefore, the theorems we present enable the use of conceptually sound valuation techniques without the significant efforts otherwise needed to obtain numerical solutions to backward recursions.

\subsection{Economic motivation}\label{sec:economicmotivation}

Current regulatory frameworks for the insurance industry prescribe so-called cost-of-capital valuation of liability cashflows.
Such approaches to valuation consider capital requirements and the costs stemming from the financing of buffer capital. In the multi-period setting, with capital requirements determined one period at the time, the randomness of future capital requirements and associated costs drive the liability valuation.

Multiperiod cost-of-capital valuation is studied in
Salzmann and W\"uthrich \cite{SalzmannWuthrich-10},
M\"ohr \cite{mhr11},
Pelsser and Salahnejhad Ghalehjooghi \cite{Pelsser-SG-16}
and
Engsner et al. \cite{engsner2017}.
A common theme is that multiperiod valuations is constructed through backward induction of one-period valuations.
Pelsser and Salahnejhad Ghalehjooghi \cite{Pelsser-SG-16} study continuous-time limits of multiperiod valuations defined in terms of one-step valuations, similar to those considered below and including cost-of-capital valuation. In \cite{Pelsser-SG-16} the convergence takes place as the length of the time periods tends to zero and the number of time periods tends to infinity. Clearly, this is a different kind of convergence than that studied here.
Multiperiod valuations are so-called time-consistent by their construction through the backward induction of one-period valuations. For mathematical properties of general multiperiod valuations, we refer to
Cheridito et al \cite{Cheridito-et-al-06}, Artzner et al. \cite{Artzner-et-al-07} and Jobert and Rogers \cite{JobertRogers-08}.

We now explain the basic ingredients of multi-period cost-of-capital valuation.
Let $V_0$ be the value at time $0$ of the liability cashflow $X$. This amount should be interpreted as the capital that needs to be transferred along with the liability to another external agent in order for the external agent (or capital provider) to accept managing the liability runoff and the associated capital costs.
Let $V_{t}$ denote the value at time $t$ of the liability cashflows $X_{>t}$. The capital $V_{t}$ is reserved at time $t$ for managing the liability. However, regulation requires the capital $R_t=\rho(-X_{t+1}-V_{t+1}\mid \mathcal{F}_t)>V_{t}$ to be set aside at time $t$. The difference $R_t-V_t$ is made available by a capital provider requiring an excess expected return $1+\eta_t$ on the provided capital between time $t$ and $t+1$. The acceptability criterion under which the capital provider accepts to provide capital gives the equation
\begin{align}\label{eq:accr_CoCnoLL}
\operatorname{E}[R_t-X_{t+1}-V_{t+1}\mid \mathcal{F}_t]=(1+\eta_t)(R_t-V_t).
\end{align}
Solving for $V_t$ yields $V_t=\phi_t(X_{t+1}+V_{t+1})$, where
\begin{align}\label{eq:phiCoCnoLL}
\phi_t(Y)=\frac{1}{1+\eta_t}\operatorname{E}[Y\mid \mathcal{F}_t]+\frac{\eta_t}{1+\eta_t}\rho(-Y\mid \mathcal{F}_t).
\end{align}
A commonly used conditional risk measure for variables with finite variance is $\rho^{\text{SD}}(Y\mid \mathcal{F}_t)=\operatorname{E}[-Y\mid \mathcal{F}_t]+c^{\text{SD}}\operatorname{var}(Y\mid \mathcal{F}_t)^{1/2}$. Note that it can be argued that $\rho^{\text{SD}}(Y\mid \mathcal{F}_t)$ is an inappropriate choice because it violates the monotonicity property: $\widetilde{Y}\geq Y$ does not imply $\rho^{\text{SD}}(Y\mid \mathcal{F}_t)\geq \rho^{\text{SD}}(\widetilde{Y}\mid \mathcal{F}_t)$. For this choice of conditional risk measure, \eqref{eq:phiCoCnoLL} takes the form
\begin{align}\label{eq:phi_esd}
\phi_t(Y)=\operatorname{E}[Y\mid \mathcal{F}_t]+c_t\operatorname{var}(Y\mid \mathcal{F}_t)^{1/2},
\end{align}
where $c_t=c^{\text{SD}}\eta_t/(1-\eta_t)$, and $\phi_t(Y)\in L^2(\mathcal{F}_t)$ if $Y\in L^2(\mathcal{F}_T)$.

It can be argued that the capital provider, seen as the share holder of the company, has limited liability and is not required to continue injecting capital if the the value $R_t$ of available asset turns out insufficient to match the value $X_{t+1}+V_{t+1}$ of the liability towards the policy holders.
For discussions on limited liability in the context of cost-of-capital valuation we refer to Albrecher et al. \cite{Albrecher-et-al-22} and M\"ohr \cite{mhr11}.
In the setting with limited liability, \eqref{eq:accr_CoCnoLL} is replaced by
\begin{align}\label{eq:accr_CoCLL}
\operatorname{E}[(R_t-X_{t+1}-V_{t+1})^{+}\mid \mathcal{F}_t]=(1+\eta_t)(R_t-V_t)
\end{align}
and \eqref{eq:phiCoCnoLL} is replaced by
\begin{align}\label{eq:phiCoCLL}
\phi_t(Y)=\rho(-Y\mid \mathcal{F}_t)-\frac{1}{1+\eta_t}\operatorname{E}[(\rho(-Y\mid\mathcal{F}_t)-Y)^{+}\mid \mathcal{F}_t].
\end{align}
Regardless of whether we consider cost-of-capital valuation with or without limited liability, it follows that $\phi_t(aY+\widetilde{Y})=a\phi_t(Y)+\widetilde{Y}$ for $a\in\mathbb{R}_+$ and $\widetilde{Y}\in L^1(\mathcal{F}_t)$ if the conditional risk measure satisfies $\rho(aY+\widetilde{Y}\mid \mathcal{F}_t)=a\rho(Y\mid \mathcal{F}_t)-\widetilde{Y}$. Moreover, $\phi_t$ inherits monotonicity from the conditional risk measure: if $\widetilde{Y}\geq Y$ implies $\rho(Y\mid \mathcal{F}_t)\geq \rho(\widetilde{Y}\mid \mathcal{F}_t)$, then $\phi_t(\widetilde{Y})\geq \phi_t(Y)$.

An alternative to an acceptability criterion based on expected excess return and a cost-of-capital rate is an acceptability criterion saying that a risk averse capital provider provides capital if the payoff resulting from providing capital is preferred, in terms of expected utility, to simply rolling this capital forward by investing it in a riskless bond (or more generally, investing it in the numeraire asset). In this setting, the acceptability criterion \eqref{eq:accr_CoCLL} is replaced by
\begin{align}\label{eq:accr_EULL}
\operatorname{E}[u_t((R_t-X_{t+1}-V_{t+1})^{+})\mid \mathcal{F}_t]=u_t(R_t-V_t),
\end{align}
where $u_t$ is an increasing and concave (utility) function.
Consequently, \eqref{eq:phiCoCLL} is replaced by
\begin{align}\label{eq:phiEULL}
\phi_t(Y)=\rho(-Y\mid \mathcal{F}_t)-u_t^{-1}\Big(\operatorname{E}\Big[u_t((\rho(-Y\mid\mathcal{F}_t)-Y)^{+})\mid\mathcal{F}_t\Big]\Big).
\end{align}
Also in this case $\phi_t$ inherits monotonicity and the property $\phi_t(Y+\widetilde{Y})=\phi_t(Y)+\widetilde{Y}$, $\widetilde{Y}\in L^1(\mathcal{F}_t)$, from a monotone and conditionally cash additive conditional risk measure. However, the property $\phi_t(aY+\widetilde{Y})=a\phi_t(Y)+\widetilde{Y}$, $a\in\mathbb{R}_+$, requires that $u_t$ is chosen a power utility function $u_t(x)=\alpha_t x^{\beta_t}$, where $\alpha_t>0$ and $\beta_t\in (0,1]$. In this case, $\phi_t$ in \eqref{eq:phiEULL} takes the form
\begin{align}\label{eq:phiPULL}
\phi_t(Y)=\rho(-Y\mid \mathcal{F}_t)-\operatorname{E}\Big[\big((\rho(-Y\mid\mathcal{F}_t)-Y)^{+}\big)^{\beta_t}\mid\mathcal{F}_t\Big]^{1/\beta_t}.
\end{align}
The mappings $\phi_t:L^1(\mathcal{F}_T)\to L^1(\mathcal{F}_t)$ given by \eqref{eq:phiCoCnoLL}, \eqref{eq:phiCoCLL} or \eqref{eq:phiPULL} all satisfy $\phi_t(Y+\widetilde{Y})=\phi_t(Y)+\widetilde{Y}$ whenever $\widetilde{Y}\in L^1(\mathcal{F}_t)$. Therefore, \eqref{eq:VtVal} can be expressed as, with $\circ$ denoting composition,
\begin{align*}
V_T(X)=0, \quad V_t(X)=\phi_t \circ \dots \circ \phi_{T-1}\bigg(\sum_{s=t+1}^{T}X_s\bigg), \quad t < T.
\end{align*}
These properties also hold for the mappings $\phi_t$ in \eqref{eq:phi_esd} with $L^1$ replaced by $L^2$.

\subsection{Main results}\label{sec:mainresults}

The main results consist of Theorems \ref{thm:convergentvaluations}, \ref{thm:CondExpVar} and \ref{thm:V0monotonicity}.
Theorems \ref{thm:convergentvaluations} and \ref{thm:CondExpVar} presents conditions under which we have convergence for multi-period values (that are typically not computable) to a computable explicit limit. Theorem \ref{thm:V0monotonicity} presents a monotonicity result for the limit expressions in terms of a partial order between filtrations.

\begin{theorem}\label{thm:convergentvaluations}
Let $X,X^1,X^2,\dots$ be random vectors in $\mathbb{R}^T$ and let $Y,Y^1,Y^2,\dots$ be random vectors in $(\mathbb{R}^d)^T$, $d\in \mathbb{N}$.
Suppose that $\mathcal{L}(X^n,Y^n)\stackrel{w}{\to} \mathcal{L}(X,Y)$ as $n\to\infty$, where $\mathcal{L}(X,Y)$ is Gaussian, and that $(\mathcal{L}(\Vert X^n\Vert))_{n\in\mathbb{N}}$ is uniformly integrable.
Suppose also that, for each $t$ and each convergent sequence $(x^n,y^n)_{\leq t}\to (x,y)_{\leq t}$ with $(x^n,y^n)_{\leq t}\in \operatorname{supp}(\mathcal{L}((X^n,Y^n)_{\leq t}))$,
\begin{align}
& \mathcal{L}\Big((X^n,Y^n) \mid (X^n,Y^n)_{\leq t}=(x^n,y^n)_{\leq t}\Big) \stackrel{w}{\to} \mathcal{L}\Big((X,Y) \mid (X,Y)_{\leq t}=(x,y)_{\leq t}\Big)
\quad \text{as } n\to\infty \label{eq:contconv}, \\
& \Big(\mathcal{L}\big(\Vert X^n\Vert \mid (X^n,Y^n)_{\leq t}=(x^n,y^n)_{\leq t}\big)\Big)_{n\in\mathbb{N}} \quad \text{is uniformly integrable}. \label{eq:condxnui}
\end{align}
For all $t$, let $V_t(X)$ and $V^n_t(X^n)$ be given by \eqref{eq:VtVal} and \eqref{eq:VntVal} with $\phi_t$ and $\phi^n_t$ given by either
\eqref{eq:phiCoCnoLL} and
\begin{align}\label{eq:phiCoCnoLL_n}
\phi^n_t(Y)=\frac{1}{1+\eta_t}\operatorname{E}[Y\mid \mathcal{F}^n_t]+\frac{\eta_t}{1+\eta_t}\rho(-Y\mid \mathcal{F}^n_t), \quad \eta_t \in \mathbb{R}_+,
\end{align}
or by \eqref{eq:phiCoCLL} and
\begin{align}\label{eq:phiCoCLL_n}
\phi^n_t(Y)=\rho(-Y\mid \mathcal{F}^n_t)-\frac{1}{1+\eta_t}\operatorname{E}\Big[\big(\rho(-Y\mid\mathcal{F}^n_t)-Y\big)^{+}\mid \mathcal{F}^n_t\Big], \quad \eta_t \in \mathbb{R}_+,
\end{align}
or by \eqref{eq:phiPULL} and
\begin{align}\label{eq:phiPULL_n}
\phi^n_t(Y)=\rho(-Y\mid \mathcal{F}^n_t)-\operatorname{E}\Big[\big((\rho(-Y\mid\mathcal{F}^n_t)-Y)^{+}\big)^{\beta_t}\mid\mathcal{F}^n_t\Big]^{1/\beta_t}, \quad \beta_t\in (0,1],
\end{align}
where $\mathcal{F}_t=\sigma((X,Y)_{\leq t})$ and $\mathcal{F}^n_t=\sigma((X^n,Y^n)_{\leq t})$, and where $\rho(-Y\mid \mathcal{F}_t)$ and $\rho(-Y\mid \mathcal{F}^n_t)$ are of the form \eqref{eq:quantriskmeas}.
Then $\lim_{n\to\infty}V^n_0(X^n)=V_0(X)$, where
\begin{align}\label{eq:V0expression}
V_0(X)=\operatorname{E}\bigg[\sum_{t=1}^T X_t\bigg]+\sum_{t=1}^T\phi_{t-1}(\varepsilon_t)\bigg(\operatorname{Var}\bigg(\sum_{u=t}^TX_u\mid (X,Y)_{\leq t-1}\bigg)
-\operatorname{Var}\bigg(\sum_{u=t}^TX_u\mid (X,Y)_{\leq t}\bigg)\bigg)^{1/2},
\end{align}
where $\varepsilon_t$ is standard normally distributed and independent of $\mathcal{F}_{t-1}$.
\end{theorem}

\begin{remark}
Notice that the conditional variances in \eqref{eq:V0expression} are nonrandom since the (joint) distribution of $(X,Y)$ is Gaussian.
Notice also $\mathcal{F}_t=\sigma((X,Y)_{\leq t})$ generates the natural filtration of $X$ if $Y$ is chosen as nonrandom (a degenerate Gaussian process). We emphasize that if the $\eta_t$ in \eqref{eq:phiCoCnoLL} and \eqref{eq:phiCoCLL} do not depend on $t$, and similarly if the $\beta_t$ in \eqref{eq:phiPULL_n} do not depend on $t$, then $\phi_{t-1}(\varepsilon_t)=\phi_0(\varepsilon_1)$ in \eqref{eq:V0expression} does not depend on $t$.
\end{remark}

\begin{remark}
Joint weak convergence $\mathcal{L}(X^n,Y^n)\stackrel{w}{\to} \mathcal{L}(X,Y)$ as $n\to\infty$ does not imply conditional weak convergence such as the convergence in \eqref{eq:contconv}. A counterexample as well as sufficient conditions for conditional weak convergence are presented in \cite{Sweeting-89}.
\end{remark}

\begin{remark}
For $\varepsilon$ is standard normally distributed and independent of $\mathcal{F}_{t-1}$,
\begin{align*}
\operatorname{V@R}_{u}(\varepsilon \mid \mathcal{F}_{t-1}) = \Phi^{-1}(1-u), \quad \operatorname{AV@R}_{u}(\varepsilon \mid \mathcal{F}_{t-1}) = \frac{1}{u}\phi\left(\Phi^{-1}(1-u)\right),
\end{align*}
where $\phi$ and $\Phi$ here denote the standard normal density and distribution function, respectively. Hence, with $\rho_0(\varepsilon)$ denoting either $\operatorname{V@R}_{u}(\varepsilon \mid \mathcal{F}_{t-1})$ or $\operatorname{AV@R}_{u}(\varepsilon \mid \mathcal{F}_{t-1})$ with values above that do not depend on $t$, we get
$$
\phi_{t-1}(\varepsilon) = \frac{\eta_{t-1}}{1 + \eta_{t-1}} \rho_0(\varepsilon)
$$
in the case of \eqref{eq:phiCoCnoLL}, and
\begin{align*}
\phi_{t-1}(\varepsilon) = \rho_0(\varepsilon) - \frac{1}{1 + \eta_{t-1}}\Big(\rho_0(\varepsilon)\Phi(\rho_0(\varepsilon)) + \phi(\rho_0(\varepsilon)) \Big)\leq \frac{\eta_{t-1}}{1 + \eta_{t-1}} \rho_0(\varepsilon)
\end{align*}
in the case of \eqref{eq:phiCoCLL}. In case of \eqref{eq:phiPULL}, the expectation $\operatorname{E}[((\rho_0(\varepsilon)-\varepsilon)^+)^{\beta_{t-1}}]$ has to be computed numerically in order to compute $\phi_{t-1}(\varepsilon)$.
\end{remark}

\begin{remark}\label{rem:noll}
As stated in Theorem \ref{thm:convergentvaluations}, the result holds when $\phi_t$ and $\phi^n_t$ are given by \eqref{eq:phiCoCnoLL} and \eqref{eq:phiCoCnoLL_n}. However, the theorem is actually proven for mappings of a more general kind, namely, for $\lambda_t\in [0,1]$ and $\mu^1_t,\mu^2_t\in \mathcal{P}([0,1])'$, the mappings $\phi_t : L^1(\mathcal{F}_{T})\to L^1(\mathcal{F}_{t})$, $\phi^n_t : L^1(\mathcal{F}^n_{T})\to L^1(\mathcal{F}^n_{t})$  given by \eqref{eq:phialtmapping} and \eqref{eq:phinaltmapping} below:
\begin{align}
\phi_t(Y)&=\lambda_t\int F^{-1}_{Y \mid \mathcal{F}_t}(p)\mu_t^1(dp)+(1-\lambda_t)\int F^{-1}_{-Y \mid \mathcal{F}_t}(p)\mu_t^2(dp), \label{eq:phialtmapping}\\
\phi^n_t(Y)&=\lambda_t\int F^{-1}_{Y \mid \mathcal{F}^n_t}(p)\mu_t^1(dp)+(1-\lambda_t)\int F^{-1}_{-Y \mid \mathcal{F}^n_t}(p)\mu_t^2(dp).   \label{eq:phinaltmapping}
\end{align}
Note that \eqref{eq:phiCoCnoLL} and \eqref{eq:phiCoCnoLL_n} correspond to $\mu_t^1(dp)=dp$ and $\lambda_t=(1+\eta_t)^{-1}$.
\end{remark}

\begin{remark}\label{rem:ll}
Both \eqref{eq:phiCoCLL}, \eqref{eq:phiCoCLL_n} and \eqref{eq:phiPULL}, \eqref{eq:phiPULL_n} are special cases of
\begin{align}
\phi_t(Y)&=\rho(-Y\mid \mathcal{F}_t)-\gamma_t\operatorname{E}\Big[\big((\rho(-Y\mid\mathcal{F}_t)-Y)^{+}\big)^{\beta_t}\mid\mathcal{F}_t\Big]^{1/\beta_t}, \label{eq:phiPULLgen}\\
\phi^n_t(Y)&=\rho(-Y\mid \mathcal{F}^n_t)-\gamma_t\operatorname{E}\Big[\big((\rho(-Y\mid\mathcal{F}^n_t)-Y)^{+}\big)^{\beta_t}\mid\mathcal{F}^n_t\Big]^{1/\beta_t}, \label{eq:phiPULLgen_n}
\end{align}
as seen by choosing $(\gamma_t,\beta_t)=((1+\eta_t)^{-1},1)$ in case of   \eqref{eq:phiCoCLL}, \eqref{eq:phiCoCLL_n}, and $\gamma_t=1$ in case of \eqref{eq:phiPULL}, \eqref{eq:phiPULL_n}.
\end{remark}

Our next result shows convergence of multi-period values in the case of where the one period valuation mapping is defined as a sum of a conditional expectation and a constant times a conditional standard deviation.
Let $c_t \in \mathbb{R}_+$ and
$\phi_t : L^2(\mathcal{F}_{T})\to L^2(\mathcal{F}_{t})$, $\phi^n_t : L^2(\mathcal{F}^n_{T})\to L^2(\mathcal{F}^n_{t})$ be given by
\begin{align}
\phi_t(Y)&=\operatorname{E} \big[Y \mid \mathcal{F}_t \big] + c_t \operatorname{var} \big[Y \mid \mathcal{F}_t \big]^{1/2}, \label{eq:phimappingII} \\
\phi^n_t(Y)&=\operatorname{E} \big[Y \mid \mathcal{F}^n_t \big] + c_t \operatorname{var} \big[Y \mid \mathcal{F}^n_t \big]^{1/2}. \label{eq:phinmappingII}
\end{align}
$(V_t(X))_{t=0}^T$ and $(V^n_t(X^n))_{t=0}^T$ are defined backward recursively by \eqref{eq:VtVal} and \eqref{eq:VntVal}, and \eqref{eq:phitab} and \eqref{eq:Vtab} hold.

\begin{theorem}\label{thm:CondExpVar}
Let $X,X^1,X^2,\dots$ be random vectors in $\mathbb{R}^T$ and let $Y,Y^1,Y^2,\dots$ be random vectors in $(\mathbb{R}^d)^T$, $d\in\mathbb{N}$. Suppose that $\mathcal{L}(X^n,Y^n)\stackrel{w}{\to} \mathcal{L}(X,Y)$ as $n\to\infty$, where $\mathcal{L}(X,Y)$ is Gaussian, and such that $(\mathcal{L}(\Vert X^n\Vert^2))_{n\in\mathbb{N}}$ is uniformly integrable. Suppose also that, for each $t$ and each convergent sequence $(x^n,y^n)_{\leq t}\to (x,y)_{\leq t}$ with $(x^n,y^n)_{\leq t}\in \operatorname{supp}(\mathcal{L}((X^n,Y^n)_{\leq t}))$,
\begin{align}
& \mathcal{L}\Big((X^n,Y^n) \mid (X^n,Y^n)_{\leq t}=(x^n,y^n)_{\leq t}\Big) \stackrel{w}{\to} \mathcal{L}\Big((X,Y) \mid (X,Y)_{\leq t}=(x,y)_{\leq t}\Big) \quad \text{as } n\to\infty, \label{eq:contconvII} \\
& \Big(\mathcal{L}(\Vert X^n\Vert^2 \mid (X^n,Y^n)_{\leq t}=(x^n,y^n)_{\leq t})\Big)_{n\in\mathbb{N}} \quad \text{is uniformly integrable}. \label{eq:condxnuiII}
\end{align}
Let $V_t(X)$ and $V^n_t(X^n)$ be given by \eqref{eq:VtVal} and \eqref{eq:VntVal} with $\phi_t$ and $\phi^n_t$ given by \eqref{eq:phimappingII} and \eqref{eq:phinmappingII} with $\mathcal{F}_t=\sigma((X,Y)_{\leq t})$ and $\mathcal{F}^n_t=\sigma((X^n,Y^n)_{\leq t})$.
Then $\lim_{n\to\infty}V^n_0(X^n)=V_0(X)$, where
\begin{align*}
V_0(X)=\operatorname{E}\bigg[\sum_{t=1}^T X_t\bigg]+\sum_{t=1}^T\phi_{t-1}(\varepsilon_t)\bigg(\operatorname{Var}\bigg(\sum_{u=t}^TX_u\mid (X,Y)_{\leq t-1}\bigg)
-\operatorname{Var}\bigg(\sum_{u=t}^TX_u\mid (X,Y)_{\leq t}\bigg)\bigg)^{1/2},
\end{align*}
where $\varepsilon_t$ is standard normally distributed and independent of $\mathcal{F}_{t-1}$.
\end{theorem}

The final theorem considers the Gaussian limit models appearing in Theorems \ref{thm:convergentvaluations} and \ref{thm:CondExpVar} and establishes a partial ordering of values $V_0(X)$ assigned to a given cashflow $X$ depending on a partial order between filtrations to which $X$ is adapted. One may guess that a larger filtration allows for more accurate predictions over time of future cashflows and thereby reduces the value $V_0(X)$. We show that this holds if the parameters of the one-step valuation mappings $\phi_t$ do not depend on $t$ and if the conditional variance of the aggregate cashflow $\sum_{s=1}^TX_s$ decays in a convex manner over time as more information becomes available.

We consider filtrations $\mathbb{F}=(\mathcal{F}_t)_{t=0}^{T}$ with $\mathcal{F}_0=\{\Omega,\emptyset\}$ and $\mathcal{F}_t=\sigma((X,Y)_{\leq t})$, for $t\geq 1$, where $(X,Y)$ is jointly Gaussian. We assume that the one-step valuation mappings $\phi_t$ do not depend on $t$ such that $\phi_{t-1}(\varepsilon_t)=\phi_0(\varepsilon_1)$ in \eqref{eq:V0expression}. We write
\begin{align*}
V_0(X,\mathbb{F})=\operatorname{E}\bigg[\sum_{t=1}^T X_t\bigg]+\phi_{0}(\varepsilon_1)\sum_{t=1}^T\bigg(&\operatorname{Var}\bigg(\sum_{u=t}^TX_u\mid \mathcal{F}_{t-1}\bigg)
-\operatorname{Var}\bigg(\sum_{u=t}^TX_u\mid \mathcal{F}_{t}\bigg)\bigg)^{1/2}.
\end{align*}
The following theorem says that for filtrations $\mathbb{F}$ and $\mathbb{G}$ of the above kind (generated by a process that is jointly Gaussian together with the cashflow process $X$), if $t\mapsto \operatorname{Var}(\sum_{s=1}^TX_s \mid \mathcal{F}_t)$ is convex, then $V_0(X,\mathbb{F})\geq V_0(X,\mathbb{G})$ whenever $\mathcal{F}_t\subseteq \mathcal{G}_t$ for every $t$. The convexity assumption is necessary.

\begin{theorem}\label{thm:V0monotonicity}
Let $(X,Y)_{t=1}^T$ and $(X,Z)_{t=1}^T$ be Gaussian processes.
Let $\mathbb{F}=(\mathcal{F}_t)_{t=0}^{T}$ and $\mathbb{G}=(\mathcal{G}_t)_{t=0}^{T}$ be filtrations with $\mathcal{F}_0=\mathcal{G}_0=\{\Omega,\emptyset\}$ and $\mathcal{F}_t=\sigma((X,Y)_{\leq t})$ and $\mathcal{G}_t=\sigma((X,Z)_{\leq t})$, for $t\geq 1$.
If $\mathcal{F}_t\subseteq \mathcal{G}_t$ for every $t$ and if $t\mapsto \operatorname{Var}(\sum_{s=1}^TX_s \mid \mathcal{F}_t)$ is convex, then $V_0(X,\mathbb{F})\geq V_0(X,\mathbb{G})$.
\end{theorem}

\begin{remark}
We emphasize that the Gaussian assumption in Theorem \ref{thm:V0monotonicity} means that the terms
\begin{align*}
c_t=\operatorname{Var}\bigg(\sum_{s=1}^TX_s \mid \mathcal{F}_{t-1}\bigg)-\operatorname{Var}\bigg(\sum_{s=1}^TX_s \mid \mathcal{F}_{t}\bigg)
\end{align*}
form a nonrandom sequence $(c_t)_{t=1}^{T}$ and that the convexity property is equivalent to
$c_1\geq c_2 \geq \dots \geq c_{T}$.
\end{remark}

\section{Proofs and auxiliary results}\label{sec:proofs}

We use the notation $Z^n_{\leq t}=(X^n,Y^n)_{\leq t}$, $Z_{\leq t}=(X,Y)_{\leq t}$ and $z^n_{\leq t}=(x^n,y^n)_{\leq t}$, $z_{\leq t}=(x,y)_{\leq t}$.

\begin{proof}[Proof of Theorem \ref{thm:convergentvaluations} assuming \eqref{eq:phialtmapping} and \eqref{eq:phinaltmapping}]
As explained in Remark \ref{rem:noll}, \eqref{eq:phiCoCnoLL} and \eqref{eq:phiCoCnoLL_n} are special cases of \eqref{eq:phialtmapping} and \eqref{eq:phinaltmapping}.

For $t=1,\dots,T$, consider mappings $\psi^n_t:\operatorname{supp}(\mathcal{L}(Z^n_{\leq t}))\to \mathbb{R}$ given by $\psi^n_T(z^n_{\leq T})=\sum_{s=1}^Tx^n_s$ and
\begin{align}\label{eq:cvfun_ta}
\psi^n_t(z^n_{\leq t})
=\lambda_t\int F^{-1}_{\psi^n_{t+1}(Z^n_{\leq t+1}) \mid Z^n_{\leq t} = z^n_{\leq t}}(p)\mu_t^1(dp)+(1-\lambda_t)\int F^{-1}_{-\psi^n_{t+1}(Z^n_{\leq t+1}) \mid Z^n_{\leq t} = z^n_{\leq t}}(p)\mu_t^2(dp),
\end{align}
and let
\begin{align}\label{eq:cvfun_0a}
\psi^n_0 =  \lambda_0\int F^{-1}_{\psi^n_{1}(Z^n_{1})}(p)\mu_0^1(dp)+(1-\lambda_0)\int F^{-1}_{-\psi^n_{1}(Z^n_{1})}(p)\mu_0^2(dp),
\end{align}
where $\mu_t^1, \mu_t^2 \in \mathcal{P}([0,1])'$. We define $\psi_t:\mathbb{R}^t\to \mathbb{R}$ and $\psi_0$ analogously without the superscript $n$.
Note that
\begin{align*}
\psi^n_t(Z^n_{\leq t})&=\sum_{s=1}^{t}X^n_s+V^n_t(X^n)=\phi^n_t(\psi^n_{t+1}(Z^n_{\leq t+1})), \; t\geq 1,
\quad \psi^n_0=V^n_0(X)=\phi^n_0(\psi^n_1(Z^n_{\leq 1})), \\
\psi_t(Z_{\leq t})&=\sum_{s=1}^{t}X_s+V_t(X)=\phi_t(\psi_{t+1}(Z_{\leq t+1})), \; t\geq 1,
\quad \psi_0=V_0(X)=\phi_0(\psi_1(Z_{\leq 1})).
\end{align*}
Therefore, the proof is complete once we show the convergence $\lim_{n\to\infty}\psi^n_0=\psi_0$.

The argument of the proof is backwards induction. We will show that for $t = T, \dots, 1$,
\begin{itemize}
\item[{\bf UI($t$):}] $(\mathcal{L}(\psi^n_t(Z^n_{\leq t}) \mid Z^n_{\leq s} = z^n_{\leq s}))_n$ is uniformly integrable for each $s = 1, \dots, t-1$ and $z^n_{\leq s} \to z_{\leq s}$,
\item[{\bf CC($t$):}] $(\psi^n_t)_n$ is continuously convergent, i.e. $\psi^n_t(z^n_{\leq t}) \to \psi_t(z_{\leq t})$, whenever $z^n_{\leq t} \to z_{\leq t}$.
\end{itemize}

\noindent
{\bf Induction base:} {\bf UI($T$)} and {\bf CC($T$)} hold.

For $s \in \{1, \dots, T-1\}$ and $z^n_{\leq s} \to z_{\leq s}$, {\bf UI($T$)} follows from \eqref{eq:condxnui} since norms on Euclidean spaces are equivalent and $|\psi^n_T(z^n_{\leq T})| = |\sum_{s=1}^T x^n_s| \leq \sum_{s=1}^T |x^n_s|$.
Moreover,
$$
\psi^n_T(z^n_{\leq T}) = \sum_{s=1}^T x^n_s \to \sum_{s=1}^T x_s = \psi_T(z_{\leq T})
$$
whenever $z^n_{\leq T} \to z_{\leq T}$, yielding {\bf CC($T$)}.

\noindent
{\bf Induction step:} {\bf UI($t+1$)}, {\bf CC($t+1$)} together imply {\bf UI($t$)}, {\bf CC($t$)}.

Fix $s\in \{1,\dots,t-1\}$ and sequence $(z^n_{\leq s})_n$ with $z^n_{\leq s}\to z_{\leq s}$.
Let $\kappa^n_{s+1,t+1}(z^n_{\leq s},\cdot)$ and $\kappa_{s+1,t+1}(z_{\leq s},\cdot)$ be regular versions of the conditional distributions
$$
\operatorname{P}((Z^n_{s+1},\dots,Z^n_{t+1})^{{\mathrm{T}}} \in \cdot \mid Z^n_{\leq s} = z^n_{\leq s})
\quad\text{and}\quad
\operatorname{P}((Z_{s+1},\dots,Z_{t+1})^{{\mathrm{T}}} \in \cdot \mid Z_{\leq s} = z_{\leq s}).
$$
By assumption \eqref{eq:contconv}, $\kappa^n_{s+1,t+1}(z^n_{\leq s},\cdot)\stackrel{w}{\to} \kappa_{s+1,t+1}(z_{\leq s},\cdot)$.
Since $(\psi^n_{t+1}(z^n_{\leq s}, \cdot ))$ is continuously convergent by the induction assumption {\bf CC($t+1$)}, by the generalized continuous mapping theorem (\cite[Theorem 4.27, Exercise 27]{Kallenberg-02}),
$$
\kappa^n_{s+1,t+1}(z^n_{\leq s}, \cdot) \circ (\psi^n_{t+1}(z^n_{\leq s}, \cdot))^{-1} \stackrel{w}{\to} \kappa_{s+1,t+1}(z_{\leq s}, \cdot) \circ (\psi_{t+1}(z_{\leq s}, \cdot))^{-1},
$$
i.e.
\begin{align}\label{conv:Weak123}
\mathcal{L}(\psi^n_{t+1}(Z_{\leq t+1}^n) \mid Z^n_{\leq s} = z^n_{\leq s}) \stackrel{w}{\to} \mathcal{L}(\psi_{t+1}(Z_{\leq t+1}) \mid Z_{\leq s} = z_{\leq s}).
\end{align}
Lemma \ref{lem:gaussian_value} shows that the limiting distribution in \eqref{conv:Weak123} is Gaussian.
By Lemma \ref{lem:BoundRec}, there exists a constant $c > 0$, such that
$$
|\psi^n_t(z)| \leq c \operatorname{E}[|\psi^n_{t+1}(Z^n_{\leq t+1})| \mid Z^n_{\leq t} = z].
$$
Hence, it is sufficient to show uniform integrability of
\begin{align}\label{eq:UI23}
\Big(\mathcal{L}\big(\operatorname{E}[|\psi^n_{t+1}(Z^{n}_{\leq t+1})| \mid Z^n_{\leq t}] \mid Z^n_{\leq s}=z^n_{\leq s}\big)\Big)_n.
\end{align}
Setting $\widetilde{W}^{n}=\operatorname{E}[|\psi^{n}_{t+1}(Z^n_{\leq t+1})| \mid Z^n_{\leq t}]$ and $W^{n}=|\psi^{n}_{t+1}(Z^n_{\leq t+1})|$, this follows from Lemma \ref{lem:TowerUI}, using {\bf UI($t+1$)} and {\bf CC($t+1$)}. Hence, {\bf UI($t$)} holds.

We proceed by showing {\bf CC($t$)}, i.e.~the convergence of
\begin{align}
\psi^n_t(z^n_{\leq t}) &=
\lambda_t\int F^{-1}_{\psi^n_{t+1}(Z^n_{\leq t+1}) \mid Z^n_{\leq t} = z_{\leq t}^n }(p) \mu_t^1(dp) \label{eq:RecConv2a} \\
& \quad + (1-\lambda_t)\int F^{-1}_{-\psi^n_{t+1}(Z^n_{\leq t+1}) \mid Z^n_{\leq t} = z^n_{\leq t} }(p) \mu_t^2(dp). \label{eq:RecConv2b}
\end{align}
We will prove the convergence assuming that $\mu^1_t$ has a bounded density and that $\operatorname{supp}(\mu^2_t)\subset [a,b]\subset (0,1)$. Other possibilities for $\mu^1_t,\mu^2_t\in \mathcal{P}([0,1])'$ are handled by the same arguments as those shown below.

We start by the integral in \eqref{eq:RecConv2a} and will show convergence by an application of Pratt's Lemma \cite{Pratt-60}.
Let $w : (0, 1) \to \mathbb{R}$, bounded by some $c' > 0$, be the density of $\mu_t^1$, i.e. $\mu_t^1(dp) = w(p)dp$. Let
$$
f_n(p) =  F^{-1}_{\psi^n_{t+1}(Z^n_{\leq t+1}) \mid Z^n_{\leq t}= z^n_{\leq t}}(p) w(p).
$$
From \eqref{conv:Weak123} follows convergence
\begin{align*}
f_n(p) \to F^{-1}_{\psi_{t+1}(Z_{\leq t+1}) \mid Z_{\leq t}= z_{\leq t}}(p)w(p)
\end{align*}
for almost every $p \in (0,1)$.
By combining \eqref{conv:Weak123} and the continuous mapping theorem applied to the absolute value function follows convergence of upper and lower bounds $l_n(p) \leq f_n(p) \leq u_n(p)$ for almost every $p\in (0,1)$:
\begin{align*}
l_n(p) &= c'F^{-1}_{-|\psi^n_{t+1}(Z^n_{\leq t+1})| \mid Z^n_{\leq t} = z^n_{\leq t} }(p) \to  c'F^{-1}_{-|\psi_{t+1}(Z_{\leq t+1})| \mid Z_{\leq t} = z_{\leq t} }(p) = l(p) \\
u_n(p) &= c'F^{-1}_{|\psi^n_{t+1}(Z^n_{\leq t+1})| \mid Z^n_{\leq t} = z^n_{\leq t} }(p) \to  c'F^{-1}_{|\psi_{t+1}(Z_{\leq t+1})| \mid Z_{\leq t} = z_{\leq t} }(p) = u(p).
\end{align*}
The induction assumption {\bf UI($t+1$)} together with
\begin{align*}
\mathcal{L}(|\psi^n_{t+1}(Z^n_{\leq t+1})| \mid Z^n_{\leq t} = z^n_{\leq t})
\stackrel{w}{\to} \mathcal{L}(|\psi_{t+1}(Z_{\leq t+1})| \mid Z_{\leq t} = z_{\leq t}),
\end{align*}
allow us to conclude from \cite[Theorem 4.11]{Kallenberg-02} that
\begin{align*}
\int l_n(p) dp = &-c'\operatorname{E}[|\psi^n_{t+1}(Z_{\leq t+1}^n)| \mid Z_{\leq t}^n = z_{\leq t}^n] \\ \to
&-c'\operatorname{E}[|\psi_{t+1}(Z_{\leq t+1})| \mid Z_{\leq t} = z_{\leq t}] = \int l(p) dp
\end{align*}
and
\begin{align*}
\int u_n(p) dp = \; &c'\operatorname{E}[|\psi^n_{t+1}(Z_{\leq t+1}^n)| \mid Z_{\leq t}^n = z_{\leq t}^n] \\ \to \;
&c'\operatorname{E}[|\psi_{t+1}(Z_{\leq t+1})| \mid Z_{\leq t} = z_{\leq t}] = \int u(p) dp.
\end{align*}
Hence, by Pratt's Lemma \cite[Theorem 1]{Pratt-60} we have convergence of the integral in \eqref{eq:RecConv2a}, i.e.
\begin{align*}
\int F^{-1}_{\psi^n_{t+1}(Z^n_{\leq t+1}) \mid Z^n_{\leq t} = z_{\leq t}^n }(p) \mu_1(dp)
&= \int f_n(p)dp \to \int f(p) dp \\
&= \int F^{-1}_{\psi_{t+1}(Z_{\leq t+1}) \mid Z_{\leq t} = z_{\leq t} }(p) \mu_1(dp).
\end{align*}
We now consider the integral in \eqref{eq:RecConv2b}. By Lemma \ref{lem:gaussian_value}, the mapping $p \mapsto F^{-1}_{-\psi_{t+1}(Z_{\leq t+1}) \mid Z_{\leq t} = z_{\leq t} }(p)$ is continuous for all $p \in (0,1)$. Therefore, from \eqref{conv:Weak123} follows pointwise convergence
$$
F^{-1}_{-\psi^n_{t+1}(Z^n_{\leq t+1}) \mid Z^n_{\leq t} = z_{\leq t}^n }(p) \to F^{-1}_{-\psi_{t+1}(Z_{\leq t+1}) \mid Z_{\leq t} = z_{\leq t} }(p)
$$
for all $p \in (0,1)$. In particular, the convergence holds at the points $a, b \in (0,1)$. Therefore, using the fact that $p \mapsto F^{-1}_{-\psi^n_{t+1}(Z^n_{\leq t+1}) \mid Z^n_{\leq t} = z^n_{\leq t} }(p)$ and $p \mapsto F^{-1}_{-\psi_{t+1}(Z_{\leq t+1}) \mid Z_{\leq t} = z_{\leq t} }(p)$ are increasing, we can find $d > 0$, such that
$$
c'' = \max \bigg(\Big|F^{-1}_{-\psi_{t+1}(Z_{\leq t+1}) \mid Z_{\leq t} = z_{\leq t} }(a)\Big|, \Big|F^{-1}_{-\psi_{t+1}(Z_{\leq t+1}) \mid Z_{\leq t} = z_{\leq t} }(b)\Big| \bigg) + d
$$
uniformly bounds the integrand, i.e.
$$
\Big|F^{-1}_{-\psi^n_{t+1}(Z^n_{\leq t+1}) \mid Z^n_{\leq t} = z^n_{\leq t} }(p)\Big| \leq c''
$$
for all $p \in [a, b]$ and for all $n \in \mathbb{N}$. Hence, the bounded convergence theorem yields
\begin{align*}
\int F^{-1}_{-\psi^n_{t+1}(Z^n_{\leq t+1}) \mid Z^n_{\leq t} = z_{\leq t}^n }(p) \mu^2_t(dp)
\to \int F^{-1}_{-\psi_{t+1}(Z_{\leq t+1}) \mid Z_{\leq t} = z_{\leq t} }(p) \mu^2_t(dp).
\end{align*}
Hence, we have shown {\bf CC($t$)} and the proof of the induction step is complete.

It remains to show convergence of $\psi^n_0$. By the continuous convergence of $(\psi^n_1)$ and the generalized continuous mapping theorem (\cite[Theorem 4.27]{Kallenberg-02}), it follows that $\mathcal{L}(\psi^n_1(Z_1^n)) \stackrel{w}{\to} \mathcal{L}(\psi_1(Z_1))$, where the limit is Gaussian (cf. Lemma \ref{lem:gaussian_value}). Therefore the convergence of
\begin{align*}
\psi^n_0=
\lambda_0\int F^{-1}_{\psi^n_{1}(Z^n_{1})}(p) \mu_0^1(dp)
 + (1-\lambda_0)\int F^{-1}_{-\psi^n_{1}(Z^n_{1}) }(p) \mu_0^2(dp)
\end{align*}
follows from arguments completely analogous to those verifying the induction step. Lemma \ref{lem:gaussian_value} together with the variance decomposition
\begin{align*}
\operatorname{Var}\bigg(\operatorname{E}\bigg[\sum_{v=u}^{T}X_v\mid Z_{\leq u}\bigg]\mid Z_{\leq u-1}\bigg)
&=\operatorname{Var}\bigg(\sum_{v=u}^{T}X_v\mid Z_{\leq u-1}\bigg)\\
&\quad-\operatorname{E}\bigg[\operatorname{Var}\bigg(\sum_{v=u}^{T}X_v\mid Z_{\leq u}\bigg)\mid Z_{\leq u-1}\bigg]\\
&=\operatorname{Var}\bigg(\sum_{v=u}^{T}X_v\mid Z_{\leq u-1}\bigg)-\operatorname{Var}\bigg(\sum_{v=u}^{T}X_v\mid Z_{\leq u}\bigg)
\end{align*}
completes the proof.
\end{proof}

\begin{proof}[Proof of Theorem \ref{thm:convergentvaluations} assuming \eqref{eq:phiPULLgen} and \eqref{eq:phiPULLgen_n}]
 The proofs is similar to the above proof of Theorem \ref{thm:convergentvaluations} assuming \eqref{eq:phialtmapping} and \eqref{eq:phinaltmapping}.

Fix $\mu_0, \dots, \mu_{T-1} \in \mathcal{P}([0,1])'$.
Let $\nu^n_T(z^n_{\leq T})=\sum_{s=1}^T x^n_s$ and $\psi^n_T(z^n_{\leq T})=\sum_{s=1}^T x^n_s$.
For $t=1,\dots,T-1$, define mappings $\nu^n_t,\psi^n_t:\operatorname{supp}(\mathcal{L}(Z^n_{\leq t}))\to \mathbb{R}$ by
\begin{align}
\nu^n_{t}(z^n_{\leq t})&=\int F^{-1}_{-\psi^n_{t+1}(Z^n_{\leq t+1})\mid Z^n_{\leq t}=z^n_{\leq t}}(p)\mu_t(dp),  \label{eq:ConvRisk3}  \\
\psi^n_{t}(z^n_{\leq t})&=\nu^n_{t}(z^n_{\leq t})-\gamma_t\operatorname{E}\Big[\Big(\big(\nu^n_{t}(z^n_{\leq t})-\psi^n_{t+1}(Z^n_{\leq t+1})\big)^+\Big)^{\beta_t}\mid Z^n_{\leq t}=z^n_{\leq t}\Big]^{1/\beta_t}. \label{eq:ConvUtility3}
\end{align}
Let
\begin{align}\label{eq:ConvRU3}
\nu^n_{0}=\int F^{-1}_{-\psi^n_{1}(Z^n_{1})}(p)\mu_0(dp),\quad
\psi^n_{0}=\nu^n_{0}-\gamma_0\operatorname{E}\Big[\Big(\big(\nu^n_{0}-\psi^n_{1}(Z^n_{1})\big)^+\Big)^{\beta_0}\Big]^{1/\beta_0}.
\end{align}
Define $\nu_t,\psi_t:\mathbb{R}^t\to\mathbb{R}$ and $\nu_0,\psi_0$ analogously without the superscript $n$.
The proof is complete once we show the convergence $\lim_{n\to\infty}\psi^n_0=\psi_0$.

The argument of the proof is backwards induction. We will show that for $t = T, \dots, 1$,
\begin{itemize}
\item[{\bf UI($t$):}] $(\mathcal{L}(\psi^n_t(Z^n_{\leq t}) \mid Z^n_{\leq s} = x^n_{\leq s}))_n$ is uniformly integrable for each $s = 1, \dots, t-1$ and $z^n_{\leq s} \to z_{\leq s}$,
\item[{\bf CC($t$):}] $(\psi^n_t)_n$ is continuously convergent, i.e. $\psi^n_t(z^n_{\leq t}) \to \psi_t(z_{\leq t})$, whenever $z^n_{\leq t} \to z_{\leq t}$.
\end{itemize}

\noindent
{\bf Induction base:} {\bf UI($T$)} and {\bf CC($T$)} hold. The argument verifying the induction base is identical to that in the proofs of Theorem \ref{thm:convergentvaluations} assuming \eqref{eq:phialtmapping} and \eqref{eq:phinaltmapping}.

\noindent
{\bf Induction step:} {\bf UI($t+1$)}, {\bf CC($t+1$)} together imply {\bf UI($t$)}, {\bf CC($t$)}.
The argument verifying that
\begin{align}\label{conv:Weak123456}
\mathcal{L}(\psi^n_{t+1}(Z^n_{\leq t+1}) \mid Z^n_{\leq s} = z^n_{\leq s}) \stackrel{w}{\to} \mathcal{L}(\psi_{t+1}(Z_{\leq t+1}) \mid Z_{\leq s} = z_{\leq s}).
\end{align}
holds is identical to that in the proofs of Theorem \ref{thm:convergentvaluations} assuming \eqref{eq:phialtmapping} and \eqref{eq:phinaltmapping}.
By Lemma \ref{lem:BoundUIUtility}, there exists a constant $c > 0$, such that
$$
|\psi^n_{t}(z)| \leq c \operatorname{E}[|\psi^n_{t+1}(Z^n_{\leq t+1})| \mid Z^n_{\leq t} = z].
$$
Hence, to verify {\bf UI($t$)} it is sufficient to show uniform integrability of
$$
\Big(\mathcal{L}\Big(\operatorname{E}[|\psi^n_{t+1}(Z^n_{\leq t+1})| \mid Z^n_{\leq t}] \mid Z^n_{\leq s} = z^n_{\leq s} \Big)\Big)_{n\in\mathbb{N}}
$$
and after setting $W^n = |\psi^n_{t+1}(Z^n_{\leq t+1})|$, $\widetilde{W}^n = \operatorname{E}[|\psi^n_{t+1}(Z_{\leq t+1}^n)| \mid Z^n_{\leq t}]$,
 {\bf UI($t$)} follows from {\bf UI($t+1$)} and \eqref{conv:Weak123456} by applying Lemma \ref{lem:TowerUI}.

We proceed by showing {\bf CC($t$)}, i.e.~the convergence of \eqref{eq:ConvUtility3}.
We first consider \eqref{eq:ConvRisk3}.
If $\mu_t$ admits a bounded density, convergence $\rho_t^n(z^n_{\leq t}) \to \rho_t(z_{\leq t})$ can be shown by Pratt's Lemma, together with \eqref{conv:Weak123456} and {\bf UI($t+1$)} (cf. Proof of Theorem \ref{thm:convergentvaluations} assuming \eqref{eq:phialtmapping} and \eqref{eq:phinaltmapping}). If $\mu_t$ satisfies $\operatorname{supp}(\mu_t) \subset [a, b]$ with $0<a<b<1$, we first note that we have pointwise convergence $F^{-1}_{-\psi^n_{t+1}(Z^n_{\leq t})\mid Z^n_{\leq t} = z^n_{\leq t}}(p) \to F^{-1}_{-\psi_{t+1}(Z_{\leq t})\mid Z_{\leq t} = z_{\leq t} }(p) $ for all $p \in (0,1)$ since the limiting distribution in \eqref{conv:Weak123456} is Gaussian (cf. Lemma \ref{lem:gaussian_value}). There exists $c > 0$ such that
\begin{align*}
\sup_n\sup_{p \in [a,b]}|F^{-1}_{-\psi^n_{t+1}(Z^n_{\leq t+1}) \mid Z^n_{\leq t} = z^n_{\leq t}}(p)|
 \leq \sup_n c \operatorname{E}[|\psi^n_{t+1}(Z^n_{\leq t+1}) | \mid Z^n_{\leq t} = z^n_{\leq t}]
\end{align*}
and the right-hand side is finite.
The inequality is due to properties of quantile functions (cf. Proof of Lemma \ref{lem:BoundRec}). The finiteness is because of the convergence of the conditional expectation which is due to \eqref{conv:Weak123456} and {\bf UI($t+1$)}. Hence, dominated convergence yields $\nu^n_t(z^n_{\leq t}) \to \nu_t(z_{\leq t})$.

We now show convergence of the conditional expectation in \eqref{eq:ConvUtility3} which can be written in terms of
\begin{align*
\operatorname{E} \Big[h_n\big(\psi^n_{t+1}(Z^n_{\leq t+1})\big) \mid Z^n_{\leq t} = z^n_{\leq t}\Big],
\end{align*}
where $h_n(w) = \big(\big(\nu^n_t(z^n_{\leq t}) - w\big)^{+}\big)^{\beta_t}$ converges continuously towards $h(w) = \big(\big(\nu_t(z_{\leq t}) - w\big)^{+}\big)^{\beta_t}$. Moreover, since
$$
|h_n(w)| \leq \bigg(\sup_n(\nu^n_t(z^n_{\leq t})) + |w|\bigg)^{\beta_t}, \quad \sup_n(\nu^n_t(z^n_{\leq t})) < \infty,
$$
and $\beta_t\in (0,1]$, uniform integrability of $\mathcal{L}\big(h_n\big(\psi^n_{t+1}(Z^n_{\leq t+1})\big) \mid Z^n_{\leq t} = z^n_{\leq t}\big)$ follows from {\bf UI($t+1$)}. This means convergence
\begin{align*
\operatorname{E}\Big[h_n\big(\psi^n_{t+1}(Z^n_{\leq t+1})\big) \mid Z^n_{\leq t} = z^n_{\leq t}\Big]
\to \operatorname{E}\Big[h\big(\psi_{t+1}(Z_{\leq t+1})\big) \mid Z_{\leq t} = z_{\leq t}\Big]
\end{align*}
which implies {\bf CC($t$)} and the proof of the induction step is complete.

Finally, the proof of the convergence $\lim_{n\to\infty}\psi^n_0=\psi_0$ follows by arguments analogous to those in the proof of Theorem \ref{thm:convergentvaluations} assuming \eqref{eq:phialtmapping} and \eqref{eq:phinaltmapping}.
\end{proof}

\begin{lemma}\label{lem:quant_int_bound}
Let $W$ be a random variable with quantile function $F^{-1}(p)=\min\{x\in\mathbb{R}:F_{W}(x)\geq p\}$. If $\mu\in\mathcal{P}([0,1])'$, then there exists $c\in (0,\infty)$ such that
\begin{align*}
\Big |\int F^{-1}_{W}(p)\mu(dp) \Big | \leq c\operatorname{E}[|W|].
\end{align*}
\end{lemma}

\begin{proof}
We have the general bounds $l,u$ given by
\begin{align*}
l=\int F^{-1}_{-|W|}(p)\mu(dp) \leq \int F^{-1}_{W}(p)\mu(dp) \leq \int F^{-1}_{|W|}(p)\mu(dp)=u.
\end{align*}
Moreover,
\begin{align*}
F^{-1}_{-|W|}(p)=-F^{-1}_{|W|}((1-p)+)=-\lim_{v\downarrow 1-p}F^{-1}_{|W|}(v).
\end{align*}
For any $p\in (0,1)$,
\begin{align*}
pF^{-1}_{|W|}(1-p)&\leq \int_{1-p}^{1}F^{-1}_{|W|}(v)dv\leq \operatorname{E}[|W|], \\
pF^{-1}_{-|W|}(p)&\geq \int_{0}^{p}F^{-1}_{-|W|}(v)dv\geq -\operatorname{E}[|W|].
\end{align*}
Consequently, if $\operatorname{supp}(\mu)\subset [a,b]\subset (0,1)$, then
\begin{align*}
l\geq -\frac{1}{a}\operatorname{E}[|W|], \quad u\leq \frac{1}{1-b}\operatorname{E}[|W|], \quad \Big | \int F^{-1}_{W}(p)\mu(dp)\Big | \leq \max\Big(\frac{1}{a},\frac{1}{1-b}\Big)\operatorname{E}[|W|].
\end{align*}
If $\mu$ has a density bounded by $c'\in (0,\infty)$, then
\begin{align*}
l\geq -c'\operatorname{E}[|W|], \quad u\leq c'\operatorname{E}[|W|], \quad \Big | \int F^{-1}_{W}(p)\mu(dp)\Big | \leq c'\operatorname{E}[|W|].
\end{align*}
\end{proof}

\begin{lemma}\label{lem:BoundRec}
If $\psi^n_t$ is given by \eqref{eq:cvfun_ta} and \eqref{eq:cvfun_0a}, then exists $c\in (0,\infty)$ such that
\begin{align*}
|\psi^n_{0}| \leq c \operatorname{E}[|\psi_{1}^n(Z^n_{1})|] \quad \text{and} \quad
|\psi^n_{t}(z)| \leq c \operatorname{E}[|\psi_{t+1}^n(Z^n_{\leq t+1})| \mid Z^n_{\leq t} = z]
\end{align*}
for any $z\in\operatorname{supp}(\mathcal{L}(Z^{n}_{\leq t}))$.
\end{lemma}

\begin{proof}
We proof the inequality for $t\geq 1$. The argument for $t=0$ is identical.
Since
\begin{align*}
|\psi^n_t(z)|
\leq \Big |\int F^{-1}_{\psi^n_{t+1}(Z^n_{\leq t+1}) \mid Z^n_{\leq t} = z}(p) \mu_t^1(dp) \Big |
+ \Big | \int F^{-1}_{-\psi^n_{t+1}(Z^n_{\leq t+1}) \mid Z^n_{\leq t} = z}(p) \mu_t^2(dp) \Big |,
\end{align*}
we can apply Lemma \ref{lem:quant_int_bound}, bounding the right-hand side by $c\operatorname{E}[|\psi^n_{t+1}(Z^n_{\leq t})| \mid Z^n_{\leq t} = z]$ for some $c$.
\end{proof}

The following result is similar to Proposition 6 in \cite{engsner2017}.

\begin{lemma}\label{lem:gaussian_value}
Let $(Z_t)_{t=1}^T=(X_t,Y_t)_{t=1}^T$ be a Gaussian process, let $\mathcal{F}_0=\{\Omega,\emptyset\}$, and let $\mathcal{F}_t=\sigma(Z_{\leq t})$ for $t\geq 1$.
Let $V_T=0$ and, for $t\in\{0,1,\dots,T-1\}$, let $V_t=\phi_t(X_{t+1}+V_{t+1})$, where $\phi_t$ is positive homogeneous and conditionally cash additive.
Suppose that for each $t$,
\begin{align*}
\text{if}\quad \operatorname{P}(W \in \cdot \mid \mathcal{F}_t)=\operatorname{P}(\widetilde{W}\in \cdot \mid \mathcal{F}_t), \quad\text{then}\quad
\phi_t(W)=\phi_t(\widetilde{W}),
\end{align*}
and if $W$ is independent of $\mathcal{F}_t$, then $\phi_t(W)$ is nonrandom.
Then
\begin{align*}
V_t=\operatorname{E}\bigg[\sum_{u=t+1}^T X_{u}\mid \mathcal{F}_t\bigg]+\sum_{u=t+1}^T \phi_{u-1}(\varepsilon_u)\operatorname{Var}\bigg(\operatorname{E}\bigg[\sum_{v=u}^T X_v\mid \mathcal{F}_u\bigg]\mid\mathcal{F}_{u-1}\bigg)^{1/2},
\end{align*}
where $\varepsilon_u$ is standard normal and independent of $\mathcal{F}_{u-1}$.
\end{lemma}

\begin{proof}
The lemma is proved by backward induction. We prove the induction step. Suppose that $X_{t+1}+V_{t+1}$ and $Z_{\leq t}$ are jointly Gaussian. Then there exists a standard normal $\varepsilon_{t+1}$ independent of $\mathcal{F}_t$ such that
\begin{align*}
\operatorname{P}(X_{t+1}+V_{t+1}\in \cdot \mid \mathcal{F}_t)=\operatorname{P}\Big(\operatorname{E}[X_{t+1}+V_{t+1}\mid \mathcal{F}_t]+\varepsilon_{t+1}\operatorname{Var}(X_{t+1}+V_{t+1}\mid \mathcal{F}_t)^{1/2}\in \cdot \mid\mathcal{F}_t\Big)
\end{align*}
and it follows from properties of the multivariate normal distribution that the conditional variance $\operatorname{Var}(X_{t+1}+V_{t+1}\mid \mathcal{F}_t)$ is nonrandom. Hence, by positive homogeneity and conditional cash additivity,
\begin{align*}
\phi_t(X_{t+1}+V_{t+1})=\operatorname{E}[X_{t+1}+V_{t+1}\mid \mathcal{F}_t]+\phi_t(\varepsilon_{t+1})\operatorname{Var}(X_{t+1}+V_{t+1}\mid \mathcal{F}_t)^{1/2}.
\end{align*}
By the induction assumption (in the nontrivial case $t+1<T$),
\begin{align*}
V_{t+1}=\operatorname{E}\bigg[\sum_{u=t+2}^T X_{u}\mid \mathcal{F}_{t+1}\bigg]+\sum_{u=t+2}^T \phi_{u-1}(\varepsilon_u)\operatorname{Var}\bigg(\operatorname{E}\bigg[\sum_{v=u}^T X_v\mid \mathcal{F}_u\bigg]\mid\mathcal{F}_{u-1}\bigg)^{1/2}.
\end{align*}
Summing up the terms and using the tower property of conditional expectations complete the proof.
\end{proof}

\begin{lemma}\label{lem:TowerUI}
Let $s \in \{1, \dots, t-1\}$ and let $(W^n)_{n\in\mathbb{N}}$ be a sequence of non-negative random variables such that
\begin{align}\label{eq:wczTUI}
\mathcal{L}(W^n \mid Z^n_{\leq s} = z^n_{\leq s}) \stackrel{w}{\to} \mathcal{L}(W \mid Z_{\leq s} = z_{\leq s}) \quad\text{as } n\to\infty
\end{align}
and
\begin{align}\label{eq:uizTUI}
\Big(\mathcal{L}(W^n \mid Z^n_{\leq s} = z^n_{\leq s})\Big)_{n\in\mathbb{N}} \quad\text{is uniformly integrable.}
\end{align}
If $\widetilde{W}^n = \operatorname{E}[W^n \mid Z^n_{\leq t}]$, then $\big(\mathcal{L}(\widetilde{W}^n \mid Z^n_{\leq s} = z^n_{\leq s}) \big)_{n\in\mathbb{N}}$ is uniformly integrable.
\end{lemma}

\begin{proof}
First note that
\begin{align*}
\operatorname{E}\Bigl[\widetilde{W}^nI\{\widetilde{W}^n \geq r\}\mid Z^n_{\leq s} \Bigr]
& =\operatorname{E}\Bigl[\operatorname{E}[W^n \mid Z^n_{\leq t}]I\{\widetilde{W}^n \geq r\}\mid Z^n_{\leq s} \Bigr] \\
& =\operatorname{E}\Bigl[\operatorname{E}[W^n I\{\widetilde{W}^n \geq r\} \mid Z^n_{\leq t}]\mid Z^n_{\leq s}\Bigr]\\
& =\operatorname{E}\Bigl[W^n I\{\widetilde{W}^n \geq r\} \mid Z^n_{\leq s}\Bigr],
\end{align*}
i.e. $\operatorname{E}[\widetilde{W}^{n}I\{\widetilde{W}^{n} \geq r\}\mid Z^n_{\leq s}=z^n_{\leq s}]=\operatorname{E}[W^{n} I\{\widetilde{W}^{n} \geq r\} \mid Z^n_{\leq s}=z^n_{\leq s}]$.
For each $\gamma_r > 0$,
\begin{align}\label{bound:23}
&\operatorname{E}[W^{n}I\{\widetilde{W}^{n} \geq r\}\mid Z^n_{\leq s}=z^n_{\leq s}]  \nonumber \\
&\quad= \operatorname{E}[W^{n}I\{\widetilde{W}^{n} \geq r\} I\{W^n \leq \gamma_r\}\mid Z^n_{\leq s}=z^n_{\leq s}] \nonumber \\
&\quad\quad + \operatorname{E}[W^{n}I\{\widetilde{W}^{n} \geq r\} I\{W^n > \gamma_r\}\mid Z^n_{\leq s}=z^n_{\leq s}] \nonumber \\
&\quad \leq \gamma_{r}\operatorname{P}(\widetilde{W}^{n} \geq r \mid Z^n_{\leq s}=z^n_{\leq s}) +\operatorname{E}[W^{n}I\{W^{n}>\gamma_{r}\}\mid Z^n_{\leq s}=z^n_{\leq s}].
\end{align}
By \eqref{eq:wczTUI} and \eqref{eq:uizTUI},
$$
\operatorname{E}[\widetilde{W}^{n}I\{\widetilde{W}^{n} \geq r\}\mid Z^n_{\leq s}=z^n_{\leq s}] < \infty \quad\text{for all } n\in\mathbb{N},
$$
and similarly with $W^{n}$ instead of $\widetilde{W}^{n}$.
Therefore, uniform integrability of $\big(\mathcal{L}(\widetilde{W}^n \mid Z^n_{\leq s} = z^n_{\leq s}) \big)_{n\in\mathbb{N}}$ is equivalent to (cf.~\cite{Kallenberg-02} p.~67)
\begin{align*}
\lim_{r\to\infty}\limsup_{n\to\infty}\operatorname{E}[\widetilde{W}^{n}I\{\widetilde{W}^n\geq r\}\mid Z^n_{\leq s}=z^n_{\leq s}]=0
\end{align*}
which follows if we show that there exists $\gamma_{r}$ such that $\lim_{r \to \infty}\limsup_{n\to\infty}[\eqref{bound:23}]=0$.
By the uniform integrability assumption \eqref{eq:uizTUI} and the finiteness of all the terms,
\begin{align*}
\lim_{r\to\infty}\limsup_{n\to\infty}\operatorname{E}[W^{n}I\{W^{n}>\gamma_{r}\}\mid Z^n_{\leq s}=z^n_{\leq s}]=0.
\end{align*}
Therefore, it only remains to show that
$$
\lim_{r \to \infty}\limsup_{n\to\infty}\gamma_{r}\operatorname{P}(\widetilde{W}^n \geq r \mid Z^n_{\leq s} = z^n_{\leq s})=0
$$
for some sequence $\gamma_{r}\to\infty$. By Markov's inequality,
\begin{align*}
\operatorname{P}(\widetilde{W}^n \geq r \mid Z^n_{\leq s} = z^n_{\leq s})
&\leq \frac{ \operatorname{E}[\widetilde{W}^n \mid Z^n_{\leq s} = z^n_{\leq s}]}{r} \\
&= \frac{\operatorname{E}[W^n \mid Z^n_{\leq s} = z^n_{\leq s}]}{r} \\
&\to \frac{\operatorname{E}[W \mid Z_{\leq s} = z_{\leq s}]}{r},
\end{align*}
where the convergence in the last step follows from Lemma 4.11 in \cite{Kallenberg-02}, using \eqref{eq:wczTUI} and \eqref{eq:uizTUI}.   Letting e.g.~$\gamma_r = \sqrt{r}$, the first term in \eqref{bound:23} can therefore be made arbitrarily small by choosing $r$ sufficiently large.
\end{proof}

\begin{lemma}\label{lem:BoundUIUtility}
If $\nu^n_t,\psi^n_t$ are given by \eqref{eq:ConvRisk3}, \eqref{eq:ConvUtility3} and \eqref{eq:ConvRU3}, then exists $c\in (0,\infty)$ such that
\begin{align*}
|\psi^n_{0}| \leq c \operatorname{E}[|\psi_{1}^n(Z^n_{1})|] \quad \text{and} \quad
|\psi^n_{t}(z)| \leq c \operatorname{E}[|\psi_{t+1}^n(Z^n_{\leq t+1})| \mid Z^n_{\leq t} = z]
\end{align*}
for any $z\in\operatorname{supp}(\mathcal{L}(Z^{n}_{\leq t}))$.
\end{lemma}

\begin{proof}[Proof of Lemma \ref{lem:BoundUIUtility}]
Consider $t\geq 1$.
Note that, by Lemma \ref{lem:quant_int_bound}, there exists $c' > 0$, such that
$$
|\nu^n_t(z)| \leq c' \operatorname{E}[|\psi^n_{t+1}(Z^n_{\leq t+1})| \mid Z^n_{\leq t} = z].
$$
Hence,
\begin{align*}
|\psi^n_t(z)| &\leq c' \operatorname{E}[ |\psi^n_{t+1}(Z^n_{\leq t+1})| \mid Z^n_{\leq t} = z] \\
& \quad +\gamma_t\operatorname{E} \Big[ \Big( \big(c' \operatorname{E}[ |\psi^n_{t+1}(Z^n_{\leq t+1})| \mid Z^n_{\leq t} = z]-\psi^n_{t+1}(Z^n_{\leq t+1})\big)^+ \Big)^{\beta_t} \mid Z^n_{\leq t} = z\Big]^{1/\beta_t}.
\end{align*}
By Jensen's inequality, the second summand above is bounded by
\begin{align*}
&\gamma_t\operatorname{E} \Big[ \big(c' \operatorname{E}[ |\psi^n_{t+1}(Z^n_{\leq t+1})| \mid Z^n_{\leq t} = z]- \psi^n_{t+1}(Z^n_{\leq t+1}) \big)^+\mid Z^n_{\leq t} = z  \Big] \\
&\quad \leq \gamma_t\operatorname{E} \Big[ c' \operatorname{E}[ |\psi^n_{t+1}(Z^n_{\leq t+1})| \mid Z^n_{\leq t} = z] + |\psi^n_{t+1}(Z^n_{\leq t+1})| \mid Z^n_{\leq t} = z\Big] \\
& \quad = (c' + 1) \gamma_t \operatorname{E}[ |\psi^n_{t+1}(Z^n_{\leq t+1})| \mid Z^n_{\leq t} = z].
\end{align*}
Hence, the statement holds with $c=c'+\gamma_t(c'+1)$.
The inequality
$|\psi^n_{0}| \leq c \operatorname{E}[|\psi_{1}^n(Z^n_{1})|]$ is shown by the same arguments.
\end{proof}

\begin{proof}[Proof of Theorem \ref{thm:CondExpVar}]
The proof proceeds in a similar manner as the proof of Theorem \ref{thm:convergentvaluations}. For $t=1,\dots,T$, consider mappings $\psi^n_t:\operatorname{supp}(\mathcal{L}(Z^n_{\leq t}))\to \mathbb{R}$ given by
\begin{align*}
\psi^n_t(z^n_{\leq t}) =
 \left\{\begin{array}{ll}
\sum_{s=1}^T x^n_s, & t=T,\\
 \operatorname{E} \big[\psi^n_{t+1}(Z^n_{\leq t+1}) \mid Z^n_{\leq t} = z^n_{\leq t}\big] + c_t \operatorname{Var}\big(\psi^n_{t+1}(Z^n_{\leq t+1}) \mid Z^n_{\leq t} = z^n_{\leq t}\big)^{1/2}, & t<T,
\end{array}\right.
\end{align*}
and
\begin{align*}
\psi^n_0 =  \operatorname{E} \big[\psi^n_{1}(Z^n_{1}) \big] + c_0 \operatorname{Var}\big(\psi^n_{1}(Z^n_{1})\big)^{1/2}.
\end{align*}
We define $\psi_t:\mathbb{R}^t\to \mathbb{R}$ and $\psi_0$ analogously without the superscript $n$.
The proof is complete once we show the convergence $\lim_{n\to\infty}\psi^n_0=\psi_0$.

The argument of the proof is backwards induction. We will show that for $t = T, \dots, 1$,
\begin{itemize}
\item[{\bf UI($t$):}] $(\mathcal{L}(\psi^n_t(Z^n_{\leq t})^2 \mid Z^n_{\leq s} = z^n_{\leq s}))_n$ is uniformly integrable for each $s = 1, \dots, t-1$ and $z^n_{\leq s} \to z_{\leq s}$,
\item[{\bf CC($t$):}] $(\psi^n_t)_n$ is continuously convergent, i.e. $\psi^n_t(z^n_{\leq t}) \to \psi_t(z_{\leq t})$, whenever $z^n_{\leq t} \to z_{\leq t}$.
\end{itemize}

\noindent
{\bf Induction base:} {\bf UI($T$)} and {\bf CC($T$)} hold.

For $s \in \{1, \dots, T-1\}$ and $z^n_{\leq s} \to z_{\leq s}$, {\bf UI($T$)} follows immediately from \eqref{eq:condxnuiII} since $\psi^n_{T}(z^n_{\leq T}) = \sum_{s=1}^T x^n_s$.
Moreover,
$$
\psi^n_{T}(z^n_{\leq T}) = \sum_{s=1}^T x^n_s \to \sum_{s=1}^T x_s = \psi_T(z_{\leq T})
$$
whenever $z_{\leq T}^n \to z_{\leq T}$, yielding {\bf CC($T$)}.

\noindent
{\bf Induction step:}{\bf UI($t+1$)}, {\bf CC($t+1$)} together imply {\bf UI($t$)}, {\bf CC($t$)}.
Fix $s\in \{1,\dots,t-1\}$ and a sequence $(z^n_{\leq s})_n$ with $z^n_{\leq s}\to z_{\leq s}$.
By the same argument as in the proof of Theorem \ref{thm:convergentvaluations},
\begin{align}\label{conv:Weak1234}
\mathcal{L}(\psi^n_{t+1}(Z^n_{\leq t+1}) \mid Z^n_{\leq s} = z^n_{\leq s}) \stackrel{w}{\to} \mathcal{L}(\psi_{t+1}(Z_{\leq t+1}) \mid Z_{\leq s} = z_{\leq s})
\end{align}
which implies
\begin{align}\label{conv:Weak12345}
\mathcal{L}((\psi^n_{t+1}(Z^n_{\leq t+1}))^2 \mid Z^n_{\leq s} = z^n_{\leq s}) \stackrel{w}{\to} \mathcal{L}((\psi_{t+1}(Z_{\leq t+1}))^2 \mid Z_{\leq s} = z_{\leq s}).
\end{align}
Lemma \ref{lem:gaussian_value} shows that the limiting distribution in \eqref{conv:Weak1234} is Gaussian. We can bound
\begin{align*}
\psi^n_t(z)^2 &\leq 2\operatorname{E}[|\psi^n_{t+1}(Z^n_{\leq t+1})| \mid Z^n_{\leq t} = z]^2 + 2c_t^2\operatorname{E}[|\psi^n_{t+1}(Z^n_{\leq t+1})|^2 \mid Z^n_{\leq t} = z] \\
&\leq (2+2c_t^2)\operatorname{E}[|\psi^n_{t+1}(Z^n_{\leq t+1})|^2 \mid Z^n_{\leq t} = z].
\end{align*}
Therefore,
\begin{align*}
\operatorname{E}[|\psi^n_t(Z^n_{\leq t})|^2 \, \mid Z^n_{\leq s} = z^n_{\leq s}]
&\leq (2+2c_t^2)\operatorname{E}[\operatorname{E}[|\psi^n_{t+1}(Z^n_{\leq t+1})|^2 \, \mid Z^n_{\leq t}]\mid Z^n_{\leq s} = z^n_{\leq s}]\\
&=(2+2c_t^2)\operatorname{E}[|\psi^n_{t+1}(Z^n_{\leq t+1})|^2 \, \mid Z^n_{\leq s} = z^n_{\leq s}].
\end{align*}
Hence, it is sufficient to show uniform integrability of
\begin{align}\label{eq:UI234}
\Big(\mathcal{L}\big(\operatorname{E}[|\psi^n_{t+1}(Z^n_{\leq t+1})|^2 \mid Z^n_{\leq t}] \mid Z^n_{\leq s}=z^n_{\leq s}\big)\Big)_n.
\end{align}
Setting $\widetilde{W}^{n}=\operatorname{E}[|\psi^n_{t+1}(Z^n_{\leq t+1})|^2 \mid Z^n_{\leq t}]$ and $W^{n}=|\psi^n_{t+1}(Z^n_{\leq t+1})|^2$, this follows from Lemma \ref{lem:TowerUI} together with \eqref{conv:Weak12345} and {\bf UI($t+1$)}. Hence, {\bf UI($t$)} holds.

We proceed by showing {\bf CC($t$)}. Let $z^n_{\leq t} \to z_{\leq t}$. We need to show convergence of
\begin{align*}
\psi^n_t(z^n_{\leq t})
&= \operatorname{E}[\psi^n_{t+1}(Z^n_{\leq t+1}) \mid Z^n_{\leq t} = z^n_{\leq t}] \\
& \quad + c_t \left(\operatorname{E}[(\psi^n_{t+1}(Z^n_{\leq t+1}))^2 \mid Z^n_{\leq t} = z^n_{\leq t}] - \operatorname{E}[\psi^n_{t+1}(Z^n_{\leq t+1}) \mid Z^n_{\leq t} = z^n_{\leq t}]^2 \right)^{1/2}.
\end{align*}
Indeed, by \eqref{conv:Weak1234} and \eqref{conv:Weak12345},
\begin{align*}
&\mathcal{L}(\psi^n_{t+1}(Z^n_{\leq t+1}) \mid Z^n_{\leq t}  = z^n_{\leq t}) \stackrel{w}{\to} \mathcal{L}(\psi_{t+1}(Z_{\leq t+1}) \mid Z_{\leq t} = z_{\leq t}), \\
&\mathcal{L}((\psi^n_{t+1}(Z^n_{\leq t+1}))^2 \mid Z^n_{\leq t}  = z^n_{\leq t}) \stackrel{w}{\to} \mathcal{L}((\psi_{t+1}(Z_{\leq t+1}))^2 \mid Z_{\leq t} = z_{\leq t}),
\end{align*}
and by \cite[Lemma 4.11]{Kallenberg-02} applied to {\bf UI($t+1$)}, it follows that
$$
\operatorname{E}[\psi^n_{t+1}(Z^n_{\leq t+1}) \mid Z^n_{\leq t} = z^n_{\leq t}] \to \operatorname{E}[\psi_{t+1}(Z_{\leq t+1}) \mid Z_{\leq t} = z_{\leq t}]
$$
and
$$
\operatorname{E}[(\psi^n_{t+1}(Z^n_{\leq t+1}))^2 \mid Z^n_{\leq t} = z^n_{\leq t}] \to \operatorname{E}[(\psi_{t+1}(Z_{\leq t+1}))^2 \mid Z_{\leq t} = z_{\leq t}],
$$
yielding {\bf CC($t$)} which completes the proof of the induction step.

As a result of the induction,
$$
\mathcal{L}(\psi^n_1(Z^n_1)) \stackrel{w}{\to} \mathcal{L}(\psi_1(Z_1))
$$
and $(\psi^n_1(Z^n_1))_n, ((\psi^n_1(Z^n_1))^2)_n$ are uniformly integrable. Therefore,
\begin{align*}
\psi^n_0 &= \operatorname{E}[\psi^n_{1}(Z^n_1)] + c_0\left(\operatorname{E}[(\psi^n_{1}(Z^n_1))^2] - \operatorname{E}[\psi^n_1(Z^n_1)]^2 \right)^{1/2} \\
 &\to \operatorname{E}[\psi_{1}(Z_1)] + c_0\left(\operatorname{E}[(\psi_{1}(Z_1))^2] - \operatorname{E}[\psi_1(Z_1)]^2 \right)^{1/2} = \psi_0.
\end{align*}
\end{proof}

\begin{lemma}\label{lem:V0monotonicity}
Let $c\in \mathbb{R}^T_+$ satisfy $\sum_{t=1}^Tc_t=1$ and $c_{t+1}\leq c_t$ for every $t$. Then $c$ solves
\begin{align}\label{optproblem}
\max_{d\in \mathcal{D}_{c,T}} \sum_{t=1}^T d_t^{1/2},
\quad \mathcal{D}_{c,T}=\bigg\{d\in\mathbb{R}^T_+:\sum_{t=1}^T d_t=1,\sum_{s=1}^t d_s \geq \sum_{s=1}^t c_s \text{ for } t=1,\dots,T\bigg\}.
\end{align}
\end{lemma}

\begin{proof}
Since $\mathcal{D}_{c,T}$ is a convex set and the objective function is concave, we have a convex optimization problem. Consequently, it is sufficient to verify that $c\in\mathcal{D}_{c,T}$ is a local optimum.
By assumption either $c_1=1$ and $c_t=0$ for $t\geq 2$ or there exists $t_0\in\{2,\dots,T\}$ such that $c_1\geq \dots \geq c_{t_0}>0$ and $c_{t_0+1}=\dots=c_{T}=0$. The first case has the trivial maximizer $d_1=1$ and $d_t=0$ for $t\geq 2$. Hence, it is sufficient to consider only the second case for a fixed $c$ and $t_0\in\{2,\dots,T\}$. In this case, \eqref{optproblem} is equivalent to, with $c=(c_1,\dots,c_{t_0})$,
\begin{align}\label{modoptproblem}
\max_{d\in \mathcal{D}_c} \sum_{t=1}^{t_0} d_t^{1/2},
\quad \mathcal{D}_c=\bigg\{d\in\mathbb{R}^{t_0}_+:\sum_{t=1}^{t_0} d_t=1,\sum_{s=1}^t d_s \geq \sum_{s=1}^t c_s \text{ for } t=1,\dots,t_0\bigg\}.
\end{align}
Since $\mathcal{D}_c$ is a convex set and the objective function is concave, we have a convex optimization problem. Consequently, it is sufficient to verify that $c\in\mathcal{D}_c$ is a local optimum.
Let
\begin{align*}
\mathcal{X}=\bigg\{x\in\mathbb{R}^{t_0}:\sum_{t=1}^{t_0} x_t=0,\sum_{s=1}^t x_s \geq 0 \text{ for } t=1,\dots,t_0\bigg\}
\end{align*}
and note that $\mathcal{D}_c\subset c+\mathcal{X}$.
Consider the set of vectors $\{b_1,\dots,b_{t_0-1}\}\subset \mathbb{R}^{t_0}$ given by
\begin{align*}
b_{i,i}=1, \; b_{i,i+1}=-1, \; b_{i,j}=0 \text{ for } j\notin \{i,i+1\}.
\end{align*}
We claim that $\mathcal{X}=\textrm{span}_+\{b_1,\dots,b_{t_0-1}\}$, where
\begin{align*}
\textrm{span}_+\{b_1,\dots,b_{t_0-1}\}=\bigg\{\sum_{k=1}^{t_0-1}\lambda_k b_k : \lambda_1,\dots,\lambda_{t_0-1}\geq 0\bigg\}.
\end{align*}
To show that $\textrm{span}_+\{b_1,\dots,b_{t_0-1}\}\subseteq \mathcal{X}$ it is sufficient to note that
\begin{align*}
&\sum_{s=1}^{t_0}\sum_{k=1}^{t_0-1}\lambda_k b_{k,s}=\sum_{k=1}^{t_0-1}\lambda_k\sum_{s=1}^{t_0}b_{k,s}=0, \\
&\sum_{s=1}^t\sum_{k=1}^{t_0-1}\lambda_k b_{k,s}=\sum_{k=1}^{t_0-1}\lambda_k\sum_{s=1}^tb_{k,s}=\lambda_t\geq 0.
\end{align*}
We now show that $\mathcal{X}\subseteq\textrm{span}_+\{b_1,\dots,b_{t_0-1}\}$. Take $x\in\mathcal{X}$ and set $\lambda_1=x_1$ (noting that $x_1\geq 0$) and $\lambda_k=\lambda_{k-1}+x_k$ for $k\geq 2$. Hence, $\lambda_k\geq 0$ for every $k$ and $\sum_{k=1}^{t_0-1}\lambda_kb_k=x$.

Let $g(d)=\sum_{t=1}^{t_0} d_t^{1/2}$ and note that $g$ is well defined and concave on $\mathcal{D}_c$, and has a well defined gradient at $c$:
\begin{align*}
\nabla g(c)^{{\mathrm{T}}}=\frac{1}{2}(c_1^{-1/2},\dots,c_{t_0}^{-1/2}).
\end{align*}
Take $d\in\mathcal{D}_c$ and notice that $d=c+x$ for $x\in\mathcal{X}$. Since $g$ is concave,
\begin{align*}
g(d)-g(c)\leq \nabla g(c)^{{\mathrm{T}}}x
=\sum_{k=1}^{t_0-1}\lambda_k\nabla g(c)^{{\mathrm{T}}}b_k
=\frac{1}{2}\sum_{k=1}^{t_0-1}\lambda_k(c_k^{-1/2}-c_{k+1}^{-1/2})\leq 0,
\end{align*}
where the last inequality holds since, for every $k$, $\lambda_k\geq 0$ and $c_k^{-1/2}-c_{k+1}^{-1/2}\leq 0$.
Hence, we have shown that $g(d)\leq g(c)$ and the proof is complete.
\end{proof}

\begin{proof}[Proof of Theorem \ref{thm:V0monotonicity}]
For $t=1,\dots,T$, set
\begin{align*}
c_t&=\operatorname{Var}\bigg(\sum_{s=1}^TX_s\bigg)^{-1}\bigg(\operatorname{Var}\bigg(\sum_{s=1}^TX_s\mid \mathcal{F}_{t-1}\bigg)-\operatorname{Var}\bigg(\sum_{s=1}^TX_s\mid \mathcal{F}_t\bigg)\bigg),\\
d_t&=\operatorname{Var}\bigg(\sum_{s=1}^TX_s\bigg)^{-1}\bigg(\operatorname{Var}\bigg(\sum_{s=1}^TX_s\mid \mathcal{G}_{t-1}\bigg)-\operatorname{Var}\bigg(\sum_{s=1}^TX_s\mid \mathcal{G}_t\bigg)\bigg).
\end{align*}
By construction, $\sum_{t=1}^Tc_t=\sum_{t=1}^Td_t=1$. Since $\mathcal{F}_t\subseteq \mathcal{G}_t$ for every $t$, $\sum_{s=t+1}^Tc_s\geq \sum_{s=t+1}^Td_s$ which is equivalent to $\sum_{s=1}^t d_s \geq \sum_{s=1}^t c_s$ for every $t$ since $\sum_{t=1}^Tc_t=\sum_{t=1}^Td_t$. Applying Lemma \ref{lem:V0monotonicity} completes the proof.
\end{proof}

\section*{Acknowledgements}
F.~Lindskog would like to acknowledge financial support from the Swedish Research Council, Project 2020-05065, and from L\"ansf\"ors\"akringars Forskningsfond, Project P9.20.


\section{Introduction}

The only candidate for the effective theory of high energy QCD is Colour Glass Condensate(CGC) approach
(see  Ref.\cite{KOLEB} for a review). Two main ideas of CGC:
the saturation of the dipole density and the new dimensional scale ($Q_s$), which increases with energy, have become  widely accepted language for discussing the high energy scattering in QCD. However the CGC approach suffers several problems. The most known of them is the power-like behaviour of the scattering amplitude at large impact parameters\cite{KW1,KW2,KW3,FIIM} that violated the Froissart theorem\cite{FROI}. We have to introduce non-perturbative corrections at large impact parameters and the embryonic stage of our understanding of confinement of quarks and gluons does not allow us to come up with a reasonable theoretical approach to  the problem.  The second well known problem is summing of the BFKL Pomeron loops\footnote{ The abbreviation BFKL Pomeron  stands  for Balitsky, Fadin, Kuraev and Lipatov Pomeron.}. This problem is a technical one, but in spite of intensive work
 \cite{BFKL,LIP,LIREV,LIFT,GLR,GLR1,MUQI,MUDI,Salam,NAPE,BART,BKP,MV,MUSA,KOLE,BRN,BRAUN,BK,KOLU,JIMWLK1,JIMWLK2,JIMWLK3,JIMWLK4,JIMWLK5,JIMWLK6,JIMWLK7,JIMWLK8,AKLL,KOLU1,KOLUD,BA05,SMITH,KLW,KLLL1,KLLL2,LEN}
 it is still far away from being solved.This situation makes the problem  one of the  principle problems,  without solving which we cannot consider the dilute-dilute and dense-dense parton densities collisions. As has been recently shown\cite{KLLL1,KLLL2}, even the Balitsky-Kovchegov (BK) equation, that governs the dilute-dense parton density  scattering (deep inelastic scattering (DIS) of electron with proton), has to be modified due to  contributions of  Pomeron loops.

  In this paper we wish to draw the attention of the reader to a different problem of the CGC approach. In the CGC approach the scattering amplitude does not exceed the unitarity limit due to shadowing corrections. The Balitsky-Kovchegov \cite{BK} non-linear equation, which sums the `fan' diagrams of interacting BFKL Pomerons (see \fig{fan}-a), generates the amplitude, which tends to unity (the unitary limit) at high energies.  However, it has been shown in Refs.\cite{LRS,BART0,BARTU,LLS,LLR}
  that  Green's function of the exchange of $n$ BFKL Pomerons does not increase as $\exp\left( n\,\Delta_{\mbox{\tiny BFKL}}Y\right)$ where $\Delta_{\mbox{\tiny BFKL}}$   is the intercept of the BFKL Pomeron and $Y$ is the rapidity of two colliding dipoles. It turns out that  Green's function grows as
   $\exp\left( \frac{n^2}{N^2_c}\Delta_{\mbox{\tiny BFKL}}Y\right)$ ($N_c$ is the number of colours)  and this increase cannot be suppressed by the shadowing corrections
   \footnote{ Actually,  in  Refs.\cite{LRS,BART0,BARTU,LLS,LLR}   the exchange of  $n$-BFKL Pomerons were considered in the double log  approximation, which leads to the anomalous dimension $\frac{\bar{\alpha}_S}{(N^2-1)^2\omega} n^4$ instead of   $\frac{\bar{\alpha}_S}{\omega} n^2$ for the exchange of $n$ BFKL Pomerons .}. The $N_c$ suppression gives rise for the hope that  such an increase can manifest itself only at very high energies, but recently  A. Kovner and  M. Li showed that for nucleus-nucleus scattering in the CGC approach there exist  large $1/N^2_c$  corrections that  are larger than the shadowing suppression\footnote{ We thank A, Kovner and M. Li for sharing with us their finding.}. In spite of these
    results for the CGC approach, the references\cite{LRS,BART0,BARTU,LLS,LLR}
    show that these large corrections can be treated in the BFKL Pomeron calculus if we introduce the vertices of interaction of four Pomerons ($I\!\!P + I\!\!P \to I\!\!P + I\!\!P$).

     \begin{figure}[ht]
    \centering
  \leavevmode
      \includegraphics[width=9cm]{fan.pdf}
      \caption{\fig{fan}-a:  the `fan' diagrams of the Balitsky-Kovchegov non-linear equation. The wavy lines denote the BFKL Pomerons. The blob shows the triple Pomeron vertex. \fig{fan}b: is the Green's function of the exchange of $n$ BFKL Pomerons.}
\label{fan}
   \end{figure}
  Bearing this in mind we wish to return to discussion of these corrections  in
 simple, but exactly solvable, two dimensional models\cite{MUSA,ACJ,AAJ,JEN,ABMC,CLR,CIAF,BIT,RS,KLremark2,SHXI,KOLEV,nestor,LEPRI,utm,utmm}.   In these models we can sum Pomeron loops and we have the same kind of corrections, which go under slang name  of many particle Regge poles      \cite{DISASTER}. In the next section we will discuss these corrections in details. Here we wish to tell that the main goal of this paper is to find the scattering amplitude taking into account the $I\!\!P + I\!\!P \to I\!\!P + I\!\!P$ vertices   in the simple two dimensional models to get experience what we can expect in QCD in the CGC approach for these  disastrous contributions.

\section{Setting the problem}

 Green's function for one Pomeron in the two dimensional model  can be viewed as a sum of the contributions:
   \begin{equation}   \label{POM1}
 G_{I\!\!P}\left( Y\right) \,\,=\,\,\sum_{k=0}^\infty\,\frac{\left( \Delta\,Y\right)^n}{n!}\,\,=\,\,e^{\Delta\,Y}
 \end{equation}
 \eq{POM1} sums the ladder diagrams in leading log(1/x) approximation in which
all produced dipoles have a strong ordering in the fractions of total momentum $x_i$:
  \begin{equation} \label{POM2}
  1\,\gg\,x_1\,\gg\,x_2  \,\gg\, \dots \,x_i\,\gg\,x_{i+1} \,\gg\,\dots\,\gg\,\,x_{n-1}\,\gg\,x_n
  \end{equation}
this leads to $1/n!$ in \eq{POM1}. The two Pomeron exchange is shown in \fig{2pom}-a and  corresponds to the Feynman diagram of \fig{2pom}-b in which all dipoles emitted by the ladder 1 are absorbed by ladder 1', and all  dipoles produced  by ladder 2 are absorbed by ladder 2'. These diagrams lead to the Green's function of the exchange of two Pomerons: $G_{2 I\!\!P} =G^2_{I\!\!P}\left( Y\right) $.  However, the diagrams in which  the dipole,  emitted from ladder 1, will be absorbed by ladder  2'  are not small in the two dimensional models. We are going to call these diagrams  ``switch diagrams" , using the terminology suggested in Refs.\cite{LRS,BART0,BARTU,LLS,LLR}.  Indeed, after first exchange ladder 1 and 2' will give the Pomeron exchange, leading to the diagram of \fig{2pom}-d. Since at given rapidity we have two `switch diagrams: dipole emitted between ladders 1 and 2' and the dipole emitted between ladders 2 and 1' we can  obtain the vertex $I\!\!P + I\!\!P \to I\!\!P + I\!\!P$ , which is equal to $ 2 \Delta$.

The diagrams of \fig{2pom}-d can be summed  in $\omega $-representation:

\begin{equation} \label{POM3}
G_2\left( \omega \right)\,\,=\,\,\sum_{k=0}^\infty \frac{1}{\omega\,-\,2\,\Delta} \left( \frac{2\,\Delta}{\omega\,-\,2\,\Delta}\right)^k\,=\,\frac{1}{\omega\,-\,4 \,\Delta}
\end{equation}
Where $1/\left( \omega - 2 \Delta\right)$ is the contribution of the exchange of two Pomerons (see \fig{2pom}-b).  Coming back to $Y$ representation one can see that $G_2\left( Y\right)\,\,=\,\,\exp\left( 4\,\Delta\,Y\right)$ instead of $G_2\left( Y\right)\,\,=\,\,\exp\left( 2\,\Delta\,Y\right)$  which is expected for two Pomeron exchange.
     \begin{figure}[ht]
    \centering
    \begin{tabular}{c c}
  \leavevmode
      \includegraphics[width=14cm]{2pom.eps}& \includegraphics[width=3.9cm]{Npom.pdf}\\
      \end{tabular}
         \caption{ Summing the switch diagrams for the exchange of two BFKL Pomerons.
      \fig{2pom}-a: the exchange of two BFKL Pomerons. \fig{2pom}-b: the exchange of two BFKL Pomeron in the two dimensional model.  \fig{2pom}-c: the switch diagram for two ladder exchange in which the produced dipoles from the Pomeron 1 are absorbed by the Pomeron 2' and vise versa.  \fig{2pom}-d: the sum of all diagrams in the two dimensional model that contribute in the Green's  function of the exchange of two Pomerons in t-channel.  \fig{2pom}-e: the Green's function of the exchange of n$n$ Pomerons.The wavy lines denote the exchange of one Pomeron with the Green's function $\exp\left( \Delta Y\right)$. The $I\!\!P + I\!\!P \to I\!\!P + I\!\!P$ vertex is denoted by blob and equals to $2 \Delta$.}
\label{2pom}
   \end{figure}
Summing the diagrams of \fig{2pom}-e we obtain the Green's function for the exchange of $n$ Pomerons: $G_n\left( Y\right)\,=\,\,\exp\left( n^2 \Delta Y\right)$. This Green's function was derived in Ref.\cite{DISASTER} for the parton approach to high energy scattering. In QCD the structure of the corrections remain to be the same and the only difference  is that the `switch' diagrams as well the intercept of $n$ Pomeron state have the smallness of the order of $1/N^2_c$. Therefore, the $I\!\!P + I\!\!P \to I\!\!P + I\!\!P$ vertex in QCD is equal to $\frac{2}{N^2_c}\,\Delta_{\mbox{\tiny BFKL}}$. leading to
$G_n\left( Y\right)\,=\,\,\exp\left( n \Delta\,Y\,+\, \frac{n(n - 1)}{N^2_c} \Delta Y\right)$.

\section{ The structure of the parton cascade}

The useful tool  for discussing  the structure of the parton cascade is the generating function, which has the following form\cite{MUDI,LELU}:
\begin{equation} \label{Z}
Z\left( Y, u\right)\,\,=\,\,\sum_{n}P_n(Y) \,u^n
\end{equation}
where $P_n(Y) $ is the probability to find $n$ dipoles with rapidity $Y$. For the scattering of one dipole  with the target we have the following initial and boundary conditions:
\begin{equation} \label{IC}
\mbox{initial condition:}\,\, Z\left( Y=0,u\right) = u; ~~~~~\mbox{boundary  condition:}\,\, Z\left( Y,u=1\right)\,=\,1;
\end{equation}
The boundary condition follows from $P_n(Y)$ being probabilities.

We need to write the evolution equation for $Z\left( Y, u\right)$ taking into account  two Pometon vertices: $\Gamma_{3 I\!\!P} = \Delta $ for $I\!\!P \to I\!\!P + I\!\!P$ and $\Gamma^{2 I\!\!P}_{2 I\!\!P}  \,=\,2 \,\Delta $ for $I\!\!P  + I\!\!P \to I\!\!P + I\!\!P$. Fortunately  the equation for $Z$ for such a cascade has been written in Ref.\cite{LELU} and it takes the form:
\begin{equation} \label{ZEQ}
\frac{ \partial\,Z\left( Y,u\right)}{\partial Y}\,\,=\,\,- \Delta u (1-u) \frac{ \partial\,Z\left( Y,u\right)}{\partial u}\,\,+\,\,\Delta\,u (u - 1)  \frac{ \partial^2\,Z\left( Y,u\right)}{\partial u^2}
\end{equation}
As has been mentioned that \eq{ZEQ} has been suggested in Ref.\cite{LELU} , but \eq{ZEQ}  has a remarkable difference in comparison with equations in this reference: the sign in front of $u^2 \partial^2/\partial\,u^2$ is plus, not minus as in our previous attempts.

 Below we discuss a bit different equation:
 \begin{equation} \label{ZEQ1}
\frac{ \partial\,Z\left( \tilde{Y},u\right)}{\partial \tilde{Y}}\,\,=\,\,- \ u (1-u) \frac{ \partial\,Z\left( \tilde{Y},u\right)}{\partial u}\,\,+\,\,\kappa\,u (u - 1)  \frac{ \partial^2\,Z\left( \tilde{Y},u\right)}{\partial u^2}
\end{equation}
where $\tilde{Y} = \Delta\,Y$ and factor $\kappa = 1/N^2_c$ takes into account a suppression for $\Gamma^{2 I\!\!P}_{2I\!\!P} = 1/N^2_c$  in QCD. For two dimensional models $\kappa = 1$.

\subsection{BFKL cascade}

Neglecting $\Gamma^{2 I\!\!P}_{2 I\!\!P}$ we obtain the well known equation for the BFKL cascade.
 \begin{equation} \label{ZEQ2}
\frac{ \partial\,Z\left( \tilde{Y},u\right)}{\partial \tilde{Y}}\,\,=\,\,- \ u (1-u) \frac{ \partial\,Z\left( \tilde{Y},u\right)}{\partial u}
\end{equation}
It is instructive to observe that \eq{ZEQ2}  leads to  a non-linear equation for $ Z\left( Y, u\right)$\cite{MUDI,LELU}.   Indeed, the general solution to \eq{ZEQ2}  is of the form  $Z (Y, u) = Z(u(Y ))$; if we substitute this function into \eq{ZEQ2} , the derivatives  $\partial Z/\partial Y$ on the l.h.s. and r.h.s. of \eq{ZEQ2} cancel, and we obtain a differential equation for the function $u(Y)$.  Using the initial condition of \eq{IC} we can re-write
\eq{ZEQ2} in the form:
\begin{equation} \label{ZEQ3}
\frac{ \partial\,Z\left( \tilde{Y},u\right)}{\partial \tilde{Y}}\,\,=\,\,- \, Z\left( Y, u\right)  \,\,+\,\,Z^2\left( Y, u\right)
\end{equation}
Note, that the scattering amplitude in our models is equal to $N\left( Y\right)\,\,=\,\,1\,\,-\,\,Z\left( Y,1-\gamma\right)$ where $\gamma$ is the amplitude of the interaction of the dipole with the target at low energy.  For $N\left( Y\right)$ we  have the nonlinear equation:
\begin{equation} \label{ZEQ4}
\frac{ \partial\,N\left( \tilde{Y}\right)}{\partial \tilde{Y}}\,\,=\,\,- \, N\left( Y\right)  \,\,+\,\,N^2\left( Y\right)
\end{equation}
which is the Balitsky-Kovchegov equation\cite{BK}  for our simple models.

The solution to \eq{ZEQ3} , which satisfies the initial and boundary conditions of \eq{IC} has the following form:
\begin{equation} \label{ZEQ5}
Z\left( Y, u\right)\,\,=\,\,\frac{u\, e^{-\Delta \,Y}}{1\,\,+\,\,u\left( e^{-\Delta \,Y}\,\,-\,\,1\right)}
\end{equation}
One can see that at large values of $Y$  $Z \,\to\,0 $  ($N \,\to\,1$). In other words, the nonlinear corrections, which stem from  triple Pomeron interactions , suppress the increase of the scattering amplitude ($N \propto\,e^{\Delta\,Y}$) and lead to the scattering amplitude which reaches the unitarity bound. We call this phenomenon the saturation of parton densities.

\subsection{Asymptotic solution}

 The general solution to \eq{ZEQ2}   has been found in Ref.\cite{KOLEV}. However, before  applying the developed technique to this particular equation we wish to see a qualitative changes in the behaviour of $Z$ at large values of $Y$ that stems from a different sign of $\Gamma^{2 I\!\!P}_{2 I\!\!P}$ than in previous attempts to develop a similar cascade.

For this purpose we are going to find the asymptotic solution at large $Y$ from the following equation:
\begin{equation} \label{5}
\frac{\partial\,Z^{\rm asymp}\left( \tilde{Y},u\right)}{\partial u}\,\,+\,\,\kappa  \frac{ \partial^2\,Z^{\rm asymp}\left( \tilde{Y},u\right)}{\partial u^2}\,\,=\,\,0
\end{equation}
It has an obvious solution
\begin{equation}\label{ZAS}
 Z^{\rm asymp}\left( u\right)  =\,\,\frac{1\,\,-\,\, e^{-\frac{ u}{\kappa}}}{1\,\,-\,\, e^{-\frac{ 1}{\kappa}}},
 \end{equation}
which satisfies the boundary condition : $Z^{\rm asymp}\left(  u=1\right)  \,\,=\,\,1$  and  $Z^{\rm asymp}\left( u\right) \propto u $ at $u \ll  1$.
To find, how the solution approaches the asymptotic one, we  are looking for the solution in the form:   $ Z\left( Y, u\right)\,\,=\,\,\left( 1\,-\,\,e^{ - \frac{u}{\kappa} -\phi(Y,u)}\right)/\left( 1\,\,-\,\, e^{-\frac{ 1}{\kappa}}\right)$, assuming that $\phi''_{uu}$ and $\phi'^2_u$ are small. The equation for $\phi$ takes the form:

\begin{equation} \label{6}
\phi'_{\tilde{Y}}\left( \tilde{Y},u\right)\,\,=\,\,u\,(1- u)\,\phi'_u\left( \tilde{Y},u\right)
\end{equation}
with the initial condition at Y=0:

\begin{equation} \label{IC}
\phi\left( \tilde{Y}=0,u\right)\,\,=\,\, \left( \frac{1}{C}\,-\,\frac{1}{\kappa} \right)\,u (1 - u) ~~~\mbox{with}~~C= \frac{1}{ 1 - \exp\left( 1/\kappa\right)}
\end{equation}

which corresponds to $Z\left( Y=0,u\right) = u$ and $Z\left( Y=0,u=1\right) = 1$.  A general solution to \eq{6} has the form:
\begin{equation} \label{7}
\phi\left( Y,u\right)\,\,=\,\,\Phi\left( \tilde{Y}  + \ln \frac{u}{1 - u}\right)
\end{equation}
where the arbitrary function $\Phi$ has to be found from the initial conditions of \eq{IC}.

Finally, the solution takes the form:
\begin{equation} \label{8}
 Z\left( Y, u\right)\,\,=\,\,\frac{1}{1 - e^{-1/\kappa}} \Bigg( 1 - \exp\left( - \frac{u}{\kappa}-
 \left( ( 1 - \exp\left( -1/\kappa\right)) -  \frac{1}{\kappa} \right)
 \frac{u(1-u)\,e^{\tilde{Y}}}{ (1 - u + u e^{\tilde{Y}})^2} \right)\Bigg)
\end{equation}

\fig{zz} shows that this solution approaches  the unitarity limit, giving a hope that the shadowing corrections could  suppress the increase of the Green's function of the $n$-Pomerons in $t$-channel. However, this figure shows  that  $Z\left( \tilde{Y},u\right)  \,<\,0$ in the limited range of $u$, where this solution violates the unitarity constraints.
Hence, we need to find exact solution for the final conclusions.

     \begin{figure}[ht]
    \centering
  \leavevmode
  \begin{tabular}{c c}
      \includegraphics[width=8cm]{ZvxUFixY.pdf} &  \includegraphics[width=8cm]{ZvsUFixY8.pdf}\\
      \fig{zz}-a & \fig{zz}-b\\
      \end{tabular}
      \caption{Graphic form of the solution of \eq{8}.\fig{zz}-a: $\kappa =1$. \fig{zz}-b: $\kappa =1/8$. $\Delta = 0.2$. }
\label{zz}
   \end{figure}
\subsection{Exact solution}

Fortunately, the general solution to \eq{ZEQ2} has been found in Ref.\cite{KOLEV}.
We discuss this solution here, repeating all steps of Ref.\cite{KOLEV} and paying special attention to the sign in the r.h.s. of \eq{ZEQ2}. First, we consider $Z\left( Y, u\right)$ in $\omega$ representation:

\begin{equation} \label{GSOL1}
Z\left( Y, u\right)\,\,=\,\,\int\limits^{\epsilon + i \infty}_{\epsilon-i \infty} \frac{d \omega}{2\,\pi\,i} e^{ \omega\,\tilde{Y}} z\left( \omega, u\right)
\end{equation}
For $z\left( \omega, u\right)$ \eq{ZEQ2} takes the form:
\begin{equation} \label{GSOL2}
\omega\,z\left( \omega, u\right)\,\,=\,\,- u(1-u) \Big(  \frac{\partial\, z\left( \omega, u\right)}{\partial\,u}\,\,+\,\,\kappa \frac{\partial^2\, z\left( \omega, u\right)}{\partial\,u^2}\Big)
\end{equation}
Plugging $z\left( \omega, u\right)\,\,=\,\,\exp\left( - \frac{u}{2 \kappa}\right) {\tilde z}\left( \omega, u\right)$ in \eq{GSOL2} we obtain:
\begin{equation} \label{GSOL3}
4\,\kappa\,\omega\,{\tilde z}\left( \omega, u\right)\,\,=\,\,-u(1-u) \Big( - {\tilde z}\left( \omega, u\right)\,\,+\,\,4\kappa^2 \,{\tilde z}''_{uu}\left( \omega, u\right)\Big)
\end{equation}
Introducing $\tilde{z}\left( \omega, u\right)= u(1-u) {\cal G}\left( \omega, u\right) =  \frac{1 - v^2}{4} {\cal G}\left( \omega, v\right)$ with $ 1- 2u = v$, we can rewrite \eq{GSOL3} in the following form:
\begin{eqnarray} \label{GSOL4}
4\,\kappa\,\omega{\cal G}\left( \omega, u\right)\,\,&=&\,\,u (1 - u){\cal G}\left( \omega, u\right)\,\,-\,\,4\kappa^2 \left( u(1-u)\, {\cal G}\left( \omega, u\right)\Rb''_{u u}\nonumber\\
&=& \left( u(1-u) + 8\kappa^2\right){\cal G}\left( \omega, u\right)\,+\,8\,\kappa^2 (2 u - 1) {\cal G}'_u\left( \omega, u\right)\,\,-\,\,4\kappa^2 u (1-u){\cal G}''_{uu}\left( \omega, u\right)\end{eqnarray}
Using $v$  we have
\begin{equation} \label{GSOL5}
(1 - v^2) \,{\cal G}''_{vv}\left( \omega, v\right)\, -\,4\, v \,{\cal G}'_v\left( \omega, v\right)\,+\,\Bigg\{-2 +\frac{1}{\kappa} \omega- \frac{1 - v^2}{4\,\kappa^2}  \Bigg\}\,{\cal G}\left( \omega, v\right)\,\,=\,\,0
\end{equation}
In Ref.\cite{KOLEV} it is noted, that  function ${\cal G}\left( \omega, u\right)$ are intimately related to prolate spheroidal wave functions  $S_{n,m}\left( c, v\right) = (1- v^2)^{\frac{m}{2}}\,{\cal G}\left( c, v\right)$\cite{ABST,SPHF}, which satisfy the following equation:

\begin{equation} \label{GSOL6}
\frac{d}{d v}\Bigg( (1- v^2)\frac{d\,S_{n,m}\left( c, v\right)}{d \,v}\Bigg)\,\,+\,\,\left( \lambda^m_n \,-\,c^2\,v^2 - \frac{m}{1 - v^2}\right) S_{n,m}\left( c, v\right)\,\,=\,\,0
\end{equation}
with $n$ and $m$ being integer numbers.

For functions  ${\cal G}\left( c, v\right)$   \eq{GSOL6} takes  the form:
\begin{equation} \label{GSOL7}
(1 - v^2)\, {\cal G}''_{v v}\left( c, v\right)\,\,-\,\,2 (m + 1) v \,{\cal G}'_v\left( c, v\right)\,\,+\,\,\Big\{ \lambda^m_n\,-\,c^2\,v^2 - m (m+1)\Big\} {\cal G}\left( c, v\right)\,\,=\,\,0
\end{equation}
 Comparing \eq{GSOL5} and \eq{GSOL7} we obtain that
 \begin{equation} \label{GSOL8}
  m =1;~~c^2= - \frac{1}{4 \,\kappa^2};~~\omega_n \,\,=\,\,\kappa \,\lambda^1_n + \frac{1}{4\,\kappa};
  \end{equation}
  Hence, the set of the eigenfunctions for the generating function $Z\left( \omega, u\right)$ has the form:
  \begin{equation} \label{GSOL9}
  Z_n\left( v\right)\,\,=\,\,\, \sqrt{1 - v^2}\,e^{ - \frac{ 1- v}{4 \,\kappa}}\,S_{n,1}\left( \frac{i}{2 \,\kappa},v\right)
  \end{equation}

  Going back to rapidity representation, one can see that

    \begin{equation} \label{GSOL10}
  Z_n \left( Y, v\right)\,\,=\,\,\,e^{ \omega_n\,\tilde{Y}} \, \sqrt{1 - v^2}\,e^{ - \frac{ 1- v}{4 \,\kappa}}\,S_{n,1}\left( \frac{i}{2 \,\kappa},v\right) \,=\,\,e^{ \kappa \lambda^1_n\,\tilde{Y}}\sqrt{1 - v^2}\,e^{ - \frac{ 1- v}{4 \,\kappa}}\,S_{n,1}\left( \frac{i}{2 \,\kappa},v\right)
  \end{equation}

  Therefore, one can see  that each eigenfunction increases as function of $\tilde{Y}$  (see Table I for the values of $\omega_n$).
  \begin{table}
  \begin{tabular}{|c|c|l|l|l|l|l|}
  \hline
  $\omega$ & $\kappa$ & 1&2&3&4&5\\
  \hline
  $\omega_n$ & 8 ($N_c=3$)&1.636&1.766&2.616&3.579&4.8\\
  \hline
  $\omega_n$ & 1 ($N_c=1$)&2.2&6.14&12.133&20.113&30.128\\
  \hline
  \end{tabular}
  \caption{First five eigenvalues $\omega_n $ of  $\omega_n = \kappa\,\lambda^1_n$.}
  \end{table}
  Generally speaking the generating function is equal to
  \begin{equation} \label{GSOL11}
  Z\left( Y, v\right)= Z^{\rm asymp} \left( u\right)\,+\,\sum_{n=1}^\infty   {\rm C_n} e^{ \omega_n\,\tilde{Y}} \, \sqrt{1 - v^2}\,e^{ - \frac{ 1- v}{4 \,\kappa}}\,S_{n,1}\left( \frac{i}{2 \,\kappa},v\right)
  \end{equation}
  where
  \begin{equation} \label{GSOL12}
  {\rm C_n}\,\,=\,\,\int^1_{-1} d v \Bigg( Z\left( Y=0, u\right) - Z^{\rm asymp} \left( u\right)\Bigg) \,\frac{e^{  \frac{ 1- v}{4 \,\kappa}}  }{ \sqrt{1 - v^2}}\,\,\frac{S_{n,1}\left( \frac{i}{2 \,\kappa},v\right)  }{  ||S_{n,1}\left( \frac{i}{2 \,\kappa},v\right)||}
  \end{equation}
   with

  \begin{equation} \label{GSOL13}
  ||S_{n,1}\left( \frac{i}{2 \,\kappa},v\right)|| \,\,=\,\,\int^1_{-1} \!\!\!d v\, | S_{n,1}\left( \frac{i}{2 \,\kappa},v\right)|^2\,\,=\,\,\frac{2\,n\,(n+1)}{2\,n\,+\,1}
  \end{equation}

  $Z^{\rm asymp} \left( u\right)$ is the solution to \eq{5} which correponds to the minimal value of $\omega =0$. It is convenient to choose it in the following form:
  \begin{equation} \label{GSOL14}
    Z^{\rm asymp} \left( u\right)\,\,=\,\,\frac{1 -e^{-\frac{u}{\kappa}}}{1 - e^{-\frac{1}{\kappa}} }
    \end{equation}
    This solution gives $ Z^{\rm asymp} \left( u=1\right) =1$ and $Z^{\rm asymp} \left( u=0\right) =0$. Hence  the difference $ Z\left( Y=0, u\right) - Z^{\rm asymp} \left( u\right)$   satisfies the following boundary condition:
      \begin{equation} \label{GSOL15}
     Z\left( Y=0, u =0\right) - Z^{\rm asymp} \left( u=0\right)  = 0;~~~       Z\left( Y=0, u =1\right) - Z^{\rm asymp} \left( u=1\right) = 0;
     \end{equation}
     This difference can be expanded in the series of \eq{GSOL11} since functions $S_{n,1}\left( \frac{i}{2 \,\kappa},v\right)$ satisfy the boundary conditions of \eq{GSOL15} being equal to zero at $v= \pm 1$\footnote{We wish to refer to Refs.\cite{KOLEV,POLY} for more detailed analysis  of solutions to  since \eq{ZEQ2} belong to this class.}

 From \eq{GSOL11} we conclude that $Z\left( Y, u\right)$ increases with energy (rapidity Y). In other words even in the simple case of the initial condition of \eq{IC} , which corresponds to the deep inelastic scattering the shadowing corrections  failed to stop the increase of the scattering amplitude.

\subsection{Energy growth from the general structure of  equation.}

  In this section we demonstrate that the energy increase    actually stems from the general structure of \eq{ZEQ2}. Indeed, \eq{GSOL2} can be rewritten in the form of the
 Sturm-Liouville equation:
 \begin{equation} \label{GS1}
  s(u) \,\omega \,Z\left( \omega, u\right) + \frac{d}{d\,u} \Bigg( p(u) Z'_u\left( \omega, u\right)\Bigg)\,\,=\,\,0 ~~~\mbox{with}~~s(u) =\frac{\kappa}{u (1 - u)}e^{\frac{u}{\kappa}}~~\mbox{and} ~~p(u) = e^{\frac{u}{\kappa}}  \end{equation}

  The Sturm- Liouville equation has the following general features\cite{KOLEV,POLY}:
  \begin{enumerate}
  \item\quad \eq{GS1} has infinite set of eigenvalues $\omega_m = \lambda_n$.  $\lambda_n$ monotonically increases with $n$ with $\lambda _n \to \infty$  at large $n$. In our case of \eq{GS1} all $\lambda_n > 0$. The least value of  $\lambda_n$ is  $\lambda_0 =0$ which corresponds the asymptotic solution of \eq{ZAS}.
  \item\quad The multiplicity of each eigenvalue is equal to 1.
   \item \quad The eigenfunctions $Z_n(u)$ are orthogonal
   \begin{equation} \label{GS2}
   \int^1_0\!\!
  \!du\,s(u) \,Z_n(u)\,Z_m(u)\,\,=\,0 ~~~\mbox{for}~~n \neq n
   \end{equation}
  \item \quad For large $n$
  \begin{equation} \label{GS3}
  \lambda_n\,\,\frac{\pi^2\,n^2}{\delta^2} ~~\mbox{with}~~\delta = \int^1_0 \!\!\!d u \sqrt{\frac{s(u)}{p(u)}}
  \end{equation}
     \item \quad  For our equation
     \begin{equation} \label{GS4}
     \delta =\sqrt{\kappa} \int^1_0 \frac{ d u}{\sqrt{u(1-u)}} \,=\,\,\sqrt{\kappa}\,\pi
     \end{equation}
     giving
     \begin{equation} \label{GS5}
     \lambda_n = \,\,\kappa\,n^2
     \end{equation}
     \end{enumerate}

   Having \eq{GS5} we can conclude, that the generating function $Z_n\left( \tilde{Y},  u \right)$
    increases with $Y$. This feature is based on the general features of the Sturm-Liouville equation  and can be  stated  without finding the exact solution.  Since the class of Sturm-Lioville equations is much  wider than our particular equation  (see \eq{ZEQ2}) , we believe that more complicate equations in the case of QCD will still have these property,

\section{Conclusions}

   The main question, that we  have answered in this paper, whether   the shadowing correction can stop the steep increase of  Green's function for the exchange of the $n$ BFKL Pomerons: $G_{n I\!\!P}\left( Y\right) \propto \exp\left( \frac{n^2}{N^2_c} \Delta_{\mbox{\tiny BFKL}} Y\right)$.  In this paper we considering the simple  Pomeron calculus in zero transverse dimension. This approach has two great advantages: (i)  it  takes into account all shadowing corrections including the summation of the Pomeron loops and (ii) it has the same as in QCD striking increase of $G_{n I\!\!P}\left( Y\right) $. Solving exactly the evolution equation we demonstrate that the shadowing corrections cannot stop the increase of scattering amplitude which violates the unitarity constraints. Hence, our answer to the question in the title is positive.

   However,
    another phenomenon could considerably increase the shadowing corrections: the saturation effects inside the parton cascade. In the two dimensional Pomeron calculus such corrections have been included  in UTM model \cite{MUSA,BIT,utm,utmm}.  On one hand,  we have demonstrated in Ref.\cite{utmm} the scattering amplitude for  this model in spite of saturation in the parton cascade coincides with the BFKL cascades at high energies.  On the other hand, we need to include $\Gamma^{2 I\!\!P}_{2 I\!\!P}$ in this model. It has not been done and we consider this as a next problem to be solved.
  If we think about the theoretical realization of the parton model (see for example Ref.\cite{WEZU}), we do not expect that the unitarity would be violated. Hence, perhaps all our problems stem from the use of the Pomeron calculus and we need to find the theoretical description beyond this approach.

    The violation of unitarity stems from the vertex $\Gamma^{2 I\!\!P}_{2I\!\!P} >0$ , which also
    appears in QCD (see Refs.\cite{LRS,BART0,BARTU,LLS,LLR} ) with the same sign.
   We do not believe that the technical complications, coming with the QCD analysis, could provide the stronger shadowing than in the simple two dimensional models. However, we are aware that we have no idea how to sum Pomeron loops in QCD and how they influence the strength of   the shadowing corrections. On the other hand, we  need to consider a possibility to go out of the Pomeron calculus to treat the high energy  amplitude in CGC (see  for example Ref.\cite{LIEFF}).

   Concluding, we believe that CGC approach  correctly describes the high energy interaction at $N_c\to \infty$ and we do not have a clue how the shadowing correction could suppress the growth of the scattering amplitude in $1/N_c$  order.

~

     {\bf Acknowledgements}

   We thank our colleagues at Tel Aviv university and UTFSM for
 encouraging discussions. Special thanks go A. Kovner, M. Li  and M. Lublinsky for stimulating and encouraging discussions on the subject of this paper. This research    was supported  by
 ANID PIA/APOYO AFB180002 (Chile), and the Tel Aviv university encouragement grant \#5731.