Antoine Jacquier

Alexandre Pannier

LARGE AND MODERATE DEVIATIONS FOR STOCHASTIC VOLTERRA SYSTEMS

Keywords: Mathematics Subject Classification. 60F10, 60G22, 91G20 stochastic Volterra equations, large deviations, moderate deviations, weak convergence, rough volatility. 1

. We show in particular how this framework encompasses most rough volatility models used in mathematical finance, yields pathwise moderate deviations for the first time and generalises many recent results in the literature.

INTRODUCTION

This paper sheds new light on the asymptotic behaviour of the class of stochastic Volterra equations (SVEs)

(1.1)

X t = X 0 + t 0 K(t, s)b(s, X s )ds + t 0 K(t, s)σ(s, X s )dW s , t ∈ [0, T ],
for some fixed time horizon T > 0, where X 0 ∈ R d , d ≥ 1, W is a multidimensional Brownian motion, K is a kernel that may be singular, and the coefficients are such that a unique pathwise solution exists. This class of models has been investigated in many fields, including nonlinear filtering [START_REF] Coutin | Abstract nonlinear filtering theory in the presence of fractional Brownian motion[END_REF] using fractional Brownian motion kernels, pharmacokinetic models [START_REF] Marie | A pathwise fractional one compartment intra-veinous bolus model[END_REF] (Langevin equation driven by fractional Brownian motion), fluid turbulence [START_REF] Chevillard | Regularized fractional Ornstein-Uhlenbeck processes and their relevance to the modeling of fluid turbulence[END_REF], and turbulence modelling in atmospheric winds or energy prices [START_REF] Barndorff-Nielsen | Assessing relative volatility/intermittency/energy dissipation[END_REF][START_REF] Corcuera | Asymptotic theory for Brownian semi-stationary processes with application to turbulence[END_REF] using Brownian semistationary processes. Mathematical finance has however been the most dynamic area by far in terms of applications of SVEs, and an in-depth study of (1.1) in the affine case with convolution kernels was recently carried out by Abi Jaber, Larsson and Pulido [START_REF] Jaber | Affine Volterra processes[END_REF]. Following previous analyses supporting non-Markovian systems [START_REF] Alòs | On the short-time behavior of the implied volatility for jump-diffusion models with stochastic volatility[END_REF][START_REF] Chronopoulou | Estimation and pricing under long-memory stochastic volatility[END_REF][START_REF] Chronopoulou | Stochastic volatility and option pricing with long-memory in discrete and continuous time[END_REF][START_REF] Comte | Long memory in continuous-time stochastic volatility models[END_REF][START_REF] Comte | Affine fractional stochastic volatility models with application to option pricing[END_REF][START_REF] Fukasawa | Asymptotic analysis for stochastic volatility: martingale expansion[END_REF], the investigation of high-frequency data in [START_REF] Gatheral | Volatility is rough[END_REF] revealed the roughness, in the sense of low Hölder regularity, of the observed time series of the instantaneous volatility of stock price processes. This suggested that fractional Brownian motion (fBm) with small Hurst parameter (H ≈ 0.1) is an accurate driver for its dynamics. Since this seminal observation, more advanced results [START_REF] Euch | The microstructural foundations of leverage effect and rough volatility[END_REF] have proposed that the drift and the diffusion coefficients should be state dependent, giving rise to the widespread development of (1.1) in quantitative finance.

For option pricing purposes, the asymptotic results in [START_REF] Alòs | On the short-time behavior of the implied volatility for jump-diffusion models with stochastic volatility[END_REF][START_REF] Bayer | Short-time near-the-money skew in rough fractional volatility models[END_REF][START_REF] Fukasawa | Short-time at-the-money skew and rough fractional volatility[END_REF] showed that the short-maturity behaviour of option prices is captured much more accurately by these rough volatility models rather than by Markovian diffusions. Reconciling the stylised facts of the markets from both the statistical and the option pricing viewpoints is the tour de force that make these models so important today. However, the loss in tractability compared to classical Itô diffusions is not negligible. The solution to (1.1) is in general not a semimartingale nor a Markov process, preventing the use of Itô calculus or Feynman-Kac type formulas. Path-dependent versions of the latter are available in some cases, in particular for affine rough volatility models [START_REF] Jaber | Affine Volterra processes[END_REF][START_REF] Cuchiero | Generalized Feller processes and Markovian lifts of stochastic Volterra processes: the affine case[END_REF][START_REF] Euch | The characteristic function of rough Heston models[END_REF][START_REF] Gatheral | Affine forward variance models[END_REF], but general results are scarce [START_REF] Viens | A martingale approach for fractional Brownian motions and related path dependent PDEs[END_REF]. Rough path theory is a natural route but is not available for H ≤ 1/4, although a regularity structure approach was recently developed [START_REF] Bayer | A regularity structure for rough volatility[END_REF]. In this context, one could turn to numerical methods to understand the dynamics or to price options but, despite new advances based on Monte-Carlo methods [START_REF] Bayer | Pricing under rough volatility[END_REF][START_REF] Bennedsen | Hybrid scheme for Brownian semistationary processes[END_REF][START_REF] Mccrickerd | Turbocharging Monte-Carlo pricing for the rough Bergomi model[END_REF], rough Donsker theorem [START_REF] Horvath | Functional central limit theorems for rough volatility[END_REF] or Fourier methods [START_REF] Euch | The characteristic function of rough Heston models[END_REF], the roughness and memory of the process seriously complicate the task.

Asymptotic methods have been used, both to provide clearer understanding of models in extreme parameter configurations and to act as proxies to numerical schemes. Large Deviations Principles (LDP), in particular, have been widely explored in mathematical finance, and we refer the interested reader to [START_REF]Large Deviations and Asymptotic Methods in Finance[END_REF] for an overview. Let {X ε } ε>0 be a sequence of random variables in some Polish space X , converging in probability to a deterministic limit X as ε goes to zero. This sequence is said to satisfy an LDP with speed ε -1 and rate function I : X → [0, +∞] if for all Borel subsets B ⊂ X , the inequalities

-inf x∈B • I(x) ≤ lim inf ε↓0 ε log P X ε ∈ B ≤ lim sup ε↓0 ε log P X ε ∈ B ≤ -inf x∈B I(x)
hold, and the level sets {x ∈ X : I(x) ≤ N } of I are compact for all N > 0. This rate function encompasses in a (relatively) concise formula first-order information about the asymptotic behaviour of complex dynamical systems. If X satisfies (1.1) and X = R d , one can consider finite-dimensional LDP for {X t } t≥0 (also called small-time LDP if the limit takes place as t goes to zero), or pathwise LDP for some rescaling of X with X = C([0, T ] : R d ). The former is easily recovered from the latter by a projection argument. Moderate deviations however are concerned with deviations of a lower order than large deviations, and thus apply to 'less rare events'. We indeed say that {X ε } ε>0 satisfies a moderate deviations principle (MDP) if {η ε } ε>0 satisfies an LDP with speed h 2 ε , where

η ε := X ε -X √ εh ε , for all ε > 0,
with lim ε↓0 h ε = +∞ and lim ε↓0 √ εh ε = 0. Since the speed of convergence of h ε is not fixed, an MDP essentially bridges the gap between the central limit regime where h ε = 1 and the LDP regime where h ε = ε -1/2 . An edifying example of the relevance of moderate deviations appears in [START_REF] Friz | Option pricing in the moderate deviations regime[END_REF] where, interested in option pricing asymptotics, the authors judiciously rescale the strikes with respect to time to expiry. Indeed, as time to expiry becomes smaller, the range of pertinent strikes naturally shrinks, and this 'moderately-out-of-the-money' regime becomes more realistic.

Large deviations for SVEs were originally studied in [START_REF] Nualart | Large deviations for stochastic Volterra equations[END_REF][START_REF] Rovira | Large deviations for stochastic Volterra equations in the plane[END_REF] with regular kernels. In the context of rough volatility, Forde and Zhang [START_REF] Forde | Asymptotics for rough stochastic volatility models[END_REF] introduced the first finite-dimensional LDP where the log-volatility is modelled by a fractional Brownian motion, and refined versions followed in [START_REF] Bayer | A regularity structure for rough volatility[END_REF][START_REF] Bayer | Short-time near-the-money skew in rough fractional volatility models[END_REF][START_REF] Friz | Precise asymptotics: robust stochastic volatility models[END_REF], while pathwise LDP for similar models were studied in [START_REF] Cellupica | Pathwise asymptotics for Volterra type rough volatility models[END_REF][START_REF] Gulisashvili | Large deviation principle for Volterra type fractional stochastic volatility models[END_REF]. Departing from regular conditions on the behaviour of the coefficients led to specific requirements, and finite-dimensional large deviations for the fractional Heston model were carried out in [START_REF] Forde | Small-time, large-time, and H → 0 asymptotics for the rough Heston model[END_REF][START_REF] Guennon | Asymptotic behavior of the fractional Heston model[END_REF], while more elaborate pathwise LDPs were derived for the rough Stein-Stein model with random starting point [START_REF] Horvath | Asymptotic behaviour of randomised fractional volatility models[END_REF], for the rough Bergomi model [START_REF] Jacquier | Pathwise large deviations for the rough Bergomi model[END_REF], and small-time LDPs for the multi-factor rough Bergomi appeared in [START_REF] Lacombe | Asymptotics for volatility derivatives in multi-factor rough volatility models[END_REF]. We emphasise at this point that no pathwise MDP was previously known in the context of rough volatility.

The Gärtner-Ellis theorem [START_REF] Dembo | Large Deviations Techniques and Applications[END_REF]Theorem 2.3.6] is the main ingredient of a finite-dimensional LDP and depends on explicit computations of certain limits of the Laplace transform. This is only available though, when the process is either Gaussian [START_REF] Forde | Asymptotics for rough stochastic volatility models[END_REF] or affine [START_REF] Forde | Small-time, large-time, and H → 0 asymptotics for the rough Heston model[END_REF]. Pathwise LDP on the other hand, have mainly been derived using the Freidlin-Wentzell approach [START_REF] Freidlin | Random Perturbations of Dynamical Systems[END_REF]: starting from known large deviations for the driving (Gaussian) process [START_REF] Deuschel | Large Deviations[END_REF]Theorem 3.4.5], they follow from a combination of approximations and continuous mapping, keeping track of the rate function. While this methodology is clear, it requires a case-by-case tailored path for each model, and in general leads to a cumbersome rate function. Furthermore, pathwise moderate deviations are so far out of reach in this approach, partially explaining the small number of related results compared to LDP.

A radically different method, introduced by Dupuis and Ellis in the monograph [START_REF] Dupuis | A Weak Convergence Approach to the Theory of Large Deviations[END_REF] and developed further by Budhiraja and Dupuis [START_REF] Budhiraja | Analysis and Approximation of Rare Events Representations and Weak Convergence Methods[END_REF], relies on the equivalence between the LDP and the Laplace principle. The family {X ε } ε>0 is said to satisfy the Laplace principle with speed ε -1 and rate function I : X → [0, +∞] if for all continuous bounded maps F : X → R,

(1.2) lim ε↓0 -ε log E exp - F (X ε ) ε = inf x∈X I(x) + F (x) .
This alternative, called the weak convergence approach, consists in proving a Laplace principle where the left-hand side pre-limit of (1.2) can be represented as a variational principle for expectations of functionals of Brownian motion [START_REF] Boué | A variational representation for certain functionals of Brownian motion[END_REF]Theorem 3.1]:

Lemma 1.1 (Boué-Dupuis). Let W be an R m -Brownian motion and F be a bounded Borel-measurable function mapping C([0, T ] : R m ) into R. Then

(1.3) -log E e -F (W ) = inf v∈A E 1 2 T 0 |v s | 2 ds + F W + • 0 v s ds ,
where

(1.4) A := v : Ω × [0, T ] → R m progressively measurable, E T 0 |v t | 2 dt < +∞ .
The representation (1.3) contains in a single formula the usual tools used in the proof of an LDP. The first term on the right-hand side comes from the relative entropy between the Wiener measure and the measure shifted by • 0 v s ds via Girsanov's theorem, under which W + • 0 v s ds is a Brownian motion. It can be interpreted as the cost of deviating from the original path and clearly indicates where the form of the rate function comes from. In essence, this representation replaces the nonlinear analysis of the Freidlin-Wentzell approach with the linear theory of weak convergence. Instead of exponential estimates, only qualitative properties of the shifted process need to be established, such as strong existence and uniqueness and tightness.

The extensive literature on the topic, summarised in [START_REF] Budhiraja | Analysis and Approximation of Rare Events Representations and Weak Convergence Methods[END_REF] and the references therein, demonstrates the strength of this generic approach which can be applied to a variety of models without appealing to their particular features. It has been used to derive LDPs, in the continuous-time case, for diffusions [START_REF] Chiarini | On large deviations for small noise Itô processes[END_REF], multiscale systems [START_REF] Budhiraja | Large deviations for small noise diffusions in a fast Markovian environment[END_REF][START_REF] Dupuis | Large deviations for multiscale diffusions via weak convergence methods[END_REF][START_REF] Spiliopoulos | Large deviations and importance sampling for systems of slow-fast motion[END_REF], SVEs with singular kernels and Lipschitz continuous coefficients [START_REF] Zhang | Euler schemes and large deviations for stochastic Volterra equations with singular kernels[END_REF], SDEs driven by infinite-dimensional Brownian motions [START_REF] Budhiraja | Large deviations for infinite dimensional stochastic dynamical systems[END_REF], by Poisson random measures [START_REF] Budhiraja | Large deviations for stochastic partial differential equations driven by a Poisson random measure[END_REF] or both [START_REF] Budhiraja | Variational representations for continuous time processes[END_REF], including stochastic PDEs. Contrary to the Freidlin-Wentzell approach, this method has also proved efficient to obtain MDPs. SVEs with Lipschitz continuous kernels [START_REF] Li | A moderate deviation principle for stochastic Volterra equation[END_REF], SDEs with jumps [START_REF] Budhiraja | Moderate deviations principles for stochastic differential equations with jumps[END_REF] and slow-fast systems [START_REF] Morse | Moderate deviations principle for systems of slow-fast diffusions[END_REF] are a few relevant examples. The latter were then tailored to the setting of stochastic volatility models in [START_REF] Jacquier | Pathwise moderate deviations for option pricing[END_REF], developing the first application of the weak convergence approach in mathematical finance, and extending the MDP results in [START_REF] Friz | Option pricing in the moderate deviations regime[END_REF] to a pathwise setting. A further appealing feature of moderate deviations is the simple form, often quadratic, of the rate function, as opposed to that provided by large deviations, thereby opening the gates to the use of importance sampling and variance reduction techniques [START_REF] Dupuis | Importance sampling for multiscale diffusions[END_REF][START_REF] Robertson | Sample path large deviations and optimal importance sampling for stochastic volatility models[END_REF][START_REF] Morse | Importance sampling for slow-fast diffusions based on moderate deviations[END_REF].

Building on this powerful approach, we provide a unified treatment of (finite-dimensional and pathwise) large and moderate deviations in the general framework (1.1) by showing the weak convergence of a perturbed system. We relax the uniqueness requirement for the limiting Volterra equation, as in [START_REF] Conforti | On small-noise equations with degenerate limiting system arising from volatility models[END_REF][START_REF] Donati-Martin | Large deviations for squares of Bessel and Ornstein-Uhlenbeck processes[END_REF] for the diffusion case, allowing us to consider coefficients that are not Lipschitz continuous and do not necessarily have sublinear growth.

The paper is organised in the following way: Section 2 introduces the framework and useful definitions. In section 3, we present abstract criteria for the validity of an LDP, extending the results by Budhiraja and Dupuis [START_REF] Budhiraja | A variational representation for positive functionals of infinite dimensional Brownian motion[END_REF]Theorem 4.1]. Our main results, Theorem 3.11 for LDP and Theorem 3.20 for MDP, are then stated in the case of convolution kernels and extended to nonconvolution kernels in Theorem 3.29. In Section 4, we show how these results apply to rough volatility models, and give precise formulae for the rough Stein-Stein, the (multi-factor) rough Bergomi and the rough Heston models. We finally gather technical proofs in the appendix.

d 1 ≥ 1, d 2 ≥ 1, M > 0, f : R d1 → R d2 , we write f M := sup{|f (x)| : |x| ≤ M }.
Unless stated otherwise, constants will be denoted by C (with possible subscript) and may be different from one proof to another. Every statement involving ε stands for all ε > 0 small enough. A family of random variables will be called tight if the corresponding measures are tight [START_REF] Dupuis | A Weak Convergence Approach to the Theory of Large Deviations[END_REF]Appendix A]. We also use the classical convention that the infimum over an empty set is equal to infinity. Finally, we recall the following definitions for clarity and notations:

