Processing math: 0%

数学物理学报, 2025, 45(1): 1-30

拟线性薛定谔方程组在有界区域上的正规化解

张倩,

清华大学数学科学系 北京 100084; 福建师范大学数学与统计学院 福州 350117

Normalized Solutions of the Quasilinear Schrödinger System in Bounded Domains

Zhang Qian,

Department of Mathematical Sciences, Tsinghua University, Beijing 100084; School of Mathematics and Statistics, Fujian Normal University, Fuzhou 350117

收稿日期: 2023-07-19   修回日期: 2024-09-9  

基金资助: 福建省高校数学学科联盟科研项目专项资金(2025SXLMQN04)

Received: 2023-07-19   Revised: 2024-09-9  

Fund supported: Fujian Alliance of Mathematics(2025SXLMQN04)

作者简介 About authors

张倩,E-mail:zhangqian9115@mail.tsinghua.edu.cn

摘要

该文关注以下非线性耦合方程组
{Δu1+ω1u112Δ(u21)u1=μ1|u1|p1u1+β|u2|p+12|u1|p32u1Δu2+ω2u212Δ(u22)u2=μ2|u2|p1u2+β|u1|p+12|u2|p32u2Ω|ui|2 dx=ρi,i=1,2,(u1,u2)H10(Ω;R2)
以及线性耦合方程组
{Δu1+ω1u112Δ(u21)u1=μ1|u1|p1u1+βu2Δu2+ω2u212Δ(u22)u2=μ2|u2|p1u2+βu1Ω|ui|2 dx=ρi,i=1,2,(u1,u2)H10(Ω;R2)
其中 ΩRN(N1) 是一个有界光滑区域,ωi, βR, μi, ρi>0, i=1,2. 而且, 若 p>1, N=1,2 且若 1<p, N\geqslant3. 应用变量替换, 一方面,证明了非线性耦合方程组正规化解的存在性和轨道稳定性, 以及当 \beta\rightarrow-\infty 时正规化解的极限行为. 另一方面, 应用极小化约束方法来获得线性耦合方程组的正规化解的存在性. 与之前的一些结果相比, 将现有结果扩展到了拟线性薛定谔方程组, 并获得了线性耦合情形下的正规化解.

关键词: 线性与非线性耦合; 有界区域; 变量替换; 正规化解; 极限行为.

Abstract

This paper is concerned with the following nonlinear coupled system
\left\{\begin{array}{l} -\Delta u_{1}+\omega_{1} u_{1}-\frac{1}{2} \Delta\left(u_{1}^{2}\right) u_{1}=\mu_{1}\left|u_{1}\right|^{p-1} u_{1}+\beta\left|u_{2}\right|^{\frac{p+1}{2}}\left|u_{1}\right|^{\frac{p-3}{2}} u_{1} \\ -\Delta u_{2}+\omega_{2} u_{2}-\frac{1}{2} \Delta\left(u_{2}^{2}\right) u_{2}=\mu_{2}\left|u_{2}\right|^{p-1} u_{2}+\beta\left|u_{1}\right|^{\frac{p+1}{2}}\left|u_{2}\right|^{\frac{p-3}{2}} u_{2} \\ \int_{\Omega}\left|u_{i}\right|^{2} \mathrm{~d} x=\rho_{i}, \quad i=1,2, \quad\left(u_{1}, u_{2}\right) \in H_{0}^{1}\left(\Omega ; \mathbb{R}^{2}\right) \end{array}\right.
and linear coupled system
\left\{\begin{array}{l} -\Delta u_{1}+\omega_{1} u_{1}-\frac{1}{2} \Delta\left(u_{1}^{2}\right) u_{1}=\mu_{1}\left|u_{1}\right|^{p-1} u_{1}+\beta u_{2} \\ -\Delta u_{2}+\omega_{2} u_{2}-\frac{1}{2} \Delta\left(u_{2}^{2}\right) u_{2}=\mu_{2}\left|u_{2}\right|^{p-1} u_{2}+\beta u_{1} \\ \int_{\Omega}\left|u_{i}\right|^{2} \mathrm{~d} x=\rho_{i}, \quad i=1,2, \quad\left(u_{1}, u_{2}\right) \in H_{0}^{1}\left(\Omega ; \mathbb{R}^{2}\right) \end{array}\right.
where \Omega\subset\mathbb R^N(N\geq1) is a bounded smooth domain, \omega_i,\ \beta\in\mathbb R, \mu_i,\ \rho_i>0,\ i=1,2. Moreover, p>1 if N=1,2 and 1<p\leqslant\frac{3N+2}{N-2} if N\geqslant3. Using change of variables, on the one hand, we prove the existence and stability of normalized solutions in nonlinear coupled system and the limiting behavior of normalized solutions as \beta\rightarrow -\infty. On the other hand, we apply the minimization constraint technique to obtain the existence of normalized solutions for linear coupled system. Compared with some previous results, we extend the existing results to the quasilinear Schrödinger system and also obtain normalized solutions for the linear coupling case.

Keywords: linear and nonlinear coupled; bounded domains; change of variables; normalized solution; limiting behavior.

PDF (779KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

张倩. 拟线性薛定谔方程组在有界区域上的正规化解[J]. 数学物理学报, 2025, 45(1): 1-30

Zhang Qian. Normalized Solutions of the Quasilinear Schrödinger System in Bounded Domains[J]. Acta Mathematica Scientia, 2025, 45(1): 1-30

1 引言

本章我们考虑以下拟线性耦合方程组

\left\{\begin{array}{l} -\Delta u_{1}+\omega_{1} u_{1}-\frac{1}{2} \Delta\left(u_{1}^{2}\right) u_{1}=\mu_{1}\left|u_{1}\right|^{p-1} u_{1}+\beta\left|u_{2}\right|^{\frac{p+1}{2}}\left|u_{1}\right|^{\frac{p-3}{2}} u_{1} \\ -\Delta u_{2}+\omega_{2} u_{2}-\frac{1}{2} \Delta\left(u_{2}^{2}\right) u_{2}=\mu_{2}\left|u_{2}\right|^{p-1} u_{2}+\beta\left|u_{1}\right|^{\frac{p+1}{2}}\left|u_{2}\right|^{\frac{p-3}{2}} u_{2} \\ \int_{\Omega}\left|u_{i}\right|^{2} \mathrm{~d} x=\rho_{i}, \quad i=1,2 \end{array}\right.
(1.1)

以及线性耦合方程组

\left\{\begin{array}{l} -\Delta u_{1}+\omega_{1} u_{1}-\frac{1}{2} \Delta\left(u_{1}^{2}\right) u_{1}=\mu_{1}\left|u_{1}\right|^{p-1} u_{1}+\beta u_{2} \\ -\Delta u_{2}+\omega_{2} u_{2}-\frac{1}{2} \Delta\left(u_{2}^{2}\right) u_{2}=\mu_{2}\left|u_{2}\right|^{p-1} u_{2}+\beta u_{1} \\ \int_{\Omega}\left|u_{i}\right|^{2} \mathrm{~d} x=\rho_{i}, \quad i=1,2, \quad\left(u_{1}, u_{2}\right) \in H_{0}^{1}\left(\Omega ; \mathbb{R}^{2}\right) \end{array}\right.

其中 (u_1,u_2)\in H_0^1(\Omega;\mathbb R^2), \Omega\subset\mathbb R^N(N\geq1) 是一个有界光滑区域, \omega_i,\ \beta\in\mathbb R, \mu_i,\ \rho_i>0, i=1,2. 而且, 当 N=1,2p>1 且当 N\geqslant31<p\leqslant\frac{3N+2}{N-2}. 方程组(1.1)与以下 Gross-Pitaevskii 方程组有关 (见文献[29]):

\left\{\begin{array}{l} \mathrm{i} \partial_{t} \Phi_{1}+\Delta \Phi_{1}+\Phi_{1}\left(\mu_{1}\left|\Phi_{1}\right|^{p-1}+\beta\left|\Phi_{2}\right|^{\frac{p+1}{2}}\left|\Phi_{1}\right|^{\frac{p-3}{2}}\right)=0,(t, x) \in \mathbb{R} \times \Omega \\ \mathrm{i} \partial_{t} \Phi_{2}+\Delta \Phi_{2}+\Phi_{2}\left(\mu_{2}\left|\Phi_{2}\right|^{p-1}+\beta\left|\Phi_{1}\right|^{\frac{p+1}{2}}\left|\Phi_{2}\right|^{\frac{p-3}{2}}\right)=0,(t, x) \in \mathbb{R} \times \Omega \end{array}\right.
(1.2)

其中 \text{i} 表示虚数单位, \Phi_i:\mathbb R^+\times\Omega\to\mathbb C 且对任意的 t>0, i=1,2, \Phi_i(t,\cdot)\in H^1_0(\Omega;\mathbb C). 这类方程组用于描述许多物理现象, 例如耗散量子力学、等离子体物理学和流体力学. 可以见文献[13,33]以及其参考文献了解更多相关的物理背景. 众所周知, 有两个量沿着 1.2)式的轨迹守恒: 质量

\mathcal{Q}(\Phi_i)=\int_\Omega|\Phi_i|^2 {\rm d}x,\ \ i=1,2

和能量

\begin{aligned} \mathcal{E}_{1}\left(\Phi_{1}, \Phi_{2}\right):= & \frac{1}{2} \int_{\Omega}\left(\left|\nabla \Phi_{1}\right|^{2}+\left|\nabla \Phi_{2}\right|^{2}+\left|\Phi_{1}\right|^{2}\left|\nabla \Phi_{1}\right|^{2}+\left|\Phi_{2}\right|^{2}\left|\nabla \Phi_{2}\right|^{2}\right) \mathrm{d} x \\ & -\frac{1}{p+1} \int_{\Omega}\left(\mu_{1}\left|\Phi_{1}\right|^{p+1}+2 \beta\left|\Phi_{1}\right|^{\frac{p+1}{2}}\left|\Phi_{2}\right|^{\frac{p+1}{2}}+\mu_{2}\left|\Phi_{2}\right|^{p+1}\right) \mathrm{d} x \end{aligned}
\begin{aligned} \mathcal{E}_{2}\left(\Phi_{1}, \Phi_{2}\right):= & \frac{1}{2} \int_{\Omega}\left(\left|\nabla \Phi_{1}\right|^{2}+\left|\nabla \Phi_{2}\right|^{2}+\left|\Phi_{1}\right|^{2}\left|\nabla \Phi_{1}\right|^{2}+\left|\Phi_{2}\right|^{2}\left|\nabla \Phi_{2}\right|^{2}\right) \mathrm{d} x \\ & -\frac{1}{p+1} \int_{\Omega}\left(\mu_{1}\left|\Phi_{1}\right|^{p+1}+\mu_{2}\left|\Phi_{2}\right|^{p+1}\right)-\beta \int_{\Omega} \Phi_{1} \Phi_{2} \mathrm{~d} x \end{aligned}

本章寻找(1.2)式的驻波解 (\Phi_1(t,x),\Phi_2(t,x)) = ({\rm e}^{{\rm i}\omega_1 t}u_1(x),{\rm e}^{{\rm i}\omega_2 t}u_2(x))

使得

(u_1,u_2)\in H^1_0(\Omega;\mathbb R^2)

且对某个预定的 \rho_1,\rho_2>0, 有

\begin{equation}\label{eq:mass_constraint} \mathcal{Q}(u_i)=\rho_i,\ \ i=1,2, \end{equation}
(1.3)

另一方面, 方程组(1.1)与拟线性 Schrödinger 方程有关, 这是在等离子体物理学中对超流体膜的研究中产生的. 见文献[9,10,21]查看更多的物理背景和相关的物理模型. 从数学上讲, 这里已有许多研究解的存在性和多重性的结果, 比如见[2,13,33]及其参考文献.

从正规化解的角度, 最简单的情形是 \mathbb R^N 上带纯幂次非线性的单个 Schrödinger 方程的情形. 已有几种经典的方法研究正规化解的存在性和多重性. 比如, 变分法见文献[5,7,37,38]及其参考文献, Jeanjean[18]引入变换 s\ast u(x)={\rm e}^{{\frac N2}s}u({\rm e}^s x) 其中 s\in\mathbb R, Bartsch, Zhong 和 Zou (见文献[6]) 提出的分歧方法, 以及文献[28]中使用的 Lyapunov-Schmidt 约化方法等.

下面介绍以下拟线性 Schrödinger 方程的研究现状

\left\{\begin{array}{l} -\Delta u+\lambda u-\frac{1}{2} \Delta\left(u^{2}\right) u=|u|^{p-1} u \text { in } \mathbb{R}^{N}, \\ \int_{\mathbb{R}^{N}} u^{2} \mathrm{~d} x=\rho>0, \end{array}\right.
(1.4)

其中 \lambda\in\mathbb R 是一个参数, 当 N\geq3p\in(1,\frac{3N+2}{N-2}) 且当 N=1,2p\in(1,+\infty). Colin, Jeanjean, Squassina [14], Jeanjean, Luo [19], Jeanjean, Luo, Wang[20]以及 Li, Zou[20,23]都考虑了单个拟线性 Schrödinger 方程(1.4) 并且应用极小化约束和扰动方法得到了正规化解的存在性. 形式上, (1.4)式的正规化解可以作为

\begin{equation}\label{eq:1.4} I(u):=\frac{1}{2}\int_{\mathbb R^N}(|\nabla u|^2+u^2 |\nabla u|^2) {\rm d}x-\frac{1}{p+1}\int_{\mathbb R^N}|u|^{p+1} {\rm d}x \end{equation}
(1.5)

限制在集合

S(\rho):=\left\{u\in H^1(\mathbb R^N):\int_{\mathbb R^N}|u|^2 {\rm d}x=\rho \right\}

下的临界点得到. 然而, 一旦我们考虑一般的有界区域, 以上方法将会失效. 这是因为文献[20]中引入的变换和文献[5]中的 Nehari-Pohožaev 流形方法在此情形下将不适用, 正如在文献[5,18]及其参考文献中所述, 我们不能得到一个有界的 Palais-Smale 序列. 正如我们所知, 只有少数文章研究有界区域上正规化解的存在性. 大概第一个结果是 Noris, Tavares 和 Verzini 在文献[27]上的结果. 他们证明了有界区域 \Omega\subset\mathbb R^N(N\geqslant1) 上具有纯幂次的方程

\left\{\begin{array}{l} -\Delta u+\lambda u=|u|^{p-1} u, \quad x \in \Omega \\ \int_{\Omega} u^{2} \mathrm{~d} x=\rho>0 \end{array}\right.
(1.6)

这里研究了球面 \Omega=B_1 上正规化解的存在性和轨道稳定性. Pierotti 和 Verzini[32]也考虑了带有特殊边界条件的一般有界区域上的正规化解. 然而, 这两篇论文都是考虑单个的 Schrödinger 方程. 因此, 方程中没有拟线性项和非线性耦合项. 确切地, 方程不包含类似 \Delta(u^2)u|u|^{p-2}u|v|^{p} 的项. 在文献[28]中, Noris, Tavares 和 Verzini 在有界区域上考虑了以下带有齐次 Dirichlet 边界条件的非线性 Schrödinger方程组的正规化解

\left\{\begin{array}{l} -\Delta u_{1}+\omega_{1} u_{1}=\mu_{1}\left|u_{1}\right|^{p-1} u_{1}+\beta\left|u_{2}\right|^{\frac{p+1}{2}}\left|u_{1}\right|^{\frac{p-3}{2}} u_{1}, \\ -\Delta u_{2}+\omega_{2} u_{2}=\mu_{2}\left|u_{2}\right|^{p-1} u_{2}+\beta\left|u_{1}\right|^{\frac{p+1}{2}}\left|u_{2}\right|^{\frac{p-3}{2}} u_{2}, \\ \int_{\Omega}\left|u_{i}\right|^{2} \mathrm{~d} x=\rho_{i}, \quad i=1,2, \quad\left(u_{1}, u_{2}\right) \in H_{0}^{1}\left(\Omega ; \mathbb{R}^{2}\right). \end{array}\right.

据我们所知, 似乎还没有有界区域上拟线性 Schrödinger 方程组正规化解的相关结果.因此,本文的第一个目的是将以上文献 [28] 的结果扩展到拟线性 Schrödinger 方程组.

对于方程组 (1.1), 我们可以在变分框架表述我们的问题, 考虑

\begin{aligned} \mathcal{J}_{1}\left(u_{1}, u_{2}\right):= & \frac{1}{2} \int_{\Omega}\left(\left|\nabla u_{1}\right|^{2}+\left|\nabla u_{2}\right|^{2}+\omega_{1} u_{1}^{2}+\omega_{2} u_{2}^{2}+\left|u_{1}\right|^{2}\left|\nabla u_{1}\right|^{2}+\left|u_{2}\right|^{2}\left|\nabla u_{2}\right|^{2}\right) \mathrm{d} x \\ & -\frac{1}{p+1} \int_{\Omega}\left(\mu_{1}\left|u_{1}\right|^{p+1}+2 \beta\left|u_{1}\right|^{\frac{p+1}{2}}\left|u_{2}\right|^{\frac{p+1}{2}}+\mu_{2}\left|u_{2}\right|^{p+1}\right) \mathrm{d} x \end{aligned}

定义在自然空间 \mathcal{X}:=X_1\times X_2, 其中

X_i:=\left\{u_i\in H_0^1(\Omega): \int_{\Omega}|u_i|^2|\nabla u_i|^2 {\rm d}x<\infty\right\},\ \ i=1,2.

与半线性方程组相比, u_1\Delta(|u_1|^2)u_2\Delta(|u_2|^2) 是不存在的, 寻找 (1.1) 式的解具有一定的困难. 这是因为当 N\geq2 时对应于\int_{\Omega}|u_1|^2|\nabla u_1|^2 {\rm d}x\ \text{和}\ \int_{\Omega}|u_2|^2|\nabla u_2|^2 {\rm d}x 的泛函在空间 \mathcal{X} 上不是可微的. 为了克服这个困难, 应用Liu, Wang 和 Wang [27] 提出的一个方法 (也可见文献 [13]), 即, 作变量替换 v_i=f^{-1}(u_i),\ i=1,2, 其中 f 定义为

\begin{equation}\label{eqf} f'(t)=\frac{1}{\sqrt{1+f^2(t)}},\ \ t\in [0,+\infty), \end{equation}
(1.7)

\begin{equation}\label{eqf1} f(t)=-f(-t),\ \ t\in (-\infty,0]. \end{equation}
(1.8)

应用变量替换, 我们将泛函重新写作 \mathcal{J}(u_1,u_2), 即

\begin{aligned} \mathcal{I}_{1}\left(v_{1}, v_{2}\right):= & \mathcal{J}_{1}\left(u_{1}, u_{2}\right) \\ = & \frac{1}{2} \int_{\Omega}\left(\left|\nabla v_{1}\right|^{2}+\left|\nabla v_{2}\right|^{2}+\omega_{1} f^{2}\left(v_{1}\right)+\omega_{2} f^{2}\left(v_{2}\right)\right) \mathrm{d} x \\ & -\frac{1}{p+1} \int_{\Omega}\left(\mu_{1}\left|f\left(v_{1}\right)\right|^{p+1}+2 \beta\left|f\left(v_{1}\right)\right|^{\frac{p+1}{2}}\left|f\left(v_{2}\right)\right|^{\frac{p+1}{2}}+\mu_{2}\left|f\left(v_{2}\right)\right|^{p+1}\right) \mathrm{d} x \end{aligned}

它在空间 H_0^1(\Omega)\times H_0^1(\Omega) 上是良定义的且是 C^1 的. 此外, 如果 (u_1, u_2)\in H_0^1(\Omega;\mathbb R^2) 且对任意的 \varphi_1,\ \varphi _2\in C_0^\infty(\Omega)

\begin{aligned} \left\langle\mathcal{I}_{1}^{\prime}\left(v_{1}, v_{2}\right),\left(\varphi_{1}, \varphi_{2}\right)\right\rangle:= & \int_{\Omega}\left(\left(1+u_{1}^{2}\right) \nabla u_{1} \nabla \varphi_{1}+\left|\nabla u_{1}\right|^{2} u_{1} \varphi_{1}+\omega_{1} u_{1} \varphi_{1}\right. \\ & \left.+\left(1+u_{2}^{2}\right) \nabla u_{2} \nabla \varphi_{2}+\left|\nabla u_{2}\right|^{2} u_{2} \varphi_{2}+\omega_{2} u_{2} \varphi_{2}\right) \mathrm{d} x \\ & -\int_{\Omega}\left(\mu_{1}\left|u_{1}\right|^{p-1} u_{1} \varphi_{1}+\mu_{2}\left|u_{2}\right|^{p-1} u_{2} \varphi_{2}\right. \\ & \left.+\beta\left|u_{2}\right|^{\frac{p+1}{2}}\left|u_{1}\right|^{\frac{p-3}{2}} u_{1} \varphi_{1}+\beta\left|u_{1}\right|^{\frac{p+1}{2}}\left|u_{2}\right|^{\frac{p-3}{2}} u_{2} \varphi_{2}\right) \mathrm{d} x \\ = & 0 \end{aligned}
(1.9)

(u_1, u_2) 是拟线性方程组(1.1)的一个正规化解, 即, 存在 (\omega_1,\omega_2)\in\mathbb R^2 使得

\left\{\begin{array}{l}-\Delta v_{1}=f^{\prime}\left(v_{1}\right)\left(\mu_{1} f\left(v_{1}\right)\left|f\left(v_{1}\right)\right|^{p-1}+\beta\left|f\left(v_{2}\right)\right|^{\frac{p+1}{2}}\left|f\left(v_{1}\right)\right|^{\frac{p-3}{2}} f\left(v_{1}\right)-\omega_{1} f\left(v_{1}\right)\right), \\-\Delta v_{2}=f^{\prime}\left(v_{2}\right)\left(\mu_{2} f\left(v_{2}\right)\left|f\left(v_{2}\right)\right|^{p-1}+\beta\left|f\left(v_{1}\right)\right|^{\frac{p+1}{2}}\left|f\left(v_{2}\right)\right|^{\frac{p-3}{2}} f\left(v_{2}\right)-\omega_{2} f\left(v_{2}\right)\right), \\\int_{\Omega} f^{2}\left(v_{i}\right) \mathrm{d} x=\rho_{i}, \quad i=1,2, \\\left(v_{1}, v_{2}\right) \in H_{0}^{1}\left(\Omega ; \mathbb{R}^{2}\right).\end{array}\right.
(1.10)

(1.10) 式的解可以看作是泛函

\begin{aligned} \mathcal{E}_{1}\left(v_{1}, v_{2}\right)= & \frac{1}{2} \int_{\Omega}\left(\left|\nabla v_{1}\right|^{2}+\left|\nabla v_{2}\right|^{2}\right) \mathrm{d} x \\ & -\frac{1}{p+1} \int_{\Omega}\left(\mu_{1}\left|f\left(v_{1}\right)\right|^{p+1}+2 \beta\left|f\left(v_{1}\right)\right|^{\frac{p+1}{2}}\left|f\left(v_{2}\right)\right|^{\frac{p+1}{2}}+\mu_{2}\left|f\left(v_{2}\right)\right|^{p+1}\right) \mathrm{d} x \end{aligned}

限制在质量约束

\mathcal{M}=\mathcal{M}_{\rho_1,\rho_2} :=\left\{(v_1,v_2)\in H^1_0(\Omega;\mathbb R^2): \int_\Omega f^2(v_i) {\rm d}x=\rho_i,\ \ i=1,2 \right\}
(1.11)

下的临界点, 这里 \omega_i 作为拉格朗日乘子出现. 类似地, (u_1, u_2) 是拟线性方程组 (1.1) 的一个正规化解, 即, 存在 (\omega_1,\omega_2)\in\mathbb R^2 使得

\left\{\begin{array}{l} -\Delta v_{1}=f^{\prime}\left(v_{1}\right)\left(\mu_{1} f\left(v_{1}\right)\left|f\left(v_{1}\right)\right|^{p-1}+\beta f\left(v_{2}\right)-\omega_{1} f\left(v_{1}\right)\right), \\ -\Delta v_{2}=f^{\prime}\left(v_{2}\right)\left(\mu_{2} f\left(v_{2}\right)\left|f\left(v_{2}\right)\right|^{p-1}+\beta f\left(v_{1}\right)-\omega_{2} f\left(v_{2}\right)\right), \\ \int_{\Omega} f^{2}\left(v_{i}\right) \mathrm{d} x=\rho_{i}, \quad i=1,2 \\ \left(v_{1}, v_{2}\right) \in H_{0}^{1}\left(\Omega ; \mathbb{R}^{2}\right) \end{array}\right.
(1.12)

(1.12) 式的解可以用相应能量泛函

\begin{aligned} \mathcal{E}_{1}\left(v_{1}, v_{2}\right)= & \frac{1}{2} \int_{\Omega}\left(\left|\nabla v_{1}\right|^{2}+\left|\nabla v_{2}\right|^{2}\right) \mathrm{d} x \\ & -\frac{1}{p+1} \int_{\Omega}\left(\mu_{1}\left|f\left(v_{1}\right)\right|^{p+1}+\mu_{2}\left|f\left(v_{2}\right)\right|^{p+1}\right) \mathrm{d} x-\beta \int_{\Omega} f\left(v_{1}\right) f\left(v_{2}\right) \mathrm{d} x \end{aligned}

限制在质量约束

\begin{equation}\label{eq:defM} \mathcal{M}=\mathcal{M}_{\rho_1,\rho_2} :=\left\{(v_1,v_2)\in H^1_0(\Omega;\mathbb R^2): \int_\Omega f^2(v_i) {\rm d}x=\rho_i,\ \ i=1,2 \right\} \end{equation}
(1.13)

下的临界点来确定, 且 \omega_i 作为拉格朗日乘子出现.

我们的主要目标是提供关于 p, \rho_1,\ \rho_2,\ \mu_1,\ \mu_2\beta 的条件, 使得 \mathcal{E}_1|_{\mathcal{M}}\mathcal{E}_2|_{\mathcal{M}} 具有一个全局或局部极小. 我们称这种解是正规化基态解. 此外, 我们考虑这些解相对于演化系统 (1.2) 的稳定性. 到目前为止我们知道, 研究正规化解的一个关键工具是 Gagliardo-Nirenberg 不等式 (见下面的 (1.14) 式), 其可用于根据二次部分来估计能量泛函中的非二次部分. 我们引入一个最佳 Gagliardo-Nirenberg 不等式[2,p458–462]

\begin{equation}\label{1.9} \int_{\mathbb R^N}|u|^{\frac{p+1}{2}} {\rm d}x\leq C_{N,p} \left(\int_{\mathbb R^N}|u| {\rm d}x\right)^{\frac{4N-(N-2)(p+1)}{2(N+2)}} \left(\int_{\mathbb R^N}|\nabla u|^2 {\rm d}x\right)^{\frac{N(p-1)}{2(N+2)}}\ \ \forall \ u\in E^1, \end{equation}
(1.14)

其中 1<p<\frac{3N+2}{N-2}E^q:=\{u\in L^q(\mathbb R^N):\nabla u \in L^2(\mathbb R^N),\ 1\leq q<+\infty\}. 在文献 [35] 中已证

C_{N,p}=\inf_{u\in H^1(\mathbb R^N)\setminus\{0\}}\frac{\|u\|_{L^{\frac{p+1}{2}}(\mathbb R^N)}^{\frac{p+1}{2}}} {\|\nabla u\|_{L^2(\mathbb R^N)}^{\frac{N(p-1)}{N+2}}\|u\|_{L^1(\mathbb R^N)}^{\frac{4N-(N-2)(p+1)}{N+2}}} =\frac{\|Z\|_{L^{\frac{p+1}{2}}(\mathbb R^N)}^{\frac{p+1}{2}}} {\|\nabla Z\|_{L^2(\mathbb R^N)}^{\frac{N(p-1)}{N+2}}\|Z\|_{L^1(\mathbb R^N)}^{\frac{4N-(N-2)(p+1)}{N+2}}},

其中 Z 是方程

\begin{equation}\label{1.10} -\Delta Z+1=Z^{\frac{p-1}{2}} \end{equation}
(1.15)

的唯一正解. 注意到不等式在 H_0^1(\Omega) 中也成立, 对任意有界区域 \Omega, 具有相同的常数 C_{N,p}. 下面, 通过把 (1.14) 式中 u 换为 u^2, 我们得到以下拟线性型 Gagliardo-Nirenberg 不等式

\begin{equation}\label{1.12} \int_{\mathbb R^N }|u|^{p+1} {\rm d}x\leq C_{N,p} \left(\int_{\mathbb R^N}|u|^2 {\rm d}x\right)^{\frac{4N-(N-2)(p+1)}{2(N+2)}} \left(\int_{\mathbb R^N}|\nabla u^2|^2 {\rm d}x\right)^{\frac{N(p-1)}{2(N+2)}}. \end{equation}
(1.16)

特别地, 在 H_0^1(\Omega) 中不等式是严格的, 除非 u 是平凡的. 如上所述, Colin, Jeanjean, Squassina [14], Jeanjean, Luo[19] 以及 Jeanjean, Luo, Wang[20]考虑极小化问题 m(a)=\inf_{u\in S(a)}I(u),I 的定义见 (1.5) 式, 1<p\leq3+\frac{4}{N}. 应用不等式 (1.16), 我们可以发现当 1<p<3+\frac{4}{N}m(a)>-\infty, 当 p>3+\frac{4}{N}m(a)=-\infty, 由于

\frac{N(p-1)}{2(N+2)}<1\ \ \text{当且仅当}\ \ p<3+\frac{4}{N}.

这表明 3+\frac{4}{N} 在拟线性 Schrödinger 方程中作用与 1+\frac{4}{N} 对于非线性 Schrödinger 方程中作用一样. 因此, 在 (1.10) 和 (1.12) 式中的指数 p 可以分为以下四种情形

(H1) 超线性, L^2-次临界: 1<p<3+\frac{4}{N};

(H2) L^2-临界: p=3+\frac{4}{N};

(H3) L^2-超临界, Sobolev-次临界:3+\frac{4}{N}<p<\frac{3N+2}{(N-2)^+};

(H4) Sobolev-临界: p=\frac{3N+2}{N-2}, N\geq3, 见文献[27,命题 5.12]关于该情形更多的性质.

本章的目标有两个: 一方面, 在 (H1)-(H4) 情形, 我们将上述在文献 [28] 中得到的结果推广到拟线性方程组; 另一方面, 我们注意到对于单个拟线性 Schrödinger方程(1.4)这里还没有关于 Sobolev-临界情形的结果, 而其他情形已得到解决. 因此, 我们处理拟线性 Schrödinger 方程组 Sobolev-临界情形 (H4). 现在我们详细描述我们的结果.

\Omega\subset\mathbb R^N(N\geq1) 是 Lipschitz 有界区域, 且令 C_N 表示 Gagliardo-Nirenberg 不等式在 L^2-临界情形下出现的常数 (见 (1.31) 式). 而且, S_N 是出现在 Sobolev 不等式中的最佳常数 (见 (1.32) 式).

首先, 我们给出 (H1) 或 (H2) 成立的情形.

定理1.1 (L^2-次临界和 L^2-临界情形: 存在性和稳定性) 设 \Omega\subset \mathbb R^N 是一个有界光滑区域. 假设以下条件其中之一成立

(i) 1<p<3+\frac{4}{N};

(ii) p=3+\frac{4}{N}\rho_1,\rho_2>0 满足

\begin{equation}\label{N=2condition2} \max\left\{\mu_1\rho_1^{\frac{2}{N}},\ \ \mu_2\rho_2^{\frac{2}{N}},\ \ \mu_1\rho_1^{\frac{2}{N}}+\mu_2\rho_2^{\frac{2}{N}}+\frac{2NC_N}{N+1} \left((\beta^+)^2-\mu_1\mu_2\right)(\rho_1\rho_2)^{\frac{2}{N}}\right\} <\frac{N+1}{2NC_{N}}, \end{equation}
(1.17)

其中 \beta^+=\max\{\beta,0\}.

(a) \inf\limits_{\mathcal{M}}\mathcal{E}_1(v_1,v_2)\in\mathcal{M} 处达到, 它是 (1.10) 式对某个 (\omega_1,\omega_2)\in\mathbb R^2 的一个正解;

(b) 基态集合

\begin{equation}\label{eqG} G=\left\{(v_1,v_2)\in H^1_0(\Omega;\mathbb C^2):\ (|v_1|,|v_2|)\in\mathcal{M},\ \ \mathcal{E}_1(v_1,v_2)=\inf_{\mathcal{M}}\mathcal{E}_1\right\} \end{equation}
(1.18)

是 (条件) 轨道稳定的.

注1.1 轨道稳定性和 (条件) 轨道稳定的定义在第4节. 见文献 [28] 及其参考文献关于更多的轨道稳定性结果. 实际上, 定理 1.1 推广文献 [27,28,32]中的主要结果到拟线性 Schrödinger 方程组.

现在我们考虑 (H3) 或 (H4) 情形成立, 即, 当 3+\frac{4}{N}<p\le\frac{3N+2}{N-2} 时 (当 N=1,2 时没有上界). 令 \alpha\geq \lambda_1(\Omega) (-\Delta\Omega 中的第一 Dirichlet 特征值)

\begin{equation} \mathcal{B}_\alpha:=\left\{(v_1,v_2)\in\mathcal{M}: \int_\Omega(|\nabla v_1|^2+|\nabla v_2|^2) {\rm d}x\leq(\rho_1+\rho_2)\alpha\right\}, \\ \mathcal{U}_\alpha:=\left\{(v_1,v_2)\in\mathcal{M}: \int_\Omega(|\nabla v_1|^2+|\nabla v_2|^2) {\rm d}x=(\rho_1+\rho_2)\alpha\right\}. \end{equation}
(1.19)

注意到 \mathcal{B}_{\alpha} \neq \varnothing , 因为它至少包含一对正规化第一特征函数, 即, (\sqrt{\rho_1} \varphi _1,\sqrt{\rho_2}\varphi _2). 定义

\begin{equation}\label{eq:calpha} c_\alpha:=\inf_{\mathcal{B}_\alpha}\mathcal{E}_1,\ \ \ \ \ \hat{c}_\alpha:=\inf_{\mathcal{U}_\alpha}\mathcal{E}_1. \end{equation}
(1.20)

由于 \mathcal{B}_\alpha\mathcal{M} 中是弱闭的, 在 Sobolev-次临界情形下, c_\alpha 限制在 \mathcal{U}_\alpha \subset\mathcal{B}_\alpha 上对任意的 \alpha\geq\lambda_1(\Omega) 可达. 但是, 为了找到 (1.1) 式的一个解, 只需要寻找 \alpha 使得 c_\alpha<\hat{c}_\alpha. 在 Sobolev-临界情形下, 由于 H^1_0(\Omega) 嵌入到 L^{\frac{4N}{N-2}}(\Omega) 不是紧的, \mathcal{E}_1|_{\mathcal{M}} 不再是弱下半连续的. 为了克服该困难, 我们首先给出方程组 (1.1) 中耦合项的 Brézis-Lieb 引理. 而且, 在 Brézis, Nirenberg [8] 和 Noris, Tavares, Verzini [28] 的启发下, 我们能够通过对质量 \rho_1,\rho_2\alpha 施加一个界恢复相对于 c_\alpha 极小化序列的紧性. 具体地讲, 我们有以下结果

定理 1.2 (Sobolev-临界情形: 存在性) 设 \Omega\subset \mathbb R^N 是一个有界光滑的区域, N\ge3p=\frac{3N+2}{N-2}. 假设 \alpha\geq\lambda_1(\Omega) 满足

\begin{equation}\label{eq:compact_intro} (\rho_1+\rho_2)\left(\alpha-\lambda_1(\Omega)\right)\le\frac{1}{\Lambda^{\frac{N-2}{2} }}, \end{equation}
(1.21)

其中

\begin{equation}\label{eq:def_Lambda} \Lambda:=\frac{2(N-2)S_{N}}{N}\max_{\{x^2+y^2=1\}}\left(\mu_1|x|^{\frac{2N}{N-2}} +\mu_2 |y|^{\frac{2N}{N-2}}+2\beta^+|xy|^{\frac{N}{N-2}}\right). \end{equation}
(1.22)

则相对于 c_\alpha 的任意极小化序列在 \mathcal{B}_\alpha 中是相对紧的. 特别地, c_\alpha 是可达的.

\rho_1=\rho_2=\rho, \mu_1=\mu_2=\mu, u_1=u_2, 我们有以下结论.

定理 1.3 假设 \Omega\subset \mathbb R^N 是一个有界光滑区域且 \mu>0. 如果

0<\rho\le \left(\frac{R(\Omega,N,\frac{3N+2}{N-2})}{\mu}\right)^{\frac{N-2}{2} },

则方程

\left\{\begin{array}{l} -\Delta u+\omega u-\frac{1}{2} u \Delta\left(u^{2}\right)=\mu|u|^{\frac{2(N+2)}{N-2}} u, \\ \int_{\Omega} u^{2} \mathrm{~d} x=\rho, \quad u \in H_{0}^{1}(\Omega), \end{array}\right.

存在一个正解 u, 它是相应能量的一个局部极小元.而且, 相应的局部基态是 (条件) 轨道稳定的.

基于以上定理, 我们引入下面的集合

A:=\left\{\left(\rho_{1}, \rho_{2}\right) \in(0, \infty)^{2}: \begin{array}{l} c_{\alpha}<\hat{c}_{\alpha} \text { 对某个 } \alpha \geq \lambda_{1}(\Omega), \\ \text { 且 } \alpha \text { 满足 (1.21) 式如果 } p=\frac{3 N+2}{N-2} \end{array}\right\} \cup\{(0,0)\}.
(1.23)

注意, A\Omega, N, p, \mu_1, \mu_2\beta 有关. 而且, 如果 (\rho_1,\rho_2)\in A, 则我们可以选择局部极小值 (v_1,v_2)\in\mathcal{M} 是 (1.10) 式对某个 (\omega_1,\omega_2)\in \mathbb R^2 的一个正解. 我们引入 ar

\begin{equation}\label{eq:newexponent} a=a(N,p):=\frac{N(p-1)}{2(N+2)},\ \ \ \ r=r(N,p):=\frac{4N-(N-2)(p+1)}{N+2}. \end{equation}
(1.24)

请注意这两个常数在本文中自然出现, 因为它们出现在 Gagliardo-Nirenberg 不等式. 我们也注意到

\begin{equation}\label{eq:sign_a-1} 0<a<1\ \ \text{当}\ \ 1<p<3+\frac{4}{N};\ \ \ \ a=1\ \ \text{当}\ \ p=3+\frac{4}{N};\ \ \ \ a>1\ \ \text{当}\ \ p>3+\frac{4}{N}. \end{equation}
(1.25)

定理 1.4 (L^2-超临界情形: 存在性) 设 \Omega\subset \mathbb R^N 是一个光滑有界区域且 3+\frac{4}{N}<p\le\frac{3N+2}{N-2}. 如果 A 定义如 (1.23) 式所示, 则 A 相对于 (0,0) 是星型的. 而且, 存在一个正常数 R=R(\Omega,N,p) 使得如果 \rho_1,\rho_2>0 满足

\begin{equation}\label{eq:assnice} \left(\max\{\mu_1\rho_1^{2r},\mu_2\rho_2^{2r}\} + \beta^+\rho_1^r\rho_2^{r} \right)\cdot(\rho_1+\rho_2)^{a-1}\le R(\Omega,N,p), \end{equation}
(1.26)

(\rho_1,\rho_2)\in A. 这里 ar 定义在 (1.24)式且 R 是明确的. (见(3.18) 式).

定理 1.5 (L^2-超临界情形: 稳定性) 设 \Omega\subset \mathbb R^N 是一个有界光滑区域, 3+\frac{4}{N}<p\le\frac{3N+2}{N-2}(\rho_1,\rho_2)\in A.\bar{\alpha}\ge\lambda_1(\Omega) 满足

c_{\bar{\alpha}}<\hat{c} _{\bar{\alpha}},\ \ \text{且}\ \bar{\alpha}\ \text{满足}\ (1.21),\ \ p=\frac{3N+2}{N-2}.

则局部基态集

\begin{equation}\label{eqGa} G_{\bar{\alpha}}:=\left\{(v_1,v_2)\in H^1_0(\Omega;\mathbb C^2):\ (|v_1|,|v_2|)\in\mathcal{B} _{\bar{\alpha}},\ \mathcal{E}_1(v_1,v_2)=c_{\bar{\alpha}}\right\}, \end{equation}
(1.27)

是轨道稳定的.

注意到在我们的结果中所有的假设都涉及 \beta 的正部 (即 \beta^+). 因此, 所有的估计关于 \beta<0 是一致的. 回顾定理 1.1 和 1.4, 当 \beta<0(\rho_1,\rho_2) 满足

\left\{\begin{array}{ll} \rho_{1}, \rho_{2}>0, & 1<p<3+\frac{4}{N}, \\ 0<\mu_{1} \rho_{1}^{\frac{2}{N}}, \mu_{2} \rho_{2}^{\frac{2}{N}}<\frac{N+1}{2 N C_{N}}, & p=3+\frac{4}{N}, \\ \left(\max \left\{\mu_{1} \rho_{1}^{2 r}, \mu_{2} \rho_{2}^{2 r}\right\}\right) \cdot\left(\rho_{1}+\rho_{2}\right)^{a-1} \leq R(\Omega, N, p), & 3+\frac{4}{N}<p \leq \frac{3 N+2}{N-2}. \end{array}\right.
(1.28)

我们陈述以下结果.

定理 1.6\Omega\subset \mathbb R^N 是一个有界光滑区域. 假设 \beta<0(\rho_1,\rho_2) 使得(1.28) 式成立. 假设 (v_{1,\beta},v_{2,\beta}) 是 (1.10) 式的基态, 且 (\omega_{1,\beta},\omega_{2,\beta})\Omega 中满足 v_{1,\beta},v_{2,\beta}>0.\{(v_{1,\beta},v_{2,\beta})\}_{\beta<0} 关于 C^{0,\alpha}(\overline \Omega) 是一致有界的, 当 \beta\to-\infty 时在 C^{0,\alpha}(\overline{\Omega})\cap H^1_0(\Omega) 中有 (v_{1,\beta},v_{2,\beta})\to(w^+,w^-), 其中 w^\pm=\pm\max\{\pm w,0\}w\in C^{0,1}(\overline \Omega) 满足

\left\{\begin{array}{l} -\Delta w-\frac{1}{2} w \Delta\left(w^{2}\right)+\omega_{1} w^{+}-\omega_{2} w^{-}=\mu_{1}\left(w^{+}\right)^{p}-\mu_{2}\left(w^{-}\right)^{p} \\ \int_{\Omega}\left(w^{+}\right)^{2} \mathrm{~d} x=\rho_{1}, \int_{\Omega}\left(w^{-}\right)^{2} \mathrm{~d} x=\rho_{2}, \quad w \in H_{0}^{1}(\Omega) \end{array}\right..

注1.2 不同于文献[14,19,20,20,23], 我们不仅证明正规化解的存在性, 而且还考虑它的轨道稳定性和 \beta\rightarrow -\infty 时的渐近行为. 然而, 当 \beta=\rho_2=0,\ \rho_1=\rho>0 时, 我们的结果对于单个拟线性方程也是新的. 而且, 我们使用的方法是将拟线性问题转化为一个半线性问题, 这与扰动方法不同.

最后, 我们考虑方程组 (1.1) 的正规化解. (1.1) 式的解将作为 C^1 函数的极小值获得

\begin{equation}\label{103} \begin{aligned} \mathcal{E}_2(v_1,v_2)&=\frac{1}{2}\int_{\Omega}(|\nabla v_1|^2+|\nabla v_2|^2) {\rm d}x-\frac{1}{p+1}\int_\Omega(\mu_1|f(v_1)|^{p+1} +\mu_2|f(v_2)|^{p+1}) {\rm d}x\\ &\ \ \ \ +\beta\int_{\Omega}f(v_1)f(v_2) {\rm d}x. \end{aligned} \end{equation}
(1.29)

定理 1.7\Omega\subset \mathbb R^N 是一个有界光滑区域. 假设以下情形之一成立

(i) 1<p<3+\frac{4}{N};

(ii) p=3+\frac{4}{N}\rho_1,\rho_2>0 满足

\begin{equation}\label{N=2condition21} \max\left\{\mu_1\rho_1^{\frac{2}{N}},\ \ \mu_2\rho_2^{\frac{2}{N}},\ \ \mu_1\rho_1^{\frac{2}{N}}+\mu_2\rho_2^{\frac{2}{N}}+ (\beta^+)^2(\rho_1\rho_2)^{\frac{1}{2}}-\frac{2NC_N}{N+1}\mu_1\mu_2(\rho_1\rho_2)^{\frac{2}{N}}\right\} <\frac{N+1}{2NC_{N}}, \end{equation}
(1.30)

其中 \beta^+=\max\{\beta,0\};

(iii) N\ge3p=\frac{3N+2}{N-2};

(iv) 3+\frac{4}{N}<p\le\frac{3N+2}{N-2}.

则对某个 (\omega_1,\omega_2)\in\mathbb R^2 存在 (1.1) 的一个正的正规化解.

注1.3 与文献[28] 相比, 我们还考虑了线性耦合情况下正规化解的存在性. 注意到当线性耦合项 \beta\int_{\Omega}u_1 u_2 {\rm d}x 替换 \beta \int_{\Omega}|u_2|^{\frac{p+1}{2}}|u_1|^{\frac{p-3}{2}}u_1 {\rm d}x, 然而, 我们在这里使用的方法不同于文献[11,12]. 他们的方法离不开平移和变换技巧, 因为我们在有界区域上考虑, 所以这些方法是失效的.

本章的结构如下. 接下来, 我们给出一些注记和记号. 定理 1.1(a) 的证明在第2节. 第 3节主要给出定理1.2-1.4的证明. 第4节关注稳定性结果的证明, 即定理1.1(b)和1.5的证明. 最后, 定理 1.6和定理 1.7的证明在第5节.

记号和预备. 本章中我们记 \lambda_1(\Omega)\Omega 中的第一 Dirichlet 拉普拉斯特征值, 且 \varphi_1 是相应的特征函数, 它在 L^2(\Omega) 中是正规化的且在 \Omega 中是正的. 为了方便起见, CC_i (i = 1, 2,\cdots) 表示正常数(可能是不同的). B_R(x) 是半径为 R 中心在 x 的球. \int_\Omega g(z) 表示积分 \int_\Omega g(z){\rm d}z. \rightarrow\rightharpoonup 分别表示强收敛和弱收敛. 我们使用如下 L^q(\Omega) (1\le q<\infty)H^1_0(\Omega)-范数

|u|_{L^q(\Omega)}^q:=\int_\Omega|u|^q, \ \ \ \ \|u\|_{H^1_0(\Omega)}^2: =\int_\Omega|\nabla u|^2.

为了不引起混淆, 我们将 |\cdot|_{L^q(\Omega)} 简记为 |\cdot|_q. 在特殊的情形 p=3+\frac{4}{N} 时记

\begin{equation}\label{eq:def_CN} C_N:=C_{N,3+\frac{4}{N}}, \end{equation}
(1.31)

然而对于 p=\frac{3N+2}{N-2}N\geq 3,

\begin{equation}\label{eq:Sobolev_constant} S_N:=C_{N,\frac{3N+2}{N-2}}. \end{equation}
(1.32)

2 L^2-次临界和 L^2-临界情形

本节我们处理条件 (H1) 和 (H2) 情形,这意味着

\begin{equation}\label{H1} 1<p\leq3+\frac{4}{N}. \end{equation}
(2.1)

我们收集了函数 f 的一些性质. 证明可见文献[25, p4–5].

引理2.1 函数 f 满足以下性质

(f_{1})f 是唯一定义的且是 C^\infty 的可逆的函数;

(f_{2})|f'(t)|\leq1,\ \ t\in \mathbb R;

(f_{3})|f(t)|\leq|t|,\ \ t\in\mathbb R;

(f_{4}) 存在一个正常数 \theta 使得

|f(t)| \geq\left\{\begin{array}{ll} \theta|t|, & |t| \leq 1, \\ \theta \sqrt{|t|}, & |t|>1 \end{array}\right.

以下引理表明 (1.10) 式的弱解就是 (1.1) 式的弱解.

引理2.2 如果 (v_1,v_2)\in H_0^1(\Omega;\mathbb R^2)\cap L^\infty_{\rm loc}(\Omega;\mathbb R^2) 是泛函 \mathcal{I} 的一个临界点, 则 (u_1,u_2)=(f(v_1),f(v_2)) 是 (1.1) 式的一个弱解.

我们将文献 [36]中研究单个拟线性 Schrödinger 方程的情形推广到方程组情形[13, p217-218]. 通过使用 (1.7) 式和引理 2.1, 我们计算并且得到

|u_i|^2=|f(v_i)|^2\leq|v_i|^2

|\nabla u_i|^2=|f^\prime(v_i)|^2|\nabla v_i|^2\leq|\nabla v_i|^2.

因此, u_1,\ u_2\in H_0^1(\Omega)\cap L_{\rm loc}^\infty(\Omega). 由于 (v_1,v_2)\mathcal{I} 的一个临界点, 对任意的 \psi _1,\ \psi _2\in H_0^1(\Omega),有

\begin{aligned} & \int_{\Omega}\left(\nabla v_{1} \nabla \psi_{1}+\nabla v_{2} \nabla \psi_{2}+\omega_{1} f\left(v_{1}\right) f^{\prime}\left(v_{1}\right) \psi_{1}+\omega_{2} f\left(v_{2}\right) f^{\prime}\left(v_{2}\right) \psi_{2}\right) \\ = & \int_{\Omega} f^{\prime}\left(v_{1}\right) \psi_{1}\left(\mu_{1} f\left(v_{1}\right)\left|f\left(v_{1}\right)\right|^{p-1}+\beta\left|f\left(v_{2}\right)\right|^{\frac{p+1}{2}}\left|f\left(v_{1}\right)\right|^{\frac{p-3}{2}} f\left(v_{1}\right)\right) \\ & +\int_{\Omega} f^{\prime}\left(v_{2}\right) \psi_{2}\left(\mu_{2} f\left(v_{2}\right)\left|f\left(v_{2}\right)\right|^{p-1}+\beta\left|f\left(v_{1}\right)\right|^{\frac{p+1}{2}}\left|f\left(v_{2}\right)\right|^{\frac{p-3}{2}} f\left(v_{2}\right)\right). \end{aligned}
(2.2)

由于 (f^{-1})^\prime(t)=\frac{1}{f^\prime(f^{-1}(t))}, 通过计算

\begin{equation}\label{eq37-2} (f^{-1})^\prime(t)=\sqrt{1+f^2(f^{-1}(t))}=\sqrt{1+t^2}, \end{equation}
(2.3)

这意味着

\begin{equation}\label{eq37-3} \nabla v_i=(f^{-1})^\prime(u_i)\nabla u_i=\sqrt{1+u_i^2}\nabla u_i. \end{equation}
(2.4)

对任意的 \varphi_1,\ \varphi_2\in C_0^\infty(\Omega), 得到

(f^\prime(v_i))^{-1}\varphi_i=\sqrt{1+u_i^2}\varphi_i\in H_0^1(\Omega)

\begin{equation}\label{eq37-4} \nabla((f^\prime(v_i))^{-1}\varphi_i)=\frac{1}{\sqrt{1+u_i^2}}u_i\varphi_i\nabla u_i +\sqrt{1+u_i^2}\nabla\varphi_i. \end{equation}
(2.5)

在 (2.2) 式中取 \psi _i=(f^\prime(v_i))^{-1}\varphi_i 且应用 (2.3), (2.4) 和 (2.5) 式, 我们得到 (1.9) 式, 这意味着 (u_1,u_2)=(f(v_1),f(v_2)) 是 (1.1) 式的一个弱解.

接下来, 我们回顾 Brézis 和 Lieb[7] 给出的一个引理: 对于 1\leq q<\infty, 如果 \{g_n\}_n\subset L^q(\Omega)L^q(\Omega) 中的一个有界序列, 满足 g_n\to g 几乎处处成立, 则

\begin{equation}\label{eq:BL} |g_n|_q^q\to|g|_q^q+|g_n-g|_q^q \ \ \text{当}\ \ n\to\infty. \end{equation}
(2.6)

下面, 我们首先给出方程组 (1.1) 中耦合项的一个 Brézis-Lieb 引理.

引理2.3\omega_{n}\rightharpoonup\omega\nu_{n}\rightharpoonup\nuH_0^1(\Omega) 中成立. 如果 \omega_{n}\rightarrow \omega\nu_{n}\rightarrow\nu\Omega 中几乎处处成立, 则

\begin{equation}\label{eq:mixed_term_conv} \int_{\Omega}|f(\omega_{n})|^{\frac{p+1}{2}} |f(\nu_{n})|^{\frac{p+1}{2}}\to\int_{\Omega}|f(\omega)|^{\frac{p+1}{2}} |f(\nu)|^{\frac{p+1}{2}}+\int_{\Omega}|f(\omega_n)-f(\omega)|^{\frac{p+1}{2}}|f(\nu_{n})-f(\nu)|^{\frac{p+1}{2}}. \end{equation}
(2.7)

通过使用引理 2.1, 得到

\int_{\Omega}|f(\omega_{n})|^{2}\leq\int_{\Omega}|\omega_{n}|^{2}, \ \ \int_{\Omega}|\nabla f(\omega_{n})|^{2} =\int_{\Omega}|f^\prime(\omega_{n})\nabla\omega_{n}|^{2}\leq \int_{\Omega}|\nabla\omega_{n}|^{2},

因此, \{f(\omega_{n})\}H_0^1(\Omega) 中有界. 根据 Sobolev 嵌入的紧性, f(\omega_{n})\in L^{\frac{4N}{N-2}}(\Omega)f(\omega_{n})\rightarrow f(\omega)\Omega 中几乎处处成立. 对于 n=1,\ 2,\cdots, 有

\begin{aligned} & \int_{\Omega}\left|f\left(\omega_{n}\right)\right|^{\frac{p+1}{2}}\left|f\left(\nu_{n}\right)\right|^{\frac{p+1}{2}}-\int_{\Omega}\left|f\left(\omega_{n}\right)-f(\omega)\right|^{\frac{p+1}{2}}\left|f\left(\nu_{n}\right)-f(\nu)\right|^{\frac{p+1}{2}} \\ = & \int_{\Omega}\left(\left|f\left(\omega_{n}\right)\right|^{\frac{p+1}{2}}-\left|f\left(\omega_{n}\right)-f(\omega)\right|^{\frac{p+1}{2}}\right)\left|f\left(\nu_{n}\right)\right|^{\frac{p+1}{2}} \\ & +\int_{\Omega}\left|f\left(\omega_{n}\right)-f(\omega)\right|^{\frac{p+1}{2}}\left(\left|f\left(\nu_{n}\right)\right|^{\frac{p+1}{2}}-\left|f\left(\nu_{n}\right)-f(\nu)\right|^{\frac{p+1}{2}}\right) \end{aligned}

由于 f(\omega_n)\rightharpoonup f(\omega)H_0^1(\Omega) 中成立, 根据文献[26,引理2.5], 得到

\int_{\Omega}\left(|f(\omega_n)|^{\frac{p+1}{2}}- |f(\omega_n)-f(\omega)|^{\frac{p+1}{2}}- |f(\omega)|^{\frac{p+1}{2}}\right)^{2}\rightarrow0\ \ \text{当}\ \ n\rightarrow\infty,

这意味着 在 L^{2}(\Omega)|f(\omega_n)|^{\frac{p+1}{2}}-|f(\omega_n)-f(\omega)|^{\frac{p+1}{2}}\rightarrow |f(\omega)|^{\frac{p+1}{2}}. 应用在 L^{2}(\Omega)|f(\nu_n)|^{\frac{p+1}{2}}\rightharpoonup |f(\nu)|^{\frac{p+1}{2}}, 这意味着

\int_{\Omega}(|f(\omega_n)|^{\frac{p+1}{2}}-|f(\omega_n)-f(\omega)|^{\frac{p+1}{2}})|f(\nu_n)| ^{\frac{p+1}{2}}\rightarrow\int_{\Omega}|f(\omega)|^{\frac{p+1}{2}}|f(\nu)|^{\frac{p+1}{2}}\ \ \text{当}\ \ n\rightarrow\infty.

类似地, 由于在 H_0^1(\Omega)f(\nu_n)\rightharpoonup f(\nu), 根据文献[26,引理2.5]可知在 L^{2}(\Omega)|f(\nu_n)|^{\frac{p+1}{2}}-|f(\nu_n) -f(\nu)|^{\frac{p+1}{2}} \rightarrow|f(\nu)|^{\frac{p+1}{2}}. 根据在 L^{2}(\Omega)|f(\omega_n)-f(\omega)|^{\frac{p+1}{2}}\rightharpoonup 0, 得到

\int_{\Omega}|f(\omega_n)-f(\omega)|^{\frac{p+1}{2}} \left(|f(\nu_n)|^{\frac{p+1}{2}}-|f(\nu_n)-f(\nu)| ^{\frac{p+1}{2}}\right)\rightarrow0\ \ \text{当}\ \ n\rightarrow\infty.

定理 1.1 (a) 的证明 应用 Hölder 不等式和拟线性型 Gagliardo-Nirenberg 不等式(1.16), 有

\begin{aligned} & \int_{\Omega}\left(\mu_{1}\left|f\left(v_{1}\right)\right|^{p+1}+2 \beta\left|f\left(v_{1}\right)\right|^{\frac{p+1}{2}}\left|f\left(v_{2}\right)\right|^{\frac{p+1}{2}}+\mu_{2}\left|f\left(v_{2}\right)\right|^{p+1}\right) \\ \leq & \mu_{1}\left|f\left(v_{1}\right)\right|_{p+1}^{p+1}+\mu_{2}\left|f\left(v_{2}\right)\right|_{p+1}^{p+1}+2 \beta^{+}\left|f\left(v_{1}\right)\right|_{p+1}^{\frac{p+1}{2}}\left|f\left(v_{2}\right)\right|_{p+1}^{\frac{p+1}{2}} \\ < & C_{N, p}\left(\mu_{1} \rho_{1}^{2 r}\left|\nabla f^{2}\left(v_{1}\right)\right|_{2}^{2 a}+\mu_{2} \rho_{2}^{2 r}\left|\nabla f^{2}\left(v_{2}\right)\right|_{2}^{2 a}+2 \beta^{+} \rho_{1}^{r} \rho_{2}^{r}\left|\nabla f^{2}\left(v_{1}\right)\right|_{2}^{a}\left|\nabla f^{2}\left(v_{2}\right)\right|_{2}^{a}\right), \end{aligned}
(2.8)

其中指数 ar 定义在 (1.24)式. 由 (1.7) 式可知

\int_{\Omega}|\nabla f^2(v)|^2 =\int_{\Omega}\Bigl|\frac{2f(v)\nabla v}{\sqrt{1+f^2(v)}} \Bigr|^2 \leq4\int_{\Omega}|\nabla v|^2.

因此, 对 (v_1,v_2)\in\mathcal{M}

\begin{aligned} \mathcal{E}_{1}\left(v_{1}, v_{2}\right)> & \frac{1}{2}\left(\left|\nabla v_{1}\right|_{2}^{2}+\left|\nabla v_{2}\right|_{2}^{2}\right) \\ & -\frac{4 C_{N, p}}{p+1}\left(\mu_{1} \rho_{1}^{2 r}\left|\nabla v_{1}\right|_{2}^{2 a}+\mu_{2} \rho_{2}^{2 r}\left|\nabla v_{2}\right|_{2}^{2 a}+2 \beta^{+} \rho_{1}^{r} \rho_{2}^{r}\left|\nabla v_{1}\right|_{2}^{a}\left|\nabla v_{2}\right|_{2}^{a}\right) \end{aligned}
(2.9)

如果 (H1) 成立, 则 0<\frac{N(p-1)}{N+2}<2, 因此, 在 (2.9) 式中, a<1, 推导出 \mathcal{E}_1 相对于 \mathcal{M} 对任意的 \rho_1,\ \rho_2>0 是强制的. 我们称 \mathcal{M} 是一类 C^1 流形. 事实上, 根据

G(v_i):=\int_\Omega|f(v_i)|^2=\rho_i>0,\ i=1,2,

由引理 2.1 可知 G 是一类 C^1 的且对 (v_1,v_2)\in\mathcal{M}, 得到

G^\prime(v_i)v_i=2\int_\Omega |f(v_i)|f^\prime(v_i)v_i =2\int_\Omega\frac{f(v_i)v_i}{\sqrt{1+v_i^2}}\neq0,\ i=1,2.

因此, 应用隐函数定理, \mathcal{M} 是一类 C^1 流形. 最后, 根据引理 2.3, Sobolev 嵌入的紧性和范数的弱下半连续性, 直接应用变分技巧 [4,注 1.5.7], \inf\limits_{\mathcal{M}}\mathcal{E}_1(v_1,v_2) 处达到, 其中 (v_1,v_2)\in\mathcal{M}, 这是因为紧嵌入 H^1_0(\Omega)\hookrightarrow L^2(\Omega). 由拉格朗日乘子原理, 对某个 (\omega_1,\omega_2)\in\mathbb R^2, (v_1,v_2) 满足 (1.10) 式. 同时, 引理2.2意味着对某个 (\omega_1,\omega_2)\in\mathbb R^2, (f(v_1),f(v_2)) 满足 (1.1) 式. 取 |v_i|, 我们假设 v_i\geq0. 根据极大值原理, 我们得到 v_i>0, 由于 \Omega 是 Lipschitz 的, 任意 v_i 关于边界是连续的.

在情形 (H2) 下, 根据 (2.9) 式, 此时 a=1, r=\frac{1}{N}, 有

\begin{aligned} \mathcal{E}_{1}\left(v_{1}, v_{2}\right)> & \frac{1}{2}\left(\left|\nabla v_{1}\right|_{2}^{2}+\left|\nabla v_{2}\right|_{2}^{2}\right) \\ & -\frac{4 C_{N, p}}{p+1}\left(\mu_{1} \rho_{1}^{2 r}\left|\nabla v_{1}\right|_{2}^{2 a}+\mu_{2} \rho_{2}^{2 r}\left|\nabla v_{2}\right|_{2}^{2 a}+2 \beta^{+} \rho_{1}^{r} \rho_{2}^{r}\left|\nabla v_{1}\right|_{2}^{a}\left|\nabla v_{2}\right|_{2}^{a}\right) \\ = & \frac{1}{2}\left|\nabla v_{1}\right|_{2}^{2}\left(1-\frac{2 N C_{N} \mu_{1} \rho_{1}^{\frac{2}{N}}}{N+1}\right)+\frac{1}{2}\left|\nabla v_{2}\right|_{2}^{2}\left(1-\frac{2 N C_{N} \mu_{2} \rho_{2}^{\frac{2}{N}}}{N+1}\right) \\ & -\frac{2 \beta^{+} N C_{N}\left(\rho_{1} \rho_{2}\right)^{\frac{1}{N}}}{N+1}\left|\nabla v_{1}\right|_{2}\left|\nabla v_{2}\right|_{2} \\ = & \frac{1}{2}\left[\left|\nabla v_{1}\right|_{2}\left|\nabla v_{2}\right|_{2}\right] \cdot A \cdot\left[\left|\nabla v_{1}\right|_{2}\left|\nabla v_{2}\right|_{2}\right]^{T}, \end{aligned}
(2.10)

其中 \begin{bmatrix}|\nabla v_1|_2 ~&|\nabla v_2|_2\end{bmatrix}^T\begin{bmatrix}|\nabla v_1|_2~ &|\nabla v_2|_2\end{bmatrix} 的转置且

A=\begin{bmatrix} 1-\frac{2NC_N\mu_1\rho_1^{\frac{2}{N}}}{N+1} & -\frac{2\beta^+ N C_N(\rho_1\rho_2)^{\frac{1}{N}}}{N+1}\\[3mm] -\frac{2\beta^+N C_N(\rho_1\rho_2)^{\frac{1}{N}}}{N+1} & 1-\frac{2NC_N\mu_2\rho_2^{\frac{2}{N}}}{N+1}\end{bmatrix}.

如果 A 是正定的, 我们得到的 \mathcal{E}_1(v_1,v_2)>0, 从而结果成立. 这里, A 正定当且仅当下面的不等式同时成立

\begin{array}{l} 1-\frac{2 N C_{N} \mu_{1} \rho_{1}^{\frac{2}{N}}}{N+1}>0, \quad 1-\frac{2 N C_{N} \mu_{2} \rho_{2}^{\frac{2}{N}}}{N+1}>0 \\ \left(1-\frac{2 N C_{N} \mu_{1} \rho_{1}^{\frac{2}{N}}}{N+1}\right)\left(1-\frac{2 N C_{N} \mu_{2} \rho_{2}^{\frac{2}{N}}}{N+1}\right)-\left(\frac{2 N C_{N}}{N+1}\right)^{2}\left(\beta^{+}\right)^{2}\left(\rho_{1} \rho_{2}\right)^{\frac{2}{N}}>0 \end{array}

即, (1.17) 式成立. 我们也推导出 \mathcal{E}_1 关于 \mathcal{M} 对任意的 \rho_1,\ \rho_2>0 是强制的.

注2.1 受文献 [28] 启发, 我们在 (\mu_1\rho_1^{\frac{2}{N}},\mu_2\rho_2^{\frac{2}{N}})-平面引入条件 (1.17). 不同于文献 [28,图 1], 这里的双曲线包含了点 \left(0,\frac{N+1}{2NC_{N}}\right)\left(\frac{N+1}{2NC_{N}},0\right).\bar x=\mu_1\rho_1^{\frac{2}{N}}, \bar y=\mu_2\rho_2^{\frac{2}{N}} 使得 (1.17) 式对应于 (\bar x,\bar y)\in G, 其中

G:=\left\{(x,y)\in\mathbb R^2:\ 0<x,y<\frac{N+1}{2NC_{N}},\ \ x+y+\frac{2NC_N}{N+1}\frac{(\beta^+)^2- \mu_1\mu_2}{\mu_1\mu_2}xy<\frac{N+1}{2NC_{N}}\right\}.

根据 \beta, 我们可以分为以下几种情形.

(i) 如果 \beta\leq 0, (1.17) 式变成矩形中的线 (\bar x,\bar y)

Q:=\left\{(x,y)\in\mathbb R^2:\ 0< x,y<\frac{N+1}{2NC_{N}} \right\}.

(ii) 如果 \beta=\sqrt{\mu_1\mu_2}, 有半个矩形

Q_1:=Q\cap\left\{(x,y)\in\mathbb R^2:\ x+y<\frac{N+1}{2NC_{N}}\right\}.

(iii) 如果 \beta>0, \beta\neq \sqrt{\mu_1\mu_2}, 曲线

\left\{(x,y)\in\mathbb R^2:\ x+y+\frac{2NC_N}{N+1} \frac{(\beta^+)^2-\mu_1\mu_2}{\mu_1\mu_2}xy =\frac{N+1}{2NC_{N}}\right\}

是一个双曲线它包含点 \left(0,\frac{N+1}{2NC_{N}}\right), \left(\frac{N+1}{2NC_{N}},0\right). 这个双曲线是

y=\left(\frac{N+1}{2NC_{N}}-x\right) \left(1+\frac{2NC_N}{N+1}\frac{(\beta^+)^2-\mu_1\mu_2}{\mu_1\mu_2}x\right)^{-1},

y=\frac{N+1}{2NC_{N}} \left(-\frac{\mu_1\mu_2}{(\beta^+)^2-\mu_1\mu_2}+\frac{N+1}{2NC_{N}} \frac{(\beta^+)^2}{(\beta^+)^2-\mu_1\mu_2} \left(\frac{N+1}{2NC_{N}}+\frac{(\beta^+)^2 -\mu_1\mu_2}{\mu_1\mu_2}x\right)^{-1}\right),

其具有垂直渐近线 x=\frac{N+1}{2NC_{N}} \frac{\mu_1\mu_2}{\mu_1\mu_2-(\beta^+)^2}. 但是, 集合 G 总是包含 Ox^+\cap\overline{Q}Oy^+\cap\overline{Q}.0<\beta<\sqrt{\mu_1\mu_2}, 它包含 Q_1, 并且当 \beta>\sqrt{\mu_1\mu_2} 包含在 Q_1 中.

注2.2 如果 \beta\leq 0, 条件 (1.17) 变成 \mu_1\rho_1^{\frac{2}{N}}, \mu_2\rho_2^{\frac{2}{N}}<\frac{N+1}{2NC_{N}}. 回到文献 [20,引理 2.1], 我们看到 Pohožaev 恒等式 (u 换成 f(v))

\frac{N+1}{2NC_{N}}=\frac{N+2}{2}|f(Z)|_2^\frac{4}{N},

Z 是拟线性 Schrödinger 方程唯一的正解[1,3,16,34]. 因此, (1.17) 式等价于

\rho_1<\left(\frac{N+2}{2}|f(Z)|_2^\frac{4}{N}\right)^{\frac{N}{2} } \mu_1^{-\frac{N}{2}} ~\text{和}~ \rho_2<\left(\frac{N+2}{2}|f(Z)|_2^\frac{4}{N}\right)^{\frac{N}{2} } \mu_2^{-\frac{N}{2}}.

这与文献 [1,3,16,34] 中的结果是一致的, 这对应于 (1.1) 式中 \beta=0 的情形.

3 L^2-超临界, Sobolev-次临界和 Sobolev-临界情形

现在开始假设 p 满足 (H3) 或 (H4), 即, p>3+\frac{4}{N}, 且 p\leq\frac{3N+2}{N-2}, N\geq 3.L^2-超临界情形 p>3+\frac{4}{N}, 由于 \mathcal{E}_1 限制在 \mathcal{M} 上对任意的 (\rho_1,\rho_2) 是强制的, 正如以下引理.

引理3.1p>3+\frac{4}{N}, 存在 (V_{1,k},V_{2,k})\in\mathcal{M}, 使得

\|(V_{1,k},V_{2,k})\|_{H^1_0(\Omega)}\to+\infty\ \ \text{且}\ \ \mathcal{E}_1(V_{1,k},V_{2,k})\to-\infty\ \ \text{当}\ \ k\to\infty.

\phi\in C^\infty_c(\Omega)\phi>0B_1 中, \int_{B_1}\phi^2=1, x_1,x_2\in\Omega 使得 x_1\neq x_2.k\in\mathbb Ni=1,2, 我们定义

V_{i,k}(x)=f^{-1}\left(\rho_i^{\frac{1}{2}}k^{\frac{N}{2}} \phi(k(x-x_i))\right),\ \ x\in\Omega.

k 充分大, 我们得到 \textrm{supp}(V_{i,k}) \subset B_{\frac{1}{k}}(x_i)\subset\Omega, i=1,2, 且 \textrm{supp}(V_{1,k}) \cap \textrm{supp}(V_{2,k})=\emptyset. 通过直接计算我们得到

\int_\Omega f^{2}(V_{i,k})=\int_\Omega \left(\rho_i^{\frac{1}{2}}k^{\frac{N}{2}} \phi(k(x-x_i))\right)^{2} =\rho_i\int_{B_1}\phi^2=\rho_i,

因此, 对充分大的 k(V_{1,k},V_{2,k})\in\mathcal{M}. 此外, 根据 (2.3) 式可知

\begin{aligned} \left\|\left(V_{1, k}, V_{2, k}\right)\right\|_{H_{0}^{1}(\Omega)} & =\left(\int_{\Omega}\left|\nabla V_{1, k}\right|^{2}+\int_{\Omega}\left|\nabla V_{2, k}\right|^{2}\right)^{\frac{1}{2}} \\ & =\left(\left(\rho_{1}+\rho_{2}\right) k^{2} \int_{\Omega}|\nabla \phi|^{2}+\left(\rho_{1}^{2}+\rho_{2}^{2}\right) k^{N+2} \int_{\Omega}|\phi|^{2}|\nabla \phi|^{2}\right)^{\frac{1}{2}} \\ & \leq k\left(\sqrt{\rho_{1}+\rho_{2}}|\nabla \phi|_{L^{2}\left(B_{1}\right)}+\frac{k^{\frac{N}{2}} \sqrt{\rho_{1}^{2}+\rho_{2}^{2}}}{2}\left|\nabla \phi^{2}\right|_{L^{2}\left(B_{1}\right)}\right) \\ & \rightarrow+\infty \text { 当 } k \rightarrow \infty \end{aligned}

因此, 当 k\to+\infty 时,

\begin{aligned} \mathcal{E}_{1}\left(V_{1, k}, V_{2, k}\right) \leq & \frac{1}{2} k^{2}\left(\left(\rho_{1}+\rho_{2}\right)|\nabla \phi|_{L^{2}\left(B_{1}\right)}+\frac{k^{\frac{N}{2}} \sqrt{\rho_{1}^{2}+\rho_{2}^{2}}}{2}\left|\nabla \phi^{2}\right|_{L^{2}\left(B_{1}\right)}\right)^{2} \\ & -\frac{\mu_{1} \rho_{1}^{\frac{p+1}{2}}+\mu_{2} \rho_{2}^{p+1}}{p+1} k^{\frac{N(p-1)}{2}}|\phi|_{L^{p+1}\left(B_{1}\right)}^{p+1} \rightarrow-\infty. \end{aligned}

这是因为 p>3+\frac{4}{N} .

在 (H3) 或 (H4) 下, 为了证明对某个 \rho_1,\rho_2 存在 (1.10) 式的一个解, 我们应用不同于第2节的方法. 对 \alpha\geq \lambda_1(\Omega), \mathcal{B} _\alpha\mathcal{U} _\alpha 定义在 (1.19)式, 然而 c_\alpha\hat c_\alpha 定义在 (1.20) 式. 注意到 \mathcal{B}_\alpha\neq \emptyset, 由于它至少包含 (\sqrt{\rho_1}\varphi _1,\sqrt{\rho_2}\varphi _2). 而且

c_{\lambda_1(\Omega)}=\hat c_{\lambda_1(\Omega)}= \mathcal{E}_1(\sqrt{\rho_1}\varphi _1,\sqrt{\rho_2}\varphi _2).

回顾 (2.9) 式和 x=|\nabla v_1|_2 以及 y=|\nabla v_2|_2, 考虑函数 \Phi:\mathbb R^2_+\to\mathbb R 定义为

\begin{equation}\label{eq:vphi_def} \Phi(x,y)=\frac{1}{2}(x^2+y^2)-\frac{4C_{N,p}}{p+1}\left(\mu_1\rho_1^{2r}x^{2a}+\mu_2\rho_2^{2r}y^{2a}+2\beta^+ \rho_1^r\rho_2^{r}x^{a}y^{a}\right), \end{equation}
(3.1)

其中 a>1. 事实上, 根据 (2.9) 式可知

\begin{equation}\label{eq:Phi_ineq} \mathcal{E}_1(v_1,v_2)>\Phi \left(|\nabla v_1|_2,|\nabla v_2|_2\right) \ \ \text{对任意的}\ \ (v_1,v_2)\in\mathcal{M}. \end{equation}
(3.2)

特别地, 这里可以从下方估计 \hat c_\alpha. 我们定义以下 \mathbb R^2 的子集合

U_\alpha=\left\{(x,y)\in\mathbb R^2_+:\ x^2+y^2=(\rho_1+\rho_2)\alpha\right\}, \ \ V_\alpha=U_\alpha\cap \left\{x\geq \sqrt{\rho_1\lambda_1(\Omega)},\ y\geq \sqrt{\rho_2\lambda_1(\Omega)}\right\}.

集合 U_\alpha 是由 \mathcal{U}_\alpha 通过定义 x=|\nabla v_1|_2y=|\nabla v_2|_2 得到的. 集合 V_\alpha 是由 (v_1,v_2)\in\mathcal{M}, |\nabla v_1|_2^2\geq\rho_1\lambda_1(\Omega) 以及 |\nabla v_2|_2^2\geq\rho_2\lambda_1(\Omega) 得到. 根据 (3.2)式, 我们推导出

\begin{equation}\label{eq:iniziostima} \hat{c}_\alpha\ge\min_{(x,y)\in V_\alpha}\Phi(x,y)= \frac{1}{2 }(\rho_1+\rho_2)\alpha-\frac{4C_{N,p}}{p+1}\max_{(x,y)\in V_\alpha} \left(\mu_1\rho_1^{2r}x^{2a}+\mu_2\rho_2^{2r}y^{2a}+2\beta^+ \rho_1^r\rho_2^{r}x^{a}y^{a}\right). \end{equation}
(3.3)

现在, 由于 V_\alpha 的定义, 最后一个最大值不能由 \alpha 精确给出 (除了特殊情形). 因此, 我们希望得到一个粗略的估计, 其中 V_\alphaU_\alpha 代替.

引理3.2 对任意的 \alpha>\lambda_1(\Omega),

\hat c_\alpha>\frac12\left((\rho_1+\rho_2)\alpha-\Lambda(\rho_1,\rho_2)(\rho_1+\rho_2)^a\alpha^a \right),

其中

\begin{equation}\label{eq:defdd} \Lambda =\Lambda(\rho_1,\rho_2):=\frac{8C_{N,p}}{p+1}\max_{t\in[\frac{\pi}{2}]}\left(\mu_1\rho_1^{2r}\cos^{2a}t+ \mu_2\rho_2^{2r}\sin^{2a}t+2\beta^+\rho_1^r\rho_2^{r}\cos^{a}t \sin^{a}t\right). \end{equation}
(3.4)

显然, 因为

V_\alpha\subset U_\alpha=\left\{(\cos t,\sin t)\sqrt{(\rho_1+\rho_2) \alpha}: t\in\left[\frac{\pi}{2}\right]\right\},

由 (3.3)式即证.

注3.1N\ge3p=\frac{3N+2}{N-2}, 我们知道 a=\frac{N}{N-2}, r=0 以及 \Lambda 不依赖于 \rho_1,\rho_2, 实际上它的定义与 (1.22)式给出的一致. 因此, 根据 S_N 在 (1.32) 式中的定义, 对任意的 (\tau_1,\tau_2)\in H^1_0(\Omega;\mathbb R^2), 得到

\begin{aligned} & \frac{2(N-2) S_{N}}{N}\left(\mu_{1}\left|\nabla \tau_{1}\right|_{2}^{\frac{2 N}{N-2}}+2 \beta^{+}\left|\nabla \tau_{1}\right|_{2}^{\frac{N}{N-2}}\left|\nabla \tau_{2}\right|_{2}^{\frac{N}{N-2}}+\mu_{2}\left|\nabla \tau_{2}\right|_{2}^{\frac{2 N}{N-2}}\right) \\ \leq & \Lambda\left(\left|\nabla \tau_{1}\right|_{2}^{2}+\left|\nabla \tau_{2}\right|_{2}^{2}\right)^{\frac{N}{N-2}}. \end{aligned}
(3.5)

为了看到这点, 对任意的 (\tau_1,\tau_2), 可以找到 t\in\left[\frac{\pi}{2}\right] 使得

|\nabla\tau_1|_2=\left(|\nabla\tau_1|_2^{2}+|\nabla\tau_{2}|_2^2\right)^{\frac{1}{2}} \cos t, \ \ |\nabla\tau_2|_2=\left(|\nabla\tau_1|_2^{2}+|\nabla\tau_{2}|_2^2\right) ^{\frac{1}{2}}\sin t,

并且代入 (3.4)式中即证.

正如之前提到的, 我们将寻找 c_\alpha 对某个 \alpha, \mathcal{E}_1\mathcal{B}_\alpha 上的极小. 第一步是证明 c_\alpha 可达. 对任意 \alpha\ge\lambda_1(\Omega), 在 Sobolev-次临界情形下很容易得到. 事实上, 此时, \mathcal{B}_\alpha 是弱紧的且 \mathcal{E}_1 是弱下半连续的, 使用变分法 [4,注 2.3.3], c_\alpha(v_1,v_2) 处达到.

接下来, 回顾(3.4)式中的 \Lambda, 在 Sobolev-临界情形下, 与 \rho_1,\rho_2 无关.

引理3.3\alpha>\lambda_1(\Omega), \rho_1,\rho_2>0 满足

(\rho_1+\rho_2)(\alpha-\lambda_1(\Omega))<\frac{1}{\Lambda^{\frac{N-2}{2}}},

且令 (v_{1,n},v_{2,n})_n 满足当 n\to\infty

\left\{\begin{array}{l} \left|v_{i, n}\right|_{2}^{2}=\rho_{i}+o(1), \quad i=1,2, \\ \left|\nabla v_{1, n}\right|_{2}^{2}+\left|\nabla v_{2, n}\right|_{2}^{2} \leq \alpha\left(\rho_{1}+\rho_{2}\right)+o(1), \\ c_{\alpha} \leq \mathcal{E}_{1}\left(v_{1, n}, v_{2, n}\right) \leq c_{\alpha}+o(1). \end{array}\right.
(3.6)

因此

(v_{1,n},v_{2,n})\to(\bar v_{1},\bar v_{2})\ \ \text{ 在}\ \ H^1_0(\Omega).

特别地, c_\alpha 可达.

根据假设条件, 存在 (\bar v_1,\bar v_2)\in H^1_0(\Omega;\mathbb R^2) 使得

\left\{\begin{array}{ll} \left|\bar{v}_{i}\right|_{2}^{2}=\rho_{i}, & i=1,2, \\ v_{i, n} \rightharpoonup \bar{v}_{i} \text { 在 } H_{0}^{1}(\Omega), & i=1,2, \\ \left|\nabla \bar{v}_{i}\right|_{2}^{2} \leq \liminf _{n \rightarrow \infty}\left|\nabla v_{i, n}\right|_{2}^{2}, & i=1,2. \end{array}\right.

注意到 (\bar{v}_1,\bar{v}_2)c_\alpha 的极小化问题, 从而

\begin{equation}\label{eq:bar_u_i_admissible} \mathcal{E}_1(\bar{v}_1,\bar{v}_2)\geq c_\alpha. \end{equation}
(3.7)

\tau_{i,n}=v_{i,n}-\bar v_i, 注意到对 i=1,2,

\begin{equation}\label{eq:v_i_n} \tau_{i,n}\rightharpoonup 0 \ \ \text{在}\ \ H^1_0(\Omega)\ \ \text{且}\ \ L^{\frac{4N}{N-2}}(\Omega), \ \ \tau_{i,n}\to0 \ \ \text{在}\ \ L^2(\Omega). \end{equation}
(3.8)

(v_{1,n},v_{2,n}) 的强收敛性

\begin{equation}\label{eq:BN_easy_case} \text{ 存在一个子列 }\ \ (\tau_{1,n_k},\tau_{2,n_k})\ \ \text{使得}\ \ |\nabla \tau_{1,n_k}|_2^2+|\nabla \tau_{2,n_k}|_2^2\to 0. \end{equation}
(3.9)

此时, 由 Sobolev 嵌入的连续性, 我们得到 c_\alpha=\mathcal{E}_1(\bar{v}_1,\bar{v}_2). 由于 c_\alpha 的一个极小化序列存在且满足 (3.6) 式, 我们得到 c_\alpha 可达.

为了完成证明, 通过矛盾, 假设 (3.9)式不成立, 因此

\begin{equation}\label{eq:BN_difficult_case} |\nabla \tau_{1,n}|_2^2+|\nabla \tau_{2,n}|_2^2\geq K>0. \end{equation}
(3.10)

\begin{aligned} & \mathcal{E}_{1}\left(v_{1, n}, v_{2, n}\right) \\ = & \frac{1}{2}\left(\left|\nabla\left(\bar{v}_{1}+\tau_{1, n}\right)\right|_{2}^{2}+\mid \nabla\left(\bar{v}_{2}+\left.\tau_{2, n}\right|_{2} ^{2}\right)\right) \\ & -\frac{N-2}{4 N}\left(\mu_{1}\left|f\left(\bar{v}_{1}+\tau_{1, n}\right)\right|_{\frac{4 N}{N-2}}^{\frac{4 N}{N-2}}+2 \beta\left|f\left(\bar{v}_{1}+\tau_{1, n}\right) f\left(\bar{v}_{2}+\tau_{2, n}\right)\right|_{\frac{2 N}{N-2}}^{\frac{2 N}{N-2}}+\mu_{2}\left|f\left(\bar{v}_{2}+\tau_{2, n}\right)\right|_{\frac{4 N}{N-2}}^{\frac{4 N}{N-2}}\right). \end{aligned}
(3.11)

根据 v_{i,n}\rightharpoonup \bar v_iH^1_0(\Omega) 中成立,对 i=1,2,

\begin{equation}\label{eq:weak_conv} |\nabla (\bar{v}_i+\tau_{i,n})|_2^2= |\nabla \bar{v}_i|_2^2+|\nabla\tau _{i,n}|_2^2+o(1)\ \ \text{当}\ \ n\to\infty. \end{equation}
(312)

根据引理 2.1 和 Sobolev 嵌入的紧性, 我们应用 (2.6) 式且 g_n=f(v_{i,n})=f(\bar v_i+\tau_{i,n}), q=\frac{4N}{N-2}, 为了得到

\begin{equation}\label{eq:L4conv} |f(\bar v_i+\tau_{i,n})|_{\frac{4N}{N-2}}^{\frac{4N}{N-2}}= |f(\bar v_i)|_{\frac{4N}{N-2}}^{\frac{4N}{N-2}}+|f(\tau_{i,n})| _{\frac{4N}{N-2}}^{\frac{4N}{N-2}} +o(1) \ \ \text{当}\ \ n\to\infty. \end{equation}
(3.13)

通过替换 (3.12), (3.13) 和 (2.7) 式且 p=\frac{3N+2}{N-2} 代入 (3.11) 式, 我们得到

\mathcal{E}_1(v_{1,n},v_{2,n})=\mathcal{E}_1(\bar v_1,\bar v_2)+\mathcal{E}_1(\tau_{1,n}, \tau _{2,n})+o(1) \ \ \text{当}\ \ n\to\infty.

最后一个式子以及 (3.6) 和 (3.7) 式意味着

\mathcal{E}_1(\tau _{1,n},\tau _{2,n})\leq o(1) \ \ \text{当}\ \ n\to\infty,

从而, 根据 (2.8), (2.9) (r=0, a=\frac{N}{N-2}C_{N,\frac{3N+2}{N-2}}=S_N), (3.5) 式可知

\begin{aligned} \left|\nabla \tau_{1, n}\right|_{2}^{2}+\left|\nabla \tau_{2, n}\right|_{2}^{2} & \leq \frac{2(N-2)}{N} S_{N}\left(\mu_{1}\left|\tau_{1, n}\right|_{2}^{\frac{2 N}{N-2}}+2 \beta\left|\tau_{1, n} \tau_{2, n}\right|_{2}^{\frac{N}{N-2}}+\mu_{2}\left|\tau_{2, n}\right|_{2}^{\frac{2 N}{N-2}}\right)+o(1) \\ & \leq \frac{2(N-2)}{N} S_{N}\left(\mu_{1}\left|\tau_{1, n}\right|_{2}^{\frac{2 N}{N-2}}+2 \beta^{+}\left|\tau_{1, n}\right|_{2}^{\frac{N}{N-2}}\left|\tau_{2, n}\right|_{2}^{\frac{N}{N-2}}+\mu_{2}\left|\tau_{2, n}\right|_{2}^{\frac{2 N}{N-2}}\right)+o(1) \\ & \leq \Lambda\left(\left|\nabla \tau_{1, n}\right|_{2}^{2}+\left|\nabla \tau_{2, n}\right|_{2}^{2}\right)^{\frac{N}{N-2}}+o(1) \end{aligned}

接下来, 我们应用 (3.10) 式重写最后一个不等式得

\left(|\nabla \tau_{1,n}|_2^2+|\nabla \tau_{2,n}|_2^2\right)^{\frac{2}{N-2}} \ge\frac{1}{\Lambda }+o(1).

结合前边的不等式以及 (3.12) 式, 我们得到当 n\to+\infty 时,

\begin{aligned} \left(\frac{1}{\Lambda}+o(1)\right)^{\frac{N-2}{2}} & \leq\left|\nabla \tau_{1, n}\right|_{2}^{2}+\left|\nabla \tau_{2, n}\right|_{2}^{2} \\ & =\left|\nabla v_{1, n}\right|_{2}^{2}+\left|\nabla v_{2, n}\right|_{2}^{2}-\left(\left|\nabla \bar{v}_{1}\right|_{2}^{2}+\left|\nabla \bar{v}_{2}\right|_{2}^{2}\right)+o(1) \\ & \leq\left(\rho_{1}+\rho_{2}\right) \alpha-\lambda_{1}(\Omega)\left(\rho_{1}+\rho_{2}\right)+o(1) \end{aligned}

这与假设矛盾.

定理 1.2 的证明 该定理是引理 3.3 当 (3.6) 式且 |v_{i,n}|_2^2=\rho_i, i=1,2, 以及 |\nabla v_{1,n}|_2^2+|\nabla v_{2,n}|_2^2\leq\alpha(\rho_1+\rho_2) 时的一个特殊情形.

引理3.4 假设 \rho_1,\rho_2>0 满足对某个 \alpha_1,\alpha_2, \lambda_1(\Omega)\le\alpha_1<\alpha_2\hat c_{\alpha_1}<\hat c_{\alpha_2}; 此外, Sobolev-临界情形 N\ge3, p=\frac{3N+2}{N-2} , 也假设

\alpha_2<\lambda_1(\Omega) +\frac{\Lambda ^{-\frac{N-2}{2}}}{\rho_1+\rho_2}.

因此 c_{\alpha_2}<\hat c_{\alpha_2}, 且 c_{\alpha_2} 达到 (1.10) 式的一个正解.

已知 c_{\alpha_2}(\bar {v}_1,\bar {v}_2) \in\mathcal{B}_{\alpha_2} 处达到. 事实上, 这在 Sobolev-次临界情形是平凡的, 然而在临界情形可由引理 3.3 得证. 注意到 c_{\alpha_2} = \min\left\{\hat c_\alpha:\lambda_1(\Omega)\le \alpha \le \alpha_2\right\}\le\hat c_{\alpha_1} < \hat c_{\alpha_2}. 因此, (\bar{v}_1,\bar {v}_2)\in\mathcal{B} _{\alpha_2}\setminus \mathcal{U}_{\alpha_2}.

我们记 \lambda_2(\Omega)-\DeltaH^1_0(\Omega) 中的第二特征值, 并且 \varphi_2 是相应地特征函数.

引理3.5

(\sqrt{\rho_1}\varphi_1,\sqrt{\rho_2} \varphi_1)\in \mathcal{U} _{\lambda_1(\Omega)}, \ \ \left(\sqrt{\rho_1}\frac{\varphi_2^+}{|\varphi_2^+|_2},\sqrt{\rho_2} \frac{\varphi_2^-}{|\varphi_2^-|_2}\right)\in \mathcal{U} _{\lambda_2(\Omega)}.

特别地,

\hat{c}_{\lambda_{j}(\Omega)} \leq \frac{\rho_{1}+\rho_{2}}{2} \lambda_{j}(\Omega), \quad \text { 其中 } \quad j=\left\{\begin{array}{ll} 1, & \text { if } \beta \geq-\sqrt{\mu_{1} \mu_{2}}, \\ 2, & \text { if } \beta<-\sqrt{\mu_{1} \mu_{2}}. \end{array}\right.

通过直接计算, 第一个断言是直接的. 我们称

\begin{aligned} \hat{c}_{\lambda_{1}(\Omega)} \leq & \mathcal{E}_{1}\left(\sqrt{\rho_{1}} \varphi_{1}, \sqrt{\rho_{2}} \varphi_{1}\right) \\ = & \frac{\rho_{1}+\rho_{2}}{2} \lambda_{1}(\Omega)-\frac{1}{p+1} \int_{\Omega}\left(\mu_{1}\left|f\left(\sqrt{\rho_{1}} \varphi_{1}\right)\right|^{p+1}+2 \beta\left|f\left(\sqrt{\rho_{1}} \varphi_{1}\right)\right|^{\frac{p+1}{2}}\left|f\left(\sqrt{\rho_{2}} \varphi_{1}\right)\right|^{\frac{p+1}{2}}\right. \\ & \left.+\mu_{2}\left|f\left(\sqrt{\rho_{2}} \varphi_{1}\right)\right|^{p+1}\right) \leq \frac{\rho_{1}+\rho_{2}}{2} \lambda_{1}(\Omega). \end{aligned}

事实上, 由 \beta\ge-\sqrt{\mu_1\mu_2}, 对任意的 \rho_1,\rho_2>0, 使用引理 2.1, 我们讨论以下情形

情形 1 |\sqrt{\rho_1}\varphi_1|\leq1|\sqrt{\rho_2}\varphi_1|\leq1.

通过直接计算, 得到

\begin{aligned} & \mu_{1}\left|f\left(\sqrt{\rho_{1}} \varphi_{1}\right)\right|^{p+1}+2 \beta\left|f\left(\sqrt{\rho_{1}} \varphi_{1}\right)\right|^{\frac{p+1}{2}}\left|f\left(\sqrt{\rho_{2}} \varphi_{1}\right)\right|^{\frac{p+1}{2}}+\mu_{2}\left|f\left(\sqrt{\rho_{2}} \varphi_{1}\right)\right|^{p+1} \\ \geq & \theta\left(\mu_{1}\left|\sqrt{\rho_{1}} \varphi_{1}\right|^{p+1}+2 \beta\left|\sqrt{\rho_{1}} \varphi_{1}\right|^{\frac{p+1}{2}}\left|\sqrt{\rho_{2}} \varphi_{1}\right|^{\frac{p+1}{2}}+\mu_{2}\left|\sqrt{\rho_{2}} \varphi_{1}\right|^{p+1}\right) \\ \geq & \theta\left(\mu_{1}\left|\sqrt{\rho_{1}} \varphi_{1}\right|^{p+1}-2 \sqrt{\mu_{1} \mu_{2}}\left|\sqrt{\rho_{1}} \varphi_{1}\right|^{\frac{p+1}{2}}\left|\sqrt{\rho_{2}} \varphi_{1}\right|^{\frac{p+1}{2}}+\mu_{2}\left|\sqrt{\rho_{2}} \varphi_{1}\right|^{p+1}\right) \\ = & \theta\left(\sqrt{\mu_{1}}\left|\sqrt{\rho_{1}} \varphi_{1}\right|^{\frac{p+1}{2}}-\sqrt{\mu_{2}}\left|\sqrt{\rho_{2}} \varphi_{1}\right|^{\frac{p+1}{2}}\right)^{2} \geq 0. \end{aligned}

情形 2 |\sqrt{\rho_1}\varphi_1|\geq1|\sqrt{\rho_2}\varphi_1|\geq1.

通过直接计算, 得到

\begin{aligned} & \mu_{1}\left|f\left(\sqrt{\rho_{1}} \varphi_{1}\right)\right|^{p+1}+2 \beta\left|f\left(\sqrt{\rho_{1}} \varphi_{1}\right)\right|^{\frac{p+1}{2}}\left|f\left(\sqrt{\rho_{2}} \varphi_{1}\right)\right|^{\frac{p+1}{2}}+\mu_{2}\left|f\left(\sqrt{\rho_{2}} \varphi_{1}\right)\right|^{p+1} \\ \geq & \theta\left(\mu_{1}\left|\sqrt{\rho_{1}} \varphi_{1}\right|^{\frac{p+1}{2}}+2 \beta\left|\sqrt{\rho_{1}} \varphi_{1}\right|^{\frac{p+1}{4}}\left|\sqrt{\rho_{2}} \varphi_{1}\right|^{\frac{p+1}{4}}+\mu_{2}\left|\sqrt{\rho_{2}} \varphi_{1}\right|^{\frac{p+1}{2}}\right) \\ \geq & \theta\left(\mu_{1}\left|\sqrt{\rho_{1}} \varphi_{1}\right|^{\frac{p+1}{2}}-2 \sqrt{\mu_{1} \mu_{2}}\left|\sqrt{\rho_{1}} \varphi_{1}\right|^{\frac{p+1}{4}}\left|\sqrt{\rho_{2}} \varphi_{1}\right|^{\frac{p+1}{4}}+\mu_{2}\left|\sqrt{\rho_{2}} \varphi_{1}\right|^{\frac{p+1}{2}}\right) \\ = & \theta\left(\sqrt{\mu_{1}}\left|\sqrt{\rho_{1}} \varphi_{1}\right|^{\frac{p+1}{4}}-\sqrt{\mu_{2}}\left|\sqrt{\rho_{2}} \varphi_{1}\right|^{\frac{p+1}{4}}\right)^{2} \geq 0. \end{aligned}

情形 3 |\sqrt{\rho_1}\varphi_1|\leq1|\sqrt{\rho_2}\varphi_1|\geq1. 通过直接计算, 得到

\begin{aligned} & \mu_{1}\left|f\left(\sqrt{\rho_{1}} \varphi_{1}\right)\right|^{p+1}+2 \beta\left|f\left(\sqrt{\rho_{1}} \varphi_{1}\right)\right|^{\frac{p+1}{2}}\left|f\left(\sqrt{\rho_{2}} \varphi_{1}\right)\right|^{\frac{p+1}{2}}+\mu_{2}\left|f\left(\sqrt{\rho_{2}} \varphi_{1}\right)\right|^{p+1} \\ \geq & \theta\left(\mu_{1}\left|\sqrt{\rho_{1}} \varphi_{1}\right|^{p+1}+2 \beta\left|\sqrt{\rho_{1}} \varphi_{1}\right|^{\frac{p+1}{2}}\left|\sqrt{\rho_{2}}\right|^{\frac{p+1}{4}}\left|\varphi_{1}\right|^{\frac{p+1}{4}}+\mu_{2}\left|\sqrt{\rho_{2}} \varphi_{1}\right|^{\frac{p+1}{2}}\right) \\ \geq & \theta\left(\mu_{1}\left|\sqrt{\rho_{1}} \varphi_{1}\right|^{p+1}-2 \sqrt{\mu_{1} \mu_{2}}\left|\sqrt{\rho_{1}} \varphi_{1}\right|^{\frac{p+1}{2}}\left|\sqrt{\rho_{2}} \varphi_{1}\right|^{\frac{p+1}{4}}+\mu_{2}\left|\sqrt{\rho_{2}} \varphi_{1}\right|^{\frac{p+1}{2}}\right) \\ = & \theta\left(\sqrt{\mu_{1}}\left|\sqrt{\rho_{1}} \varphi_{1}\right|^{\frac{p+1}{2}}-\sqrt{\mu_{2}}\left|\sqrt{\rho_{2}} \varphi_{1}\right|^{\frac{p+1}{4}}\right)^{2} \geq 0. \end{aligned}

情形 4 |\sqrt{\rho_1}\varphi_1|\geq1|\sqrt{\rho_2}\varphi_1|\leq1. 通过直接计算, 得到

\begin{aligned} & \mu_{1}\left|f\left(\sqrt{\rho_{1}} \varphi_{1}\right)\right|^{p+1}+2 \beta\left|f\left(\sqrt{\rho_{1}} \varphi_{1}\right)\right|^{\frac{p+1}{2}}\left|f\left(\sqrt{\rho_{2}} \varphi_{1}\right)\right|^{p+1}+\mu_{2}\left|f\left(\sqrt{\rho_{2}} \varphi_{1}\right)\right|^{p+1} \\ \geq & \theta\left(\mu_{1}\left|\sqrt{\rho_{1}} \varphi_{1}\right|^{\frac{p+1}{2}}+2 \beta\left|\sqrt{\rho_{1}}\right|^{\frac{p+1}{4}}\left|\varphi_{1}\right|^{\frac{p+1}{4}}\left|\sqrt{\rho_{2}} \varphi_{1}\right|^{\frac{p+1}{2}}+\mu_{2}\left|\sqrt{\rho_{2}} \varphi_{1}\right|^{p+1}\right) \\ \geq & \theta\left(\mu_{1}\left|\sqrt{\rho_{1}} \varphi_{1}\right|^{\frac{p+1}{2}}-2 \sqrt{\mu_{1} \mu_{2}}\left|\sqrt{\rho_{1}} \varphi_{1}\right|^{\frac{p+1}{4}}\left|\sqrt{\rho_{2}} \varphi_{1}\right|^{\frac{p+1}{2}}+\mu_{2}\left|\sqrt{\rho_{2}} \varphi_{1}\right|^{p+1}\right) \\ = & \theta\left(\sqrt{\mu_{1}}\left|\sqrt{\rho_{1}} \varphi_{1}\right|^{\frac{p+1}{4}}-\sqrt{\mu_{2}}\left|\sqrt{\rho_{2}} \varphi_{1}\right|^{\frac{p+1}{2}}\right)^{2} \geq 0. \end{aligned}

另一方面, 由以上四种情形可知

\begin{aligned} \hat{c}_{\lambda_{2}(\Omega)} & \leq \mathcal{E}_{1}\left(\sqrt{\rho_{1}} \frac{\varphi_{2}^{+}}{\left|\varphi_{2}^{+}\right|_{2}}, \sqrt{\rho_{2}} \frac{\varphi_{2}^{-}}{\left|\varphi_{2}^{-}\right|_{2}}\right) \\ & =\frac{\rho_{1}+\rho_{2}}{2} \lambda_{2}(\Omega)-\frac{1}{p+1} \int_{\Omega}\left(\mu_{1}\left|f\left(\sqrt{\rho_{1}} \frac{\varphi_{2}^{+}}{\left|\varphi_{2}^{+}\right|_{2}}\right)\right|^{p+1}+\mu_{2}\left|f\left(\sqrt{\rho_{2}} \frac{\varphi_{2}^{-}}{\left|\varphi_{2}^{-}\right|_{2}}\right)\right|^{p+1}\right) \\ & \leq \frac{\rho_{1}+\rho_{2}}{2} \lambda_{2}(\Omega) \end{aligned}

引理3.6\rho_1,\rho_2>0 使得

\begin{equation}\label{eq:mainassL} \Lambda (\rho_1,\rho_2) \cdot (\rho_1+\rho_2)^{a-1} \le \frac{(a-1)^{a-1}}{a^a} \lambda_j(\Omega)^{-(a-1)}, \end{equation}
(3.14)

其中 j=1, \beta\ge-\sqrt{\mu_1\mu_2}, 否则, j=2.\bar \alpha= \frac{a}{a-1}\lambda_i(\Omega).c_{\bar{\alpha}}(\bar{v}_1,\bar{v}_2)\in\mathcal{B}_{\bar{\alpha}}\setminus \mathcal{U}_{\bar{\alpha}} 处达到, 满足 \mathcal{E}_1(\bar{v}_1,\bar{v}_2)=c_{\bar{\alpha}}, 这意味着 (\bar{v}_1,\bar{v}_2)\mathcal{E}_1 |_{\mathcal{M} } 的一个局部极小, 对应于 (1.10) 式的某个正解 (\omega_1,\omega_2)\in \mathbb R^2. 等价地, (\rho_1,\rho_2)\in A.

根据引理 3.4, 我们的目标是寻找 \bar{\alpha}>\lambda_j(\Omega) 使得

\begin{equation}\label{eq:c<hat_c_app} \hat{c}_{\lambda_j(\Omega)} < \hat{c}_{\bar{\alpha}}. \end{equation}
(3.15)

我们寻找满足 (3.15) 式的充分条件. 使用引理 3.2 和 3.5, 只需要寻找 \bar{\alpha}>\lambda_j(\Omega) 使得

\frac{\rho_1+\rho_2}{2} \lambda_j(\Omega) \le \frac12\left(\bar{\alpha}(\rho_1+\rho_2) -\Lambda (\rho_1,\rho_2) (\rho_1+\rho_2)^a\bar{\alpha}^a \right)

(由引理 3.5, 左边大于等于 \hat{c}_{\lambda_j(\Omega)}. 而且, 根据引理 3.2, 右边严格小于 \hat c_{\bar{\alpha}}).

\begin{equation}\label{eq:passa} \Lambda(\rho_1,\rho_2)(\rho_1+\rho_2)^{a-1}\le\frac{\bar{\alpha}- \lambda_j(\Omega)}{\bar{\alpha}^a}. \end{equation}
(3.16)

通过直接计算, 根据 p>3+\frac{4}{N} 可知 a=\frac{N(p-1)}{2(N+2)}>1, 对于右边最可能的选择是

\max_{\alpha\ge\lambda_j(\Omega)}\frac{\alpha-\lambda_j(\Omega) }{\alpha^a} = \frac{(a-1)^{a-1}}{a^a}\lambda_j(\Omega)^{-(a-1)}, \ \ \ \ \bar{\alpha}=\frac{a}{a-1}\lambda_j(\Omega).

\bar{\alpha} 的选取上是可能的, 因为 (3.16) 式等价于 (3.14) 式, 这是引理的假设. 显然 \bar{\alpha}>\lambda_j(\Omega). 因此, 为了应用引理 3.4 得出结论, 我们只需要证, 当 N\ge3p=\frac{3N+2}{N-2}

\begin{equation}\label{eq:passa2} \bar\alpha<\lambda_1(\Omega)+\frac{1}{\Lambda ^{\frac{N-2}{2}}(\rho_1+\rho_2)} \end{equation}
(3.17)

成立即可. 这是直接的, 因为 a=\frac{N}{N-2}, 根据 (3.16) 式可知

\Lambda(\rho_1+\rho_2)^{\frac{2}{N-2}}\le\frac{\bar{\alpha}- \lambda_j(\Omega)}{\bar{\alpha}^{\frac{N}{N-2}}} <\frac{\bar{\alpha}-\lambda_1(\Omega)}{(\bar{\alpha}-\lambda_1(\Omega))^{\frac{N}{N-2}}} =\frac{1}{(\bar{\alpha}- \lambda_1(\Omega))^{\frac{2}{N-2}}},

这等价于 (3.17) 式.

注3.2 如果 (\bar v_1,\bar v_2) 满足 (1.10) 式, 解 (\bar v_1,\bar v_2)(\sqrt{\rho_1} \varphi_1,\sqrt{\rho_2} \varphi_1) 不一致. 事实上, 它满足 (1.10) 式当且仅当对任意的 i=1,2

\begin{aligned} \lambda_{1}(\Omega) \sqrt{\rho_{i}} \varphi_{1} & =f^{\prime}\left(\sqrt{\rho_{i}} \varphi_{1}\right)\left(\left(\mu_{i}+\beta\right) f^{p-1}\left(\sqrt{\rho_{i}} \varphi_{1}\right)-\omega_{i}\right) f\left(\sqrt{\rho_{i}} \varphi_{1}\right) \\ & \Longleftrightarrow \lambda_{1}(\Omega)=f^{\prime}\left(\sqrt{\rho_{i}} \varphi_{1}\right)\left(\left(\mu_{i}+\beta\right) f^{p-1}\left(\sqrt{\rho_{i}} \varphi_{1}\right)-\omega_{i}\right) f\left(\sqrt{\rho_{i}} \varphi_{1}\right)\left(\sqrt{\rho_{i}} \varphi_{1}\right)^{-1} \\ & \Longleftrightarrow \lambda_{1}(\Omega)=\omega_{i}=0, \beta=-\mu_{1}=-\mu_{2} \end{aligned}

受文献[28,29]启发, 我们将研究非线性 Schrödinger 方程在 Sobolev-临界情形的方法推广到拟线性 Schrödinger 方程组. 下面可以明确一些量与 \rho_1,\rho_2 的相关性: 为了简化记号, 定义

G(v_1,v_2):=\int_\Omega\left(\mu_1|f(v_1)|^{p+1}+2\beta|f(v_1)|^{\frac{p+1}{2}} |f(v_2)|^{\frac{p+1}{2}}+\mu_2 |f(v_2)|^{p+1}\right),

并且引入最优化问题 M_\alpha(\rho_1,\rho_2):=\sup\limits_{\mathcal{U}_\alpha(\rho_1,\rho_2)}G. 注意到 \hat c_\alpha(\rho_1,\rho_2)\!=\! \frac{1}{2}\alpha(\rho_1+\rho_2)-\frac{1}{p+1}M_\alpha(\rho_1,\rho_2), 以及 \hat c_\alpha(\rho_1,\rho_2)(v_1,v_2)\in \mathcal{U}_\alpha(\rho_1,\rho_2) 处可达当且仅当 M_\alpha(\rho_1,\rho_2)(v_1,v_2) 处达到. 对任意的 (\rho_1,\rho_2)\in A\setminus\{(0,0)\}, 根据 A 的定义, 存在 \alpha>\lambda_1(\Omega) 以及 (\bar{v}_1,\bar{v}_2) \in\mathcal{B}_\alpha, (1.10) 式的一个解满足 \mathcal{E}_1 (\bar v_1,\bar v_2)=c_{\alpha}< \hat c_{\alpha}, 且当 N\ge3, p=\frac{3N+2}{N-2}\alpha 满足 (1.21) 式. 注意到假设 c_{\alpha}<\hat c_{\alpha} 意味着 \int_\Omega(|\nabla \bar{v}_1|^2+|\nabla\bar{v}_2|)^2<(\rho_1+\rho_2)\alpha, 因此

\bar{\alpha}:=\frac{1}{\rho_1+\rho_2}\int_\Omega(|\nabla \bar{v}_1|^2+|\nabla \bar{v}_2|^2)<\alpha.

从而, (\bar v_1,\bar v_2)\in\mathcal{U}_{\bar{\alpha}}(\rho_1,\rho_2) 达到 \hat c_{\bar{\alpha}}= c_\alpha.

引理3.7 如果 s>0, 则 (s \bar v_1,s \bar v_2)\in\mathcal{U}_{\bar{\alpha}} (s^2\rho_1,s^2\rho_2) 达到

\hat c_{\bar \alpha}(s^2 \rho_1,s^2\rho_2)=\frac{s^2}2\bar{\alpha} (\rho_1+\rho_2)-\frac{1}{p+1}G(s\bar v_1,s\bar v_2 ).

通过直接计算, (v_1,v_2)\in\mathcal{U}_{\bar\alpha}(\rho_1,\rho_2)\Leftrightarrow (sv_1,sv_2)\in\mathcal{U}_{\bar\alpha}(s^2\rho_1,s^2\rho_2),

G(sv_1,sv_2)=\int_\Omega\left(\mu_1|f(sv_1)|^{p+1}+2\beta|f(sv_1)|^{\frac{p+1}{2}} |f(sv_2)|^{\frac{p+1}{2}}+\mu_2 |f(sv_2)|^{p+1}\right).

引理3.8s\in(0,1), |t| 充分小,

(U_1(t),U_2(t)):=\bigg(f^{-1}\left(s\sqrt{\rho_1} \frac{f(\bar v_1+t\tau _1)}{|f(\bar v_1+t\tau _1|_2)}\right), f^{-1}\left(s\sqrt{\rho_2} \frac{f(\bar v_2+t\tau _2)}{|f(\bar v_2+t\tau _2)|_2}\right)\bigg),

其中 (\tau_1,\tau_2)\in H^1_0(\Omega,\mathbb R^2) 使得

\begin{array}{l} \int_{\Omega} \frac{1}{1+f^{2}\left(\bar{v}_{1}\right)}\left(\nabla \bar{v}_{1} \nabla \tau_{1}\left(1+s^{2} f^{2}\left(\bar{v}_{1}\right)\right)-\frac{\left|\nabla \bar{v}_{1}\right|^{2}\left(1+2 s^{2} f^{2}\left(\bar{v}_{1}\right)\right)}{\rho_{1}} \int_{\Omega} \frac{f\left(\bar{v}_{1}\right) \tau_{1}}{\sqrt{1+\left|\bar{v}_{1}\right|^{2}}}\right. \\ \left.-\frac{f\left(\bar{v}_{1}\right) \tau_{2}\left|\nabla \bar{v}_{1}\right|^{2}\left(1+s^{2} f^{2}\left(\bar{v}_{1}\right)\right)}{\left(1+f^{2}\left(\bar{v}_{1}\right)\right) \sqrt{1+\left|\bar{v}_{1}\right|^{2}}}\right) \\ +\int_{\Omega} \frac{1}{1+f^{2}\left(\bar{v}_{2}\right)}\left(\nabla \bar{v}_{2} \nabla \tau_{2}\left(1+s^{2} f^{2}\left(\bar{v}_{2}\right)\right)-\frac{\left|\nabla \bar{v}_{2}\right|^{2}\left(1+2 s^{2} f^{2}\left(\bar{v}_{2}\right)\right)}{\rho_{2}} \int_{\Omega} \frac{f\left(\bar{v}_{2}\right) \tau_{2}}{\sqrt{1+\left|\bar{v}_{2}\right|^{2}}}\right. \\ \left.-\frac{f\left(\bar{v}_{2}\right) \tau_{i}\left|\nabla \bar{v}_{2}\right|^{2}\left(1+s^{2} f^{2}\left(\bar{v}_{2}\right)\right)}{\left(1+f^{2}\left(\bar{v}_{2}\right)\right) \sqrt{1+\left|\bar{v}_{2}\right|^{2}}}\right)<0, \quad i=1,2 \end{array}
\int_{\Omega}\left(\frac{\mu_{1}\left|f\left(\bar{v}_{1}\right)\right|^{p} \tau_{1}}{\sqrt{1+f^{2}\left(\bar{v}_{1}\right)}}+\frac{\beta\left|f\left(\bar{v}_{1}\right)\right|^{\frac{p-1}{2}} \tau_{1}}{\sqrt{1+f^{2}\left(\bar{v}_{1}\right)}}\left|f\left(\bar{v}_{2}\right)\right|^{\frac{p+1}{2}}+\frac{\beta\left|f\left(\bar{v}_{2}\right)\right|^{\frac{p-1}{2}} \tau_{2}}{\sqrt{1+f^{2}\left(\bar{v}_{2}\right)}}\left|f\left(\bar{v}_{1}\right)\right|^{\frac{p+1}{2}}+\frac{\mu_{2}\left|f\left(\bar{v}_{2}\right)\right|^{p}}{\sqrt{1+f^{2}\left(\bar{v}_{2}\right)}} \tau_{2}\right)>0.

则对任意的 t,

\frac{\rm d}{{\rm d}t}\|(U_1(t),U_2(t))\|^2_{H^1_0(\Omega)}\bigg|_{t=0}<0, \ \ \ \ \frac{\rm d}{{\rm d}t} \mathcal{E}_1 (U_1(t),U_2(t))\bigg|_{t=0}<0,

(U_1(t),U_2(t))\in\mathcal{M}_{s^2\rho_1,s^2\rho_2}.

经过简单计算,

\int_{\Omega}f^2\left(f^{-1}\left(s\sqrt{\rho_i} \frac{f(\bar v_i+t\tau_i)}{|f(\bar v_i+t\tau_i|_2)}\right)\right) =\int_{\Omega}\left(s\sqrt{\rho_i} \frac{f(\bar v_i+t\tau_i)}{|f(\bar v_i+t\tau_i)|_2}\right)^2=s^2\rho_i,\ \ i=1,2,

对任意的 t, (U_1(t),U_2(t))\in\mathcal{M} _{s^2\rho_1,s^2\rho_2}. 根据 (2.3) 式, 计算

\begin{array}{l} \left.\frac{\mathrm{d}}{\mathrm{~d} t}\left(\int_{\Omega}\left|\nabla\left(f^{-1}\left(s \sqrt{\rho_{i}} \frac{f\left(\bar{v}_{i}+t \tau_{i}\right)}{\mid f\left(\bar{v}_{i}+\left.t \tau_{i}\right|_{2}\right)}\right)\right)\right|^{2}\right)\right|_{t=0} \\ =\left.\frac{\mathrm{d}}{\mathrm{~d} t}\left(\int_{\Omega}\left|\sqrt{1+s^{2} \rho_{i} \frac{f^{2}\left(f\left(\bar{v}_{i}+t \tau_{i}\right)\right)}{\left|f\left(\bar{v}_{i}+t \tau_{i}\right)\right|_{2}^{2}}} \frac{s \sqrt{\rho_{i}} f^{\prime}\left(\bar{v}_{i}+t \tau_{i}\right) \nabla\left(\bar{v}_{i}+t \tau_{i}\right)}{\left|f\left(\bar{v}_{i}+t \tau_{i}\right)\right|_{2}}\right|^{2}\right)\right|_{t=0} \\ =2 s^{2}\left(\int _ { \Omega } \frac { 1 } { 1 + f ^ { 2 } ( \overline { v } _ { i } ) } \left(\nabla \bar{v}_{i} \nabla \tau_{i}\left(1+s^{2} f^{2}\left(\bar{v}_{i}\right)\right)-\frac{\left|\nabla \bar{v}_{i}\right|^{2}\left(1+2 s^{2} f^{2}\left(\bar{v}_{i}\right)\right)}{\rho_{i}} \int_{\Omega} \frac{f\left(\bar{v}_{i}\right) \tau_{i}}{\sqrt{1+\left|\bar{v}_{i}\right|^{2}}}\right.\right. \\ \left.\left.-\frac{f\left(\bar{v}_{i}\right) \tau_{i}\left|\nabla \bar{v}_{i}\right|^{2}\left(1+s^{2} f^{2}\left(\bar{v}_{i}\right)\right)}{\left(1+f^{2}\left(\bar{v}_{i}\right)\right) \sqrt{1+\left|\bar{v}_{i}\right|^{2}}}\right)\right), \quad i=1,2. \end{array}

\begin{aligned} & \left.\frac{\mathrm{d}}{\mathrm{~d} t}\left\|\left(U_{1}(t), U_{2}(t)\right)\right\|_{H_{0}^{1}(\Omega)}^{2}\right|_{t=0} \\ = & 2 s^{2}\left(\int _ { \Omega } \frac { 1 } { 1 + f ^ { 2 } ( \overline { v } _ { 1 } ) } \left(\nabla \bar{v}_{1} \nabla \tau_{1}\left(1+s^{2} f^{2}\left(\bar{v}_{1}\right)\right)-\frac{\left|\nabla \bar{v}_{1}\right|^{2}\left(1+2 s^{2} f^{2}\left(\bar{v}_{1}\right)\right)}{\rho_{1}} \int_{\Omega} \frac{f\left(\bar{v}_{1}\right) \tau_{1}}{\sqrt{1+\left|\bar{v}_{1}\right|^{2}}}\right.\right. \\ & \left.-\frac{f\left(\bar{v}_{1}\right) \tau_{1}\left|\nabla \bar{v}_{1}\right|^{2}\left(1+s^{2} f^{2}\left(\bar{v}_{1}\right)\right)}{\left(1+f^{2}\left(\bar{v}_{1}\right)\right) \sqrt{1+\left|\bar{v}_{1}\right|^{2}}}\right)+\int_{\Omega} \frac{1}{1+f^{2}\left(\bar{v}_{2}\right)}\left(\nabla \bar{v}_{2} \nabla \tau_{2}\left(1+s^{2} f^{2}\left(\bar{v}_{2}\right)\right)\right. \\ & \left.\left.-\frac{\left|\nabla \bar{v}_{2}\right|^{2}\left(1+2 s^{2} f^{2}\left(\bar{v}_{2}\right)\right)}{\rho_{i}} \int_{\Omega} \frac{f\left(\bar{v}_{2}\right) \tau_{2}}{\sqrt{1+\left|\bar{v}_{2}\right|^{2}}}-\frac{f\left(\bar{v}_{2}\right) \tau_{i}\left|\nabla \bar{v}_{2}\right|^{2}\left(1+s^{2} f^{2}\left(\bar{v}_{2}\right)\right)}{\left(1+f^{2}\left(\bar{v}_{2}\right)\right) \sqrt{1+\left|\bar{v}_{2}\right|^{2}}}\right)\right)<0. \end{aligned}

而且, 根据

\begin{aligned} \left.\frac{\mathrm{d}}{\mathrm{~d} t}\left(\int_{\Omega} \mu_{i}\left|f\left(U_{i}(t)\right)\right|^{p+1}\right)\right|_{t=0} & =\left.\frac{\mathrm{d}}{\mathrm{~d} t}\left(\int_{\Omega} \mu_{i}\left|s \sqrt{\rho_{i}} \frac{f\left(\bar{v}_{i}+t \tau_{i}\right)}{\left|f\left(\bar{v}_{i}+t \tau_{i}\right)\right|_{2}}\right|^{p+1}\right)\right|_{t=0} \\ & =(p+1) s^{p+1} \int_{\Omega} \mu_{i} \frac{\left|f\left(\bar{v}_{i}\right)\right|^{p}}{\sqrt{1+f^{2}\left(\bar{v}_{i}\right)}} \tau_{i}, \quad i=1,2 \end{aligned}

可知

\begin{aligned} & \left.\frac{\mathrm{d}}{\mathrm{~d} t} G\left(U_{1}(t), U_{2}(t)\right)\right|_{t=0} \\ = & \left.\frac{\mathrm{d}}{\mathrm{~d} t}\left(\int_{\Omega}\left(\mu_{1}\left|f\left(U_{1}(t)\right)\right|^{p+1}+2 \beta\left|f\left(U_{1}(t)\right)\right|^{\frac{p+1}{2}}\left|f\left(U_{2}(t)\right)\right|^{\frac{p+1}{2}}+\mu_{2}\left|f\left(U_{2}(t)\right)\right|^{p+1}\right)\right)\right|_{t=0} \\ = & (p+1) s^{p+1} \int_{\Omega}\left(\mu_{1} \frac{\left|f\left(\bar{v}_{1}\right)\right|^{p} \tau_{1}}{\sqrt{1+f^{2}\left(\bar{v}_{1}\right)}}+\beta \frac{\left|f\left(\bar{v}_{1}\right)\right|^{\frac{p-1}{2}} \tau_{1}}{\sqrt{1+f^{2}\left(\bar{v}_{1}\right)}\left|f\left(\bar{v}_{2}\right)\right|^{\frac{p+1}{2}}}\right. \\ & \left.+\beta \frac{\left|f\left(\bar{v}_{2}\right)\right|^{\frac{p-1}{2}} \tau_{2}}{\sqrt{1+f^{2}\left(\bar{v}_{2}\right)}}\left|f\left(\bar{v}_{1}\right)\right|^{\frac{p+1}{2}}+\mu_{2} \frac{\left|f\left(\bar{v}_{2}\right)\right|^{p}}{\sqrt{1+f^{2}\left(\bar{v}_{2}\right)}} \tau_{2}\right)>0 \end{aligned}

因此, 对任意的 0<s<1p>1,

\begin{align*} &~~~~ \frac{\rm d}{{\rm d}t}\mathcal{E}_1(U_1(t),U_2(t))\bigg|_{t=0}\\ &=s^2\bigg(\int_\Omega \frac{1}{1+f^2(\bar v_1)}\bigg(\nabla \bar v_1\nabla \tau_1(1+s^2f^2(\bar v_1)) -\frac{|\nabla \bar v_1|^2(1+2s^2f^2(\bar v_1))}{\rho_1} \int_\Omega \frac{f(\bar v_1)\tau_1}{\sqrt{1+|\bar v_1|^2}}\\ &\ \ \ \ -\frac{f(\bar v_1)\tau_1|\nabla \bar v_1|^2(1+s^2f^2(\bar v_1)) }{(1+f^2(\bar v_1)) \sqrt{1+|\bar v_1|^2}} \bigg)+\int_\Omega \frac{1}{1+f^2(\bar v_2)}\bigg(\nabla \bar v_2\nabla \tau_2(1+s^2f^2(\bar v_2))\\ &\ \ \ \ -\frac{|\nabla \bar v_2|^2(1+2s^2f^2(\bar v_2))}{\rho_i} \int_\Omega \frac{f(\bar v_2)\tau_2}{\sqrt{1+|\bar v_2|^2}} -\frac{f(\bar v_2)\tau_i|\nabla \bar v_2|^2(1+s^2f^2(\bar v_2)) }{(1+f^2(\bar v_2)) \sqrt{1+|\bar v_2|^2}} \bigg)\bigg)\\ &\ \ \ \ -s^{p+1}\int_\Omega\bigg(\mu_1\frac{|f(\bar v_1)|^p \tau_1}{\sqrt{1+f^2(\bar v_1)} } +\beta\frac{|f(\bar v_1)|^{\frac{p-1}{2} }\tau_1}{\sqrt{1+f^2(\bar v_1)}} |f(\bar v_2)|^{\frac{p+1}{2} }\\ &\ \ \ \ +\beta\frac{|f(\bar v_2)|^{\frac{p-1}{2} }\tau_2}{\sqrt{1+f^2(\bar v_2)}} |f(\bar v_1)|^{\frac{p+1}{2} } +\mu_2\frac{|f(\bar v_2)|^p}{\sqrt{1+f^2(\bar v_2)} }\tau_2\bigg) <0. \end{align*}

引理3.9A 如 (1.23) 定义所示. 则 A 对于 (0,0) 是星型的.

根据引理 3.8

(U_1(0),U_2(0))\!=\!\bigg(f^{-1}\left(s\sqrt{\rho_1} \frac{f(\bar v_1)}{|f(\bar v_1|_2)}\right), f^{-1}\left(s\sqrt{\rho_2} \frac{f(\bar v_2)}{|f(\bar v_2)|_2}\right)\bigg) \!=\!(s\bar v_1,s\bar v_2)\!\in\mathcal{U}_{\bar\alpha} (s^2\rho_1,s^2\rho_2),

存在正的和小的常数 \varepsilon,\tau 使得 (U_1(\tau ),U_2(\tau))\in\mathcal{U} _{\bar\alpha-\varepsilon}(s^2\rho_1,s^2\rho_2). 结合引理 3.7 和 3.8 我们得到

\hat{c}_{\bar{\alpha}-\varepsilon}(s^2\rho_1,s^2\rho_2) \le \mathcal{E}_1 (U_1(\tau),U_2(\tau))<\mathcal{E}_1(U_1(0),U_2(0))=\hat{c}_{\bar{\alpha}} (s^2\rho_1,s^2\rho_2)

应用引理 3.4, \alpha_1=\bar\alpha-\varepsilon\alpha_2=\bar\alpha 可知对任意的 s\in(0,1), (s^2\rho_1,s^2\rho_2)\in A.

接下来, 我们估计 \Lambda . 在这点上, 引理 3.6 的主要假设对于函数 \Lambda(\rho_1,\rho_2) 定义在 (3.4)式. 为了读者方便, 回顾

\Lambda (\rho_1,\rho_2)=\frac{8C_{N,p}}{p+1} \max_{t\in \left[\frac{\pi}{2}\right]} \left(\mu_1 \rho_1^{2r}\cos^{2a}t +\mu_2 \rho_2^{2r}\sin^{2a}t +2\beta^+ \rho_1^r\rho_2^{r} \cos^{a}t \sin^{a}t\right),

其中

a=\frac{N(p-1)}{2(N+2)}\in\left(1,\frac{N}{N-2}\right], \ \ r = \frac{4N-(N-2)(p+1)}{N+2}\in\left[0,\frac{1}{N}\right).

显然, \Lambda 是一个关于 (\rho_1,\rho_2)r-齐次多项式, 但它的显式表达只能在少数特定情况下推导出来. 接下来证明定理 1.4 中的条件 (1.26), R=R(\Omega,N,p) 的定义为

\begin{equation}\label{eq:def_R} R(\Omega,N,p)=\frac{p+1}{8C_{N,p}} \frac{(a-1)^{a-1}}{a^a} \lambda_j(\Omega)^{-(a-1)}, \end{equation}
(3.18)

意味着引理 3.6 中的假设 (3.14) 式成立.

定理 1.4 的证明 引理 3.9 可知 A 对于 (0,0) 是星型的. 我们从上面估计 \Lambda , 由 a>1

\begin{aligned} \Lambda\left(\rho_{1}, \rho_{2}\right) & \leq \Lambda^{\prime}\left(\rho_{1}, \rho_{2}\right) \\ : & =\frac{8 C_{N, p}}{p+1} \max _{t \in\left[\frac{\pi}{2}\right]}\left(\mu_{1} \rho_{1}^{2 r} \cos ^{2} t+\mu_{2} \rho_{2}^{2 r} \sin ^{2} t+2 \beta^{+} \rho_{1}^{r} \rho_{2}^{r} \cos t \sin t\right) \\ & =\frac{4 C_{N, p}}{p+1} \max _{t \in\left[\frac{\pi}{2}\right]}\left(\left(\mu_{1} \rho_{1}^{2 r}+\mu_{2} \rho_{2}^{2 r}\right)+\left(\mu_{1} \rho_{1}^{2 r}-\mu_{2} \rho_{2}^{2 r}\right) \cos 2 t+2 \beta^{+} \rho_{1}^{r} \rho_{2}^{r} \sin 2 t\right) \\ & =\frac{4 C_{N, p}}{p+1} \max _{x^{2}+y^{2}=1}\left(\left(\mu_{1} \rho_{1}^{2 r}+\mu_{2} \rho_{2}^{2 r}\right)+\left(\mu_{1} \rho_{1}^{2 r}-\mu_{2} \rho_{2}^{2 r}\right) x+2 \beta^{+} \rho_{1}^{r} \rho_{2}^{r} y\right) \end{aligned}

通过计算得出

\begin{aligned} \Lambda^{\prime}\left(\rho_{1}, \rho_{2}\right) & =\frac{4 C_{N, p}}{p+1}\left(\mu_{1} \rho_{1}^{2 r}+\mu_{2} \rho_{2}^{2 r}+\sqrt{\left(\mu_{1} \rho_{1}^{2 r}-\mu_{2} \rho_{2}^{2 r}\right)^{2}+4\left(\beta^{+} \rho_{1}^{r} \rho_{2}^{r}\right)^{2}}\right) \\ & \leq \frac{4 C_{N, p}}{p+1}\left(\mu_{1} \rho_{1}^{2 r}+\mu_{2} \rho_{2}^{2 r}+\left|\mu_{1} \rho_{1}^{2 r}-\mu_{2} \rho_{2}^{2 r}\right|+2 \beta^{+} \rho_{1}^{r} \rho_{2}^{r}\right) \\ & =\frac{8 C_{N, p}}{p+1}\left(\max \left\{\mu_{1} \rho_{1}^{2 r}, \mu_{2} \rho_{2}^{2 r}\right\}+\beta^{+} \rho_{1}^{r} \rho_{2}^{r}\right). \end{aligned}

因此, 根据 (1.26) 和 (3.18) 式, 我们得到 (3.14) 式, 从而应用引理 3.6 得出结论.

注3.3 与 (3.14) 式相比, (1.26) 式更精确; 而且, 两个条件与 \beta\leq0 一致. 事实上, 当 t=0t=\frac{\pi}{2}, 最大值原理知 \Lambda' 是可达的且

\Lambda(\rho_1,\rho_2)=\Lambda '(\rho_1,\rho_2)= \frac{8C_{N,p}}{p+1}\max\{\mu_1\rho_1^{2r},\mu_2\rho_2^{2r}\}\ \ \text{其中}\ \ \beta\leq0.

4 正规化基态解的轨道稳定性

本节中, 我们主要证明正规化基态解的轨道稳定性, 即定理 1.1(b) 和 1.5. 具体地, 我们的目标是证明分别定义在 (1.18) 和 (1.27) 式的 GG_{\bar\alpha} 的稳定性.

注意到全局极小值也是局部极小值, 我们为所有情况提供了统一的证明. 回顾 (1.19)式中的定义 \mathcal{B}_\alpha,\ \mathcal{U}_\alpha, (1.20) 式中的定义 c_\alpha,\ \hat{c}_\alpha. 一方面, 对于 p\leq3+\frac{4}{N}(\rho_1,\rho_2) 满足定理的假设 1.1, 根据 (2.9) 以及 (2.10) 式, 我们可以推导出存在 \bar\alpha>\lambda_1(\Omega) 使得

\{(|v_1|,|v_2|)\in H^1_0(\Omega;\mathbb R^2):(v_1,v_2)\in G\}\subseteq\mathcal{B}_{\bar\alpha},
\inf_{\mathcal{M}}\mathcal{E}_1=c_{\bar{\alpha}}<\hat {c}_{\bar{\alpha}}.

特别地, G=G_{\bar{\alpha}}. 另一方面, 对 p>3+\frac{4}{N}(\rho_1,\rho_2)\in A, 取 \bar{\alpha}\ge\lambda_1(\Omega) 使得 c_{\bar\alpha}<\hat c_{\bar\alpha}, 根据定理 1.5, 在 p=\frac{3N+2}{N-2} 时满足 (1.21) 式. 因此, 通过取 \bar{\alpha}>\lambda_1(\Omega) 证明 G_{\bar\alpha} 的稳定性.

受文献[[28],第 4 节] 启发, 我们称集合 \mathcal{G}\subset H^1_0(\Omega;\mathbb C^2) 是轨道稳定的, 如果对任意的 \varepsilon >0 存在 \delta>0 使得 (\psi_1,\psi_2)\in H^1_0(\Omega;\mathbb C^2) 满足 \text{dist}_{H^1_0}((\psi_1,\psi_2),\mathcal{G})<\delta, \text{dist}_{H^1_0} 表示 {H^1_0}-距离, 因此

\left\{\begin{array}{l} \mathrm{i} \partial_{t} \Psi_{1}+\Delta \Psi_{1}=f^{\prime}\left(\Psi_{1}\right)\left(\mu_{1}\left|f\left(\Psi_{1}\right)\right|^{p-1} f\left(\Psi_{1}\right)+\beta\left|f\left(\Psi_{2}\right)\right|^{\frac{p+1}{2}}\left|f\left(\Psi_{1}\right)\right|^{\frac{p-3}{2}} f\left(\Psi_{1}\right)\right)=0 \\ \mathrm{i} \partial_{t} \Psi_{2}+\Delta \Psi_{2}=f^{\prime}\left(\Psi_{2}\right)\left(\mu_{2}\left|f\left(\Psi_{2}\right)\right|^{p-1} f\left(\Psi_{2}\right)+\beta\left|f\left(\Psi_{1}\right)\right|^{\frac{p+1}{2}}\left|f\left(\Psi_{2}\right)\right|^{\frac{p-3}{2}} f\left(\Psi_{2}\right)\right)=0 \\ \Psi_{i}(0, \cdot)=\psi_{i}(\cdot), \quad \Psi_{i}(t, \cdot) \in H_{0}^{1}\left(\Omega ; \mathbb{C}^{2}\right) \end{array}\right.

的解 (\Psi_{1},\Psi_{2}) 满足对 0\leq t<+\infty,

\begin{equation}\label{eq:os1} (\Psi_1(t,\cdot),\Psi_2(t,\cdot))\ \ \text{ 可以连续到一个解} \end{equation}
(4.1)

\begin{equation}\label{eq:os2} \sup\limits_{t>0}\text{dist}_{H^1_0}((\Psi_1(t,\cdot),\Psi_2(t,\cdot)),\mathcal{G})<\varepsilon, \end{equation}
(4.2)

其中 \Psi_i=f^{-1}(\Phi _i),\ i=1,2.

引理4.1 对于 (v_1,v_2)\in G_{\bar\alpha}, 存在 \theta_1,\theta_2\in\mathbb R 使得 (v_1,v_2)=({\rm e}^{{\rm i}\omega \theta_1}|v_1|,{\rm e}^{{\rm i}\omega \theta_2}|v_2|). 特别地,

\inf \left\{\mathcal{E}_1 (v_1,v_2): (v_1,v_2) \in H^1_0(\Omega;\mathbb C^2), (|v_1|,|v_2|) \in\mathcal{B} _{\bar\alpha}\right\}=c_{\bar\alpha},

然而

\inf \left\{\mathcal{E}_1 (v_1,v_2): (v_1,v_2) \in H^1_0(\Omega;\mathbb C^2), (|v_1|,|v_2|) \in\mathcal{U} _{\bar\alpha}\right\} =: \tilde c_{\bar \alpha}\le \hat c_{\bar \alpha}.

给定 (\tau_1,\tau_2)\in \mathcal{U}_{\bar\alpha}, 根据 \mathcal{U}_{\bar\alpha} 的定义, 我们知道

\int_\Omega(|\nabla \tau_1|^2+|\nabla \tau_2|^2)x=(\rho_1+\rho_2)\bar\alpha,

易知 (\tau_1,\tau_2) \in H^1_0(\Omega;\mathbb C^2)(|\tau_1|,|\tau_2|) \in\mathcal{U}_{\bar\alpha}, 因此 \tilde c_{\bar\alpha}\le\hat c_{\bar\alpha}. 对于 (v_1,v_2)\in G_{\bar\alpha}, 通过使用 diamagnetic 不等式[24,定理 7.21], 得到

\int_\Omega|\nabla|v_i||^2\le\int_\Omega|\nabla v_i|^2,\ \ i=1,2.

结合 (1.8) 式, 可得到 f(|v_i|)=f(v_i),\ i=1,2

c_{\bar\alpha}\le\mathcal{E}_1(|v_1|,|v_2|)\le\mathcal{E}_1(v_1,v_2)=c_{\bar\alpha}.

因此, 对 i=1,2,

\int_\Omega|\nabla|v_i||^2=\int_\Omega|\nabla v_i|^2,

因此不等式成立, 从而 v_i|v_i| 的复数, 即 v_i={\rm e}^{{\rm i}\omega\theta_i}|v_i| 对某个 \theta_i\in\mathbb R 成立, 从而引理得证.

引理4.2 假设 \bar \alpha 如上定义. 对 \{(\psi_{1,n},\psi_{2,n})\}\subset H^1_0(\Omega;\mathbb C^2) 满足

\begin{equation}\label{eq:stability_compact101} \int_\Omega|\psi_{i,n}|^2\to\rho_i,\ \ \ i=1,2,\ \ \mathcal{E}_1(\psi_{1,n},\psi_{2,n})\to c_{\bar\alpha}\ \ \text{当}\ \ n\to\infty, \end{equation}
(4.3)

并且对 n 充分大

\begin{equation}\label{eq:stability_compact2} \int_\Omega(|\nabla\psi_{1,n}|^2+|\nabla \psi_{2,n}|)^2\leq (\rho_1+\rho_2)\bar \alpha +\text{o}(1). \end{equation}
(4.4)

则存在 (v_1,v_2)\in G_{\bar \alpha} 使得在 H^1_0(\Omega;\mathbb C^2)(\psi_{1,n},\psi_{2,n})\to (v_1,v_2).

通过使用 (4.4) 式, 知道存在 (\bar\psi_1,\bar\psi_2)\in H^1_0(\Omega;\mathbb C^2) 使得在 H^1_0(\Omega;\mathbb C)\psi_{i,n}\rightharpoonup v_i 且在 L^2(\Omega;\mathbb C)\psi_{i,n}\to v_i, i=1,2. 因此, 根据 (4.3), (4.4) 式和引理 4.1 可知, 对 i=1,2, 有

\int_\Omega|v_i|^2 =\rho_i,\ \ \int_\Omega(|\nabla v_1|^2+|\nabla v_2|^2)\leq(\rho_1+\rho_2)\bar\alpha,\ \ \mathcal{E}_1(v_{1},v_{2})\ge c_{\bar\alpha}.

一方面, 当 p<\frac{3N+2}{N-2} , Sobolev 嵌入的紧性意味着 (\psi_{1,n},\psi_{2,n})\to(v_1,v_2)L^{p+1}(\Omega;\mathbb C^2) 中成立. 因此, 根据引理 2.1 和 2.3 可知

\int_\Omega|f(\psi_{i,n})|^{p+1}=|f(v_{i})|^{p+1},\ \ \int_\Omega|f(\psi_{1,n})|^{\frac{p+1}{2} }|f(\psi_{2,n})|^{\frac{p+1}{2} } =|f(v_{1})|^{\frac{p+1}{2}}|f(v_{2})|^{\frac{p+1}{2}},\ \ i=1,2.

由范数的弱下半连续性

c_{\bar\alpha}\le\mathcal{E}_1(v_{1},v_{2})\le\liminf_{n\to+\infty} \mathcal{E}_1(\psi_{1,n},\psi_{2,n})=c_{\bar\alpha},

且结合 (v_1,v_2)\in G_{\bar\alpha} 知在 H^1_0 中强收敛,.

另一方面, 当 p=\frac{3N+2}{N-2}, 根据引理 3.3 可知该结果成立: 事实上, 引理对实值函数成立, 但在引理 4.1 之后, 可以直接验证其证明也适用于复值函数.

接下来, 我们给出稳定性的证明.

引理4.3 \bar \alpha 如上所述. 如果 c_{\bar \alpha}<\tilde c_{\bar \alpha}, 则 G_{\bar \alpha} 是(条件)轨道稳定的.

用反证法, 假设 \{(\psi_{1,n},\psi_{2,n})\}\subset H^1_0 (\Omega;\mathbb C^2), (v_{1,n},v_{2,n})\in G_{\bar \alpha}\bar \varepsilon >0 使得

\begin{equation}\label{eq:psi_to_u} \lim_{n\to\infty} \|(\psi_{1,n},\psi_{2,n})-(v_{1,n},v_{2,n})\|_{H^1_0(\Omega;\mathbb C^2)}=0, \end{equation}
(4.5)

\begin{equation} \sup\limits_{t>0}\text{dist}_{H^1_0}((\Psi_{1,n}(t,\cdot),\Psi_{2,n}(t,\cdot)),G_{\bar \alpha}) \geq 2\bar\varepsilon, \end{equation}
(4.6)

其中 (\Psi_{1,n},\Psi_{2,n}) 是 (1.2) 式带初始条件 (\psi_{1,n},\psi_{2,n}) 的一个解. 则存在 \{t_n\} 使得, 取 \phi_{i,n}(x):=\Psi_{i,n}(t_n,x), i=1,2

\begin{equation}\label{eq:stability_phi} \text{dist}_{H^1_0}((\phi_{1,n},\phi_{2,n}),G_{\bar \alpha}) \geq\bar\varepsilon. \end{equation}
(4.7)

接下来, 我们证明 \{(\phi_{1,n},\phi_{2,n})\} 满足 (4.3) 和 (4.4) 式. 则引理 4.2 与 (4.7) 式矛盾.

引理 4.2 意味着 G_{\bar \alpha} 是紧的. 因此, 通过使用 (4.5) 式, 存在 (v_1,v_2)\in G_{\bar \alpha} 使得

\begin{equation}\label{eq:strongconvergence_stability} (\psi_{1,n},\psi_{2,n})\to (v_{1},v_{2}) \ \ \text{在}\ \ H^1_0(\Omega;\mathbb C^2). \end{equation}
(4.8)

结合 Sobolev 嵌入的连续性, 根据引理 2.1 和 2.3 可知

\int_\Omega|f(\psi_{i,n})|^{p+1}=|f(v_{i})|^{p+1},\ \ \int_\Omega|f(\psi_{1,n})|^{\frac{p+1}{2} }|f(\psi_{2,n})|^{\frac{p+1}{2} } =|f(v_{1})|^{\frac{p+1}{2}}|f(v_{2})|^{\frac{p+1}{2}},\ \ i=1,2,

这意味着 (\psi_{1,n},\psi_{2,n}) 满足 (4.3) 式. 因此质量和能量守恒可知

\int_\Omega|\phi_{i,n}|^2=\int_\Omega|\psi_{i,n}|^2\to\rho_i,\ \ i=1,2, \ \ \text{和}\ \ \mathcal{E}_1(\phi_{1,n},\phi_{2,n}) =\mathcal{E}_1(\psi_{1,n},\psi_{2,n})\to c_{\bar\alpha},

因此 (\phi_{1,n},\phi_{2,n}) 也满足 (4.3) 式. 接下来, 我们验证, 至少有一个序列 (\phi_{1,n}, \phi_{2,n}) 满足 (4.4) 式, 即, 对充分大的 n,

\begin{equation}\label{eq:stability_contr2} \int_\Omega (|\nabla \phi_{1,n}|^2 + |\nabla \phi_{2,n}|^2) \leq(\rho_1+\rho_2) \bar\alpha + \text{o}(1). \end{equation}
(4.9)

事实上, 通过矛盾, 假设存在 \bar n\in\mathbb N\bar\varepsilon>0 使得

\int_\Omega(|\nabla\phi_{1,n}|^2+|\nabla\phi_{2,n}|^2) \geq(\rho_1+\rho_2)\bar\alpha+\bar\varepsilon.

由于对充分大的 n,

\begin{aligned} \int_{\Omega}\left(\left|\nabla \Psi_{1, n}(0, \cdot)\right|^{2}+\left|\nabla \Psi_{2, n}(0, \cdot)\right|^{2}\right) & =\int_{\Omega}\left(\left|\nabla \psi_{1, n}\right|^{2}+\left|\nabla \psi_{2, n}\right|^{2}\right) \\ & \leq \int_{\Omega}\left(\left|\nabla v_{1, n}\right|^{2}+\left|\nabla v_{2, n}\right|^{2}\right)+\mathrm{o}(1) \\ & \leq\left(\rho_{1}+\rho_{2}\right) \bar{\alpha}+\mathrm{o}(1) \end{aligned}

则存在 \bar t_n\in (0,t_n) 使得 (\Psi_{1,n} (\bar t_n,\cdot), \Psi_{2,n}(\bar t_n,\cdot)) 满足 (4.3) 式且

\int_\Omega (|\nabla \Psi_{1,n}(\bar t_n,\cdot)|^2 + |\nabla \Psi_{2,n} (\bar t_n,\cdot)|^2) = (\rho_1+\rho_2)\bar \alpha + \text{o}(1),

特别地, (4.4) 式. 根据引理 4.2, 存在 (\bar v_1,\bar v_2)\in G_{\bar \alpha} 使得

\int_\Omega (|\nabla \bar v_{1}|^2 + |\nabla\bar v_2|^2) =(\rho_1+\rho_2)\bar \alpha,

这与假设 c_{\bar\alpha}<\tilde c_{\bar \alpha} 矛盾.

假设 c_{\bar \alpha} <\tilde c_{\bar \alpha}, 我们证明引理 4.3. 现在验证由于 c_{\bar \alpha}<\hat c_{\bar \alpha}, 这个假设是满足的.

引理4.4\bar \alpha 如上所述. 则 c_{\bar \alpha}<\tilde c_{\bar\alpha}.

通过矛盾, 假设 \tilde c=c, 则存在 \varepsilon _n\to 0(\tau_{1,n},\tau_{2,n})\in H^1_0(\Omega;\mathbb C^2) 使得

\|(\tau_{1,n},\tau_{2,n})\|^2_{H^1_0(\Omega;\mathbb C^2)}=\bar\alpha(\rho_1+\rho_2),\ \ \int_\Omega|\tau_{i,n}|^2=\rho_i,

且对任意的 n,

c_{\bar \alpha}\leq \mathcal{E}_1(\tau_{1,n},\tau_{2,n})\leq c_{\bar\alpha}+\varepsilon_n, \ \ i=1,2

v_{i,n}:=|\tau _{i,n}|, i=1,2, diamagnetic 不等式意味着

\begin{equation}\label{eq:c_tilde1} \|(v_{1,n},v_{2,n})\|^2_{H^1_0(\Omega;\mathbb R^2)} \leq \|(\tau_{1,n},\tau_{2,n})\|^2_{H^1_0(\Omega;\mathbb C^2)}=(\rho_1+\rho_2)\bar\alpha, \end{equation}
(4.10)

因此 (v_{1,n},v_{2,n})c_{\bar\alpha} 的极小化问题且

\begin{equation}\label{eq:c_tilde2} c_{\bar\alpha}\leq\mathcal{E}_1(v_{1,n},v_{2,n})\leq\mathcal{E}_1(\tau_{1,n},\tau _{2,n}) \leq c_{\bar\alpha}+\varepsilon_n. \end{equation}
(4.11)

特别地,

\begin{equation}\label{eq:c_tilde3} \frac{1}{2}\left(\|(\tau _{1,n},\tau _{2,n})\|^2_{H^1_0(\Omega;\mathbb C^2)}-\|(v_{1,n},v_{2,n})\|^2_{H^1_0(\Omega; \mathbb R^2)}\right)=\mathcal{E}_1 (\tau _{1,n},\tau _{2,n})-\mathcal{E}_1 (v_{1,n},v_{2,n})\leq \varepsilon _n. \end{equation}
(4.12)

由引理 4.2, 在 H^1_0(\tau _{1,n},\tau _{2,n})\to (\tau _{1,\infty},\tau _{2,\infty})(v_{1,n},v_{2,n})\to (v_{1,\infty},v_{2,\infty}). 由 (4.11), (4.12) 式, 结合 Sobolev 嵌入的连续性和引理 2.1 与 2.3, 推导出

\int_\Omega|f(v_{i,n})|^{p+1}=|f(v_{i,\infty})|^{p+1},\ \ i=1,2,
\int_\Omega|f(v_{1,n})|^{\frac{p+1}{2} }|f(v_{2,n})|^{\frac{p+1}{2} } =|f(v_{1,\infty})|^{\frac{p+1}{2}}|f(v_{2,\infty})|^{\frac{p+1}{2}},

\mathcal{E}_1(v_{1,\infty},v_{2,\infty})= c_{\bar \alpha}\ \ \text{且}\ \ \|(v_{1,\infty},v_{2,\infty})\|^2_{H^1_0(\Omega;\mathbb R^2)}= \|(\tau _{1,\infty},\tau _{2,\infty})\|^2_{H^1_0(\Omega;\mathbb C^2)} =\bar\alpha(\rho_1+\rho_2).

(v_{1,\infty},v_{2,\infty})\in\mathcal{U}_{\bar\alpha}, 这与 c_{\bar\alpha}<\hat c_{\bar \alpha} 矛盾.

定理 1.1(b) 和 1.5 的证明 回顾本节第一段, 我们必须证明集合 G_{\bar\alpha} 是 (条件) 轨道稳定的. 证明直接由引理 4.3 和引理4.4 得证.

5 \beta\to-\infty 时的渐近性

本节的目标是完成定理 1.6 的证明. 设 \mu_1,\mu_2>0, 且取 \rho_1,\ \rho_2>0 满足

\left\{\begin{array}{ll} \rho_{1}, \rho_{2}>0, & 1<p<3+\frac{4}{N}, \\ 0<\mu_{1} \rho_{1}^{\frac{2}{N}}, \mu_{2} \rho_{2}^{\frac{2}{N}}<\frac{N+1}{2 N C_{N}}, & p=3+\frac{4}{N}, \\ \max \left\{\mu_{1} \rho_{1}^{2 r}, \mu_{2} \rho_{2}^{2 r}\right\} \cdot\left(\rho_{1}+\rho_{2}\right)^{a-1} \leq \frac{(a-1)^{a-1}}{a^{a}} \lambda_{2}(\Omega)^{-(a-1)}, & 3+\frac{4}{N}<p \leq \frac{3 N+2}{N-2}. \end{array}\right.
(5.1)

注意到所有的条件与 \beta 无关. 通过应用定理 1.1 和 1.4, 引理 3.6 和 (3.18) 式, 推导出对给定的 \beta<0, 存在正函数 v_{1,\beta},v_{2,\beta} 以及 \omega_{1,\beta},\omega_{2,\beta}\in \mathbb R 使得

\left\{\begin{array}{l} -\Delta v_{1, \beta}=f^{\prime}\left(v_{1, \beta}\right)\left(\mu_{1} f^{p}\left(v_{1, \beta}\right)+\beta f^{\frac{p-1}{2}}\left(v_{1, \beta}\right) f^{\frac{p+1}{2}}\left(v_{2, \beta}\right)-\omega_{1, \beta} f\left(v_{1, \beta}\right)\right) \\ -\Delta v_{2, \beta}=f^{\prime}\left(v_{2, \beta}\right)\left(\mu_{2} f^{p}\left(v_{2, \beta}\right)+\beta f^{\frac{p+1}{2}}\left(v_{2, \beta}\right) f^{\frac{p+1}{2}}\left(v_{1, \beta}\right)-\omega_{2, \beta} f\left(v_{2, \beta}\right)\right), \\ \int_{\Omega} f^{2}\left(v_{i}\right)=\rho_{i}, \quad i=1,2 \\ \left(v_{1}, v_{2}\right) \in H_{0}^{1}\left(\Omega ; \mathbb{R}^{2}\right), \end{array}\right.
(5.2)

\begin{equation}\label{eq_uniformalphabar} \mathcal{E}_1(v_{1,\beta},v_{2,\beta})=\inf_\mathcal{M} \mathcal{E}_1,\ \ \ \ \text{当}\ \ 1<p\leq 3+\frac{4}{N}, \end{equation}
(5.3)
\begin{equation}\label{eq_uniformalphabar2} \mathcal{E}_1(v_{1,\beta},v_{2,\beta}) =\inf_{\mathcal{B}_{\bar\alpha}}\mathcal{E}_1,\ \ (v_{1,\beta},v_{2,\beta})\in\mathcal{B}_{\bar\alpha}\setminus \mathcal{U}_{\bar\alpha},\ \ \ \ \text{当}\ \ 3+\frac{4}{N}<p\leq\frac{3N+2}{N-2}, \end{equation}
(5.4)

其中在 (5.4) 式中 \bar{\alpha}:=\frac{a}{a-1}\lambda_2(\Omega).

引理5.1 在以上假设下, 存在一个常数 C>0\beta 无关, 使得

\|v_{i,\beta}\|_{H^1_0(\Omega)} + |v_{i,\beta} |_{L^\infty(\Omega)} +|\omega_{i,\beta}| \leq C\ \ \ \ \text{对每个}\ \ \beta<0, \ \ i=1,2.

(\xi_1,\xi_2)\in \mathcal{M} (如果 p\leq 3+\frac{4}{N})(\xi_1,\xi_2)\in \mathcal{B} _{\bar{\alpha}} (如果 p>3+\frac{4}{N}), 且在其余情形有 \xi_1\cdot\xi_2\equiv 0. 根据(5.3)-(5.4) 式, 有

\begin{matrix}\label{111} \mathcal{E}_1(v_{1,\beta},v_{2,\beta})\leq\mathcal{E}_1(\xi_1,\xi_2)&= \frac{1}{2}\int_\Omega(|\nabla \xi_{1}|^2+|\nabla \xi_{2}|^2)-\frac{1}{p+1}\int_\Omega (\mu_1|f(\xi_1)|^{p+1}+\mu_2 |f(\xi_2)|^{p+1})\\ &=:C_1, \end{matrix}
(5.5)

其中 C_1\beta<0 无关.

根据 (5.4) 式的第一段, 我们得到 \{(v_{1,\beta},v_{2,\beta})\}_{\beta<0}H^1_0(\Omega) 中对 p>3+\frac{4}{N} 是一致有界的. 当 1<p<3+\frac{4}{N}, H^1_0-有界性可根据 (5.5) 式以及估计

\begin{aligned} \mathcal{E}_{1}\left(v_{1, \beta}, v_{2, \beta}\right) \geq & \left|\nabla v_{1, \beta}\right|_{2}^{2 a}\left(\frac{1}{2}\left|\nabla v_{1, \beta}\right|_{2}^{2-2 a}-\frac{C_{N, p}}{p+1} \mu_{1} \rho_{1}^{2 r}\right) \\ & +\left|\nabla v_{2, \beta}\right|_{2}^{2 a}\left(\frac{1}{2}\left|\nabla v_{2, \beta}\right|_{2}^{2-2 a}-\frac{C_{N, p}}{p+1} \mu_{2} \rho_{2}^{2 r}\right) \end{aligned}

这与 (2.9) 式关于 \beta<0 相对应, 而对 p=3+\frac{4}{N}, 根据 (5.5) 式和

\mathcal{E}_1(v_{1,\beta},v_{2,\beta})\geq\frac{1}{2}\left(1- \frac{2NC_N\mu_1\rho_1^{\frac{2}{N}}}{N+1}\right)|\nabla v_{1,\beta}|_2^2 +\frac{1}{2}\left(1-\frac{2NC_N\mu_2\rho_2^{\frac{2}{N}}}{N+1}\right) |\nabla v_{2,\beta}|_2^2,

(见 (2.10) 式和 \beta<0). 根据 Sobolev 嵌入 H^1_0(\Omega)\hookrightarrow L^{p+1}(\Omega), 我们得到 \{(v_{1,\beta},v_{2,\beta})\}_{\beta<0}L^{p+1}-范数下是一致有界的. 特别地

\begin{align*} 0&\leq\frac{2(-\beta)}{p+1}\int_\Omega (f(v_{1,\beta})f(v_{2,\beta}))^{\frac{p+1}{2} }\\ &\leq \mathcal{E}_1(v_{1,\beta},v_{2,\beta})+\frac{1}{p+1} \int_\Omega (\mu_1f^{p+1}(v_{1,\beta}) +\mu_2 f^{p+1}(v_{2,\beta}))\leq C_2. \end{align*}

通过用 f(v_{1,\beta}) 试验 (5.2) 式中的第一个方程以及用 f(v_{2,\beta}) 试验 (5.2) 式中的第二个方程, 应用前边的估计得到 i,j\in \{1,2\}, i\neq j

\begin{aligned} \rho_{i}\left|\omega_{i, \beta}\right| & =\left|\int_{\Omega}\left(\mu_{i} f\left(v_{i, \beta}\right)^{p+1}+\beta\left(f\left(v_{1, \beta}\right) f\left(v_{2, \beta}\right)\right)^{\frac{p+1}{2}}-\left|\nabla v_{i, \beta}\right|^{2}\right)\right| \\ & \leq \int_{\Omega}\left(\mu_{i} f\left(v_{i, \beta}\right)^{p+1}+|\beta|\left(f\left(v_{1, \beta}\right) f\left(v_{2, \beta}\right)\right)^{\frac{p+1}{2}}+\left|\nabla v_{i, \beta}\right|^{2}\right) \leq C_{3}. \end{aligned}

现在, 我们应用 Brézis-Kato-Moser 型理论[p1264-1265], 得到 \{(v_{1,\beta},v_{2,\beta})\}_{\beta<0}L^\infty-一致有界性.

定理 1.6 的证明 应用引理 5.1, \{(v_{1,\beta},v_{2,\beta})\}_{\beta<0} 满足文献[39,定理 1.3 和 1.5]的条件. 因此该序列在 C^{0,\alpha}(\overline \Omega) 中对任意的 0<\alpha<1 是一致有界的, 且存在 v_1,v_2\in C^{0,1}(\overline \Omega), 在 \Omegav_1,v_2\geq 0, 并且 (\omega_1,\omega_2)\in\mathbb R^2 使得当 \beta\to -\infty

v_{i,\beta}\to v_i\ \ \text{在}\ \ C^{0,\alpha}(\overline\Omega)\cap H^1_0(\Omega),\ \ \omega_i\to\omega_i.

根据文献 [15,定理 1.1], 得到

-\Delta(v_1-v_2)\geq \mu_1f^{p}(v_1)f^\prime(v_1)-\mu_2f^p(v_2)f^\prime(v_2) -\omega_1f(v_1)f^\prime(v_1)+\omega_2f(v_2)f^\prime(v_2)\ \ \text{在}\ \ \Omega,

-\Delta(v_2-v_1)\geq \mu_2f^{p}(v_2)f^\prime(v_2)-\mu_1f^p(v_1)f^\prime(v_1) -\omega_2f(v_2)f^\prime(v_2)+\omega_1f(v_1)f^\prime(v_1)\ \ \text{在}\ \ \Omega. 注意到 f^{\frac{p+1}{2}}(v_{1,\beta})f^{\frac{p+1}{2}}(v_{2,\beta})\leq |v_{1,\beta}v_{2,\beta}|^{\frac{p+1}{2}}=0, 根据引理 2.1, 取 w:=v_1-v_2 即证.

参考文献

Adachi S, Watanabe T.

Uniqueness of the ground state solutions for quasilinear Schrödinger equations

Nonlinear Anal, 2012, 75(2): 819-833

[本文引用: 2]

Agueh M.

Sharp Gagliardo-Nirenberg inequalities via p-Laplacian type equations

NoDEA Nonlinear Differential Equations Appl, 2008, 15: 457-472

[本文引用: 2]

Ambrosetti A, Wang Z Q.

Positive solutions to a class of quasilinear elliptic equations on \mathbb R

Discrete Contin Dyn Syst, 2003, 9(1): 55-68

[本文引用: 2]

Badiale M, Serra E. Semilinear Elliptic Equations for Beginners: Existence Results via the Variational Approach. London: Springer-Verlag, 2011

[本文引用: 2]

Bartsch T, Jeanjean L, Soave N.

Normalized solutions for a system of coupled cubic Schrödinger equations on \mathbb R^3

J Math Pures Appl, 2016, 106(4): 583-614

[本文引用: 3]

Bartsch T, Zhong X, Zou W.

Normalized solutions for a coupled Schrödinger system

Math Ann, 2021, 380: 1713-1740

[本文引用: 1]

Brézis H, Lieb E.

A relation between pointwise convergence of functions and convergence of functionals

Proc Amer Math Soc, 1983, 88(3): 486-490

[本文引用: 2]

Brézis H, Nirenberg L.

Positive solutions of nonlinear elliptic equations involving critical Sobolev exponents

Commun Pure Appl Math, 1983, 36(4): 437-77

[本文引用: 1]

Brizhik L, Eremko A, Piette B, Zahrzewski W J.

Electron self-trapping in a discrete two-dimensional lattice

Phys D, 2001, 159(1/2): 71-90

[本文引用: 1]

Brizhik L, Eremko A, Piette B, Zahkrzewski W J.

Static solutions of a D-dimensional modified nonlinear Schrödinger equation

Nonlinearity, 2003, 16(4): 1481-1497

[本文引用: 1]

Chen Z, Zou W.

Normalized solutions for nonlinear Schrödinger systems with linear couples

J Math Anal Appl, 2021, 499(1): 125013

[本文引用: 1]

Chen Z, Zhong X, Zou W.

Normalized solutions for nonlinear Schrödinger systems with special mass-mixed terms: The linear couple case

arXiv preprint arXiv:2107.12564, 2021

[本文引用: 1]

Colin M, Jeanjean L.

Solutions for a quasilinear Schrödinger equation: a dual approach

Nonlinear Anal, 2004, 56(2): 213-226

[本文引用: 4]

Colin M, Jeanjean L, Squassina M.

Stability and instability results for standing waves of quasi-linear Schrödinger equations

Nonlinearity, 2010, 23(6): 1353-1385

[本文引用: 3]

Dancer N, Wang K, Zhang Z.

The limit equation for the Gross-Pitaevskii equations and S Terracini's conjecture

J Funct Anal, 2012, 262(3): 1087-1131

[本文引用: 1]

Gladiali F, Squassina M.

Uniqueness of ground states for a class of quasi-linear elliptic equations

Adv Nonlinear Anal, 2012, 1: 159-179

[本文引用: 2]

Gou T, Jeanjean L.

Multiple positive normalized solutions for nonlinear Schrödinger systems

Nonlinearity, 2018, 31(5): 2319-2345

Jeanjean L.

On the existence of bounded Palais-Smale sequences and application to a Landesman-Lazer-type problem set on \mathbb R^N

Proc Roy Soc Edinburgh Sect A, 1999, 129(4): 787-809

[本文引用: 2]

Jeanjean L, Luo T.

Sharp non-existence results of prescribed L^2-norm solutions for some class of Schrödinger-Poisson and quasi-linear equations

Z Angew Math Phys, 2013, 64(4): 937-954

[本文引用: 3]

Jeanjean L, Luo T, Wang Z Q.

Multiple normalized solutions for quasi-linear Schrödinger equations

J Differential Equations, 2015, 259(8): 3894-3928

[本文引用: 7]

Kurihara S.

Large-amplitude quasi-solitons in superfluid films

J Phys Soc Japan, 1981, 50(10): 3262-3267

[本文引用: 1]

Li H, Zou W.

Normalized ground states for semilinear elliptic systems with critical and subcritical nonlinearities

J Fixed Point Theory Appl, 2021, 23: 1-30

Li H, Zou W.

Quasilinear Schrödinger equations: ground state and infinitely many normalized solutions

Pacific J Math, 2023, 322(1): 99-138

[本文引用: 2]

Lieb E, Loss M.

Analysis, Graduate Studies in Mathematics

Rhode Island: American Mathematical Society, 2001

[本文引用: 1]

Liu J, Wang Y, Wang Z Q.

Solutions for quasilinear Schrödinger equations via the Nehari method

Comm Partial Differential Equations, 2004, 29(5/6): 879-901

[本文引用: 1]

Moroz V, Schaftingen J.

Groundstates of nonlinear Choquard equations: existence, qualitative properties and decay asymptotics

J Funct Anal, 2013, 265(2): 153-184

[本文引用: 2]

Noris B, Tavares H, Verzini G.

Existence and orbital stability of the ground states with prescribed mass for the L^2-critical and supercritical NLS on bounded domains

Anal PDE, 2015, 7(8): 1807-1838

[本文引用: 4]

Noris B, Tavares H, Verzini G.

Normalized solutions for nonlinear Schrödinger systems on bounded domains

Nonlinearity, 2019, 32(3): Article 1044

[本文引用: 12]

Noris B, Tavares H, Verzini G.

Stable solitary waves with prescribed L^2-mass for the cubic Schrödinger system with trapping potentials

Discrete Contin Dyn Syst, 2015, 35: 6085-6112

[本文引用: 2]

Pellacci B, Pistoia A, Vaira G, Verzini G.

Normalized concentrating solutions to nonlinear elliptic problems

J Differential Equations, 2021, 275: 882-919

Perez-Garcia V, Michinel H, Cirac J, et al.

Dynamics of Bose-Einstein condensates: Variational solutions of the Gross-Pitaevskii equations

Phys Rev A, 1997, 56(2): Article 1424

Pierotti D, Verzini G.

Normalized bound states for the nonlinear Schrödinger equation in bounded domains

Calc Var Partial Differential Equations, 2017, 56: 1-27

[本文引用: 2]

Poppenberg M, Schmitt K, Wang Z Q.

On the existence of soliton solutions to quasilinear Schrödinger equations

Calc Var Partial Differential Equations, 2002, 14(3): 329-344

[本文引用: 2]

Selvitella A.

Uniqueness and nondegeneracy of the ground state for a quasilinear Schrödinger equation with a small parameter

Nonlinear Anal, 2011, 74(5): 1731-1737

[本文引用: 2]

Serrin J, Tang M.

Uniqueness of ground states for quasilinear elliptic equations

Indiana Univ Math J, 2000, 49: 897-923

[本文引用: 1]

Severo U.

Symmetric and nonsymmetric solutions for a class of quasilinear Schrödinger equations

Adv Nonlinear Stud, 2008, 8(2): 375-389

[本文引用: 1]

Soave N.

Normalized ground states for the NLS equation with combined nonlinearities

J Differential Equations, 2020, 269(9): 6941-6987

[本文引用: 1]

Soave N.

Normalized ground states for the NLS equation with combined nonlinearities: The Sobolev critical case

J Funct Anal, 2020, 279(6): 108610

[本文引用: 1]

Soave N, Tavares H, Terracini S, Zilio A.

Hölder bounds and regularity of emerging free boundaries for strongly competing Schrödinger equations with nontrivial grouping

Nonlinear Anal, 2016, 138: 388-427

[本文引用: 1]

/