\section{Introduction}

One of the key challenges in machine learning research concerns the ``generalization'' behavior of learning algorithms. That is: if a learning algorithm performs well on the training set, what guarantees can one provide on its performance on new samples?

While the question of generalization is understood in many settings (\citep{introToLearning,learningBook}), existing bounds and techniques provide vacuous expressions when employed to show the generalization capabilities of deep neural networks (DNNs) (\citep{normsNN,vcDimNN, genErrDL1,genErrDL2}). In general,
classical measures of model expressivity (such as Vapnik-Chervonenkis
(VC) dimension (\citep{vcDim}), Rademacher complexity (\citep{rademacherComplexity}), etc.), fail to explain the generalization abilities
of DNNs due to the fact that they are typically over-parameterized models with less
training data than model parameters.
A novel approach was introduced by \citep{explBiasMI}, and \citep{infoThGenAn} (further developed by \citep{conditionalMI,tighteningMI,fullVersionGeneralization,tpcGenerr} and many others), where information-theoretic techniques are used to link the generalization capabilities of a learning algorithm to information measures. These quantities are algorithm-dependent and can be used to analyze the generalization capabilities of general classes of updates and models, \textit{e.g.}, noisy iterative algorithms like Stochastic Gradient Langevin Dynamics (SGLD) (\citep{genErrMISGLD, optimalNoiseSGLD}), which can thus be applied to deep learning. Moreover, it has been shown that information-theoretic bounds can be non-vacuous and reflect the real generalization behavior even in deep
learning settings (\citep{genErrNoisyIterative,nonvacuousPAC, genErrSGDLDataDependent,genErrNoisyIterativeCMI}).

In this work we adopt and expand the framework introduced by \citep{genErrMISGLD}, but instead of focusing on the mutual information between the input and output of an iterative algorithm, we compute the maximal leakage (\citep{leakageLong}).
Maximal leakage, together with other information measures of the Sibson/R\'enyi family (maximal leakage can be shown to be Sibson Mutual information of order infinity (\citep{leakageLong})), have been linked to high-probability bounds on the generalization error (\citep{fullVersionGeneralization}). In particular, given a learning algorithm $\mathcal{A}$ trained on data-set $S$ (made of $n$ samples) one can provide the following guarantee in case of the $0-1$ loss:
\begin{equation}
    \mathbb{P}(|\text{gen-err}(\mathcal{A},S)|\geq \eta) \leq 2\exp(-2n\eta^2 + \ml{S}{\mathcal{A}(S)}).
\end{equation}
This deviates from much of the literature in which the focus is on bounding the \textbf{expected} generalization error instead (\citep{infoThGenAn,conditionalMI}).
Consequently, if one can guarantee that for a class of algorithms, the maximal leakage between the input and the output is bounded, then one can provide an \textbf{exponentially decaying} (in the number of samples $n$) bound on the probability of having a large generalization error. This is in general not true for mutual information, which can typically only guarantee a linearly decaying bound on the probability of the same event (\citep{learningMI}). Moreover, a bound on maximal leakage implies a bound on mutual information (cf.~\Cref{eq:sibsonMI-increasing}) and, consequently, a bound on the expected generalization error of $\mathcal{A}$. The main advantage of maximal leakage lies in the fact that it depends on the distribution of the samples only through its support. It is thus naturally independent from the distribution over the samples and particularly amenable to analysis, especially in additive noise settings.

The contributions of this work can be summarized as follows:
\begin{itemize}
    \item we derive novel bounds on $\ml{S}{\mathcal{A}(S)}$ whenever $\mathcal{A}$ is a noisy, iterative algorithm (SGLD-like), which then implies generalization with high-probability;
    \item we show that the bounds provided on maximal leakage strictly improve the bounds provided by \citep{genErrMISGLD}, and we thus provide a tighter bound on the expected generalization error of said algorithms as well;
    \item we show that, under certain assumptions, adding Gaussian noise is asymptotically optimal in the number of dimensions $d$. In particular, we prove that the maximal leakage (and, consequently, the mutual information) between the input and output of this family of algorithms goes to $0$ with the number of dimensions. This implies that the input and output are \textit{asymptotically independent} which is consistent with practical observations
 - larger neural networks generalize better;
  \item we leverage the analysis to extrapolate the optimal type of noise to be added (in the sense that it minimizes the induced maximal leakage), based on the assumptions imposed on the algorithm. In particular,
  \begin{itemize}
      \item if one assumes the $L_p$ norm of the gradient to be bounded, with $p\leq 2$, our analysis shows that adding Gaussian noise is asymptotically optimal;
      \item if one assumes the $L_\infty$ norm of the gradient to be bounded, then adding uniform noise is optimal;
  \end{itemize}
  Hence, the analysis and computation of maximal leakage can be used to inform the design of novel noisy, iterative algorithms.
\end{itemize}