Definition 2.1. Let g be a function from R d to R n .
• It has linear growth if there exists

C L > 0 such that |g(x)| ≤ C L (1 + |x|), for all x ∈ R d ;
• if it is uniformly continuous, it admits a continuous and increasing modulus of continuity ρ g : R + → R + , with ρ g (0) = 0 and |g(x) -g(y)| ≤ ρ g (|x -y|), for all x, y ∈ R d ; • it is locally δ-Hölder continuous with δ ∈ (0, 1) if, for all M > 0, there exists

C M > 0 such that |g(x) -g(y)| ≤ C M |x -y| δ , for all |x| ∨ |y| ≤ M .
2.2. Framework. We consider small-noise convolution stochastic Volterra equations (SVE)

(2.1)

X ε t = X ε 0 + t 0 K(t -s)b ε (s, X ε s )ds + ϑ ε t 0 K(t -s)σ ε (s, X ε s )dW s , t ∈ T,
taking values in R d with d ≥ 1, where ε > 0, and ϑ ε > 0 tends to zero as ε goes to zero. For each

ε > 0, X ε 0 ∈ R d , b ε : T × R d → R d and σ ε : T × R d → R d×m
are Borel-measurable functions, and W is an m-dimensional Brownian motion on the filtered probability space (Ω, F , {F t } t∈T , P) satisfying the usual conditions. The kernel function K : T → R d × R d , of convolution type, is allowed to be singular, thus encompassing fractional processes, in particular the recent literature on rough volatility [START_REF] Alòs | On the short-time behavior of the implied volatility for jump-diffusion models with stochastic volatility[END_REF][START_REF] Bayer | Pricing under rough volatility[END_REF][START_REF] Euch | The microstructural foundations of leverage effect and rough volatility[END_REF][START_REF] Gatheral | Volatility is rough[END_REF]. Components of the system are in general correlated, the correlation matrix being implicitly encoded in the diffusion coefficient σ ε . General existence and uniqueness results for such stochastic Volterra equations are so far out of reach, and our conditions below are sufficient and general enough for most applications. In order to state them precisely, we first introduce several definitions and concepts: Definition 2.2. For any ε > 0, a solution to (2.1) is an R d -valued progressively measurable stochastic process X ε satisfying (2.1) almost surely and such that

P t 0 |K(t -s)b ε (s, X ε s )| + |K(t -s)σ ε (s, X ε s )| 2 ds < ∞, for all t ∈ T = 1.
We shall assume that the (singular) convolution kernel satisfies the following condition, which is essentially a multivariate version of the one given in [1, Condition (2.5)]: Assumption 2.3. The kernel K : T → R d×d is an upper triangular matrix satisfying the following conditions: K ∈ L 2 (T : R d×d ) and there exists γ ∈ (0, 2] such that, for h small enough,

h 0 |K(t)| 2 dt + T 0 |K(t + h) -K(t)| 2 dt = O(h γ ).
We refer to [1, Example 2.3] for a broad range of kernels that satisfy this assumption. Of particular interest in mathematical finance is the Riemann-Liouville kernel K(t) = t H-1 2 , for H ∈ (0, 1 2 ) implying γ = 2H. Moreover if K is locally Lipschitz and K satisfies Assumption 2.3 then so does the product K K; this includes the gamma and power-law kernels which are related to the class of Brownian semistationary processes [START_REF] Barndorff-Nielsen | Brownian semistationary processes and volatility/intermittency[END_REF].

Remark 2.4. This setup covers in particular the following two useful forms for the kernel:

• K = Diag(k 1 , • • • , k d ) is a diagonal matrix, where each k i : T → R satisfies Assumption 2.3.
• The drift and diffusion coefficients of any sub-system of (2.1) can be convoluted with different kernels. As an example, the one-dimensional SVE

Y ε t = Y ε 0 + t 0 K 1 (t -s)b Y ε (s, Y ε s ) ds + t 0 K 2 (t -s)σ Y ε (s, Y ε s ) dW (1)
s is the first component of (2.1) with d = 2 and

K = K 1 K 2 0 0 , b ε = b Y ε 0 , σ ε = 0 0 σ Y ε 0 .
One could in fact use general matrices to eliminate the auxiliary state, such as

K = (K 1 K 2 ) : T → R 1,2
in the example above. We stick to square matrices for consistency.

Volterra systems appearing in the literature, and in particular in the mathematical finance one, have a specific structure in the sense that only one component satisfies an SVE with (singular) kernel, and can be dealt with independently of the other component. This particular structure allows us to relax some conditions on the coefficients, and we shall leverage on it whenever needed. We make this more specific through the following two definitions: Definition 2.5. Let Υ ⊂ 1, d and Γ : R |Υ| → [0, ∞). We define S Γ Υ as the set of functions f for which there exists a strictly positive constant

C Υ such that, for all x ∈ R d , (2.2) |f (x)| ≤ C Υ 1 + |x| Υ c + Γ x (Υ) ,
where |x| Υ c := i∈Υ c x (i) and x (Υ) := (x (i) ) i∈Υ .

Definition 2.6. The process X ε admits an autonomous S Γ Υ -subsystem {X ε,(l) } l∈Υ if for all 1 ≤ i, j ≤ d such that K ij = 0, b (j) ε (t, •) and all the components of the row σ (j) ε (t, •) satisfy the following for small enough ε and uniformly in t ∈ T:

• if i ∈ Υ, they have linear growth and do not depend on X ε,(k) for k ∈ Υ c ;

• if i ∈ Υ c , they belong to S Γ Υ . Example 2.7. The motivation for Definition 2.6 is to be able to handle (rough) stochastic volatility models, ubiquitous in mathematical finance, where linear growth of all the coefficients may not hold. Consider for example the rough Bergomi model [START_REF] Bayer | Pricing under rough volatility[END_REF] 

      X (1) t = - 1 2 t 0 exp X (2) s ds + t 0 exp 1 2 X (2) s dB s , X (2)
t = y 0 -at 2H + t 0 (t -s) H-1 2 dW s ,
where B and W are Brownian motions with correlation ρ ∈ (-1, 1). After dropping the dependence in ε, this fits into the setup of (2.1)

with d = 3, y 0 ∈ R, a > 0, H ∈ (0, 1 2 ), ρ = 1 -ρ 2 and K(t) =   1 0 0 0 t 2H-1 t H-1 2 0 0 0   , b t, (x 1 , x 2 ) =   -e x2 /2 -a/(2H) 0   , σ t, (x 1 , x 2 ) =   ρ e x2/2 ρ e x2/2 0 0 1 0   ,
where the third component is meaningless but allows us to handle the two different kernels. Here X admits X (2) as autonomous subsystem with Υ = {2} and Γ(x 2 ) = 1 + e x2 .

The following set of assumptions, inspired by [START_REF] Chiarini | On large deviations for small noise Itô processes[END_REF], completes our framework:

H1. X ε 0 converges to x 0 ∈ R d as ε tends to zero.

H2. For all ε > 0 small enough, the coefficients b ε and σ ε are measurable maps on T × R b) The process X ε admits an autonomous S Γ Υ -subsystem. H4. The SVE (2.1) is pathwise unique for small enough ε > 0.

H2 ensures that, on compact subsets of T × R d , the convergence of b ε and σ ε is uniform and that b and σ are uniformly continuous. H1, H2, H3a are standard and easily verifiable. H3b is unusual but includes a large number of functions; Assumption 3.9 will complete it to indicate the role of Γ such as to include Example 2.7. Moreover, the growth conditions from H3 are uniform in ε and therefore apply to the limits b and σ. The main restrictions arise from H4, although the latter is satisfied if, for instance, the coefficients b ε and σ ε are locally Lipschitz continuous for small enough ε > 0. This condition was relaxed in [START_REF] Mytnik | Uniqueness for Volterra-type stochastic integral equations[END_REF] to the one-dimensional case where K(t) = t -α , for α ∈ (0, 1 2 ) and σ(x) = x γ , for γ ∈ ( 1 2(1-α) , 1], which is clearly not Lispchitz continuous. Furthermore, to the best of our knowledge there currently exists no pathwise LDP for stochastic equations where pathwise uniqueness fails.

One can compare our setup with the SVE considered by Zhang in [START_REF] Zhang | Euler schemes and large deviations for stochastic Volterra equations with singular kernels[END_REF], where a pathwise LDP was derived under the assumptions of Lipschitz continuity and linear growth of the coefficients. We relax both these assumptions. Indeed, our framework covers Hölder-continuous diffusion coefficients, as mentioned in the previous paragraph, and functions with non-linear growth through the concept of autonomous S Γ Υ -subsystem, as in Example 2.7.

LARGE AND MODERATE DEVIATIONS

As discussed in the introduction, our goal is to provide pathwise large and moderate deviations for the general convolution stochastic Volterra system (2.1), and then extend these to nonconvolution kernels. The classical Freidlin-Wentzell approach, used in [START_REF] Freidlin | Random Perturbations of Dynamical Systems[END_REF], has limitations regarding the behaviour of the coefficients, and the rate function is often rather cumbersome to write. We follow here instead the weak convergence approach developed by Dupuis and Ellis [START_REF] Dupuis | A Weak Convergence Approach to the Theory of Large Deviations[END_REF]. We first introduce the reader to their abstract setting, and refine the large deviations result by Budhiraja and Dupuis [START_REF] Budhiraja | A variational representation for positive functionals of infinite dimensional Brownian motion[END_REF] to our general setup. We then show how this abstract framework applies to the small-noise stochastic Volterra system (2.1), first proving pathwise large deviations, and then the moderate deviations counterpart.

3.1. Weak convergence approach: the abstract setting. Given a family of Borel-measurable functions {G ε } ε>0 from W m to W d , we enquire about the large deviations behaviour of the family of random variables {G ε (W )} ε>0 as ε tends to zero, where W is a standard Brownian motion on the filtered probability space above. For each N > 0, the spaces of bounded deterministic and stochastic controls

S N := v ∈ L 2 : T 0 |v s | 2 ds ≤ N and A N := v ∈ A : v ∈ S N almost surely , (3.1)
with A introduced in (1.4), are equipped with the weak topology on L 2 (T × Ω) such that they are closed and even compact (by Banach-Alaoglu-Bourbaki theorem). Budhiraja and Dupuis [START_REF] Budhiraja | A variational representation for positive functionals of infinite dimensional Brownian motion[END_REF] assume, for any sequence {v ε } ε>0 in A N converging weakly to v ∈ A N , the existence of a limit in distribution of

G ε W + 1 ε • 0 v ε
s ds which is uniquely characterised by v. However, such uniqueness may fail when the coefficients of the system (1.1) (in particular the diffusion coefficient σ) are not locally Lipschitz, as is the case for the Feller diffusion for example (in this case without singular kernel, a dedicated analysis was carried out in [START_REF] Conforti | On small-noise equations with degenerate limiting system arising from volatility models[END_REF][START_REF] Donati-Martin | Large deviations for squares of Bessel and Ornstein-Uhlenbeck processes[END_REF] using the Freidlin-Wentzell approach). We relax here this uniqueness assumption by replacing the limiting trajectory by a perturbed version.

φ = lim n↑∞ G εn W + ε -1 n • 0 v εn s ds in distribution , and G 0 v,N is empty if v is not in A N . For all v ∈ A, we also denote G 0 v := N ∈N G 0 v,N .
Then we define the functional I :

W d → R + given by (3.2) I(φ) := inf 1 2 T 0 |v s | 2 ds : v ∈ L 2 such that φ ∈ G 0 v .
Definition 3.2. We say that φ ∈ W d is uniquely characterised if there exists a sequence {v n } n∈N ⊂ L 2 such that

(3.3) G 0 v n = {φ} and 1 2 T 0 |v n s | 2 ds ≤ I(φ) + 1 n , for all n ∈ N.
In particular, if there exists v ∈ L 2 which attains the infimum in (3.2) and G 0 v = {φ} then φ is uniquely characterised, because one can choose v n = v for all n ∈ N.

In the display (3.2) and Definition 3.2, φ, v, v n and v are all deterministic. Assumption 3.3. For any δ > 0 and any φ ∈ W d such that I(φ) < +∞, there exists φ δ uniquely characterised such that φ -φ δ T ≤ δ and I(φ) -I(φ δ ) ≤ δ. Remark 3.4. This assumption is reminiscent of [START_REF] Donati-Martin | Large deviations for squares of Bessel and Ornstein-Uhlenbeck processes[END_REF]Proposition 3.3], where the authors resolve the non-uniqueness issue in the diffusion case. A similar problem is also at the core of [START_REF] Budhiraja | Many-server asymptotics for join-the-shortest-queue: Large deviations and rare events[END_REF]Lemma 5.1] in an infinite-dimensional setting.

Our abstract large deviations result is the following, extending [START_REF] Budhiraja | A variational representation for positive functionals of infinite dimensional Brownian motion[END_REF]Theorem 4.4], at least when the underlying Hilbert space is L 2 (T : R m ), to the non-uniqueness case. Theorem 3.5. Assume that (i) For all N > 0, all v ∈ A N and all families {v ε } ε>0 in A N converging in distribution to v as ε tends to zero, We defer the proof to Appendix A.1; the lower bound can be tackled as in [START_REF] Budhiraja | A variational representation for positive functionals of infinite dimensional Brownian motion[END_REF], and we therefore concentrate on the upper bound. The idea is that the Laplace principle (1.2) upper bound involves an infimum so deriving it only requires a δ-optimal path. Hence a perturbation will also do the trick, provided one knows how to handle the control associated to it. In [START_REF] Budhiraja | A variational representation for positive functionals of infinite dimensional Brownian motion[END_REF]Theorem 4.4], unique characterisation of the limiting element in (i) is granted, and the set G 0 v is a singleton that takes the form G 0 ( • 0 v s ds), where they view G 0 as a map. In that case Assumption 3.3 is clearly satisfied since φ δ can be taken as φ itself.

G ε W + ε -1 • 0 v ε s ds ε>0 is tight. (ii)
3.2. Application to stochastic Volterra systems. We now show how the abstract setting developed above in Section 3.1 applies to the small-noise stochastic Volterra system (2.1) and why pathwise uniqueness is so fundamental. If H4 holds, define the functional G ε as the Borel-measurable map associating the multidimensional Brownian motion W to the solution of the stochastic Volterra equation (2.1), that is: G ε (W ) = X ε . For any control v ∈ A N , N > 0 (introduced in (3.1)) and any ε > 0, the process W := W + ϑ -1 ε • 0 v s ds is a P-Brownian motion by Girsanov's theorem, where

d P dP := exp - 1 ϑ ε m i=1 T 0 v (i) s dW (i) s - 1 2ϑ 2 ε T 0 |v s | 2 ds .
Hence the shifted version X ε,v := G ε ( W ) appearing in Theorem 3.5(i) is the strong unique solution of (2.1) under P, with X ε and W replaced by X ε,v and W . Because P and P are equivalent, X ε,v is also the unique strong solution, under P, of the controlled equation

(3.4) X ε,v t = X ε 0 + t 0 K(t -s) b ε (s, X ε,v s ) + σ ε (s, X ε,v s )v s ds + ϑ ε t 0 K(t -s)σ ε (s, X ε,v s )dW s .
Under appropriate conditions, and using the notations set in H1, H2, we heuristically observe that taking ε to zero, the system (3.4) reduces to the deterministic Volterra equation

(3.5) φ t = x 0 + t 0 K(t -s) b(s, φ s ) + σ(s, φ s )v s ds.
We will show later that the set G 0 v corresponds to the set of solutions of (3.5). Example 3.8. To illustrate the need for a set G 0 v rather than a singleton, consider the Feller diffusion

X t = x 0 + κ t 0 (θ -X s ) ds + t 0 X s dW s ,
for t ∈ T, with x 0 , κ, θ > 0. Letting t → εt and denoting X ε t := X εt yields, by scaling,

X ε t = x 0 + κε t 0 (θ -X ε s ) ds + √ ε t 0 X ε s dW s , which is exactly (2.1) with d = 1, K ≡ 1 , ϑ ε = √ ε , b ε (x) = κε(θ -x), σ ε (x) = √ x.
For v ∈ L 2 , taking limits as ε tends to zero in the corresponding controlled equation (3.4) yields (3.5), or

φ t = x 0 + t 0 φ s v s ds, t ∈ T.

Uniqueness of this Volterra equation does not hold in general because of the non-Lipschitz coefficient, and thus G 0

v corresponds to the set of non-negative solutions. Consider for example x 0 = 1, T = [0, 4] and the control

