This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

The vanishing discount problem for monotone systems of Hamilton-Jacobi equations:
a counterexample to the full convergence

Hitoshi Ishii Institute for Mathematics and Computer Science
Tsuda University
2-1-1 Tsuda, Kodaira, Tokyo, 187-8577 Japan.
hitoshi.ishii@waseda.jp
Abstract.

In recent years there has been intense interest in the vanishing discount problem for Hamilton-Jacobi equations. In the case of the scalar equation, B. Ziliotto has recently given an example of the Hamilton-Jacobi equation having non-convex Hamiltonian in the gradient variable, for which the full convergence of the solutions does not hold as the discount factor tends to zero. We give here an explicit example of nonlinear monotone systems of Hamilton-Jacobi equations having convex Hamiltonians in the gradient variable, for which the full convergence of the solutions fails as the discount factor goes to zero.

Key words and phrases:
systems of Hamilton-Jacobi equations, vanishing discount, full convergence
2010 Mathematics Subject Classification:
35B40, 35D40, 35F50, 49L25
This paper is dedicated to Neil Trudinger on the occasion of his 80th birthday.

1. Introduction

We consider the system of Hamilton-Jacobi equations

(1) {λu1(x)+H1(Du1(x))+B1(u1(x),u2(x))=0 in 𝕋n,λu2(x)+H2(Du2(x))+B2(u1(x),u2(x))=0 in 𝕋n,\begin{cases}\lambda u_{1}(x)+H_{1}(Du_{1}(x))+B_{1}(u_{1}(x),u_{2}(x))=0\ &\text{ in }\mathbb{T}^{n},\\ \lambda u_{2}(x)+H_{2}(Du_{2}(x))+B_{2}(u_{1}(x),u_{2}(x))=0\ &\text{ in }\mathbb{T}^{n},\end{cases}

where λ>0\lambda>0 is a given constant, the functions Hi:nH_{i}:\mathbb{R}^{n}\to\mathbb{R} and Bi:2RB_{i}:\mathbb{R}^{2}\to R, with i=1,2i=1,2, are given continuous functions, and 𝕋n\mathbb{T}^{n} denotes the nn-dimensional flat torus n/n\mathbb{R}^{n}/\mathbb{Z}^{n}.

In a recent paper [IJ], the authors have investigated the vanishing discount problem for a nonlinear monotone system of Hamilton-Jacobi equations

(2) {λu1(x)+G1(x,Du1(x),u1(x),u2(x),,um(x))=0 in 𝕋n,λum(x)+Gm(x,Dum(x),u1(x),u2(x),,um(x))=0 in 𝕋n,\begin{cases}\lambda u_{1}(x)+G_{1}(x,Du_{1}(x),u_{1}(x),u_{2}(x),\ldots,u_{m}(x))=0\ &\text{ in }\mathbb{T}^{n},\\ \phantom{\lambda u_{1}(x)+G_{1}(x,Du_{1}(x),u_{1}(x),u_{2}(x)}\vdots&\\ \lambda u_{m}(x)+G_{m}(x,Du_{m}(x),u_{1}(x),u_{2}(x),\ldots,u_{m}(x))=0\ &\text{ in }\mathbb{T}^{n},\end{cases}

and established under some hypotheses on the GiC(𝕋n×n×m)G_{i}\in C(\mathbb{T}^{n}\times\mathbb{R}^{n}\times\mathbb{R}^{m}) that, when uλ=(u1λ,,umλ)C(𝕋n)mu^{\lambda}=(u_{1}^{\lambda},\ldots,u_{m}^{\lambda})\in C(\mathbb{T}^{n})^{m} denoting the (viscosity) solution of (2), the whole family {uλ}λ>0\{u^{\lambda}\}_{\lambda>0} converges in C(𝕋n)mC(\mathbb{T}^{n})^{m} to some u0C(𝕋n)mu^{0}\in C(\mathbb{T}^{n})^{m} as λ0+\lambda\to 0+. The constant λ>0\lambda>0 in the above system is the so-called discount factor.

The hypotheses on the system are the convexity, coercivity, and monotonicity of the GiG_{i} as well as the solvability of (2), with λ=0\lambda=0. Here the convexity of GiG_{i} is meant that the functions n×m(p,u)Gi(x,p,u)\mathbb{R}^{n}\times\mathbb{R}^{m}\ni(p,u)\mapsto G_{i}(x,p,u) are convex. We refer to [IJ] for the precise statement of the hypotheses.

Prior to work [IJ], there have been many contributions to the question about the whole family convergence (in other words, the full convergence) under the vanishing discount, which we refer to [IJ, DZ2, DFIZ, IMT1, IMT2, IS, CCIZ] and the references therein.

In the case of the scalar equation, B. Ziliotto [Zi] has recently shown an example of the Hamilton-Jacobi equation having non-convex Hamiltonian in the gradient variable for which the full convergence does not hold. In Ziliotto’s approach, the first step is to find a system of two algebraic equations

(3) {λu+f(uv)=0,λv+g(vu)=0,\left\{\begin{aligned} &\lambda u+f(u-v)=0,\\ &\lambda v+g(v-u)=0,\end{aligned}\right.

with two unknowns u,vu,v\in\mathbb{R} and with a parameter λ>0\lambda>0 as the discount factor, for which the solutions (uλ,vλ)(u^{\lambda},v^{\lambda}) stay bounded and fail to fully converge as λ0+\lambda\to 0+. Here, an “algebraic” equation means not to be a functional equation. The second step is to interporate the two values uλu^{\lambda} and vλv^{\lambda} to get a function of x𝕋1x\in\mathbb{T}^{1} which satisfies a scalar non-convex Hamilton-Jacobi equation in 𝕋1\mathbb{T}^{1}.

In the first step above, Ziliotto constructs f,gf,g based on a game-theoretical and computational argument, and the formula for f,gf,g is of the minimax type and not quite explicit. In [IH], the author has reexamined the system given by Ziliotto, with a slight generality, as a counterexample for the full convergence in the vanishing discount.

Our purpose in this paper is to present a system (3), with an explicit formula for f,gf,g, for which the solution (uλ,vλ)(u^{\lambda},v^{\lambda}) does not fully converge to a single point in 2\mathbb{R}^{2}. A straightforward consequence is that (1), with B1(u1,u2)=f(u1u2)B_{1}(u_{1},u_{2})=f(u_{1}-u_{2}) and B2(u1,u2)=g(u2u1)B_{2}(u_{1},u_{2})=g(u_{2}-u_{1}), has a solution given by

(u1λ(x),u2λ(x))=(u,v) for x𝕋n,(u_{1}^{\lambda}(x),u_{2}^{\lambda}(x))=(u,v)\ \ \text{ for }x\in\mathbb{T}^{n},

under the assumption that Hi(x,0)=0H_{i}(x,0)=0 for all x𝕋nx\in\mathbb{T}^{n}, and therefore, gives an example of a discoutned system of Hamilton-Jacobi equations, the solution of which fails to satisfy the full convergence as the discount factor goes to zero.

The paper consists of two sections. This introduction is followed by Section 2, the final section, which is divided into three subsections. The main results are stated in the first subsection of Section 2, the functions f,gf,g, the key elements of (3), are contstructed in the second subsection, and the final subsection provides the proof of the Main results.

2. A system of algebraic equations and the main results

Our main focus is now the system

{λu+f(uv)=0λv+g(vu)=0,\begin{cases}\lambda u+f(u-v)=0&\\ \lambda v+g(v-u)=0,\end{cases}

where f,gC(,)f,g\in C(\mathbb{R},\mathbb{R}) are nondecreasing functions, to be constructed, and λ>0\lambda>0 is a constant, to be sent to zero.

We remark that, due to the monotonicity assumption on f,gf,g, the mapping (u,v)(f(uv),g(vu)),22(u,v)\mapsto(f(u-v),g(v-u)),\,\mathbb{R}^{2}\to\mathbb{R}^{2} is monotone. Recall that, by definition, a mapping (u,v)(B1(u,v),B2(u,v)),22(u,v)\mapsto(B_{1}(u,v),B_{2}(u,v)),\,\mathbb{R}^{2}\to\mathbb{R}^{2} is monotone if, whenever (u1,v1),(u2,v2)2(u_{1},v_{1}),(u_{2},v_{2})\in\mathbb{R}^{2} satisfy u1u2v1v2u_{1}-u_{2}\geq v_{1}-v_{2} (resp., v1v2u1u2v_{1}-v_{2}\geq u_{1}-u_{2}), we have B1(u1,v1)B1(u2,v2)B_{1}(u_{1},v_{1})\geq B_{1}(u_{2},v_{2}) (resp., B2(u1,v1)B2(u2,v2)B_{2}(u_{1},v_{1})\geq B_{2}(u_{2},v_{2}))

2.1. Main results

Our main results are stated as follows.

Theorem 1.

There exist two increasing functions f,gC(,)f,g\in C(\mathbb{R},\mathbb{R}) having the properties (a)–(c):

  1. (a)

    For any λ>0\lambda>0 there exists a unique solution (uλ,vλ)2(u_{\lambda},v_{\lambda})\in\mathbb{R}^{2} to (3),

  2. (b)

    the family of the solutions (uλ,vλ)(u_{\lambda},v_{\lambda}) to (3), with λ>0\lambda>0, is bounded in 2\mathbb{R}^{2},

  3. (c)

    the family {(uλ,vλ)}λ>0\{(u_{\lambda},v_{\lambda})\}_{\lambda>0} does not converge as λ0+\lambda\to 0+.

It should be noted that, as mentioned in the introduction, the above theorem has been somewhat implicitly established by Ziliotto [Zi]. In this note, we are interested in a simple and easy approach to finding functions f,gf,g having the properties (a)–(c) in Theorem 1.

The following is an immediate consequence of the above theorem.

Corollary 2.

Let HiC(n,)H_{i}\in C(\mathbb{R}^{n},\mathbb{R}), i=1,2i=1,2, satisfy H1(0)=H2(0)=0H_{1}(0)=H_{2}(0)=0. Let f,gC(,)f,g\in C(\mathbb{R},\mathbb{R}) be the functions given by Theorem 1, and set B1(u1,u2)=f(u1u2)B_{1}(u_{1},u_{2})=f(u_{1}-u_{2}) and B2(u1,u2)=g(u1u1)B_{2}(u_{1},u_{2})=g(u_{1}-u_{1}) for all (u1,u2)2(u_{1},u_{2})\in\mathbb{R}^{2}. For any λ>0\lambda>0, let (uλ,1,uλ,2)(u_{\lambda,1},u_{\lambda,2}) be the (viscosity) solution of (1). Then, the functions uλ,iu_{\lambda,i} are constants, the family of the points (uλ,1,uλ,2)(u_{\lambda,1},u_{\lambda,2}) in 2\mathbb{R}^{2} is bounded, and it does not converge as λ0+\lambda\to 0+.

Notice that the convexity of HiH_{i} in the above corollary is irrelevant, and, for example, one may take Hi(p)=|p|2H_{i}(p)=|p|^{2} for i𝕀i\in\mathbb{I}, which are convex functions.

We remark that a claim similar to Corollary 2 is valid when one replaces Hi(p)H_{i}(p) by degenerate elliptic operators Fi(x,p,M)F_{i}(x,p,M) as far as Fi(x,0,0)=0F_{i}(x,0,0)=0, where MM is the variable corresponding to the Hessian matrices of unknown functions. (See [CIL] for an overview on the viscosity solution approach to fully nonlinear degenerate elliptic equations.)

2.2. The functions f,gf,g

If f,gf,g are given and (u,v)2(u,v)\in\mathbb{R}^{2} is a solution of (2), then w:=uvw:=u-v satisfies

(3) λw+f(w)g(w)=0.\lambda w+f(w)-g(-w)=0.

Set

(4) h(r)=f(r)g(r) for r,h(r)=f(r)-g(-r)\ \ \text{ for }r\in\mathbb{R},

which defines a continuous and nondecreasing function on \mathbb{R}.

To build a triple of functions f,g,hf,g,h, we need to find two of them in view of the relation (4). We begin by defining function hh.

For this, we discuss a simple geometry on xyxy-plane as depicted in Fig. 1 below. Fix 0<k1<k20<k_{1}<k_{2}. The line y=12k2+k1(x+12)y=-\frac{1}{2}k_{2}+k_{1}(x+\frac{1}{2}) has slope k1k_{1} and crosses the lines x=1x=-1 and y=k2xy=k_{2}x at P:=(1,12(k1+k2))\mathrm{P}:=(-1,-\frac{1}{2}(k_{1}+k_{2})) and Q:=(12,12k2)\mathrm{Q}:=(-\frac{1}{2},-\frac{1}{2}k_{2}), respectively, while the line y=k2xy=k_{2}x meets the lines x=1x=-1 and x=12x=-\frac{1}{2} at R:=(1,k2)\mathrm{R}:=(-1,-k_{2}) and Q=(12,12k2)\mathrm{Q}=(-\frac{1}{2},-\frac{1}{2}k_{2}), respectively.

Choose k>0k^{*}>0 so that 12(k1+k2)<k<k2\frac{1}{2}(k_{1}+k_{2})<k^{*}<k_{2}. The line y=kxy=k^{*}x crosses the line y=12k2+k1(x+12)y=-\frac{1}{2}k_{2}+k_{1}(x+\frac{1}{2}) at a point S:=(x,y)\mathrm{S}:=(x^{*},y^{*}) in the open line segment between the points P=(12,12(k1+k2))\mathrm{P}=(-\frac{1}{2},-\frac{1}{2}(k_{1}+k_{2})) and Q=(12,12k2)\mathrm{Q}=(-\frac{1}{2},-\frac{1}{2}k_{2}). The line connecting R=(1,k2)\mathrm{R}=(-1,-k_{2}) and S=(x,y)\mathrm{S}=(x^{*},y^{*}) can be represented by y=k2+k+(x+1)y=-k_{2}+k^{+}(x+1), with k+:=y+k2x+1>k2k^{+}:=\frac{y^{*}+k_{2}}{x^{*}+1}>k_{2}.

Oyyxx1-11/2-1/2y=k2xy=k_{2}xy=12k2+k1(x+12)y=-\frac{1}{2}k_{2}+k_{1}(x+\frac{1}{2})y=kxy=k^{*}xPQSR
Figure 1. Graph of ψ\psi.

We set

ψ(x)={k2x for x(,1][1/2,),min{k2+k+(x+1),12k2+k1(x+12)} for x(1,12).\psi(x)=\begin{cases}k_{2}x\qquad\qquad\text{ for }x\in(-\infty,-1]\cup[-1/2,\infty),&\\ \min\{-k_{2}+k^{+}(x+1),-\frac{1}{2}k_{2}+k_{1}(x+\frac{1}{2})\}\ \ \text{ for }x\in(-1,-\frac{1}{2}).&\end{cases}

It is clear that ψC()\psi\in C(\mathbb{R}) and increasing on \mathbb{R}. The building blocks of the graph y=ψ(x)y=\psi(x) are three lines whose slpoes are k1<k2<k+k_{1}<k_{2}<k^{+}. Hence, if x1>x2x_{1}>x_{2}, then ψ(x1)ψ(x2)k1(x1x2)\psi(x_{1})-\psi(x_{2})\geq k_{1}(x_{1}-x_{2}), that is, the function xψ(x)k1xx\mapsto\psi(x)-k_{1}x is nondecreasing on \mathbb{R}.

Next, we set for jj\in\mathbb{N},

ψj(x)=2jψ(2jx) for x.\psi_{j}(x)=2^{-j}\psi(2^{j}x)\ \ \text{ for }x\in\mathbb{R}.

It is clear that for all jj\in\mathbb{N}, ψjC()\psi_{j}\in C(\mathbb{R}), the function xψj(x)k1xx\mapsto\psi_{j}(x)-k_{1}x is nondecreasing on \mathbb{R}, and

ψj(x){>k2x for all x(2j,2j1),=k2x otherwise.\psi_{j}(x)\begin{cases}>k_{2}x\ \ &\text{ for all }x\in(-2^{-j},-2^{-j-1}),\\ =k_{2}x\ \ &\text{ otherwise}.\end{cases}

We set

η(x)=maxjψj(x) for x.\eta(x)=\max_{j\in\mathbb{N}}\psi_{j}(x)\ \ \text{ for }x\in\mathbb{R}.

It is clear that ηC()\eta\in C(\mathbb{R}) and xη(x)k1xx\mapsto\eta(x)-k_{1}x is nondecreasing on \mathbb{R}. Moreover, we see that

η(x)=k2x for all x(,12][0,),\eta(x)=k_{2}x\ \ \text{ for all }x\in(-\infty,-\tfrac{1}{2}]\cup[0,\infty),

and that if 2j<x<2j1-2^{-j}<x<-2^{-j-1} and jj\in\mathbb{N},

η(x)=ψj(x)>k2x.\eta(x)=\psi_{j}(x)>k_{2}x.

Note that the point S=(x,y)\mathrm{S}=(x^{*},y^{*}) is on the graph y=ψ(x)y=\psi(x) and, hence, that for any jj\in\mathbb{N}, the point (2jx,2jy)(2^{-j}x^{*},2^{-j}y^{*}) is on the graph y=η(x)y=\eta(x). Similarly, since the point S=(x,y)\mathrm{S}=(x^{*},y^{*}) is on the graph y=kxy=k^{*}x and for any jj\in\mathbb{N}, the point (2jx,2jy)(2^{-j}x^{*},2^{-j}y^{*}) is on the graph y=kxy=k^{*}x. Also, for any jj\in\mathbb{N}, the point (2j,k22j)(-2^{-j},-k_{2}2^{-j}) lies on the graphs y=η(x)y=\eta(x) and y=k2xy=k_{2}x.

Fix any d1d\geq 1 and define hC()h\in C(\mathbb{R}) by

h(x)=η(xd).h(x)=\eta(x-d).

For the function hh defined above, we consider the problem

(5) λz+h(z)=0.\lambda z+h(z)=0.
Lemma 3.

For any λ0\lambda\geq 0, there exists a unique solution zλz_{\lambda}\in\mathbb{R} of (5).

Proof.

Fix λ0\lambda\geq 0. The function xh(x)+λxx\mapsto h(x)+\lambda x is increasing on \mathbb{R} and satisfies

limx(h(x)+λx)= and limx(h(x)+λx)=.\lim_{x\to\infty}(h(x)+\lambda x)=\infty\ \ \text{ and }\ \ \lim_{x\to-\infty}(h(x)+\lambda x)=-\infty.

Hence, there is a unique solution of (5). ∎

For any λ0\lambda\geq 0, we denote by zλz_{\lambda} the unique solution of (5). Since h(d)=0h(d)=0, it is clear that z0=dz_{0}=d.

Lemma 4.

Let λ>0\lambda>0 and k>0k>0. Let (z,w)2(z,w)\in\mathbb{R}^{2} be the point of the intersection of two lines y=λxy=-\lambda x and y=k(xd)y=k(x-d). Then

z=kdk+λ.z=\frac{kd}{k+\lambda}.
Proof.

By the assumption, we have

w=λz=k(zd),w=-\lambda z=k(z-d),

and hence,

z=kdk+λ.z=\frac{kd}{k+\lambda}.\qed
Lemma 5.

There are sequences {μj}\{\mu_{j}\} and {νj}\{\nu_{j}\} of positive numbers converging to zero such that

zμj=k2dk2+μj and zνj=kdk+νj.z_{\mu_{j}}=\frac{k_{2}d}{k_{2}+\mu_{j}}\ \ \text{ and }\ \ z_{\nu_{j}}=\frac{k^{*}d}{k^{*}+\nu_{j}}.
Proof.

Let jj\in\mathbb{N}. Since (2j,k22j)(-2^{-j},-k_{2}2^{-j}) is on the intersection of the graphs y=k2xy=k_{2}x and y=η(x)y=\eta(x), it follows that (2j+d,k22j)(-2^{-j}+d,-k_{2}2^{-j}) is on the intersection of the graphs y=k2(xd)y=k_{2}(x-d) and y=h(x)y=h(x). Set

(6) μj=k22jd2j,\mu_{j}=\frac{k_{2}2^{-j}}{d-2^{-j}},

and note that μj>0\mu_{j}>0 and that

μj(d2j)=k22j,-\mu_{j}(d-2^{-j})=-k_{2}2^{-j},

which says that the point (d2j,k22j)(d-2^{-j},-k_{2}2^{-j}) is on the line y=μjxy=-\mu_{j}x. Combining the above with

k22j=h(d2j)-k_{2}2^{-j}=h(d-2^{-j})

shows that d2jd-2^{-j} is the unique solution of (5). Also, since (d2j,μj(d2j))=(d2j,k22j)(d-2^{-j},-\mu_{j}(d-2^{-j}))=(d-2^{-j},-k_{2}2^{-j}) is on the line y=k2(xd)y=k_{2}(x-d), we find by Lemma 4 that

zμj=k2dk2+μj.z_{\mu_{j}}=\frac{k_{2}d}{k_{2}+\mu_{j}}.

Similarly, since (2jx,2jy)(2^{-j}x^{*},2^{-j}y^{*}) is on the intersection of the graphs y=kxy=k^{*}x and y=η(x)y=\eta(x), we deduce that if we set

(7) νj:=2jyd+2jx=2j|y|d2j|x|,\nu_{j}:=-\frac{2^{-j}y^{*}}{d+2^{-j}x^{*}}=\frac{2^{-j}|y^{*}|}{d-2^{-j}|x^{*}|},

then

zνj=kdk+νj.z_{\nu_{j}}=\frac{k^{*}d}{k^{*}+\nu_{j}}.

It is ovbvious by (6) and (7) that the sequences {μj}j\{\mu_{j}\}_{j\in\mathbb{N}} and {νj}j\{\nu_{j}\}_{j\in\mathbb{N}} are decreasing and converge to zero. ∎

We fix k0(0,k1)k_{0}\in(0,k_{1}) and define f,gC()f,g\in C(\mathbb{R}) by f(x)=k0(xd)f(x)=k_{0}(x-d) and

g(x)=f(x)h(x).g(x)=f(-x)-h(-x).

It is easily checked that g(x)(k1k0)xg(x)-(k_{1}-k_{0})x is nondecreasing on \mathbb{R}, which implies that gg is increasing on \mathbb{R}, and that h(x)=f(x)g(x)h(x)=f(x)-g(-x) for all xx\in\mathbb{R}. We note that

(8) f(d)=h(d)=g(d)=0.f(d)=h(d)=g(-d)=0.

2.3. Proof of the main results

We fix f,g,hf,g,h as above, and consider the system (2).

Lemma 6.

Let λ>0\lambda>0. There exists a unique solution of (3).

The validity of the above lemma is well-known, but for the reader’s convenience, we provide a proof of the lemma above.

Proof.

By the choice of f,gf,g, the functions f,gf,g are nondecreasing on \mathbb{R}. We show first the comparison claim: if (u1,v1),(u2v2)2(u_{1},v_{1}),(u_{2}v_{2})\in\mathbb{R}^{2} satisfy

(9) λu1+f(u1v1)0,λv1+g(v1u1)0,\displaystyle\lambda u_{1}+f(u_{1}-v_{1})\leq 0,\quad\lambda v_{1}+g(v_{1}-u_{1})\leq 0,
(10) λu2+f(u2v2)0,λv2+g(v2u2)0,\displaystyle\lambda u_{2}+f(u_{2}-v_{2})\geq 0,\quad\lambda v_{2}+g(v_{2}-u_{2})\geq 0,

then u1u2u_{1}\leq u_{2} and v1v2v_{1}\leq v_{2}. Indeed, contrary to this, we suppose that max{u1u2,v1v2}>0\max\{u_{1}-u_{2},v_{1}-v_{2}\}>0. For instance, if max{u1u2,v1v2}=u1u2\max\{u_{1}-u_{2},v_{1}-v_{2}\}=u_{1}-u_{2}, then we have u1v1u2v2u_{1}-v_{1}\geq u_{2}-v_{2} and u1>u2u_{1}>u_{2}, and moreover

0λu1+f(u1v1)λu1+f(u2v2)>λu2+f(u2v2),0\geq\lambda u_{1}+f(u_{1}-v_{1})\geq\lambda u_{1}+f(u_{2}-v_{2})>\lambda u_{2}+f(u_{2}-v_{2}),

yielding a contradiction. The other case when max{u1u2,v1v2}=v1v2\max\{u_{1}-u_{2},v_{1}-v_{2}\}=v_{1}-v_{2}, we find a contradiction, 0>λv2+g(v2u2)0>\lambda v_{2}+g(v_{2}-u_{2}), proving the comparison.

From the comparison claim, the uniqueness of the solutions of (3) follows readily.

Next, we may choose a constant C>0C>0 so large that (u1,v1)=(C,C)(u_{1},v_{1})=(-C,-C) and (u2,v2)=(C,C)(u_{2},v_{2})=(C,C) satisfy (9) and (10), respectively. We write SS for the set of all (u1,u2)2(u_{1},u_{2})\in\mathbb{R}^{2} such that (9) hold. Note that (C,C)S(-C,-C)\in S and that for any (u,v)S(u,v)\in S, uCu\leq C and vCv\leq C. We set

u\displaystyle u^{*} =sup{u:(u,v)S for some v},\displaystyle=\sup\{u\,:\,(u,v)\in S\ \text{ for some }v\},
v\displaystyle v^{*} =sup{v:(u,v)S for some u}.\displaystyle=\sup\{v\,:\,(u,v)\in S\ \text{ for some }u\}.

It follows that Cu,vC-C\leq u^{*},v^{*}\leq C. We can choose sequences

{(un1,vn1)}n,{(un2,vn2)}nS\{(u_{n}^{1},v_{n}^{1})\}_{n\in\mathbb{N}},\,\{(u_{n}^{2},v_{n}^{2})\}_{n\in\mathbb{N}}\subset S

such that {un1},{vn2}\{u_{n}^{1}\},\{v_{n}^{2}\} are nondecreasing,

limnun1=u and limnvn2=v.\lim_{n\to\infty}u_{n}^{1}=u^{*}\ \ \text{ and }\ \ \lim_{n\to\infty}v_{n}^{2}=v^{*}.

Observe that for all nn\in\mathbb{N}, un2uu_{n}^{2}\leq u^{*}, vn1vv_{n}^{1}\leq v^{*}, and

0λun1+f(un1vn1)λun1+f(un1v),0\geq\lambda u_{n}^{1}+f(u_{n}^{1}-v_{n}^{1})\geq\lambda u_{n}^{1}+f(u_{n}^{1}-v^{*}),

which yields, in the limit as nn\to\infty,

0λu+f(uv).0\geq\lambda u^{*}+f(u^{*}-v^{*}).

Similarly, we obtain 0λv+g(vu)0\geq\lambda v^{*}+g(v^{*}-u^{*}). Hence, we find that (u,v)S(u^{*},v^{*})\in S.

We claim that (u,v)(u^{*},v^{*}) is a solution of (3). Otherwise, we have

0>λu+f(uv) or  0>λv+g(vu).0>\lambda u^{*}+f(u^{*}-v^{*})\ \ \text{ or }\ \ 0>\lambda v^{*}+g(v^{*}-u^{*}).

For instance, if the former of the above inequalities holds, we can choose ε>0\varepsilon>0, by the continuity of ff, so that

0>λ(u+ε)+f(u+εv).0>\lambda(u^{*}+\varepsilon)+f(u^{*}+\varepsilon-v^{*}).

Since (u,v)S(u^{*},v^{*})\in S, we have

0λv+g(vu)λv+g(vuε).0\geq\lambda v^{*}+g(v^{*}-u^{*})\geq\lambda v^{*}+g(v^{*}-u^{*}-\varepsilon).

Accordingly, we find that (u+ε,v)S(u^{*}+\varepsilon,v^{*})\in S, which contradicts the definition of uu^{*}. Similarly, if 0>λv+g(vu)0>\lambda v^{*}+g(v^{*}-u^{*}), then we can choose δ>0\delta>0 so that (u,v+δ)S(u^{*},v^{*}+\delta)\in S, which is a contradiction. Thus, we conclude that (u,v)(u^{*},v^{*}) is a solution of (3). ∎

Theorem 7.

For any λ>0\lambda>0, let (uλ,vλ)(u_{\lambda},v_{\lambda}) denote the unique solution of (3). Let {μj},{νj}\{\mu_{j}\},\{\nu_{j}\} be the sequences of positive numbers from Lemma 5. Then

limjuμj=k0dk2 and limjuνj=k0dk.\lim_{j\to\infty}u_{\mu_{j}}=\frac{k_{0}d}{k_{2}}\ \ \text{ and }\ \ \lim_{j\to\infty}u_{\nu_{j}}=\frac{k_{0}d}{k^{*}}.

In particular,

lim infλ0uλk0dk2<k0dklim supλ0uλ.\liminf_{\lambda\to 0}u_{\lambda}\leq\frac{k_{0}d}{k_{2}}<\frac{k_{0}d}{k^{*}}\leq\limsup_{\lambda\to 0}u_{\lambda}.

With our choice of f,gf,g, the family of solutions (uλ,vλ)(u_{\lambda},v_{\lambda}) of (2), with λ>0\lambda>0, does not converge as λ0\lambda\to 0.

Proof.

If we set zλ=uλvλz_{\lambda}=u_{\lambda}-v_{\lambda}, then zλz_{\lambda} satisfies (5). By Lemma 5, we find that

zμj=k2dk2+μj and zνj=kdk+νj.z_{\mu_{j}}=\frac{k_{2}d}{k_{2}+\mu_{j}}\ \ \text{ and }\ \ z_{\nu_{j}}=\frac{k^{*}d}{k^{*}+\nu_{j}}.

Since uλu_{\lambda} satisfies

0=λuλ+f(zλ)=λuλ+k0(zλd),0=\lambda u_{\lambda}+f(z_{\lambda})=\lambda u_{\lambda}+k_{0}(z_{\lambda}-d),

we find that

uμj=k0(zμjd)μj=k0dμj(k2k2+μj1)=k0dμjμjk2+μj=k0dk2+μj,u_{\mu_{j}}=-\frac{k_{0}(z_{\mu_{j}}-d)}{\mu_{j}}=-\frac{k_{0}d}{\mu_{j}}\left(\frac{k_{2}}{k_{2}+\mu_{j}}-1\right)=-\frac{k_{0}d}{\mu_{j}}\frac{-\mu_{j}}{k_{2}+\mu_{j}}=\frac{k_{0}d}{k_{2}+\mu_{j}},

which shows that

limjuμj=k0dk2.\lim_{j\to\infty}u_{\mu_{j}}=\frac{k_{0}d}{k_{2}}.

A parallel computation shows that

limjuνj=k0dk.\lim_{j\to\infty}u_{\nu_{j}}=\frac{k_{0}d}{k^{*}}.

Recalling that 0<k<k20<k^{*}<k_{2}, we conclude that

lim infλ0uλk0dk2<k0dklim supλ0uλ.\liminf_{\lambda\to 0}u_{\lambda}\leq\frac{k_{0}d}{k_{2}}<\frac{k_{0}d}{k^{*}}\leq\limsup_{\lambda\to 0}u_{\lambda}.\qed

We remark that, since

limλ0zλ=d and vλ=uλzλ,\lim_{\lambda\to 0}z_{\lambda}=d\ \ \text{ and }\ \ v_{\lambda}=u_{\lambda}-z_{\lambda},
limjvμj=k0dk2d and limjvνj=k0dkd.\lim_{j\to\infty}v_{\mu_{j}}=\frac{k_{0}d}{k_{2}}-d\ \ \text{ and }\ \ \lim_{j\to\infty}v_{\nu_{j}}=\frac{k_{0}d}{k^{*}}-d.

We give the proof of Theorem 1.

Proof of Theorem 1.

Assertions (a) and (c) are consequences of Lemma 6 and Theorem 7, respectively.

Recall (8). That is, we have f(d)=h(d)=g(d)=0f(d)=h(d)=g(-d)=0. Setting (u2,v2)=(d,0)(u_{2},v_{2})=(d,0), we compute that for any λ>0\lambda>0,

λu2+f(u2v2)>f(d)=0 and λv2+g(v2u2)=g(d)=0.\lambda u_{2}+f(u_{2}-v_{2})>f(d)=0\ \ \text{ and }\ \ \lambda v_{2}+g(v_{2}-u_{2})=g(-d)=0.

By the comparison claim, proved in the proof of Lemma 6, we find that uλdu_{\lambda}\leq d and vλ0v_{\lambda}\leq 0 for any λ>0\lambda>0. Simlarly, setting (u1,v1)=(0,d)(u_{1},v_{1})=(0,-d), we find that for any λ>0\lambda>0,

λu1+f(u1v1)=f(d)=0 and λv1+g(v1u1)g(v1u1)=g(d)=0,\lambda u_{1}+f(u_{1}-v_{1})=f(d)=0\ \ \text{ and }\ \ \lambda v_{1}+g(v_{1}-u_{1})\leq g(v_{1}-u_{1})=g(-d)=0,

which shows by the comparison claim that uλ0u_{\lambda}\geq 0 and vλdv_{\lambda}\geq-d for any λ>0\lambda>0. Thus, the sequence {(uλ,vλ)}λ>0\{(u_{\lambda},v_{\lambda})\}_{\lambda>0} is bounded in 2\mathbb{R}^{2}, which proves assertion (b). ∎

Proof of Corollary 2.

For any λ>0\lambda>0, we set (uλ,vλ)2(u_{\lambda},v_{\lambda})\in\mathbb{R}^{2} be the unique solution of (3). Since H1(0)=H2(0)=0H_{1}(0)=H_{2}(0)=0, it is clear that the constant function (uλ,1(x),uλ,2(x)):=(uλ,vλ)(u_{\lambda,1}(x),u_{\lambda,2}(x)):=(u_{\lambda},v_{\lambda}) is a classical solution of (1). By a classical uniqueness result (see, for instance, [IK, Theorem 4.7]), (uλ,1,uλ,2)(u_{\lambda,1},u_{\lambda,2}) is a unique viscosity solution of (1). The rest of the claims in Corollary 2 is an immediate consequence of Theorem 1. ∎

Some remarks are in order. (i) Following [Zi], we may use Theorem 7 as the primary cornerstone for building a scalar Hamilton-Jacobi equation, for which the vanishing discount problem fails to have the full convergence as the discount factor goes to zero.

(ii) In the construction of the functions f,gC(,)f,g\in C(\mathbb{R},\mathbb{R}) in Theorem 7, the author has chosen dd to satisfy d1d\geq 1, but one may choose any d>0d>0. In the process, the crucial step is to find the function h(x)=f(x)g(x)h(x)=f(x)-g(-x), with the properties: (a) the function xh(x)εxx\mapsto h(x)-\varepsilon x is nondecreasing on \mathbb{R} for some ε>0\varepsilon>0 and (b) the curve y=h(x)y=h(x), with x<dx<d, meets the lines y=p(xd)y=p(x-d) and y=q(xd)y=q(x-d), respectively, at PnP_{n} and QnQ_{n} for all nn\in\mathbb{N}, where p,q,dp,q,d are positive constants such that ε<p<q\varepsilon<p<q, and the sequences {Pn}n,{Qn}n\{P_{n}\}_{n\in\mathbb{N}},\,\{Q_{n}\}_{n\in\mathbb{N}} converge to the point (d,0)(d,0). Thus, a possible choice of hh among many other ways is the following. Define first η:\eta\,:\,\mathbb{R}\to\mathbb{R} by η(x)=x(sin(log|x|)+2)\eta(x)=x(\sin(\log|x|)+2) if x0x\not=0, and η(0)=0\eta(0)=0 (see Fig. 2). Fix d>0d>0 and set h(x)=η(xd)h(x)=\eta(x-d) for xx\in\mathbb{R}. Note that if x0x\not=0,

η(x)=sin(log|x|)+cos(log|x|)+2[22,2+2],\eta^{\prime}(x)=\sin(\log|x|)+\cos(\log|x|)+2\in[2-\sqrt{2},2+\sqrt{2}],

and that if we set xn=exp(2πn)x_{n}=-\exp(-2\pi n) and ξn=exp(2πn+π2)\xi_{n}=-\exp\left(-2\pi n+\frac{\pi}{2}\right), nn\in\mathbb{N}, then

η(xn)=2xn and η(ξn)=3ξn.\eta(x_{n})=2x_{n}\ \ \text{ and }\ \ \eta(\xi_{n})=3\xi_{n}.

The points Pn:=(xn+d,2xn)P_{n}:=(x_{n}+d,2x_{n}) are on the intersection of two curves y=h(x)y=h(x) and y=2(xd)y=2(x-d), while the points Qn:=(d+ξn,3ξn)Q_{n}:=(d+\xi_{n},3\xi_{n}) are on the intersection of y=h(x)y=h(x) and y=3(xd)y=3(x-d). Moreover, limPn=limQn=(d,0)\lim P_{n}=\lim Q_{n}=(d,0).

Oyyxx
Figure 2. Graph of η\eta (slightly deformed).

Acknowledgments

The author was supported in part by the JSPS Grants KAKENHI No. 16H03948, No. 20K03688, No. 20H01817, and No. 21H00717.

References