\subsection{Related Work}
The line of work exploiting information measures to bound the expected generalization started in~\citep{explBiasMI,infoThGenAn} and was then refined with a variety of approaches considering Conditional Mutual Information~\citep{conditionalMI,genErrNoisyIterativeCMI}, the Mutual Information between individual samples and the hypothesis~\citep{ismi} or improved versions of the original bounds~\citep{ISIT2019,improvingMICondProc}. Other approaches employed the Kullback-Leibler Divergence with a PAC-Bayesian approach~\citep{pacbayes,nonvacuousPAC}.
Moreover, said bounds were then characterized for specific SGLD-like algorithms, denoted as ``noisy, iterative algorithms'' and used to provide novel, non-vacuous bounds for Neural Networks~\citep{genErrMISGLD,genErrSGDLDataDependent,genErrNoisyIterativeCMI, wang2023generalization} as well as for SGD algorithms~\citep{sgdBoundInfoTh}. Recent efforts tried to provide the optimal type of noise to add in said algorithms and reduce the (empirical) gap in performance between SGLD and SGD~\citep{optimalNoiseSGLD}.
All of these approaches considered the KL-Divergence or (variants of) Shannon's Mutual Information. General bounds on the expected generalization error leveraging arbitrary divergences were given in~\citep{tpcGenerr,convexAnalysisGenErr}.
Another line of work considered instead bounds on the probability of having a large generalization error~\citep{learningMI,fullVersionGeneralization,conditionalSibsonErrorBounds} and focused on large families of divergences and generalizations of the Mutual Information (in particular of the Sibson/R\'enyi-family, including conditional versions).

\section{Preliminaries, Setup, and a General Bound}
\subsection{Preliminaries}
\subsubsection{Information Measures}
The main building block of the information measures considered in this work is the R\'enyi's $\alpha$-divergence between two measures $P$ and $Q$, $D_\alpha(P\|Q)$ (which can be seen as a parametrized generalization of the Kullback Leibler-divergence)~\cite[Definition 2]{RenyiKLDiv}. Starting from R\'enyi's Divergence and the geometric averaging that it involves, Sibson built the notion of Information Radius~\citep{infoRadius} which can be seen as a special case of the following quantity~\citep{verduAlpha}: $
    I_\alpha(X,Y) = \min_{Q_Y} D_\alpha(P_{XY}\|P_{X}Q_Y). $
Sibson's $I_\alpha(X,Y)$ represents a generalization of Shannon's mutual information, indeed one has that: $
\lim_{\alpha\to 1}I_\alpha(X,Y)=I(X;Y) = \mathbb{E}_{P_{XY}}\left[\log\left(\frac{dP_{XY}}{dP_XP_Y}\right)\right].$
Differently, when $\alpha\to\infty$, one gets: \begin{equation}
I_\infty(X,Y)=\log\mathbb{E}_{P_Y}\left[ \esssup_{P_X} \frac{dP_{XY}}{dP_XP_Y} \right] =\ml{X}{Y},
\end{equation} where $\ml{X}{Y}$ denotes the maximal leakage from $X$ to $Y$, a recently defined information measure with an operational meaning in the context of privacy and security~\citep{leakageLong}. Maximal leakage represents the main quantity of interest for the scope of this paper, as  it is amenable to analysis and has been used to bound the generalization error~\citep{fullVersionGeneralization}. As such, we will bound the maximal leakage between the input and output of generic noisy iterative algorithms.

To that end, we mention a few useful properties of $\ml{X}{Y}$. If $X$ and $Y$ are jointly continuous random variables, then~\cite[Corollary 4]{leakageLong}
\begin{align}
    \ml{X}{Y} = \log \int \esssup_{P_X} f_{Y|X} (y|x) dy,
\end{align}
where $f_{Y|X}$ is the conditional pdf of $Y$ given $X$. Moreover, maximal leakage satisfies the following chain rule (the proof of which is given in Appendix~\ref{app:lem-chain-rule}):
\begin{lemma} \label{lem:chain-rule-maxLeak}
    Given a triple of random variables $(X,Y_1,Y_2)$, then
    \begin{align}
    \label{eq:lem-chain-rule}
    \ml{X}{Y_1, Y_2} \leq \ml{X}{Y_1} + \cml{X}{Y_2}{Y_1},
\end{align}
where the conditional maximal leakage $\cml{X}{Y_2}{Y_1}  = \esssup_{P_{Y_1}} \cml{X}{Y_2}{Y_1=y_1}$, where the latter term is interpreted as the maximal leakage from $X$ to $Y_2$ with respect to the distribution $P_{XY_2 | Y_1 = y_1}$. Consequently, for random variables $(X, (Y_i)_{i=1}^n)$,
\begin{align} \label{eq:lem:chain-rule-general}
    \ml{X}{Y^n} \leq \sum_{i=1}^n \cml{X}{Y_i}{Y^{i-1}}.
\end{align}
\end{lemma}
Moreover, one can relate $\ml{X}{Y}$ to $I(X;Y)$ through $I_\alpha$. Indeed, an important property of $I_\alpha$ is that it is non-decreasing in $\alpha$, hence for every $\infty>\alpha>1$:
\begin{equation} \label{eq:sibsonMI-increasing}
    I(X;Y) = I_1(X,Y) \leq I_\alpha(X,Y) \leq I_\infty(X,Y) = \ml{X}{Y}.
\end{equation}
For more details on Sibson's $\alpha$-MI we refer the reader to~\citep{verduAlpha}, as for maximal leakage the reader is referred to~\citep{leakageLong}.

\subsubsection{Learning Setting}

Let $\mathcal{Z}$ be the sample space, $\mathcal{W}$ be the hypothesis space, and $\ell : \mathcal{W} \times \mathcal{Z} \rightarrow \mathbb{R}_+$ be a loss function. Say $\mathcal{W} \subseteq \mathbb{R}^d$. Let $S=(Z_1, Z_2, \ldots, Z_n)$ consist of $n$ i.i.d samples, where $Z_i \sim P$, with $P$ unknown.
A learning algorithm $\mathcal{A}$ is a mapping $\mathcal{A}:\mathcal{Z}^n\to \mathcal{W}$ that given a sample $S$ provides a hypothesis $W=\mathcal{A}(S)$. $\mathcal{A}$ can be either a deterministic or a randomized mapping and undertaking a probabilistic (and information-theoretic) approach one can then equivalently consider $\mathcal{A}$ as a family of conditional probability distributions $P_{W|S=s}$ for $s\in \mathcal{Z}^n$ \textit{i.e.}, an information channel. Given a hypothesis $w\in\mathcal{W}$ the true risk of $w$ is denoted as follows:
\begin{equation}
    L_{P_Z}(w) = \mathbb{E}_{P}[\ell(w,Z)]
\end{equation}
while the empirical risk of $w$ on $S$ is denoted as follows:
\begin{equation}
    L_{S}(w) = \frac1n\sum_{i=1}^n \ell(w,Z_i).
\end{equation}
Given a learning algorithm $\mathcal{A}$, one can then define its generalization error as follows:
\begin{equation}
\text{gen-err}_\mathcal{P}(\mathcal{A},S)=L_{\mathcal{P}}(\mathcal{A}(S))-L_{S}(\mathcal{A}(S)).
\end{equation}
Since both $S$ and $\mathcal{A}$ can be random, $\text{gen-err}_\mathcal{P}(\mathcal{A},S)$ is a random variable and one can then study its expected value or its behavior in probability. Bounds on the expected value of the generalization error in terms of information measures are given in~\citep{infoThGenAn,ISIT2019,ismi,conditionalMI} stating different variants of the following bound (\cite[Theorem 1]{infoThGenAn}):
if $\ell(w,Z)$ is $\sigma^2$-sub-Gaussian\footnote{A $0$-mean random variable $X$ is said to be $\sigma^2$-sub-Gaussian if $\log\mathbb{E}[\exp(\lambda X)]\leq \sigma^2\lambda^2/2$ for every $\lambda\in\mathbb{R}$.} then
\begin{equation} \label{eq:gen-bound-MI}
    \left|\mathbb{E}[\text{gen-err}_\mathcal{P}(\mathcal{A},S)] \right| \leq \sqrt{\frac{2\sigma^2 I(S;\mathcal{A}(S))}{n}}.
\end{equation}
Thus, if one can prove that the mutual information between the input and output of a learning algorithm $\mathcal{A}$ trained on $S$ is bounded (ideally, growing less than linearly in $n$) then the expected generalization error of $\mathcal{A}$ will vanish with the number of samples.
Alternatively,~\cite{fullVersionGeneralization} demonstrate high-probability bounds, involving different families of information measures. One such bound, which is relevant to the scope of this paper is the following
\cite[Corollary 2]{fullVersionGeneralization}: assume $\ell(w,Z)$ is $\sigma^2$-sub-Gaussian and let $\alpha>1$
\begin{equation}
    \mathbb{P}(|\text{gen-err}_P(\mathcal{A},S)|\geq t) \leq 2\exp\left(-\frac{\alpha-1}{\alpha}\left(\frac{nt^2}{2\sigma^2}- I_\alpha(S,\mathcal{A}(S))\right) \right) \label{eq:genErrProbBound},
\end{equation}
taking the limit of $\alpha\to\infty$ in~\eqref{eq:genErrProbBound} leads to the following~\cite[Corollary 4]{fullVersionGeneralization}:
\begin{equation}
    \mathbb{P}(|\text{gen-err}_P(\mathcal{A},S)|\geq t) \leq 2\exp\left(-\left(\frac{nt^2}{2\sigma^2}- \ml{S}{\mathcal{A}(S)}\right)\right) \label{eq:genErrProbBoundMaxLeak}.
\end{equation}
Thus, in this case, if one can prove that the maximal leakage between the input and output of a learning algorithm $\mathcal{A}$ trained on $S$ is bounded, then the \textbf{probability} of the generalization error of $\mathcal{A}$ being larger than any constant $t$ will decay \textbf{exponentially fast} in the number of samples $n$.
\subsection{Problem Setup}
We consider iterative algorithms, where each update is of the following form:
\begin{align} \label{eq:iteration}
W_t = g(W_{t-1}) - \eta_t F(W_{t-1},Z_t) + \xi_t, ~ \forall ~ t \geq 1,
\end{align}
where $Z_t \subseteq S$ (sampled according to some distribution), $g: \mathbb{R}^{d} \rightarrow \mathbb{R}^{d}$ is a deterministic function, $F(W_{t-1},Z_t)$ computes a direction (e.g., gradient), $\eta_t$ is a constant step-size, and $\xi_t =( \xi_{t1}, \ldots, \xi_{td})$ is noise. We will assume for the remainder of this paper that $\eta_t$ has an absolutely continuous distribution.
Let $T$ denote the total number of iterations, $W^t = (W_1, W_2, \ldots W_t)$, and $Z^t = (Z_1, Z_2, \ldots, Z_t)$. The algorithms under consideration further satisfy the following two assumptions
\begin{itemize}
\item {\bf Assumption 1 (Sampling):} The sampling strategy is agnostic to parameter vectors:
\begin{align}
P(Z_{t+1}|Z^t, W^t, S) = P(Z_{t+1} |Z^t, S).
\end{align}
\item {\bf Assumption 2 ($\mathbf{L_p}$-Boundedness):} For some $p>0$ and $L>0$,  $\sup_{w,z} \| F(w,z) \|_p \leq L$.
\end{itemize}
As a consequence of the first assumption and the structure of the iterates, we get:
\begin{align} \label{eq:markov}
P(W_{t+1} | W^t, Z^T, S) = P(W_{t+1} | W_t, Z_{t+1}).
\end{align}
The above setup was proposed by Pensia \emph{et al.}~\citep{genErrMISGLD}, who specifically studied the case $p=2$.
Denoting by $W$ the final output of the algorithm (some function of $W^T$), they show that
\begin{theorem}[{\cite[Theorem 1]{genErrMISGLD}}] \label{thm:Pensia}
If the boundedness assumption holds for $p=2$ and
$\xi_t \sim \mathcal{N}(0, \sigma^2_t I_d)$, then
\begin{align} \label{eq:thm-Pensia}
I(S;W) \leq \frac{d}{2} \sum_{t=1}^T \log \left(1 + \frac{\eta^2_t L^2}{d \sigma^2_t} \right).
\end{align}
\end{theorem}
By virtue of inequality~\eqref{eq:gen-bound-MI}, this yields a bound on the expected generalization error.

In this work, we derive bounds on the maximal leakage between $\ml{S}{W}$ for iterative noisy algorithms, which leads to high-probability bounds on the generalization error (cf. equation~\eqref{eq:genErrProbBoundMaxLeak}). We consider different scenarios in which $F$ is bounded in $L_1$, $L_2$, or $L_\infty$ norm, and the added noise is Laplace, Gaussian, or Uniform.

\subsection{Notation} \label{sec:notation}
Given $d \in \mathbb{N}$, $w \in \mathbb{R}^d$, and $r >0$, let $\mathcal{B}_p^{d}(w,r) = \{ x \in \mathbb{R}^d: \| x - w \|_p \leq r \}$ denote the $L_p$-ball of radius $r$ and center $w$, and let $V_p (d,r)$ denote its corresponding volume. When the dimension $d$ is clear from the context, we may drop the superscript and write $\mathcal{B}_p(w,r)$. Given a set $S$, we denote its complement by $\overline{S}$. The $i$-th component of $w_t$ will be denoted by $w_{ti}$.

We denote the pdf of the noise $\xi_t$ by $f_t: \mathbb{R}^d \rightarrow \mathbb{R}$. The following functional will be useful for our study: given $d \in \mathbb{N}$, $p > 0$, a pdf $f: \mathbb{R}^d \rightarrow R$, and an $r \geq 0$, define
\begin{align}
    \label{eq:def-h}
    h(d,p,f,r) \vcentcolon=  \int \limits_{\overline{\mathcal{B}_p^d}(0, r)} \sup_{x \in \mathcal{B}_p^d(0,r)} f(w - x) \mathrm{d} w.
\end{align}
We denote the ``positive octant'' by $A_d$, i.e.,
\begin{align}
    \label{eq:def-positive-octant}
    A_d \vcentcolon= \{ w \in \mathbb{R}^d: w_i \geq 0, \text{ for all } i \in \{1,2,\ldots,d\} \}.
\end{align}
Since we will mainly consider pdfs that are symmetric (Gaussian, Laplace, uniform), the $h$ functional ``restricted'' to $A_d$ will be useful:
\begin{align}
    \label{eq:def-h-positive}
    h_+ (d,p,f,r) \vcentcolon=  \int \limits_{\overline{\mathcal{B}_p^d}(0, r) \cap A_d } \sup_{x \in \mathcal{B}_p^d(0,r)} f(w - x) \mathrm{d} w.
\end{align}

\subsection{General Bound}

\begin{proposition} \label{prop:general}
    Suppose $f_{t}: \mathbb{R}^d \rightarrow \mathbb{R}$ is maximized for $x = 0$. If Assumptions 1 and 2 hold for some $p > 0$, then
\begin{align}
    \label{eq:prop-general}
    \ml{S}{W} \leq \sum_{t=1}^T \log \left( f_{t} (0) V_p(d,\eta_t L) + h(d,p,f_t, \eta_t L) \right),
\end{align}
where $h$ is defined in equation~\eqref{eq:def-h}.
\end{proposition}

The above bound is appealing as it implicitly poses an optimization problem: given a constraint on the noise pdf $f_t$ (say, a bounded variance), one may choose $f_t$ as to minimize the upper bound in equation~\eqref{eq:prop-general}. Moreover, despite its generality, we show that it is tight in several interesting cases, including when $p=2$ and $f_t$ is the Gaussian pdf.
The series of steps leading to the upper bound include only one inequality (the source of the ``looseness'' of the bound), which could be viewed as due to replacing Assumption 2 by the following statement: for all $w$, $\cup_z F(w,z) = \mathcal{B}_p^d(0,L)$, i.e., in addition to assuming $F(w,z) \in \mathcal{B}_p^d(0,L)$, we assume that, for every $w$, every point in the ball $\mathcal{B}_p^d(0,L)$ is attained for some $z$.

In the next section, we consider several scenarios for different values of $p$ and different noise distributions. As a testament to the tractability of maximal leakage, we derive exact semi-closed form expressions for the bound of Proposition~\ref{prop:general}. Finally, it is worth noting that the form of the bound allows to choose different noise distributions at different time steps, but these examples are outside the scope of this paper. \\

\begin{proof}
We proceed as in the the work of \citep{genErrMISGLD}:
\begin{align}
\ml{S}{W} \leq \ml{Z^T}{W^T} \leq \sum_{t=1}^T \cml{Z^T}{W_t}{W^{t-1}} = \sum_{t=1}^T \cml{Z_t}{W_t}{W_{t-1}},
\end{align}
where the first inequality follows from Lemma 2 of \citep{genErrMISGLD} and the data processing inequality for maximal leakage~\cite[Lemma 1]{leakageLong}, the second inequality follows Lemma~\ref{lem:chain-rule-maxLeak}, and the equality follows from~\eqref{eq:markov}. Now,
\begin{align}
\exp\left\lbrace \cml{Z_t}{W_t}{W_{t-1}=w_{t-1}} \right\rbrace
    & =     \int_{\mathbb{R}^d} \esssup_{P_{z_t}} p(w_t | Z_t) \mathrm{d} w_t \\
    & = \int_{\mathbb{R}^d} \esssup_{P_{z_t}} f_t \left( w_t - g(w_{t-1}) + \eta_t F(w_{t-1},Z_t) \right) \mathrm{d} w_t, \\
    & = \int_{\mathbb{R}^d} \esssup_{P_{z_t}} f_t \left( w_t + \eta_t F(w_{t-1},Z_t) \right) \mathrm{d} w_t,
\end{align}
where the last equality follows from a change of a variable $\tilde{w_t} = w_t - g(w_{t-1})$. Finally, since $\eta_t F(w_{t-1},z_t) \in \mathcal{B}_p(0, \eta_t L)$ by assumption, we can further upper-bound the above by:
\begin{align}
     & \exp\left\lbrace  \cml{Z_t}{W_t}{W_{t-1}=w_{t-1}} \right\rbrace \\
     & \leq \int_{\mathbb{R}^d} \sup\limits_{x_t \in \mathcal{B}_p(0,\eta_t L)} f_t \left( w_t + x_t \right) \mathrm{d} w_t  \label{eq:prop-ball-Lp} \\
     & = \int\limits_{\mathcal{B}_p(0,\eta_t L) } \sup\limits_{x_t \in \mathcal{B}_p(0,\eta_t L)}  f_{t} \left( w_{ti} + x_{ti} \right) \mathrm{d} w_t + \int\limits_{\mathcal{B}_p^\mathsf{c}(0,\eta_t L)} \sup\limits_{ x_t \in \mathcal{B}_p(0,\eta_t L)} f_t \left( w_t + x_t \right) \mathrm{d} w_t \\
     &  = f_{t} (0) V_p(d, \eta_t L) + \int\limits_{ \overline{\mathcal{B}_p} (0,\eta_t L)} \sup\limits_{ x_t \in \mathcal{B}_p(0,\eta_t L)} f_t \left( w_t - x_t \right) \mathrm{d} w_t,
\end{align}
where the last equality follows from the assumptions on $f_{t}$.
\end{proof}

\section{Boundedness in $L_2$-Norm}

Considering the case where $F$ computes a gradient, then boundedness in $L_2$-norm is a common assumption. It is commonly enforced, for instance, using gradient clipping \citep{abadi2016tensorflow,abadi2016deep,chen2020understanding}.

The case in which $p=2$ and the noise is Gaussian leads to the strongest result in this paper:

\begin{theorem} \label{thm:gaussian-ml}
If the boundedness assumption holds for $p \leq 2$ and
$\xi_t \sim \mathcal{N}(0, \sigma^2_t I_d)$, then
\begin{align} \label{eq:thm-gaussian-ml}
\ml{S}{W} \leq  \sum_{t=1}^T \log \left( \frac{V_2 (d,\eta_t L)}{(2 \pi \sigma_t^2)^{d/2}} +  \frac{1}{\Gamma \left( \frac{d}{2} \right)  } \sum_{i=0}^{d-1} \Gamma \left( \frac{i+1}{2} \right)
\left(  \frac{\eta_t L}{ \sigma_t \sqrt{2}}  \right)^{d-1-i} \right),
\end{align}
where $\displaystyle V_2(d,r) = \frac{\pi^{d/2}}{\Gamma \left(\frac{d}{2}+1 \right)} r^d$. Consequently, for fixed $T$,
\begin{align} \label{eq:thm-gaussian-limit}
\lim_{d \rightarrow \infty} \ml{S}{W}=0.
\end{align}
\end{theorem}

Remarkably, equation~\eqref{eq:thm-gaussian-limit} states that, as $d$ grows, $S$ and $W$ are asymptotically independent. The bound is asymptotically optimal for $\ml{S}{W}$ (indeed it yields an equality). More importantly, the induced high probability bound by equation~\eqref{eq:genErrProbBoundMaxLeak} is also optimal. Indeed, at the limit when $S$ and $W$ are independent, the bound~\eqref{eq:genErrProbBoundMaxLeak} recovers (the order optimal) McDiarmid's inequality \textit{i.e.}, under the assumptions of~\Cref{thm:gaussian-ml} and considering the $0-1$ loss:
\begin{equation}
    \mathbb{P}(|\text{gen-err}(\mathcal{A},S)|\geq t) \leq 2\exp(-2nt^2).
\end{equation} This can be seen as an explanation of the (arguably un-intuitive) phenomenon that deeper networks often generalize better (also analyzed by~\citep{wang2023generalization}).

By contrast, this is not captured in the bound by~\citep{genErrMISGLD} given in equation~\eqref{eq:thm-Pensia}. Indeed, it is growing as a function of $d$, and tends to a non-zero value:
\begin{align}
    \lim_{ d \rightarrow \infty} \frac{d}{2} \sum_{t=1}^T \log \left(1 + \frac{\eta^2_t L^2}{d \sigma^2_t} \right) & = \sum_{t=1}^T \frac{\eta^2_t L^2}{2 \sigma^2_t},
\end{align}
where the equality follows from the fact that $\displaystyle \lim_{n \rightarrow \infty} (1+c/n)^c = e^c$ for all $c \in \mathbb{R}$. Notably, however, $I(S;W) \leq \ml{S}{W}$ (cf. equation~\eqref{eq:sibsonMI-increasing}), so that for large $d$, we have
\begin{align}
    I(S;W) \leq \ml{S}{W}\leq \text{right-hand side of~\eqref{eq:thm-gaussian-ml}} \leq \text{ right-hand side of~\eqref{eq:thm-Pensia}}.
\end{align}
As such, the bound in Theorem~\ref{thm:gaussian-ml} is also a tighter upper bound on $I(S;W)$.

Moreover, note that even if the parameter $L$ is large (e.g., Lipschitz constant of a neural network \citep{genErrSGDLDataDependent}), it appears in~\eqref{eq:thm-gaussian-ml} normalized by $\Gamma(d/2)$ so its effect is significantly dampened (as $d$ is also typically very large).

Finally, note that the bound in Proposition~\ref{prop:general} is increasing in $p$: this can be seen from line~\eqref{eq:prop-ball-Lp}, where the supremum over $\mathcal{B}_p$ can be further upper-bounded by a supremum over $\mathcal{B}_{p'}$ for $p' > p$. Therefore for $q \leq p$, the bound induced by Proposition~\ref{prop:general} is smaller. The bound in~\Cref{thm:gaussian-ml} corresponds to $p=2$ and goes to 0 (as $d$ grows), hence the bound induced by Proposition~\ref{prop:general} goes to 0 for all $q \leq p$. In particular, adding Gaussian noise is asymptotically optimal (in the sense discussed above) when Assumption 2 holds for any $p \leq 2$. \\

\begin{proof}
To show that the right hand side of Equation~\eqref{eq:thm-gaussian-ml} goes to zero as as $d\rightarrow\infty$ is equal to zero, we use Stirling's approximation of the Gamma function: for all $ x >0$,
\begin{align}
    \sqrt{2 \pi} x^{x-\frac{1}{2}} e^{-x} \leq \Gamma(x) \leq \sqrt{2 \pi} x^{x-\frac{1}{2}} e^{-x} e^{\frac{1}{12x}}.
\end{align}
The details of the computation can be found in Appendix~\ref{app:stirling}.
We now turn to the proof of inequality~\eqref{eq:thm-gaussian-ml}. The conditions of Proposition~\ref{prop:general} are satisfied, thus it is sufficient to prove the bound for $p=2$:
\begin{align}
    \ml{S}{W} & \leq \sum_{t=1}^T \log \left( f_{t} (0) V_2(d,\eta_t L) + \int \limits_{\overline{\mathcal{B}_2}(0, \eta_t L)} \sup_{x_t \in \mathcal{B}_2(0,\eta_{t} L)} f_t(w_t - x_t) \mathrm{d} w_t \right) \\
    & = \sum_{t=1}^T \log \!\! \left( \!\! \frac{V_2(d,\eta_t L)}{(2 \pi \sigma_t^2)^{\frac{d}{2}}} + \!\!\!\! \int  \limits_{\overline{\mathcal{B}_2} (0, \eta_t L)} \!\! \sup_{x_t \in \mathcal{B}_2(0,\eta_{t} L)} \!\! \frac{1}{(2 \pi \sigma_t^2)^{\frac{d}{2}}} \exp \left\lbrace -\frac{\|w_t-x_t\|^2_2}{ 2 \sigma_t^2} \right\rbrace \mathrm{d} w_t \!\! \right) \! . \label{eq:prop-for-gaussian}
\end{align}
Hence, it remains to show that the second term inside the $\log$ matches that of equation~\eqref{eq:thm-gaussian-ml}. To that end, note that the point in $\mathcal{B}_2(0,\eta_t L)$ that minimizes the distance to $w_t$ is given $\frac{\eta_t L}{\| w_t \|} w_t$. So we get
\begin{align}
\| w_t - x_t \| \geq \| w_t - \frac{\eta_t L}{\| w_t \|} w_t \| = \| w_t\| - \eta_t L.
\end{align}
 Then,
\begin{align}
h(d,2,f_t,\eta_t L) & =     \int \limits_{\overline{\mathcal{B}_2} (0, \eta_t L)} \sup_{x_t \in \mathcal{B}_2(0,\eta_{t} L)} \frac{1}{(2 \pi \sigma_t^2)^{\frac{d}{2}}} \exp \left\lbrace -\frac{\|w_t-x_t\|^2_2}{ 2 \sigma_t^2} \right\rbrace \mathrm{d} w_t \\
&  = \int\limits_{\overline{\mathcal{B}_2} (0, \eta_t L)}  \frac{1}{(2 \pi \sigma_t^2)^{\frac{d}{2}}} \exp \left\lbrace - \frac{(\|w_t\|_2 - \eta_t L )^2}{2 \sigma_t^2} \right\rbrace \mathrm{d} w_t.
\label{eq:almost-there}
\end{align}
To evaluate this integral, we use spherical coordinates (details in Appendix~\ref{app:spherical}). Then,
\begin{align} \label{eq:MlGaussInt}
h(d,2,f_t,\eta_t L) = \left ( \frac{\eta_t L}{ \sigma_t \sqrt{2}} \right)^{d-1} \frac{1}{\Gamma \left( \frac{d}{2} \right)  } \sum_{i=0}^{d-1} \left( \frac{ \sigma_t \sqrt{2}}{\eta_t L}  \right)^i \Gamma \left( \frac{i+1}{2} \right).
\end{align}
Combining equations~\eqref{eq:prop-for-gaussian} and~\eqref{eq:MlGaussInt} yields~\eqref{eq:thm-gaussian-ml}.
\end{proof}

\begin{remark}
One could also derive a semi-closed form bound for the case in which the added noise is uniform. However, in that case $\ml{S}{W}$ goes to infinity as $d$ goes to infinity. The same behavior holds if the added noise is Laplace. Since the Gaussian noise leads to an asymptotically optimal bound, we skip the analysis of uniform and Laplace noise.
\end{remark}

\section{Boundedness in $L_\infty$-Norm}

The bound in Proposition~\ref{prop:general} makes minimal assumptions about the pdf $f_t$. In many practical scenarios we have more structure we could leverage. In particular, we make the following standard assumptions in this section:
\begin{itemize}
    \item $\xi_t$ is composed of i.i.d components. Letting $f_{t0}$ be the pdf of a component, $\displaystyle f_t(x_t) \!\! = \!\! \prod_{i=1}^d f_{t0}(x_{ti})$.
    \item $f_{t0}$ is symmetric around 0 and non-increasing over $[0,\infty).$
\end{itemize}
In this setting, Proposition~\ref{prop:general} reduces to a very simple form for $p = \infty$:

\begin{theorem} \label{thm:Linf}
Suppose $f_{t}$ satisfies the above assumptions. If Assumptions 1 and 2 hold for $p = \infty$, then
\begin{align} \label{eq:thm-Linf}
\ml{S}{W} \leq \sum_{t=1}^T d \log \left( 1 + 2 \eta_t L f_{t0}(0) \right).
\end{align}
\end{theorem}

Unlike the bound of Theorem~\ref{thm:gaussian-ml}, the limit of $d$ to infinity here is infinite. However, the bounded-$L_\infty$ assumption is \textit{weaker} than assuming a bounded $L_2$-norm. Moreover, the assumption of having a bounded $L_\infty$-norm is satisfied in~\citep{pichapati2019adaclip} where the authors clipped the gradient in terms of the $L_\infty$-norm, thus ``enforcing'' the assumption. On the other hand, the theorem has an intriguing form as, under standard assumptions, the bound depends on $f_{t0}$ only through $f_{t0}(0)$. This naturally leads to an optimization problem: given a certain constraint on the noise, which distribution $f^\star$ minimizes $f(0)$? The following theorem shows that, if the noise is required to have a bounded variance, then $f^\star$ corresponds to the uniform distribution:

\begin{theorem} \label{thm:opt-noise}
 Let $\mathcal{F}$ be the family of probability densities (over $\mathbb{R}$) satisfying for each $f \in \mathcal{F}$:
 \begin{enumerate}
 \item $f$ is symmetric around 0.
 \item $f$ is non-increasing over $[0, \infty)$.
 \item ${\rm \mathbf  E}_f [X^2] \leq \sigma^2$.
 \end{enumerate}
Then, the distribution minimizing $f(0)$ over $\mathcal{F}$ is the uniform distribution $\mathcal{U}( -\sigma \sqrt{3}, \sigma \sqrt{3})$.
 \end{theorem}

That is, uniform noise is optimal in the sense that it minimizes the upper bound in Theorem~\ref{thm:Linf} under bounded variance constraints. The proof of Theorem~\ref{thm:opt-noise} is deferred to Appendix~\ref{app:opt-noise}.

\subsection{Proof of Theorem~\ref{thm:Linf}}
    Since the assumptions of Proposition~\ref{prop:general} hold, then
    \begin{align}
        \ml{S}{W} & \leq \sum_{t=1}^T \log \left( f_t(0) V_\infty(d,\eta_t L) + \int \limits_{\overline{\mathcal{B}_\infty}(0, \eta_t L)} \sup_{x_t \in \mathcal{B}_\infty(0,\eta_{t} L)} f_t(w_t - x_t) \mathrm{d} w_t \right) \\
        & = \sum_{t=1}^T \log \left( (2\eta_t L f_{t0}(0))^d  +
        \int \limits_{\overline{\mathcal{B}_\infty}(0, \eta_t L)} \prod_{i=1}^d \sup_{ x_{ti}: |x_{ti}| \leq \eta_t L} f_{t0} (w_{ti}-x_{ti}) \mathrm{d} w_t
        \right).
    \end{align}
It remains to show that $h(d,\infty,f_t,\eta_t L)$ (i.e., the second term inside the $\log$ in~\Cref{eq:def-h}) is equal to $(1+2\eta_t L f_{t0}(0))^d -(2\eta_t L f_{t0}(0))^d $. We will derive a recurrence relation for $h$ in terms of $d$. To simplify the notation, we drop the subscript $t$ and ignore the dependence of $h$ on $p=\infty$, $f_t$, and $\eta_t L$, so that we simply write $h(d)$ (and correspondingly, $h_+(d)$, cf.~\Cref{eq:def-h-positive}).

By symmetry, $h(d) = 2^d h_+(d)$. Letting $w^{d-1}:=(w_1,\ldots,w_{d-1})$, we will decompose the integral over $\overline{\mathcal{B}_\infty^d}(0,\eta_t L)$ into two disjoint subsets: 1) $ w^{d-1} \notin \mathcal{B}_\infty^{d-1}(0,\eta_t L)$, in which case $w_d$ can take any value in $\mathbb{R}$, and 2) $w^{d-1} \in \mathcal{B}_\infty^{d-1}(0,\eta_t L)$, in which case $w_d$ must satisfy $|w_d| > \eta_t L$.
\begin{align}
    h_+(d) &  =
    \int \limits_{\overline{\mathcal{B}^{d-1}_\infty}(0, \eta_t L) \cap A_{d-1} } \prod_{i=1}^{d-1} \sup_{ x_{i}: |x_{i}| \leq \eta_t L} f (w_{i}-x_{i})  \int_{0}^\infty \sup_{x_d: |x_d| \leq \eta_t L} f(w_d - x_d) \mathrm{d} w_d \mathrm{d} w^{d-1}  \label{eq:Linf-messy-1}\\
    & ~ + \int \limits_{{\mathcal{B}^{d-1}_\infty}(0, \eta_t L) \cap A_{d-1} } \prod_{i=1}^{d-1} \sup_{ x_{i}: |x_{i}| \leq \eta_t L} f (w_{i}-x_{i})  \int_{\eta_t L}^\infty \sup_{x_d: |x_d| \leq \eta_t L} f(w_d - x_d) \mathrm{d} w_d \mathrm{d} w^{d-1} \label{eq:Linf-messy-2}
