\section{Introduction}

One of the tacit rules adopted in the modern portfolio optimization models is the expected utility hypothesis on the psychology of decision-making. It dates back to the pioneering work of Merton \cite{M69,M75} in which an investor chooses a portfolio built upon a simplified market and seeks to maximize the expected utility from consumption. For models that do not consider consumption behavior, investors' targets are usually to maximize the expected utility of terminal wealth or/plus either exogenous or endogenous penalty terms. Such goals seem quite natural since the majority of investors are concerned about how much wealth they will accumulate when they decide to liquidate their portfolios. However, in practice, practitioners and institutions, especially those with long-term investment plans, do not always keep their eyes on the capital values of their portfolios. Instead, they care more about their portfolios' growth or return rates (such as seasonal or annual return, internal rate of return), and often adopt relative return rates (such as excess return rate over a stock index or interest rate) to measure their portfolios' performance.

Taking growth rate rather than capital value as a performance evaluation metric has many advantages in practice. It allows for a comparison between various portfolios because it removes the impact of the initial endowment and investment horizon. For fund managers, achieving a higher return rate represents excellent investment skills and brings fame and fortune.
Moreover, it is difficult or almost impossible to do comparisons on portfolios' capital values between worldwide practitioners and institutions who trade different currencies, because exchange rates have so significant impacts that they may change the comparison results. By contrast, taking growth rates (such as the well-known gross domestic products of different regions) as measurement objectives will avoid the setback of using capital values.
In this light, it seems more natural for practitioners to measure their portfolios' performances in terms of growth rates rather than capital values.

Taking growth rate as a measurement objective is not new in the financial economics literature. It at least goes back to the seminal Markowitz's \cite{M52} single-period mean-variance model.
In Markowitz's original formulation, the mean-variance analysis was imposed on the return rate right from the beginning. Ever since most models on mean-variance analysis focus on capital values since they are equivalent in the single-period framework.
Recently, Dai et al. \cite{DJKX20} investigated a dynamic mean-variance portfolio choice problem on growth rate. The optimal policies under specific settings are found to be consistent with several conventional investment wisdoms that are usually contradicted by models based on capital values.

Kelly \cite{K56} introduced the growth-optimal portfolio (GOP) that is to maximize the expected growth rate. In theory and practice, the GOP has enjoyed wide applications in a large number of areas including portfolio theory, utility theory, game theory, information theory, asset pricing theory, and insurance theory (e.g. \cite{ABHSV00,T00}). However, these models ignore the relevant risk when selecting a desirable portfolio. This research intends to partly fill the gap in a continuous-time setting by investigating a GOP model with behavioral risk control.

Another motivation for our paper is to take the behavioral elements into our model.
It is well known that the classical expected utility theory (EUT) and Markowitz's mean-variance theory (MVT) have failed to explain many phenomena and paradoxes (such as the equity premium puzzle \cite{CC99}) in practice. The Nobel prize-winning work Kahneman and Tversky's \cite{KT79, TK92} cumulative prospect theory (CPT) has been widely regarded as a competing alternative to EUT and MVT. There are three key ingredients in the CPT: a reference point, an $S$-shaped utility\footnote{It is convex below the reference point and concave above it.} and a probability weighting (or distortion) function. In our model, the risk is measured by the CPT preference, and the measurement objective is the relative growth rate to some given random or deterministic benchmark rate. We use the well-known Kahneman and Tversky's piece-wise power utility and general smooth probability weighting function. Incorporating behavioral criteria into continuous-time portfolio choice problems has received much attention in recent years (see, e.g., Berkelaar et al. \cite{BKP04}, Jin and Zhou \cite{JZ08}, He and Zhou \cite{HZ11}, Zhang et al. \cite{ZJZ11}, Xu \cite{X16}, Hou and Xu \cite{HX16}, Mi and Xu \cite{MX21}). These works focus on capital value maximization.
To our best knowledge, less attentions are paid to growth rate oriented portfolio choice problems under behavioral criteria. A recent progress comes from Wei and Xu \cite{WX21} which focuses on the mean-risk problem.

Our model leads to a non-concave utility maximization problem under nonlinear Choquet expectation. Due to the presence of probability distortion, the classical stochastic control methods are no more applicable. We adopt the martingale method to transform the problem into a quantile optimization problem, which turns out to be a new type of $M$-shaped\footnote{It is first concave, then convex, and finally concave from left to the right on the real line.} utility maximization problem. We then use the concavification and quantile optimization methods to tackle it, and
derive the optimal growth rate in closed form.
We find that the benchmark growth rate has a significant impact on investment behaviors.

We also present a numerical example to compare our results with Zhang et al.'s \cite{ZJZ11}.
Our terminal relative growth rate is higher than or at least equal to Zhang et al.'s in the intermediate and bad market states. By contrast, ours is lower than theirs in the good market states. From an economic perspective, our strategy is more suitable for loss-sensitive investors than Zhang et al.'s, whereas their strategy is more suitable for those betting for extremely good market states.

The rest of this paper is organized as follows. We present our financial market and model in Section 2 and then turn it into a static random variable optimization problem by the martingale method in Section 3. The latter is transformed into a quantile optimization problem in Section 4 and solved in Section 5. Section 6 presents a numerical example to compare our results to Zhang et al.'s \cite{ZJZ11}. Some proofs and additional analyses are placed in the Appendices.