v t := -1, if t ∈ [0, 2) 1, if t ∈ [2, 4].
The function

φ t := (t-2) 2 4
is clearly a solution, but so is ϕ equal to φ on [0, 2] and null on [START_REF] Alòs | On the short-time behavior of the implied volatility for jump-diffusion models with stochastic volatility[END_REF][START_REF] Barndorff-Nielsen | Brownian semistationary processes and volatility/intermittency[END_REF]. The square root function is indeed locally Lipschitz away from zero, and uniqueness can thus be guaranteed as long as the solution remains positive. The perturbation φ δ := φ + δt is now the unique solution to [START_REF] Donati-Martin | Large deviations for squares of Bessel and Ornstein-Uhlenbeck processes[END_REF]Proposition 3.3] shows that φ δ satisfies Assumption 3.3.

φ δ t = 1 + t 0 φ δ s v δ s ds, for all t ∈ [0, 4], where v δ := φδ / φ δ . The infimum in (3.2) is attained by v δ and G 0 v δ = {φ δ }, thus φ δ is uniquely characterised. Furthermore,
In [START_REF] Bhamidi | Rare event asymptotics for exploration processes for random graphs[END_REF], the authors were also confronted to a limiting equation with multiple solutions. Instead of perturbing the path φ, they perturb the control in a way that the resulting equation has a unique solution which is precisely φ, i.e. G 0 v δ = {φ}. This approach may seem more natural; however, it is not always obvious how to perturb the control ensuring uniqueness of the ODE, while our formulation makes it more straightforward. Before stating the main large and moderate deviations results for small-noise stochastic Volterra equations, we introduce the following assumption, monitoring the moments of the controlled equation: Assumption 3.9. Let X ε,v be the pathwise unique solution to (3.4). If H3a holds then the present assumption is satisfied. If instead H3b holds, then there exists ε 0 > 0 such that, for any p ≥ 1 and N > 0,

(3.6) sup E Γ (X ε,v t ) (Υ) p : t ∈ T, v ∈ A N , ε ∈ (0, ε 0 ) < ∞, (3.7) sup Γ (φ t ) (Υ) : t ∈ T, v ∈ S N , φ ∈ G 0 v < ∞.
Remark 3.10. In the following, H3b will always be complemented by Assumption 3.9.

Large Deviations.

Armed with the abstract setting in Section 3.1, and its application to the stochastic Volterra system (2.1) in Section 3.2, we can at last show large deviations for the latter:

Theorem 3.11 (Large Deviations). Under H1 -H4, Assumptions 2.3, 3.3 and 3.9, the family {X ε } ε>0 , unique solution of (2.1), satisfies a Large Deviations Principle with rate function (3.2) and speed ϑ -2 ε , where G 0 v is the set of solutions of the limiting equation (3.5). Remark 3.12. We recall that Assumption 3.3 is automatically satisfied if the limiting equation (3.5) has a unique solution. Also, it is only necessary to check Assumption 3.9 if H3a does not hold.

Technical preliminary results

. The proof will rely on the following results: Lemma 3.13 (proved in Section A.2) shows the moment bound of the controlled process defined by (3.4), Lemma 3.15 (proved in Section A.3) demonstrates the tightness and Lemma 3.18 (proved in Section A.4) deals with the compactness of the level sets of the rate function. They will then allow the use of Theorem 3.5.

Lemma 3.13 (LDP Moment bound)

. Under H1 -H4, Assumptions 2.3 and 3.9, for all p ≥ 2, N > 0, v ∈ A N and ε > 0 small enough, there exists a constant c > 0 independent of ε, v, t such that

(3.8) sup t∈T E |X ε,v t | p ≤ c.
Remark 3.14. This bound also holds for any solution φ of (3.5) under the same assumptions, therefore there also exists

c 0 > 0 such that sup{ φ T : v ∈ S N such that φ ∈ G 0 v } ≤ c 0 .
The following lemma deals with 3.11(i) by showing tightness of

X ε,v ε = G ε W + ϑ -1 ε • 0 v ε s ds . Lemma 3.15 (LDP Tightness). Consider H1 -H4, Assumptions 2.3 and 3.9. If p > 2 ∨ 2/γ, N > 0 and {v ε } ε>0 is a family in A N , then X ε,v ε
admits a version which is Hölder continuous on T of any order α < γ/2-1/p, uniformly for all ε > 0. Denoting again this version by X ε,v ε , one has for all ε > 0 small enough

(3.9) E     sup 0≤s<t≤T X ε,v ε t -X ε,v ε s |t -s| α   p  ≤ C, for all α ∈ [0, γ/2 -1/p), where C is a constant independent of ε, v ε , t. Moreover, the family of random variables {X ε,v ε } ε>0 is tight in W d .
Remark 3.16. This lemma entails that for all N > 0, v ∈ S N , any solution to (3.5) also has Hölder continuous paths of the same order.

Lemma 3.17. The set G 0 v from Definition 3.1 is characterised by

G 0 v = φ : Ω → W d φ t = x 0 + t 0 K(t -s) b(s, φ s ) + σ(s, φ s )v s ds, for all t ∈ T .
The following lemma proves Theorem 3.5(ii) and its proof can be found in Appendix A.4.

Lemma 3.18 (LDP Compactness

). Under H2, H3, Assumptions 2.3 and 3.9, the functional I in (3.2) has compact level sets.

Leveraging on the above lemmas, the Large Deviations Principle (Theorem 3.11) is a direct consequence of Theorem 3.5.

3.3.2.

Proof of Lemma 3.17. For N ∈ N and v ∈ A N , we first need to identify G 0 v,N , defined in (3.1). Consider a subsequence {ε n } n∈N ⊂ R + with lim n↑∞ ε n = 0 and a sequence {v εn } n∈N ⊂ A N such that lim n↑∞ v εn = v in distribution, and assume that X n := X εn,v εn converges in distribution to some random variable φ with values in W d . We also denote v n , ε n , X n 0 , b n , σ n along this subsequence.

By Skorohod representation theorem we can work with almost sure convergence for the purpose of identifying the limit. Hence X n , v n n≥0 converges almost surely in the product topology on W d × S N , and the limit is the W d × S N -valued random variable (φ, v). The convergence of the couple also takes place in distribution, so that we can follow the technique in [START_REF] Chiarini | On large deviations for small noise Itô processes[END_REF] to identify the limit. For t ∈ T, define Φ t :

S N × W d → R as Φ t (f, ω) := ω t -x 0 - t 0 K(t -s) b(s, ω s ) + σ(s, ω s )f s ds ∧ 1.
Clearly, Φ t is bounded and we show that it is also continuous. Indeed, let ω n → ω in W d and f n → f in S N with respect to the weak topology. H2 implies the existence of continuous moduli of continuity ρ b and ρ σ for both coefficients on compact subsets (see Definition 2.1). Since the paths ω n , n ≥ 1 and ω are continuous, they are also uniformly bounded and hence these moduli are available. Then, using Cauchy-Schwarz inequality and the fact that |x ∧ 1 -y ∧ 1| ≤ |x -y| for all x, y > 0,

|Φ t (f, ω) -Φ t (f n , ω n )| ≤ |ω t -ω n t | + t 0 |K(t -s)| |b(s, ω s ) -b(s, ω n s )| ds + t 0 |K(t -s)| σ(s, ω s ) -σ(s, ω n s ) f n s + σ(s, ω s ) f s -f n s ds ≤ ω -ω n T + ρ b (|ω -ω n |) T K 1 + ρ σ (|ω -ω n |) T K 2 f n 2 + σ(•, ω) T t 0 |K(t -s)| |f s -f n s | ds.
Since K(t -•) ∈ L 2 and f n tends to f weakly in L 2 then the last integral converges to zero as n goes to infinity. Moreover

lim n↑∞ ω -ω n T = 0, f n 2 ≤ √ N for all n ≥ 0 and K 2 + σ(•, ω) T < ∞, which proves that Φ t is continuous, and therefore lim n↑∞ E [Φ t (v n , X n )] = E [Φ t (v, φ)] .
We now prove that the left-hand side is actually equal to zero. We start with the observation that, using BDG inequality,

E [Φ t (v n , X n )] ≤ |X n 0 -x 0 | + t 0 |K(t -s)| E |b n (s, X n s ) -b(s, X n s )| ds + t 0 |K(t -s)| E [|σ n (s, X n s ) -σ(s, X n s )| |v n s |] ds + ϑ εn E t 0 |K(t -s)σ n (s, X n s )| 2 ds 1 2 . (3.10)
The bounds (A.2) and (A.4) show how to control the last term under H3a and H3b respectively, hence there exists C 1 > 0 independent of t and n such that

E t 0 |K(t -s)σ n (s, X n s )| 2 ds 1 2 ≤ C 1 .
However the convergence of b n , σ n only occurs on compact subsets so we use a localisation argument. For all n ≥ 0, M > 0 we introduce

A M n := ω ∈ Ω : X n (ω) T > M .
The uniform (in n ∈ N) Hölder regularity of X n , encompassed by (3.9), entails the existence, for all

p > 2 ∨ 2/γ, of C 2 (p), C 3 (p) > 0 independent of n such that E sup t∈T |X n t | p ≤ C 2 (p) |X n 0 | p + T pα ≤ C 3 (p),
for some 0 < α < γ/2 -1/p and where X n 0 is uniformly bounded by 2 |x 0 | for n large enough. Markov's inequality then implies that (3.11) lim

M↑∞ sup n∈N P A M n ≤ lim M↑∞ sup n∈N C 3 (p) M p = 0.
Moreover, for all n ∈ N, ω ∈ Ω \ A M n , and t ∈ T, |X n t (ω)| is bounded by M , which means

|b n (t, X n t (ω)) -b(t, X n t (ω))| ≤ b n (t, •) -b(t, •) M ,
which tends to zero uniformly on T as n goes to infinity (and likewise for σ n ) from H2. Define now

I n := t 0 |K(t -s)| |b n (s, X n s ) -b(s, X n s )| + |σ n (s, X n s ) -σ(s, X n s )| |v n s | ds
and observe that, using Jensen and Cauchy-Schwarz inequalities, the growth condition on the coefficients from H3 and the moment bounds on X n from (3.8), there exists C 4 > 0 independent of n such that

E |I n | 2 ≤ 2t t 0 |K(t -s)| 2 E |b n (s, X n s ) -b(s, X n s )| 2 ds + 2N t 0 |K(t -s)| 2 E |σ n (s, X n s ) -σ(s, X n s )| 2 ds ≤ 2 K 2 2 sup s≤t tE |b n (s, X n s ) -b(s, X n s )| 2 + N E |σ n (s, X n s ) -σ(s, X n s )| 2 ≤ C 4 . (3.12)
Let us fix ǫ > 0 and choose M ǫ > 0 large enough such that sup n∈N P A Mǫ n ≤ ǫ 2 /C 4 ; this choice is possible because of (3.11). Therefore, using the bound (3.12) and Cauchy-Schwarz inequality to separate I n and ½ A Mε n , one obtains lim sup n↑∞

E[I n ] = lim sup n↑∞ E I n ½ A Mε n + ½ Ω\A Mε n ≤ lim sup n↑∞ C 4 P A Mε n + K 1 b n -b Mε T + √ N K 2 σ n -σ Mε T ≤ ǫ.
It follows from (3.10) that

lim n↑∞ E Φ t (v n , X n ) ≤ lim n↑∞ |X n 0 -x 0 | + E[I n ] + ϑ εn C 1 ≤ ǫ, hence lim n↑∞ E [Φ t (v n , X n )] = 0 since ǫ > 0 was chosen arbitrarily. The equality E [Φ t (v, φ)] = 0
implies that φ satisfies (3.5) almost surely, for all t ∈ T. Since φ has continuous paths, it satisfies (3.5) for all t ∈ T, almost surely, which means G 0 v,N consists of all the solutions of (3.5). Since this definition is independent of N , it extends to G 0 v , which yields the claim.

Moderate Deviations.

Let h ε tend to infinity such that ϑ ε h ε tends to zero as ε goes to zero and define X to be the limit in law of X ε , which we identified in the previous subsection as a solution of the Volterra equation (3.13)

X t = x 0 + t 0 K(t -s)b(s, X s )ds.
Then the MDP for {X ε } ε>0 is equivalent to the LDP for the family {η ε } ε>0 defined as

η ε := X ε -X ϑ ε h ε = G ε (W ) -X ϑ ε h ε =: T ε (W ),
where T ε : W m → W d is a Borel-measurable map for each ε > 0. Therefore η ε satisfies the following SVE for all ε > 0, and is its unique solution if H4 holds.

η ε t = X ε 0 -x 0 ϑ ε h ε + t 0 K(t -s) b ε s, X s + ϑ ε h ε η ε s -b(s, X s ) ϑ ε h ε ds + t 0 K(t -s) σ ε s, X s + ϑ ε h ε η ε s h ε dW s .
Similarly to the LDP case we are interested in a certain shift of the driving Brownian motion, controlled by v ∈ A. For all ε > 0, let

(3.14) η ε,v := T ε W + h ε • 0 v s ds = G ε W + h ε • 0 v s ds -X ϑ ε h ε .
For convenience we introduce the sequence {Θ ε,v } ε>0 defined, for all ε > 0, v ∈ A, t ∈ T, by

Θ ε,v t := G ε W + h ε • 0 v s ds (t) = X ε 0 + t 0 K(t -s) b ε (s, Θ ε,v s ) + ϑ ε h ε σ ε (s, Θ ε,v s )v s ds + ϑ ε t 0 K(t -s)σ ε (s, Θ ε,v s )dW s .
This sequence satisfies the bound (3.8) and converges weakly towards X since ϑ ε h ε σ ε tends to zero as ε tends to zero. Finally the process defined by (3.14) satisfies

η ε,v t = X ε 0 -x 0 ϑ ε h ε + t 0 K(t -s) b ε s, Θ ε,v s -b(s, X s ) ϑ ε h ε ds + t 0 K(t -s)σ ε s, Θ ε,v s v s ds + 1 h ε t 0 K(t -s)σ ε s, Θ ε,v s dW s . (3.15)
For all v ∈ A, we define T 0 v to be the solution of the limiting equation

(3.16) ψ t = t 0 K(t -s) ∇b(s, X s )ψ s + σ(s, X s )v s ds.
The form of the limit equation is dramatically simpler than for the LDP and much easier to compute. Moreover T 0 v is well defined because the linearity of the equation and Assumption 2.3 grant uniqueness for free, provided ∇b exists. Hence we will need the following assumptions: H5. For each t ∈ T, the function b(t, •) is continuously differentiable and b is Lipschitz continuous. H6. There exists δ > 0 such that σ(t, •) is locally δ-Hölder continuous, uniformly for all t ∈ T. H7.

lim ε↓0 ϑ ε h ε -1 |X ε 0 -x 0 | = 0. H8. There exist ε 0 > 0, a sequence {ν ε } ε>0 with lim ε↓0 ν ε (ϑ ε h ε ) -1 = 0 and a function Ξ : R d → R such that |b ε (t, x) -b(t, x)| ≤ ν ε Ξ(x) for all t ∈ T, ε ∈ (0, ε 0 )
, where for all p ≥ 1, N > 0,

(3.17) sup E Ξ Θ ε,v t p , ε ∈ (0, ε 0 ), v ∈ A N , t ∈ T < ∞.
(3.18) Λ(ψ) := inf 1 2 T 0 |v t | 2 dt : v ∈ L 2 , ψ = T 0 v , and
Λ(ψ) = +∞ if this set is empty.
The proof of the moderate deviations theorem follows a similar structure to that of Theorem 3.11, making use of Theorem 3.5. It will rely on moment bounds in Lemma 3.21 (proved in Section B.1), tightness in Lemma 3.22 (proved in Section B.2), weak convergence in Lemma 3.23 (proved in Section B.3), and finally compactness of the level sets in Lemma 3.24.

Lemma 3.21 (MDP Moment bound)

. Under H2 -H5, H7, H8, Assumptions 2.3 and 3.9, for all p ≥ 2, N > 0, v ∈ A N and ε > 0 small enough, there exists c > 0 independent of ε, v, t such that (3.19) sup

t∈T E |η ε,v t | p ≤ c.
Lemma 3.22 (MDP tightness). Let p > 2 ∨ 2/γ, N > 0 and a family {v ε } ε>0 in A N . Under H2 -H5, H7, H8, Assumptions 2.3 and 3.9, η ε,v ε admits a version which is Hölder continuous on T of any order α < γ/2 -1/p, uniformly for all ε > 0. Denoting again this version by η ε,v ε , one has for all ε > 0 small enough,