\end{align}
The innermost integral of line~\eqref{eq:Linf-messy-2} is independent of $w^{d-1}$ so that the outer integral is equal to $h_{+}(d-1)$. Similarly, the innermost integral of line~\eqref{eq:Linf-messy-1} is independent of $w^{d-1}$, and the supremum in the outer integral yields $f(0)$ for every $i$. Hence, we get
\begin{align}
    h(d) = \left(1+2\eta_t L f(0) \right) h (d-1)+ (2 \eta_t L f(0))^{d-1}, \label{eq:recurrence-Linf}
\end{align}
the detailed proof of which is deferred to Appendix~\ref{app:recurrence-Linf}.
Finally, it is straightforward to check that $h(1) =1$, hence $h(d) = (1+2\eta_t L f(0))^d-(2\eta_t L f(0))^d$.

\section{Boundedness in $L_1$-Norm}

In this section, we consider the setting where Assumption 2 holds for $p=1$. By Proposition~\ref{prop:general}, any bound derived for $p=2$ holds for $p=1$ as well. In particular, Theorem~\ref{thm:gaussian-ml} applies so that $\ml{S}{W}$ goes to zero when the noise is Gaussian. Nevertheless, it is possible to compute a semi-closed form directly for $p=1$ (cf. Theorem~\ref{thm:l1-gaussian} below), which would be inherently tighter.

