Time-inconsistency with rough volatility
Bingyan Han Hoi Ying Wong
q-fin.MF math.OC 91G80, 91A80, 60G22, 60H20






Time-Inconsistency with Rough Volatility[An earlier version of this paper was circulated and cited under the title “Time-consistent feedback strategies with Volterra processes". The authors would like to thank two anonymous referees and the editors for their careful reading and valuable comments, which have greatly improved the manuscript. Bingyan Han is supported by UIC Start-up Research Fund (Reference No: R72021109).]
    Bingyan Han Division of Science and Technology, BNU-HKBU United International College, Zhuhai, Guangdong, China, bingyanhan@uic.edu.cn
	Hoi Ying WongDepartment of Statistics, The Chinese University of Hong Kong, Hong Kong SAR, China, hywong@cuhk.edu.hk

    December 21, 2021
========================================================================================================================================================================================================================================================================================================================================================================================================================================



Time-Inconsistency with Rough Volatility[An earlier version of this paper was circulated and cited under the title “Time-consistent feedback strategies with Volterra processes". The authors would like to thank two anonymous referees and the editors for their careful reading and valuable comments, which have greatly improved the manuscript. Bingyan Han is supported by UIC Start-up Research Fund (Reference No: R72021109).]
    Bingyan Han Division of Science and Technology, BNU-HKBU United International College, Zhuhai, Guangdong, China, bingyanhan@uic.edu.cn
	Hoi Ying WongDepartment of Statistics, The Chinese University of Hong Kong, Hong Kong SAR, China, hywong@cuhk.edu.hk

    December 21, 2021
========================================================================================================================================================================================================================================================================================================================================================================================================================================




	 In this paper, we consider equilibrium strategies under Volterra processes and time-inconsistent preferences embracing mean-variance portfolio selection (MVP). Using a functional Itô calculus approach, we overcome the non-Markovian and non-semimartingale difficulty in Volterra processes. The equilibrium strategy is then characterized by an extended path-dependent Hamilton-Jacobi-Bellman equation system under a game-theoretic framework. A verification theorem is provided. We derive explicit solutions to three problems, including MVP with constant risk aversion, MVP for log returns, and a mean-variance objective with a linear controlled Volterra process. We also thoroughly examine the effect of volatility roughness on equilibrium strategies. Numerical experiments demonstrate that trading rules with rough volatility outperform the classic counterparts.
	
 
	Keywords: Time-inconsistency, rough volatility, functional Itô calculus, mean-variance portfolio selection, Volterra Heston model.
	

	Mathematics Subject Classification: 91G80, 91A80, 60G22, 60H20. 

	JEL Classification: C72, C73, G11. 




§ INTRODUCTION

Portfolio selection is a fundamental and leading problem in mathematical finance. Pioneered by <cit.>, mean-variance portfolio selection (MVP) is well recognized as a cornerstone of modern portfolio theory. Its intuitive and flexible formulation has attracted the attention of numerous researchers, who have sought to strengthen the original framework. Examples include, but are not limited to, <cit.>. In contrast to expected utility theory <cit.>, MVP suffers from time-inconsistency induced by the variance operator. When an investor with an initial value of (t, x) perceives that the derived strategy is no longer optimal at a later state (s, X_s) for s > t, time-inconsistency occurs. As noted in <cit.>, time-consistency is a fundamental requirement for any reasonable strategic decision-making and optimization. An agent should only choose strategies from which he/she will not deviate. <cit.> identifies that MVP has an adjustment term that provides “an incentive for the investor to deviate from his optimal strategy at a later time."  However, studies such as <cit.> have neglected time-inconsistency and have only provided the pre-committed strategy. 


Time-inconsistency has generated an astonishing amount of controversial opinions toward the notion of optimality. A remedy is to consider an equilibrium strategy. The intent is to formulate a game between the current agent and his/her future selves and to then derive an equilibrium of the game as the strategy. Several treatments are available based on different methodologies. <cit.> follow the classic dynamic programming framework and derive an extended Hamilton-Jacobi-Bellman (HJB) equation to characterize the equilibrium. <cit.> formulates the game in an open-loop optimization and derives the equilibrium via the stochastic maximum principle. <cit.> consider a partition on the whole planning time horizon and obtain the equilibrium by taking limits. With a fixed-point argument, <cit.> further distinguishes between strong and weak equilibria. Related discussions include <cit.> and the references therein. In this paper, we exploit the extended HJB methodology for its wider applications, including MVP. 

<cit.> document that hedging demand from stochastic volatility can comprise a substantial percentage of total equilibrium stock exposure. Realistic modeling for volatility then plays an indispensable role in investment decision-making. We use rough volatility models, recently proposed by <cit.>. In this seminal work, <cit.> uses the fractional Brownian motion (fBm) to demonstrate the roughness of volatility. By developing a rigorous statistical estimation and inference, <cit.> further confirms that volatility is even rougher than reported in <cit.>. These models are consistent with some stylized facts of financial time series and have several desired theoretical properties. They capture the term structure of the implied volatility (IV) surface, especially for the explosion of the at-the-money (ATM) skew when the maturity is close to zero <cit.>, which smooth volatility models fail to do. Examples of rough volatility models include the fBm <cit.>, the fractional Ornstein-Uhlenbeck (fOU) process <cit.>, the rough Bergomi (rBergomi) model <cit.>, the fractional Heston model <cit.>, and the rough Heston model <cit.>. The rough Heston model has received particular attention and has been extended to the Volterra Heston model <cit.> and the affine forward variance (AFV) model <cit.>. The economic interpretation of rough volatility is explained by <cit.> via metaorders in high-frequency trading, by <cit.> via the no-arbitrage property, by <cit.> via heterogeneity in near-term downside risk, and by <cit.> via long memory behavior. However, like the debate on the short-range or long-range dependence in volatility, the understanding of rough volatility is still in development.

We are interested in the equilibrium strategies under a more realistic stochastic financial environment delineated by rough volatility. We are the first to explore time-inconsistency with rough volatility, although related works, such as <cit.>, and the references therein are available for alternative portfolio problems under rough volatility. Despite empirical evidence of rough volatility, its non-Markovian and non-semimartingale nature challenges the classic methodology of equilibria. We also consider Volterra processes for generality. Previous results in the literature cannot be directly adopted to tackle equilibrium strategies under Volterra processes. <cit.> accounts for the linear non-Markovian systems, but the application is limited to linear-quadratic control problems. The locally mean-variance efficient (LMVE) approach in <cit.> can deal with semimartingales but the Volterra process is not typically a semimartingale.

To tackle these difficulties, we adopt a general methodology called functional Itô calculus with applications far beyond time-inconsistency and rough volatility. <cit.> first develops a pathwise calculus for non-anticipative functionals, motivated by pricing and hedging path-dependent derivatives. By defining the time and spatial derivatives, the classic Itô formula is extended to the functional Itô formula for path-dependent functionals in <cit.>. The functional Itô calculus is useful for a wide class of optimization and decision-making problems. To better incorporate financial and insurance risks, <cit.> considers the applications of the functional Itô calculus in convex risk measures with non-Markovian jump-diffusion processes. <cit.> uses Dupire's functional Itô calculus to motivate their definition of contracts in principal-agent problems. Under the framework of <cit.>, <cit.> presents two pathwise versions of the master formula in Fernholz's stochastic portfolio theory and elucidates their performance with empirical data. <cit.> applies the functional Itô calculus to stochastic differential games and optimal control problems with delay. 

However, the aforementioned literature relies on a semimartingale assumption, whereas the Volterra Heston model and general Volterra processes are non-semimartingales. Recently, <cit.> has developed a powerful toolkit for the functional Itô formula to analyze functionals of Volterra processes. Heuristically, their approach aims to “recover" the flow property of Volterra processes by incorporating an auxiliary non-anticipative process Θ^t (<ref>) into the path ω. Their elegant results enable us to derive the extended path-dependent HJB (PHJB) equation system in Theorem <ref>, which extends the results of <cit.> and has potential applications in addition to portfolio selection under rough volatility. Nevertheless, we stress that the development of the PHJB system is non-trivial even given the existing results. We also offer an example of the unsolved future problem suggested in <cit.>:
bjork2017time
	“The present theory depends critically on the Markovian structure. It would be interesting to see what can be done without this assumption.”


In Section <ref>, we apply the general framework for Volterra processes to time-consistent (TC) MVP under rough volatility. We refer to the agent under a rough stochastic environment as a rough investor. Should a rough investor buy more or less when volatility is rougher? When should he/she change his/her preference to rough? Furthermore, how profitable are rough strategies according to empirical data? We observe that volatility roughness dramatically increases the demand for hedging. Additionally, rough investors' attitude toward roughness depends on the payoff functional in mind. By deriving explicit solutions to classic problems in <cit.>, we present a thorough analysis and disentangle the complicated relationship between time-inconsistency and roughness. Our findings advocate rough volatility models as a promising alternative for the classic models adopted in <cit.>.

	
  * Section <ref> considers TC-MVP with constant risk aversion <cit.> under the Volterra Heston model. It is referred to as the const-MV case. In our sensitivity analysis, we find that when the investment horizon is long, the const-MV strategy suggests high demand on stocks with smoother volatility. For a short investment horizon, the const-MV strategy increases exposure to stocks with rougher volatility. We refer to this phenomenon as the investment horizon effect. In the const-MV case, the change point in this effect is irrelevant to the investor's risk aversion. In a simulation study, we find that a rough investor demands up to 40% more than a smooth investor.
	
  * Section <ref> investigates TC-MVP for log returns <cit.> under a rough stochastic environment. It is referred to as the log-MV case. The investment horizon effect in the const-MV case is retained in the log-MV case. However, the main difference is that heterogeneity in risk aversion changes when investors start to prefer rough. A more risk-averse investor prefers rough much earlier. In general, the log-MV case prohibits bankruptcy and is thus more conservative than the const-MV case. Rough investors increase their stock demand by at most 9% compared with <cit.>, as detailed in the simulation study. In addition, the log-MV criterion implies different roughness-related behaviors, compared with constant relative risk aversion (CRRA) utility.
	
  * Section <ref> studies a problem with the mean-variance (MV) objective and a linear controlled Volterra process. Distinct from portfolio selection, the convolution and the control appear together in the state process (<ref>). We manage to derive an explicit strategy, which demonstrates the potential application of our framework to problems with controlled Volterra state processes.
	
	
  * Via the empirical study in Section <ref>, we highlight that trading strategies with rough volatility dominate all classic counterparts in <cit.>. They gain greater terminal wealth with a better Sharpe ratio, even under the volatile U.S. equity markets from May 2018 to April 2019. The empirical performance substantiates the claim that our proposed strategies are more effective in capturing market movements.



The reminder of this paper is organized as follows. Section <ref> describes the general framework with Volterra processes. Section <ref> reviews the Volterra Heston model as a main example of rough volatility. Section <ref> derives the extended PHJB equation system. Section <ref> discusses the solutions to three problems using the MV objective. Section <ref> presents the numerical study. Finally, Section <ref> concludes the paper. The functional Itô calculus in <cit.> is summarized in Appendix <ref> for a self-contained article. All mathematical proofs are deferred to Appendix <ref>.



§ PROBLEM FORMULATION



 §.§ Volterra processes

Let T>0 be a deterministic finite investment horizon. We first present a general model with applications beyond rough volatility. Consider a controlled n-dimensional stochastic Volterra integral equation (SVIE) on [0, T]:

    X^_t = x + ∫^t_0 μ(t; r, X^_, (r, X^_)) dr + ∫^t_0 σ(t; r, X^_, (r, X^_)) dW_r,

where X^_ refers to the whole past path of the process (X^_s)_0 ≤ s ≤ r, W is a d-dimensional standard Brownian motion, and μ, σ are adapted with suitable dimensions. The feedback strategy  is a k-dimensional deterministic measurable function. We provide a rigorous definition of admissible strategies in Definition <ref>. The main example of (<ref>) in this paper is the two-dimensional process (M^, ν) in (<ref>) and (<ref>).  It is also worth mentioning again that SVIE (<ref>) is non-Markovian and non-semimartingale in general.

We consider feedback strategies (r, X^_) that depend on the whole path X^_ instead of solely on the current value X^_r of the process as in <cit.>. This setting is more reasonable because investors can always base their decisions on the observed history of the process. 

Before formally defining the equilibrium feedback strategies, we impose the following standing assumption throughout this paper.

	Controlled SVIE (<ref>) admits a unique in law continuous weak solution (X^, W) on some complete probability space (Ω, , ) with a filtration = {_t }_ 0 ≤ t ≤ T that satisfies the usual conditions. (X^, W) is -adapted and
	
    [ sup_ 0 ≤ t ≤ T |X^_t |^p ] < ∞,

	for any p ≥ 1.


As noted in the Appendix of <cit.>, a sufficiently large moment constant p is enough. However, this p is not explicit. As our primary focus is time-inconsistency, we do not pursue potentially more general conditions validating Assumption <ref> in this paper. We verify Assumption <ref> for some examples in Section <ref> and refer interested readers to <cit.> for further results. Indeed, as  is a feedback strategy, we can regard μ^ and σ^ in (<ref>) as the drift and the diffusion, respectively. The results without controls in <cit.> then become applicable to our cases. Compared with Assumption 3.1 of <cit.>, Assumption <ref> further requires (<ref>) to admit a unique in law solution. For a given feedback strategy , it is natural for our problem to attain a unique reward functional (<ref>) under , which requires the law of SVIE (<ref>) to be unique. We also need a continuous solution to (<ref>). This condition is relatively mild when a feedback strategy is considered. The concatenated path (<ref>) is later justified as continuous under this condition. We fix a weak solution (X^, W) to (<ref>) once the feedback strategy  is given. Although the probability space and Brownian motion are also parts of the solution to SVIE (<ref>), the distribution of X^ is unique. We denote [·] and [·|_t] as the expectation and conditional expectation, respectively, under the probability measure , which is a part of the weak solution under a generic control . As in <cit.>, we use = ^X^, W.



 §.§ Example: Volterra Heston model

Consider a two-dimensional standard Brownian motion W ≜ (W_1, W_2). An important example of (<ref>) is a rough version of the classic Heston model <cit.>:

    ν_t= ν_0     + 1/Γ(H + 1/2)∫_0^t (t - r)^H - 1/2κ (ϕ - ν_r ) dr  + 1/Γ(H + 1/2)∫_0^t (t - r)^H-1/2σ√(ν_r) d B_r,

where Γ(·) is the Gamma function and H is the Hurst parameter. dB_r = ρ dW_1r + √(1 - ρ^2) dW_2r and ν_0, κ, ϕ, and σ are positive constants. The correlation ρ between stock price and volatility is also constant. When H = 1/2, the model is reduced to the classic Heston model adopted in <cit.>. The volatility trajectories of (<ref>) have almost surely Hölder regularity H - ε, for all ε > 0, as shown in <cit.>. Therefore, (<ref>) is called the rough Heston model and the Hurst parameter H is an index of volatility roughness. The smaller H is, the rougher the volatility is. With H of order 0.1, <cit.> shows that the rough Heston model provides remarkable fits for volatility skews, including cases of extreme short maturity. Thus, it captures the near-term risks implied by ATM skew explosion.

Extending the rough Heston model (<ref>), the Volterra Heston model in <cit.> reads as follows:

    ν_t= ν_0 + κ∫_0^t K(t - r)(ϕ - ν_r ) dr + ∫_0^t K(t - r) σ√(ν_r) d B_r,

where K(·) is the kernel function. By setting K(t) = t^H-1/2/Γ(H+1/2), namely the fractional kernel, (<ref>) recovers (<ref>). In line with <cit.>, we impose the following assumption on the kernel function.


	The kernel K is strictly positive and completely monotone. There exists τ∈(0,2] such that ∫_0^h K(t)^2 d t=O(h^τ) and ∫_0^T(K(t+h)-K(t))^2 d t=O(h^τ) for every T<∞. 


Like <cit.>, the risky asset (stock) price S_t is postulated as follows:

    dS_t = S_t (_t + θν_t) dt + S_t √(ν_t) dW_1t,    S_0 > 0,

with a deterministic bounded risk-free rate _t>0 and constant θ≠ 0. The market price of risk, or risk premium, is then given by θ√(ν_t). Such a risk premium specification is widely used in the literature, such as in <cit.> and <cit.>. The risk-free rate _t >0 is the return of a risk-free asset available in the market. Indeed, a general Heston specification <cit.> is also tractable, as indicated in Remark <ref>. We adopt (<ref>) to simplify the presentation.

We quote the following result from <cit.>, which guarantees the existence and weak uniqueness of the Volterra Heston model.

	Under Assumption <ref>, the stochastic Volterra equation (<ref>)-(<ref>) has a unique in law _≥ 0×_≥ 0-valued continuous weak solution for any initial condition (S_0, ν_0) ∈_≥ 0×_≥ 0.



	Pathwise uniqueness for (<ref>)-(<ref>) remains an open problem in general. We mention <cit.> as a related result with kernel K ∈ C^1([0, T], ) and <cit.> for certain smooth kernels. However, the strong uniqueness of (<ref>)-(<ref>) is left open for singular kernels. For weak solutions, Brownian motion is free to construct as needed. In the sequel, we fix a solution (S, ν, W_1, W_2) to (<ref>)-(<ref>) as other solutions share the same law. Furthermore, the boundary point 0 may be reachable for the Volterra Heston model and the property of the boundary point 0 is left open.


Multiply the dollar amount of wealth in the stock by √(ν_t) and denote it as the investment strategy . Let M^_t be the wealth process. Then, M^_t satisfies

    d M^_t = (_t M^_t + θ√(ν_t) u_t ) dt + u_t dW_1t,    M^_0 > 0.

It is clear that (M^, ν) in (<ref>) and (<ref>) is a special case of the Volterra process (<ref>). We handle the general time-inconsistent problems in (<ref>) in a unified way, whereas the applications focus on (<ref>) and (<ref>).



 §.§ Equilibrium under time-inconsistent preferences


For time-inconsistent problems, we must consider the state process starting from time t ∈ [0, T). For s ≥ t, the general state process (<ref>) can be decomposed as follows:

    X^_s = x    + ∫^t_0 μ(s; r, X^_, (r, X^_)) dr + ∫^t_0 σ(s; r, X^_, (r, X^_)) dW_r 
       + ∫^s_t μ(s; r, X^_, (r, X^_)) dr + ∫^s_t σ(s; r, X^_, (r, X^_)) dW_r.

Following <cit.>, we define

    Θ^t, _s ≜ x     + ∫^t_0 μ(s; r, X^_, (r, X^_)) dr + ∫^t_0 σ(s; r, X^_, (r, X^_)) dW_r,    t ≤ s ≤ T.

t ↦Θ^t, _s is a semimartingale for 0 ≤ t ≤ s. Using Θ^t, _s, we concatenate a path ω as

    ω_s =  (X^⊗_t Θ^t, )_s ≜ X^_s 𝕀_{0 ≤ s < t} + Θ^t, _s 𝕀_{t ≤ s ≤ T}.

Although ω is defined on [0, T], it is adapted to _t. ω is - continuous.

An interpretation of Θ^t, _s is that it can be written as follows:

    Θ^t, _s = [ X^_s - ∫^s_t μ(s; r, X^_, (r, X^_)) dr | _t ],

which is related to the modified forward processes <cit.>. It represents the current view of process distributions in the future.

Particularly, if we consider the Volterra Heston model (<ref>), then

    Θ^t_s = ν_0 + κ∫_0^t K(s - r)(ϕ - ν_r ) d r + ∫_0^t K(s - r) σ√(ν_r) d B_r,

which corresponds to the variance part of Θ^t, in (<ref>). As  does not appear in the variance process, we exclude it from the notation, which becomes Θ^t. We further denote the concatenated path ω for the variance process as follows: 

    ω^ν_s =  (ν⊗_t Θ^t)_s ≜ν_s 𝕀_{0 ≤ s < t} + Θ^t_s 𝕀_{t ≤ s ≤ T}.

For Θ^t_s in (<ref>), <cit.> shows that Θ^t_s can be represented by the forward variance curve [ν_s|_t]. Therefore, Θ^t_s can be roughly replicated with financial products, such as variance swaps.  

At time t, for a realized path ω, we have the following:

    X^t, ω, _s = ω_s     + ∫^s_t μ(s; r, X^t, ω, _, (r, X^t, ω, _)) dr 
        + ∫^s_t σ(s; r, X^t, ω, _, (r, X^t, ω, _)) dW_r,   t ≤ s ≤ T, 
    
    X^t, ω, _s = ω_s,       0 ≤ s < t,

where the notation X^_s is replaced with X^t, ω, _s to highlight its dependence on t and the path ω. For t ≤ s ≤ T, Θ^t, _s is then interpreted as follows:

    Θ^t, _s = ω_s = x + ∫^t_0 μ(s; r, X^t, ω, _, (r, X^t, ω, _)) dr + ∫^t_0 σ(s; r, X^t, ω, _, (r, X^t, ω, _)) dW_r.
 

For a given feedback strategy , let

    μ^(t; r, ω) ≜μ(t; r, ω_, (r, ω_)),   σ^(t; r, ω) ≜σ(t; r, ω_, (r, ω_)).

As it is enough for us to consider μ^ and σ^ with the same singularities, we encounter two cases only. If lim_r→ tμ^(t; r, ·) = ∞ and lim_r→ tσ^(t; r, ·) = ∞, it is called a singular case; otherwise, if lim_r→ tμ^(t; r, ·) < ∞ and lim_r→ tσ^(t; r, ·) < ∞, it is called a regular case <cit.>. 

We introduce the reward functional as follows:

    J(t, ω; ) ≜   [ ∫^T_t C(t, ω_t, r,  X^t, ω, _, (r, X^t, ω, _)) dr + F(t, ω_t, X^t, ω, _T ∧·) | _t ]     + G(t, ω_t, [ X^t, ω, _T | _t]),

where X^t, ω, is given by (<ref>).

Functional (<ref>) has nested MV criterion (<ref>) as a special case. SVIE (<ref>) is not time-consistent due to the absence of the flow property <cit.>. However, we focus on the time-inconsistency issue from the objective function J, which originates from its dependence on the current time t, the current state ω_t, and the nonlinear function G. The reward functional (<ref>) does not satisfy the Bellman optimality principle. A strategy that maximizes (<ref>) at the current time t may no longer be optimal at a certain future time. We refer readers to <cit.> for motivations for and examples of (<ref>). 

Conceptually, the non-Markov property implies that it is not enough to record the current state X^_t only. More information from _t is needed. For Volterra processes, the concatenated path ω is sufficient. Furthermore, by writing the reward J as a functional of ω = X^⊗_t Θ^t, rather than X^_t∧· only, J preserves some nice regularity properties, such as continuity, under mild conditions. To clarify it further, although Θ^t, is a functional of X^_t∧·, the dependence is usually discontinuous under uniform convergence due to the stochastic integrals involved. Readers may refer to <cit.> for a specific example. <cit.> discovers that the flow property can be recovered by including Θ^t,, resulting in the functional Itô formula. Section <ref> uses specific examples to clarify the rationale more concretely. 

Let m be a generic positive value for the polynomial growth rate, which may vary from line to line. By the supremum norm ||ω||_T ≜sup_ 0 ≤ t ≤ T |ω_t| defined in Appendix <ref>, we introduce continuity in ω under ||·||_T.

	Properties for F and G:
	
		
  * For any fixed s and y, F is of polynomial growth in ω. That is,
		
    |F(s, y, ω)| ≤ C_0[ 1 + ||ω||^m_T],

		for some constants C_0, m > 0. 
		
		
  * For any fixed s and y, G(s, y, z) is continuously differentiable in z.
	

Similarly, for a given feedback strategy , let

    C^(t, ω, s, y)  ≜ C(s, y, t, ω_t ∧·,  (t, ω_t ∧·)).



	 is said to be an admissible strategy, denoted by ∈, if 
	
		
  * Assumption <ref> holds.
		
  * (a) If μ^ and σ^ are regular, then for a fixed t ∈ [0, T], assume that μ^(t; r, ω) and σ^(t; r, ω) are right-continuous in r ∈ [0, t] and continuous in ω∈Ω.  ∂_t μ^(t; r, ·) and ∂_t σ^(t; r, ·) exist for t ∈ [r, T]. For φ = μ^, σ^, ∂_t μ^, ∂_t σ^, 
		
    |φ(t; r, ω)| ≤ C_0[ 1 + ||ω||^m_T],

		for some constants C_0, m > 0.
		
		(b) If μ^ and σ^ are singular, then for a fixed t ∈ [0, T], assume that μ^(t; r, ω) and σ^(t; r, ω) are right-continuous in r ∈ [0, t) and continuous in ω∈Ω.  For φ = μ^, σ^, suppose that ∂_t φ (t; r, ·) exists for t ∈ (r, T], and there exists 0 < H < 1/2 such that for any 0 ≤ r < t ≤ T,
		
    |φ(t; r, ω)|    ≤ C_0[ 1 + ||ω||^m_T](t - r)^H-1/2, 
    
    		|∂_t φ(t; r, ω)|    ≤ C_0[ 1 + ||ω||^m_T](t - r)^H - 3/2,

		for some constants C_0, m > 0.
		
		
  * For any fixed s and y, C^ is continuous in (t, ω). C^ is of polynomial growth in ω, uniformly in t. That is,
		
    |C^(t, ω, s, y)| ≤ C_0[ 1 + ||ω||^m_T],

		for some constants C_0, m > 0.
		
  * For any fixed s and y, 
		
    [    sup_ t ≤ r ≤ T| C(s, y, r,  X^t, ω, _, (r, X^t, ω, _)) | + | F(s, y, X^t, ω, _T ∧·) | | _t ]  
        + | G(s, y, [ X^t, ω, _T | _t]) | ≤ C_0[ 1 + ||ω||^m_T],

		for some constants C_0, m > 0, which are independent of t.
	


As a strategy that is optimal at time t may no longer be optimal afterward, an agent is motivated to deviate from it. <cit.> argue that any reasonable agent should only choose strategies from which she will not deviate. Informally, the agent who is aware of time-inconsistency can consider her selves at different future times as different agents. The agent at time t controls the state X^ exactly at time t by choosing a control function (t, ·). A candidate admissible equilibrium strategy  should have the following property: if for each r > t, the agent at time r chooses (r, ·), then it is optimal for the agent at time t to choose (t, ·). Such an equilibrium strategy is formulated backwardly by induction. Therefore, the derived  will not be deviated from.

The informal arguments above work well in discrete time; but in continuous time, controlling on a time set of Lebesgue measure zero has no effect on the state process. Thus, the agent at time t is allowed to act on an infinitesimally small interval [t, t+h) and then send h to zero. Formally, let (r, ω_) be a deterministic map that is also admissible. Perturbing  in the same way as <cit.> yields the following:

    _h(r, ω_) = {[ (r, ω_),    t ≤ r < t+h,; (r, ω_),    t+h ≤ r ≤ T. ].

If we denote the solution to SVIE (<ref>) with _h as X^_h, the feedback strategy reads as follows:

    _h(r, X^_h_) = {[ (r, X^_h_),    t ≤ r < t+h,; (r, X^_h_),    t+h ≤ r ≤ T. ].

A crucial characteristic of the feedback (closed-loop) formulation is that perturbing  on [t, t+h) implicitly affects the strategies on [t+h, T] through X^_h. It is different with open-loop strategies whose value on [t+h, T] is unchanged <cit.>.

Loosely speaking, the candidate equilibrium  should satisfy the following property: If all agents at [t+h, T] agree to use , then it is optimal for the agent at time t to adopt . Mathematically, we have the following definition.


	Consider a candidate equilibrium law . For any t ∈ [0, T) and ∈, where  is defined in Definition <ref>, define _h as in (<ref>),  is an (weak) equilibrium strategy if 
	
    lim inf_h ↓ 0J(t, ω; ) - J(t, ω; _h)/h≥ 0,

	for any ω∈Λ̃(, t).


In Definition <ref>, we consider a path-dependent counterpart of the concept of support. Λ̃(, t) is the support of paths for X^⊗_t Θ^t, conditional on _0. The support is the set of ω∈Ω such that any neighborhood of ω has a positive measure under the distribution of X^⊗_t Θ^t,. The metric is induced by norm ||·||_T. Roughly speaking, the support contains all possible situations for the paths. We refer to <cit.> for the rationale for considering the support rather than the whole space Ω.


	Similar to <cit.>, the definition of support differs from the standard definition in the literature. We refer readers to the footnote under <cit.> for further details. Characterizing the support Λ̃(, t) under SVIE (<ref>) remains an open problem. Related work of which we are aware includes <cit.>. However, in the examples we consider, the support is clear and relatively straightforward to obtain.


	As noted in <cit.>,  under (<ref>) may merely be a stationary point. Recent works have also considered equilibria under the following condition:
	
    J(t, ω; _h) ≤ J(t, ω; ),

	where _h is selected in certain sets. <cit.> clarifies three notions, namely strong, regular, and weak equilibria. <cit.> considers a stochastic control problem in which the generator of certain Markov chains can be controlled, with a definition such as that in (<ref>). However, weak equilibria should be considered first, as other types of equilibrium strategies are under stronger conditions that may be too restrictive.


To emphasize that probability is also part of the weak solution, denote the expectation and conditional expectation under equilibrium control  by [·] and [· | _t], respectively. [·] and [· | _t] are under a perturbed control _h in (<ref>). Therefore, in Definition <ref>, J(t, ω; ) is under [· | _t] and J(t, ω; _h) is under [· | _t].



 §.§ Extended path-dependent HJB equation

The following notations are useful. Define

    f^(t, ω, s, y)    ≜[ F(s, y, X^t, ω, _T ∧·) | _t], 
    
    g^(t, ω)    ≜[X^t, ω, _T | _t ], 
    
    c^r, (t, ω, s, y)    ≜[ C(s, y, r,  X^t, ω, _, (r,  X^t, ω, _)) | _t ].

When =, denote

    f(t, ω, s, y)    ≜[ F(s, y, X^t, ω, _T ∧·) | _t], 
    
    g(t, ω)    ≜[X^t, ω, _T | _t ], 
    
    c^r(t, ω, s, y)    ≜[ C(s, y, r,  X^t, ω, _, (r,  X^t, ω, _)) | _t ].

Our convention is that the last two arguments (s, y) are reserved for state-dependence. These auxiliary functions are reduced to their counterparts in <cit.> when there is no path dependence. 

First, we derive a recursive relationship, which extends <cit.> to the non-Markovian case applicable to our problem. To do so, we investigate the problem at time t+ h. Denote the path

    ω^t+h_s = {[         ω_s,   0 ≤ s < t,;  X^t, ω, _s, t ≤ s < t+h,;   Θ^t+h, _s, t+h ≤ s ≤ T, ].

where

    Θ^t+h, _s = x     + ∫^t+h_0 μ(s; r, X^t, ω, _, (r, X^t, ω, _)) dr + ∫^t+h_0 σ(s; r, X^t, ω, _, (r, X^t, ω, _)) dW_r.

Note that ω^t+h is adapted to _t+h but not _t. To make the notation compact, we write the following:

    ω^t+h_s = (X^t, ω, ⊗_t+hΘ^t + h, )_s ≜ X^t, ω, _s 𝕀_{0 ≤ s < t+h} + Θ^t + h, _s 𝕀_{t+h ≤ s ≤ T}.

Θ^t + h, is only defined on [t+h, T]. X^t, ω, _s ≠ω_s for  s ∈ (t, t+h) and Θ^t + h, _s ≠ω_s for s ∈ [t+h, T].


	For a general admissible feedback strategy , the reward functional J satisfies the following recursion:
	
    J(t, ω; ) =    [ J(t+h, ω^t+h; ) | _t ]  
        - {∫^T_t+h[ c^r, (t+h, ω^t+h, t+h, ω^t+h_t+h) | _t ] dr  
            - ∫^T_t [ c^r, (t+h, ω^t + h, t, ω_t)| _t ] dr }
        - {[ f^(t+h, ω^t+h, t+h, ω^t+h_t+h) | _t ] - [ f^(t+h, ω^t+h, t, ω_t) | _t ] }
         -{[ G(t+h, ω^t+h_t+h, g^(t+h, ω^t+h)) | _t ] - G(t, ω_t, [ g^(t+h, ω^t+h) | _t] ) }.


The proof of Lemma <ref> solely applies the tower property of conditional expectation and does not rely on the functional Itô formula. However, the verification theorem does need the functional Itô formula in <cit.>, quoted as Theorem <ref>. The derivatives and spaces C^1, 2_+( Λ) and C^1, 2_+, α(Λ) in Theorem <ref> are defined in <cit.>, which are also briefly reviewed in the Appendix <ref>.


	Suppose that (1)-(2) of Definition <ref> hold. Let f ∈ C^1, 2_+( Λ) for the regular case or f ∈ C^1, 2_+, α( Λ) for the singular case with β≜α + H - 1/2 > 0. The constant H is defined in Definition <ref> (2) for the singular case. Then,
	
    df(t, X^⊗_t Θ^t, ) =   ∂_t f(t, X^⊗_t Θ^t, ) dt + ∂_ω f(t, X^⊗_t Θ^t, ), μ^t,  dt     + 1/2∂^2_ωω f(t, X^⊗_t Θ^t, ), (σ^t, , σ^t, ) dt      + ∂_ω f(t, X^⊗_t Θ^t, ), σ^t,  d W_t,   -,

	where for φ = μ, σ, the notation φ^t, _s ≜φ^(s; t, ·) emphasizes the dependence on s ∈ [t, T].
	
	For the singular case, the derivatives related to ω are defined as follows:
	
    ∂_ω f(t, ω), φ^t,    ≜lim_δ↓ 0∂_ω f(t, ω), φ^δ, t, ,    φ = μ, σ, 
    ∂^2_ωω f(t, ω), (σ^t, , σ^t, )   ≜lim_δ↓ 0∂^2_ωω f(t, ω), (σ^δ, t, , σ^δ, t, ),

	where φ^δ, t, _s ≜φ^( s ∨ (t + δ); t, ·) for 0 < δ≤ T - t, is the truncated function. It also emphasizes the dependence on s ∈ [t, T].


Define the value function as follows:

    V(t, ω) = J(t, ω; ).

For a general admissible control  and a functional f(t, ω) that satisfies Assumption <ref>, as specified later, denote the operator ^ as follows:

    (^ f)(t, ω) ≜∂_t f(t, ω) + ∂_ω f(t, ω), μ^t,  + 1/2∂^2_ωω f(t, ω), (σ^t, , σ^t, ),

where we omit the arguments in μ^t, and σ^t,. The derivatives in (<ref>) are defined in (<ref>), (<ref>), and (<ref>) for regular cases, while (<ref>) and (<ref>) are for singular cases. The operator ^ only applies to variables within parentheses. For instance, (^ f)(t, ω, t, ω_t) operates on t, ω, t, ω_t, whereas (^ f^s, y)(t, ω) operates on t, ω only.


	The extended PHJB equation system is defined as follows:
	
		
  * The function V satisfies
		
    sup_∈{    (^ V)(t, ω) + C(t, ω_t, t, ω_t ∧·,  (t, ω_t ∧·)) - ∫^T_t (^ c^r)(t, ω, t, ω_t) dr     + ∫^T_t (^ c^t, ω_t, r)(t, ω) dr - (^ f)(t, ω, t, ω_t) + (^ f^t, ω_t)(t, ω)     - ^(G ♢ g)(t, ω) + ∂_y G(t, ω_t, g(t, ω)) (^ g)(t, ω) } = 0,    0 ≤ t ≤ T, 
        V(T, ω) = F(T, ω_T, ω) + G(T, ω_T, ω_T).
 
		Let  be the strategy that attains the supremum.
		
  * For each fixed s and y, f^s,y(t, ω) is defined as follows: 
		
    (^ f^s, y)(t ,ω) = 0,    0 ≤ t ≤ T,    f^s, y(T, ω) = F(s, y, ω).

		
  * The function g satisfies 
		
    (^ g)(t, ω) = 0,    0 ≤ t ≤ T,    g(T, ω) = ω_T.

		
  * For each fixed s, r, and y, c^s, y, r is defined by
		
    (^ c^s, y, r) (t, ω)     = 0,    0 ≤ t ≤ r, 
    
    		c^s, y, r(r, ω)     = C(s, y, r, ω_, (r, ω_)).

		
  * The notations have the following meanings:
		
    f(t, ω, s, y)    = f^s,y(t, ω),      c^r(t, ω, s, y) = c^s, y, r(t, ω), 
    
    		(G ♢ g)(t, ω)    = G(t, ω_t, g(t, ω)),     ∂_y G(t, ω_t, y) = ∂ G/∂ y (t, ω_t, y).

		
  * The probabilistic interpretations are as follows:
		
    f^s, y(t, ω)     = [ F(s, y, X^t, ω, _T ∧·) | _t ],      g(t, ω) = [X^t, ω, _T | _t ], 
    
    		c^s, y, r(t, ω)    = [ C(s, y, r,  X^t, ω, _, (r,  X^t, ω, _)) | _t ],    0 ≤ t ≤ r.

	
	Equations (<ref>)-(<ref>) above hold for ω∈Λ̃(, t), t ∈ [0, T].


	The spatial region in (<ref>)-(<ref>) is expressed as Λ̃(, t),   t ∈ [0, T], which is consistent with Definition <ref>. For example, for MVP with state-dependent risk aversion in <cit.>, the assumption is that wealth stays positive implicitly. This implies that the system in <cit.> holds for region M > 0 instead of M ∈.


We must emphasize the dependence on ω and ω_t ∧· in (<ref>)-(<ref>). Although the functionals V, c^r, c^s, y, r, f, f^s, y, and g generally depend on the whole path ω, the strategy  only depends on ω_t ∧·, paths up to time t. C(t, ω_t, t, ω_t ∧·,  (t, ω_t ∧·)) also only depends on ω_t ∧·. This follows from the definition of paths in (<ref>) and the fact that  and C only depend on X^ but not Θ^t, directly. If there is no path dependence, the system in (<ref>)-(<ref>) reduces to the one in <cit.>.

We impose the regularity condition, Assumption <ref>, on the functionals appearing in the extended PHJB system in Definition <ref>. This condition validates that all of the derivatives are well defined, although it is not the mildest condition. Indeed, we require the functionals to have spatial derivatives on Ω rather than merely on the Λ̃(, t). 

	For the regular case, 
	
		
  * V, f, G♢ g, g ∈ C^1, 2_+ (Λ);
		
  * For any fixed s and y, f^s, y∈ C^1, 2_+ (Λ);
		
  * For any fixed r ∈ [0, T], c^r ∈ C^1, 2_+ ([0, r] ×Ω); and
		
  * For any fixed s, y, and fixed r ∈ [0, T], c^s, y, r∈ C^1, 2_+ ([0, r] ×Ω).
	
	For the singular case, let α∈ (0, 1).
	
		
  * V, f, G♢ g, g ∈ C^1, 2_+,α (Λ).
		
  * For any fixed s and y, f^s, y∈ C^1, 2_+,α (Λ);
		
  * For any fixed r ∈ [0, T], c^r ∈ C^1, 2_+,α ([0, r] ×Ω); and
		
  * For any fixed s, y, and any fixed r ∈ [0, T], c^s, y, r∈ C^1, 2_+,α ([0, r] ×Ω).
	
	In addition, β≜α + H - 1/2 > 0, where the constant H is defined in Definition <ref> (2) for the singular case.


We frequently encounter stochastic integrals that are required to be true martingales. Lemma <ref> is useful for the related justification. For ease of notation, we denote ∂_ω f, σ^t, · W  ≜∫^·_t ∂_ω f(r, X^⊗_r Θ^r, ), σ^r,  dW_r for later use.

	Suppose that  is admissible. Let f be a general functional and f ∈ C^1, 2_+ (Λ) for the regular case or f ∈ C^1, 2_+, α (Λ) for the singular case, with β≜α + H - 1/2 > 0. Then,
	
    [∫^T_t | ∂_ω f(r, X^⊗_r Θ^r, ), σ^r, |^2 dr | _t ] < ∞,

	which implies that ∂_ω f, σ^t, · W is a true martingale.


We can now provide the verification theorem, which is one of the main results of this paper. The proof is in the same spirit of <cit.> but invokes Lemmas <ref> and <ref> and the functional Itô formula in <cit.>. The proof consists of two steps. First, we show that the interpretations in Definition <ref> (6) hold and that V(t, ω) = J(t, ω; ). The functional Itô formula and Lemma <ref> are applied to prove the martingale property of three functions in Definition <ref> (6). V(t, ω) = J(t, ω; ) is verified in a similar way, under the conditions in Definition <ref>. Second, we prove that  is indeed an equilibrium strategy under Definition <ref>. The recursive relationship in Lemma <ref> with _h gives one representation of J(t, ω; _h). The PHJB equation (<ref>) with _h implies an inequality related to V(t, ω), which equals J(t, ω; ) according to the first step. A comparison of these two sides deduces (<ref>) as desired. In the whole proof, including Θ^t, recovers the flow property and overcomes the non-Markovian and non-semimartingale difficulty.


	Suppose that the extended PHJB system in (<ref>)-(<ref>) in Definition <ref> admits a solution (V, f, g, f^s, y, c^r, c^s, y, r) that satisfies Assumption <ref>. If  realizes the supremum in (<ref>) for V and  is admissible, then  is an equilibrium law in the sense of Definition <ref> and V is the corresponding value function.




§ EXAMPLES

In this section, we investigate the impact of volatility roughness. We apply the general framework to some specific decision-making situations with explicit or semi-closed form solutions. They include TC-MVP with constant risk aversion <cit.>, TC-MVP for log returns <cit.>, and the MV objective with a linear controlled Volterra process. We refer to TC-MVP with constant risk aversion as the const-MV case and to TC-MVP for log returns as the log-MV case. We focus on the Volterra Heston model, which is a specific form of SVIE (<ref>). 

The concept of resolvent is used frequently. Kernel R on [0, ∞) is referred to as the resolvent, or the resolvent of the second kind, of K if

    K*R(t) = R*K(t) = K(t) - R(t),   ∀  t ≥ 0,

where * denotes the convolution operation:

    K * R(t) = ∫^t_0 K(t - s) R(s) ds,   ∀  t > 0.

The integral is extended to t = 0 by right-continuity if possible. Further properties of these definitions can be found in <cit.>. Examples of kernels are available in Table <ref>.

Let R_λ be the resolvent of λ K such that 

    λ K * R_λ = R_λ * ( λ K) = λ K - R_λ.

If λ = 0, R_λ/λ = K and R_λ = 0.



Unlike variance, the wealth process (<ref>) does not have a convolution feature. Roughly speaking, certain Markov property is thus maintained. The dependence on wealth does not involve the entire trajectories. The following examples demonstrate this point.














 §.§ Const-MV: TC-MVP under constant risk aversion


Consider the TC-MVP in <cit.> under the Volterra Heston model (<ref>) and wealth (<ref>):

    _t [ M^_T ] - γ/2Var_t [M^_T] = _t [ M^_T - γ/2 (M^_T)^2] + γ/2( _t[M^_T] )^2,

where the constant γ > 0 reflects the risk aversion level. The general reward functional in (<ref>) then becomes

    F(t, ω_t, M^t, ω, _T ∧·) = M^_T - γ/2 (M^_T)^2,    G(t, ω_t, y) = γ/2 y^2.



To solve the PHJB equation system in Definition <ref>, we highlight that c^r = 0 and that f is not state-dependent. Consider the following Ansatz for V in (<ref>) and g in (<ref>). Denote the current wealth at time t by M and recall Θ^t_s as defined in (<ref>):

    V(t, ω)    =    V(t, M, Θ^t_[t, T])  = V_1(t) M + ∫^T_t V_2(s) Θ^t_s ds + V_0(t), 
    
    g(t, ω)    =    g (t, M, Θ^t_[t, T]) = g_1(t) M + ∫^T_t g_2(s) Θ^t_s ds + g_0(t),

where V_1, V_0, g_1, and g_0 are deterministic continuously differentiable functions and V_2 and g_2 satisfy suitable integrability conditions. This Ansatz implies that the functions V and g depend on current wealth M and Θ^t_[t, T] only. 

As V and g are linear functionals of Θ^t_[t, T], the direct calculation proceeds as follows:

    (^ V)(t, ω) =   V̇_1(t) M - V_2(t) ν + V̇_0(t) 
        + ( M + θ√(ν) u) V_1(t) + (κϕ - κν) ∫^T_t V_2(s)K(s-t)ds, 
    
    (^ g)(t, ω) =   ġ_1(t) M - g_2(t) ν + ġ_0(t) 
         + ( M + θ√(ν) u) g_1(t) + (κϕ - κν) ∫^T_t g_2(s)K(s-t)ds, 
    ^(G ♢ g)(t, ω) =   γ g (t, M, Θ^t_[t, T]) (^ g)(t, ω) + γ/2 g^2_1(t) u^2 
        + σρ u √(ν)γ g_1(t) ∫^T_t g_2(s) K(s-t) ds + γ/2( ∫^T_t g_2(s) K(s-t) ds )^2 σ^2 ν, 
    ∂_y G(t, ω_t, g(t,    ω)) = γ g (t, M, Θ^t_[t, T]).

We use the fact that ω^ν is continuous at time t and Θ^t_t = ν_t ≜ν.

Equation (<ref>) in Definition <ref> becomes

    sup_∈{   V̇_1(t) M - V_2(t) ν + V̇_0(t) + ( M + θ√(ν) u) V_1(t) + (κϕ - κν) ∫^T_t V_2(s)K(s-t)ds     - γ/2 g^2_1(t) u^2 -  σρ u √(ν)γ g_1(t) ∫^T_t g_2(s) K(s-t) ds     - γ/2( ∫^T_t g_2(s) K(s-t) ds )^2 σ^2 ν} = 0,    V_1(T) = 1,   V_0(T) = 0.

Therefore,

    (t, ν_t) = θ V_1(t) - γσρ g_1(t) ∫^T_t g_2(s)K(s-t)ds/γ g^2_1(t)√(ν_t).

Furthermore, we have g_1(T) = 1, g_0(T) = 0, and

    (^ g)(t, ω) =   ġ_1(t) M - g_2(t) ν + ġ_0(t) +  M g_1(t) + θ^2 ν V_1(t)/γ g_1(t)  - ρσθν∫^T_t g_2(s)K(s-t)ds 
        + (κϕ - κν ) ∫^T_t g_2(s)K(s-t)ds = 0.


By separation of variables and recognizing that g_1(t) = V_1(t) from (<ref>) and (<ref>), we obtain the following:

    ġ_1(t) + _t g_1(t) = 0,    g_1(T) = 1, 
        g_2(t) + (κ + ρσθ) ∫^T_t g_2(s)K(s-t)ds - θ^2/γ = 0, 
       ġ_0(t) + κϕ∫^T_t g_2(s)K(s-t)ds = 0,    g_0(T) = 0

and

    V̇_1(t) + _t V_1(t) = 0,    V_1(T) = 1, 
        V_2(t) + κ∫^T_t V_2(s)K(s-t)ds + γσ^2/2( ∫^T_t g_2(s) K(s-t) ds )^2         - (θ - γσρ∫^T_t g_2(s)K(s-t)ds)^2/2 γ= 0, 
       V̇_0(t) + κϕ∫^T_t V_2(s)K(s-t)ds = 0,    V_0(T) = 0.


The system in (<ref>)–(<ref>) can be solved explicitly. First, g_1(t) = V_1(t) = e^∫^T_t _sds. (<ref>) is a linear Volterra integral equation (VIE). Existence and uniqueness can be determined using <cit.>. Let λ = κ + ρσθ and recall that R_λ is the resolvent of λ K:

    g_2(t) = θ^2/γ - θ^2/γ∫^T_t R_λ(s-t)ds.

Additionally, a useful result is

    ∫^T_t g_2(s)K(s-t)ds = θ^2/γ∫^T_t R_λ(s-t)/λ ds.

g_0 can then be solved directly. V_2 in (<ref>) is also a linear VIE that can be solved in the same way as g_2. V_0 is solved after V_2. Although the calculation is straightforward, it is lengthy. As such, we omit it here. Ultimately,

    (t, ν_t) = θ/γ e^-∫^T_t _sds√(ν_t) -  ρσθ^2/γ e^-∫^T_t _sds∫^T_t R_λ(s-t)/λ ds √(ν_t).


Then the support for the wealth process is . The first term in  (<ref>) is the same as the constant volatility case. The second term can be interpreted as a hedge for the randomness from stochastic volatility. Roughness alters the hedge through the resolvent R_λ.

Verifying that  in (<ref>) is admissible in the sense of Definition <ref> is straightforward. Indeed, Assumption <ref> holds in view of the moment estimates in <cit.> for ν. Other requirements in Definition <ref> are direct. We summarize the analysis above in the following lemma.
 
	Problem (<ref>) under the Volterra Heston model (<ref>) has an equilibrium strategy  given by (<ref>), which is admissible in the sense of Definition <ref>. The value function is given by (<ref>), with V_1, V_2, and V_0 given by (<ref>), (<ref>), and (<ref>), respectively. g in (<ref>) is given by (<ref>), with g_1, g_2, and g_0 given by (<ref>), (<ref>), and (<ref>), respectively.




 §.§ Log-MV: TC-MVP for log returns

Instead of considering preferences for terminal wealth, <cit.> argues that an analysis based on log returns is more plausible. The derived equilibrium strategy is wealth-dependent and will not short sell risky assets with positive excess returns over the long-term horizon. Suppose that the proportional amount of wealth in the stock is . M^ is the wealth process corresponding to . Denote L_t^≜ln M^_t. For ease of notation, we write L_t = L_t^. It follows that 

    d L_t = [_t + θν_t π_t - 1/2π^2_t ν_t ] dt + √(ν_t)π_t dW_1t.


Consider the TC-MVP in <cit.> under the Volterra Heston model (<ref>) and log return (<ref>):

    _t [ L_T ] - γ/2Var_t [L_T] = _t [ L_T - γ/2 (L_T)^2] + γ/2( _t[L_T] )^2.


With a slight abuse of notation, we try the following Ansatz for V in (<ref>) and for g in (<ref>):

    V(t, ω)    =    V(t, L, Θ^t_[t, T])  = L + ∫^T_t V_2(s) Θ^t_s ds + V_0(t), 
    
    g(t, ω)    =    g (t, L, Θ^t_[t, T]) = L + ∫^T_t g_2(s) Θ^t_s ds + g_0(t),

where V_0 and g_0 are deterministic continuously differentiable functions and V_2 and g_2 satisfy the suitable integrability conditions.

Equation (<ref>) in Definition <ref> becomes

    sup_∈{    - V_2(t) ν + V̇_0(t) +  + θνπ - 1/2νπ^2 + (κϕ - κν) ∫^T_t V_2(s)K(s-t)ds     - γ/2νπ^2 -   γρσνπ∫^T_t g_2(s) K(s-t) ds     - γ/2( ∫^T_t g_2(s) K(s-t) ds )^2 σ^2 ν} = 0,    V_0(T) = 0.

Therefore,

    π̂_t = θ - γρσ∫^T_t g_2(s)K(s-t)ds/1 + γ.

Furthermore, (^ g)(t, ω)=0 gives rise to

    - g_2(t) ν + ġ_0(t) +  +  θν/1 + γ[θ - γρσ∫^T_t g_2(s)K(s-t)ds ]   
         -ν/2(1+γ)^2[θ - γρσ∫^T_t g_2(s)K(s-t)ds]^2 + (κϕ - κν ) ∫^T_t g_2(s)K(s-t)ds = 0, 
          g_0(T) = 0.


By separation of variables, we obtain

    g_2(t) =    (1+2γ)θ^2/2(1+γ)^2 - [κ + γ^2 ρσθ/(1+γ)^2] ∫^T_t g_2(s)K(s-t)ds 
        - γ^2ρ^2σ^2/2(1+γ)^2(∫^T_t g_2(s)K(s-t)ds )^2.

Let ψ(T-t) = ∫^T_t g_2(s)K(s-t)ds, then convolve both sides of (<ref>) with kernel K(·) and change T-t to t. This yields the following Riccati-Volterra equation:

    ψ(t) = ∫^t_0 K(t - s) [ - γ^2ρ^2σ^2/2(1+γ)^2ψ^2(s) - (κ + γ^2 ρσθ/(1+γ)^2) ψ(s) + (1+2γ)θ^2/2(1+γ)^2] ds.


Furthermore,

    V_2(t) + κ∫^T_t V_2(s)K(s-t)ds + γσ^2/2ψ^2(T-t) - 1/2(1+γ)[θ - γρσψ(T-t)]^2 = 0, 
       V̇_0(t) + _t + κϕ∫^T_t V_2(s)K(s-t)ds = 0,    V_0(T) = 0, 
       ġ_0(t) + _t + κϕψ(T-t) = 0,    g_0(T) = 0.



	Suppose that κ + γ^2 ρσθ/(1+γ)^2 > 0. Then, (<ref>) has a unique global continuous solution on [0, T]. Define 
	
    H(w) ≜γ^2ρ^2σ^2/2(1+γ)^2  w^2 - (κ + γ^2 ρσθ/(1+γ)^2) w - (1+2γ)θ^2/2(1+γ)^2≜ H_2 w^2 +  H_1 w + H_0.

	Then,
	
    0 < ψ(t) ≤ -r_1(t) < - w_* ,   ∀  t > 0,

	with w_* = -H_1 - √(H_1^2 - 4H_2H_0)/2H_2 < 0 and r_1(t) ≜ Q^-1_1 ( ∫^t_0 K(s) ds ), where Q_1(w) = -∫^0_w du/H(u).
	
		In addition, system (<ref>)-(<ref>) has a unique continuous solution (V_2, V_0, g_0) on [0, T].
	



Ultimately, an equilibrium strategy is expressed as follows:

    π̂_t = θ/1 + γ - γρσ/1 + γψ(T-t).


For the admissibility of , we have the following result about Assumption <ref>.

	Assume that (<ref>) has a unique continuous solution on [0, T]. Suppose that
	
    [ e^c ∫^T_0 ν_s ds] < ∞,

	with constant c expressed as follows:
	
    c = max{ 2 p |θ| sup_t ∈ [0, T] | π̂_t |, (8p^2 - 2p) sup_t ∈ [0, T]π̂^2_t },   for certain  p > 1.

	Then, wealth M^* under  satisfies 
	
    [ sup_ t ∈ [0, T] | M^*_t |^p ] < ∞.


Thus, we have the following lemma.
 
	Problem (<ref>) under the Volterra Heston model (<ref>) has an equilibrium strategy  given by (<ref>). If Assumption (<ref>) holds for a large enough constant c>0, then the equilibrium strategy (<ref>) is admissible in the sense of Definition <ref>. The value function is given by (<ref>), with V_2 and V_0 given by (<ref>) and (<ref>), respectively. g in (<ref>) is given by (<ref>), with g_2 and g_0 given by (<ref>) and (<ref>), respectively.


	For the specific fractional kernel K(t) = t^H-1/2/Γ(H+1/2), Assumption (<ref>) holds for a large enough constant c>0 when the time horizon T is sufficiently small. See <cit.>.


	For the general Heston specification considered in <cit.>, 
	
    dS_t = S_t (_t + θν^1+δ/2δ_t) dt + S_t ν^1/2δ_t dW_1t.

	An equilibrium strategy is
	
    π̂_t = [ θ/1 + γ - γρσ/1 + γψ(T-t) ] ν^δ - 1/2δ_t.

	The related proof and verification of admissibility are identical.




 §.§ TC-MV under linear controlled Volterra processes

	In this section, we consider the same objective (<ref>) as in Section <ref>: 
	
    _t [ X^_T - γ/2 (X^_T)^2] + γ/2( _t[X^_T] )^2.

	However, the state process X^ is a one-dimensional linear controlled Volterra process given by (<ref>). B, b, D, and σ are one-dimensional essentially bounded and deterministic measurable functions:
	
    X^_t = x + ∫^t_0 K(t-r)(B_r u_r + b_r)dr + ∫^t_0 K(t-r)(D_r u_r + σ_r) dW_r.

	
	Consider the following Ansatz for V and g. Θ^t, is defined as in (<ref>) but with =.
	
    V(t, ω)    =    V_1(t) Θ^t, _T + V_0(t), 
    
    	g(t, ω)    =    g_1(t) Θ^t, _T + g_0(t),

	where V_1, V_0, g_1, and g_0 are deterministic continuously differentiable functions.
	
	Similarly, we have
	
    (^ V)(t, ω) =   V̇_1(t) Θ^t, _T + V̇_0(t) + V_1(t) K(T-t) (B_t u_t + b_t), 
    
    	(^ g)(t, ω) =   ġ_1(t) Θ^t, _T + ġ_0(t) + g_1(t) K(T-t) (B_t u_t + b_t), 
    ^(G ♢ g)(t, ω) =   γ g(t, ω) (^ g)(t, ω) + γ/2 g^2_1(t) K^2(T - t) (D_t u_t + σ_t)^2.

	
	Equation (<ref>) in Definition <ref> is
	
    sup_∈{   V̇_1(t) Θ^t, _T + V̇_0(t) + V_1(t) K(T-t) (B_t u_t + b_t)     - γ/2 g^2_1(t) K^2(T - t) (D_t u_t + σ_t)^2  } = 0,    V_1(T) = 1,   V_0(T) = 0.

	Along with (^ g)(t, ω) =0, we obtain the following:
	
    V_1(t)     = g_1(t) = 1, 
    V̇_0(t)     + (B_t - γ K(T-t) D_t σ_t)^2/2 γ D^2_t + K(T-t)b_t - γ/2 K^2(T-t)σ^2_t = 0,    V_0(T) = 0, 
    ġ_0(t)     + B^2_t - γ K(T-t) D_t B_t σ_t/γ D^2_t + K(T-t) b_t = 0,    g_0(T) = 0.

	An equilibrium control is
	
    _t = B_t - γ K(T- t) D_t σ_t/γ K(T- t) D^2_t.

	When K= id, the solution reduces to that found in <cit.>. Verifying the admissibility of (<ref>) is straightforward, as follows:
	
		Suppose that Assumption <ref> holds and D^2_t ≥δ >0. Then, (<ref>) and (<ref>) admit unique solutions on [0, T]. (<ref>) is an admissible equilibrium strategy in the sense of Definition <ref> for the MV objective (<ref>) under the state (<ref>). 
	



§ NUMERICAL ANALYSIS

In this section, we numerically analyze the effects of roughness on stock demand and improvements on profitability. We concentrate on the rough Heston model <cit.>, namely kernel K(t) = t^H-1/2/Γ(H+1/2), H ∈ (0, 1/2]. We conduct this numerical study from three perspectives. First, we perform a sensitivity analysis by varying the roughness only and fixing the other parameters. However, roughness can have sophisticated interactions with other variables. Thus, second, we further consider a calibration situation. Third, we conduct an empirical study based on recent market data. Strategies with rough volatility outperform their classic counterparts in terms of profit and the Sharpe ratio. We mainly compare four strategies: the const-MV case (<ref>), the log-MV case (<ref>), the pre-committed MVP <cit.>, and Merton's portfolio problem with power utility (CRRA) <cit.>.

The interaction between time-inconsistency and volatility roughness can be complicated. We answer three questions. First, when volatility is rougher, should investors increase or decrease their stock demand? Second, do investors with different levels of risk aversion behave differently with respect to volatility roughness? Third, is it more profitable to incorporate rough volatility?



 §.§ Sensitivity analysis



  §.§.§ Const-MV

The equilibrium strategy in Lemma <ref> can be further simplified for the fractional kernel K(t) = t^α-1/Γ(α) with α = H + 1/2. For λ  = 0,

    ∫^T_t R_λ(s-t)/λ ds = ∫^T_t K(s - t) ds = (T - t)^α/Γ(α + 1).

For λ≠ 0, the property of Mittag-Leffler functions in <cit.> shows that

    ∫^T_t R_λ(s-t)/λ ds = 1 - E_α, 1 ( - λ (T - t)^α)/λ≜F^α, λ(T - t)/λ.

This enables us to compute the equilibrium strategy using an explicit formula. 

Figure (<ref>) illustrates the plot for the hedging term,

    -  ρσθ^2/γ e^-∫^T_t _sds∫^T_t R_λ(s-t)/λ ds.

Note that we set a negative correlation between stock price and volatility in the plot to reflect the leverage effect. Figure <ref> numerically shows the following phenomenon. When the investment horizon is long (i.e., t is small), the const-MV strategy advocates investing more if the stock is smoother. However, near the end of the investment horizon, the const-MV strategy prefers to invest more if the stock is rougher. We refer to this phenomenon as the investment horizon effect. It is distinct from the pre-committed MVP strategy obtained in <cit.>. When stock volatility is smooth, the equilibrium strategy reduces the stock position gradually until the end of the investment period. In contrast, when stock volatility is rough, the equilibrium strategy suggests a relatively steady holding of the stock for a sufficiently long investment horizon but rapidly slashes the holding near the end of the investment horizon. The latter phenomenon also occurs in optimal investment problems with position limits, but our problem has no constraints on the position. In fact, the investment horizon effect can be shown mathematically by deriving asymptotic estimates of (<ref>) and (<ref>). We refer to the following corollary.


	Suppose that α∈ (1/2, 1). Then, for sufficiently large values of T - t, integral ∫^T_t R_λ(s-t)/λ ds is increasing on α. For sufficiently small values of T - t, ∫^T_t R_λ(s-t)/λ ds is decreasing on α.




The behavior in Figure <ref> can be further enhanced using different investment horizons. In Figure <ref>, an investor with a relatively short horizon (e.g., 1 year) would simply have more stock demands if volatility is rougher. If the investment horizon is long enough (e.g., 10 years) and volatility is rougher, the investor purchases less for the first 8 years and purchases more for the remaining time. The diminished amount of wealth invested in the stock can be up to 40% of the classic Heston counterpart. If other variables do not change, increments in volatility roughness reduce the willingness to retain more stock exposure in the long term. Roughness does have a material impact on the investment decisions of the investors with long planning horizons.






  §.§.§ Log-MV

For the fractional kernel, the Riccati-Volterra equation (<ref>) can be solved numerically using the fractional Adams method, which is detailed in <cit.>. The assumption in Corollary <ref> is validated for settings in Figures <ref> and <ref>. Assumption (<ref>) is satisfied by Remark <ref>.

The investment horizon effect is also observed in log-MV equilibrium strategies. When volatility is rougher, investors purchase less at first and more later on, as shown in Figure <ref>. This is consistent with the const-MV counterpart. However, there is an important difference between the two criteria. For the const-MV case, the moment to prefer rough, shown as the interactions between the curves in Figure <ref>, is not affected by heterogeneity in risk aversion γ. This is not the case for the log-MV case. As shown in Figure <ref>, when log-MV investors are more risk averse (i.e., the risk aversion parameter γ is larger), they prefer rough earlier. This indicates that if only roughness varies and increases, then preferring rough earlier can minimize risk.





  §.§.§ Comparison with pre-committed MVP and CRRA utility

We first summarize the comparison of parameter dependence for the const-MV, log-MV, pre-committed MVP, and CRRA utility cases. All four strategies depend on the risk premium θ, risk aversion[Pre-committed MVP is replaced by target terminal wealth.] γ, the correlation ρ between stock and volatility, volatility-of-volatility σ, the mean reversion speed κ, the investment horizon T, and the Hurst parameter H. For simplicity, we refer to these seven parameters as the primary parameters. The CRRA and log-MV strategies have the simplest parameter dependence. Only the primary parameters affect stock demand. The const-MV strategy relies on an additional parameter, namely the risk-free rate. The pre-committed MVP strategy depends on all parameters. Interestingly, the pre-committed MVP strategy is the only strategy that depends on the long-term mean level of volatility. Table <ref> summarizes the comparison along with wealth dependence.

<cit.> observes that volatility-of-volatility has a material impact on the pre-committed MVP in the sensitivity analysis. However, it is not observed for const-MV and log-MV investors. Alternately, the investment horizon becomes vital for the time-consistent alternatives. 

<cit.> argues that log-MV and CRRA criteria are analogous. Their structure of strategies is almost the same. This similarity is maintained when volatility is rough, but the effect of roughness is disparate. In short investment horizons (e.g., 1 year), CRRA investors allocate less to the stock when volatility is rougher <cit.>. Log-MV investors do the opposite. In long investment horizons (e.g., 10 years), CRRA and log-MV investors still have distinct preferences regarding roughness. 





  §.§.§ Linear controlled Volterra processes
 



To further enrich the analysis, we include a numerical study on the objective (<ref>) under the state process (<ref>). For simplicity, consider constant parameters B, D, σ and the fractional kernel. The equilibrium control (<ref>) then reduces to

    _t = B/γ D^2 (T - t)^1/2 - HΓ(H+1/2) - σ/D.

When H = 1/2,  _t is a constant. If the state process is rough with H<1/2, then _t is larger at the beginning and smaller at the end, as shown in Figure <ref>. Moreover, this phenomenon is the same for different time horizons. As the fractional kernel appears together with the control, it exhibits distinct effects compared with the rough Heston model. When the agent is more risk-averse (a larger γ), _t is smaller. In the extreme case with γ→∞, the effect of roughness is eliminated and _t → - σ/D.



 §.§ Simulation study

In the sensitivity analysis, we only vary the roughness level. In general, this is unrealistic as other parameters are also likely to change. For example, <cit.> document the connection between volatility roughness and components with fast mean reversion. It is also observed that to capture the deep near-term volatility skew, calibration with the classic Heston model usually results in a greater mean reversion speed κ and in a greater volatility-of-volatility σ.

In this section, we leverage the information from implied volatility (IV). Theoretically, the roughness estimates based on realized and implied volatility should coincide. However, this does not hold in reality. The IV surface represents the current view of the future. In particular, ATM skew explosion indicates near-term downside risks, such as earnings announcements <cit.>. As the real parameters underlying the IV surface are unknown, we first adopt the simulated IV surface in <cit.>. Given the simulated IV data, investors calibrate two sets of parameters for the Heston and rough Heston models. We contrast the two strategies induced from the calibrated parameters. The investor under the Heston model (hereafter, Heston investor) uses the calibrated parameters in <cit.>. The investor under the rough Heston model (hereafter, rough investor) uses <cit.>. The variance process is simulated using the lifted Heston method in <cit.>. The parameters for simulation are given in <cit.>. We set M_0 = 1, = 0.01, θ = 1.5, T = 10, and γ = 0.5. Additionally, we implement the Euler scheme for the stock process. The simulation is run with 250 time steps for 1 year, corresponding to the 250 trading days per year.  



Figure (<ref>) plots the dollar amount in the stock for the const-MV investor, and Figure (<ref>) depicts the proportional amount of wealth in the stock for the log-MV investor. In both cases, the investors demand more if the rough Heston model is adopted. Volatility roughness dramatically changes the investment decisions. The rough const-MV investor almost doubles the Heston counterpart. The rough log-MV investor increases stock demand by nearly 10%. These results are interesting. The advantage of the rough Heston model is that it better captures the volatility smiles with short maturities. In other words, the rough Heston model captures the near-term downside risk, whereas the classic Heston model fails to do so. Surprisingly, the downside risk alters the investment over almost the entire time horizon. We interpret these adjustments as hedging against risk.





Figures <ref> and <ref> demonstrate the distribution of terminal wealth for the const-MV and log-MV investors with the rough Heston model or the classic Heston model. Figures <ref> and <ref> have 5,000 simulation paths. Rough investors tend to obtain higher terminal wealth but with higher variance. The reason for this may be that they bear the near-term risk represented by roughness. The log-MV case in Figure <ref> looks conservative compared with the const-MV case. However, the const-MV case in Figure <ref> and the log-MV case in Figure <ref> are not fairly comparable. Risk aversion γ is set to 0.5, but one is for wealth and another is for log return. The scale makes the actual risk aversion level different.



 §.§ Empirical study

To further elucidate the need to account for rough volatility, we evaluate the performance of the strategies on real financial data. We download the CBOE S&P 500 options data from OptionMetrics via Wharton Research Data Services. The most recent 1-year time period we can obtain is from 2018/05/01 to 2019/05/01. The U.S. equity market was relatively volatile during this period, plunging in October 2018 and December 2018. Due to the computational burden involved in calibrating the rough Heston model, we adopt the following procedure to implement the trading strategies. Consider the investment horizon as 1 month. The models are calibrated on European call options data observed on the first trading day of that month. We then use the calibrated parameters for trading in the whole month. The portfolios are rebalanced at a daily frequency (i.e., the step size equals 1/250). For simplicity, we ignore transaction costs and price impact. After this month, we redo the calibration and implement the trading rules for the next 1-month investment horizon. Overall, for each model, we have 12 sets of calibrated parameters in a year. Tables <ref> and <ref> report the sample means and sample standard deviations (std) of the calibration results. In Table <ref>, the calibrated Hurst parameter is around 0.15 with a slight variance. Therefore, volatility roughness is a persistent fact in the equity market. In Table <ref>, it is further confirmed that the classic Heston model tends to yield a considerably greater mean reversion speed κ and volatility-of-volatility σ.





In Tables <ref> and <ref>, we present the terminal wealth and Sharpe ratio at several risk aversion levels. In the const-MV and log-MV cases, trading strategies under rough volatility dominate the classic Heston counterparts <cit.> with both higher terminal wealth and a better Sharpe ratio. This confirms that the rough Heston model captures profit opportunities missed by the classic Heston model. Rough volatility models are therefore an attractive alternative for volatility modeling. Furthermore, trading strategies under the log-MV criterion are less sensitive to the choice of risk aversion. Figure <ref> illustrates the special case in which risk aversion equals 2. Const-MV strategies yield higher profits, whereas log-MV strategies are more stable, even during the market downturn in October and December 2018. Log-MV strategies yield lower Sharpe ratios mainly due to lower excess returns. In addition, we emphasize that directly comparing the const-MV and log-MV cases is not fair as the actual risk aversion preferences are distinct.











§ CONCLUDING REMARKS

In this paper, we use the functional Itô calculus to examine in depth the equilibrium strategies under time-inconsistent preferences and in a rough stochastic environment. Volatility roughness significantly modifies investment demands. Our general framework also embraces Volterra processes with potential applications beyond rough volatility. Several interesting problems are left for future research. The first is the existence and uniqueness of solutions to the extended PHJB equation system in Definition <ref>. The second is the time-inconsistent open-loop control problem under Volterra processes.




§ BRIEF SUMMARY OF THE FUNCTIONAL ITÔ CALCULUS (VIENS AND ZHANG, 2019)


Let Ω≜ C^0([0, T], ^n) be the sample space with continuous paths, Ω̅≜ D^0([0, T],  ^n) be the sample space with càdlàg (right-continuous with left limits) paths, and

    Ω_t  ≜ C^0([t, T], ^n),   Λ≜ [0, T] ×Ω,   Λ̅≜{ (t, ω) ∈ [0, T] ×Ω̅: ω|_[t, T]∈Ω_t }, 
        ||ω||_T ≜sup_ 0 ≤ t ≤ T |ω_t|,   ((t, ω), (t', ω')) ≜ |t - t'| + ||ω - ω'||_T.


Let C^0(Λ̅) be the space of functions f: Λ̅→, which are continuous under . For f ∈ C^0(Λ̅), define the time derivative as follows:

    ∂_t f(t, ω) ≜lim_δ↓ 0f(t + δ, ω) - f(t, ω)/δ,    for all  (t, ω) ∈Λ̅,

whenever the limit exists.

Given (t, ω) ∈Λ̅, the spatial derivative with respect to ω, denoted by ∂_ω f(t, ω), is a linear operator on Ω_t and defined as the Fréchet derivative with respect to ω𝕀_[t, T]. That is,

    f(t, ω + η𝕀_[t, T]) - f(t, ω) = ∂_ω f(t, ω), η + o(||η𝕀_[t, T]||_T),    for any η∈Ω_t.

Furthermore, ∂_ω f(t, ω) satisfies the definition of the Gateaux derivative:

    ∂_ω f(t, ω), η = lim_ε→ 0f(t, ω + εη𝕀_[t, T]) - f(t, ω)/ε,    for any η∈Ω_t.

The perturbation is on [t, T], not on [0, t). If η∈Ω_s for certain s < t, the derivative is understood as follows:

    ∂_ω f(t, ω), η≜∂_ω f(t, ω), η𝕀_[t, T].


The second-order derivative ∂^2_ωω f(t, ω) is defined as a bilinear operator on Ω_t ×Ω_t:

    ∂_ω f(t, ω + η_1 𝕀_[t, T]), η_2  - ∂_ω f(t, ω), η_2 = ∂^2_ωω f(t, ω), (η_1, η_2) + o(||η_1 𝕀_[t, T]||_T),
 
for any η_1, η_2 ∈Ω_t. If η_1, η_2 ∈Ω_s for certain s < t, the derivative is understood in the same way as in (<ref>).

For the well-posedness of these derivatives, we refer readers to <cit.>.

We introduce two spaces, namely C^1, 2_+ (Λ) and C^1, 2_+, α (Λ), from <cit.>, under which the functional Itô formula in <cit.> holds.

	Suppose that f ∈ C^0(Λ̅) and ∂_ω f exists for all (t, ω) ∈Λ̅. 
	
		
  * ∂_ω f is said to have polynomial growth if there exist constants C_0, m > 0 such that
		
    |∂_ω f(t, ω), η| ≤ C_0 [1 + ||ω||^m_T] ||η𝕀_[t, T]||_T,   ∀  (t, ω) ∈Λ̅,  η∈Ω.
 
		
  * ∂_ω f is said to be continuous if, for all η∈Ω, the mapping (t, ω) ∈Λ̅↦∂_ω f(t, ω), η is continuous under .
		
  * ∂^2_ωω f is said to have polynomial growth if there exist constants C_0, m > 0 such that
		
    |∂^2_ωω    f(t, ω), (η_1, η_2)| 
       ≤ C_0 [1 + ||ω||^m_T] ||η_1𝕀_[t, T]||_T ||η_2 𝕀_[t, T]||_T,  ∀  (t, ω) ∈Λ̅,  η_1, η_2 ∈Ω.
 
		
  * ∂^2_ωω f is said to be continuous if, for all η_1, η_2 ∈Ω, the mapping (t, η_1, η_2) ∈Λ̅_2 ↦∂^2_ωω f(t, ω), (η_1, η_2) is continuous under '((t, ω_1, ω_2), (t', ω'_1, ω'_2)) ≜ |t - t'| + ||ω_1 - ω'_1||_T + ||ω_2 - ω'_2||_T, where Λ̅_2 ≜{ (t, ω_1, ω_2) ∈ [0, T] ×Ω̅×Ω̅: ω_1|_[t, T],  ω_2 |_[t, T]∈Ω_t }.
	


	Let C^1, 2(Λ̅) ⊂ C^0(Λ̅) be the set of all f with continuous derivatives ∂_t f, ∂_ω f, ∂^2_ωω f on Λ̅. Let C^1, 2_+(Λ̅) be the set of all f ∈ C^1, 2(Λ̅) such that all derivatives have polynomial growth and ∂^2_ωω f(t, ω), (η, η) is locally uniformly continuous in ω with polynomial growth, namely there exists a constant m > 0 and a bounded modulus of continuity function ϱ. For all (t, ω), (t, ω') ∈Λ̅ and η∈Ω_t, we have
	
    |∂^2_ωω f(t, ω) - ∂^2_ωω f(t, ω'), (η, η)| ≤[1 + ||ω||^m_T + ||ω'||^m_T ] || η𝕀_[t, T] ||^2_T ϱ(||ω - ω'||_T).

	C^1, 2_+(Λ) is defined in the same spirit of C^1, 2_+(Λ̅), with Λ̅ replaced by Λ.


	f ∈ C^1, 2_+ (Λ) is said to vanish diagonally at a rate of α∈ (0, 1), denoted by f ∈ C^1, 2_+, α (Λ), if there exists an extension of f in C^1, 2_+ (Λ̅), still denoted by f, such that for every 0 ≤ t < T, 0 < δ≤ T -t, and η, η_1, η_2 ∈Ω_t with supports contained in [t, t + δ]:
	
		
  * ∀ ω∈Ω̅ satisfying ω𝕀_[t, T]∈Ω_t,
		
    |∂_ω f(t, ω), η|    ≤ C_0 [1 + ||ω||^m_T] ||η||_T δ^α, 
    |∂^2_ωω f(t, ω), (η_1, η_2)|    ≤ C_0 [1 + ||ω||^m_T] |||η_1||η_2|||_T δ^2α;

		
  * For any other ω' ∈Ω̅ satisfying ω'𝕀_[t, T]∈Ω_t,
		
    |∂_ω f(t, ω)     - ∂_ω f(t, ω'), η| 
       ≤[1 + ||ω||^m_T + ||ω'||^m_T ] || η||_T ϱ(||ω - ω'||_T) δ^α, 
    |∂^2_ωω f(t, ω)     - ∂^2_ωω f(t, ω'), (η_1, η_2)| 
       ≤[1 + ||ω||^m_T + ||ω'||^m_T ] |||η_1||η_2|||_T ϱ(||ω - ω'||_T) δ^2α.

	
	Constant m > 0 denotes the polynomial growth rate and ϱ is a bounded modulus of the continuity function.

α characterizes the level of singularity in the diagonal of time. Finally, the functional Itô formula is quoted in Theorem <ref>.




















§ PROOFS OF RESULTS



 §.§ Proof of Lemma <ref>


	By the tower property of the conditional expectation and the definitions of c^r,, f^, g^, and ω^t+h,  
	
    J(t, ω; ) =    ∫^T_t [ [C(t, ω_t, r,  X^t, ω, _, (r, X^t, ω, _)) | _t+h] | _t ] dr      + [ [ F(t, ω_t, X^t, ω, _T ∧·) | _t+h] | _t ] + G(t, ω_t, [ [ X^t, ω, _T | _t+h] | _t] ) 
    	=    ∫^T_t [ c^r, (t+h, ω^t + h, t, ω_t)| _t ] dr  
        + [ f^(t+h, ω^t+h, t, ω_t) | _t ] + G(t, ω_t, [ g^(t+h, ω^t+h) | _t] ).

	Meanwhile, the definition of the reward functional in (<ref>) indicates the following:
	
    J(t+h, ω^t+h; ) =   [ ∫^T_t+h C(t+h, ω^t+h_t+h, r,  X^t+h, ω^t+h, _, (r, X^t+h, ω^t+h, _)) dr  | _t+h]    + [ F(t+h, ω^t+h_t+h, X^t+h, ω^t+h, _T ∧·) | _t+h]     + G(t+h, ω^t+h_t+h, [ X^t+h, ω^t+h, _T | _t+h]) 
    	=   ∫^T_t+h c^r, (t+h, ω^t+h, t+h, ω^t+h_t+h) dr + f^(t+h, ω^t+h, t+h, ω^t+h_t+h)     + G(t+h, ω^t+h_t+h, g^(t+h, ω^t+h)),

	where
	
    X^t+h, ω^t+h, _s    =   ω^t+h_s + ∫^s_t+hμ(s; r, X^t+h, ω^t+h, _, (r, X^t+h, ω^t+h, _)) dr        + ∫^s_t+hσ(s; r, X^t+h, ω^t+h, _, (r, X^t+h, ω^t+h, _)) dW_r,   t+h ≤ s ≤ T, 
    	X^t+h, ω^t+h, _s    =   ω^t+h_s,   0 ≤ s < t+h.

	
	Taking the conditional expectation at _t on both sides of (<ref>) yields the following:
	
    [ J(t+h, ω^t+h; ) | _t ] =   ∫^T_t+h[ c^r, (t+h, ω^t+h, t+h, ω^t+h_t+h) | _t ] dr 
        + [ f^(t+h, ω^t+h, t+h, ω^t+h_t+h) | _t ] 
        + [ G(t+h, ω^t+h_t+h, g^(t+h, ω^t+h)) | _t ].

	The result follows by combining (<ref>) and (<ref>).



 §.§ Proof of Lemma <ref>


	Let m be a generic positive value that may vary from line to line. We first present the proof for the regular case. By (1)-(2) in Definition <ref> and the assumption that ∂_ω f has polynomial growth,
	
    [∫^T_t | ∂_ω f(r, X^⊗_r Θ^r, ), σ^r, |^2 dr | _t ] 
       ≤ C_0 [ sup_ t ≤ r ≤ T| ∂_ω f(r, X^⊗_r Θ^r, ), σ^r, |^2 | _t ]    ≤ C_0 [( 1 + sup_ t ≤ r ≤ Tsup_0 ≤ s ≤ T |(X^⊗_r Θ^r, )_s|^m )^2  sup_ t ≤ r ≤ Tsup_ r ≤ s ≤ T | σ^(s; r, X^_)|^2 | _t ]    ≤ C_0 [ (1 + sup_ t ≤ r ≤ Tsup_0 ≤ s ≤ T |(X^⊗_r Θ^r, )_s|^m ) | _t ]    ≤ C_0 [ (1 + sup_ 0 ≤ s ≤ T |X^_s|^p )| _t ] < ∞.

	For the singular case, for [r, T], consider the partition r = r_∞ < ... < r_k < ... < r_0 = T, where r_k = r + T-r/2^k. Then,
	
    [∫^T_t | ∂_ω f(r, X^⊗_r Θ^r, ), σ^r, |^2 dr | _t ] 
        = [∫^T_t | lim_δ↓ 0∂_ω f(r, X^⊗_r Θ^r, ), σ^δ, r, |^2 dr | _t ]     = [∫^T_t | lim_δ↓ 0∑^∞_k = 0∂_ω f(r, X^⊗_r Θ^r, ), σ^δ, r, _s 𝕀_s ∈ [r_k+1, r_k)|^2 dr | _t ]    ≤ C_0 [∫^T_t (1 + sup_0 ≤ s ≤ T |(X^⊗_r Θ^r, )_s|^m )^2 | lim_δ↓ 0∑^∞_k = 0 ||σ^δ, r, _s 𝕀_s ∈ [r_k+1, r_k)||_T (r_k - r_k+1)^α|^2 dr | _t ] ,

	where we use the assumption that f vanishes diagonally at a rate of α∈ (0, 1) and the fact that ∂_ω f is a linear operator in the last inequality. By (2) in Definition <ref> for the singular case,
	
    ∑^∞_k = 0 ||σ^δ, r, _s 𝕀_s ∈ [r_k+1, r_k)||_T (r_k - r_k+1)^α
       ≤  C_0 (1 + sup_0 ≤ s ≤ T |(X^⊗_r Θ^r, )_s|^m ) ∑^∞_k = 0( r_k+1∨ (r+δ) - r )^H - 1/2( T-r/2^k+1)^α.

	For any 0 < δ≤ T - r, there exists an integer z such that T-r/2^z+1 < δ≤T-r/2^z. Then, 
	
    ∑^∞_k = 0(r_k+1∨ (r+δ) - r )^H - 1/2( T-r/2^k+1)^α    = ∑^z-1_k=0(T-r/2^k+1)^H - 1/2( T-r/2^k+1)^α + ∑^∞_k = zδ^H - 1/2( T-r/2^k+1)^α    = ∑^z-1_k=0(T-r/2^k+1)^β + δ^H - 1/2(T-r/2^z+1)^α∑^∞_k = 0( 1/2^α)^k.

	Note that T-r/2^z+1 < δ implies that (T-r/2^z+1)^α < δ^α. We obtain the following:
	
    ∑^∞_k = 0 ||σ^δ, r, _s 𝕀_s ∈ [r_k+1, r_k)||_T (r_k - r_k+1)^α   ≤ C_0 (1 + sup_0 ≤ s ≤ T |(X^⊗_r Θ^r, )_s|^m ) [ (T-r)^β + δ^β].

	Finally,
	
    [∫^T_t | ∂_ω f(r, X^⊗_r Θ^r, ), σ^r, |^2 dr | _t ] 
       ≤ C_0 [∫^T_t (1 + sup_0 ≤ s ≤ T |(X^⊗_r Θ^r, )_s|^m )^4 (T - r)^2 β dr | _t ]    ≤ C_0 [ (1 + sup_ t ≤ r ≤ Tsup_0 ≤ s ≤ T |(X^⊗_r Θ^r, )_s|^m ) | _t ]    ≤ C_0 [ (1 + sup_ 0 ≤ s ≤ T |X^_s|^p )| _t ] < ∞,

	as desired.




 §.§ Proof of Theorem <ref>


	First, we show that the interpretations in Definition <ref> (6) hold and V(t, ω) = J(t, ω; ).
	
	By (<ref>), (<ref>), (<ref>) and Lemma <ref>, f^s,y(t, X^⊗_t Θ^t, ), g(t, X^⊗_t Θ^t, ), and c^s, y, r(t, X^⊗_t Θ^t, ) are martingales. By the boundary conditions in (<ref>), (<ref>), and (<ref>) and note that ω = X^t, ω, ⊗_t Θ^t,, we derive the following:
	
    f^s, y(t, ω)     = [ F(s, y, X^t, ω, _T ∧·) | _t ],      g(t, ω) = [X^t, ω, _T | _t ], 
    
    	c^s, y, r(t, ω)    = [ C(s, y, r,  X^t, ω, _, (r,  X^t, ω, _)) | _t ],    0 ≤ t ≤ r.

	
	By (1)-(4) of Definition <ref>, 
	
    (^ V)(t, ω) + C(t, ω_t, t, ω_t ∧·,  (t, ω_t ∧·)) - ∫^T_t (^ c^r)(t, ω, t, ω_t) dr     - (^ f)(t, ω, t, ω_t)  - ^(G ♢ g)(t, ω) = 0.
 
	
	As  is admissible and V satisfies Assumption <ref>, we apply the functional Itô formula in Theorem <ref> to V and then claim that ∂_ω V, σ^t, ·Ŵ is a true martingale, where the Brownian motion Ŵ is a part of the weak solution under . Along with (<ref>), we obtain the following:
	
    [ V(T, X^t, ω, _T ∧·)| _t ] = [ V(T, X^t, ω, ⊗_T Θ^T, )| _t ]     = V(t, X^t, ω, ⊗_t Θ^t, ) + [ ∫^T_t (^ V)(s, X^t, ω, ⊗_s Θ^s, ) ds | _t ]     =  V(t, ω)  - [ ∫^T_t C(s, X^t, ω, ⊗_s Θ^s, , s, X^t, ω, _s ∧·, ( s, X^t, ω, _s ∧·)) ds | _t ] 
          +  [ ∫^T_t ∫^T_s (^ c^r)(s, X^t, ω, ⊗_s Θ^s, , s, X^t, ω, _s ) dr ds | _t ]       +  [ ∫^T_t (^ f)(s, X^t, ω, ⊗_s Θ^s, , s, X^t, ω, _s ) ds | _t ]       +  [ ∫^T_t ^(G ♢ g)(s, X^t, ω, ⊗_s Θ^s, ) ds | _t ].

	
	For the third term, Fubini's theorem holds under the polynomial growth rate condition on the derivatives of c^r by Assumption <ref> and the conditions in (1)-(2) of Definition <ref>. Lemma <ref> shows that ∂_ω c^r, σ^t, ·Ŵ is a true martingale. Hence, the definition of c^r leads to
	
    [ ∫^T_t ∫^T_s (^ c^r)(s, X^t, ω, ⊗_s Θ^s, , s, X^t, ω, _s ) dr ds | _t ] 
        = [ ∫^T_t ∫^r_t (^ c^r)(s, X^t, ω, ⊗_s Θ^s, , s, X^t, ω, _s ) ds dr | _t ] 
        = ∫^T_t [ ∫^r_t (^ c^r)(s, X^t, ω, ⊗_s Θ^s, , s, X^t, ω, _s ) ds  | _t ] dr
        = ∫^T_t {[ c^r(r, X^t, ω, ⊗_r Θ^r, , r, X^t, ω, _r )| _t ] - c^r(t, X^t, ω, ⊗_t Θ^t, , t, X^t, ω, _t ) }dr
       = ∫^T_t {[ C(r, X^t, ω, _r, r,  X^t, ω, _, (r,  X^t, ω, _)) | _t ] 
                -  [ C(t, X^t, ω, _t, r,  X^t, ω, _, (r,  X^t, ω, _)) | _t ] } dr.

	
	For the fourth term, we use the same arguments:
	
    [ ∫^T_t (^ f)(s, X^t, ω, ⊗_s Θ^s, , s, X^t, ω, _s ) ds | _t ] 
        = [ f (T, X^t, ω, ⊗_T Θ^T, , T, X^t, ω, _T )  | _t ] -  f (t, X^t, ω, ⊗_t Θ^t, , t, X^t, ω, _t) 
        = [ F (T, X^t, ω, _T, X^t, ω, _T ∧·)  | _t ] - [ F(t, ω_t, X^t, ω, _T ∧·) | _t].

	
	Similarly, for the fifth term,
	
    [ ∫^T_t ^(G ♢ g)(s, X^t, ω, ⊗_s Θ^s, ) ds | _t ] 
        = [ (G♢ g) (T, X^t, ω, ⊗_T Θ^T, )  | _t ] -  (G♢ g) (t, X^t, ω, ⊗_t Θ^t, ) 
        = [ G(T, X^t, ω, _T, X^t, ω, _T)  | _t ] - G(t, ω_t, [ X^t, ω, _T | _t]).

	
	By the boundary condition in (<ref>), we obtain the following:
	
    V(t, ω) =    [ V(T, X^t, ω, _T ∧·)| _t ] + ∫^T_t [ C(t, X^t, ω, _t, r,  X^t, ω, _, (r,  X^t, ω, _)) | _t ] dr 
        - [ F (T, X^t, ω, _T, X^t, ω, _T ∧·)  | _t ] + [ F(t, ω_t, X^t, ω, _T ∧·) | _t] 
       - [ G (T, X^t, ω, _T, X^t, ω, _T)  | _t ] + G(t, ω_t, [ X^t, ω, _T | _t]) 
    
    	=   ∫^T_t [ C(t, X^t, ω, _t, r,  X^t, ω, _, (r,  X^t, ω, _)) | _t ] dr 
        + [ F(t, ω_t, X^t, ω, _T ∧·) | _t] + G(t, ω_t, [ X^t, ω, _T | _t]) 
    
    	=    J(t, ω; ).

	In other words, we verify that V is the value function with .
	
	Next, we show that  is indeed an equilibrium strategy under Definition <ref>. We apply the recursive relationship in Lemma <ref> with _h. Note that ω^t+h =X^t, ω, _h⊗_t+hΘ^t + h, _h, then
	
    J(t, ω; _h)     =  [ J(t+h, X^t, ω, _h⊗_t+hΘ^t + h, _h; _h) | _t ]
          - {∫^T_t+h[ c^r, _h(t+h, X^t, ω, _h⊗_t+hΘ^t + h, _h, t+h, X^t, ω, _h_t+h) | _t ] dr          - ∫^T_t [ c^r, _h(t+h, X^t, ω, _h⊗_t+hΘ^t + h, _h, t, ω_t)| _t ] dr }      - {[ f^_h(t+h, X^t, ω, _h⊗_t+hΘ^t + h, _h, t+h, X^t, ω, _h_t+h) | _t ]           - [ f^_h(t+h, X^t, ω, _h⊗_t+hΘ^t + h, _h, t, ω_t) | _t ] }      -{[ G(t+h, X^t, ω, _h_t+h, g^_h(t+h, X^t, ω, _h⊗_t+hΘ^t + h, _h)) | _t ]          - G(t, ω_t, [ g^_h(t+h, X^t, ω, _h⊗_t+hΘ^t + h, _h) | _t] ) }.

	As _h = on [t+h, T], conditional on _t+h, (X^t, ω, _h_s)_s ∈ [t+h, T] has the same distribution as (X^t, ω, _s)_s ∈ [t+h, T]. Then,
	
    J(t+h, X^t, ω, _h⊗_t+hΘ^t + h, _h; _h) = V(t+h, X^t, ω, _h⊗_t+hΘ^t + h, _h).

	For t+h ≤ r ≤ T,
	
    c^r, _h(t+h, X^t, ω, _h⊗_t+hΘ^t + h, _h, t+h, X^t, ω, _h_t+h) 
    	=    [ C(t+h, X^t, ω, _h_t+h, r, X^t, ω, _h_, ( r, X^t, ω, _h_)) | _t+h] 
    
    	=     c^r(t+h, X^t, ω, _h⊗_t+hΘ^t + h, _h, t+h, X^t, ω, _h_t+h).

	
	When t ≤ r ≤ t+h,
	
    c^r, _h(t+h, X^t, ω, _h⊗_t+hΘ^t + h, _h, t, ω_t) = C(t, ω_t, r, X^t, ω, _h_, ( r, X^t, ω, _h_)).

	
	When t+h ≤ r ≤ T,
	
    c^r, _h(t+h, X^t, ω, _h⊗_t+hΘ^t + h, _h, t, ω_t)     = [ C(t, ω_t, r, X^t, ω, _h_, ( r, X^t, ω, _h_)) | _t+h]     = c^r(t+h, X^t, ω, _h⊗_t+hΘ^t + h, _h, t, ω_t).

	
	Similarly,
	
    f^_h(t+h, X^t, ω, _h⊗_t+hΘ^t + h, _h, t+h, X^t, ω, _h_t+h) 
        = [ F(t+h, X^t, ω, _h_t+h, X^t, ω, _h_T ∧·) | _t+h]     =  f(t+h, X^t, ω, _h⊗_t+hΘ^t + h, _h, t+h, X^t, ω, _h_t+h),

	
    f^_h(t+h, X^t, ω, _h⊗_t+hΘ^t + h, _h, t, ω_t) =  f(t+h, X^t, ω, _h⊗_t+hΘ^t + h, _h, t, ω_t),

	and
	
    g^_h(t+h, X^t, ω, _h⊗_t+hΘ^t + h, _h)     = [X^t, ω, _h_T | _t+h] 
        =  g(t+h, X^t, ω, _h⊗_t+hΘ^t + h, _h).

	
	Therefore, (<ref>) is reduced to
	
    J(t, ω; _h)     =  [ V(t+h, X^t, ω, _h⊗_t+hΘ^t + h, _h) | _t ] 
          - {∫^T_t+h[ c^r(t+h, X^t, ω, _h⊗_t+hΘ^t + h, _h, t+h, X^t, ω, _h_t+h) | _t ] dr          - ∫^t+h_t [  C(t, ω_t, r, X^t, ω, _h_, ( r, X^t, ω, _h_)) | _t ] dr          - ∫^T_t+h[ c^r(t+h, X^t, ω, _h⊗_t+hΘ^t + h, _h, t, ω_t) | _t ] dr }      - {[ f(t+h, X^t, ω, _h⊗_t+hΘ^t + h, _h, t+h, X^t, ω, _h_t+h) | _t ]           - [  f(t+h, X^t, ω, _h⊗_t+hΘ^t + h, _h, t, ω_t) | _t ] }      -{[ G(t+h, X^t, ω, _h_t+h, g(t+h, X^t, ω, _h⊗_t+hΘ^t + h, _h)) | _t ]           - G(t, ω_t, [ g(t+h, X^t, ω, _h⊗_t+hΘ^t + h, _h) | _t] ) }.

	
	Meanwhile, from the PHJB equation (<ref>) for V, we apply the functional Itô formula and Lemma <ref>. Note the right-continuity in time and continuity in ω assumption from (2)-(3) of Definition <ref>:
	
    [ V(t+h, X^t, ω, _h⊗_t+hΘ^t + h, _h) | _t ] - V(t, ω)     + [ ∫^t+h_t C(s, X^t, ω, _h_s, s, X^t, ω, _h_s ∧·, ( s, X^t, ω, _h_s ∧·) ) ds | _t ]     - [ ∫^t+h_t ∫^T_s (^_h c^r)(s, X^t, ω, _h⊗_s Θ^s, _h, s, X^t, ω, _h_s) dr ds | _t ]     + [ ∫^t+h_t ∫^T_s (^_h c^t, ω_t, r)(s, X^t, ω, _h⊗_s Θ^s, _h) dr ds | _t ] 
        - {[ f(t+h, X^t, ω, _h⊗_t+hΘ^t + h, _h, t+h, X^t, ω, _h_t+h) | _t ] - f(t, ω, t, ω_t) }    + {[ f(t+h, X^t, ω, _h⊗_t+hΘ^t + h, _h, t, ω_t) | _t ] - f(t, ω, t, ω_t) }    - {[ G(t+h, X^t, ω, _h_t+h, g(t+h, X^t, ω, _h⊗_t+hΘ^t + h, _h)) | _t ] - G(t, ω_t, g(t, ω)) }    + { G(t, ω_t, [ g(t+h, X^t, ω, _h⊗_t+hΘ^t + h, _h) | _t ] ) - G(t, ω_t, g(t, ω)) }≤ o(h).

	
	We further simplify the C, c^r, and c^t, ω_t, r terms in (<ref>). By Fubini's theorem, we obtain the following:
	
    [∫^t+h_t ∫^T_s (^_h c^r)(s, X^t, ω, _h⊗_s Θ^s, _h, s, X^t, ω, _h_s) dr ds | _t ]    = [ ∫^T_t ∫^r ∧ (t+h)_t  (^_h c^r)(s, X^t, ω, _h⊗_s Θ^s, _h, s, X^t, ω, _h_s) ds dr | _t ]    = ∫^t+h_t [ ∫^r_t (^_h c^r)(s, X^t, ω, _h⊗_s Θ^s, _h, s, X^t, ω, _h_s) ds | _t ] dr 
          +  ∫^T_t+h[ ∫^t+h_t (^_h c^r)(s, X^t, ω, _h⊗_s Θ^s, _h, s, X^t, ω, _h_s) ds | _t ] dr    =∫^t+h_t {[ c^r(r, X^t, ω, _h⊗_r Θ^r, _h, r, X^t, ω, _h_r) | _t ] - c^r(t, ω, t, ω_t) }dr       + ∫^T_t+h{[ c^r(t+h, X^t, ω, _h⊗_t+hΘ^t+h, _h, t+h, X^t, ω, _h_t+h) | _t ] - c^r(t, ω, t, ω_t) }dr.

	Similarly,
	
    [ ∫^t+h_t ∫^T_s (^_h c^t, ω_t, r)(s, X^t, ω, _h⊗_s Θ^s, _h) dr ds | _t ]    =∫^t+h_t {[ c^r(r, X^t, ω, _h⊗_r Θ^r, _h, t, ω_t) | _t ] - c^r(t, ω, t, ω_t) }dr 
          + ∫^T_t+h{[ c^r(t+h, X^t, ω, _h⊗_t+hΘ^t+h, _h, t, ω_t) | _t ] - c^r(t, ω, t, ω_t) }dr

	and
	
    ∫^t+h_t [ c^r(r, X^t, ω, _h⊗_r Θ^r, _h, r, X^t, ω, _h_r) | _t ]dr     = ∫^t+h_t [ [ C(r, X^t, ω, _h_r, r, X^t, ω, _h_, (r, X^t, ω, _h_)) | _r ] | _t ]dr 
        = ∫^t+h_t [ C(r, X^t, ω, _h_r, r, X^t, ω, _h_, (r, X^t, ω, _h_)) | _t ]dr.

	
    ∫^t+h_t [ c^r(r, X^t, ω, _h⊗_r Θ^r, _h, t, ω_t) | _t ]dr 
        = ∫^t+h_t [ C(t, ω_t, r, X^t, ω, _h_, (r, X^t, ω, _h_)) | _t ]dr.

	
	By the Lebesgue differentiation theorem held under (3)-(4) of Definition <ref>, we obtain the following:
	
    ∫^t+h_t [ C(r, X^t, ω, _h_r, r, X^t, ω, _h_, (r, X^t, ω, _h_)) | _t ]dr     = ∫^t+h_t [ C(t, ω_t, r, X^t, ω, _h_, (r, X^t, ω, _h_)) | _t ]dr + o(h).

	
	Thus, (<ref>) and (<ref>) are reduced to
	
    - [ ∫^t+h_t ∫^T_s (^_h c^r)(s, X^t, ω, _h⊗_s Θ^s, _h, s, X^t, ω, _h_s) dr ds | _t ]     + [ ∫^t+h_t ∫^T_s (^_h c^t, ω_t, r)(s, X^t, ω, _h⊗_s Θ^s, _h) dr ds | _t ] 
    	=     -∫^T_t+h[ c^r(t+h, X^t, ω, _h⊗_t+hΘ^t+h, _h, t+h, X^t, ω, _h_t+h) | _t ]dr 
        + ∫^T_t+h[ c^r(t+h, X^t, ω, _h⊗_t+hΘ^t+h, _h, t, ω_t) | _t ] dr + o(h).

	Using the same argument in (<ref>) yields the following:
	
    [ ∫^t+h_t C(s, X^t, ω, _h_s, s, X^t, ω, _h_s ∧·, ( s, X^t, ω, _h_s ∧·) ) ds | _t ] 
        = [ ∫^t+h_t  C(t, ω_t, r, X^t, ω, _h_, ( r, X^t, ω, _h_)) dr | _t ] + o(h).

	
	Combining (<ref>), (<ref>), and (<ref>) yields
	
    [ V(t+h, X^t, ω, _h⊗_t+hΘ^t + h, _h) | _t ] - V(t, ω)     + [ ∫^t+h_t  C(t, ω_t, r, X^t, ω, _h_, ( r, X^t, ω, _h_)) dr | _t ]    -∫^T_t+h[ c^r(t+h, X^t, ω, _h⊗_t+hΘ^t+h, _h, t+h, X^t, ω, _h_t+h) | _t ]dr     + ∫^T_t+h[ c^r(t+h, X^t, ω, _h⊗_t+hΘ^t+h, _h, t, ω_t) | _t ] dr 
        - [ f(t+h, X^t, ω, _h⊗_t+hΘ^t + h, _h, t+h, X^t, ω, _h_t+h) | _t ]     + [ f(t+h, X^t, ω, _h⊗_t+hΘ^t + h, _h, t, ω_t) | _t ]     - [ G(t+h, X^t, ω, _h_t+h, g(t+h, X^t, ω, _h⊗_t+hΘ^t + h, _h)) | _t ]     + G(t, ω_t, [ g(t+h, X^t, ω, _h⊗_t+hΘ^t + h, _h) | _t ] ) ≤ o(h).

	
	Compared with (<ref>), we ensure that
	
    J(t, ω; _h) - V(t, ω) ≤ o(h).

	As V(t, ω) = J(t, ω; ) is shown previously,
	
    J(t, ω; _h) - J(t, ω; ) ≤ o(h),

	as desired.




  §.§.§ Proof of Corollary <ref>


	We first prove the results for ψ. Consider ψ̃= - ψ. Then, ψ̃ satisfies
	
    ψ̃= K * H(ψ).
 
	w_* is the unique root of H(w) = 0 on (-∞, w_max] with w_max≜ - H_1/2H_2. H(w) satisfies Assumption A.1 in <cit.>. Therefore, <cit.> with a(t) ≡ 0 implies that (<ref>) has a unique global continuous solution and that
	
    w_* < r_1(t) ≤ψ̃(t) < 0,   ∀  t > 0.

	This gives the desired result for ψ.

	(<ref>) is a linear VIE. The existence and uniqueness results of V_2 are thus determined using <cit.> or <cit.>. (<ref>) and (<ref>) are linear ordinary differential equations (ODEs).




  §.§.§ Proof of Corollary <ref>


	M^* under  is given by
	
    M^*_t = M_0 exp[ ∫^t_0 (_s + θπ̂_s ν_s - 1/2π̂^2_s ν_s)ds + ∫^t_0 √(ν_s)π̂_s d Ŵ_1s],

	where Brownian motion Ŵ_1 is a part of the weak solution under .
 	
	By Doob's maximal inequality and <cit.>,
	
    [ sup_ t ∈ [0, T] | M^*_t |^p ] 
       ≤ C [ sup_ t ∈ [0, T]| e^- ∫^t_0 θπ̂_s ν_s ds|^2p] + C [ sup_ t ∈ [0, T]|  exp(- 1/2∫^t_0  π̂^2_s ν_s ds + ∫^t_0 √(ν_s)π̂_s d Ŵ_1s) |^2p] 
       ≤ C [ e^2p ∫^T_0 |θπ̂_s| ν_s ds] + C [ exp(- ∫^T_0 p π̂^2_s  ν_s ds + ∫^T_0 2p π̂_s √(ν_s) d Ŵ_1s) ].

	The first term is finite by assumption (<ref>) with constant 2 p |θ| sup_t ∈ [0, T] | π̂_t |. The second term is also finite. In fact, by Hölder's inequality and assumption (<ref>) with a constant (8p^2 - 2p) sup_t ∈ [0, T]π̂^2_t,
	
    [ exp(- ∫^T_0 p π̂^2_s  ν_s ds + ∫^T_0 2p π̂_s √(ν_s) d Ŵ_1s) ] 
       ≤{[ e^(8p^2 - 2p) ∫^T_0 π̂^2_s  ν_s ds] }^1/2{[ exp(- 8 p^2 ∫^T_0 π̂^2_s  ν_s ds + 4 p ∫^T_0 π̂_s √(ν_s) d Ŵ_1s) ] }^1/2
        < ∞.


















































  §.§.§ Proof of Corollary <ref>


	When λ = 0, the claim is directly verifiable. We suppose that λ≠ 0. By <cit.> or <cit.>,
	
    F^α, λ(T - t)/λT-t →∞∼1/λ - 1/λ^2 Γ(1 - α) (T - t)^α.

	Note that
	
    ∂_α[Γ(1 - α) (T - t)^α] = [ln(T - t) - ψ( 1 - α)] Γ(1 - α) (T - t)^α > 0,

	where ψ(z) = Γ'(z)/Γ(z) is the polygamma function and ψ( 1 - α) < 0. We get the first part of the claim.
	
	Furthermore,
	
    F^α, λ(T - t)/λT-t → 0^+∼(T - t)^α/Γ(α + 1).

	When T - t is small, (T - t)^α is decreasing on α. Finally, (T - t)^α/Γ(α + 1) decreases with α.







	
	[Abi JaberAbi Jaber2019]abi2019lifting
	Abi Jaber, E. (2019).
	Lifting the Heston model.
	Quantitative Finance, 19(12), 1995–2013.
	
	[Abi Jaber and El EuchAbi Jaber and
		El Euch2019]abi2019multifactor
	Abi Jaber, E. and El Euch, O. (2019).
	Multifactor approximation of rough volatility models.
	SIAM Journal on Financial Mathematics, 10(2), 309–349.
	
	[Abi Jaber, Larsson and PulidoAbi Jaber
		et al.2019]abi2017affine
	Abi Jaber, E., Larsson, M., and Pulido, S. (2019).
	Affine Volterra processes.
	The Annals of Applied Probability, 29(5), 3155–3200.
	
	[Alòs, León and VivesAlòs
		et al.2007]alos2007short
	Alòs, E., León, J. A., and Vives, J. (2007).
	On the short-time behavior of the implied volatility for
	jump-diffusion models with stochastic volatility.
	Finance and Stochastics, 11(4), 571–589.
	
	[Basak and ChabakauriBasak and
		Chabakauri2010]basak2010dynamic
	Basak, S. and Chabakauri, G. (2010).
	Dynamic mean-variance asset allocation.
	The Review of Financial Studies, 23(8), 2970–3016.
	
	[Bayer, Friz and GatheralBayer
		et al.2016]bayer2016pricing
	Bayer, C., Friz, P., and Gatheral, J. (2016).
	Pricing under rough volatility.
	Quantitative Finance, 16(6), 887–904.
	
	[Bäuerle and DesmettreBäuerle and
		Desmettre2020]bauerle2020portfolio
	Bäuerle, N. and Desmettre, S. (2020).
	Portfolio optimization in fractional and rough Heston models.
	SIAM Journal on Financial Mathematics, 11(1), 240–273.
	
	[Björk, Khapko and MurgociBjörk
		et al.2017]bjork2017time
	Björk, T., Khapko, M., and Murgoci, A. (2017).
	On time-inconsistent stochastic control in continuous time.
	Finance and Stochastics, 21(2), 331–360.
	
	[Björk and MurgociBjörk and
		Murgoci2014]bjork2014theory
	Björk, T. and Murgoci, A. (2014).
	A theory of Markovian time-inconsistent stochastic control in
	discrete time.
	Finance and Stochastics, 18(3), 545–592.
	
	[Björk, Murgoci and ZhouBjörk
		et al.2014]bjork2014mean
	Björk, T., Murgoci, A., and Zhou, X. Y. (2014).
	Mean-variance portfolio optimization with state-dependent risk
	aversion.
	Mathematical Finance, 24(1), 1–24.
	
	[BrunnerBrunner2017]brunner2017volterra
	Brunner, H. (2017).
	Volterra integral equations: an introduction to theory and
		applications, volume 30.
	Cambridge University Press.
	
	[Cont and FourniéCont and
		Fournié2013]cont2013functional
	Cont, R. and Fournié, D.-A. (2013).
	Functional Itô calculus and stochastic integral representation of
	martingales.
	The Annals of Probability, 41(1), 109–133.
	
	[Cvitanić, Possamaï and
		TouziCvitanić et al.2017]cvitanic2017moral
	Cvitanić, J., Possamaï, D., and Touzi, N. (2017).
	Moral hazard in dynamic risk management.
	Management Science, 63(10), 3328–3346.
	
	[CzichowskyCzichowsky2013]czichowsky2013time
	Czichowsky, C. (2013).
	Time-consistent mean-variance portfolio selection in discrete and
	continuous time.
	Finance and Stochastics, 17(2), 227–271.
	
	[Dai, Jin, Kou and XuDai
		et al.2020]dai2020log
	Dai, M., Jin, H., Kou, S., and Xu, Y. (2020).
	A dynamic mean-variance analysis for log returns.
	Management Science.
	<https://doi.org/10.1287/mnsc.2019.3493>.
	
	[DupireDupire2019]dupire2019functional
	Dupire, B. (2019).
	Functional Itô calculus.
	Quantitative Finance, 19(5), 721–729.
	
	[Ekeland and PirvuEkeland and
		Pirvu2008]ekeland2008investment
	Ekeland, I. and Pirvu, T. A. (2008).
	Investment and consumption without commitment.
	Mathematics and Financial Economics, 2(1), 57–86.
	
	[El Euch and RosenbaumEl Euch and
		Rosenbaum2019]el2019characteristic
	El Euch, O. and Rosenbaum, M. (2019).
	The characteristic function of rough Heston models.
	Mathematical Finance, 29(1), 3–38.
	
	[Fouque and HuFouque and
		Hu2018]fouque2018optimalSIFN
	Fouque, J.-P. and Hu, R. (2018).
	Optimal portfolio under fast mean-reverting fractional stochastic
	environment.
	SIAM Journal on Financial Mathematics, 9(2), 564–601.
	
	[Fouque and HuFouque and
		Hu2019]fouque2019optimalMF
	Fouque, J.-P. and Hu, R. (2019).
	Optimal portfolio under fractional stochastic environment.
	Mathematical Finance, 29(3), 697–734.
	
	[Fukasawa, Takabatake and WestphalFukasawa
		et al.2019]fukasawa2019volatility
	Fukasawa, M., Takabatake, T., and Westphal, R. (2019).
	Is volatility rough?
	arXiv preprint arXiv:1905.04852.
	
	[Gatheral, Jaisson and RosenbaumGatheral
		et al.2018]gatheral2018volatility
	Gatheral, J., Jaisson, T., and Rosenbaum, M. (2018).
	Volatility is rough.
	Quantitative Finance, 18(6), 933–949.
	
	[Gatheral and Keller-ResselGatheral and
		Keller-Ressel2019]gatheral2019affine
	Gatheral, J. and Keller-Ressel, M. (2019).
	Affine forward variance models.
	Finance and Stochastics, 23(3), 501–533.
	
	[Gerhold, Gerstenecker and PinterGerhold
		et al.2019]gerhold2019moment
	Gerhold, S., Gerstenecker, C., and Pinter, A. (2019).
	Moment explosions in the rough Heston model.
	Decisions in Economics and Finance, 42(2), 575–608.
	
	[Glasserman and HeGlasserman and
		He2019]glasserman2019buy
	Glasserman, P. and He, P. (2019).
	Buy rough, sell smooth.
	Quantitative Finance, 1–16.
	
	[Gripenberg, Londen and StaffansGripenberg
		et al.1990]gripenberg1990volterra
	Gripenberg, G., Londen, S.-O., and Staffans, O. (1990).
	Volterra integral and functional equations, volume 34.
	Cambridge University Press.
	
	[Guennoun, Jacquier, Roome and ShiGuennoun
		et al.2018]guennoun2018asymptotic
	Guennoun, H., Jacquier, A., Roome, P., and Shi, F. (2018).
	Asymptotic behavior of the fractional Heston model.
	SIAM Journal on Financial Mathematics, 9(3),
	1017–1045.
	
	[Han and WongHan and Wong2020a]han2019mean
	Han, B. and Wong, H. Y. (2020a).
	Mean-variance portfolio selection under Volterra Heston model.
	Applied Mathematics and Optimization.
	<https://doi.org/10.1007/s00245-020-09658-3>.
	
	[Han and WongHan and
		Wong2020b]han2019merton
	Han, B. and Wong, H. Y. (2020b).
	Merton's portfolio problem under Volterra Heston model.
	Finance Research Letters, 101580.
	
	[He and JiangHe and
		Jiang2019]he2018equilibrium
	He, X. D. and Jiang, Z. (2019).
	Equilibrium strategies in continuous-time time-inconsistent problems.
	Available at SSRN 3308274.
	
	[Hu, Jin and ZhouHu et al.2012]hu2012time
	Hu, Y., Jin, H., and Zhou, X. Y. (2012).
	Time-inconsistent stochastic linear–quadratic control.
	SIAM Journal on Control and Optimization, 50(3),
	1548–1572.
	
	[Huang and ZhouHuang and
		Zhou2018]huang2018strong
	Huang, Y.-J. and Zhou, Z. (2018).
	Strong and weak equilibria for time-inconsistent stochastic control
	in continuous time.
	arXiv preprint arXiv:1809.09243. Forthcoming in Mathematics of
		Operations Research.
	
	[Jusselin and RosenbaumJusselin and
		Rosenbaum2020]jusselin2018no
	Jusselin, P. and Rosenbaum, M. (2020).
	No-arbitrage implies power-law market impact and rough volatility.
	Mathematical Finance, 1–28.
	<https://doi.org/10.1111/mafi.12254>.
	
	[KalininKalinin2019]kalinin2019support
	Kalinin, A. (2019).
	Support characterization for regular path-dependent stochastic Volterra integral equations.
	arXiv preprint arXiv:1908.10786.
	
	[Keller-Ressel, Larsson and
		PulidoKeller-Ressel et al.2018]keller2018affine
	Keller-Ressel, M., Larsson, M., and Pulido, S. (2018).
	Affine rough models.
	arXiv preprint arXiv:1812.08486.
	
	[Li and NgLi and Ng2000]li2000optimal
	Li, D. and Ng, W.-L. (2000).
	Optimal dynamic portfolio selection: Multiperiod mean-variance
	formulation.
	Mathematical Finance, 10(3), 387–406.
	
	[LiuLiu2007]liu2007portfolio
	Liu, J. (2007).
	Portfolio selection in stochastic environments.
	The Review of Financial Studies, 20(1), 1–39.
	
	[MainardiMainardi2014]mainardi2014some
	Mainardi, F. (2014).
	On some properties of the Mittag-Leffler function E_α
		(-t^α), completely monotone for t> 0 with 0<α< 1.
	Discrete and Continuous Dynamical Systems-B, 19(7),
	2267–2278.
	
	[MarkowitzMarkowitz1952]markowitz1952portfolio
	Markowitz, H. (1952).
	Portfolio selection.
	The Journal of Finance, 7(1), 77–91.
	
	[MertonMerton1969]merton1969lifetime
	Merton, R. C. (1969).
	Lifetime portfolio selection under uncertainty: The continuous-time
	case.
	The Review of Economics and Statistics, 247–257.
	
	[Mytnik and SalisburyMytnik and
		Salisbury2015]mytnik2015uniqueness
	Mytnik, L. and Salisbury, T. S. (2015).
	Uniqueness for Volterra-type stochastic integral equations.
	arXiv preprint arXiv:1502.05513.
	
	[SaporitoSaporito2019]saporito2019stochastic
	Saporito, Y. F. (2019).
	Stochastic control and differential games with path-dependent
	influence of controls on dynamics and running cost.
	SIAM Journal on Control and Optimization, 57(2),
	1312–1327.
	
	[Schied, Speiser and VoloshchenkoSchied
		et al.2018]schied2018model
	Schied, A., Speiser, L., and Voloshchenko, I. (2018).
	Model-free portfolio theory and its functional master formula.
	SIAM Journal on Financial Mathematics, 9(3),
	1074–1101.
	
	[SiuSiu2016]siu2016functional
	Siu, T. K. (2016).
	A functional Itô's calculus approach to convex risk measures with
	jump diffusion.
	European Journal of Operational Research, 250(3),
	874–883.
	
	[StrotzStrotz1955]strotz1955myopia
	Strotz, R. H. (1955).
	Myopia and inconsistency in dynamic utility maximization.
	The Review of Economic Studies, 23(3), 165–180.
	
	[Viens and ZhangViens and
		Zhang2019]viens2017martingale
	Viens, F. and Zhang, J. (2019).
	A martingale approach for fractional Brownian motions and related
	path dependent PDEs.
	The Annals of Applied Probability, 29(6), 3489–3540.
	
	[YongYong2012]yong2012time
	Yong, J. (2012).
	Time-inconsistent optimal control problems and the equilibrium HJB
	equation.
	Mathematical Control and Related Fields, 2(3), 271–329.
	
	[Zhao, Shen and WeiZhao
		et al.2014]zhao2014consumption
	Zhao, Q., Shen, Y., and Wei, J. (2014).
	Consumption–investment strategies with non-exponential discounting
	and logarithmic utility.
	European Journal of Operational Research, 238(3),
	824–835.
	
	[Zhou and LiZhou and
		Li2000]zhou2000continuous
	Zhou, X. Y. and Li, D. (2000).
	Continuous-time mean-variance portfolio selection: A stochastic
	LQ framework.
	Applied Mathematics and Optimization, 42(1), 19–33.
	