(3.20) E sup 0≤s<t≤T η ε,v ε t -η ε,v ε s |t -s| α p ≤ C, for all α ∈ [0, γ/2 -1/p), where C is a constant independent of ε, v ε , s, t. Moreover, the fam- ily {η ε,v ε } ε>0 is tight in W d .
We recall that η ε,v ε = T ε W + h ε • 0 v s ds and ψ = T 0 v , hence the lemma above deals with Theorem 3.5 (i). The following one identifies the limit set as the unique solution to (3.16). It is thus more precise than in the LDP case, and justifies the form of the rate function (3.18).

Lemma 3.23 (MDP weak convergence).

Let N > 0, a family {v ε } ε>0 such that, for all ε > 0, v ε ∈ A N and v ε converges in distribution to v ∈ A N , and ψ the unique solution of (3.16). Under H2 -H8, Assumptions 2.3 and 3.9, η ε,v ε converges in distribution to ψ as ε goes to zero.

Item (ii) is dealt with in the following lemma. Lemma 3.24 (MDP compactness). Under H2, H3, H5, Assumptions 2.3 and 3.9, the functional Λ defined by (3.18) has compact level sets.

Proof. Noticing that ∇b(t, X t ) and σ(t, X t ) are uniformly bounded on T by continuity, this lemma boils down to a particular case of Lemma 3.18. Theorem 3.5(iii) is immediate by uniqueness of (3.16), therefore all the conditions are met and Theorem 3.20 follows as a direct application of Theorem 3.5.

3.5. Extension to non-convolution kernels. The analysis undertaken in this paper is based, both for notational convenience and with a view towards application, on convolution kernels. Different assumptions were studied in the literature, in particular Decreusefond [START_REF] Decreusefond | Regularity properties of some stochastic Volterra integrals with singular kernels[END_REF] considered the properties of the map f → • 0 K(•, s)f (s) ds in order to include the fractional Brownian motion in his setting.

3.5.1. Setting. We call a kernel a map K : T 2 → R for which both t 0 K(t, s) 2 ds and K(t, s) are finite for all t ∈ T and s = t. The associated space is defined as

K := u : T → R, {F t }-progressively measurable, such that E t 0 K(t, s)u(s) 2 ds < ∞, for all t ∈ T .
Hence, for all u ∈ K the stochastic integral

M K t (u) := t 0 K(t, s)u(s)dW s
is well defined for all t ∈ T in the Itô sense. For any α ∈ (0, 1), we denote the Riemann-Liouville integral I α and derivative D α as

(3.21) (I α f )(t) := 1 Γ(α) t 0 (t -s) α-1 f (s)ds, (D α f )(t) := d dt (I 1-α f )(t), for f ∈ L 1 , t ∈ T. Define I α,p := I α (L p ) equipped with the norm f Iα,p := D α f L p . If α > 1 p , then I α,p ⊂ C α-1 p 0
, the space of (α -1 p )-Hölder continuous functions null at time 0. Let K denote the linear map associated to K(t, s) by

(3.22) Kf (t) := t 0 K(t, s)f (s)ds,
and introduce, for x ∈ (0, ∞) \ {2},

(3.23) θ(x) := 2x 2 -x .
Given the space inclusions above, the following assumption implies precise Hölder regularity for the integral (3.22): Assumption 3.25. There exist χ ∈ (1, 2) and γ > 1/θ(χ) for which K is continuous from L 2 to I γ+ 1 2 ,2 and from L χ to I γ,θ(χ) . Example 3.26. The operators associated to the following kernels satisfy Assumption 3.25:

• The Riemann-Liouville kernel J H (t, s) = (t -s) H-1 2 + Γ(H + 1 2 )
, with H ∈ (0, 1), satisfies this assumption with γ = H and any χ < 2 [32, Theorem 4.1].

• The fractional Brownian motion kernel