Considering the optimality of Gaussian noise for the $p=2$ case, and the optimality of uniform noise (in the sense discussed above) for $p = \infty$ case, one might wonder if Laplace noise is optimal for the $p=1$ case. We answer this question in the negative, as the limit of the leakage in this case is a non-zero constant (cf. Theorem~\ref{thm:l1-laplace}), as opposed to the zero limit when the noise is Gaussian.

\subsection{Bound for Laplace noise}
We say $X$ has a Laplace distribution, denoted by $X \sim \mathrm{Lap}(\mu,1/\lambda)$, if its pdf is given by $f(x)= \frac{\lambda}{2}e^{-\lambda |x-\mu|}$ for $x \in \mathbb{R}$, for some $\mu \in \mathbb{R}$ and $\lambda >0$. The corresponding variance is given by $2/\lambda^2$.
\begin{theorem} \label{thm:l1-laplace}
If the boundedness assumptions holds for $p=1$ and
$\xi_t$ is composed of i.i.d components, each of which is $\sim \mathrm{Lap}(0, \frac{\sigma_t}{\sqrt{2}})$, then
\begin{align} \label{eq:thm-l1-laplace}
\ml{S}{W} \leq  \sum_{t=1}^T \log \left( \frac{V_1 (d,\eta_t L)}{( \sigma_t \sqrt{2})^{d}} +  \sum_{i=0}^{d-1} \frac{(\sigma_t \eta_t L/ \sqrt{2})^i}{ i ! }  \right),
\end{align}
where $\displaystyle V_1(d,r)  = \frac{(2r)^d}{ d!}$. Consequently, for fixed $T$,
\begin{align} \label{eq:thm-l1-laplace-limit}
\lim_{d \rightarrow \infty} \ml{S}{W} \leq \sum_{t=1}^T \frac{\sigma_t \eta_t L}{ \sqrt{2}}.
\end{align}
\end{theorem}