\section{Problem formulation}\label{Problem formulation}
Let $T>0$ denote a fixed finite investment horizon throughout this paper. Consider a complete filtered probability space $(\Omega, \mathcal{F}, \ensuremath{\opfont{P}})$, on which a standard $n$-dimensional Brownian motion $W(\cdot) \equiv (W^1(\cdot),\ldots, W^n(\cdot))'$ is defined. We assume the uncertainty of the market entirely comes from the Brownian motion and let the information filtration $\{\ensuremath{\mathcal{F}}_t\}_{0\leq t\leq T}$ be generated by the Brownian motion, augmented by all the $\ensuremath{\opfont{P}}$-null sets. Let $\mathcal{F}={\ensuremath{\mathcal{F}}}_T$. For any vector or matrix $M$, we use $M'$ to denote its transpose and $\Vert M\Vert=\sqrt{\text{trace}(MM')}$ to denote the Euclidean norm.
All random variables (r.v.s) considered in this paper are assumed to be $\ensuremath{\mathcal{F}}_T$-measurable.

\subsection{Market and portfolio}

The financial market consists of $n+1$ assets which are traded continuously over the investment horizon $[0, T]$. The market is frictionless, i.e., there are no trading constraints, transaction costs or taxes etc. The first asset is a bond whose price $S_0(\cdot)$ evolves according to an ordinary differential equation (ODE):
\begin{align*}
\begin{cases}
\ensuremath{\operatorname{d}\! } S_0(t) = r(t)S_0(t)\ensuremath{\operatorname{d}\! t}, \quad t\in[0,T], \\
S_0(0) = s_0 > 0,
\end{cases}
\end{align*}
where $r(t)$ is the instantaneous interest rate of the bond at time $t$. The remaining $n$ assets are {stocks} (also called risky assets), and their prices $S_{i}(\cdot)$, $i=1,2,\ldots,n$, are modeled by stochastic differential equations (SDEs):
\begin{align}\label{stockprice}
\begin{cases}
\ensuremath{\operatorname{d}\! } S_i(t) = S_i(t)\big\{\mu_i(t)\ensuremath{\operatorname{d}\! t} + \sum_{j=1}^n\sigma_{ij}(t)\ensuremath{\operatorname{d}\! } W^j(t)\big\}, \quad t\in[0,T], \\
S_i(0) = s_i > 0,
\end{cases}
\end{align}
where $\mu_i : [0,T]\times\Omega\rightarrow \mathbb{R}$ is called the \emph{appreciation rate} of stock $i$, and $\sigma_{ij}: [0,T]\times\Omega\rightarrow \mathbb{R}$ is called the \emph{volatility coefficient} of stock $i$ with respect to $W^j$ at time $t$. Define the appreciation rate vector process $\mu=(\mu_1, \mu_2,\cdots, \mu_n)':[0,T]\times\Omega\rightarrow \mathbb{R}^n$ and
the volatility matrix process $\sigma:=(\sigma_{ij})_{n\times n}:[0,T]\times\Omega\rightarrow \mathbb{R}^{n\times n}$.
We impose the following standard assumption which ensures the market is free of arbitrage opportunities.
\begin{assmp}\label{ASS2.1}
The processes $r$, $\mu$, $\sigma$ are $\{\mathcal{F}_t\}$-progressively measurable and satisfy
\[\int_{0}^{T}|r(s)| d s<+\infty, \text { a.s. }\] and \[\int_{0}^{T}\left[\Vert \mu(t)\Vert+\Vert\sigma(t)\Vert^2\right] \ensuremath{\operatorname{d}\! } t<+\infty, \text { a.s. }\]
Moreover, there exists a unique $\{\mathcal{F}_t\}$-progressively measurable, uniformly bounded risk premium process $\theta:[0,T]\times\Omega\rightarrow \mathbb{R}^{m}$ such that \[\sigma(t)\theta(t)= (\mu_{1}(t)-r(t),\ldots, \mu_n(t)-r(t))^{\prime} \; a.s.,\; a.e. \; t\in[0,T].\]
\end{assmp}
The last statement in \citeassmp{ASS2.1} implies that the matrix $\sigma(t)$ is invertible a.s. for a.e. $t\in[0,T]$.

Consider an investor (``She'') with an initial endowment $x_0>0$. She is a small investor so her actions cannot affect asset prices in the market. Let $\pi_i(t)$ denote the \emph{proportion} of her total wealth invested in stock $i$ at time $t$, $ i = 1,\ldots, n$. The proportion invested in the bond at time $t$ is $1-\sum_{i=1}^{n}\pi_i(t)$.
We call the vector process $\pi:=(\pi_1,\ldots,\pi_n)^{\prime}:[0,T]\times\Omega\rightarrow \mathbb{R}^n$ a \emph{portfolio} (also called a \emph{strategy}) and denote $X^{\pi}$ the agent's \emph{wealth process} with the implementation of the portfolio $\pi$. Assume that the trading of assets takes place continuously in a self-financing fashion (i.e., there is no consumption or income). Then the wealth process $X^{\pi}$ evolves according to the following SDE (see Karatazas and Shreve \cite{KS98}):
\begin{align}\label{wealthprocess}
\begin{cases}
\ensuremath{\operatorname{d}\! } X^{\pi}(t) =X^{\pi}(t) \big[(r(t)+ \pi(t)^{\prime}\sigma(t)\theta(t))\ensuremath{\operatorname{d}\! t} + \pi(t)^{\prime}\sigma(t) \ensuremath{\operatorname{d}\! } W(t)\big],\quad t\geqslant 0, \\
X^{\pi}(0) =x_0>0.
\end{cases}
\end{align}
A portfolio $\pi$ is called admissible if it is an $\{\mathcal{F}_t\}$-progressively measurable process such that
\[\int_{0}^{T} \Vert\sigma(s)^{\prime}\pi(s)\Vert^2 \ensuremath{\operatorname{d}\! s}<\infty, \text { a.s., }\]
and \eqref{wealthprocess} admits a unique strong solution $X^{\pi}$.
From now on, we focus on admissible portfolios, the set of which is denoted by $\mathscr{A}$.

Applying It\^{o}'s lemma to $\log X^{\pi}(t)$, we obtain
\begin{align*}
\qquad X^{\pi}(t)= x_0 \exp\bigg\{\int_{0}^{t}\big(r(s) + \pi(s)^{\prime}\sigma(s)\theta(s)
&-\frac{1}{2}\Vert\sigma(s)^{\prime}\pi(s)\Vert^2 \big)\ensuremath{\operatorname{d}\! s}\\
&+\int_{0}^{t}\pi(s)^{\prime}\sigma(s) \ensuremath{\operatorname{d}\! } W(s)\bigg\}>0. \qquad
\end{align*}
Therefore, we have a no-bankruptcy condition inherent in the wealth process.

\par
A particular feature of our model is that the performance is measured on the \emph{terminal relative growth rate} (also called \emph{terminal relative return rate}) rather than on the terminal wealth as in most financial economics literature. The investor's objective is to minimize the behavioral criterion of the terminal relative growth rate.
To this end, we need to introduce a \emph{growth rate process} and a \emph{benchmark growth rate} (also called the \emph{target growth rate} or simply the \emph{benchmark}).

The growth rate process $R^{\pi}$ associated with a portfolio $\pi$ is defined as
$$R^{\pi}(t)=\log (X^{\pi}(t)/x_{0}),\quad t\geq 0.$$
The normalized growth rate over $[0,t]$ is then given by $t^{-1}R^{\pi}(t)$ for $t>0$.
By It\^{o}'s Lemma and \eqref{wealthprocess}, we have
\begin{align}\label{returnprocess}
\begin{cases}
\ensuremath{\operatorname{d}\! } R^{\pi}(t)= \big(r(t)+ \pi(t)^{\prime}\sigma(t)\theta(t)-\frac{1}{2}\pi(t)^{\prime}\sigma(t) \sigma(t)^{\prime} \pi(t)\big)\ensuremath{\operatorname{d}\! t} + \pi(t)^{\prime}\sigma(t) \ensuremath{\operatorname{d}\! } W(t),\quad t\geqslant 0,\\
R^{\pi}(0)=0.
\end{cases}
\end{align}
Note that $R^{\pi}(t)$ may take negative values, meaning that the position of the investor slides into a loss.

\par
The benchmark growth rate, denoted by $\wp$, is a given ${\ensuremath{\mathcal{F}}}_T$-measurable random variable. The corresponding terminal benchmark wealth is $x_0e^{\wp}$.
In practice, $\wp$ can be chosen as the return of a stock index or structured product, the deposit or loan interest rate, or a constant target rate, etc.

We assume that the investor will not use any strategy that underperforms the benchmark significantly. Mathematically, we impose the following constraint:
\begin{align}\label{lowerbound}
R^{\pi}(T)-\wp \geq -c,
\end{align}
where $c>0$ is a constant that represents the \emph{risk tolerance level}.
To exclude trivial cases, we assume
\begin{align}\label{lowerbound1}
\text{ there exists at least one admissible portfolio $\pi_{0}$ such that } R^{\pi_{0}}(T)-\wp > -c.
\end{align}
 We will show that this requirement is equivalent to \citeassmp{feasibility} below, which can be easily verified in terms of $\wp$, $c$, and market parameters.

Now we present our optimal portfolio selection problem with behavioral criterion on the terminal relative growth rate:
\begin{align}\label{Optimalcontrol0}
\max_{\pi\in\mathscr{A}} & \quad \BV{u(R^{\pi}(T)-\wp)}\\
\mathrm{s.t.} &\quad \text{$(\pi, R^{\pi})$ satisfies \eqref{returnprocess} and \;}R^{\pi}(T)-\wp \geq -c. \nonumber
\end{align}
In this model, the investor adopts the CPT preference $\BV{u(\cdot)}$ to evaluate the terminal relative growth rate. Here, $\BV{\cdot}$ is a nonlinear expectation (indeed a Choquet expectation) and $u$ is an $S$-shaped utility function, which will be specified in the next section.

\blue{It is important to point out that neither the constraints nor the objective of Problem \eqref{Optimalcontrol0} depend on the initial endowment $x_{0}$, as we consider proportional investment strategies and logarithmic returns.
Consequently, the optimal portfolio $\pi^*$ and the corresponding growth rate $R^{\pi^*}$ must be independent of $x_0$, meaning that they are the same for all investors regardless of initial endowments.}

\subsection{CPT preference}
\noindent

In our model, the investor uses a CPT preference $\BV{u(\cdot)}$. It depends on a nonlinear expectation $\BV{\cdot}$ and a non-concave utility function $u$.

\par
First, we introduce the nonlinear expectation $\BV{\cdot}$ induced by a probability weighting (or distortion) function. A function $w\colon [0,1]\mapsto[0,1] $ is called a probability weighting function if it is strictly increasing and continuously differentiable with
 $w(0)=0$ and $w(1)=1$ (three types of probability weighting functions are shown in Figure \ref{figure:distortion}).
The expectation $\BV{\cdot}$ for a r.v. $\xi$ is a Choquet expectation defined by
\begin{equation}\label{Choquet integral}
\BV{\xi}=\int_0^{\infty}w(1-F_{\xi}(x))\ensuremath{\operatorname{d}\! x}-\int_{-\infty}^0 (1-w(1-F_{\xi}(x)))\ensuremath{\operatorname{d}\! x},
\end{equation}
provided that at least one of the integrals is finite. The Choquet expectation $\BV{\cdot}$ is in general nonlinear, unless there is no probability weighting (i.e. $w$ is the identity function) in which case $\BV{\cdot}$ becomes the standard linear mathematical expectation $\BE{\cdot}$ and the CPT preference $\BV{u(\cdot)}$ becomes the classical EUT preference $\BE{u(\cdot)}$. Choquet expectations have many applications in statistics, economics, and finance; see Schmeidler \cite{S89}, Wakker \cite{W01} and references therein.
\begin{remark}
Tversky and Kahneman \cite{TK92} used different probability weighting functions for the gain part (i.e., $\xi> 0$) and loss part (i.e., $\xi<0$). In this paper, we use the same probability weighting function for them. It is possible to consider different probability weighting functions for the two parts using the ``divide and conquer'' machinery developed by Jin and Zhou \cite{JZ08} to study the corresponding model. We leave this to interested readers.
\end{remark}
\begin{remark}
When there is no probability weighting, our model reduces to the classical EUT model with a non-concave utility function. The latter has been studied in Bian et al. \cite{BCX19} and Guan et al. \cite{GLXY17}, among others.
\end{remark}
\par
Second, the CPT involves an $S$-shaped utility function $u$ which is convex on $(-\infty, 0]$ and concave on $[0,\infty)$. Following the seminal work of Tversky and Kahneman \cite{TK92},
we adopt the piece-wise power utility function:
\begin{align}\label{valuefunction}
u(x)=\begin{cases}
x^{\alpha}, &\quad x\geq 0;\\
-\kappa(-x)^{\beta}, &\quad x<0,
\end{cases}
\end{align}
where $0<\alpha$, $\beta<1$ are risk aversion parameters, and $\kappa>0$ represents the degree of {loss aversion}.\footnote{The specific parameters studied in \cite{TK92} are $\alpha=\beta=0.88$ and $\kappa=2.5$.
} The utility function $u$ is continuous and strictly increasing on $(-\infty,\infty)$, $C^\infty$ smooth except at the point $x=0$, convex on $(-\infty, 0]$ and concave on $[0,\infty)$ (thus $S$-shaped).
One can generalize our model to more general utility functions at the cost of less explicit solutions.

From the CPT preference $\BV{u(\cdot)}$ introduced above, we see that, on one hand, it is fully nonlinear; on the other hand, it is non-concave. Due to the presence of nonlinear expectation $\BV{\cdot}$ in the target, the problem \eqref{Optimalcontrol0} is not a standard stochastic control problem (namely unlike those in Yong and Zhou \cite{YZ99}). The nonlinear expectation leads to the failure of dynamic programming and maximum principle, so the standard results and tools in the classical control theory cannot be applied to tackle it immediately. Although the delicate control theory of mean-field type may be adopted to study our problem, it seems very hard to get explicit solutions by that approach. One should notice that the Choquet expectation $\BV{\cdot}$ is essentially different from the nonlinear $g$-expectation introduced by Peng \cite{P97, P99}; see Chen et al. \cite{CCD05}. So the results and tools from the $g$-expectation theory can not be applied to our problem as well. We will provide closed-form solutions by relatively simple arguments.

\section{Martingale method}
As standard tools such as the dynamic programming in the control theory cannot be applied to tackle Problem \eqref{Optimalcontrol0}, we use the martingale method instead. The latter allows us to reduce the dynamic stochastic control problem to a static (but still non-concave) r.v. optimization problem. We need to assume the market is complete in order to use the martingale method.\footnote{We can deal with the problem in some special incomplete markets such as when all the market parameters are deterministic. In this case, the market is essentially equivalent to a complete one; see \cite{HZ11} for a more detailed discussion.}

We introduce the \emph{state price density process}
\[\rho(t)=\exp\left(-\int_{0}^{t}\Big(r(s)+\tfrac{1}{2}\Vert\theta(s)\Vert^{2}\Big)\ensuremath{\operatorname{d}\! s}-\int_{0}^{t}\theta(s)'\ensuremath{\operatorname{d}\! } W(s)\right),\quad t\in[0,T].\]
Then It\^{o}'s lemma gives
\[\ensuremath{\operatorname{d}\! } \rho(t)=-\rho(t)\big(r(t)\ensuremath{\operatorname{d}\! t}+\theta(t)'\ensuremath{\operatorname{d}\! } W(t)\big).\]
We denote $\rho=\rho(T)$\footnote{In the financial economics literature, $\rho(T)$ is often called the pricing kernel.} from now on.
For any admissible portfolio $\pi$, we have, by It\^{o}'s lemma and \eqref{returnprocess},
\begin{align}
\ensuremath{\operatorname{d}\! }\; (\rho(t)e^{R^{\pi}(t)})=\rho(t)e^{R^{\pi}(t)}\big(\pi(t)^{\prime}\sigma(t)-\theta(t)^{\prime}\big)\ensuremath{\operatorname{d}\! } W(t).
\end{align}
Hence, the process $\rho(t)e^{R^{\pi}(t)}$ is a local martingale. Because it is a positive process, it is a supermartingale. Thus
\begin{align}\label{budget1}
\BE{\rho e^{R^{\pi}(T)}}\leq 1.
\end{align}
This is called the \emph{budget constraint} for the investor.\footnote{In the financial economics literature, the budget constraint often refers to $\BE{\rho X^{\pi}(T)}\leq x_{0}$ which is equivalent to \eqref{budget1}.} For any admissible portfolio, the corresponding terminal growth rate must obey the budget constraint.
We can add the budget constraint into Problem \eqref{Optimalcontrol0} without changing its nature:
\begin{align}\label{Optimalcontrol}
\max_{\pi\in\mathscr{A}} & \quad \BV{u(R^{\pi}(T)-\wp)}\\
\mathrm{s.t.} &\quad \text{$(\pi, R^{\pi})$ satisfies \eqref{returnprocess}}, \quad \BE{\rho e^{R^{\pi}(T)}}\leq 1,\quad R^{\pi}(T)-\wp \geq -c.\nonumber
\end{align}

The following result plays a key role in transforming the dynamic stochastic optimal control problem \eqref{Optimalcontrol} to a static r.v. optimization problem.

\begin{lemma}\label{martingalemethod}
For any random variable $\xi$ with $\BE{\rho e^{\xi}}= 1$, there exists an admissible portfolio $\pi$ such that $R^{\pi}(T) = \xi$.
\end{lemma}
All proofs are placed in \citeapp{proofs}.

With the help of this result, solving the problem \eqref{Optimalcontrol} reduces to solving the following static r.v. optimization problem:
\begin{align}\label{Optimalcontrol1}
\max_{\xi} & \quad \BV{u(\xi-\wp)}\\
\mathrm{s.t.} &\quad \BE{\rho e^{\xi}}\leq 1,\quad\xi-\wp \geq -c.\nonumber
\end{align}

\begin{lemma}\label{martingalemethod2}
An admissible portfolio $\pi$ is optimal to Problem \eqref{Optimalcontrol} if and only if $R^{\pi}(T)$
is optimal to Problem \eqref{Optimalcontrol1}.
\end{lemma}

We can now rewrite \eqref{Optimalcontrol1} in a more compact form:
\begin{align}\label{Obj3}
\max_{\zeta} &\quad \BV{v(\zeta)}\\
\mathrm{s.t.} &\quad \BE{\zeta \eta}\leq 1,\quad \zeta\geq \hat{c}, \nonumber
\end{align}
where
\[v(x)=u\left( \log x \right),\quad \zeta=e^{\xi-\wp},\quad \eta=\rho e^{\wp},\quad \hat{c}=e^{-c}.\]
Here $\eta$ is termed the \emph{benchmark adjusted pricing kernel} under the benchmark $\wp$.
Moreover, $\zeta$ represents the terminal wealth relative to the benchmark, $X^{\pi}(T)/(x_{0}e^{\wp})$.

In \eqref{Obj3}, we can view $\BV{v(\cdot)}$ as a preference for $\zeta$, which is neither a CPT nor a rank-dependent utility theory (RDUT) preference due to the specific shape of $v$. In fact, $v$ is $M$-shaped as shown below, as opposed to the $S$-shaped utility in CPT or the concave utility in RDUT.
Recall that $\zeta$ stands for the terminal relative wealth, we term $v$ the \emph{relative utility function} and $\BV{v(\cdot)}$ the \emph{relative CPT preference}.

Finally, we present an equivalent condition to fulfill the requirement \eqref{lowerbound1}.
\begin{assmp}[Feasibility conditioin]\label{feasibility}
The trio of the pricing kernel $\rho$, the benchmark $\wp$ and the risk tolerance level $c>0$ satisfies
\[ \BE{\rho e^{\wp}} < e^{c}.\]
\end{assmp}
In fact, the requirement \eqref{lowerbound1} is fulfilled if and only if \citeassmp{feasibility} holds.
First, \eqref{lowerbound1} together with the budget constraint \eqref{budget1} clearly implies \citeassmp{feasibility} holds. On the other hand, by
\citelem{martingalemethod}, we can find $\pi_{0}$ such that
\[R^{\pi_{0}}(T)=\wp-\log \left(\BE{\rho e^{\wp}}\right),\]
which together with \citeassmp{feasibility} gives \eqref{lowerbound1} immediately. Therefore, \eqref{lowerbound1} and \citeassmp{feasibility} are equivalent.
Henceforth, we assume \citeassmp{feasibility} holds. \citeassmp{feasibility} is also equivalent to $\hat{c}\BE{\eta}<1$, so we see that there are infinitely many feasible solutions to Problem \eqref{Obj3}.

\begin{remark}
If the benchmark is constant and the interest rate process is deterministic, then \citeassmp{feasibility} is equivalent to $\wp<c +\int_{0}^{T}r(s)\ensuremath{\operatorname{d}\! } s$. This implies that the investor can not set up an extremely high target of $\wp$ for a fixed level of risk tolerance $c$; otherwise the problem may be ill-formulated, which is consistent with the saying that ``a higher return comes with a higher risk''.
\end{remark}

\section{Quantile formulation}

In this section, we focus on Problem \eqref{Obj3}. The Choquet expectation \eqref{Choquet integral} involves two elements: the probability weighting function $w(\cdot)$ and the probability distribution function $F_{\xi}$. The former imposes great complexity if we attempt to evaluate the nonlinear expectation directly. Expressing the objective in terms of the quantile of the decision variable will enable us to treat the two elements separately. This is the motivation and advantage of the so-called \emph{quantile formulation}, which has been used by Carlier and Dana \cite{CD08} to study risk-sharing problems, Jin and Zhou \cite{JZ08} and Xia and Zhou \cite{XZ16} to solve portfolio selection problems under behavioral criterion, Hou and Xu \cite{HX16} to study mean-variance portfolio optimization with intractable claims, and Wei \cite{W18,W21} to solve utility maximization with risk constraints. Notably, the third author \cite{X16} provided a simple form of this technique. Here we apply Xu's method to solve Problem \eqref{Obj3}.

For a r.v. $\xi$ with the probability distribution $F_{\xi}$, its quantile function is defined as
\[Q_{\xi}(p):=\inf\{z\in \mathbb{R}\;|\;F_{\xi}(z)> p\},\quad p \in(0,1), \]
with the convention that $Q_{\xi}(0)=\lim_{p\to 0}Q_{\xi}(p)$
and $Q_{\xi}(1)=\lim_{p\to 1}Q_{\xi}(p)$.
It is easy to see that $Q_{\xi}$ is increasing and right-continuous on $(0,1)$.
On the other hand, if $Q$ is increasing and right-continuous on $(0,1)$, then $Q=Q_{\xi}$, where $\xi=Q(U)$ with $U$ being uniformly distributed on $(0,1)$.
Hence,
\[\setq=\big\{Q:(0,1)\to\ensuremath{\operatorname{\mathbb{R}}} \;\big|\; \text{$Q$ is increasing and right-continuous} \big\},\]
denotes the set of quantile functions for all r.v.s.

In terms of the quantile function, we can express the Choquet expectation $\BV{\xi}$ as
\[\BV{\xi}=\int_{0}^{1}Q_{\xi}(p)w'(1-p)\ensuremath{\operatorname{d}\! p},\]
which leads to
\begin{align} \label{Obj4}
\BV{v(\zeta)}&=\int_{0}^{1}Q_{v(\zeta)}(p)w'(1-p)\ensuremath{\operatorname{d}\! p} =\int_{0}^{1}v(Q_{\zeta}(p))w'(1-p)\ensuremath{\operatorname{d}\! p}.
\end{align}
By definition, we have $\zeta\geq \hat{c}$ if and only if
\begin{align}\label{quantileconstraint1}
Q_{\zeta}(0)\geq \hat{c}.
\end{align}
Therefore, the objective functional and constraint $\zeta\geq \hat{c}$ in Problem \eqref{Obj3} can be expressed in terms of the quantile function of $\zeta$. We then need to express the other constraint $ \BE{\zeta \eta}\leq 1$ in terms of the quantile function of $\zeta$. This is impossible for general $\zeta$ since $\BE{\zeta \eta}$ depends not only on the quantile of $\zeta$, but also on the joint law of $\zeta$ and $\eta$. Nevertheless, we do not need such expressions for all $\zeta$ and it suffices to obtain them for the optimal solution.
\par
In the neoclassical models (such as EUT) and behavioral finance models for terminal wealth (such as Jin and Zhou \cite{JZ08} and He and Zhou \cite{HZ11}), the optimal solution $\zeta$ must be anti-comonotonic with the pricing kernel $\rho$, meaning that the optimal terminal wealth will be higher if the market status is getting better (i.e., $\rho$ is getting smaller).
This may not be the case in our model since Problem \eqref{Obj3} relies on $\eta$ instead of $\rho$. In the same spirit of Xu \cite{X14}, we can show that any optimal candidate $\zeta$ for Problem \eqref{Obj3} must be anti-comonotonic with $\eta$. Hence, $\zeta$ is not anti-comonotonic with $\rho$ unless $\rho$ and $\eta$ are comonotonic. Intuitively, the benchmark should be higher if the market status is getting better, e.g., when the benchmark is the return of a stock index. Mathematically, this means $\wp$ is decreasing in $\rho$. In particular, if the benchmark is chosen as $\wp=k\log\rho^{-1}$ with $k>0$, then $\zeta$ is indeed anti-comonotonic with $\rho$ when $k<1$ and comonotonic when $k>1$. Therefore, the benchmark may significantly influence people's investment strategy. When one setup an aggressive target (namely $k>1$), the output may be opposite to the classical observation that one should have a higher relative return in a good market scenario.

The above analysis eventually leads to the following quantile formulation of Problem \eqref{Obj3}:
\begin{align}\label{quantile1}
\max_{Q\in\setq} &\quad \int_{0}^{1}v(Q(p))w'(1-p)\ensuremath{\operatorname{d}\! p} \\
\mathrm{s.t.} &\quad \int_{0}^{1}Q(p)Q_{\eta}(1-p)\ensuremath{\operatorname{d}\! p}=1,\quad Q(0)\geq \hat{c}.\nonumber
\end{align}

We now establish the relationship between Problem \eqref{quantile1} and Problem \eqref{Obj3}.
Let $\setu$ be the set of random variables that are uniformly distributed on $(0,1)$ and comonotonic with $\eta$.
Then, $\setu$ is not empty, and $\eta=Q_{\eta}(U)$ a.s. for any $U\in\setu$. Moreover, the distribution function of $\eta$, $F_{\eta}$, is continuous if and only if $\setu=\{F_{\eta}(\eta)\}$ (see, e.g. Xu \cite{X14}).
\begin{prop}\label{quantileformulation}
A random variable $\zeta^{*}$ is an optimal solution to Problem \eqref{Obj3}, if and only if
$\zeta^{*}=Q^{*}(1-U)$ a.s., where $U\in\setu$ and $Q^{*}$ is an optimal solution to Problem \eqref{quantile1}.
\end{prop}

Following Xu \cite{X16}, we introduce the following notations to simplify the expressions in Problem \eqref{quantile1},
\begin{align
\nu(p)&=1-w^{-1}(1-p),\quad p\in[0, 1],\nonumber \\
G(p) &=Q(\nu(p)),\quad p\in(0, 1), \nonumber\\
\varphi(p)&=-\int_{p}^{1}Q_{\eta}(1-\nu(s))\nu'(s)\ensuremath{\operatorname{d}\! s},\quad p\in[0,1]. \label{varphi}
\end{align}
Then, Problem \eqref{quantile1} can be rewritten as (see \citeapp{changeofvariable})
\begin{align}\label{quantile2}
\max_{G\in\setq} &\quad \int_{0}^{1}v(G(p))\ensuremath{\operatorname{d}\! p} \\
\mathrm{s.t.} &\quad \int_{0}^{1}G(p)\varphi'(p)\ensuremath{\operatorname{d}\! p}=1,\quad G(0)\geq \hat{c}. \nonumber
\end{align}

As will be clear shortly, $v$ is not concave and thus Problem \eqref{quantile2} is not a concave optimization problem. Therefore, the Lagrange method can not be applied directly to tackle it. This is different from the problem in Xu \cite{X16} wherein the constraint $G(0)\geq \hat{c}$ is absent and the $v$ is concave. Hence, Xu's \cite{X16} result is not readily applicable and a new method is called for.

\section{Optimal solution}
We intend to use the concavification method to solve the problem \eqref{quantile2}.
Before doing so, let us first study the properties of the function $v$.
\subsection{Shape of the relative utility $v$}\label{UOR}
\noindent
We start with the shape of $v$.
Note that
\begin{align*}
v(x)=u(\log x)=\begin{cases}
(\log x)^{\alpha}, &\quad x\geq 1;\\
-\kappa(-\log x)^{\beta}, &\quad 0< x<1.
\end{cases}
\end{align*}
Hence we have
\[v'(x)=u'(\log x)x^{-1}>0, \quad x\neq 1,\]
and
\[v''(x)=[u''(\log x)-u'(\log x)]x^{-2},\quad x>0.
\]
We see from the above relations that $v$ would be increasing and global concave if so was $u$, which is
the case in EUT and RDUT. By contrast, the shape of $v$ in our model is essentially different from theirs. Indeed,
\begin{itemize}
\item When $x>1$, we have $u''(\log x)<0$ and $u'(\log x)>0$, so $v''(x)<0$ and $v$ is strictly concave.
\item When $0<x<1$, we have
\begin{align*}
v''(x)&=[u''(\log x)-u'(\log x)]x^{-2}\\
&=-\kappa\beta\big((\beta-1)+(- \log x)\big)(-\log x)^{\beta-2} x^{-2}\\
&\begin{cases}
>0, &\quad e^{\beta-1}<x<1;\\
<0, &\quad 0<x< e^{\beta-1}.
\end{cases}
\end{align*}
So $v$ is strictly concave on $(0, e^{\beta-1}]$ and strictly convex on $[e^{\beta-1},1)$.
\end{itemize}
To summarize, the function $v$ is continuous, strictly increasing on $(-\infty,\infty)$, strictly concave, respectively, on $(0, e^{\beta-1}]$ and $[1,\infty)$, and strictly convex on $[e^{\beta-1},1]$. We call such \emph{concave-convex-concave-shaped} utility an \emph{$M$-shaped} utility function.
Since $\hat{c}$ and $\beta$ are entirely decided by the investor's preference, both $\hat{c} \leq e^{\beta-1}$ (\citefig{fig3.1}) and $\hat{c}> e^{\beta-1}$ (\citefig{fig3.2}) are possible.
\begin{figure}[H]
\begin{center}
\begin{picture}(280,160)
\linethickness{0.6pt}
\put(40,89){\vector(1,0){190}}
\put(50,0){\vector(0,1){150}}
\put(215,78){$x$}
\qbezier(54,14)(57,39)(96,51)
\qbezier(96,51)(136,60)(144,88)
\qbezier(144,88)(149,120)(209,138)
\put(200,122){$v(x)$}
\put(145,78){$1$}

\put(75, 92){$e^{\beta-1}$}
\qbezier[30](80, 45)(80, 66)(80, 88)
\put(62, 92){$\hat{c}$}
\qbezier[38](65, 35)(65, 62)(65, 88)
\end{picture}
\caption{The function $v$ with $\hat{c} \leq e^{\beta-1}$.}\label{fig3.1}
\label{fig4}
\end{center}
\end{figure}

\begin{figure}[H]
\begin{center}
\begin{picture}(280,160)
\linethickness{0.6pt}
\put(40,89){\vector(1,0){190}}
\put(50,0){\vector(0,1){150}}
\put(215,78){$x$}
\qbezier(54,14)(57,39)(96,51)
\qbezier(96,51)(136,60)(144,88)
\qbezier(144,88)(149,120)(209,138)
\put(200,122){$v(x)$}
\put(145,78){$1$}

\put(75, 92){$e^{\beta-1}$}
\qbezier[30](80, 45)(80, 66)(80, 88)
\put(102, 92){$\hat{c}$}
\qbezier[28](105, 54)(105, 76)(105, 88)
\end{picture}
\caption{The function $v$ with $\hat{c} > e^{\beta-1}$.}\label{fig3.2}
\label{fig4}
\end{center}
\end{figure}

\par
Because $v$ is $M$-shaped, Problem \eqref{quantile2} is not a concave or $S$-shaped utility optimization problem like those in Jin and Zhou \cite{JZ08} and Xu \cite{X16}. The $M$-shaped utility has been rarely studied in the financial economics literature. The earliest one we can find dates back to Friedman and Savage \cite{FS48} that initially proposed such shape of utility to explain the purchasing of both insurance and lottery tickets (the ``F-S hypothesis''). As demonstrated in \cite{FS48}, the concave segment indicates a preference for certainty (insurance) while the convex segment indicates a preference for risk (gambling). In their paper, the $M$-shaped utility was endowed with a reasonable interpretation that a lower socioeconomic class consumer whose income placed corresponding to the first concave segment wishes to shift himself up to a higher socioeconomic class whose income placed corresponding to the last concave segment, the convex segment with increasing marginal utility represented a transitional stage between the two classes.

Although the shape of the utility $v$ in our model is also $M$-shaped, there are some differences we would like to highlight. First, in terms of the meaning of utility, it refers to the utility of ``income'' in the context of \cite{FS48}. Whereas the ``utility'' investigated underline refers to the utility of the relative growth rate $\zeta=X^{\pi}(T)/x_{0}e^{\wp}$. As one may note, the form of the relative utility underline relies heavily on the hypothesis of the piece-wise utility of wealth $u$. The shape of the relative utility $v$ is entirely determined by the parameters $\alpha,\beta,\kappa$. It is possible to have other kinds of non-concave, non-$M$-shaped shaped relative utilities $v$ if one would use general $S$-shaped or other shaped utilities $u$.

More interestingly, as explained in \cite{FS48}, $M$-shape utility is only plausible for consumer units whose level of income is in the first concave segment of the utility function. If that was so, then, in turn, the ``F-S'' hypothesis implies that the investor's ability to invest, measured by the relative return variable $\zeta$, should take values in the first concave segment, namely the region $(0,e^{\beta-1})$. In other words,, the gross return of investor should subject to $0<{X^{\pi}(T)}/{x_0}<e^{\wp+\beta-1}$. If we set $\beta=0.88$ and $\wp$ to the common value of risk-free interest rate (normally less than $0.05$), then these investors would typically refer to those who suffer losses in the market as $e^{\wp+\beta-1}<1$.

\subsection{Concavified problem}
\noindent
One common approach to overcome the non-concavity in Problem \eqref{quantile2} is to consider its \emph{concavified} problem. This approach has been used to solve non-concave utility maximization problems by Bian et al. \cite{BCX19} and Guan et al. \cite{GLXY17} for expected utility models, \blue{by Bi et al. \cite{BYCF21} for behavioral models.}

Let us define the concave envelope of $v$ on $(0,\infty)$, that is, the smallest concave function dominating $v$ on $(0,\infty)$, denoted by $\hat{v}_{0}$. Mathematically, it is given by
\[\hat{v}_{0}(x)=\sup_{\substack{0<y\leq x\leq z\\ y\neq z}}\left(\frac{z-x}{z-y}v(y)+\frac{x-y}{z-y}v(z)\right),\quad x>0.\]
We have the following explicit expression for $\hat{v}_{0}$.
\begin{lemma}\label{twopoints}
There are two scalars $0<a<1<b$ such that the function $\hat{v}_{0}$ is concave on $(0,\infty)$, coincides with $v$ on $(0, a]\cup[b,\infty)$, and is linear on $[a,b]$. Moreover,
\begin{align}\label{abequation}
v'(a)=v'(b)=\frac{v(b)-v(a)}{b-a}.
\end{align}
\end{lemma}

The locations of $a$ and $b$ are demonstrated in \citefig{fig4}. They are determined by the utility parameters $\alpha,\beta$, and $\kappa$.
The global concave envelope function $\hat{v}_{0}$ is different from $v$ on the interval $(a,b)$, shown in the red dot line.
\begin{figure}[H]
\begin{center}
\begin{picture}(280,160)
\linethickness{0.6pt}
\put(40,89){\vector(1,0){190}}
\put(50,0){\vector(0,1){150}}
\put(215,78){$x$}
\qbezier(54,14)(57,39)(96,51)
\qbezier(96,51)(136,60)(144,88)
\qbezier(144,88)(149,120)(209,138)
\put(200,122){$v(x)$}
\put(145,78){$1$}

\linethickness{0.9pt}
{{\color{red}\qbezier[80](69, 41)(113, 79)(158, 116)}}
\qbezier[25](68,41)(68, 65)(68,88)
\put(65,92){$a$}
\qbezier[14](158,116)(158,103)(158,90)
\put(156,78){$b$}
\end{picture}
\caption{The function $\hat{v}_{0}$ in red dot line on $(a,b)$.}
\label{fig4}
\end{center}
\end{figure}

A naive way to solve Problem \eqref{quantile2} is described as follows. We first try to solve the following \emph{globally concavified problem}:
\begin{align}\label{quantile3}
\max_{G\in\setq } &\quad \int_{0}^{1}\hat{v}_{0}(G(p))\ensuremath{\operatorname{d}\! p} \\
\mathrm{s.t.} &\quad \int_{0}^{1}G(p)\varphi'(p)\ensuremath{\operatorname{d}\! p}=1,\quad G(0)\geq \hat{c}. \nonumber
\end{align}
We then try to show that the optimal solution to Problem \eqref{quantile3} is also optimal to Problem \eqref{quantile2}. Unfortunately, this approach fails due to the lower bound constraint $Q(0)\geq \hat{c}$.
In fact, one can easily see from the formulation of Problem \eqref{quantile2} that its optimal solution depends on the utility function $v$ only on $[\hat{c}, \infty)$. However, in Problem \eqref{quantile3}, $\hat{v}_{0}$ depends on $v$ on $(0,\infty)$, so these two problems are unlikely to share the same optimal solution. It is natural to consider the local concave envelope of $v$ on $[\hat{c},\infty)$ rather than the global one $\hat{v}_{0}$ on $(0,\infty)$.
\par
The above discussion motivates us to consider the concave envelope function of $v$ on $[\hat{c},\infty)$ that is denoted by $\hat{v}_{\qlow}$. It is the smallest concave function dominating $v$ on $[\hat{c},\infty)$, called the local concave envelope of $v$ and given by
\begin{align}\label{hatvdef}
\hat{v}_{\qlow}(x)=\sup_{\substack{\hat{c}\leq y\leq x\leq z\\ y\neq z}}\left(\frac{z-x}{z-y}v(y)+\frac{x-y}{z-y}v(z)\right),\quad x\geq \hat{c}.
\end{align}
We now introduce the \emph{locally concavified problem}:
\begin{align}\label{q1}
\max_{G\in\setq } &\quad \int_{0}^{1}\hat{v}_{\qlow}(G(p))\ensuremath{\operatorname{d}\! p} \\
\mathrm{s.t.} &\quad \int_{0}^{1}G(p)\varphi'(p)\ensuremath{\operatorname{d}\! p}=1,\quad G(0)\geq \hat{c}. \nonumber
\end{align}
This is a concave optimization problem
which is easier to solve compared to the non-concave optimization Problem \eqref{quantile2}.
The two problems look different but we will show that they share the optimal solution and optimal value,
\blue{
provided that they are well-posed (that is, their optimal values are finite); see \citeprop{wellposedness}.
}
\par
Before proving the above claim, let us first study the properties of the local concave envelope function $\hat{v}_{\qlow}$.
Clearly, $\hat{v}_{\qlow}\leq \hat{v}_{0}$ on $[\hat{c},\infty)$. It is also easy to show that they are identical on $[\hat{c},\infty)$ if and only if $\hat{c}\leq a$; see \citefig{fig4.5}.

\begin{figure}[H]
\begin{center}
\begin{picture}(280,160)
\linethickness{0.6pt}
\put(40,89){\vector(1,0){190}}
\put(50,0){\vector(0,1){150}}
\put(215,78){$x$}
\qbezier(54,14)(57,39)(96,51)
\qbezier(96,51)(136,60)(144,88)
\qbezier(144,88)(149,120)(209,138)
\put(200,122){$v(x)$}

\linethickness{0.9pt}
{{\color{red}\qbezier[80](69, 41)(113, 79)(158, 116)}}
\qbezier[25](68,41)(68, 65)(68,88)
\put(65,92){$a$}
\qbezier[14](158,116)(158,103)(158,90)
\put(156,78){$b$}

\qbezier[30](55,28)(55,58)(55,88)
\put(52,92){$\hat{c}$}
\end{picture}
\caption{The functions $\hat{v}_{0}$ and $\hat{v}_{\qlow}$ coincide, when $\hat{c}\leq a$.}
\label{fig4.5}
\end{center}
\end{figure}

If $a<\hat{c}$, then the function $\hat{v}_{\qlow}$ coincides with $v$ on $[d,\infty)$ and is linear on $[\hat{c},d]$ for some $1<d<b$. Moreover, $\hat{v}_{\qlow}< \hat{v}_{0}$ on $[\hat{c}, b]$ and $\hat{v}_{\qlow}>v$ on $(\hat{c}, d)$; see \citefig{fig5}. Similar as before, we can find the value of $d$ via an algebraic equation
\[v'(d)=\frac{v(d)-v(\hat{c})}{d-\hat{c}}, \quad\hat{c}>a.\]

\begin{figure}[H]
\begin{center}
\begin{picture}(280,160)
\linethickness{0.6pt}
\put(40,89){\vector(1,0){190}}
\put(50,0){\vector(0,1){150}}
\put(215,78){$x$}
\qbezier(54,14)(57,39)(96,51)
\qbezier(96,51)(136,60)(144,88)
\qbezier(144,88)(149,120)(209,138)
\put(200,122){$v(x)$}

\linethickness{0.9pt}
{{\color{red}\qbezier[80](69, 41)(113, 79)(158, 116)}}
\qbezier[25](68,41)(68, 65)(68,88)
\put(65,92){$a$}
\qbezier[14](158,116)(158,103)(158,90)
\put(156,78){$b$}

\qbezier[20](100,54)(100,71)(100,88)
\put(97,92){$\hat{c}$}
{{\color{blue}\qbezier[47](100,53)(125,81)(150, 109)}}
\qbezier[12](150,90)(150,99)(150,108)
\put(146,78){$d$}
\end{picture}
\caption{The function $\hat{v}_{0}$ in red dot line and $\hat{v}_{\qlow}$ in blue dot line, when $\hat{c}>a$.}
\label{fig5}
\end{center}
\end{figure}

From an economic perspective, the larger the risk tolerance level $c$ (i.e., the smaller the value of $\hat{c}$), the smaller the local concave envelope function $\hat{v}_{\qlow}$ as well as the optimal value of Problem \eqref{quantile2}.
\par
Because $\hat{v}_{\qlow}'$ is continuous and decreasing on $[\hat{c},\infty)$, we may define its right-continuous inverse function as
\begin{align}\label{Leftinverse}
I(x):=\inf\{y\geq \hat{c}\mid \hat{v}_{\qlow}'(y)\leq x\}\vee \hat{c},\quad x>0.
\end{align}
Its properties are summarized in the following lemma.
\begin{lemma}\label{Iprop}
The function $I$ is decreasing and right-continuous, and satisfies the following properties.
\begin{enumerate}
\item If $a\leq \hat{c}$, then
\begin{align}\label{I1}
I(x)=\begin{cases}
(v')^{-1}(x), &\quad \text{ if $0<x<v'(d)$};\\
\hat{c},&\quad \text{ if $x\geq v'(d)$}.
\end{cases}
\end{align}

\item
If $a>\hat{c}$, then
\begin{align}\label{I2}
I(x)=\begin{cases}
(v')^{-1}(x), &\quad \text{ if $0<x<v'(a)$};\\
a,&\quad \text{ if $x=v'(a)$};\\
(v')^{-1}(x), &\quad \text{ if $v'(a)<x<v'(\hat{c})$};\\
\hat{c},&\quad \text{ if $x\geq v'(\hat{c})$}.
\end{cases}
\end{align}

\item For any $x>0$, \[\hat{v}_{\qlow}(I(x))=v(I(x)). \]
\item For any $x>0$, \[\max_{y\geq \hat{c}}(v(y)-xy)=\max_{y\geq \hat{c}}(\hat{v}_{\qlow}(y)-xy)=\hat{v}_{\qlow}(I(x))-xI(x).\]
\end{enumerate}
\end{lemma}
In particular, $I$ has a unique jump at $x=v'(d)$ when $a\leq \hat{c}$ or at $x=v'(a)$ when $a> \hat{c}$, and it does not take values in $(\hat{c},d)$ when $a\leq \hat{c}$ or in $(a,b)$ when $a> \hat{c}$.

\subsection{Optimal solution}

We now use employ the method of Xu \cite{X16} to solve all the problems formulated in the previous sections, including Problems \eqref{Optimalcontrol0}, \eqref{quantile2} and \eqref{q1}.

Let $\widehat{\varphi}$ be the concave envelope of $\varphi$ on $[0,1]$, that is,
\begin{align}\label{deltadef}
\widehat{\varphi}(p):=\sup_{\substack{0\leq x\leq p\leq y\leq 1\\ x\neq y}}\left(\frac{y-p}{y-x}\varphi(x)+\frac{p-x}{y-x}\varphi(y)\right),\quad p\in[0,1].
\end{align}
Because $\varphi(1)>\varphi(p)$ for all $p\in(0,1)$, we conclude that $\widehat{\varphi}$ is strictly increasing on $(0,1)$.
Let $\widehat{\varphi}'$ denote the right-continuous derivative of $\widehat{\varphi}$, that is,
\[\widehat{\varphi}'(p):=\lim_{\varepsilon\to 0+}\frac{\widehat{\varphi}(p+\varepsilon)-\widehat{\varphi}(p)}{\varepsilon},\quad p \in(0,1).\]
Then $\widehat{\varphi}'>0$ on $(0,1)$.

\blue{
The following proposition characterizes the well-posedness of Problems \eqref{Optimalcontrol0}, \eqref{quantile2} and \eqref{q1}.
\begin{prop}\label{wellposedness}
The following statements are equivalent.
\begin{enumerate}
\item Problem \eqref{Optimalcontrol0} is well-posed.
\item Problem \eqref{quantile2} is well-posed.
\item Problem \eqref{q1} is well-posed.
\item It holds that
\begin{align}\label{well1}
\int_0^1 \big(\max\{-\log \widehat{\varphi}'(p),\;0\}\big)^{\alpha} \ensuremath{\operatorname{d}\! p}<\infty.
\end{align}
\end{enumerate}
\end{prop}
}

The following theorem presents the optimal solution to our model.

\begin{thm}[Optimal solution]\label{vt}
Suppose \eqref{well1} holds.
If there exists a constant $\lambda>0$ such that
\begin{align}\label{lambda1}
\int_{0}^{1}I(\lambda \widehat{\varphi}'(p))\widehat{\varphi}'(p)\ensuremath{\operatorname{d}\! p}=1,
\end{align}
then $$G^{*}(p)=I(\lambda \widehat{\varphi}'(p)),\quad p \in(0,1),$$ is optimal to Problems \eqref{q1} and \eqref{quantile2}. Moreover,
$$Q^{*}(p)=I(\lambda \widehat{\varphi}'(1-w(1-p))),\quad p \in(0,1),$$ is optimal to Problem \eqref{quantile1}.

Furthermore, let $U\in\setu$. Then there exists an $\{\mathcal{F}_t\}$-progressively measurable process $Z$ such that
\[\BE{\rho(T) e^{\wp} G^{*}(1-w(U))\;\big|\;\mathcal{F}_t}=1+\int_0^t Z(s)'\ensuremath{\operatorname{d}\! } W(s), \quad t\in[0,T].\]
Furthermore,
\[\pi^*(t)=(\sigma'(t))^{-1}\frac{\rho(t) Z(t)}{\BE{\rho(T) e^{\wp} G^{*}(1-w(U))\;\big|\;\mathcal{F}_t}},\quad t\in[0,T],\]
is optimal to Problem \eqref{Optimalcontrol0} and $$R^{\pi^{*}}(T)=\wp+\log (G^{*}(1-w(U)))$$ is the corresponding optimal terminal growth rate.
\end{thm}

\blue{
\begin{remark}
Theorem \ref{vt} characterizes the optimal solution assuming the existence of a Lagrange multiplier $\lambda$ that solves \eqref{lambda1}. In Appendix \ref{existenceofLagrangemultiplier} we discuss when this assumption holds.
\end{remark}
}

Because $I$ does not take values in $(\hat{c},d)$ when $a\leq \hat{c}$ or in $(a,b)$ when $a> \hat{c}$, so does $G^{*}$. This indicates that the optimal terminal relative growth rate and wealth will not take values close to the benchmark. From an economic point of view, the investor shall take relatively large or small returns rather than those close to the benchmark. In other words, she would seek relatively large returns at the risk of getting relatively small returns.

In Berkelaar et al. \cite{BKP04} and Jin and Zhou \cite{JZ08}, the optimal terminal relative wealth turns out to consist of a pricing kernel-dependent gain part and a constant loss part with a strict gap between them.
In our model, when the agent's risk tolerance level $c$ is no more than the threshold $-\log a$ (i.e., $a\leq\hat{c}$), the same phenomenon happens by \eqref{I1}. However, when the risk tolerance level exceeds this threshold,
we see from \eqref{I2} that the loss part is no more a constant. Instead, it is also benchmark-adjusted pricing kernel dependent and takes values from $a$ to $\hat{c}$. This is a new financial phenomenon observed in our model. It shows that risk tolerance level plays a critical role in determining the optimal strategy.

Zhang et al. \cite{ZJZ11} considered a CPT model with a lower bound constraint on the terminal relative wealth. Their model is similar to ours, except that they maximize the utility of terminal relative wealth and use separate probability weighting functions for gain and loss parts. Different from our case, the loss part in their model consists of a constant moderate loss and a constant maximum loss due to different probability weighting functions for gain and loss parts. These two constants become the same when the probability weighting functions for gain and loss parts are the same (see \eqref{parawealth1} below).

In the next section, we make a comparison between their results and ours.

\section{A comparison}

We rewrite Problem \eqref{Optimalcontrol} as
\begin{align}\label{paramodel}
\max_{\pi\in\mathscr{A}} & \quad \BV{u(\log (X^{\pi}(T)/x_{0})-\wp)}.\\
\mathrm{s.t.} &\quad\BE{\rho X^{\pi}(T)}\leq x_0,\quad X^{\pi}(T)\geq x_{0}e^{\wp -c}. \nonumber
\end{align}
In Problem \eqref{paramodel}, the terminal return $R^{\pi}(T)$ in \eqref{Optimalcontrol} has been replaced by the traditional terminal wealth $X^{\pi}(T)$.

Zhang et al. \cite{ZJZ11} considered a CPT model with a lower bound constraint on the terminal relative wealth.
For comparison, we assume the agent takes the same terminal benchmark wealth at time $T$ as ours, given by $x_0e^{\wp}$; also, she uses the same lower bound.
Let $L$ denote the constant maximum loss given by
\[L=x_{0}e^{\wp}-x_{0}e^{\wp -c}.\]
Their problem can be formulated as
\begin{align}\label{paramodel1}
\max_{\pi\in\mathscr{A}'} & \quad \BV{u(X^{\pi}(T)-x_0e^{\wp})}.\\
\mathrm{s.t.} &\quad\BE{\rho X^{\pi}(T)}\leq x_{0},\quad X^{\pi}(T)-x_0e^{\wp}\geq -L, \nonumber
\end{align}
when the probability weighting functions for gain and loss parts are identical as in Problem \eqref{paramodel}.
Problems \eqref{paramodel1} and \eqref{paramodel} share the same constraint on the terminal wealth $X^{\pi}(T)$, but they have different objective functionals. The admissible portfolio sets of these two problems are also slightly different, but this will not affect our comparison.

After some tedious calculations, we can show that under our numerical setting to be specified below, the optimal terminal wealth to Problem \eqref{paramodel1} is
\begin{align}\label{parawealth1}
X^{\pi}(T)-x_0e^{\wp}=\begin{cases}
\;(u')^{-1}\Big( \frac{\lambda_1\rho}{w'(F_{\rho}(\rho))}\Big),&\quad \mbox{if } \frac{\lambda_1\rho}{w'(F_{\rho}(\rho))} \leq \alpha l^{\alpha-1};\\[10pt]
\; -L,&\quad \mbox{if } \frac{\lambda_1\rho}{w'(F_{\rho}(\rho))}>\alpha l^{\alpha-1},
\end{cases}
\end{align}
where $l$ solves
$$\alpha l^{\alpha-1} =\frac{l^\alpha + \kappa L^\beta }{l+L},$$
and $\lambda_1>0$ is the Lagrange multiplier that solves the budget constraint $\BE{\rho X^{\pi}(T)}= x_{0}$.

The optimal terminal wealth to Problem \eqref{paramodel} is given by
\begin{enumerate}
\item If $a\leq \hat{c}$, then
\begin{align}\label{parawealth2}
X^{\pi}(T)/x_0e^{\wp}=\begin{cases}
\;(v')^{-1} \Big( \frac{\lambda_2 e^{\wp} \rho}{w'(F_{\rho}(\rho))} \Big),&\quad \mbox{if} \frac{\lambda_2 e^{\wp} \rho}{w'(F_{\rho}(\rho))} \leq \frac{\alpha (\log d)^{\alpha-1}}{d};\\
\; \hat{c},&\quad \mbox{if } \frac{\lambda_2 e^{\wp} \rho}{w'(F_{\rho}(\rho))} > \frac{\alpha (\log d)^{\alpha-1}}{d}.
\end{cases}
\end{align}
\item If $a>\hat{c}$, then
\begin{align}\label{parawealth3}
X^{\pi}(T)/x_0e^{\wp}=\begin{cases}
(v')^{-1} \Big( \frac{\lambda_2 e^{\wp} \rho}{w'(F_{\rho}(\rho))} \Big) , &\quad \text{if } \frac{\lambda_2 e^{\wp} \rho}{w'(F_{\rho}(\rho))} \leq \frac{\kappa \beta (- \log a)^{\beta-1}}{ a};\\
(v')^{-1} \Big( \frac{\lambda_2 e^{\wp} \rho}{w'(F_{\rho}(\rho))} \Big), &\quad \text{if } \frac{\kappa \beta (- \log a)^{\beta-1}}{a} < \frac{\lambda_2 e^{\wp} \rho}{w'(F_{\rho}(\rho))} \leq \frac{\kappa \beta (- \log \hat{c})^{\beta-1}}{ \hat{c}};\\
\hat{c},&\quad \text{if } \frac{\lambda_2 e^{\wp} \rho}{w'(F_{\rho}(\rho))} > \frac{\kappa \beta (- \log \hat{c})^{\beta-1}}{ \hat{c}}.
\end{cases}
\end{align}
\end{enumerate}
where $\lambda_2>0$ is the Lagrange multiplier that solves the budget constraint.

The investor in Problem \eqref{paramodel1} classifies the future states into two scenarios: the good states when $\frac{\lambda_1\rho}{w'(F_{\rho}(\rho))} \leq \alpha l^{\alpha-1}$, in which case the payoff is decreasing in the state price density $\rho$; the bad states when $\frac{\lambda_1\rho}{w'(F_{\rho}(\rho))}>\alpha l^{\alpha-1}$, in which case the investor accepts a fixed loss. Moreover, there is a jump discontinuity in the terminal wealth at $\frac{\lambda_1\rho}{w'(F_{\rho}(\rho))}=\alpha l^{\alpha-1}$.

By contrast, the optimal terminal wealth to Problem \eqref{paramodel} exhibits two patterns depending on the relationship between $a$ and $\hat{c}$. If $a \leq \hat{c}$ (e.g., $c=0.1$ or $c=0.2$), the optimal terminal wealth resembles that to \eqref{paramodel1}: the good states when $\frac{\lambda_2 e^{\wp} \rho}{w'(F_{\rho}(\rho))} \leq \frac{\alpha (\log d)^{\alpha-1}}{d}$, in which case the payoff is decreasing in the state price density $\rho$; the bad states when $\frac{\lambda_2 e^{\wp} \rho}{w'(F_{\rho}(\rho))} > \frac{\alpha (\log d)^{\alpha-1}}{d}$, in which case the investor accepts a fixed loss; and there is a jump discontinuity in the terminal wealth at $\frac{\lambda_2 e^{\wp} \rho}{w'(F_{\rho}(\rho))} = \frac{\alpha (\log d)^{\alpha-1}}{d}$.

If $a > \hat{c}$ (e.g., $c=0.3$), the investor in our model classifies the future states into three scenarios: the good states when $\frac{\lambda_2 e^{\wp} \rho}{w'(F_{\rho}(\rho))} \leq \frac{\kappa \beta (- \log a)^{\beta-1}}{ a}$ and intermediate states when $\frac{\kappa \beta (- \log a)^{\beta-1}}{a} < \frac{\lambda_2 e^{\wp} \rho}{w'(F_{\rho}(\rho))} \leq \frac{\kappa \beta (- \log \hat{c})^{\beta-1}}{ \hat{c}}$, in which cases the payoff is decreasing in the state price density $\rho$ but there is a jump discontinuity at $\frac{\lambda_2 e^{\wp} \rho}{w'(F_{\rho}(\rho))} = \frac{\kappa \beta (- \log a)^{\beta-1}}{ a}$; the bad states when $\frac{\lambda_2 e^{\wp} \rho}{w'(F_{\rho}(\rho))} > \frac{\kappa \beta (- \log \hat{c})^{\beta-1}}{ \hat{c}}$, in which case the investor accepts a fixed loss.

For ease of illustration and comparison, we make the following assumptions in our numerical examples. The investment horizon is 1 year, i.e., $T=1$, and the initial wealth is $x_0=1$. The investment opportunity set is deterministic with $r=0.02$ and $|\theta|=0.2$, so $\rho$ is lognormally distributed, i.e., $\ln \rho \sim N(\mu_{\rho},\sigma_{\rho}^2)$ with $\mu_\rho = -(r+\frac{|\theta|^2}{2})T=-0.04$ and $\sigma_\rho = \sqrt{|\theta|^2T}=0.2$. We follow \cite{TK92} to set $\alpha=\beta=0.88$ and $\kappa=2.5$. In addition, we assume the benchmark growth rate $\wp$ is given by $\wp = (r + g)T$, where $g$ represents the growth rate in excess of the risk-free rate. In our numerical examples, $c$ will take values in $\{0.1,0.2,0.3\}$ and $g$ will take values in $\{0, 0.05, -0.05\}$.

We consider three types of probability weighting functions:
\begin{enumerate}
\item The identity function $w(p)=p, ~ p \in [0,1]$, i.e., there is no probability distortion;

\item The power function $w(p)=\sqrt{p}, ~ p \in [0,1]$;

\item The inverse $S$-shaped probability weighting function proposed in Jin and Zhou \cite{JZ08}:
\begin{equation*}
w(p) = \left\{
\begin{aligned}
k e^{(\bar{a}+\bar{b})\Phi ^{-1}(\bar{p})+\frac{\bar{a}^2}{2}} \Phi \left(\Phi ^{-1}(p)+\bar{a}\right), \quad & p\le \bar{p},\\
A + ke^{\frac{\bar{b}^2}{2}} \Phi \left(\Phi ^{-1}(p)-\bar{b} \right) , \quad & p> \bar{p},\\
\end{aligned}
\right.
\end{equation*}
where $\Phi$ is the standard normal probability distribution, $\bar{a}$ and $\bar{b}$ are two nonnegative parameters, and
$$k = \frac{1}{e^{\frac{\bar{b}^2}{2}}\Phi \left(-\Phi ^{-1}(\bar{p})+\bar{b} \right) + e^{(\bar{a}+\bar{b})\Phi ^{-1}(\bar{p})+\frac{\bar{a}^2}{2}} \Phi \left(\Phi ^{-1}(\bar{p})+\bar{a}\right)}, \quad A = 1 - ke^{\frac{\bar{b}^2}{2}}.$$
We choose the following parameterization: $\bar{p}=0.3$, $\bar{a}=1.6\sigma_{\rho}$, $\bar{b}=0.8 \sigma_{\rho}$.

\end{enumerate}

Figure \ref{figure:distortion} displays the three types of probability weighting functions under our parameterization.

\begin{figure}[htbp]
\centering

\includegraphics[width=5in]{./Figures/distortion.eps}

\caption{Three types of probability weighting functions}\label{figure:distortion}
\end{figure}

Figures \ref{figure:no_weighting}, \ref{figure:power}, and \ref{figure:JZ} display the optimal terminal wealth to \eqref{paramodel} and \eqref{paramodel1} as functions of $\rho$ when the weighting function is given by the identity weighting function, the power weighting function, and the Jin-Zhou weighting function, respectively, in different scenarios.

\bigskip

\begin{figure}[htbp]
\centering
\begin{minipage}[t]{2in}
\centering
\centerline{$c=0.1$ and $g=0$}
\includegraphics[width=2.2in]{./Figures/no_distortion/case1.eps}
\end{minipage}
\begin{minipage}[t]{2in}
\centering
\centerline{$c=0.1$ and $g=0.05$}
\includegraphics[width=2.2in]{./Figures/no_distortion/case2.eps}
\end{minipage}
\begin{minipage}[t]{2in}
\centering
\centerline{$c=0.1$ and $g=-0.05$}
\includegraphics[width=2.2in]{./Figures/no_distortion/case3.eps}
\end{minipage}

\quad

\quad

\begin{minipage}[t]{2in}
\centering
\centerline{$c=0.2$ and $g=0$}
\includegraphics[width=2.2in]{./Figures/no_distortion/case4.eps}
\end{minipage}
\begin{minipage}[t]{2in}
\centering
\centerline{$c=0.2$ and $g=0.05$}
\includegraphics[width=2.2in]{./Figures/no_distortion/case5.eps}
\end{minipage}
\begin{minipage}[t]{2in}
\centering
\centerline{$c=0.2$ and $g=-0.05$}
\includegraphics[width=2.2in]{./Figures/no_distortion/case6.eps}
\end{minipage}

\quad

\quad

\begin{minipage}[t]{2in}
\centering
\centerline{$c=0.3$ and $g=0$}
\includegraphics[width=2.2in]{./Figures/no_distortion/case7.eps}
\end{minipage}
\begin{minipage}[t]{2in}
\centering
\centerline{$c=0.3$ and $g=0.05$}
\includegraphics[width=2.2in]{./Figures/no_distortion/case8.eps}
\end{minipage}
\begin{minipage}[t]{2in}
\centering
\centerline{$c=0.3$ and $g=-0.05$}
\includegraphics[width=2.2in]{./Figures/no_distortion/case9.eps}
\end{minipage}
\caption{Optimal terminal wealth as functions of $\rho$ under the identity weighting function. }\label{figure:no_weighting}
\end{figure}

\begin{figure}[htbp]
\centering
\begin{minipage}[t]{2in}
\centering
\centerline{$c=0.1$ and $g=0$}
\includegraphics[width=2.2in]{./Figures/power/case1.eps}
\end{minipage}
\begin{minipage}[t]{2in}
\centering
\centerline{$c=0.1$ and $g=0.05$}
\includegraphics[width=2.2in]{./Figures/power/case2.eps}
\end{minipage}
\begin{minipage}[t]{2in}
\centering
\centerline{$c=0.1$ and $g=-0.05$}
\includegraphics[width=2.2in]{./Figures/power/case3.eps}
\end{minipage}

\quad

\quad

\begin{minipage}[t]{2in}
\centering
\centerline{$c=0.2$ and $g=0$}
\includegraphics[width=2.2in]{./Figures/power/case4.eps}
\end{minipage}
\begin{minipage}[t]{2in}
\centering
\centerline{$c=0.2$ and $g=0.05$}
\includegraphics[width=2.2in]{./Figures/power/case5.eps}
\end{minipage}
\begin{minipage}[t]{2in}
\centering
\centerline{$c=0.2$ and $g=-0.05$}
\includegraphics[width=2.2in]{./Figures/power/case6.eps}
\end{minipage}

\quad

\quad

\begin{minipage}[t]{2in}
\centering
\centerline{$c=0.3$ and $g=0$}
\includegraphics[width=2.2in]{./Figures/power/case7.eps}
\end{minipage}
\begin{minipage}[t]{2in}
\centering
\centerline{$c=0.3$ and $g=0.05$}
\includegraphics[width=2.2in]{./Figures/power/case8.eps}
\end{minipage}
\begin{minipage}[t]{2in}
\centering
\centerline{$c=0.3$ and $g=-0.05$}
\includegraphics[width=2.2in]{./Figures/power/case9.eps}
\end{minipage}
\caption{Optimal terminal wealth as functions of $\rho$ under the power weighting function. }\label{figure:power}
\end{figure}

\begin{figure}[htbp]
\centering
\begin{minipage}[t]{2in}
\centering
\centerline{$c=0.1$ and $g=0$}
\includegraphics[width=2.2in]{./Figures/JinZhou/case1.eps}
\end{minipage}
\begin{minipage}[t]{2in}
\centering
\centerline{$c=0.1$ and $g=0.05$}
\includegraphics[width=2.2in]{./Figures/JinZhou/case2.eps}
\end{minipage}
\begin{minipage}[t]{2in}
\centering
\centerline{$c=0.1$ and $g=-0.05$}
\includegraphics[width=2.2in]{./Figures/JinZhou/case3.eps}
\end{minipage}

\quad

\quad

\begin{minipage}[t]{2in}
\centering
\centerline{$c=0.2$ and $g=0$}
\includegraphics[width=2.2in]{./Figures/JinZhou/case4.eps}
\end{minipage}
\begin{minipage}[t]{2in}
\centering
\centerline{$c=0.2$ and $g=0.05$}
\includegraphics[width=2.2in]{./Figures/JinZhou/case5.eps}
\end{minipage}
\begin{minipage}[t]{2in}
\centering
\centerline{$c=0.2$ and $g=-0.05$}
\includegraphics[width=2.2in]{./Figures/JinZhou/case6.eps}
\end{minipage}

\quad

\quad

\begin{minipage}[t]{2in}
\centering
\centerline{$c=0.3$ and $g=0$}
\includegraphics[width=2.2in]{./Figures/JinZhou/case7.eps}
\end{minipage}
\begin{minipage}[t]{2in}
\centering
\centerline{$c=0.3$ and $g=0.05$}
\includegraphics[width=2.2in]{./Figures/JinZhou/case8.eps}
\end{minipage}
\begin{minipage}[t]{2in}
\centering
\centerline{$c=0.3$ and $g=-0.05$}
\includegraphics[width=2.2in]{./Figures/JinZhou/case9.eps}
\end{minipage}
\caption{Optimal terminal wealth as functions of $\rho$ under the Jin-Zhou weighting function. }\label{figure:JZ}
\end{figure}

In all scenarios under the three types of probability weighting functions, Figures \ref{figure:no_weighting}, \ref{figure:power} and \ref{figure:JZ} show that our terminal wealth is higher than or at least the same as Zhang et al.'s in the intermediate and bad states (i.e. $\rho$ is not very small). By contrast, our one is lower than Zhang et al.'s in the good states (i.e. $\rho$ is very small). From an economic viewpoint, the performance of our terminal wealth is more stable than Zhang et al.'s, and it is less sensitive to the pricing kernel (i.e. the curve is flatter). It suggests that our strategy is more suitable for loss-sensitive investors than Zhang et al.'s, whereas their strategy is more suitable for those investors betting for extremely good states.

\blue{
Once we have the expression for the optimal terminal wealth, we can solve for the time-$t$ wealth via $$X(t) = \frac{1}{\rho(t)}\BE{\rho(T) X(T)\;\big|\;\mathcal{F}_t},$$
which can be expressed as a deterministic function of $t$ and $\rho (t)$, say $\Psi(t,\rho(t))$. It is then straightforward to verify that the optimal portfolio weight is given by
$$\pi (t) = - \Psi_{x}(t,\rho(t)) \cdot \frac{\rho (t)}{X(t)} \left( \sigma (t)' \right)^{-1} \theta. $$
It is well-known that $\left( \sigma (t)'
 \right)^{-1} \theta$ is the (instantaneous) mean-variance portfolio, or the GOP, i.e., it is optimal for an investor maximizing the logarithm utility of terminal wealth. Therefore, $-\Psi_{x}(t,\rho(t)) \cdot {\rho (t)}/{X(t)}$ represents the risk exposure relative to the GOP.

Figure \ref{figure:no_weighting_strategy} presents the optimal risk exposures as functions of $\rho (t)$ when $t=0.5$ and the weighting function is given by the identity weighting function.\footnote{Optimal risk exposures under other distortion functions can be obtained in a similar fashion.} As $\rho (t)$ increases, the optimal risk exposures in both our model and Zhang et al.'s first decrease, then increase and finally decrease. In all scenarios, when $\rho (t)$ is small, the risk exposure in our model is significantly lower than that in Zhang et al.'s. When $g$ is small, the risk exposure in our model is slightly higher than Zhang et al.'s when $\rho (t)$ takes intermediate values.

\begin{figure}[htbp]
\centering
\begin{minipage}[t]{2in}
	\centering
	\centerline{$c=0.1$ and $g=0$}
	\includegraphics[width=2.2in]{./Figures/no_distortion/portfolio/case1.eps}
\end{minipage}
\begin{minipage}[t]{2in}
	\centering
	\centerline{$c=0.1$ and $g=0.05$}
	\includegraphics[width=2.2in]{./Figures/no_distortion/portfolio/case2.eps}
\end{minipage}
\begin{minipage}[t]{2in}
	\centering
	\centerline{$c=0.1$ and $g=-0.05$}
	\includegraphics[width=2.2in]{./Figures/no_distortion/portfolio/case3.eps}
\end{minipage}

\quad

\quad

\begin{minipage}[t]{2in}
	\centering
	\centerline{$c=0.2$ and $g=0$}
	\includegraphics[width=2.2in]{./Figures/no_distortion/portfolio/case4.eps}
\end{minipage}
\begin{minipage}[t]{2in}
	\centering
	\centerline{$c=0.2$ and $g=0.05$}
	\includegraphics[width=2.2in]{./Figures/no_distortion/portfolio/case5.eps}
\end{minipage}
\begin{minipage}[t]{2in}
	\centering
	\centerline{$c=0.2$ and $g=-0.05$}
	\includegraphics[width=2.2in]{./Figures/no_distortion/portfolio/case6.eps}
\end{minipage}

\quad

\quad

\begin{minipage}[t]{2in}
	\centering
	\centerline{$c=0.3$ and $g=0$}
	\includegraphics[width=2.2in]{./Figures/no_distortion/portfolio/case7.eps}
\end{minipage}
\begin{minipage}[t]{2in}
	\centering
	\centerline{$c=0.3$ and $g=0.05$}
	\includegraphics[width=2.2in]{./Figures/no_distortion/portfolio/case8.eps}
\end{minipage}
\begin{minipage}[t]{2in}
	\centering
	\centerline{$c=0.3$ and $g=-0.05$}
	\includegraphics[width=2.2in]{./Figures/no_distortion/portfolio/case9.eps}
\end{minipage}
\caption{Optimal exposures as functions of $\rho (t)$ under the identity weighting function; $t=0.5$.}\label{figure:no_weighting_strategy}
\end{figure}

}

\newpage
\section*{Appendix}