K H (t, s) = (t -s) H-1 2 + Γ(H + 1 2 ) F H - 1 2 , 1 2 -H, H + 1 2 , 1 - t s ,
where F is the Gauss hypergeometric function, also satisfies this assumption with the same parameters as above [

∈ K ∩ L θ(χ) (Ω × T). Then M K (u) has a measurable version M K (u) which is α-Hölder continuous for all α < γ -1/θ(χ).
From now on, we only consider the measurable version of the stochastic integral. Although this theorem was proved in a one-dimensional setting, it also covers multi-dimensional stochastic Volterra integrals by considering their components individually and summing them.

X ε t = X ε 0 + t 0 K(t, s)b ε (s, X ε s ) ds + ϑ ε t 0 K(t, s)σ ε (s, X ε s )dW s , t ∈ T,
which was studied in [START_REF] Coutin | Volterra differential equations with singular kernels[END_REF] without the ε-dependence, and where the coefficients live in the same spaces as those from (2.1). To complete the non-convolution setup we also need the following condition.

Assumption 3.28. There exists q > 2 such that (3.25)

sup t∈T t 0 |K(t, s)| -θ(q) ds < ∞,
with θ(q) introduced in (3.23).

Let p > q > 2, then 2 < -θ(p) < -θ(q) and Hölder's and Jensen's inequalities yield

t 0 |K(t, s)f (s)| 2 ds p 2 ≤ t 0 |K(t, s)| -θ(p) ds p-2 2 t 0 |f (s)| p ds ≤ t p-q q t 0 |K(t, s)| -θ(q) ds p -θ(q) t 0 |f (s)| p ds.
This replaces the Gronwall-type inequality derived for convolution kernels in Lemma A.1. Hence replacing Assumption 2.3 by the condition (3.25) one recovers the moments bounds of Lemmata 3.13 and 3.21 for the processes {X ε,v } ε>0 and {η ε,v } ε>0 and for any p ≥ 1. Setting in particular p = θ(χ) from Theorem 3.27, then for any v ∈ A N , N > 0 and ξ ∈ {b ε , σ ε } we have ξ(X ε,v ) ∈ K ∩ L p (Ω × T) thanks to the growth conditions H3. Therefore Assumption 3.25, Theorem 3.27 and Assumption 3.28 yield the almost sure Hölder regularity of the following processes defined on T:

• 0 K(•, s)ξ(X ε,v s ) ds and • 0 K(•, s)ξ(X ε,v s ) dW s .
Hence we recover the tightness of Lemmata 3.15 and 3.22 under this new set of assumptions. Notice that we can consider kernels consisting of both convolution and non-convolution components.

Finally we can extend the LDP and MDP results without further modifications:

APPLICATION TO ROUGH VOLATILITY

We now show how our results (Theorems 3.11, 3.20 and 3.29) apply to a large class of models recently developed in mathematical finance. Originally proposed by Comte and Renault [START_REF] Comte | Long memory in continuous-time stochastic volatility models[END_REF] with financial econometrics applications in mind, rough volatility models were rediscovered later in the context of option pricing in [START_REF] Alòs | On the short-time behavior of the implied volatility for jump-diffusion models with stochastic volatility[END_REF][START_REF] Bayer | Pricing under rough volatility[END_REF][START_REF] Fukasawa | Asymptotic analysis for stochastic volatility: martingale expansion[END_REF][START_REF] Gatheral | Volatility is rough[END_REF], developed and extended widely, and have now become the new standards of volatility modelling. They usually take the following form:

(4.1)        X t = - 1 2 t 0 Σ(Y s ) ds + t 0 Σ(Y s )dB s , Y t = y 0 + t 0 K 1 (t -s)b(Y s ) ds + t 0 K 2 (t -s)ζ(Y s )dW s ,
where both X and Y are one-dimensional, K 1 , K 2 ∈ L 2 (T : R + ) and B and W are two standard Brownian motions with d B, W t = ρ dt, for some correlation parameter ρ ∈ (-1, 1). We further define ρ := 1 -ρ 2 , and set X 0 = 0 without loss of generality. Here X denotes the logarithm of a stock price process, and Σ(Y ) its instantaneous volatility. We adopt a slight abuse of notation, as X previously denoted the multidimensional system, but writing now X as the log-stock price is consistent with the mathematical finance literature and should not create any confusion. We summarise in Table 1 the most common rough volatility models used in mathematical finance, indicating where their asymptotic behaviours were covered, and where our framework not only encompasses those, but fills the gaps so far missing. As discussed below, our application to the rough Heston model is conditional on the latter to have a unique pathwise solution, a problem that remains open so far. The detailed analysis of these cases is then provided in Section 4.2 in the small-time case, and in Section 4.3 for their tail behaviours.

Models

Small-time rescaling (general).

In the small-time case, we need to assume some scaling behaviour for the kernel functions. We say that a function f :

R → R is homogeneous of degree α ∈ R if f (λx) = λ α f (x) holds for all x, λ ∈ R. Assumption 4.1. K 1 and K 2 are homogeneous of degrees ̟ ∈ (-1 2 , 1 2 ] and H -1 2 ∈ (-1 2 , 1 2 ].
Since K 1 is homogeneous of degree ̟, then

h 0 K 1 (t) 2 dt = h 0 K 1 (1) 2 t 2̟ dt ≤ K 1 (1) 2 1 + 2̟ h 1+2̟ , for any h > 0, T 0 K 1 (t + h) -K 1 (t) 2 dt = K 1 (1) 2 T 0 (t + h) ̟ -t ̟ 2 dt = O h 1+2̟ , for h small enough,
and so Assumption 2.3 is satisfied with γ = 1 + 2̟ ∈ (0, 2], and likewise for K 2 with γ = 2H. Under this assumption, the rescalings

X ε t := ε H-1 2 X εt and Y ε t := Y εt turn (4.1) into        X ε t = - ε H+ 1 2 2 t 0 Σ(Y ε s ) ds + ε H t 0 Σ(Y ε s )dB s , Y ε t = y 0 + ε 1+̟ t 0 K 1 (t -s)b(Y ε s ) ds + ε H t 0 K 2 (t -s)ζ(Y ε s )dW s , (4.2)
so that we are precisely in the framework of (2.1)

with d = 3, ϑ ε = ε H , K(t) =   1 0 0 0 K 1 (t) K 2 (t) 0 0 0   , b ε (t, (x, y)) =   -1 2 ε H+ 1 2 Σ(y) ε 1+̟ b(y) 0   , σ ε (t, (x, y)) =   ρ Σ(y) ρ Σ(y) 0 0 ζ(y) 0   ,
where, similarly to Example 2.7, the additional dimension allows to handle the two different kernels.

Note that σ ε does not depend on ε but encodes the correlation. The controlled equation (3.4) for the second component reads

Y ε,v t = y 0 + ε 1+̟ t 0 K 1 (t -s)b(Y ε,v s ) ds + ε H t 0 K 2 (t -s)ζ(Y ε,v s )dW s + t 0 K 2 (t -s)ζ(Y ε,v s )v s ds,
for each t ∈ T, ε > 0 and v ∈ A. Note that the dynamics of X ε do not feed back into Y ε and that Σ ∈ S |Σ| {2} in the sense of Definition 2.5. The following assumption stands throughout this section: Assumption 4.2 (Small-time assumptions).

• K 2 (t) := t H-1 2 /Γ(H + 1 2 ); • Σ,
E Σ Y ε,v t p < ∞;
• the equation for Y ε in (4.2) is pathwise unique for small enough ε > 0.

The choice of kernel K 2 is a common setup in rough volatility models and allow for more explicit results. These conditions ensure that H2 holds with limit coefficients b ≡ (0, 0, 0) ⊤ and σ = σ ε . Furthermore, Y ε is an autonomous subsystem in the sense of Definition 2.6 and H3 and the bound (3.6) hold. An pathwise unique solution of the system (4.1) exists since X ε is explicit from Y ε , and H4 is satisfied.

4.1.1. Large deviations. For each control v ∈ S N with N > 0, the limit equation (3.5) of the volatility in the large deviations regime reads (4.4)

ϕ t = y 0 + t 0 (t -s) H-1 2 Γ(H + 1 2 ) ζ(ϕ s )v s ds, for t ∈ T.
From the uniform bound on ϕ derived in Remark 3.14 and the continuity of Σ, we obtain that |Σ(ϕ t )| is uniformly bounded in t ∈ T and in v ∈ S N , hence (3.7) holds. Therefore Assumption 3.9 and H1 -H4 follow from Assumption 4.2. Mimicking the fractional integral notation from Section 3.5, we introduce for convenience the notations

I H+ 1 2 x (f ) := x + t 0 (t -s) H-1 2 Γ(H + 1 2 )
f s ds and I

H+ 1 2 x (L 1 ) := I H+ 1 2 x (f ), f ∈ L 1 ,
and the fractional derivative D is defined in (3.21). From now on, to simplify the statements, we write Z ε ∼ LDP(I, ε -1 ) to express that the family of random variables {Z ε } ε>0 satisfies an LDP with rate function I and speed ε -1 , as ε tends to zero.

(L1) (X ε , Y ε ) ∼ LDP I, ε -2H
, where I : W 2 → R + is given by

I(φ, ϕ) = inf 1 2 T 0 u 2 t + v 2 t dt : u, v ∈ L 2 , φ t = t 0 Σ(ϕ s ) ρu s + ρv s ds, ϕ = I H+ 1 2 y0 (ζ(ϕ)v) ; (L2) X ε ∼ LDP I X , ε -2H
, where

I X (φ) = inf I(φ, ϕ) : ϕ ∈ I H+ 1 2 y0 (L 1 ) , if φ ∈ AC 0 and infinity otherwise; (L3) ε H-1 2 X ε ∼ LDP I X 1 , ε -2H where I X 1 (x) = inf I X (φ) : φ 1 = x for all x ∈ R; (L4) Y ε ∼ LDP I Y , ε -2H
, where

I Y (ϕ) = 1 2 T 0 D H+ 1 2 (ϕ -y 0 )(t) ζ(ϕ t ) 2 ½ ζ(ϕt) =0 dt, if ϕ ∈ I H+ 1 2 y0
(L 1 ) and infinity otherwise;

(L5) Y ε ∼ LDP I Y 1 , ε -2H , with I Y 1 (y) = inf I Y (ϕ) : ϕ 1 = y for y ∈ R. While (
L1), (L2) and (L4) deal with pathwise large deviations, (L3) and (L5) are one-dimensional large deviations statements, about the marginal distributions of X and Y . In this small-time behaviour case, we recover the same scaling as in [START_REF] Forde | Small-time, large-time, and H → 0 asymptotics for the rough Heston model[END_REF][START_REF] Forde | Asymptotics for rough stochastic volatility models[END_REF].

Proof. (L1)

As discussed above, the assumptions of Theorem 3.11 are satisfied, so that the three-dimensional process (X ε , Y ε , Z ε ), where Z ε ≡ 0 for all ε > 0, satisfies an LDP with rate function

J(φ, ϕ, ψ) = inf 1 2 T 0 u 2 t + v 2 t + w 2 t dt : u, v, w ∈ L 2 , φ t = t 0 Σ(ϕ s ) ρu s + ρv s ds, ϕ = I H+ 1 2 y0 ζ(ϕ)v , ψ ≡ 0 . The map (X ε , Y ε , Z ε ) → (X ε , Y ε
) is continuous so that the contraction principle [33, Theorem 4.2.1] yields an LDP for (X ε , Y ε ) with rate function inf J(φ, ϕ, ψ) : ψ ≡ 0 = J(φ, ϕ, 0), which corresponds to I. (L2) Since the map (X ε , Y ε ) → X ε is continuous, the claim follows from the contraction principle. (L3) Projecting the pathwise large deviations (L2) onto the last coordinate point t = 1 is equivalent to applying the contraction principle, and the claim follows immediately. (L4) A direct application of Theorem 3.11 yields an LDP with rate function

I Y (ϕ) = inf 1 2 T 0 v 2 t dt : v ∈ L 2 , ϕ = I H+ 1 2 y0 (ζ(ϕ)v) ;
Inverting it as above ends the proof and (L5) follows from the contraction principle.

We observe that in some special cases one can reach a more explicit expression for I.

I(φ, ϕ) = 1 2ρ 2 T 0   φ2 t Σ(ϕ t ) ½ Σ(ϕt) =0 - 2ρ φt D H+ 1 2 (ϕ -y 0 )(t) ζ(ϕ t ) Σ(ϕ t ) ½ ζ(ϕt)Σ(ϕt) =0 + D H+ 1 2 (ϕ -y 0 )(t) ζ(ϕ t ) 2 ½ ζ(ϕt) =0   dt, if φ ∈ AC 0 and ϕ ∈ I H+ 1 2 y0
(L 1 ), and infinity otherwise.

Proof. We start from the definition of I given in Proposition 4.3(L1). For each φ ∈ AC 0 and ϕ ∈ I

H+ 1 2 y0
(L 1 ), reverting the integral which defines φ gives

u t = 1 ρ φt Σ(ϕ t ) -ρv t ½ Σ(ϕt) =0
, for all t ∈ T, because whenever Σ(ϕ t ) = 0, although u is not uniquely determined by φ, the optimal choice of control (the one minimising the cost) is u t = 0, see [27, Remark 2.3] for more details. In the uncorrelated case ρ = 0, the same reasoning for the equation that ϕ solves yields for all t ∈ T

v t = 1 ζ(ϕ t ) D H+ 1 2 (ϕ -y 0 )(t)½ ζ(ϕt) =0 .
Furthermore, in the special case where ζ(ϕ t ) = 0 almost everywhere the equality above holds almost everywhere, which is sufficient for the optimisation problem (because of the correlation, v may not be equal to zero even if ζ(ϕ) is). Plugging these into the rate function yields the claim. The last

condition stands because if φ / ∈ AC 0 or ϕ / ∈ I H+ 1 2 y0
(L 1 ) then they cannot satisfy the equations and therefore the infimum takes place over an empty set. 4.1.2. Moderate deviations. We now show how our moderate deviations results apply to the rough volatility model (4.1). Let h ε = ε -β for any β ∈ (0, H), and define the two-dimensional process

η ε := 1 ϑ ε h ε X ε , Y ε -y 0 = 1 ε H-β X ε , Y ε -y 0 .
The case β = 0 corresponds to the Central Limit Theorem, whereas β = H is the LDP regime, so that MDP precisely corresponds to some interpolation between the two. Regarding the assumptions note that y ε 0 = y 0 and b ε clearly tends to zero as ε goes to zero, hence it is trivial that b ≡ 0 is continuously differentiable and Lipschitz continuous, and thus H5 and H7 hold. Now let Assumption 4.2 hold. Denoting the i-th component of b with b (i) , one notices that b

ε (t, (x, y)) -b (2) (t, (x, y)) ≤ ε 1+̟ C L (1 + |y|) by linear growth of b. For H8 to hold, one then requires that ε 1+̟-(H-β) and ε H+ 1 2 -(H-β) both tend to zero as ε goes to zero. Moreover the bound (4.3) implies (3.17).

Assumption 4.5 (Moderate deviations assumptions).

• The parameters H, ̟ and β are such that (1 + ̟) ∧ (H + 1 2 ) > H -β; • There exists δ > 0 such that Σ and ζ are locally δ-Hölder continuous.

Notice that the first inequality is always satisfied if H ≤ 1 2 . Therefore, Assumptions 4.1, 4.2, 4.5 imply H1 -H8 and Assumptions 2.3 and 3.9. Similarly to the LDP case, and recalling the definition of MDP from the introduction, we write Z ε ∼ MDP β (Λ, l ε ) if in fact ε β-H (Z ε -Z) ∼ LDP(Λ, l ε ), for any l ε > 0 converging to zero as ε tends to zero, where Z is the limit in distribution of Z ε . We also denote the subset of W d of absolutely continuous functions by AC, and AC 0 := {φ ∈ AC, φ 0 = 0}, and refer to (3.21) for the definition of the Riemann-Liouville fractional derivative.

(M1) (X ε , Y ε ) ∼ MDP β Λ, ε -2β , where Λ : W 2 → R + is given by Λ(φ, ϕ) = 1 2ρ 2 T 0   φ2 t Σ(y 0 ) - 2ρ φt D H+ 1 2 (ϕ -y 0 )(t) ζ(y 0 ) Σ(y 0 ) + D H+ 1 2 (ϕ -y 0 )(t) ζ(y 0 ) 2   dt, if φ ∈ AC 0 and ϕ ∈ I H+ 1 2 y0
(L 1 ), and infinity otherwise;

(M2) X ε ∼ MDP β Λ X , ε -2β , where Λ X (φ) = 1 2Σ(y 0 ) T 0 φ2
t dt, if φ ∈ AC 0 and infinity otherwise;

(M3) ε H-1 2 X ε ∼ MDP β Λ X 1 , ε -2β , with Λ X 1 (x) = x 2 2Σ(y 0 ) , for x ∈ R; (M4) Y ε ∼ MDP β Λ Y , ε -2β , where Λ Y : W → R + is given by Λ Y (ϕ) = 1 2 T 0 D H+ 1 2 (ϕ -y 0 )(t) ζ(y 0 ) 2 ½ ζ(y0) =0 dt, if ϕ ∈ I H+ 1 2 y0
(L 1 ) and infinity otherwise;

(M5) Y ε ∼ MDP β (Λ Y 1 , ε -2β ), where Λ Y 1 (y) = 1 2 y 2 for y ∈ R.

Proof. (M1)

As discussed above, the assumptions of Theorem 3.20 are satisfied, thus it yields an MDP with rate function

Λ(φ, ϕ) = 1 2 T 0 u 2 t + v 2 t dt : u, v ∈ L 2 , φ t = t 0 Σ(y 0 ) ρu s + ρv s ds, ϕ = I H+ 1 2 y0 (ζ(y 0 )v) ,
and inverting it as in the LDP case gives the claim. (M2) The contraction principle implies that an MDP for X ε holds with rate function

Λ X (φ) = inf Λ(φ, ϕ) : ϕ ∈ I H+ 1 2 y0 (L 1 ) . Let ψ ∈ AC 0 such that ψ := D H+ 1 2 (ϕ-y0) ζ(y0)
, then the rate function translates to

Λ X (φ) = inf 1 2ρ 2
Λ X 1 (x) = inf Λ X (φ) : φ ∈ AC 0 , φ 1 = x = inf 1 2Σ(y 0 ) T 0 φ2 t dt : φ ∈ AC 0 , φ 1 = x .
Setting T = 1 the optimal path under the constraint φ 1 = x is φ t = xt by the Euler-Lagrange equation. Again, plugging it into the rate function ends the proof. (M4) Theorem 3.20 gives an MDP for Y ε with rate function

Λ Y (ϕ) = inf 1 2 T 0 v 2 t dt : v ∈ L 2 , ϕ = I H+ 1 2 y0 (ζ(y 0 )v) .

Inverting it yields (M4). (M5) By contraction principle one obtains Λ

Y 1 (y) = inf Λ Y (ϕ), ϕ 1 = y .
Then setting ψ as in (M2) it boils down to the same optimisation problem as for (M3).

As in the large deviations results, (M1), (M2) and (M4) correspond to pathwise statements, whereas (M3) and (M5) are finite-dimensional results about the marginal distributions. For the log-stock price, (M3) corresponds precisely to the moderately-out-of-the-money regime presented and justified in [START_REF] Friz | Option pricing in the moderate deviations regime[END_REF] (for diffusion volatility models), based on the observation that the range of observable strikes grows with maturity. Furthermore, one can always apply Theorem 3.20 in the degenerate case Σ(y 0 )ζ(y 0 ) = 0 although the rate functions take slightly different forms. 4.1.3. Implied volatility asymptotics. We can easily deduce from the above results the asymptotic behaviour of the implied volatility, a standard norm for quoting option prices. For each maturity t ≥ 0 and log-moneyness k ∈ R, the implied volatility σ(t, k) is the unique non-negative solution to C BS t, k, σ(t, k) = C(t, k), where C BS corresponds to the price of a European Call option under the Black-Scholes model, and C a given Call option price (for example in a rough volatility model). This notion is only well defined if the underlying stock price is a true martingale, which we have not assumed so far, and may require additional conditions on the coefficients. This will be the case though in all our examples below, but for now, with the current level of generality, we assume it: Assumption 4.7. The process exp(X ε ) in (4.2) is a true martingale for all small enough ε > 0.

Small-time implied volatility asymptotics can be derived from Properties 4.

lim t↓0 σ t, kt 1 2 -H 2 =        k 2 2 inf x≥k I X 1 (x) , if k > 0, k 2 2 inf x≤k I X 1 (x) , if k < 0.
(MDP) Under the same assumptions as Proposition 4.6 and for any β ∈ (0, H), k = 0, lim t↓0 σ t, kt

lim t↓0 t 2H log P(t H-1 2 X t ≥ k) = -inf x≥k I X 1 (x).
Meanwhile in the Black-Scholes model with constant volatility σ > 0 the log-price process satisfies X t = -σ 2 t 2 + σB t for all t ∈ T, and simple Gaussian computations yield the large deviations behaviour

lim t↓0 t 2H log P(X t ≥ kt 1 2 -H ) = - k 2 2σ 2 .
The claim then follows directly from [49, Corollary 7.1], and by symmetry for the case k < 0.

(MDP) Following the same arguments as above, we obtain

lim t↓0 σ t, kt 1 2 -β 2 =        k 2 2 inf x≥k Λ X 1 (x) , if k > 0, k 2 2 inf x≤k Λ X 1 (x) , if k < 0.
Plugging in the expression of Λ X 1 from (M3) finishes the proof.

This concludes the presentation of the general results for rough volatility models. The next sections display the diversity of the models found in the literature and how large and moderate deviations principles apply to them. [START_REF] Horvath | Asymptotic behaviour of randomised fractional volatility models[END_REF] is an extension of the classical Stein-Stein volatility model [START_REF] Stein | Stock price distributions with stochastic volatility -an analytic approach[END_REF] to the fractional setting. It corresponds to (4.1) with K 1 ≡ 1 (hence

Small-time rescaling (examples).

Rough Stein-Stein. The rough Stein-Stein, suggested in

̟ = 0), K 2 (t) = t H-1 2 /Γ(H + 1 2 ), H ∈ (0, 1 2 ), y 0 > 0, Σ(y) = y 2 , b(y) = κ(θ -y), κ, θ > 0 and ζ(y) ≡ ξ > 0.
The coefficients are Lipschitz continuous and well-behaved, hence Assumptions 4.2 and 4.5 are easily checkable and the limit equation ( 4.4) has a unique solution, hence Propositions 4.3 and 4.6 apply. Note that because ζ is a positive constant, Corollary 4.4 gives the rate function I in integral form and one can solve (L5) in closed-form using the Euler-Lagrange equation in a similar way as in the proof of Proposition 4.6. Furthermore, since Y is Gaussian its exponential moments are finite and Novikov's condition [60, Section 3.5.D] ensures that Assumption 4.7 holds. Therefore, Corollary 4.8 yields the small-time behaviour of the implied volatility. Notice that the LDP and MDP for this model still hold when replacing the Riemann-Liouville kernel with the standard fractional Brownian motion by virtue of Theorem 3.29. The pathwise LDP for this model was first derived in [START_REF] Horvath | Asymptotic behaviour of randomised fractional volatility models[END_REF] albeit with the different scaling

X ε t := ε H-1 2 +2β X εt and Y ε t := ε β Y εt , for β > 0.
4.2.2. Rough Bergomi. The rough Bergomi model as presented in [START_REF] Bayer | Pricing under rough volatility[END_REF] reads

         X t = - 1 2 t 0 V s ds + t 0 V s dB s , V t = V 0 exp t 0 (t -s) H-1 2 Γ(H + 1 2 ) dW s -at 2H ,
with V 0 > 0 and a ∈ R. A pathwise LDP for this model first appeared in [START_REF] Jacquier | Pathwise large deviations for the rough Bergomi model[END_REF] using the Freidlin-Wentzell approach and a tailored proof. This case is quite intricate because the exponential does not satisfy the linear growth bound but we circumvented this issue by introducing the notion of autonomous system, illustrated in Example 2.7 and completed by Assumption 3.9 and H3b. Not only does this framework unifies the result of [START_REF] Jacquier | Pathwise large deviations for the rough Bergomi model[END_REF] with other rough volatility models, but it also leads to a pathwise MDP.

With Y := log(V ), the system (X, Y ) fits into (4.1) where

K 1 (t) = t 2H-1 , K 2 (t) = t H-1 2 /Γ(H + 1
2 ), for H ∈ (0, 1 2 ), Σ(y) = exp(y), b(y) = -a/(2H), a > 0, ζ(y) ≡ 1 and y 0 = log(V 0 ). The bound (4.3) is then satisfied since t 0 (t -s) H-1 2 dW s is a Gaussian process with exponentional moments bounded uniformly in t ∈ T, and for each v ∈ A N , N > 0:

t 0 (t -s) H-1 2 v s ds ≤ N t 2H 2H ,
almost surely, by Cauchy-Schwarz inequality. Therefore sup t∈T E exp Y ε,v t is finite, yielding the claim. Moreover, the volatility equation is explicit so we shall not be concerned with uniqueness and the rest of Assumptions 4.2 and 4.5 is straightforward to check. This implies that Propositions 4.3 and 4.6 apply, and so does Corollary 4.4. Again, Theorem 3.29 guarantees that the LDP and MDP still hold when K 2 is replaced with the non-convolution fractional Brownian motion kernel. Gassiat [START_REF] Gassiat | On the martingale property in the rough Bergomi model[END_REF] showed that, if ρ ≤ 0, then the stock price process is a true martingale, ensuring that Assumption 4.7 holds, and implied volatility asymptotics thus follow from Corollary 4.8. [START_REF] Euch | The characteristic function of rough Heston models[END_REF] the rough Heston model fits into the framework of (4.1) with

Rough Heston. As introduced in

K 1 (t) = K 2 (t) = t H-1 2 /Γ(H+ 1 2 ), for H ∈ (0, 1 2 ), y 0 > 0, Σ(y) = y, b(y) = κ θ-y , κ > 0, θ ≥ 0 and ζ(y) = ξ √ y, ξ ∈ R d .
Linear growth and local Hölder continuity of the coefficients clearly hold.

The weak existence and uniqueness was proved in [START_REF] Jaber | Affine Volterra processes[END_REF], however the square-root coefficient brings an issue for pathwise uniqueness of the SVE. We assume here that there exists a set U of coefficients (H, κ, θ, ξ, ρ, y 0 ) such that pathwise uniqueness indeed stands. The only known result so far is due to [START_REF] Yamada | On the uniqueness of solutions of stochastic differential equations[END_REF] in the smooth case H = 1 2 . We also recall that pathwise uniqueness was proved for ζ(y) = y γ where γ > 1 2H+1 in [START_REF] Mytnik | Uniqueness for Volterra-type stochastic integral equations[END_REF], but does not encompass the square root case. Therefore Assumption 4.2 holds in those two cases. On a heuristic note remark that, even if pathwise uniqueness fails, there is a unique candidate for G ε since there exists a unique strong solution until the first hitting time of zero. The issue is it may not satisfy the SVE anymore after that time, but should be consistent for small-time LDP. Moreover, uniqueness of the limit equation (4.4) only holds up to first hitting time of zero. Hence we will make use of the uniqueness relaxation presented in Section 3.1 and similar arguments as in Example 3.8 to prove that Assumption 3.3 holds. The suggested rate function (3.2) reads now

I(φ, ϕ) = inf 1 2 t 0 u 2 t + v 2 t dt : u, v ∈ L 2 , φ t = t 0 √ ϕ s ρu s + ρv s ds, ϕ t = I H+ 1 2 y0 (ξ √ ϕv) .
We emphasise that ϕ above solves the Volterra equation (4.5)

ϕ t = y 0 + ξ t 0 (t -s) H-1 2 Γ(H + 1 2 )
√ ϕ s v s ds, for all t ∈ T.

Lemma 4.9. Let v ∈ L 2 and ϕ satisfying the Volterra equation (4.5) with H ∈ (0, 1 2 ). Then the set

D := {t ∈ [0, T ] : ϕ t > 0} has Lebesgue measure T .
Proof. We follow some arguments in the proof of [START_REF] Jaber | Affine Volterra processes[END_REF]Theorem 3.6]. Let us drop the subscript in the kernel and write it K for clarity, and introduce its resolvent of the first kind L( dt) := t -H-1 2 Γ( 12 -H) dt [53, Definition 5.5.1]. Moreover, for h > 0, define ∆ h K(t) := K(t+h) and for every measurable function f on R + and measure g on R + , (f * g)(t) := t 0 f (t-s)dg(s). It is proved, in [1, Equation (3.9)], that ∆ h K * L is non-decreasing but in fact in this special (rough) case, it is strictly increasing. Indeed, the authors show that in the general case, for all 0 ≤ s ≤ t ≤ T ,

(∆ h K * L)(t) -(∆ h K * L)(s) = h 0 K(h -u) L(s + du) -L(t + du) ,
which is positive because K > 0 and L is decreasing. Furthermore K is decreasing and L > 0 thus

0 < (∆ h K * L)(t) < (K * L)(t) = 1,
where the equality holds by definition. Let