\begin{proof}
    We give a high-level description of the proof (as similar techniques have been used in proofs of earlier theorems) and defer the details to Appendix~\ref{app:l1-laplace}.  Since the multivariate Laplace distribution (for i.i.d variables) depends on the $L_1$-norm of the corresponding vector of variables, we need to solve the following problem: given $R >0$ and $w \notin \mathcal{B}_1(0,R)$, compute
    \begin{align} \label{eq:to-solve-B1-L1}
        \inf_{x \in \mathcal{B}_1(0,R)} \| w- x \|_1.
    \end{align}
The closest element in $\mathcal{B}_1(0,R)$ will lie on the hyperplane defining $\mathcal{B}_1$ that is in the same octant as $w$, so the problem reduces to projecting a point on a hyperplane in $L_1$-distance (the proof in the appendix does not follow this argument but arrives at the same conclusion). Then, we need to compute $h(d,1,f_t,\eta_t L)$. We use a similar approach as in the proof of Theorem~\ref{thm:Linf}, that is, we split the integral and derive a recurrence relation.
\end{proof}

\subsection{Bound for Gaussian noise}

Finally, we derive a bound on the induced leakage when the added noise is Gaussian:

\begin{theorem} \label{thm:l1-gaussian}
If the boundedness assumptions holds for $p=1$ and
$\xi_t \sim \mathcal{N}(0, \sigma^2_t I_d)$, then
\begin{align} \label{eq:thm-l1-gaussian}
     \ml{S}{W} \leq  \sum_{t=1}^T \log \left( \frac{V_1(d,R_t)}{  (2 \pi \sigma^2)^{\frac{d}{2}}}+\frac{(2\eta_t L)^{d-1} (\sigma_t \sqrt{2 d}) }{ (2 \pi \sigma_t^2)^{\frac{d}{2}} ((d-1)!) }   \sum_{i=0}^{d-1} \left( \frac{\sigma_t \sqrt{2 d}}{\eta_t L} \right)^i \Gamma\left(\frac{i+1}{2}\right) \right).
\end{align}
\end{theorem}
~\Cref{thm:l1-gaussian} is tighter than~\Cref{thm:gaussian-ml} for any given $d$ moreover one has, again, that:
\begin{align}
    \lim_{d \rightarrow \infty} \ml{S}{W} = 0. \label{eq:limitToInfinityL1}
\end{align}
\Cref{eq:limitToInfinityL1} follows from~\Cref{thm:gaussian-ml} and the fact that the bound in Proposition~\ref{prop:general} is increasing in $p$ (cf. discussion below~\Cref{thm:gaussian-ml}).

In order to prove~\Cref{thm:l1-gaussian} one has to solve a problem similar to the one introduced in~\Cref{thm:l1-laplace} (cf. equation~\eqref{eq:to-solve-B1-L1}). However, in this case a different norm is involved: i.e., given $R >0$ and $w \notin \mathcal{B}_1(0,R)$, one has to compute
\begin{align}
    \inf_{x \in \mathcal{B}_1(0,R)} \|w-x\|_2.
\end{align}
Again, one can argue that the point achieving the infimum lies on the hyperplane defining $\mathcal{B}_1$ that is in the same octant as $w$. In other words, the minimizer $x^\star$ is such that the sign of each component is the same sign as the corresponding component of $w$ (and lies on the boundary of $\mathcal{B}_1$). Thus, we are simply projecting a point on a hyperplane. The induced integral is solved by an opportune choice of change of variables.  The details of the proof are given in Appendix~\ref{app:l1-gaussian}.

\newpage


\section{Introduction}\label{sec1}

In the current energy crisis, surging prices for fossil fuels add to the debate on alternatives to gas, coal, and oil for energy generation. In addition to ongoing and planned efforts to expand renewable energy sources, many governments have recently declared or emphasized their intention to invest in new nuclear power plants: Japan is considering to invest in new reactors for the first time since the meltdown of the Fukushima Daiichi plant, France recently pledged to construct up to 14 new plants, although specifics remain unclear, and the United Kingdom announced funding for a Sizewell C station to pursue their target of 24 GW nuclear capacity by 2050  \citep{ft, nyt, bb}.

Currently surging prices for fossil fuels only encourage existing plans for nuclear power that pre-date the crisis and originate from government strategies to mitigate carbon emissions and combat climate change \citep{Goldstein2019}. One the one hand, from a technical perspective, nuclear power offers two advantages: First, nuclear plants are dispatchable. While today's plants do not operate flexibly and thus provide baseload supply, research suggests new reactor concepts or even established light-water reactors are capable of flexible operation \citep{Jenkins2018, Lynch2022, mit_future_2018}. As a result, nuclear power could complement fluctuating supply from renewables and ensure security of supply in decarbonized energy systems. Second, nuclear offers a great energy potential. Compared to renewables, nuclear power is far less constrained by land availability or meteorological conditions \citep{landUse}. At the same time, decarbonization in heating, industry, and transport requires electricity as a primary source of energy, either directly or indirectly using synthetic fuels produced from electricity \citep{Luderer2021,Bogdanov2021}. Therefore, nuclear power could substantially contribute to fulfilling the increasing demand for electricity.

On the other hand, there are economic arguments against investing into nuclear power. New plants are capital-intensive; especially since actual construction cost and time frequently exceed forecasts \citep{Lovins2022,Rothwell2022}. For instance, in France, offical construction costs of Flamanville 3 so far quadrupled to 12,600 €/kW while project completion is currently delayed by more than a decade \citep{edf_universal_2022,Rothwell2022}. Similarly, costs of Olkiluoto 3 in Finland tripled over a construction time of 17 years \citep{dw}. Costs of Hinkley Point C in the United Kingdom have doubled to date, and the reactor is currently scheduled to start operation in 2027---11 years after construction started \citep{bbc}. Whether this trend is reversible, is subject to debate: Drawing on the historic expansion of pressurized water reactors in France, \citet{Berthelemy2015} suggest reinforced investments and standardization could lower construction costs and time. \citet{Grubler2010} however states that the expansion of pressurized water reactors never achieved positive learning effects in the first place.

With the high costs of nuclear power, other technologies might be more efficient to achieve decarbonization. Instead of potentially flexible nuclear plants, imports from other regions, batteries, or seasonal storage using synthetic fuels can cover demand when supply from wind and photovoltaic (PV) is low. In addition, electrification in heating, industry, and transport adds a substantial share of flexible demand. For example, electricity consumption of electrolyzers can adjust to wind and solar generation, thereby supporting the integration of renewables \citep{Wang2018}. In the same way, flexible electric heating or charging of electric vehicles can adapt to renewable supply \citep{Schill2020, Schuller2015}. Beyond flexibility, high investment costs can render nuclear power more expensive even compared to renewables in unfavorable sites, like offshore wind turbines in deep waters or PV in locations with little sunshine \citep{mckenna2014}.

Against this background, we investigate the economic efficiency of nuclear power in two steps: First, we review projected and actual costs of new nuclear capacity to compute the conceivable range of levelized costs of energy (LCOE), discuss plausibility and compare the outcome to renewables. Second, we compute the efficient share of nuclear power for previously obtained cost ranges in a decarbonized energy system with a comprehensive energy planning model \citep{Goeke2020a}. Due to its extensive scope, the model can consider all alternatives to nuclear power for the provision of flexibility, most importantly cross-border exchange of energy, flexible demand, and storage technologies. Previous research limited the scope to single regions omitting the exchange of energy as a source of flexibility, neglected the demand-side flexibility induced by electrification, or only featured short-term storage systems \citep{Duan2022,Baik2021}. The cost review focuses on OECD countries due to discrepancy of nuclear development between OECD and non-OECD countries observed in \citet{Rothwell2022}; the techno-economic analysis covers the European continent.

Inevitable simplifications of the techno-economic model are made in a way that favors investment in nuclear power. Therefore, the computed share of nuclear power in the energy mix is not an accurate estimate, but rather an upper bound of what might be if significant nuclear cost reductions were to be achieved. This applies in particular because our techno-economic analysis omits external social costs of nuclear power associated with the risk of accidents, waste management, and proliferation \citep{leveque2014,lordan2021, wheatley_reassessing_2016, lehtveer_nuclear_2015}. Although their robust quantification is difficult, external costs for nuclear typically exceed estimates for renewable energies \citep{stirling_limits_1997}.

The remainder of this paper is structured as follows: The next section \ref{methods} describes the methodology to compute LCOEs and introduces the applied energy planning model. Section \ref{main} presents the results for LCOE analysis and the planning model. In Section \ref{diss} these results are discussed and Section \ref{conclusion} concludes.

\section{Methods} \label{methods}

The first subsection \ref{lcoe} describes the computation of LCOEs for nuclear with an emphasis on the review of construction costs and financing costs during construction. Then, we introduce the energy planning model applied for a more detailed evaluation of nuclear power. The applied model uses the open AnyMOD.jl modeling framework and is described with greater detail in previous publications \citep{Goeke2023,Goeke2020a}.

\subsection{Computation of LCOEs} \label{lcoe}

