REVIEW 5 major objections 4 minor 23 references
Portfolio Optimization Managing Value at Risk under Heavy Tail Return, using Stochastic Maximum Principle
T0 review · 5 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read This paper claims to provide a closed-form approximate optimal portfolio strategy for heavy-tailed stock returns by maximizing median wealth subject to a Value-at-Risk constraint, and shows empirically that higher risk-free rates reduce…
desk verdict Good idea for continuous-time quantile-based portfolio choice, but the central SDE is misspecified and the median/VaR labels are swapped, so the derived optimal strategy solves a different problem. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the asymptotic covariance matrix $\Lambda$ of two sample quantiles, with entries $\Lambda_{ij}=p_i(1-p_j)/(f_if_j)$ for $i,j\in\{1,2\}$, where $p_1=0.05$, $p_2=0.5$, and $f_1,f_2$ are the densities at those quantiles. The paper uses $\Lambda$ as the diffusion coefficient in the SDE governing the quantile processes, turning a statistical fact about sample estimates into a stochastic control model. The second piece of machinery is the sigmoid approximation of the indicator function in the VaR constraint, which makes the objective smooth; the stochastic maximum principle then supplies adjoint equations whose reduced solution leads to the first-order condition and, after a further approximation, to formula (7).
What would settle it
Simulate a heavy-tailed return process with known parameters, form rolling 5% and 50% sample quantiles over time, and compare the empirical quadratic covariation of those quantile paths with the value implied by $\Lambda$; a mismatch means the quantile-SDE assumption is false, so formula (7) would not be the optimum for the simulated process.
Extended reading notes
Core claim
On its own terms, the paper's result is a closed-form approximate solution to a constrained stochastic control problem, whose state is the pair of wealth processes associated with the two quantiles: the median process $\bar L_1$ carries the utility, and the 5% quantile process $\bar L_2$ carries the Value-at-Risk constraint. After forming the Hamiltonian, solving the two adjoint equations with zero matrix-valued terms, and setting the derivative with respect to the portfolio weight to zero, the paper obtains the optimal portfolio weight $$\bar\pi_t=\frac{$e^{{-\beta t}}$\psi_0 \bar L_1(t)^\gamma $r^{2}$-\psi k_1\bar L_2(t)(b_2-r)(r-\$\beta$)}{\left($e^{{-\beta t}}$\psi_0 \bar L_1(t)^\gamma-\psi k_1\bar L_2(t)\right)(b_2-r)r}.$$ The paper presents this as the optimal strategy for the approximated problem and, via Ekeland's theorem, as an approximate optimum of the original non-smooth problem; it then supports the formula with simulations showing that higher risk-free rates lower the optimal stock fraction and that changes in $\gamma$ and $\beta$ leave the strategy largely unchanged.
Load-bearing premise
The load-bearing assumption is that the asymptotic covariance matrix $\Lambda$ of quantile estimates, a statement about how estimates vary across samples, can be used directly as the instantaneous random-shock matrix of the quantile processes' own time-series SDE; if the quantile paths do not follow that SDE, formula (7) has not been derived for the actual heavy-tailed return process.
Editorial extensions
If this is right
- If formula (7) is the optimal strategy, a higher risk-free rate $r$ reduces the fraction held in the risky asset, matching the paper's empirical finding.
- The strategy is implementable without specifying the return distribution, because calibration needs only sample quantiles and kernel-density estimates of the quantile densities.
- Putting more weight on the utility objective and less on the VaR constraint increases both the stock allocation and the final wealth, consistent with accepting more tail risk for higher median gains.
- The framework extends to logarithmic and exponential utility, and other constant-risk-aversion forms, at the cost of algebraic complexity, so the median/VaR formulation is not tied to power utility.
- Because the constraint is smoothed, the resulting strategy is an approximate optimum rather than an exact one; Ekeland's theorem is invoked to justify the approximation.
Reading between the lines
- A testable extension is to replace the 5% quantile process by an expected-shortfall or conditional-VaR functional in the same Hamiltonian machinery, since the calibration only requires quantile asymptotics and a smooth approximation of the constraint.
- A cautious reading treats formula (7) as a calibrated heuristic rather than an exact optimum for the original heavy-tailed process, because the derivation replaces the return process with a quantile diffusion and then smooths the constraint; simulated fat-tailed data would show how large the gap is.
- The empirical insensitivity to $\gamma$ and $\beta$ suggests the allocation is governed mainly by the ratio of median to tail quantile levels and by the risk-free rate; if confirmed, the rule could be condensed into a simple tail-thickness-dependent trading formula.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper considers a continuous-time portfolio problem with one risky asset and one risk-free asset when the risky return is heavy-tailed. Since moments may not exist, the authors propose to work with two return quantiles (the 5% and the 50% quantiles), use a multivariate normal approximation for their joint distribution, and model the quantile dynamics by an Itô SDE. They then formulate a constrained stochastic control problem that maximizes discounted power utility of a median-wealth process subject to a Value-at-Risk constraint, solve it by the stochastic maximum principle, and derive a closed-form approximate optimal strategy in equation (7). The model is calibrated nonparametrically with a kernel density estimator on daily Entergy Corporation returns, and the resulting strategy and wealth paths are simulated for various parameter choices. The paper claims that the results accord with financial intuition.
Significance. If the derivation were correct, the paper would offer a genuinely useful construction: a quantile-based, nonparametric dynamic allocation rule that bypasses nonexistent moments under heavy tails. The choice of quantiles and the use of asymptotic normality of sample quantiles are appropriate building blocks, and the empirical setting is relevant. However, the central analytical chain contains load-bearing errors: the quantile SDE is misspecified, the objective and constraint labels are reversed, the adjoint BSDEs are solved by setting the martingale integrands to zero without justification, and the closed-form formula (7) rests on an unvalidated exponential approximation. These are not local presentation issues; they affect the claimed optimal strategy and therefore the entire contribution. The paper provides no machine-checked proofs, no code, and no out-of-sample validation, so the empirical conclusions are also not independently verifiable.
major comments (5)
- [Section 2, quantile SDE after the Proposition] The diffusion matrix in the SDE for (dX1, dX2) is set equal to the asymptotic covariance matrix Lambda, but with independent Brownian motions the instantaneous covariance of the increments is Lambda Lambda^T, not Lambda. For example, Var(dX1) would be (p1(1-p1)/f1^2)^2 + (p1(1-p2)/(f1 f2))^2, not p1(1-p1)/f1^2. A correct diffusion model requires a matrix square root A with A A^T = Lambda. Moreover, the displayed Lambda also omits the 1/N factor from the standard asymptotic covariance of sample quantiles. Since equations (1), (2), the Hamiltonian, the adjoint equations, and formula (7) are all derived from this misspecified SDE, the closed-form strategy is not the optimal strategy for the process the authors intend to model.
- [Section 2, wealth processes L(1) and L(2)] The paper says L(1) is the median wealth, but its dynamics are driven by dX1 = dX(p1) with p1 = 0.05, the lower 5% quantile; it then says L(2) is the lower-quantile (VaR) process, but its dynamics are driven by dX2 = dX(p2) with p2 = 0.5, the median. Consequently the objective in equation (2) maximizes a process driven by the 5% return quantile, while the constraint is applied to a process driven by the median. This reverses the roles of the objective and the risk constraint, so the optimization problem actually solved is not the stated one.
- [Section 2, adjoint equations (3) and (4)] The adjoint BSDEs are not solved; the martingale integrands q(t) and Q(t) are simply set to zero. In the stochastic maximum principle, q and Q are part of the solution and are determined by the BSDE, and q appears in the Hamiltonian through a trace term involving sigma. Because sigma depends on the control pi, dropping q changes the first-order condition used to derive the strategy. No argument is given that q = 0 and Q = 0 are the actual integrands for this problem, so the resulting s(t), S(t), and ultimately equation (7) are not established as the SMP solution.
- [Section 3, derivation of equation (7) from equation (6)] The step from equation (6) to equation (7) replaces the ratio of exponential terms by 1 without any error bound or numerical justification. This approximation is the step that produces the claimed closed-form optimal strategy. Since the resulting formula is the central output of the paper, the approximation must be validated with explicit error estimates or at least a careful numerical check over the relevant parameter range; as written, equation (7) is at best a heuristic.
- [Section 3, calibration and empirical illustration] The parameters b2, Q0.05, f0.05, f0.5, and the KDE bandwidth are calibrated on the same Entergy Corporation series that is then used to generate the optimal-strategy and wealth curves. The reported regularities, such as higher risk-free rates reducing stock investment, are therefore in-sample properties of the fitted formula rather than out-of-sample evidence. A holdout period or a proper cross-validation scheme is needed before the empirical claims can support the analytical conclusions.
minor comments (4)
- [Section 3 and List of Figures] The figures are not actually included; the text and the figure list contain placeholders such as 'Figure 1 and Figure 2 should be placed here.' This makes it impossible to verify the quantitative claims about confidence intervals and wealth paths.
- [Section 2, utility discussion] The text refers to negative exponential utility as 'Constant Absolute Relative Risk Aversion (CARA)' and says RRA = eta x; this appears to conflate CARA and CRRA. Negative exponential utility has constant absolute risk aversion, not constant relative risk aversion.
- [Section 2, Remark after equation (6)] The Remark invokes 'Ekeland's theorem' but cites reference [6], which is Clarke's book, and it does not state the metric or the sense in which the solution is approximate. As written, the remark cannot justify the combined effect of the indicator smoothing, the q = 0 and Q = 0 simplifications, and the exponential approximation.
- [Section 3, simulated wealth magnitudes] For r = 0.0004 the text reports portfolio wealth reaching an order of 10^15 starting from L1 = 1 over a 795-day horizon; such magnitudes suggest possible numerical instability in the discretization of equation (7) and should be diagnosed rather than reported as a financial outcome.
Circularity Check
No circularity: the derivation is self-contained; suspected issues are modeling/calibration concerns, not equivalence by construction.
full rationale
The paper's claimed derivation does not reduce to its inputs. The asymptotic joint normality of sample quantiles is quoted from Beach and Davidson [4], an external source, and the quantile SDE is then written as a modeling step rather than as a consequence of the optimality conditions. The maximum-principle framework is imported from Yong and Zhou [23] and applied to the stated Hamiltonian; the optimality condition (6) and the closed-form approximation (7) are obtained by solving the resulting first-order conditions, not by assuming the conclusion. The only self-citation, Biswas and Mukherjee [5], appears in the literature review as discrete-time background and is not load-bearing for any theorem or formula in the paper. The numerical section calibrates Q0.05, f0.05, f0.5, b2, and bandwidth on the Entergy series, but it presents this as an empirical illustration rather than as an out-of-sample prediction; the reported comparative static that higher r lowers stock investment follows from evaluating formula (7) at different values of r, so it is not a fitted parameter renamed as a prediction. Concerns such as using Lambda directly as a diffusion matrix and the labeling of L(1) and L(2) relative to the 5% and 50% quantiles are correctness/modeling issues, not circularity: no equation in the chain is equivalent to its own input by construction.
Assumptions & free parameters
free parameters (9)
- b2 =
0.00599
- Q0.05 =
0.00077
- f0.05 =
47.63579
- f0.5 =
68.43975
- alpha =
10
- epsilon =
0.00001
- KDE bandwidth =
0.00271447
- psi, psi0 =
(0.6,0.8), (0.8,0.6), (0.95,0.312), (0.312,0.95)
- rho =
0.22941
assumptions (8)
- domain assumption The return distribution F is continuous, differentiable to second order, and strictly monotonic at the quantiles of interest.
- standard math Sample quantiles from the return series are asymptotically jointly normal with covariance matrix Lambda from Beach and Davidson.
- ad hoc to paper The asymptotic quantile covariance matrix Lambda can serve directly as the diffusion coefficient in a continuous-time SDE for the quantile processes.
- domain assumption The market is frictionless and the portfolio is self-financing.
- standard math Yong and Zhou's stochastic maximum principle applies to the sigmoid-approximated constrained problem and the adjoint BSDEs have unique solutions.
- ad hoc to paper The martingale integrand q(t) and Q(t) in the adjoint BSDEs can be set to zero.
- ad hoc to paper The ratio of exponential terms in equation (6) is approximately 1 for all t in the horizon.
- ad hoc to paper The indicator function in the VaR constraint can be replaced by a sigmoid with parameters alpha=10 and epsilon=0.00001 without materially changing the optimum.
Cite this review
Pith. "Pith review of Portfolio Optimization Managing Value at Risk under Heavy Tail Return, using Stochastic Maximum Principle." pith.science (2026). https://pith.science/paper/DRYYFVO2
@misc{pith2026190803905,
author = {Pith},
title = {Pith review of: Portfolio Optimization Managing Value at Risk under Heavy Tail Return, using Stochastic Maximum Principle},
year = {2026},
howpublished = {\url{https://pith.science/paper/DRYYFVO2}},
note = {Machine review of arXiv:1908.03905}
}
read the original abstract
We consider an investor, whose portfolio consists of a single risky asset and a risk free asset, who wants to maximize his expected utility of the portfolio subject to managing the Value at Risk (VaR) assuming a heavy tailed distribution of the stock prices return. We use a stochastic maximum principle to formulate the dynamic optimisation problem. The equations which we obtain does not have any explicit analytical solution, so we look for accurate approximations to estimate the value function and optimal strategy. As our calibration strategy is non-parametric in nature, no prior knowledge on the form of the distribution function is needed. We also provide detailed empirical illustration using real life data. Our results show close concordance with financial intuition.We expect that our results will add to the arsenal of the high frequency traders.
Figures
Figures from the paper (5 more)
Reference graph
Works this paper leans on
-
[1]
A. Agarwal & R. Sircar, Portfolio Benchmarking under Drawdown Constraint and Stochastic Sharpe Ratio, SIAM J. Financial Math, 9-2, 435–464 (2018)
work page 2018
-
[2]
Y . Aït-Sahalia & A. W. Lo, Nonparametric Estimation of State-Price Densities Implicit in Financial Asset Prices, The Journal of Finance, 53 (2), 499 - 547 (1998)
work page 1998
-
[3]
T. W. Archibald & E. Possani, Investment and operational decisions for start-up companies: a game theory and Markov decision process approach, Annals of Operations Research, https://doi.org/10.1007/s10479-019- 03426-5 (2019)
-
[4]
C. M. Beach & R. Davidson, Distribution -Free Statistical Inference with Lorenz Curves and Income Shares, Review of Economic Studies, L 723-735 (1983)
work page 1983
-
[5]
S. Biswas & D. Mukherjee, Discrete Time Portfolio Optimization managing Value at Risk under heavy tail return distribution, International Journal of Mathematical Modelling and Numerical Optimization, Accepted (2020)
work page 2020
-
[6]
F. H. Clarke, Optimization and Nonsmooth Analysis , Classics in Applied Mathematics, SIAM, Philadelphia, 1990
work page 1990
- [7]
-
[8]
G. B. Folland, Real analysis, Pure and Applied Mathematics, John Wiley & Sons Inc., New york (1999)
work page 1999
Show all 23 references
-
[9]
Fu, Variance Optimization for Continuous-Time Markov Decision Processes, Scientific Research Open Access, V olume 09(02),10.4236/ojs.2019.92014 (2019)
Y . Fu, Variance Optimization for Continuous-Time Markov Decision Processes, Scientific Research Open Access, V olume 09(02),10.4236/ojs.2019.92014 (2019)
2019
-
[10]
S. J. Grossman & Z. Zhou, Optimal investment strategies for controlling drawdowns, Math. Finance 3(3), 241–276 (1993)
1993
-
[11]
F. C. Harmantzis, L. Miao & Y . Chien, Empirical study of value-at-risk and expected shortfall models with heavy tails, The Journal of Risk Finance, V ol. 7 Issue: 2, pp.117-135 (2006) 22
2006
-
[12]
Y . S. Kim, R. Giacometti, S. T. Rachev, Frank J. Fabozzi & Domenico Mignacca, Measuring financial risk and portfolio optimization with a non-Gaussian multivariate model, Annals of Operations Research, V olume 201, Issue 1, pp 325–343 (2012)
2012
-
[13]
Y . Li, A. Li & Z. Liu, Two Ways of Calculating VaR in Risk Management — An Empirical Study Based on CSI300 Index, Procedia Computer Science, V olume 139, Pages 432-439 (2018)
2018
-
[14]
Markovich, Non parametric analysis of univariate heavy-tailed distribution, Wiley (2007)
N. Markovich, Non parametric analysis of univariate heavy-tailed distribution, Wiley (2007)
2007
-
[15]
M. Á. M. Mato, Classic and modern measures of risk in fixed-income portfolio optimization, The Journal of Risk Finance, V ol. 6 Issue: 5, 416-423 (2005)
2005
-
[16]
https://www.quantopian.com, Accessed: 2018-10-10
2018
-
[17]
Perez, D
I. Perez, D. Hodge & H. Le, Markov decision process algorithms for wealth allocation problems with de- faultable bonds, Advances in Applied Probability, V olume 48(02), pp. 392-405 (2016)
2016
-
[18]
M. L. Puterman, Markov Decision Processes: Discrete stochastic dynamic programming, 79-82, John Wiley New Jersey (2005)
2005
-
[19]
Rao, Linear Statistical Inference and Its Applications, John Wiley & Sons, New York (1965)
C.R. Rao, Linear Statistical Inference and Its Applications, John Wiley & Sons, New York (1965)
1965
-
[20]
D. E. Régisa & R. Artesa, Using multi-state markov models to identify credit card risk, SciELO, V olume 26(2), http://dx.doi.org/10.1590/0103-6513.160814 (2016)
2016
-
[21]
P. A. Samuelson, Lifetime portfolio selection by dynamic stochastic programming, Stochastic Optimization Models in Finance, Pages 517-524 (1975)
1975
-
[22]
L. Tang, An actor-critic-based portfolio investment method inspired by benefit-risk optimization, Journal of Algorithms and Computational Technology, https://doi.org/10.1177/1748301818779059 (2018)
2018 doi
-
[23]
Yong & X
J. Yong & X. Y . Zhou, Stochastic Controls Hamiltonian Systems and HJB Equations, Springer, New York (1999) 23 List of Figures 1 Average optimal strategy and its corresponding 95% interval for different rates of interest (r). . . . 25 2 Portfolio wealth and its corresponding 9...
1999
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.