ϕ t = y 0 + t 0 K(t -s)ξ √ ϕ s v s ds =: y 0 + (K * z)(t)
, where z is trivially a semimartingale, hence from [1, Equation (2.15)]:

y 0 + (∆ h K * dz)(t) = 1 -(∆ h K * L)(t) y 0 + ∆ h K * L (0)ϕ t + d(∆ h K * L) * ϕ (t),
which is strictly positive because because y 0 , ϕ ≥ 0 and the two lines before. Now let us suppose there exists an interval [t, t + h] ⊂ T on which ϕ = 0. Then

ϕ t+h = y 0 + t 0 K(t + h -s) √ ϕ s v s ds = y 0 + (∆ h K * dz)(t) > 0,
which is a contradiction. Hence no such interval exists and the claim follows.

Remark 4.10. This argument works for any rough kernel but not for the diffusion case H = 1 2 . We refer to [START_REF] Donati-Martin | Large deviations for squares of Bessel and Ornstein-Uhlenbeck processes[END_REF]Proposition 3.3] for the latter.

The previous lemma allows to invert the integrals as showed in the proof of Corollary 4.4 and yields a more explicit form for I:

(4.6) I(φ, ϕ) = T 0 ½ ϕt>0 2ρ 2 ϕ t   φ2 t -2ρ φt D H+ 1 2 (ϕ -y 0 )(t) ξ + D H+ 1 2 (ϕ -y 0 )(t) ξ 2   dt, if the integral is well defined, φ ∈ AC 0 , ϕ ∈ I H+ 1 2 y0
(L 1 ), and I = +∞ otherwise. We can now prove the following: Lemma 4.11. Let H ∈ (0, 1 2 ), then the functional I satisfies Assumption 3.3.

Proof. Note that any solution ϕ of (4.5) is non-negative. Let (φ, ϕ) be such that I(φ, ϕ) is finite. Then, for each δ > 0, define ϕ δ t := ϕ t + δt H+ 1 2 such that ϕ δ is strictly positive. Therefore from definition (3.21) we have:

D H+ 1 2 (ϕ δ -y 0 )(t) -D H+ 1 2 (ϕ -y 0 )(t) = 1 Γ( 1 2 -H) d dt t 0 (t -s) -H-1 2 (ϕ δ s -ϕ s ) ds = δ Γ( 1 2 -H) d dt t 0 (t -s) -H-1 2 s H+ 1 2 ds = δ π H + 1 2 Γ( 1 2 -H) cos(πH)
, which entails convergence as δ goes to zero, uniformly on T. Now define the control

v t := D H+ 1 2 (ϕ -y 0 )(t) √ ϕ t ½ ϕt>0 ,
and v belongs to L 2 since I(φ, ϕ) is finite. Then for each δ > 0 the control v δ defined as

v δ t := D H+ 1 2 (ϕ δ -y 0 )(t) ϕ δ t ≤ D H+ 1 2 (ϕ δ -y 0 )(t) √ ϕ t
is also in L 2 because ½ ϕt>0 = 1 almost everywhere. Furthermore, for all t ∈ D, lim δ↓0 ϕ δ t -1 = (ϕ t ) -1 and therefore lim δ↓0 v δ t = v t . Let P (ϕ) denote the term between brackets in (4.6) divided by 2ρ 2 , which is non-negative by design since it corresponds to ϕ(u 2 +v 2 ). Therefore, by Lemma 4.9,

I(φ, ϕ) -I(φ, ϕ δ ) = D 1 ϕ t - 1 ϕ δ t P (ϕ) t dt + D 1 ϕ δ t P (ϕ) t -P (ϕ δ ) t dt,
where the first integrand is smaller than P (ϕ) t /ϕ t for all t ∈ D and this upper bound belongs to L 1 by assumptions. Hence the dominated convergence theorem implies that the first integral goes to zero. From the calculations above we deduce that P (ϕ) t -P (ϕ δ ) t tends to zero uniformly as δ goes to zero, hence the second integrand converges pointwise and, for δ small enough, is dominated by P (ϕ)/ϕ. A second application of DCT yields convergence of the integral, and the claim follows.

Therefore the large and moderate deviations from Propositions 4.3 and 4.6 apply if the coefficients belong to U and H ∈ (0, 1 2 ). Observe that Proposition 4.6(M3) agrees with [42, Section 3.5], although the routes taken differ significantly. El Euch and Rosenbaum [START_REF] Euch | Perfect hedging in rough Heston models[END_REF]Appendix B] showed that Assumption 4.7 is satisfied, and the implied volatility behaviour thus follows from Corollary 4.8.

Multi-factor rough

Bergomi. Let W be an R m+1 -Brownian motion, Z := Z (1) , • • • , Z (m) † where Z (j) t := t 0 K (j) (t -s)dW (j)
s , and we allow

K (j) ∈ L 2 (T, R + ), j ∈ 1, m , to be homogeneous of different degrees H j -1 2 with H j ∈ (0, 1]. Therefore, the variance of Z is proportional to A t := t 2H1 /(2H 1 ), . . . , t 2Hm /(2H m ) † ,
for all t ∈ T. Assume without loss of generality that the H j are ordered by increasing values, then we will design the rescaling at the speed ε -2H1 . Denote m ⋆ := max{j ∈ 1, m : H j = H 1 }. Let U and V be m-dimensional square matrices, and Y an m-dimensional process defined for all t ∈ T by

Y t := y 0 + UZ t - 1 2 VA t , y 0 ∈ R m .
The log-price reads

X t = - 1 2m t 0 m i=1 exp Y (i) s ds + t 0 1 m m i=1 exp Y (i) s dB s , where B = ρW (m+1) + m j=1 ρ j W (j) , ρ 2 + m i=1 ρ 2 i = 1. The rescaling X ε t = ε H1-1 2 X εt , Y ε t := Y εt yields              X ε t = - ε H1+ 1 2 2m t 0 m i=1 exp Y ε,(i) s ds + ε H1 t 0 1 m m i=1 exp Y ε,(i) s dB s , Y ε,(i) t = y (i) 0 + m j=1 ε Hj U ij Z (j) t -V ij (εt) 2Hj 2H j , for all i ∈ 1, m .
As we will shift each BM by ε -H1 v (j) , we notice that ε Hj -H1 U ij goes to zero if H j > H 1 , i.e. if j > m ⋆ . It means that the roughest component(s) (the one(s) with H 1 ) will outweigh the others, and only the former will make a contribution to the rate function.

Although similar to its one-dimensional counterpart, this model does not fit into the framework of (4.1). Regarding the assumptions of Theorems 3.11, we only check H3b and Assumption 3.9 because the others are standard and similar to the one-dimensional case. Clearly (Y ε, (1)

, • • • , Y ε,(m)
) is an autonomous subsystem. As a Gaussian process, UZ has exponential moments of all orders and for all N > 0, j ∈ 1, m and v (j) ∈ A N :

t 0 K (j) (t -s)v (j) s ds ≤ √ N K (j) 2 almost surely, thus exp Y ε,(i),v ∈ L p (Ω) for all p ≥ 1.
Therefore the bound (3.6) and Assumption 3.9 are satisfied. This estimate also checks that (3.17) and thus H8 stand. Since ϑ ε = ε H1 , we define for the moderate deviations regime

h ε = ε β , β ∈ (0, H 1 ).
Corollary 4.12. The pathwise LDP and MDP hold.

• (X ε , Y ε ) ∼ LDP I, ε -2H1
where for all φ ∈ W and ϕ ∈ W m :

I(φ, ϕ) = inf 1 2 T 0 u 2 t + |v t | 2 dt : u ∈ L 2 (T, R), v ∈ L 2 (T, R m ), φ t = t 0 1 m m i=1 e ϕ (i) s 1 2   ρu s + m j=1 ρ j v (j) s   ds, ϕ (i)
t = y (i) 0 + m ⋆ j=1 U ij t 0 K (j) (t -s)v (j) s ds . • (X ε , Y ε ) ∼ MDP β Λ, ε -2β
where for all φ ∈ W and ϕ ∈ W m :

Λ(φ, ϕ) = inf 1 2 T 0 u 2 t + |v t | 2 dt : u ∈ L 2 (T, R), v ∈ L 2 (T, R m ), (4.7)
φ t = t 0 1 m m i=1 e y (i) 0 1 2   ρu s + m j=1 ρ j v (j) s   ds, ϕ (i) t = m ⋆ j=1 U ij t 0 K (j) (t -s)v (j)
s ds .

Proof. The LDP is a direct application of Theorem 3.11 and the MDP of Theorem 3.20.

One can also recover the LDP and MDP for (X ε ) ε>0 as well as the small-time LDP and MDP by contraction principle, as in Propositions 4.3 and 4.6.

If U is lower triangular (i.e. U ij = 0 for all i < j), for instance if it arises from the Cholesky decomposition of a covariance matrix, then for all φ ∈ W, ϕ ∈ W m , one can derive the vector v recursively, followed by u. Note that if m ⋆ < m, ϕ may not be attainable by the restrained number of controls {v (j) , j ∈ 1, m ⋆ }. Example 4.13. Consider the case m = 2. Let K (1)

(t) = K (2) (t) = t H-1 2 /Γ(H + 1 2 ) for H ∈ (0, 1 2 )
, hence m ⋆ = 2, and U be lower triangular (i.e. U 12 = 0). Then in the moderate deviations setting v and u are explicit from (4.7): 2) .