Whether nuclear power is an economically efficient option for decarbonization strongly depends on its costs. Operational costs for fuel or maintenance only constitute a minor share while construction and financing account for 80\% of the total costs of nuclear power \citep{haas_historical_2019, mackerron_nuclear_1992}. At the same time, costs of construction and financing are also subject to the greatest uncertainty. Therefore, we review construction and financing costs in greater depth and vary them across the observed range in the subsequent analysis. For all other parameters, section \ref{cost_analysis} of the appendix provides a review. Table \ref{tab:assumptions} lists the assumptions used for the analysis. For these parameters, we use optimistic assumptions that correspond to at least the 25th percentile of the reviewed data. All monetary values in this paper are inflation-adjusted and given in US-\$\textsubscript{2018}. Section \ref{lcoeSup} of the appendix provides the corresponding conversion factors.

To assess the range of overnight construction costs, we reviewed different academic publications and industry reports providing 88 individual future projections or reported costs for different reactor concepts and countries \citep{nrel_2021_2021,oecd_full_2018,shirvan_overnight_2022, Rothwell2022, barkatullah_current_2017, lazard2021, wealer_investing_2021, international_atomic_energy_agency_advances_2020, green_smr_2019, rothwell_economics_2016, Duan2022, Grubler2010, stein_advancing_2022, mit_future_2018, boldon_small_2014, stewart_capital_2021, nrel_2021_2021, iea_projected_2020, tolley_economic_2004, dixon_advanced_2017}. We put a focus on light-water reactors (LWR), the most common reactor concept\footnote{Table \ref{tab:OCC_results} in the appendix shows that considering a broader range of reactor types does not give lower cost projections.}, and OECD countries since \citep{Rothwell2022} observes non-OECD estimates are not comparable. Cost projections are often provided in the means of n-th-of-a-kind cost, presenting assumptions of sometimes significant cost reductions through standardization and learning effects, e.g., \citep{shirvan_overnight_2022}. These projections however are not given in relation to a specific year or cumulated investment, complicating the comparison with cost projections for renewables.

The left boxplot in Fig. \ref{fig:OCC_boxplot} summarizes the distribution of the 45 remaining estimates. Overall, projections range from 1,914 US-\$  per kW projected by \citet{oecd_full_2018} to 12,600 US-\$ per kW, the actual costs of the Flamanville-3 project \citep{Rothwell2022}. The median is at 5,430 US-\$ per kW; 50\% of estimates are between 4,328 and 9,966 US-\$ per kW. Future projections are generally much lower than recent observations of actual costs for the European Pressurized Reactor (EPR) and the Westinghouse AP1000. Observations of actual costs are still preliminary since of the 4 relevant projects, Olkiluoto 3 has not yet started commercial operation and Flamanville 3, Hinkley Point C and Voglte Station are still under construction.\footnote{For the Shin Hanul 1 reactor that recently started operation in South Korea, final costs were not publicly disclosed. KEPCO's financial statements suggest costs as low as 3,120 US-\$ per kW \citep{kepco_interim_2022}. This seems very low when taking a construction time of approx. 12 years into account.} Some literature suggests that historically, learning rates of 5 to 10\% were achieved in the nuclear industy which would coincide with the mentioned projections \citep{Berthelemy2015, mignacca_economics_2020, lloyd_expanding_2020}. However, the discrepancy between projected and actual costs is in line with previous analyses on the actual development of nuclear construction costs in Western countries since the 1970s \citep{Grubler2010, davis_prospects_2012, koomey_reactor-level_2007, escobar_rangel_revisiting_2015, harding_economics_2007}, placing doubt on supposed learning rates.
\begin{figure}
    \centering
    \includegraphics[scale=0.1]{figures/boxplotMain.pdf}
    \caption{Boxplots for overnight construction costs and time}
    \label{fig:OCC_boxplot}
\end{figure}

The right boxplot in Fig. \ref{fig:OCC_boxplot} presents a similar image for construction time. Future projections range from 5 to 9 years, but recent and ongoing projects are currently taking 10 to 17 years. These overruns drive up financing costs and can also threaten energy security.
\begin{table}
\caption{Parameters for nuclear power} \label{tab:assumptions}
\centering
\begin{tabular}{p{3cm}|p{2cm}|p{2cm}}
    \hline
    Parameter                           &Unit       &Value / Range\\
    \hline
    Overnight construction cost   &US-\$/kW     &1,914 - 12,600\\
    Annual fixed operational cost       &US-\$/kW     &88.81\\
    Variable operational cost  (incl. fuel)          &US-\$/MWh    &10.96\\
    Capacity factor                     &\%         &95\\
    Construction time             &years      &4 - 10\\
    Depreciation period        &years      &40\\
    \hline
\end{tabular}
\end{table}

Together, construction and financing costs constitute total construction costs (TCC). Financing costs reflect the interest paid and are typically calculated from the length of depreciation and the interest rate. Following the methodology for nuclear power investments introduced in \citep{rothwell_economics_2016}, our analysis additionally accounts for costs during construction. We compute the financing costs during construction $f$ using Equation \ref{idc} introduced by \citet{rothwell_economics_2016}. In the equation $i$ is the interest rate, $t$ the construction time, and $c$ the overnight construction costs.

\begin{equation} \label{idc}
f = (\frac{i}{2 \cdot t} + \frac{i^2}{6 \cdot t^2})*c
\end{equation}

For the comparison of LCOEs and analysis with the energy system model, the annuity $a$ is computed for the sum of financing and overnight construction costs according to Equation \ref{ann}. In the formula, $d$ is the depreciation period.

\begin{equation} \label{ann}
a = (f + c) \cdot \frac{(1+i)^d \cdot i}{(1+i)^d -1}
\end{equation}

Finally, the LCOE $l$ for a single year corresponds to the sum of annuity and variable costs $v$ divided by the yearly generation. The yearly generation corresponds to the product of the utilization rate $u$ and the numbers of hours in a year.

\begin{equation}
l = \frac{a + v \cdot u \cdot 8760}{u \cdot 8760}
\end{equation}

\subsection{Energy system model}

The deployed planning model decides on the resources and technologies to supply an exogenous demand for energy services. Fig. \ref{fig:overMod} provides an overview of the considered primary energy sources, secondary energy carriers, and the categories for final demand. Besides nuclear power, supply of primary energy is restricted to renewable sources, namely wind, hydro, photovoltaic, and biomass. To capture intermediate steps between final demand and primary supply, the model includes secondary carriers, like hydrogen or synthetic gas and oil. Finally, the model must meet the demand for residential and industrial heat, passenger and freight transport, synthetic fuels for aviation and navigation, and electricity.

\begin{figure}
	\centering
		\includegraphics[scale=0.618]{figures/modelOverview.pdf}
	\caption{Overview of model structure}
	\label{fig:overMod}
\end{figure}

To cover final demand, the model decides on the expansion and operation of predefined technologies. Technologies in the model generate, convert, store, or transport energy carriers or provide an energy service. For instance, wind turbines generate electricity that can either satisfy the final electricity demand, be temporarily stored by a battery or be converted to hydrogen by an electrolyzer or to heat by a heat pump to satisfy final demand. Overall, the model includes 22 energy carriers and 125 distinct technologies. Except for nuclear, the choice and parametrization of technologies is based on reports of the \citet{DEA} and on \citet{Robinius2020} for the transport sector. The graph in Fig. \ref{fig:powerAll2} shows all technoloiges supplying electricity. In the graph, vertices either represent energy carriers, depicted as colored squares, or technologies, depicted as gray circles. Entering edges of technologies refer to input carriers; outgoing edges refer to outputs. The appendix \ref{graphSec} provides a comprehensive graphs with all considered technologies.

\begin{figure}
	\centering
		\includegraphics[scale=0.25]{figures/powerAll_wNuclear_cropped.png}
	\caption{Technologies for electricity supply}
	\label{fig:powerAll2}
\end{figure}

The model's objective is to minimize the total system costs which consist of annualized expansion and operational costs for technologies and costs of energy imports from outside the system. Accordingly, the model takes a strictly techno-economic social planner approach and does not consider the different agents within the system, such as producers, consumers, or governments, and their behavior. The model is formulated as a linear optimization problem without constraints on unit commitment. Therefore, operational restrictions, e.g., ramping rates or start-up times, for nuclear power, or any other technology, are neglected. Previous studies agree that operational detail has little impact on results if models include options for short-term flexibility, like batteries or demand-side response \citep{Poncelet2016,Poncelet2020}. In addition, most technologies subject to considerable operational restrictions are excluded from our analysis, for instance coal power plants. As a result, the approach only creates a positive bias in favor of nuclear power, which is in line with our intention to compute an upper bound on its deployment.

The temporal scope of the model consists of a single year. Using a brownfield approach, today's transmission infrastructure is available without further investment due to its long technical lifetime. In addition, the model can decide to invest into additional infrastructure. The same applies to all kinds of hydro power plants. Additional transmission capacities and all other technologies require investment to be available. Technology and demand data is based on projections for the year 2040. The model uses an hourly temporal resolution for electricity to capture fluctuations of renewables and the sensitivity of the power grid to imbalances. For other energy carriers, the model uses a coarser but at least daily resolution to reduce computational load and capture the inherent flexibility \citep{Goeke2020b}. For instance, hydrogen uses a daily resolution to capture its physical inertia and the resulting storage capabilities of pipelines. The resolutions used for other carriers are discussed in the section below.

For timeseries of electricity demand, heat demand, and renewable supply, the model uses historic data from 2008. Since the scope only includes a single climatic year, planning does not consider extreme events with low renewable supply but high demand. On the other hand, we neither consider the risk of coincident nuclear outages threatening security of supply, for instance during the Winter of 2022/2023 in France \citep{nyt2022}.

\subsubsection{Spatial resolution and starting grid} \label{map}

The map in Fig. \ref{fig:startGrid} shows the spatial resolution of the model. The spatial scope includes the European Union plus the United Kingdom, Switzerland, Norway, and the entire Balkan region. The spatial resolution of the model divides this scope into 96 clusters with distinct energy balances and capacities.

\begin{figure}
	\centering
		\includegraphics[scale=0.65]{figures/startGrid.pdf}
	\caption{Pre-existing grid infrastructure (solid) and new potential connections (dashed) in the model}
	\label{fig:startGrid}
\end{figure}

Representation of the electricity grid distinguishes HVAC and HVDC grids and aggregates clusters according to the balancing zones of the European power market. On the one hand, this resolution neglects interzonal congestions, but on the other hand, it allows to build on specific ENTSO-E data on limits and costs for the expansion of net-transfer capacities in the European power grid \citet{entsoe2020}. More detailed resolutions must rely on unspecific cost assumptions that range from 400 to 1,000 EUR per MW and km \citep{Fuersch2013,Gerbaulet2017}. Neglecting interzonal congestions does not impose a bias in favor of renewables and against nuclear power since nuclear capacities are typically centralized in remote areas and require substantial interzonal connection with demand centers. For renewables this is only partially true: While onshore and especially offshore wind depend on transmission grids, too, PV can be built very close to demand. In addition, we use a transport instead of flow formulation for power grid operation since previous research found this simplification to be sufficiently accurate \citep{Neumann2020a}.

To exchange hydrogen between clusters, the model can build transmission infrastructure between neighboring clusters, or overseas in plausible cases. These pipelines are subject to costs of 0.4 million EUR per GW and km and energy losses of 2.44\% per 1,000 km \citep{DEA}. The distance between the geographic center of clusters serves as an estimate for pipeline length.

\subsubsection{Flexibility and electricity demand}

Nuclear power has a key advantage over renewable wind or PV---it is dispatchable. In our analysis, nuclear plants are considered to be fully flexible and are not subject to operational restrictions. To what extent these characteristic drives investment strongly depends on if and how a model includes alternatives to nuclear for the provision of flexibility. One flexibility option, the interconnection of regions, is covered by the large spatial scope and consideration of transmission infrastructure in our model. This section details the representation of two additional options: energy storage and flexible demand.

First, the model includes several technologies for electricity storage. It can invest in batteries or utilize the pre-existing pumped storage and reservoir plants. Second, and more importantly, the model considers investment into storage technologies for other carriers, too. Tanks and caverns for hydrogen storage decouple electricity consumption by electrolyzers from hydrogen demand adding substantial flexibility since the modeled system requires a great amount of hydrogen for high temperature process heat and synthetic aviation fuels. In addition, using stored hydrogen for electricity generation provides an option for long-term storage that batteries or pumped storage cannot. Finally, the model includes different technologies for heat storage. At additional investment costs, electric heating systems on the consumer level can be equipped with short-term heat storages to mitigate load peaks. In heating networks storage options include water tanks for short and thermal pits for long storage durations.

Beyond storage, the model captures to what extent final demand is elastic and can adapt to supply. To represent this, it varies temporal resolution by energy carrier. For the final demand for electricity, we assume no elasticity and use an hourly resolution.

The demand for residential and industrial heat uses a 4-hour resolution, but the generation profile of each individual heating system on the consumer level is fixed to match demand. Technologies within heating networks can instead operate flexible but must jointly meet the final demand of the connected consumers, similar to the power grid. In residential heating, the 4-hour resolution captures the thermal inertia of buildings \citep{Heinen2017}; in heating networks the inertia of the network  itself \citep{Triebs2022}; and in industrial heating the possibility to reschedule processes. It should be noted that despite the 4-hour resolution, industrial heating is mostly inflexible in this approach, because demand is almost constant, and systems must continuously operate close to full capacity.

Demand from BEVs is flexible within limits: First, an hourly profile restricts charging and reflects when BEVs are connected to the grid. Second, for each day the sum of charged electricity must match consumption. To implement the latter, transport services use a daily resolution in the model. This approach is a compromise regarding the flexibility of BEVs in future power systems which is still subject to significant uncertainty. On the one hand, we assume charging can be flexible to some degree and can adapt to supply \citep{Strobel2022}. On the other, we do not consider bidirectional loading, also known as vehicle-to-grid, where BEVs can operate like regular batteries and supply energy back to the grid \citep{Hannan2022}. Apart from BEVs, other electricity consumption in the transport sector, for instance for rail transport, is inflexible and uses an hourly resolution.

\subsubsection{Renewable supply}

In the model, the level of investment into nuclear power depends on the potential and costs of its alternative renewable energy. Therefore, this section provides further details on the used assumptions regarding renewable supply.

Assumptions on the technical potential of wind and solar can vary substantially. In our study, we use data from \citet{Auer2020} for the potential of wind and solar. To put these numbers into context, we compare them against results of a meta-analysis on renewable potential in Germany in the appendix \ref{supEE} \citep{Risch2022}. The material demonstrates that the assumed potential is in line with results of the meta-analysis, except for PV, for which our assumptions are more conservative.

To reflect that sites for renewables differ in quality, wind and PV technologies are further subdivided. Onshore wind and openspace PV each have three sub-technologies; rooftop PV has two. Energy yield and capacity potential for each sub-technology are computed for each of the 96 clusters based on the respective distribution of full-load hours. Offshore wind is subdivided into two technologies with identical full-load hours, but different investment costs representing shallow and deep waters. Figures in section \ref{supEE} of the appendix sort the entire potential of different renewables by full-load hours and demonstrate how increased deployment results in smaller yields.

The use of biomass in each country is subject to an upper energy limit that sums to 1,081 TWh for the entire model \citep{jrc2015}. In addition to domestic production, the model can import renewable hydrogen by ship at costs of 131.8 US-\$\textsubscript{2018} per MWh and by pipeline from Marocco or Egypt at 90.7 and 86.8 US-\$\textsubscript{2018} per MWh, respectively \citep{Hampp2021}.

\section{Results} \label{main}

This section \ref{main} first compares the computed range of LCOEs for nuclear power with renewable sources. Afterwards, we present the results of the more elaborate analysis applying the energy planning model.

\subsection{Levelized cost comparison}

Based on the review of construction projections, we compute a conceivable range for LCOE for nuclear power using the formula presented in section \ref{lcoe}. Fig. \ref{fig:lcoe} compares this range against historic LCOEs on the left side (a) and against cost projections for renewable energy on the right (b). Historic LCOEs for all energy sources build on the corresponding reports by \citeauthor{lazard2021} from 2009 to 2021.\footnote{To exclude subsidies and harmonize interest rates, we re-computed LCOEs based on the reported techno-economic assumptions instead of directly copying the values in the Lazard reports \citep{lazard2009,lazard2010,lazard2011,lazard2012,lazard2013,lazard2014,lazard2015,lazard2016,lazard2017,lazard2018,lazard2019,lazard2020,lazard2021}.} The projection for nuclear results from the parameters in Table \ref{tab:assumptions}, projections for renewable sources use 2040 data from \citet{DEA} for cost assumptions and the 10th to 90th percentile of respective energy yield in Europe for full-load hours. For a full documentation of assumptions on costs and renewable potential, see the respective sections in the appendix \ref{cost_analysis} and \ref{lcoeSup}.

\begin{figure}[!htbp]
	\centering
		\includegraphics[scale=0.1]{figures/lcoeComparison.pdf}
	\caption{Comparison of LCOEs for generation technologies}
	\label{fig:lcoe}
\end{figure}

The comparison of historic LCOEs shows how, in the past, costs for nuclear power increased but dropped for renewable energy, in particular for PV. Current costs for nuclear power reported by Lazard are at the upper end of our projections for the future. For wind, current costs are already within the project range; for PV, current costs are above projections implying a continuous degression of costs. Overall, the comparison of projected LCOEs suggests that nuclear can only compete with renewables if the industry can reverse the historic trend and achieve substantial learning effects and economies of scale to reduce project costs.

\subsection{Planning model analysis}

However, comparing LCOEs for nuclear and renewables is of limited insight because nuclear and renewables ultimately provide different goods: nuclear plants are dispatchable and can, within limits, adapt supply to demand; renewable energy sources like wind and PV cannot. Therefore, we elaborate our comparison and deploy a techno-economic planning model of the European energy system for further analysis. The model captures fluctuations of wind and PV and assumes fully flexible operation of nuclear plants. Its objective is to minimize the sum of investment and operational costs for meeting energy demand in a fully decarbonized system.

Fig. \ref{fig:costShare} presents the key results of the energy system model analysis. On the left axis, the graph shows the share of nuclear power in total electricity generation depending on the overnight construction costs and time. The right axis provides the corresponding distribution of cost projections in the literature and actual costs reported by recent projects.

\begin{figure}[!htbp]
	\centering
		\includegraphics[scale=0.1]{figures/costShare.pdf}
	\caption{Generation share of nuclear and distribution of overnight cost projections}
	\label{fig:costShare}
\end{figure}

The results confirm the tentative conclusion from the LCOE analysis. Only if construction costs and time drop substantially, nuclear power provides a significant share of electricity. At 7 years construction time, overnight costs must drop to 4,120 US-\$ per kW to achieve a 10\% share in generation, at 4 years construction time this threshold increases to 4,450 US-\$ per kW. At the same time, only 20\% of overnight construction cost projections are below 4,120 US-\$ per kW and actual costs of recent projects lie above 8,000 US-\$ per kW. At the lowest construction costs, projected by \citet{oecd_full_2018}, nuclear power provides up to 52 \% of total electricity generation. Non-nuclear electricity generation largely corresponds to wind and PV. Other sources are hydro which is kept at today's level supplying 540 TWh. Biomass supplies at most 13.6 TWh since the model predominantly deploys the limited potential to cover the demand for synthetic fuels instead.

Despite its assumed flexibility, nuclear power rather substitutes than complements fluctuating supply from wind and PV generation. In the most extreme case with a 52 \% generation share, nuclear power plants are operating at a utilization rate of 94.4\%; close to the 95\% limit. This changes only slightly at smaller nuclear and consequently higher renewable shares. At a nuclear share of 7.9\%, utilization still amounts to 88.7\%, which is above the average of European nuclear plants operating today \citep{wealer_ten_2021}. This implies that even if technically capable, nuclear power is economically not well suited to provide system flexibility due to its high investment costs.

To highlight the system's renewable integration, Fig. \ref{fig:durCur} exemplarily compares the residual load curves of Germany for a system with the highest nuclear share of 52\% and no nuclear capacity at all. Residual load is the remaining electricity demand after subtracting supply of fluctuating renewables from total demand. Plotted in descending order, the curve shows the residual peak load on the y-axis, the amount of residual demand as the area above the x-axis, and excess energy as the area below the x-axis. Appendix \ref{timeSeries} shows an extract of the hourly model results from which the duration curves were computed.

\begin{figure}[!htbp]
	\centering
		\includegraphics[scale=0.1]{figures/durationCurve.pdf}
	\caption{Residual load curve for Germany in a system with a nuclear share of 52\% (left) and no nuclear at all (right)}
	\label{fig:durCur}
\end{figure}

In the graphs, there are two residual load curves: 'Flexible demand' shows the actual residual load in the model; 'inflexible demand' is computed ex-post assuming demand cannot adapt to supply. To what extent demand is flexible depends on its use. The conventional final demand for electricity is inflexible, but only constitutes around 39\% of total demand, exact shares differing from case to case. Additionally, the creation of synthetic fuels takes up 24\%, electrified transportation 20\% and heating 17\% of total electricity demand. This demand can be flexible within specific restrictions detailed in Section \ref{methods}. Some of this flexibility comes without costs, like flexible charging of battery electric vehicles, but most require the model to invest, for instance into heat or hydrogen storage.

In the 52\% nuclear system, wind and PV only account for 10\% of generation in Germany and the need for flexibility is small. The residual load is flat and does not substantially adapt to supply. For the system without nuclear, results are contrary. To reduce the amount of residual demand and excess energy, demand substantially adapts to supply. For instance, electrolyzers operate more flexibly, thereby causing their full-load hours to decrease to 3,343 compared to 5,979 in the 52\% nuclear system while more than doubling their capacity from 35.9 to 77.1 GW. Investments in heat and hydrogen storage provide additional flexibility. Storage capacities for district heating increase to 40.9 TWh compared to 18.2 Twh in the 52\% nuclear system; for hydrogen even to 12.4 from 1.7 TWh.

Coloured areas in the graph indicate how the system covers the residual load. In the 52\% nuclear system, 132 GW of nuclear operate at full capacity and cover the predominant share of residual load. Additional imports of 86.9 TWh from neighboring countries cover peak loads. Vice versa, Germany exports 78.5 TWh of excess supply. In the system without nuclear, the residual load is smaller, but imports are higher and cover 244 TWh of demand---the predominant share of residual load. Dispatchable thermal plants fueled by hydrogen and biogas only provide 7.3 TWh of generation. In both cases, there is no additional investment in electricity storage and both systems only deploy the existing pumped storage and hydro reservoirs.

Although the integration of renewables requires additional investments, nuclear power is only a cost-efficient alternative at very low construction costs as shown previously. Otherwise, investment costs of nuclear significantly outweigh the costs of renewables and corresponding flexibility options. For instance, the system without nuclear shows 9.7 billion US-\$ of additional annualized investment into demand side flexibility compared to the 52\% nuclear system. Adding 8.2 billion US-\$ for additional grid infrastructure and 3.0 billion US-\$ for thermal plants, total costs of renewable integration add up to 20.9 billion US-\$. But on the other hand, even if overnight construction costs of nuclear drop to 5,430 US-\$ per kW, the median of projections in Fig. \ref{fig:costShare}, the 52\% nuclear system is already 46.9 billion US-\$ more expensive than the system without nuclear. This 46.9 billion US-\$ difference is composed of a 67.8 billion US-\$ increase in generation costs only partially offset by 20.9 billion US-\$ of additional integration costs in the system without nuclear exclusively deploying renewables. For reference, the total costs of renewable generation amount to 150.5 billion US-\$ in this case.

Fig. \ref{fig:mapRes} further investigates the role of import and export for both systems. It shows electricity generation and net exchange of hydrogen and electricity by high-voltage alternating current (HVAC) or direct current (HVDC) grids for each country. Exchange of hydrogen and electricity increases substantially in the system without nuclear reflecting the 8.2 billion US-\$ rise of costs for grid infrastructure. The maps show that this infrastructure not only serves renewable integration but reduces generation costs as well because it allows to utilize renewables where full-load hours are highest, for instance in southern Europe for PV or in the North Sea for offshore wind.

\begin{figure}[!htbp]
	\centering
		\includegraphics[scale=0.65]{figures/mapResults.pdf}
	\caption{Electricity generation and exchange the 52\% nuclear system (top) and the system without nuclear (bottom)}
	\label{fig:mapRes}
\end{figure}

For a comprehensive overview of energy flows and conversion processes in the compared energy systems, see appendix \ref{sankey}.

\section{Discussion} \label{diss}

Although building on an extensive review of cost projections and a detailed planning model, our analysis faces some inevitable limitations. Overall, these limitations impose a positive bias on nuclear power, so results are rather an upper bound than reliable estimates for the cost-efficient share of nuclear.

First, there are simplifications of the techno-economic model itself. We consider construction cost of nuclear (and all other technologies) as exogenous, although costs can depend on the level of investment due to learning effects. Therefore, cases in Fig. \ref{fig:costShare} with low nuclear shares, but construction costs dropping below 5,000 US-\$ per kW, are implausible. If learning can in fact achieve such a drop in costs, it will require capacity investments that will already raise the share of nuclear in power generation above 10\%. To achieve this share, 41 GW of nuclear would have to be built which approximately corresponds to 25 EPR-sized reactors. A nuclear share of 50\% requires 132.5 GW of capacity, i.e. 83 EPR-sized reactors. Given the current state of ongoing projects in France, Finland and the UK, this seems rather challenging \citep{Rothwell2022}. Considering learning effects endogenously is not an option, due to the uncertainty of their existence and computational limits \citep{Berthelemy2015,Grubler2010, harding_economics_2007}.

The model focuses on the single year 2040, but decarbonization of the energy system is a process over several decades. Computational limits however require to limit the analysis to one year. On the one hand, this year should be sufficiently far in the future to allow for learning effects of nuclear. On the other hand, years further into the future will increasingly favor renewables since learning effects are expected to keep reducing their costs, too \citep{DEA}. Against this background, we choose 2040 as the earliest date until that substantial reductions of costs are conceivable, an optimistic assumption considering the construction times of recent European nuclear projects.

On the more technical side, we omit unit commitment constraints, as already mentioned. Of all considered technologies, these operational restrictions are most pronounced for nuclear power. The method section discusses further technical limitations like the representation of extreme events and grid constraints.

Second, some adverse characteristics of nuclear power escape the techno-economic perspective of the analysis. For instance, costs of managing nuclear waste and decommissioning old plants are substantial but not included in our analysis due to unavailability of reliable data and the long time frame \citep{lordan2021, Lovins2022}. Furthermore, nuclear power carries the risk of severe accidents, proliferation of nuclear material and radioactive contamination \citep{leveque2014, wheatley_reassessing_2016, lehtveer_nuclear_2015}. Although robust quantification of these externalities is difficult and therefore not part of the analysis, their social costs typically exceed estimates for renewable energy \citep{stirling_limits_1997}.

\section{Conclusion} \label{conclusion}

Our least-cost analysis indicates that nuclear power is only a cost efficient option for full decarbonization of Europe's energy system if construction costs at least half and drop below 4,000 US-\$\textsubscript{2018} per kW. In contrast, costs and construction time for nuclear power in the last ten years have been increasing rather than decreasing. This assessment uses optimistic assumptions for other characteristics of nuclear, like availability and operational costs, and still omits social costs, like the risk of accidents or waste management, since these externalities escape robust quantification.

Results show that nuclear power reduces investments into system flexibility. Since nuclear plants are, in contrast to wind or PV, dispatchable, their system integration induces lower costs for grid infrastructure, energy storage, demand-side flexibility, or thermal backup plants. However, if generation costs of nuclear do not drop substantially, they outweigh this benefit. For instance, assuming overnight construction costs of 5,430 US-\$ per kW, the median of compiled projections, a 52\% generation share of nuclear will reduce investments into system flexibility by 20.9 bn US-\$ per year but increase generation costs by 67.8 bn US-\$ or 44\% compared to a fully renewable system.

A crucial detail to understand these results is that apart from grid infrastructure, most investments into system flexibility occur outside of the power sector. Instead of battery storage or conventional demand-side management, new electricity demand induced by the decarbonization of heating, transport, and industry adds substantial flexibility. For instance, hydrogen production, accounting for 25\% of electricity demand, supports renewable integration by adapting to renewable supply and using storage to match production with demand.

Results do not indicate that fluctuating renewables and flexibly operated nuclear plants are economic complements. Even when neglecting operational restrictions, nuclear power runs at utilization rates close to 90\%. From a cost perspective, it is preferable to invest into system flexibility rather than running capital-intensive nuclear plants at low utilization rates.

Our analysis puts previous least-cost analyses on nuclear power that were limited to the power sector and a single country into perspective. As we showed, this setup omits two key options for system integration of fluctuating renewables: First, grid infrastructure connecting regions with different generation profiles, and second, the potential flexibility of electricity demand added to decarbonize the heating, transport, and industry sector. As a result, previous analyses overestimated the value of dispatchable nuclear power plants. Not considering other types of firm generation, like hydrogen or biogas fueled plants, adds to this bias.

Our results suggest that efforts to advance nuclear technology---if undertaken at all---should focus on reducing costs and construction time. If plans to extend nuclear capacity are followed through, and renewable investments consequently neglected, project construction delays might put European energy security at risk. Concepts for more flexible operation are not promising because even fully flexible plants almost exclusively operate at full capacity for economic reasons. Since the share of nuclear also affects investments into grid infrastructure and energy storage, path dependencies in system development arise. The risk of nuclear power not achieving substantial cost reductions should thus be accounted for in long-term planning. Potential non-electrical applications for nuclear power, such as high temperature provision for industrial processes, are not considered in this paper but should be addressed in future research.

\section*{Supplementary material}

All data and running scripts used in this paper are openly available here: \url{https://github.com/leonardgoeke/EuSysMod/tree/greenfield_nu} Due to their size, model outputs are not part of the repository but are available upon request.

The applied model uses the open AnyMOD.jl modeling framework \citep{Goeke2020a}. The applied version of the tool is openly available under this link. \url{https://github.com/leonardgoeke/AnyMOD.jl/releases/tag/flexibleElectrificationWorkingPaper}

\section{Acknowledgments}

The research leading to these results has received funding from the Deutsche Forschungsgemeinschaft (DFG) under project number 423336886.

\printcredits

\bibliographystyle{elsarticle-num-names}