v (1) = 1 U 11 D H+ 1 2 ϕ (1) , v (2) = 1 U 22 D H+ 1 2 ϕ (2) - U 21 U 22 v (1) , u = 1 ρ √ 2 φ exp y (1) 0 + exp y (2) 0 -1 2 -ρ 1 v (1) -ρ 2 v (
Remark 4.14. We can similarly consider multidimensional versions of the other models presented in this chapter and derive large and moderate deviation principles. We only work out the computations for the multi-factor rough Bergomi model because it is the most relevant in the literature.

4.3. Tail rescaling. We now investigate tail rescalings, which generally have the form X ε = εX, such that an LDP provides asymptotic estimates on P(X ε ≥ 1) = P(X ≥ ε -1 ). The MDP for the whole system is not available in this case because Y := lim ε↓0 Y ε ≡ 0 hence the limit equation for X ε , arising from (3.16), would be independent of the control. Note that the theory does not break down but the rate function is trivial (equals zero at zero and +∞ everywhere else). Furthermore, the exponential function prevents the study of such a rescaling in the rough Bergomi model.

:= ε 2 X t , the system becomes        X ε t = - 1 2 t 0 (Y ε s ) 2 ds + ε t 0 Y ε s dB s , Y ε t = εy 0 + t 0 κ (εθ -Y ε s ) ds + ε t 0 ξ (t -s) H-1 2 Γ(H + 1 2 ) dW s , (4.8)
where the coefficients are identical to the small-time case. Although the rescaling is different, Assumption 2.3 and H1 -H4 are easily satisfied in a similar way, the limit equation (3.5) has a unique solution, and therefore Theorem 3.11 applies.

Corollary 4.15. The following hold:

(L1) (X ε , Y ε ) ∼ LDP I, ε -2 with I(φ, ϕ) = 1 2 T 0 u 2 t + v 2 t dt, where        u = 1 ρ φ ϕ + 1 2 ϕ -ρv ½ ϕ =0 , v = 1 ξ D H+ 1 2 (ϕ) + κI 1 2 -H (ϕ) , if φ ∈ AC 0 , ϕ ∈ I H+ 1 2 0
and infinity otherwise.

(L2) X ε ∼ LDP I X , ε -2 with I X (φ) = inf I(φ, ϕ) : ϕ ∈ I H+ 1 2 0 if φ ∈ AC 0 and infinity otherwise. (L3) For each t ∈ T, ε 2 X t ∼ LDP I X t , ε -2 , where I X t (x) = inf I X (φ) : φ t = x .
Proof. For (L1), Theorem 3.11 entails that the rate function is

I(φ, ϕ) = inf 1 2 T 0 u 2 t + v 2 t dt : u, v ∈ L 2 , φ t = - 1 2 t 0 ϕ 2 s ds + t 0 ϕ s ρu s + ρv s ds, ϕ t = - t 0 κ ϕ s ds + t 0 ξ (t -s) H-1 2 Γ(H + 1 2 ) v s ds .
Inverting the integrals as in Corollary 4.4 to obtain the unique controls and using D H+ 1 2 I 1 = I 1 2 -H yields the claim. Similarly to the small-time case, (L2) follows from the contraction principle, and one only needs to fix t ∈ T to prove (L3).

One can prove an LDP for Y ε in a similar way; a more interesting problem is the moderate deviations setting. Recall that an MDP for the couple (X ε , Y ε ) would have a trivial rate function because the limit equation of X ε is independent of the control. However, since the diffusion coefficient of Y ε is constant equal to ξ, one can obtain an MDP for Y ε . More surprisingly, the limit equations in the large deviations (3.5) and moderate deviations (3.16) regimes coincide, which leads to identical rate functions. Notice that ϑ ε = ε in this example, and therefore let h ε = ε -β where β ∈ (0, 1).

Corollary 4.16. Y ε ∼ MDP β Λ Y , ε -2β where Λ Y (ϕ) = 1 2ξ 2 T 0 D H+ 1 2 (ϕ)(t) + κI 1 2 -H (ϕ)(t) 2 dt, if ϕ ∈ I
Λ Y (ϕ) = inf 1 2 T 0 v 2 t dt : v ∈ L 2 , ϕ t = - t 0 κϕ s ds + ξ t 0 (t -s) H-1 2 Γ(H + 1 2 ) v s ds .
Inverting the integral yields the claim.

Rough Heston. After the rescaling

Y ε t := ε 2 Y t and X ε t := ε 2 X t , this model introduced in Section 4.2.3 takes the form        X ε t = - 1 2 t 0 Y ε s ds + ε t 0 Y ε s dB s , Y ε t = ε 2 y 0 + t 0 κ (t -s) H-1 2 Γ(H + 1 2 ) ε 2 θ -Y ε s ds + ε t 0 ξ (t -s) H-1 2 Γ(H + 1 2 )
Y ε s dW s .

Clearly H1 -H3 hold and we recall that U is the set of coefficients such that pathwise uniqueness, and hence H4, hold. We appeal to the uniqueness relaxation in the same way as the small-time case to prove the following result, which extends [27, Theorem 1.1] to the rough case.

Corollary 4.17. If the rough Heston coefficients belong to U and H ∈ (0, 1 2 ), then the following hold

(L1) (X ε , Y ε ) ∼ LDP I, ε -2 where I(φ, ϕ) = 1 2 T 0 u 2 t + v 2 t dt, where          u = 1 ρ φ √ ϕ + 1 2 √ ϕ -ρv ½ ϕ>0 , v = 1 ξ √ ϕ D H+ 1 2 (ϕ) + κϕ ½ ϕ>0 ,
where φ ∈ AC 0 and ϕ ∈ I (L 1 ) and infinity otherwise.

(L2) X ε ∼ LDP I X , ε -2 with I X (φ) = inf I(φ, ϕ) : ϕ ∈ I H+ 1 2 0 if φ ∈ AC 0 and infinity otherwise. (L3) For each t ∈ T, ε 2 X t ∼ LDP I X t , ε -2
, where I X t (x) = inf I X (φ) : φ t = x . Proof. The proof is similar to the small-time case. The potential rate function for the couple is

I(φ, ϕ) = inf 1 2 T 0 u 2 t + v 2 t dt : u, v ∈ L 2 , φ t = - 1 2 t 0 ϕ s ds + t 0 √ ϕ s ρu s + ρv s ds, ϕ t = t 0 (t -s) H-1 2 Γ(H + 1 2 ) -κϕ s + ξ √ ϕ s v s ds .
The same arguments that were used to prove Lemmata 4.9 and 4.11 in the small-time case can be applied again here. They entail that Assumption 3.

lim k↑∞ 1 k log P(X t ≥ k) = -inf y≥1 I X t (y).
In the Black-Scholes model with constant volatility σ > 0, we can directly compute

lim k↑∞ 1 k 2 log P(X t ≥ k) = -
= E F • G ε W + ε -1 • 0 v ε t dt
. Theorem 3.5(i) entails the existence of φ ∈ G 0 v,N such that the subsequence y : R + → R defined as y(ε) := inf α<ε x(α) has a subsequence converging to E [F (φ)]. By definition lim inf ε↓0 x ε = lim ε↓0 y ε , which implies that y has a limit in [-∞, +∞] and by uniqueness this limit must be E [F (φ)]. Therefore we deduce:

lim inf ε↓0 E 1 2 T 0 |v ε t | 2 dt + F • G ε W + ε -1 • 0 v ε t dt ≥ E 1 2 T 0 |v t | 2 dt + F (φ) ≥ inf 1 2 T 0 |v t | 2 dt + F (φ) : v ∈ L 2 such that φ ∈ G 0 v ∩ W d ,
which suggests the potential rate function I defined in (3.2) and concludes the proof of the lower bound.

Then we prove the Laplace principle upper bound, for all F ∈ C b (W d : R):

lim sup ε↓0 -ε 2 log E e -F •G ε (W )/ε 2 ≤ inf ψ∈W d {I(ψ) + F (ψ)}.
Assume that the right-hand side is finite otherwise there is nothing to prove. Fix ǫ > 0 and let φ ∈ W d such that

I(φ) + F (φ) ≤ inf ψ∈W d {I(ψ) + F (ψ)} + ǫ.
Since F is continuous at φ, there exists δ ∈ (0, ǫ) such that |F (φ) -F (ϕ)| ≤ ǫ for all ϕ ∈ W d such that φ -ϕ T ≤ δ. If φ is uniquely characterised then the proof is the same as in [START_REF] Budhiraja | A variational representation for positive functionals of infinite dimensional Brownian motion[END_REF]. Otherwise, by Theorem 3.5(iii), we can choose φ δ uniquely characterised such that φ -φ δ T ≤ δ and I(φ) -I(φ δ ) ≤ δ, which implies I(φ) + F (φ) -I(φ δ ) -F (φ δ ) ≤ 2ǫ. Hence, combining inequalities we obtain

I(φ δ ) + F (φ δ ) ≤ inf ψ∈W d {I(ψ) + F (ψ)} + 3ǫ.
Moreover, there exist {v n } n∈N in L 2 such that (3.3) is satisfied with φ δ and m ≥ 1/ǫ such that

G 0 v m = {φ δ } and 1 2 T 0 |v m t | 2 dt ≤ I(φ δ ) + 1 m ≤ I(φ δ ) + ǫ,
and therefore the remainder of the upper bound proof unfolds identically.

Along the subsequence

{ε n } n≥0 , G εn (W + ε -1 n • 0 v m t dt
) converges in distribution to φ δ by item (i). Using the variational representation formula (1.3) and the convergence we obtain

lim sup n↑∞ -ε 2 n log E exp - F • G εn (W ) ε 2 n = lim sup n↑∞ inf v∈A E 1 2 T 0 |v t | 2 dt + F • G εn W + ε -1 n • 0 v t dt ≤ lim sup n↑∞ E 1 2 T 0 |v m t | 2 dt + F • G εn W + ε -1 n • 0 v m t dt = 1 2 T 0 |v m t | 2 dt + F (φ δ ) ≤ I(φ δ ) + F (φ δ ) + ǫ ≤ inf ψ∈Ω {I(ψ) + F (ψ)} + 4ǫ.
Since ǫ > 0 is arbitrary this concludes the proof.

:= b ε (s, X ε,v s ½ s≤τn )
and σ n s := σ ε (s, X ε,v s ½ s≤τn ). We start by assuming that all the coefficients satisfy the linear growth condition H3a. We fix n ∈ N and observe that, almost surely:

|X ε,v t | p ½ t≤τn ≤4 p-1 |X ε 0 | p + t 0 K(t -s)b n s ds p + t 0 K(t -s)σ n s v s ds p + ϑ p ε t 0 K(t -s)σ n s dW s p =:4 p-1 |X ε 0 | p + I n + II n + III n ,
because if t > τ n then the left-hand side is zero while the right-hand side is non-negative, and if t ≤ τ n then s ≤ τ n for all s ∈ [0, t] and the τ n dependence vanishes on both sides of the inequality.

For ε small enough we can bound |X ε 0 | by 2 |X 0 | and ϑ ε by 1 and we will do so repetitively in the sequel. Using Hölder's and Jensen's inequalities, we obtain the following estimates almost surely: (A.1)

I n ≤ t 0 |K(t -s)| 4 p |b n s | 2 ds p 2 t 0 |K(t -s)| 2-4 p ds p 2 ≤ t p 2 +4 2(p-2) K p-2 2 t 0 |K(t -s)| 2 |b n s | p ds, and
(A.2) II n ≤ N p 2 t 0 |K(t -s)| 2-4/p |K(t -s)| 4 p |σ n s | 2 ds p 2 ≤ N p 2 K p-2 2 t 0 |K(t -s)| 2 |σ n s | p ds,
where we also used that

T 0 |v s | 2 ds ≤ N almost surely. Notice that for fixed t ∈ T, { u 0 K(t - s)σ n s dW s : u ∈ [0, t]
} is a continuous local martingale and is bounded in L 2 (Ω). Hence, using Burkholder-Davis-Gundy (BDG) inequality and similar calculations as (A.2) there exists C p > 0 such that

E[III n ] ≤ C p K p-2 2 t 0 |K(t -s)| 2 E |σ n s | p ds.
From the linear growth condition on b ε and σ ε (uniform in ε > 0) we deduce that there exists C 1 > 0 independent of ε, v, n, t such that, for all n ∈ N,

f n t := E |X ε,v t | p ½ t≤τn satisfies the inequality (A.3) f n t ≤ C 1 + C 1 t 0 |K(t -s)| 2 f n s ds.
The following lemma (Lemma A.1) yields a uniform bound in both n ∈ N and t ∈ T for f n t . Taking the limit as n goes to infinity and using Fatou's lemma concludes the first part of the proof.

Lemma

f (t) ≤ c 1 -c 1 t 0 R(s) ds ≤ c 1 -c 1 T 0 R(s) ds,
where R is the (non-positive) resolvent of second kind of -K [1, Equation (2.11)], proving the lemma.

If only H3b and Assumption 3.9 hold with an autonomous sub-system Υ (see Definition 2.6), then by the previous calculations for all l ∈ Υ, the components (X ε,v ) (l) satisfy the bound (3.8) because their coefficients have linear growth. Then we turn our attention to the components (X ε,v ) (i) , i / ∈ Υ. Using (3.6) and Hölder's inequality as in (A.2), we obtain that for all 1 ≤ j ≤ d such that K ij = 0 and for all 1 ≤ k ≤ m:

E t 0 K ij (t -s)σ (jk) ε s, X ε,v s ½ s≤τn 2 ds p/2 ≤ C p Υ E t 0 |K ij (t -s)| 2 1 + X ε,v s ½ s≤τn Υ c + Γ (X ε,v s ½ s≤τn ) (Υ) 2 ds p/2 ≤ C p Υ E t 0 |K ij (t -s)| 2 1 + |X ε,v s ½ s≤τn | Υ c + Γ (X ε,v s ½ s≤τn ) (Υ) p ds t 0 |K ij (t -s)| 2 ds p/2-1 ≤ 3 p-1 C p Υ K p-2 2 t 0 |K ij (t -s)| 2 E 1 + X ε,v s ½ s≤τn p Υ c + Γ (X ε,v s ½ s≤τn ) (Υ) p ds ≤ C 2 + C 2 t 0 |K ij (t -s)| 2 E |X ε,v s ½ s≤τn | p Υ c ds, (A.4)
for some C 2 > 0. Applying the same calculations to the other terms and summing all the coefficients we fall back on (A.3). Taking the limit and applying Fatou's lemma again conclude the proof.

A.3. LDP tightness: Proof of Lemma 3.15. Let us fix p > 2 ∨ 2/γ, N > 0, a family {v ε } ε>0 in A N and ε > 0. For clarity we will write b u := b ε (u, X ε,v u ) and σ u := σ ε (u, X ε,v u ) for all u ∈ T. Then, for all 0 ≤ s < t ≤ T , using Cauchy-Schwarz and BDG inequalities as in the previous proof we obtain:

E X ε,v ε t -X ε,v ε s p ≤ 6 p-1 E s 0 K(t -u) -K(s -u) b u du p + 6 p-1 E t s K(t -u)b u du p + 6 p-1 N p/2 E s 0 K(t -u) -K(s -u) σ u 2 du p/2 + 6 p-1 N p/2 E t s |K(t -u)σ u | 2 du p/2 + 6 p-1 ϑ p ε C p E s 0 K(t -u) -K(s -u) σ u 2 du p/2 + 6 p-1 ϑ p ε C p E t s |K(t -u)σ u | 2 du p/2
.

In a first step we assume that all the coefficients satisfy the linear growth condition H3a. Analogous calculations to the proof of Lemma 3.13, bounds on sup t∈T,ε>0 E X ε,v ε t p , linear growth from H3a and Assumption 2.3 lead to

E X ε,v ε t -X ε,v ε s p ≤ C 1 s 0 K(t -u) -K(s -u) 2 du p/2 + t s K(t -u) 2 du p/2 = C 1 s 0 K(u + t -s) -K(u) 2 du p/2 + t-s 0 K(u) 2 du p/2 ≤ C 2 (t -s) γp/2 ,
for some C 1 , C 2 > 0 independent of ε, v ε , s, t. Again, if there are components such that only H3b holds with Assumption 3.9 then following the example of (A.4) yields the same result. Then Kolmogorov continuity theorem asserts that X ε,v ε admits a version which is Hölder continuous on T of any order α < γ/2 -1/p, uniformly in ε > 0 because C 2 does not depend on ε, and which satisfies (3.9). Furthermore, Aldous theorem [START_REF] Billingsey | Convergence of probability measures[END_REF]Theorem 16.10] states that the sequence {X ε,v ε } ε>0 is tight.

A.4. LDP compactness: Proof of Lemma 3.18. We prove that for all N > 0, the sublevel sets

L N := {φ ∈ W d : I(φ) ≤ N } of the map I : W d → R given by (3.
2) or more precisely by (A.5)

I(φ) = inf 1 2 T 0 |v s | 2 ds : v ∈ L 2 , φ t = x 0 + t 0 K(t -s) b(s, φ s ) + σ(s, φ s )v s ds
are compact. Fix N > 0 and consider an arbitrary sequence J := {φ n } n∈N ⊂ L N ; we will show that there exists a converging subsequence the limit of which belongs to L N . Interestingly enough, the proof parallels, in a deterministic context, the proofs of bound, Hölder continuity and convergence of X ε,v . Relative compactness. According to Arzelà-Ascoli's theorem, the family J is relatively compact in W d if and only if {φ n t } is bounded uniformly in n ∈ N and in t ∈ T and J is equicontinuous.

Moreover, for all n ∈ N and all t ∈ T, there exists

v n ∈ L 2 such that 1 2 T 0 |v n t | 2 dt ≤ N and φ n ∈ G 0 v n , which means v n ∈ S 2N and φ n t = x 0 + t 0 K(t -s) b(s, φ n s ) + σ(s, φ n s )v n s ds.
Hence Remarks 3.14 and 3.16 grant the uniform bound and equicontinuity respectively. Therefore J is relatively compact which entails that L N is relatively compact for any N > 0.

Closure. Let {φ n } n∈N be a converging sequence of L N and denote its limit by φ ∈ W d . The controls v n associated to φ n through (A.5) belong to S 2N which is a compact space with respect to the weak topology. Hence there exists a subsequence {n k } k∈N such that v n k converges weakly in L 2 to a limit v ∈ S 2N and lim k↑∞ φ n k = φ. Now let us prove that φ ∈ L N . For clarity we replace n k by n from now on. The convergence as n goes to +∞ and the continuity of the paths entail

sup n∈N sup t∈T |φ n t | + |φ t | < +∞,
such that the paths lie in compact subsets of R d and H2 asserts that uniform continuity of the coefficients b and σ hold. Therefore they admit continuous moduli of continuity that we respectively name ρ b and ρ σ . Using Cauchy-Schwarz inequality and H3 we get that for all t ∈ T:

t 0 K(t -s)b(s, φ n s )ds - t 0 K(t -s)b(s, φ s )ds ≤ K L 1 ρ b |φ n • -φ • | T , t 0 K(t -s)σ(s, φ n s )v n s ds - t 0 K(t -s)σ(s, φ s )v s ds ≤ t 0 K(t -s) σ(s, φ n s ) -σ(s, φ s ) v n s ds + t 0 |K(t -s)σ(s, φ s )(v n s -v s )| ds ≤ K 2 v n 2 ρ σ (|φ n • -φ • |) T + σ(φ) T t 0 |K(t -s)(v n s -v s )| ds,
and both estimates converge towards zero as n tends to infinity. Therefore, for all t ∈ T

φ t = lim n↑∞ φ n t = lim n↑∞ x 0 + t 0 K(t -s) b(s, φ n s ) + σ(s, φ n s )v n s ds = x 0 + t 0 K(t -s) b(s, φ s ) + σ(s, φ s )v s ds, so that φ ∈ L N since v ∈ S 2N ,
E |η ε,v t | p ≤ 5 p-1 |X ε 0 -x 0 | p ϑ ε h ε p (B.1) + 5 p-1 E t 0 K(t -s) b ε s, Θ ε,v s -b s, Θ ε,v s ϑ ε h ε ds p + 5 p-1 E t 0 K(t -s) b s, Θ ε,v s -b(s, X s ) ϑ ε h ε ds p + 5 p-1 N p/2 E t 0 K(t -s)σ ε s, Θ ε,v s 2 ds p/2 + 5 p-1 C p h p ε E t 0 K(t -s)σ ε s, Θ ε,v s 2 ds p/2
.

The first term converges by H7 and is thus bounded. Notice that H8 entails

E t 0 K(t -s) b ε s, Θ ε,v s -b s, Θ ε,v s ϑ ε h ε ds p ≤ ν ε ϑ ε h ε p E t 0 K(t -s)Ξ Θ ε,v s ds p ,
which is bounded because ν ε (ϑ ε h ε ) -1 tends to zero and, using Cauchy-Schwarz and Jensen's inequalities in the same way as (A.1) and the bound (3.17), In the first four terms, Cauchy-Schwarz inequality allows to separate the kernels from the random variables. For the last four terms, analogous calculations to (A.2) achieve a similar separation of kernels and random variables. Then linear growth or (3.9) and bounds on sup t∈T,ε>0 E Ξ Θ ε,v ε t p and sup t∈T,ε>0 E η ε,v ε t p lead to the existence of C 1 > 0 independent of t and ε such that

E η ε,v ε t -η ε,v ε s p ≤8 p-1 E ν ε ϑ ε h ε s 0 K(t -u) -K(s -u) Ξ Θ ε,v ε u du p + 8 p-1 E ν ε ϑ ε h ε t s K(t -u)Ξ Θ ε,v ε u du p + 8 p-1 E s 0 K(t -u) -K(s -u) C b η ε,v ε u du p + 8 p-1 E t s K(t -u)C b η ε,v ε u du p + 8 p-1 N p/2 E s 0 K(t -u) -K(s -u) σ ε u, Θ ε,v ε
E η ε,v ε t -η ε,v ε s p ≤ C 1 s 0 K(t -u) -K(s -u)
E η ε,v ε t -η ε,v ε s p ≤ C 2 (t -s) γp/2 .
Then Kolmogorov continuity theorem asserts that η ε,v ε admits a version which is Hölder continuous on T of any order α < γ/2 -1/p, uniformly in ε > 0 and which satisfies (3.20). Furthermore, Aldous theorem [START_REF] Billingsey | Convergence of probability measures[END_REF]Theorem 16.10] states that the sequence {η ε,v ε } ε>0 is tight.

B.3. MDP weak convergence: Proof of Lemma 3.23. We have shown in Lemma 3.22 that for any subsequence {ε k } k∈N , {η ε k ,v ε k } k∈N and {v ε k } k∈N are tight as families of random variables with values in W d and S N respectively. By Skorohod representation theorem we can work with almost sure convergence for the purpose of identifying the limit. Hence there exists a subsubsequence, denoted hereafter η k , v k , that converges almost surely in the product topology on W d × S N to some W d × S N -valued limit (η 0 , v) in a possibly different probability space (Ω 0 , F 0 , P 0 ) as n tends to +∞. We also denote ε k , b k , σ k , X k 0 , Θ k along this subsequence. The convergence of the couple also takes place in distribution, and we follow the same method as in the LDP case which comes from [START_REF] Chiarini | On large deviations for small noise Itô processes[END_REF]. For all t ∈ [0, T ], let Ψ t : S N × W d → R such that Ψ t (f, ω) := ω t -t 0 K(t -s) ∇b(s, X s )ω s + σ(s, X s )f s ds ∧ 1.

Clearly, Ψ t is bounded and one can show its continuity along the same lines as in the LDP proof but in a simpler way because of the linearity. Therefore

lim k↑∞ E Ψ t (v k , η k ) = E 0 Ψ t (v, η 0 ) ,
and we prove that the left-hand side is actually equal to zero. By H5 and Taylor's formula there exists a sequence of R d -valued stochastic processes {R ε } ε>0 such that b s, X s + ϑ ε h ε η ε,v s -b(s, X s ) ϑ ε h ε = ∇b s, X s η ε,v s + R ε (s), (B.4) and a constant C R > 0 such that

|R ε (s)| ≤ C R ϑ ε h ε |η ε,v s | 2 .
We recall that ∇b(•, X) T is finite by H5 and observe that, by (3.19):

(B.5)

E |R ε (u)| p ≤ C R ϑ ε h ε p E η ε,v ε s 2p < ∞.
Again starting from (3.15), we use H8, the Taylor estimate (B.4) and Itô isometry to get

E Ψ t (η k , v k ) 2 ≤ 5 X k 0 -x 0 2 ϑ 2 ε k h 2 ε k + 5 ν ε k ϑ ε k h ε k 2 E t 0 K(t -s)Ξ Θ k s ds 2 + 5E t 0 K(t -s)R ε (s) ds 2 + 5E t 0 K(t -s) σ k s, Θ k s -σ s, X s v k s ds 2 + 5 h 2 ε k E t 0 K(t -s)σ k s, Θ k s 2 ds =:5 I k + II k + III k + IV k + h -2 ε k V k .
H7 and H8 tell us that I k +II k tends to zero as ε goes to zero while an application of Cauchy-Schwarz inequality and the bound (B.5) yields the same conclusion for III k .

To deal with IV k , recall that Θ k converges in distribution towards X as k tends to infinity. Since X is deterministic, the convergence actually takes place in probability, with respect to the topology of uniform convergence. Moreover Θ k is uniformly bounded in L r (Ω) for all r > p, thus the family { Θ k p } k≥0 is uniformly integrable. Therefore the convergence also occurs with respect to the L p (Ω)-norm.

The modulus of continuity of σ is only available on compact sets of T × R d so we define a constant M > X T and introduce the following sets, for each k ∈ N:

E k := ω ∈ Ω : Θ k (ω) -X T ≤ M ,
with the observation that lim k↑∞ P(E k ) = 1 thanks to the previous argument. Since X is uniformly bounded, Θ k t (ω) ≤ 2M for all t ∈ T, ω ∈ E k , k ≥ 0. Therefore using Cauchy-Schwarz inequality,

E [IV k ] ≤ 2N E t 0 |K(t -s)| 2 σ k s, Θ k s -σ s, Θ k s 2 ds(½ E k + ½ E c k ) + 2N E t 0 |K(t -s)| 2 σ s, Θ k s -σ(s, X s ) 2 ds(½ E k + ½ E c k ) ,
where we will use the localisation to obtain convergence in the first term and Hölder continuity in the second. Let us assume for the moment that linear growth H3a holds. We use that Θ k is uniformly bounded by 2M in E k and linear growth for both σ k and σ to obtain

sup s∈T E σ k s, Θ k s -σ(s, Θ k s ) 2 = sup s∈T E σ k s, Θ k s -σ(s, Θ k s ) 2 ½ E k + sup s∈T E σ k s, Θ k s -σ(s, Θ k s ) 2 ½ E c k ≤ σ k -σ 2M 2 T + sup s∈T E ½ E c k C 2 L 2 + 2|Θ k s | 2 , (B.6)
which tends to zero as k goes to infinity because of H2 for the first term and because P(Ω \ E k ) tends to zero for the second. Moreover, by H6, there exists δ > 0 such that σ is locally δ-Hölder continuous thus there exist C 2M > 0 such that

sup s∈T E ½ E k σ s, Θ k s -σ(s, X s ) 2 ≤ sup s∈T E ½ E k C 2M Θ k s -X s 2δ ,
which tends to zero. Finally, linear growth leads to (B.7)

sup s∈T E ½ E c k σ s, Θ k s -σ(s, X s ) 2 ≤ E ½ E c k C 2 L 2 + |Θ k s |+ X s 2 ,
which also tends to zero because P(Ω\E k ) tends to zero. If only H3b with Assumption 3.9 hold then a different bound depending on (2.2) would replace those in (B.6) and (B.7), by noticing that X = G 0 (0). In both cases the above estimates tend to zero as k tends to infinity, hence E [IV k ] converges towards zero. Finally, {V k } k∈N is uniformly bounded across k ≥ 0 as (B.3) shows, thus h -p ε k V k tends to zero. We have proved that lim k↑∞ E Ψ t (v k , η k ) = 0, and this entails that the limit η 0 satisfies (3.16) P 0 -almost surely, for all t ∈ T. Since η 0 has continuous paths, this holds for all t ∈ T, P 0 -almost surely and the solution is unique therefore we conclude that η 0 = ψ. Every subsequence has a subsequence for which this convergence holds therefore η ε,v ε converges weakly towards ψ as ε goes to zero.

2. GENERAL FRAMEWORK 2 . 1 .

 21 Notations. We consider a fixed time horizon T > 0, and denote T := [0, T ], R + := [0, +∞) and R + := [0, +∞]. For d 1 ≥ 1, d 2 ≥ 1, |•| denotes the Euclidean norm in R d1 and the Frobenius norm in R d1×d2 , and d 1 , d 2 := {d 1 , • • • , d 2 }. For p ≥ 1, L p stands short for L p (T), and • 2 is the usual L 2 norm. Furthermore, for d ≥ 1, W d := C(T : R d ) represents the space of continuous functions from T to R d , equipped with the supremum norm ϕ T := sup t∈T |ϕ t | for any ϕ ∈ W d . Finally, for any

  d and converge pointwise to b and σ as ε goes to zero. Moreover, b(t, •) and σ(t, •) are continuous on R d , uniformly in t ∈ T. H3. Either a) or b) holds: a) For all ε > 0 small enough, b ε and σ ε have linear growth uniformly in ε and in t ∈ T.

Definition 3 . 1 .

 31 For all N ∈ N and v ∈ A N we define G 0 v,N := φ : Ω → W d there exist {ε n } n∈N ⊂ R + with lim n↑∞ ε n = 0 and a sequence {v εn } n∈N ⊂ A N with lim n↑∞ v εn = v in distribution such that

32 ,Theorem 3 . 27 .

 32327 Theorem 4.2]. Decreusefond's main result yields the Hölder regularity of the stochastic Volterra integral [32, Theorem 3.1]: Let Assumption 3.25 hold and u

3. 5 . 2 .

 52 Large and moderate deviations. For each ε > 0 consider the stochastic Volterra equation(3.24)

Theorem 3 . 29 .

 329 The conclusions of Theorems 3.11 and 3.20 stand for {X ε } ε>0 defined by(3.24) when each component of the kernel K satisfies either Assumption 2.3 or Assumptions 3.25 and 3.28.

  ζ and b are continuous on R; • b and ζ are of linear growth; • Σ is either of linear growth or such that for all p ≥ 1, N > 0 and ε > 0 small enough, (4.3) sup t∈T,v∈AN

Proposition 4 . 3 (

 43 Large deviations). Under Assumptions 3.3, 4.1 and 4.2, the following hold:

Corollary 4 . 4 .

 44 Under the same assumptions as Proposition 4.3, if ρ = 0 or ζ(ϕ t ) = 0 almost everywhere, the rate function I can be written

( 1 )

 1 ε (t, (x, y)) -b (1) (t, (x, y)) ≤ ε H+ 1 2 |Σ(y)| and b

Proposition 4 . 6 .

 46 Under Assumptions 4.1, 4.2, 4.5 and the condition ζ(y 0 )Σ(y 0 ) = 0, the following moderate deviations hold:

  ∈ AC , which can be solved as a variational problem as in [59, Corollary 2.4]. The corresponding Euler-Lagrange equation reads ψ = ρ φ/ Σ(y 0 ) hence ψ = ρ φ/ Σ(y 0 ) because ψ0 = 0 by definition. Plugging into the above equation finishes the proof. (M3) The rate function is given by contraction principle as

  3 and 4.6 in a similar fashion. The explicit form of the MDP rate function allows a closed form expression. Corollary 4.8. Let Assumption 4.7 hold. (LDP) Under the same assumptions as Proposition 4.3,

1 2

 1 -β 2 = Σ(y 0 ).Proof.(LDP) Consider the case k > 0. Proposition 4.3(L3) translates into

H+ 1 2 0

 2 and infinity otherwise. Proof. From (4.8), b ε (y) = κ(εθ -y) converges to b(y) = -κy and the diffusion coefficient is constant, hence H2 -H6 are easily satisfied. Moreover, b ε -b ≡ εκθ and ε 1-(H-β) tends to zero therefore H7 and H8 also hold. Theorem 3.20 thus yields an MDP with rate function

1 2σ 2 1 . 5 .

 215 t , and, similarly to the small-time case, the proof follows from [49, Corollary 7.1]. APPENDIX A. TECHNICAL LARGE DEVIATIONS PROOFS A.Abstract relaxation: Proof of Theorem 3.The proof follows [13, Theorem 4.4]. The lower bound proof stands as it is until the last series of inequalities. For any sequence {v ε } ε>0 in A N , N > 0 converging weakly to v and any F ∈ C b (W d : R), define the function x : R + → R by x(ε) :

A. 2 .

 2 LDP moment bounds: Proof of Lemma 3.13. Let us fix p ≥ 2, N > 0, v ∈ A N , ε > 0 and t ∈ T. Let τ n := inf{t ≥ 0 : |X ε,v t | ≥ n} ∧ T for all n ∈ N. For clarity we write b n s

C 1 ,

 1 where C 1 is a positive constant that does not depend on ε. Since b is globally Lipschitz continuous, there exists C b > 0 such that for all s ∈ T:(B.2) b s, X s + ϑ ε h ε η ε,v s -b(s, X s ) ≤ C b ϑ ε h ε |η ε,v s | .Therefore, using Cauchy-Schwarz and Jensen's inequalities again

E t 0 KE t 0 K 2 E

 002 (t -s) b s, Θ ε,v s -b(s, X s ) ϑ ε h ε ds p ≤ (t -s)C b |η ε,v s | ds p |η ε,v s | p ds.The last two terms of (B.1) are also uniformly bounded in n and ε because, similarly to (A.2), Hölder's inequality and the linear growth condition H3a. If the latter fails we rely on H3b, Assumption 3.9 and the same calculations as in (A.4) to obtain a similar bound.

u 2 du p/ 2 + 8 2 du 2 du 2 du p/ 2 .

 2282222 p-1 N p/2 E t s K(t -u)σ ε u, Θ ε,v ε u -u) -K(s -u) σ ε u, Θ ε,v ε u -u)σ ε u, Θ ε,v ε u

Remark 3.6. Item

  The functional I defined by(3.2) has compact level sets. (iii) Assumption 3.3 holds. Then the family {G ε (W )} ε>0 satisfies the Laplace principle and, by equivalence, the Large Deviations Principle with rate function I and speed ε -2 . (i) entails that, for all N > 0, all v ∈ A N and all families {v ε } ε>0 in A N converging in distribution to v, there exists a subsequence such that lim n G εn W + ε -1

	definition belongs to G 0 v,N .	n	• 0 v εn s ds exists and by
	Remark 3.7. In the large deviations literature, a rate function is sometimes called 'good' if it has
	compact level sets. All the rate functions in the present paper satisfy this requirement (item (ii) above
	takes care of that) therefore we drop the adjective 'good'.

  Remark 3.19. H5 entails that(3.13) has a unique solution and yields the bound ∇b(•, X) T < ∞ by continuity. H7 implies H1. We have already proved in Lemma 3.8 that the moments of all orders of Θ ε,v are bounded hence (3.17) is automatically satisfied if Ξ is of polynomial growth. This is however not sufficient for the applications we have in mind where Ξ is of exponential growth.

The main theorem of this section is the following. Theorem 3.20 (Moderate Deviations). Under H2 -H8, Assumptions 2.3 and 3.9, the family {η ε } ε>0 satisfies a Large Deviations Principle (equivalently {X ε } ε>0 satisfies a Moderate Deviations Principle) with speed h 2 ε and rate function

TABLE 1 .

 1 Summary of rough volatility results and form of the rate functions (CF=closed-form; IF=integral form; OP=optimisation problem; shadowed cells are new contributions from this paper). σ corresponds to the implied volatility, defined precisely in Section 4.1.3.

						Rough	multi-factor
		Rough Stein-Stein	Bergomi	rough Bergomi	Rough Heston
		Small-time	Tail	Small-time	Small-time	Small-time [42]	Tail
		LDP [56]	MDP	LDP [56]	MDP	LDP [58]	MDP	LDP [61]	MDP	LDP MDP LDP
	(X ε , Y ε )	IF	IF	IF	-	IF	IF	IF	IF	IF	IF	IF
	X ε	OP	IF	OP	-	OP	IF	OP	OP	OP	IF	OP
	X ε	OP	CF	OP	-	OP	CF	OP	OP	OP	CF	OP
	σ	OP	CF	OP	-	OP	CF	OP	OP	OP	CF	OP
	Y ε	IF	IF	IF	IF	IF	IF	IF	IF	IF	IF	IF
	Y ε	CF	CF	IF	IF	CF	CF	CF	CF	IF	CF	IF

  4.3.1. Rough Stein-Stein. This model was defined in Section 4.2.1, but with the rescaling Y ε t := εY t and X ε t

  3 holds and hence Theorem 3.11 applies, and the form of the rate function in (L1) follows by inverting the relationships between (u, v) and (φ, ϕ). (L2) and (L3) follow from the same steps as in Corollary 4.15.4.3.3. Implied volatility asymptotics.We can also obtain implied volatility asymptotics since, by the same arguments as before, exp(X ε ) is a martingale in both the rough Stein-Stein and rough Heston models. Mapping ε 2 to 1/k we have from Corollaries 4.15 and 4.17 respectively that, for each t ∈ T,

	Corollary 4.18. In both the rough Stein-Stein and rough Heston models, for each t ∈ T, the implied volatility σ satisfies lim k↑∞ σ(t, k) 2 t k = -1 1 2 inf y≥1 I X t (y) ,
	where I X t is the respective rate function, given in Corollaries 4.15(L3) and 4.17(L3).
	Proof.

  A.1. Let f : T → R + and K a kernel satisfying Assumption 2.3. If there exists c 1 , c 2 ≥ 0 such that then f is uniformly bounded on T by a constant depending only on c 1 , c 2 , K 2 , T . If c 1 = 0 then f = 0.

	f (t) ≤ c 1 + c 2	0	t	|K(t -s)| 2 f (s) ds,	for all t ∈ T,
	Proof of Lemma A.1. By definition |K(t -s)| is a Volterra kernel in the sense of [53, Definition 9.2.1]. Following similar arguments as in the proof 2 = d i,j=1 |K ij (t -s)| 2 ½ s≤t 2 , and K(t, s) := c 2 |K(t -s)|
	of [1, Lemma 3.1], the generalised Gronwall lemma [53, Theorem 9.8.2] yields the bound

MDP moment bounds: Proof of Lemma 3.21.

  which concludes the proof of the closure and therefore of the compactness of L N . Let p ≥ 2, N > 0, v ∈ A N , ε > 0 and t ∈ T. Starting from (3.15), we use Cauchy-Schwarz and BDG inequalities to obtain

	APPENDIX B. TECHNICAL MODERATE DEVIATIONS PROOFS
	B.1.

  Overall this results in the existence of a constant C 2 > 0 independent of ε, v, t such thatE |η ε,v t | p ≤ C 2 + C 2 ]ds, for all t ∈ T, and Lemma A.1 yields the bound uniform in ε and t. B.2. MDP tightness: Proof of Lemma 3.22. Let p > 2 ∨ 2/γ, N > 0, a sequence {v ε } ε>0 in A N , ε > 0, and 0 ≤ s < t ≤ T . We proceed as in Lemma 3.15; starting from (3.15), applying consecutively H8, then (B.2), Cauchy-Schwarz and BDG inequalities we obtain

	0	t	|K(t -s)|	2 E [|η ε,v s |

p

  Hence Assumption 2.3 yields the existence of a constant C 2 > 0 such that

					2 du	p/2	+	s	t	K(t -u) 2 du	p/2
	= C 1	0	s	K(u + t -s) -K(u)	2 du	p/2	+	0	t-s	K(u) 2 du	p/2	.

DEPARTMENT OF MATHEMATICS, IMPERIAL COLLEGE LONDON, AND ALAN TURING INSTITUTE

Email address: a.jacquier@imperial.ac.uk

DEPARTMENT OF MATHEMATICS, IMPERIAL COLLEGE LONDON

Email address: a.pannier17@imperial.ac.uk