コンテンツにスキップ

作用素不等式と積分作用素:Löwner–Heinz・Heinz–Kato・Schur テスト・Schatten クラス・Daleckii–Krein・トレース不等式

対象読者: 本テーマのコンパクト作用素・トレースクラスの章とスペクトル定理の章を学び、複素関数論の最大値原理を知っている人。補間理論は仮定しません。

前提となる記事:コンパクト作用素・Hilbert–Schmidt 作用素・トレースクラス:特異値で測る作用素の大きさ、スペクトル定理と関数計算:非有界自己共役作用素を射影値測度で対角化する

生 Markdown

数のあいだには大小関係があり、0≤a≤b0\le a\le b なら a2≤b2a^2\le b^2 のように、関数で写しても大小は保たれます。行列(作用素)にも「AA は BB 以下」という大小の付け方がありますが、関数で写すと大小が保たれないことがあります。この記事では、どんな関数なら大小が保たれるのか、また作用素の大きさを不等式でどう測るのかを、証明つきで順に調べます。読み終えると、行列の平方根や対数の単調性、積分作用素の有界性の判定、トレースの不等式を、次元によらない形で使えるようになります。

  • 自己共役作用素の順序 A≤BA\le B(Löwner 順序。すべてのベクトル xx で ⟨x,Ax⟩≤⟨x,Bx⟩\langle x,Ax\rangle\le\langle x,Bx\rangle という大小関係)は、関数で写すと壊れやすい順序です。A≤BA\le B から A2≤B2A^2\le B^2 は出ませんが、0≤α≤10\le\alpha\le1 なら Aα≤BαA^\alpha\le B^\alpha です(Löwner–Heinz の不等式)。証明は tαt^\alpha を「t/(s+t)t/(s+t) の重ね合わせ」と見る積分表示で、log⁡\log も同じく単調です。
  • 複素補間(三線定理。帯状の領域の両端での大きさから、内側での大きさを抑える定理)を使うと、Heinz–Kato の不等式 ∣⟨y,Tx⟩∣≤∥Aαx∥∥B1−αy∥|\langle y,Tx\rangle|\le\|A^\alpha x\|\|B^{1-\alpha}y\| と Cordes の不等式 ∥AαBα∥≤∥AB∥α\|A^\alpha B^\alpha\|\le\|AB\|^\alpha が得られます。
  • 積分作用素の L2L^2 有界性(核 K(x,y)K(x,y) の積分作用素のノルムが有限であること)は、正の試験関数を 1 組見つければ判定できます(Schur テスト)。Hilbert の不等式(定数 π\pi)と Hardy の不等式(定数 22)がすぐに出ます。
  • Schatten ノルム(特異値の pp 乗和の 1/p1/p 乗)には正規直交系による変分表示があり、そこから三角不等式・完備性・Hölder 不等式 ∥KL∥1≤∥K∥p∥L∥q\|KL\|_1\le\|K\|_p\|L\|_q が証明できます。f(x)g(−i∇)f(x)g(-i\nabla) の Hilbert–Schmidt ノルムは (2π)−d/2∥f∥2∥g∥2(2\pi)^{-d/2}\|f\|_2\|g\|_2 で、三線定理で p=∞p=\infty とつなぐと Kato–Seiler–Simon の不等式 ∥f(x)g(−i∇)∥p≤(2π)−d/p∥f∥p∥g∥p\|f(x)g(-i\nabla)\|_p\le(2\pi)^{-d/p}\|f\|_p\|g\|_p(p≥2p\ge2)が、2 進分解で gg を弱 LpL^p に広げると Cwikel の評価 sn≤Cn−1/p∥f∥p[g]p,∞s_n\le Cn^{-1/p}\|f\|_p[g]_{p,\infty}(p>2p\gt2)が出ます。
  • f(A)−f(B)f(A)-f(B) はスペクトル射影で挟んだ「二重作用素積分」∑f[1](λi,μj)Pi(A−B)Qj\sum f^{[1]}(\lambda_i,\mu_j)P_i(A-B)Q_j で書け、Hilbert–Schmidt ノルムでは ∥f(A)−f(B)∥2≤∥f′∥∞∥A−B∥2\|f(A)-f(B)\|_2\le\|f'\|_\infty\|A-B\|_2 が次元によらず成り立ちます。微分版が Daleckii–Krein 公式です。Klein の不等式から Gibbs の変分原理と Peierls–Bogoliubov の不等式を導き、Weyl の優越不等式と Lie–Trotter の積公式から Golden–Thompson の不等式 tr eA+B≤tr(eAeB)\mathrm{tr}\,e^{A+B}\le\mathrm{tr}(e^Ae^B) を証明します。
  • 並列和 (P−1+Q−1)−1(P^{-1}+Q^{-1})^{-1} の変分表示から重み付き幾何平均の同時凹性を示し、それを左右の掛け算作用素に使って Lieb の凹性定理(Ando の証明)と相対エントロピーの同時凸性を証明します。

実数では 0≤a≤b0\le a\le b なら a2≤b2a^2\le b^2 ですが、行列ではそうなりません。順序 A≤BA\le B は「すべてのベクトルで ⟨x,Ax⟩≤⟨x,Bx⟩\langle x,Ax\rangle\le\langle x,Bx\rangle」という条件で、AA と BB が可換でなければ、固有ベクトルの向きのずれが単調性を壊すからです。にもかかわらず、平方根や対数は順序を保ちます。この非対称な現象を正確に述べるのが作用素単調関数の理論で、量子情報の相対エントロピーの単調性や、Powers–Størmer の不等式(状態の距離の比較)の基礎になっています。

もう一つの主題は、作用素を「数のように」評価する技術です。非有界作用素の定義域の比較(∥Tx∥≤∥Ax∥\|Tx\|\le\|Ax\| から AαA^\alpha についての評価を引き出す)、積分作用素のノルム、Schatten クラスの積、関数 ff を作用素に入れたときのリプシッツ性、統計力学の分配関数の評価がそれにあたります。どれも、有限次元の不等式を次元によらない形で証明することが要点です。

H\mathcal{H} は複素 Hilbert 空間です。自己共役な A,B∈B(H)A,B\in\mathcal{B}(\mathcal{H}) について、全 xx で ⟨x,Ax⟩≤⟨x,Bx⟩\langle x,Ax\rangle\le\langle x,Bx\rangle のとき A≤BA\le B と書きます(Löwner 順序)。A≥0A\ge0 を正、A≥ε>0A\ge\varepsilon>0 となるとき正定値と呼びます。有界自己共役作用素の関数 f(A)f(A) はスペクトル定理による関数計算(定理 6.2[スペクトル定理と関数計算])で定め、σ(A)\sigma(A) 上で連続な ff については ∥f(A)∥=max⁡σ(A)∣f∣\|f(A)\|=\max_{\sigma(A)}|f| です。

補題 2.1(Löwner 順序の基本性質)✓

自己共役な有界作用素について、

  1. A≤BA\le B なら任意の C∈B(H)C\in\mathcal{B}(\mathcal{H}) で C∗AC≤C∗BCC^*AC\le C^*BC。
  2. 0≤A≤B0\le A\le B なら ∥A∥≤∥B∥\|A\|\le\|B\|。
  3. A,BA,B が正定値なら、A≤BA\le B と ∥A1/2B−1/2∥≤1\|A^{1/2}B^{-1/2}\|\le1 は同値。
  4. A,BA,B が正定値で A≤BA\le B なら B−1≤A−1B^{-1}\le A^{-1}。
証明(補題 2.1)

1: ⟨x,C∗ACx⟩=⟨Cx,ACx⟩≤⟨Cx,BCx⟩\langle x,C^*ACx\rangle=\langle Cx,ACx\rangle\le\langle Cx,BCx\rangle。2: 正作用素では ∥A∥=sup⁡∥x∥=1⟨x,Ax⟩\|A\|=\sup_{\|x\|=1}\langle x,Ax\rangle(定理 6.10[有界作用素と共役作用素])。3: x=B−1/2yx=B^{-1/2}y と置けば、⟨x,Ax⟩≤⟨x,Bx⟩\langle x,Ax\rangle\le\langle x,Bx\rangle は ∥A1/2B−1/2y∥2≤∥y∥2\|A^{1/2}B^{-1/2}y\|^2\le\|y\|^2 と同じです。4: 3 から ∥A1/2B−1/2∥≤1\|A^{1/2}B^{-1/2}\|\le1 で、共役を取って ∥B−1/2A1/2∥≤1\|B^{-1/2}A^{1/2}\|\le1。これは 3 を正定値作用素 B−1B^{-1}, A−1A^{-1} に使った形((B−1)1/2((A−1)−1)1/2(B^{-1})^{1/2}((A^{-1})^{-1})^{1/2})なので、B−1≤A−1B^{-1}\le A^{-1} です。

3. 作用素単調関数と Löwner–Heinz の不等式

Section titled “3. 作用素単調関数と Löwner–Heinz の不等式”

定義 3.1(作用素単調関数)↗

区間 II 上の連続関数 ff が作用素単調であるとは、任意の Hilbert 空間と、スペクトルが II に含まれる自己共役な有界作用素 A≤BA\le B について f(A)≤f(B)f(A)\le f(B) となることをいう。

次の補題は、tαt^\alpha を t/(s+t)t/(s+t) という簡単な関数を ss について足し合わせたものとして書く公式です。t/(s+t)t/(s+t) の単調性なら ss を固定して確かめられるので、これを使って A≤BA\le B から Aα≤BαA^\alpha\le B^\alpha を引き出します(定理 3.3 の証明)。

補題 3.2(tαt^α の積分表示)✓

0<α<10<\alpha<1, t≥0t\ge0 について

tα=sin⁡παπ∫0∞sα−1ts+t ds.t^\alpha=\frac{\sin\pi\alpha}{\pi}\int_0^\infty s^{\alpha-1}\frac{t}{s+t}\,ds .
証明(補題 3.2)

t=0t=0 では両辺 00 です。t>0t>0 なら s=tus=tu と置くと右辺の積分は tα∫0∞uα−11+udut^\alpha\int_0^\infty\frac{u^{\alpha-1}}{1+u}du です。11+u=∫0∞e−(1+u)rdr\frac1{1+u}=\int_0^\infty e^{-(1+u)r}dr と Tonelli の定理から

∫0∞uα−11+udu=∫0∞e−r(∫0∞uα−1e−urdu)dr=Γ(α)∫0∞e−rr−αdr=Γ(α)Γ(1−α)=πsin⁡πα\int_0^\infty\frac{u^{\alpha-1}}{1+u}du=\int_0^\infty e^{-r}\Bigl(\int_0^\infty u^{\alpha-1}e^{-ur}du\Bigr)dr=\Gamma(\alpha)\int_0^\infty e^{-r}r^{-\alpha}dr=\Gamma(\alpha)\Gamma(1-\alpha)=\frac{\pi}{\sin\pi\alpha}

で、最後は相反公式(定理 6.5[整関数とガンマ関数])です。

定理 3.3(Löwner–Heinz の不等式)✓

0≤α≤10\le\alpha\le1 なら t↦tαt\mapsto t^\alpha は [0,∞)[0,\infty) 上の作用素単調関数である。すなわち 0≤A≤B0\le A\le B なら Aα≤BαA^\alpha\le B^\alpha。

証明(定理 3.3)

α=0,1\alpha=0,1 は自明なので 0<α<10<\alpha<1 とします。s>0s>0 について s+A≤s+Bs+A\le s+B はどちらも正定値なので、補題 2.1 の 4 から (s+B)−1≤(s+A)−1(s+B)^{-1}\le(s+A)^{-1}、したがって

A(s+A)−1=1−s(s+A)−1≤1−s(s+B)−1=B(s+B)−1.A(s+A)^{-1}=1-s(s+A)^{-1}\le1-s(s+B)^{-1}=B(s+B)^{-1}.

x∈Hx\in\mathcal{H} と AA のスペクトル測度 μx\mu_x について、補題 3.2 と Tonelli の定理(被積分関数は非負)から

⟨x,Aαx⟩=∫λαdμx(λ)=sin⁡παπ∫0∞sα−1⟨x,A(s+A)−1x⟩ ds\langle x,A^\alpha x\rangle=\int\lambda^\alpha d\mu_x(\lambda)=\frac{\sin\pi\alpha}{\pi}\int_0^\infty s^{\alpha-1}\langle x,A(s+A)^{-1}x\rangle\,ds

で、BB についても同じです。被積分関数を比べれば ⟨x,Aαx⟩≤⟨x,Bαx⟩\langle x,A^\alpha x\rangle\le\langle x,B^\alpha x\rangle です。

系 3.4(対数は作用素単調)✓

A,BA,B が正定値で A≤BA\le B なら log⁡A≤log⁡B\log A\le\log B。

証明(系 3.4)

σ(A)∪σ(B)⊂[m,M]⊂(0,∞)\sigma(A)\cup\sigma(B)\subset[m,M]\subset(0,\infty) とすると、(tα−1)/α→log⁡t(t^\alpha-1)/\alpha\to\log t(α↓0\alpha\downarrow0)は [m,M][m,M] 上一様なので、関数計算のノルムの等式から (Aα−1)/α→log⁡A(A^\alpha-1)/\alpha\to\log A がノルムで成り立ちます(BB も同様)。定理 3.3 から (Aα−1)/α≤(Bα−1)/α(A^\alpha-1)/\alpha\le(B^\alpha-1)/\alpha で、極限に移ります。

例 3.5(t² は作用素単調でない)

A=(1111)A=\begin{pmatrix}1&1\\1&1\end{pmatrix}, B=(2111)B=\begin{pmatrix}2&1\\1&1\end{pmatrix} とすると B−A=(1000)≥0B-A=\begin{pmatrix}1&0\\0&0\end{pmatrix}\ge0 で 0≤A≤B0\le A\le B ですが、

B2−A2=(5332)−(2222)=(3110)B^2-A^2=\begin{pmatrix}5&3\\3&2\end{pmatrix}-\begin{pmatrix}2&2\\2&2\end{pmatrix}=\begin{pmatrix}3&1\\1&0\end{pmatrix}

は行列式が −1<0-1<0 なので正ではありません。AA と BB が可換でないことが原因で、可換なら A≤BA\le B から f(A)≤f(B)f(A)\le f(B) は任意の単調関数で成り立ちます。

注意 3.6

(0,∞)(0,\infty) 上の関数が作用素単調であるための必要十分条件は、上半平面に正則に延びて虚部が非負になる(Pick 関数である)ことです(Löwner の定理)。tαt^\alpha が作用素単調なのは 0≤α≤10\le\alpha\le1 に限られます。この記事ではこの特徴づけは証明しません。補題 3.2 の積分表示は、その「t/(s+t)t/(s+t) の正の重ね合わせ」という構造の具体例です。

定理 4.1(Heinz–Kato の不等式)✓

A,B≥0A,B\ge0 を有界作用素、T∈B(H)T\in\mathcal{B}(\mathcal{H}) が全 x,yx,y で ∥Tx∥≤∥Ax∥\|Tx\|\le\|Ax\|, ∥T∗y∥≤∥By∥\|T^*y\|\le\|By\| を満たすとする。このとき 0≤α≤10\le\alpha\le1 について

∣⟨y,Tx⟩∣≤∥Aαx∥ ∥B1−αy∥(x,y∈H).|\langle y,Tx\rangle|\le\|A^\alpha x\|\,\|B^{1-\alpha}y\|\qquad(x,y\in\mathcal{H}).

証明の方針: α=0\alpha=0 と α=1\alpha=1 のときは不等式が仮定そのものです。その間の α\alpha を扱うために、α\alpha を複素数 zz に置き換えた関数 F(z)F(z) を作り、帯の両端での評価から内側の評価を三線定理で引き出します。

証明(定理 4.1)

正定値の場合. A,B≥ε>0A,B\ge\varepsilon>0 とします。log⁡A\log A, log⁡B\log B は有界自己共役で、Az=ezlog⁡A=∑nzn(log⁡A)n/n!A^z=e^{z\log A}=\sum_nz^n(\log A)^n/n! はノルムで正則な整関数、∥Ax+iy∥=max⁡σ(A)λx\|A^{x+iy}\|=\max_{\sigma(A)}\lambda^x は帯 0≤x≤10\le x\le1 で有界、AiyA^{iy} はユニタリです。u,v∈Hu,v\in\mathcal{H} について F(z)=⟨u,Bz−1TA−zv⟩F(z)=\langle u,B^{z-1}TA^{-z}v\rangle は帯の上で有界連続、内部で正則です。∥T∗y∥≤∥By∥\|T^*y\|\le\|By\| に y=B−1wy=B^{-1}w を入れると ∥T∗B−1∥≤1\|T^*B^{-1}\|\le1、同様に ∥TA−1∥≤1\|TA^{-1}\|\le1 なので

∣F(iy)∣≤∥B−1T∥∥u∥∥v∥=∥T∗B−1∥∥u∥∥v∥≤∥u∥∥v∥,∣F(1+iy)∣≤∥TA−1∥∥u∥∥v∥≤∥u∥∥v∥.|F(iy)|\le\|B^{-1}T\|\|u\|\|v\|=\|T^*B^{-1}\|\|u\|\|v\|\le\|u\|\|v\|,\qquad|F(1+iy)|\le\|TA^{-1}\|\|u\|\|v\|\le\|u\|\|v\|.

Hadamard の三線定理(定理 5.4[最大値原理の一般化と正則関数族])から ∣F(α)∣≤∥u∥∥v∥|F(\alpha)|\le\|u\|\|v\|、すなわち ∥Bα−1TA−α∥≤1\|B^{\alpha-1}TA^{-\alpha}\|\le1 です。よって

∣⟨y,Tx⟩∣=∣⟨B1−αy, Bα−1TA−α Aαx⟩∣≤∥B1−αy∥∥Aαx∥.|\langle y,Tx\rangle|=|\langle B^{1-\alpha}y,\,B^{\alpha-1}TA^{-\alpha}\,A^\alpha x\rangle|\le\|B^{1-\alpha}y\|\|A^\alpha x\|.

一般の場合. Aε=A+εA_\varepsilon=A+\varepsilon は ∥(A+ε)x∥2=∥Ax∥2+2ε⟨x,Ax⟩+ε2∥x∥2≥∥Ax∥2\|(A+\varepsilon)x\|^2=\|Ax\|^2+2\varepsilon\langle x,Ax\rangle+\varepsilon^2\|x\|^2\ge\|Ax\|^2 から仮定を満たし、BεB_\varepsilon も同様です。(t+ε)α→tα(t+\varepsilon)^\alpha\to t^\alpha は [0,∥A∥][0,\|A\|] 上一様なので Aεαx→AαxA_\varepsilon^\alpha x\to A^\alpha x で、極限に移れます。

∥Tx∥≤∥Ax∥\|Tx\|\le\|Ax\| は T∗T≤A2T^*T\le A^2 と同じです。T=1T=1, B=AB=A と取れば仮定は自明に成り立ち、結論は ∣⟨y,x⟩∣≤∥Aαx∥∥A1−αy∥|\langle y,x\rangle|\le\|A^\alpha x\|\|A^{1-\alpha}y\|(Cauchy–Schwarz の補間)です。非有界作用素への応用では、∥Tx∥≤∥Ax∥\|Tx\|\le\|Ax\| という定義域の比較から AαA^\alpha に関する評価を引き出すのに使います。

定理 4.2(Cordes の不等式)✓

A,B≥0A,B\ge0 を有界作用素、0≤α≤10\le\alpha\le1 とすると ∥AαBα∥≤∥AB∥α\|A^\alpha B^\alpha\|\le\|AB\|^\alpha。

証明(定理 4.2)

まず BB を正定値とします。c=∥AB∥c=\|AB\| と置くと ∥ABx∥≤c∥x∥\|ABx\|\le c\|x\| なので BA2B≤c2BA^2B\le c^2、補題 2.1 の 1 を C=B−1C=B^{-1} で使って A2≤c2B−2A^2\le c^2B^{-2}。定理 3.3 を指数 α\alpha で使うと A2α≤c2αB−2αA^{2\alpha}\le c^{2\alpha}B^{-2\alpha}、再び 1 を C=BαC=B^\alpha で使って BαA2αBα≤c2αB^\alpha A^{2\alpha}B^\alpha\le c^{2\alpha}。C*-恒等式(定理 3.4[有界作用素と共役作用素])から ∥AαBα∥2=∥BαA2αBα∥≤c2α\|A^\alpha B^\alpha\|^2=\|B^\alpha A^{2\alpha}B^\alpha\|\le c^{2\alpha} です。一般の BB は B+εB+\varepsilon で近似し、ε→0\varepsilon\to0 とします。

定理 5.1(Schur テスト)✓

(X,μ)(X,\mu), (Y,ν)(Y,\nu) を σ\sigma-有限な測度空間、K:X×Y→CK:X\times Y\to\mathbb{C} を可測関数とする。正の可測関数 pp on XX, qq on YY と定数 C1,C2C_1,C_2 があって、

∫Y∣K(x,y)∣q(y) dν(y)≤C1p(x)  (μ-a.e. x),∫X∣K(x,y)∣p(x) dμ(x)≤C2q(y)  (ν-a.e. y)\int_Y|K(x,y)|q(y)\,d\nu(y)\le C_1p(x)\ \ (\mu\text{-a.e. }x),\qquad\int_X|K(x,y)|p(x)\,d\mu(x)\le C_2q(y)\ \ (\nu\text{-a.e. }y)

ならば、(Kf)(x)=∫YK(x,y)f(y) dν(y)(Kf)(x)=\int_YK(x,y)f(y)\,d\nu(y) は f∈L2(Y)f\in L^2(Y) について a.e. xx で絶対収束し、∥Kf∥L2(X)≤C1C2 ∥f∥L2(Y)\|Kf\|_{L^2(X)}\le\sqrt{C_1C_2}\,\|f\|_{L^2(Y)} を満たす。

証明(定理 5.1)

∣K∣∣f∣=(∣K∣q)1/2⋅(∣K∣/q)1/2∣f∣|K||f|=(|K|q)^{1/2}\cdot(|K|/q)^{1/2}|f| と分けて Cauchy–Schwarz の不等式を使うと

∫Y∣K(x,y)∣∣f(y)∣ dν≤(C1p(x))1/2(∫Y∣K(x,y)∣∣f(y)∣2q(y)dν)1/2.\int_Y|K(x,y)||f(y)|\,d\nu\le\Bigl(C_1p(x)\Bigr)^{1/2}\Bigl(\int_Y|K(x,y)|\frac{|f(y)|^2}{q(y)}d\nu\Bigr)^{1/2}.

2 乗して xx で積分し、Tonelli の定理(定理 3.4[直積測度と Fubini–Tonelli の定理])で積分の順序を入れ替えると

∫X(∫Y∣K∣∣f∣ dν)2dμ≤C1∫Y∣f(y)∣2q(y)(∫X∣K(x,y)∣p(x) dμ)dν≤C1C2∥f∥22.\int_X\Bigl(\int_Y|K||f|\,d\nu\Bigr)^2d\mu\le C_1\int_Y\frac{|f(y)|^2}{q(y)}\Bigl(\int_X|K(x,y)|p(x)\,d\mu\Bigr)d\nu\le C_1C_2\|f\|_2^2 .

とくに左辺の内側の積分は a.e. xx で有限です。

p=q=1p=q=1 と取れば、「行和と列和が有界なら L2L^2 有界」という Young 型の判定になります。重み p,qp,q を選べるのが Schur テストの強みで、斉次な核では冪関数がうまくはまります。

例 5.2(Hilbert の不等式)

L2(0,∞)L^2(0,\infty) 上の K(x,y)=1/(x+y)K(x,y)=1/(x+y) に p(x)=q(x)=x−1/2p(x)=q(x)=x^{-1/2} を使うと、y=xuy=xu と置いて

∫0∞y−1/2x+ydy=x−1/2∫0∞u−1/21+udu=πx−1/2\int_0^\infty\frac{y^{-1/2}}{x+y}dy=x^{-1/2}\int_0^\infty\frac{u^{-1/2}}{1+u}du=\pi x^{-1/2}

(補題 3.2 の証明の積分で α=12\alpha=\frac12)、対称性から C1=C2=πC_1=C_2=\pi です。よって

∣∫0∞ ⁣ ⁣∫0∞f(x)g(y)‾x+ydx dy∣≤π∥f∥2∥g∥2.\Bigl|\int_0^\infty\!\!\int_0^\infty\frac{f(x)\overline{g(y)}}{x+y}dx\,dy\Bigr|\le\pi\|f\|_2\|g\|_2 .

定数 π\pi は最良です(Mellin 変換で対角化すると、この作用素のスペクトルが [0,π][0,\pi] であることがわかります)。

6. Schatten クラスの変分表示と Hölder 不等式

Section titled “6. Schatten クラスの変分表示と Hölder 不等式”

特異値 sns_n と Schatten クラス Sp\mathcal{S}_p は 定義 7.3[コンパクト作用素・Hilbert–Schmidt 作用素・トレースクラス] のとおりです。この節では H\mathcal{H} を無限次元の可分 Hilbert 空間とします(有限次元では、以下で使うミニマックス原理を Courant–Fischer の定理 定理 3.1[行列の固有値評価と不等式] に置き換えれば同じです)。

補題 6.1(特異値の比較)✓

K∈K(H)K\in\mathcal{K}(\mathcal{H}), A,B∈B(H)A,B\in\mathcal{B}(\mathcal{H}) について sn(K∗)=sn(K)s_n(K^*)=s_n(K) かつ sn(AKB)≤∥A∥sn(K)∥B∥s_n(AKB)\le\|A\|s_n(K)\|B\|。

証明(補題 6.1)

コンパクトな正作用素 CC の nn 番目に大きい固有値(足りなければ 00)を λn(C)\lambda_n(C) とすると、ミニマックス原理(定理 6.3[作用素の収束と摂動])を −C-C に使って λn(C)=−μn(−C)\lambda_n(C)=-\mu_n(-C) です(−C-C の本質スペクトルは {0}\{0\})。系 6.5[作用素の収束と摂動] から 0≤C1≤C20\le C_1\le C_2 なら λn(C1)≤λn(C2)\lambda_n(C_1)\le\lambda_n(C_2)。sn(K)2=λn(K∗K)s_n(K)^2=\lambda_n(K^*K) で、K∗A∗AK≤∥A∥2K∗KK^*A^*AK\le\|A\|^2K^*K なので sn(AK)≤∥A∥sn(K)s_n(AK)\le\|A\|s_n(K)。極分解 K=V∣K∣K=V|K| から KK∗=V∣K∣2V∗KK^*=V|K|^2V^* で、VV は Ran∣K∣‾\overline{\mathrm{Ran}|K|} から RanK‾\overline{\mathrm{Ran}K} へのユニタリなので、KK∗KK^* と K∗KK^*K の正の固有値は重複度も込めて一致し、sn(K∗)=sn(K)s_n(K^*)=s_n(K)。最後に sn(KB)=sn(B∗K∗)≤∥B∥sn(K)s_n(KB)=s_n(B^*K^*)\le\|B\|s_n(K) です。

定理 6.2(Schatten ノルムの変分表示)✓

1≤p<∞1\le p<\infty, K∈K(H)K\in\mathcal{K}(\mathcal{H}) とする。

∑nsn(K)p=sup⁡{∑n∣⟨gn,Khn⟩∣p:(gn),(hn) は正規直交系}\sum_ns_n(K)^p=\sup\Bigl\{\sum_n|\langle g_n,Kh_n\rangle|^p:(g_n),(h_n)\text{ は正規直交系}\Bigr\}

(両辺とも ∞\infty でありうる)。したがって ∥⋅∥p\|\cdot\|_p は Sp\mathcal{S}_p 上のノルムで、Sp\mathcal{S}_p は Banach 空間であり、∥K∗∥p=∥K∥p\|K^*\|_p=\|K\|_p, ∥AKB∥p≤∥A∥∥K∥p∥B∥\|AKB\|_p\le\|A\|\|K\|_p\|B\| が成り立つ。

証明の方針: 右辺の和を、特異値分解の係数 anka_{nk} で行・列の両方向から押さえます。Bessel の不等式で行和・列和が 11 以下になることが鍵です。

証明(定理 6.2)

特異値分解 K=∑ksk∣fk⟩⟨ek∣K=\sum_ks_k|f_k\rangle\langle e_k|(定理 7.2[コンパクト作用素・Hilbert–Schmidt 作用素・トレースクラス])を使い、ank=∣⟨gn,fk⟩∣∣⟨ek,hn⟩∣a_{nk}=|\langle g_n,f_k\rangle||\langle e_k,h_n\rangle| と置きます。Cauchy–Schwarz と Bessel の不等式から ∑nank≤1\sum_na_{nk}\le1, ∑kank≤1\sum_ka_{nk}\le1 です。∣⟨gn,Khn⟩∣≤∑kskank|\langle g_n,Kh_n\rangle|\le\sum_ks_ka_{nk} で、重み anka_{nk} についての Hölder の不等式から

(∑kskank)p≤(∑kank)p−1∑kankskp≤∑kankskp.\Bigl(\sum_ks_ka_{nk}\Bigr)^p\le\Bigl(\sum_ka_{nk}\Bigr)^{p-1}\sum_ka_{nk}s_k^p\le\sum_ka_{nk}s_k^p .

nn について足すと ∑kskp∑nank≤∑kskp\sum_ks_k^p\sum_na_{nk}\le\sum_ks_k^p。gn=fng_n=f_n, hn=enh_n=e_n で等号です。

右辺の形から、ℓp\ell^p の Minkowski の不等式(定理 4.5[L^p 空間と関数解析への導入])で三角不等式が、⟨g,K∗h⟩=⟨h,Kg⟩‾\langle g,K^*h\rangle=\overline{\langle h,Kg\rangle} から ∥K∗∥p=∥K∥p\|K^*\|_p=\|K\|_p が出ます。∥AKB∥p\|AKB\|_p の評価は 補題 6.1 です。完備性: (Km)(K_m) が ∥⋅∥p\|\cdot\|_p のコーシー列なら、∥⋅∥≤∥⋅∥p\|\cdot\|\le\|\cdot\|_p から作用素ノルムで KK に収束し、KK はコンパクトです。有限個の正規直交ベクトルについて

(∑n≤N∣⟨gn,(K−Km)hn⟩∣p)1/p=lim⁡l→∞(∑n≤N∣⟨gn,(Kl−Km)hn⟩∣p)1/p≤sup⁡l≥m∥Kl−Km∥p\Bigl(\sum_{n\le N}|\langle g_n,(K-K_m)h_n\rangle|^p\Bigr)^{1/p}=\lim_{l\to\infty}\Bigl(\sum_{n\le N}|\langle g_n,(K_l-K_m)h_n\rangle|^p\Bigr)^{1/p}\le\sup_{l\ge m}\|K_l-K_m\|_p

で、上限を取れば ∥K−Km∥p≤sup⁡l≥m∥Kl−Km∥p→0\|K-K_m\|_p\le\sup_{l\ge m}\|K_l-K_m\|_p\to0 です。

定理 6.3(Schatten クラスの Hölder 不等式)✓

1≤p,q≤∞1\le p,q\le\infty, 1p+1q=1\frac1p+\frac1q=1 とする(S∞\mathcal{S}_\infty はコンパクト作用素に作用素ノルムを入れたもの)。K∈SpK\in\mathcal{S}_p, L∈SqL\in\mathcal{S}_q なら KL∈S1KL\in\mathcal{S}_1 で ∥KL∥1≤∥K∥p∥L∥q\|KL\|_1\le\|K\|_p\|L\|_q。

証明(定理 6.3)

p=∞p=\infty は 定理 6.2 の ∥KL∥1≤∥K∥∥L∥1\|KL\|_1\le\|K\|\|L\|_1 です。∥KL∥1=∥L∗K∗∥1\|KL\|_1=\|L^*K^*\|_1 で、p,qp,q の一方は 22 以上なので、2≤p<∞2\le p<\infty としてよいです。

補題. p≥2p\ge2 と正規直交系 (vk)(v_k) について ∑k∥Kvk∥p≤∥K∥pp\sum_k\|Kv_k\|^p\le\|K\|_p^p。実際、∣K∣2=∑jsj2∣ej⟩⟨ej∣|K|^2=\sum_js_j^2|e_j\rangle\langle e_j| と書くと ∥Kv∥2=∑jsj2wj\|Kv\|^2=\sum_js_j^2w_j, wj=∣⟨ej,v⟩∣2w_j=|\langle e_j,v\rangle|^2, ∑jwj≤1\sum_jw_j\le1 で、t↦tp/2t\mapsto t^{p/2} は凸で 00 で 00 なので(Jensen の不等式)∥Kv∥p≤∑jwjsjp\|Kv\|^p\le\sum_jw_js_j^p。kk について足し、Bessel の不等式 ∑k∣⟨ej,vk⟩∣2≤1\sum_k|\langle e_j,v_k\rangle|^2\le1 を使います。

本体. KLKL と LL の特異値分解を KL=∑nσn∣φn⟩⟨ψn∣KL=\sum_n\sigma_n|\varphi_n\rangle\langle\psi_n|, L=∑ktk∣vk⟩⟨uk∣L=\sum_kt_k|v_k\rangle\langle u_k| とすると

∥KL∥1=∑n⟨φn,KLψn⟩=∑n∑ktk⟨uk,ψn⟩⟨φn,Kvk⟩≤∑ktk(∑n∣⟨uk,ψn⟩∣2)12(∑n∣⟨φn,Kvk⟩∣2)12,\|KL\|_1=\sum_n\langle\varphi_n,KL\psi_n\rangle=\sum_n\sum_kt_k\langle u_k,\psi_n\rangle\langle\varphi_n,Kv_k\rangle\le\sum_kt_k\Bigl(\sum_n|\langle u_k,\psi_n\rangle|^2\Bigr)^{\frac12}\Bigl(\sum_n|\langle\varphi_n,Kv_k\rangle|^2\Bigr)^{\frac12},

Bessel の不等式から右辺は ∑ktk∥Kvk∥\sum_kt_k\|Kv_k\| 以下で、数列の Hölder の不等式と補題から ≤∥L∥q(∑k∥Kvk∥p)1/p≤∥L∥q∥K∥p\le\|L\|_q\bigl(\sum_k\|Kv_k\|^p\bigr)^{1/p}\le\|L\|_q\|K\|_p です。KLKL はコンパクトなので特異値分解をもち、上の計算は n≤Nn\le N の有限和でも同じ上界を与えるので、∑nσn<∞\sum_n\sigma_n<\infty、すなわち KL∈S1KL\in\mathcal{S}_1 もこれで示されています。

これで 定理 7.4[コンパクト作用素・Hilbert–Schmidt 作用素・トレースクラス] がすべて証明されました。より一般に 1r=1p+1q\frac1r=\frac1p+\frac1q で ∥KL∥r≤∥K∥p∥L∥q\|KL\|_r\le\|K\|_p\|L\|_q も成り立ちますが、これには特異値の積についての Horn の不等式が必要で、ここでは扱いません。

次は散乱理論や Lieb–Thirring 型不等式で使う、「位置の関数 × 運動量の関数」の評価です。U\mathcal{U} を L2(Rd)L^2(\mathbb{R}^d) 上のユニタリな Fourier 変換、W=U−1\mathcal{W}=\mathcal{U}^{-1} とします(定理 4.1[ℝⁿ 上のフーリエ変換])。

命題 6.4(Kato–Seiler–Simon の不等式(p = 2))✓

f,g∈L2(Rd)f,g\in L^2(\mathbb{R}^d) とし、φ∈L1∩L2\varphi\in L^1\cap L^2 に対し f(x)g(−i∇)φ=f⋅W(g⋅Uφ)f(x)g(-i\nabla)\varphi=f\cdot\mathcal{W}(g\cdot\mathcal{U}\varphi) と置く。この作用素は Hilbert–Schmidt 作用素に拡張され、

∥f(x)g(−i∇)∥2=(2π)−d/2∥f∥2∥g∥2.\|f(x)g(-i\nabla)\|_2=(2\pi)^{-d/2}\|f\|_2\|g\|_2 .
証明(命題 6.4)

g∈L1∩L2g\in L^1\cap L^2 なら、Fubini の定理から

W(g Uφ)(x)=(2π)−d∬g(ξ)φ(y)ei(x−y)⋅ξ dy dξ=(2π)−d/2∫(Wg)(x−y)φ(y) dy.\mathcal{W}(g\,\mathcal{U}\varphi)(x)=(2\pi)^{-d}\iint g(\xi)\varphi(y)e^{i(x-y)\cdot\xi}\,dy\,d\xi=(2\pi)^{-d/2}\int(\mathcal{W}g)(x-y)\varphi(y)\,dy .

両辺は g∈L2g\in L^2 について連続です(左辺は ∥g Uφ∥2≤∥Uφ∥∞∥g∥2\|g\,\mathcal{U}\varphi\|_2\le\|\mathcal{U}\varphi\|_\infty\|g\|_2、右辺は Young の不等式 定理 5.1[ℝⁿ 上のフーリエ変換] の ∥h∗φ∥2≤∥h∥2∥φ∥1\|h*\varphi\|_2\le\|h\|_2\|\varphi\|_1)ので、等式は全 g∈L2g\in L^2 で成り立ちます。よって作用素は核 k(x,y)=(2π)−d/2f(x)(Wg)(x−y)k(x,y)=(2\pi)^{-d/2}f(x)(\mathcal{W}g)(x-y) の積分作用素で、Plancherel の定理から

∥k∥L2(R2d)2=(2π)−d∫∣f(x)∣2∫∣(Wg)(x−y)∣2dy dx=(2π)−d∥f∥22∥g∥22.\|k\|_{L^2(\mathbb{R}^{2d})}^2=(2\pi)^{-d}\int|f(x)|^2\int|(\mathcal{W}g)(x-y)|^2dy\,dx=(2\pi)^{-d}\|f\|_2^2\|g\|_2^2 .

定理 6.4[コンパクト作用素・Hilbert–Schmidt 作用素・トレースクラス] から、核が L2L^2 の積分作用素は Hilbert–Schmidt で、ノルムは核の L2L^2 ノルムです。

p=2p=2 の等式と、p=∞p=\infty の自明な評価 ∥f(x)g(−i∇)∥≤∥f∥∞∥g∥∞\|f(x)g(-i\nabla)\|\le\|f\|_\infty\|g\|_\infty の間を補間すると、2≤p<∞2\le p<\infty の評価が得られます。Schatten クラス全体の補間理論を用意する代わりに、必要な形だけを三線定理(定理 5.4[最大値原理の一般化と正則関数族])から直接示します。帯 {0≤Re z≤1}\{0\le\mathrm{Re}\,z\le1\} を S‾\overline{S} と書きます。

補題 6.5(Schatten ノルムの三線補題)✓

2≤p<∞2\le p<\infty, θ=2/p\theta=2/p とする。F:S‾→B(H)F:\overline{S}\to\mathcal{B}(\mathcal{H}) は、任意の u,v∈Hu,v\in\mathcal{H} について z↦⟨u,F(z)v⟩z\mapsto\langle u,F(z)v\rangle が S‾\overline{S} 上で有界かつ連続、内部で正則であるとし、

sup⁡y∈R∥F(iy)∥≤M0,F(1+iy)∈S2,sup⁡y∈R∥F(1+iy)∥2≤M1\sup_{y\in\mathbb{R}}\|F(iy)\|\le M_0,\qquad F(1+iy)\in\mathcal{S}_2,\quad \sup_{y\in\mathbb{R}}\|F(1+iy)\|_2\le M_1

を満たすとする。F(θ)F(\theta) がコンパクトなら F(θ)∈SpF(\theta)\in\mathcal{S}_p で、∥F(θ)∥p≤M01−θM1θ\|F(\theta)\|_p\le M_0^{1-\theta}M_1^{\theta}。

証明の方針: ℓp\ell^p ノルムを、双対的な数列 cnc_n との内積として測ります。z=0z=0 で ℓ∞\ell^\infty 側、z=1z=1 で ℓ2\ell^2 側に対応するよう cn(z)c_n(z) の指数を選び、z=θz=\theta で元の pp に戻ります。

証明(補題 6.5)

定理 6.2 により、有限個の正規直交系 (gn)n≤N(g_n)_{n\le N}, (hn)n≤N(h_n)_{n\le N} について an=⟨gn,F(θ)hn⟩a_n=\langle g_n,F(\theta)h_n\rangle と置き、∑n≤N∣an∣p≤(M01−θM1θ)p\sum_{n\le N}|a_n|^p\le(M_0^{1-\theta}M_1^\theta)^p を示せば十分です(変分表示の右辺の和は有限部分和の上限です)。∥a∥p=(∑∣an∣p)1/p>0\|a\|_p=(\sum|a_n|^p)^{1/p}\gt0 としてよく、p′=p/(p−1)p'=p/(p-1) と置いて

cn=sgn an‾ ∣an∣p−1∥a∥pp−1c_n=\overline{\mathrm{sgn}\,a_n}\,\frac{|a_n|^{p-1}}{\|a\|_p^{p-1}}

とすると ∑n∣cn∣p′=1\sum_n|c_n|^{p'}=1, ∑ncnan=∥a∥p\sum_nc_na_n=\|a\|_p です(sgn a=a/∣a∣\mathrm{sgn}\,a=a/|a|、sgn 0=0\mathrm{sgn}\,0=0)。cn=∣cn∣eiαnc_n=|c_n|e^{i\alpha_n} と書き、

cn(z)=eiαn∣cn∣p′(1−z/2)(cn=0 なら cn(z)=0),Φ(z)=∑n≤Ncn(z)⟨gn,F(z)hn⟩c_n(z)=e^{i\alpha_n}|c_n|^{p'(1-z/2)}\quad(c_n=0\text{ なら }c_n(z)=0),\qquad\Phi(z)=\sum_{n\le N}c_n(z)\langle g_n,F(z)h_n\rangle

と置きます。∣cn(z)∣=∣cn∣p′(1−Re z/2)|c_n(z)|=|c_n|^{p'(1-\mathrm{Re}\,z/2)} は帯の上で有界なので、Φ\Phi は S‾\overline{S} 上で有界連続、内部で正則です。p′(1−θ/2)=p′(1−1/p)=1p'(1-\theta/2)=p'(1-1/p)=1 なので cn(θ)=cnc_n(\theta)=c_n, Φ(θ)=∥a∥p\Phi(\theta)=\|a\|_p です。

Re z=0\mathrm{Re}\,z=0 では ∣cn(iy)∣=∣cn∣p′|c_n(iy)|=|c_n|^{p'} なので ∣Φ(iy)∣≤∑n∣cn∣p′∥F(iy)∥≤M0|\Phi(iy)|\le\sum_n|c_n|^{p'}\|F(iy)\|\le M_0。Re z=1\mathrm{Re}\,z=1 では ∣cn(1+iy)∣2=∣cn∣p′|c_n(1+iy)|^2=|c_n|^{p'} なので、Cauchy–Schwarz の不等式と 定理 6.2(p=2p=2)から

∣Φ(1+iy)∣≤(∑n∣cn(1+iy)∣2)1/2(∑n∣⟨gn,F(1+iy)hn⟩∣2)1/2≤∥F(1+iy)∥2≤M1.|\Phi(1+iy)|\le\Bigl(\sum_n|c_n(1+iy)|^2\Bigr)^{1/2}\Bigl(\sum_n|\langle g_n,F(1+iy)h_n\rangle|^2\Bigr)^{1/2}\le\|F(1+iy)\|_2\le M_1 .

三線定理から ∥a∥p=∣Φ(θ)∣≤M01−θM1θ\|a\|_p=|\Phi(\theta)|\le M_0^{1-\theta}M_1^\theta です。

f(x)g(−i∇)f(x)g(-i\nabla) を一般の LpL^p の関数について定めるため、稠密な部分空間 D={φ∈L2(Rd):Uφ∈Cc∞(Rd)}\mathcal{D}=\{\varphi\in L^2(\mathbb{R}^d):\mathcal{U}\varphi\in C_c^\infty(\mathbb{R}^d)\} を使います。D\mathcal{D} は Schwartz 空間に含まれるので D⊂L1∩L2\mathcal{D}\subset L^1\cap L^2 で、命題 6.4 の定義域にも入っています。gg が局所的に L1∩L2L^1\cap L^2 に属すれば、φ∈D\varphi\in\mathcal{D} について g Uφ∈L1∩L2g\,\mathcal{U}\varphi\in L^1\cap L^2 なので W(g Uφ)∈L2∩L∞\mathcal{W}(g\,\mathcal{U}\varphi)\in L^2\cap L^\infty(∥Wh∥∞≤(2π)−d/2∥h∥1\|\mathcal{W}h\|_\infty\le(2\pi)^{-d/2}\|h\|_1)です。さらに f∈Lpf\in L^p(2≤p≤∞2\le p\le\infty)なら、12=1p+1r\frac12=\frac1p+\frac1r と置いて Hölder の不等式(定理 4.2[L^p 空間と関数解析への導入])と ∥h∥r≤∥h∥22/r∥h∥∞1−2/r\|h\|_r\le\|h\|_2^{2/r}\|h\|_\infty^{1-2/r} から

f(x)g(−i∇)φ:=f⋅W(g Uφ)∈L2f(x)g(-i\nabla)\varphi:=f\cdot\mathcal{W}(g\,\mathcal{U}\varphi)\in L^2

が意味を持ちます。

定理 6.6(Kato–Seiler–Simon の不等式)✓

2≤p<∞2\le p<\infty, f,g∈Lp(Rd)f,g\in L^p(\mathbb{R}^d) とする。D\mathcal{D} 上の作用素 f(x)g(−i∇)f(x)g(-i\nabla) は Sp\mathcal{S}_p の元に一意に拡張され、

∥f(x)g(−i∇)∥p≤(2π)−d/p∥f∥p∥g∥p.\|f(x)g(-i\nabla)\|_p\le(2\pi)^{-d/p}\|f\|_p\|g\|_p .
証明(定理 6.6)

単関数の場合. f=∑iαi1Aif=\sum_i\alpha_i\mathbf{1}_{A_i}, g=∑jβj1Bjg=\sum_j\beta_j\mathbf{1}_{B_j}(有限和、AiA_i は互いに素で測度有限、αi≠0\alpha_i\ne0。Bj,βjB_j,\beta_j も同様)とします。f,g∈L2∩L∞f,g\in L^2\cap L^\infty なので、命題 6.4 から f(x)g(−i∇)f(x)g(-i\nabla) は Hilbert–Schmidt 作用素、とくにコンパクトです。sgn αi=αi/∣αi∣\mathrm{sgn}\,\alpha_i=\alpha_i/|\alpha_i| として

F(z)=∑i,j∣αi∣pz/2∣βj∣pz/2 sgn αi sgn βj  1Ai(x) 1Bj(−i∇)F(z)=\sum_{i,j}|\alpha_i|^{pz/2}|\beta_j|^{pz/2}\,\mathrm{sgn}\,\alpha_i\,\mathrm{sgn}\,\beta_j\;\mathbf{1}_{A_i}(x)\,\mathbf{1}_{B_j}(-i\nabla)

と置きます(1B(−i∇)=W1B U\mathbf{1}_{B}(-i\nabla)=\mathcal{W}\mathbf{1}_B\,\mathcal{U} は射影)。⟨u,F(z)v⟩\langle u,F(z)v\rangle は ecze^{cz} 型の項の有限和なので、帯の上で有界連続、内部で正則です。

  • F(iy)=uy(x)vy(−i∇)F(iy)=u_y(x)v_y(-i\nabla) で、∣uy∣=1∪Ai|u_y|=\mathbf{1}_{\cup A_i}, ∣vy∣=1∪Bj|v_y|=\mathbf{1}_{\cup B_j} です。掛け算作用素と WvyU\mathcal{W}v_y\mathcal{U} のノルムは 11 以下なので ∥F(iy)∥≤1\|F(iy)\|\le1。
  • F(1+iy)=uy′(x)vy′(−i∇)F(1+iy)=u'_y(x)v'_y(-i\nabla) で、∣uy′∣=∣f∣p/2|u'_y|=|f|^{p/2}, ∣vy′∣=∣g∣p/2|v'_y|=|g|^{p/2} です。命題 6.4 から ∥F(1+iy)∥2=(2π)−d/2∥∣f∣p/2∥2∥∣g∣p/2∥2=(2π)−d/2∥f∥pp/2∥g∥pp/2\|F(1+iy)\|_2=(2\pi)^{-d/2}\||f|^{p/2}\|_2\||g|^{p/2}\|_2=(2\pi)^{-d/2}\|f\|_p^{p/2}\|g\|_p^{p/2}。
  • θ=2/p\theta=2/p で ∣αi∣pθ/2sgn αi=αi|\alpha_i|^{p\theta/2}\mathrm{sgn}\,\alpha_i=\alpha_i なので F(θ)=f(x)g(−i∇)F(\theta)=f(x)g(-i\nabla) で、これはコンパクトです。

補題 6.5 から

∥f(x)g(−i∇)∥p≤11−θ((2π)−d/2∥f∥pp/2∥g∥pp/2)2/p=(2π)−d/p∥f∥p∥g∥p.\|f(x)g(-i\nabla)\|_p\le1^{1-\theta}\Bigl((2\pi)^{-d/2}\|f\|_p^{p/2}\|g\|_p^{p/2}\Bigr)^{2/p}=(2\pi)^{-d/p}\|f\|_p\|g\|_p .

一般の場合. 単関数 fk→ff_k\to f, gk→gg_k\to g を、∣fk∣≤∣f∣|f_k|\le|f|, ∣gk∣≤∣g∣|g_k|\le|g| かつほとんどいたるところ収束するように取ると、優収束定理で LpL^p でも収束します。D\mathcal{D} 上で

fk(x)gk(−i∇)−fl(x)gl(−i∇)=(fk−fl)(x)gk(−i∇)+fl(x)(gk−gl)(−i∇)f_k(x)g_k(-i\nabla)-f_l(x)g_l(-i\nabla)=(f_k-f_l)(x)g_k(-i\nabla)+f_l(x)(g_k-g_l)(-i\nabla)

で、右辺の各項に単関数の場合を使うと、Xk=fk(x)gk(−i∇)X_k=f_k(x)g_k(-i\nabla) は Sp\mathcal{S}_p のコーシー列です。定理 6.2 の完備性から Xk→X∈SpX_k\to X\in\mathcal{S}_p で、∥X∥p=lim⁡∥Xk∥p≤(2π)−d/p∥f∥p∥g∥p\|X\|_p=\lim\|X_k\|_p\le(2\pi)^{-d/p}\|f\|_p\|g\|_p。

残るのは φ∈D\varphi\in\mathcal{D} で Xφ=f⋅W(g Uφ)X\varphi=f\cdot\mathcal{W}(g\,\mathcal{U}\varphi) となることです。∥⋅∥≤∥⋅∥p\|\cdot\|\le\|\cdot\|_p から Xkφ→XφX_k\varphi\to X\varphi です。一方

Xkφ−f W(g Uφ)=(fk−f) W(gk Uφ)+f W((gk−g) Uφ)X_k\varphi-f\,\mathcal{W}(g\,\mathcal{U}\varphi)=(f_k-f)\,\mathcal{W}(g_k\,\mathcal{U}\varphi)+f\,\mathcal{W}\bigl((g_k-g)\,\mathcal{U}\varphi\bigr)

で、Uφ\mathcal{U}\varphi は有界で台がコンパクトなので、優収束定理から (gk−g) Uφ→0(g_k-g)\,\mathcal{U}\varphi\to0 は L1L^1 と L2L^2 の両方で成り立ち、W(⋅)\mathcal{W}(\cdot) は L2∩L∞L^2\cap L^\infty で 00 に収束します。W(gk Uφ)\mathcal{W}(g_k\,\mathcal{U}\varphi) の L2L^2 ノルムと L∞L^\infty ノルムは ∣gk∣≤∣g∣|g_k|\le|g| から kk によらず有界です。上の Hölder の評価 ∥fh∥2≤∥f∥p∥h∥r\|fh\|_2\le\|f\|_p\|h\|_r を両項に使えば、右辺は L2L^2 で 00 に収束します。よって XX は D\mathcal{D} 上で f(x)g(−i∇)f(x)g(-i\nabla) に一致し、D\mathcal{D} は稠密なので拡張は一意です。

p→∞p\to\infty の極限では右辺は ∥f∥∞∥g∥∞\|f\|_\infty\|g\|_\infty になり、自明な作用素ノルムの評価に一致します。p<2p\lt2 では同じ形の評価は成り立ちません(注意 6.9)。

gg が LpL^p に入らないが弱 LpL^p に入る場合(典型例は ∣ξ∣−1|\xi|^{-1} で、∣{∣ξ∣−1>t}∣=ωdt−d|\{|\xi|^{-1}\gt t\}|=\omega_dt^{-d}。ωd\omega_d は単位球の体積)には、Sp\mathcal{S}_p に入ることは期待できませんが、特異値の減衰の速さ sn=O(n−1/p)s_n=O(n^{-1/p}) は保たれます。これが Cwikel の評価で、Schrödinger 作用素の束縛状態の個数の評価(CLR 不等式)の核心です。弱 LpL^p の擬ノルム [g]p,∞=sup⁡t>0t ∣{∣g∣>t}∣1/p[g]_{p,\infty}=\sup_{t\gt0}t\,|\{|g|\gt t\}|^{1/p} は 定義 2.1[補間と古典的不等式] のとおりです。まず特異値についての簡単な補題を用意します。

補題 6.7(有界な摂動と特異値)✓

B∈B(H)B\in\mathcal{B}(\mathcal{H}), R∈K(H)R\in\mathcal{K}(\mathcal{H}) で B+RB+R がコンパクトなら、すべての nn で sn(B+R)≤∥B∥+sn(R)s_n(B+R)\le\|B\|+s_n(R)。また R∈S2R\in\mathcal{S}_2 なら sn(R)≤n−1/2∥R∥2s_n(R)\le n^{-1/2}\|R\|_2。

証明(補題 6.7)

K=B+RK=B+R とし、sn(K)>0s_n(K)\gt0 としてよいです。特異値分解 K=∑ksk(K)∣fk⟩⟨ek∣K=\sum_ks_k(K)|f_k\rangle\langle e_k|, R=∑ktk∣fk′⟩⟨ek′∣R=\sum_kt_k|f'_k\rangle\langle e'_k|(定理 7.2[コンパクト作用素・Hilbert–Schmidt 作用素・トレースクラス]。tk=sk(R)t_k=s_k(R))を取ります。nn 次元の空間 L=span{e1,…,en}L=\mathrm{span}\{e_1,\ldots,e_n\} には e1′,…,en−1′e'_1,\ldots,e'_{n-1} のすべて(RR の特異値が n−1n-1 個未満ならあるものすべて)に直交する単位ベクトル xx があります(高々 n−1n-1 個の一次条件)。x∈Lx\in L なので ∥Kx∥2=∑k≤nsk(K)2∣⟨ek,x⟩∣2≥sn(K)2\|Kx\|^2=\sum_{k\le n}s_k(K)^2|\langle e_k,x\rangle|^2\ge s_n(K)^2、また ∥Rx∥2=∑k≥ntk2∣⟨ek′,x⟩∣2≤tn2\|Rx\|^2=\sum_{k\ge n}t_k^2|\langle e'_k,x\rangle|^2\le t_n^2 です。よって sn(K)≤∥Kx∥≤∥Bx∥+∥Rx∥≤∥B∥+sn(R)s_n(K)\le\|Kx\|\le\|Bx\|+\|Rx\|\le\|B\|+s_n(R)。後半は特異値の単調性から nsn(R)2≤∑k≤nsk(R)2≤∥R∥22ns_n(R)^2\le\sum_{k\le n}s_k(R)^2\le\|R\|_2^2 です。

定理 6.8(Cwikel の評価)✓

p>2p\gt2 とし、f∈Lp(Rd)f\in L^p(\mathbb{R}^d)、gg は [g]p,∞<∞[g]_{p,\infty}\lt\infty を満たす可測関数とする。このとき D\mathcal{D} 上の作用素 f(x)g(−i∇)f(x)g(-i\nabla) はコンパクト作用素 AA に一意に拡張され、すべての n≥1n\ge1 で

sn(A)≤Cd,p n−1/p ∥f∥p [g]p,∞,Cd,p=5+8(2π)−d/2(2p2p−4)1/2.s_n(A)\le C_{d,p}\,n^{-1/p}\,\|f\|_p\,[g]_{p,\infty},\qquad C_{d,p}=5+8(2\pi)^{-d/2}\Bigl(\frac{2^p}{2^p-4}\Bigr)^{1/2}.

証明の方針: ff と gg を大きさ 2j2^j, 2k2^k ごとの層に分け、積の大きさ 2j+k2^{j+k} が小さい側は作用素ノルムで、大きい側は Hilbert–Schmidt ノルムで評価します(要点は証明のあとにまとめています)。

証明(定理 6.8)

準備. f⋅W(g Uφ)f\cdot\mathcal{W}(g\,\mathcal{U}\varphi) が φ∈D\varphi\in\mathcal{D} で意味を持つには gg が局所的に L1∩L2L^1\cap L^2 であればよく(上の注意)、それは 補題 2.2[補間と古典的不等式] と p>2p\gt2 から出ます: 半径 ρ\rho の球 BρB_\rho 上で ∫Bρ∣g∣2=∫0∞2t ∣{∣g∣>t}∩Bρ∣ dt≤∫0∞2tmin⁡(∣Bρ∣,[g]p,∞pt−p) dt<∞\int_{B_\rho}|g|^2=\int_0^\infty2t\,|\{|g|\gt t\}\cap B_\rho|\,dt\le\int_0^\infty2t\min(|B_\rho|,[g]_{p,\infty}^pt^{-p})\,dt\lt\infty。両辺は ff と gg について斉次なので、∥f∥p=[g]p,∞=1\|f\|_p=[g]_{p,\infty}=1 としてよいです。

j,k∈Zj,k\in\mathbb{Z} について

Fj={2j≤∣f∣<2j+1},fj=f1Fj,aj=∣Fj∣;Gk={2k≤∣g∣<2k+1},gk=g1GkF_j=\{2^j\le|f|\lt2^{j+1}\},\quad f_j=f\mathbf{1}_{F_j},\quad a_j=|F_j|;\qquad G_k=\{2^k\le|g|\lt2^{k+1}\},\quad g_k=g\mathbf{1}_{G_k}

と置きます。∑j2jpaj≤∫∣f∣p=1\sum_j2^{jp}a_j\le\int|f|^p=1 です。また {∣g∣≥s}=⋂t<s{∣g∣>t}\{|g|\ge s\}=\bigcap_{t\lt s}\{|g|\gt t\} なので ∣{∣g∣≥s}∣≤s−p|\{|g|\ge s\}|\le s^{-p}、とくに ∣Gk∣≤2−kp|G_k|\le2^{-kp} です。整数 mm を固定し、fj≠0f_j\ne0 となる各 jj で gg を「j+k≤mj+k\le m の部分」と「j+k>mj+k\gt m の部分」に分けます:

ej=g 1{∣g∣<2m−j+1}=∑k≤m−jgk,hj=g 1{∣g∣≥2m−j+1}=∑k≥m−j+1gk.e_j=g\,\mathbf{1}_{\{|g|\lt2^{m-j+1}\}}=\sum_{k\le m-j}g_k,\qquad h_j=g\,\mathbf{1}_{\{|g|\ge2^{m-j+1}\}}=\sum_{k\ge m-j+1}g_k .

f=∑jfjf=\sum_jf_j で FjF_j は互いに素なので、φ∈D\varphi\in\mathcal{D} について FjF_j の上で f W(g Uφ)=fjW(ej Uφ)+fjW(hj Uφ)f\,\mathcal{W}(g\,\mathcal{U}\varphi)=f_j\mathcal{W}(e_j\,\mathcal{U}\varphi)+f_j\mathcal{W}(h_j\,\mathcal{U}\varphi) です。

有界な部分. ∣fj∣<2j+1|f_j|\lt2^{j+1} と FjF_j の互いに素であることから

∥∑jfjW(ej Uφ)∥2=∑j∥fjW(ej Uφ)∥2≤∫(∑j4j+1∣ej(ξ)∣2)∣Uφ(ξ)∣2 dξ.\Bigl\|\sum_jf_j\mathcal{W}(e_j\,\mathcal{U}\varphi)\Bigr\|^2=\sum_j\|f_j\mathcal{W}(e_j\,\mathcal{U}\varphi)\|^2\le\int\Bigl(\sum_j4^{j+1}|e_j(\xi)|^2\Bigr)|\mathcal{U}\varphi(\xi)|^2\,d\xi .

g(ξ)≠0g(\xi)\ne0 を固定すると、ej(ξ)≠0e_j(\xi)\ne0 となるのは 2j∣g(ξ)∣<2m+12^j|g(\xi)|\lt2^{m+1} を満たす jj、すなわちある j0j_0 以下のすべての整数で、等比級数から

∑j4j+1∣ej(ξ)∣2=4∣g(ξ)∣2∑j≤j04j=163(2j0∣g(ξ)∣)2<643 4m.\sum_j4^{j+1}|e_j(\xi)|^2=4|g(\xi)|^2\sum_{j\le j_0}4^j=\frac{16}{3}\bigl(2^{j_0}|g(\xi)|\bigr)^2\lt\frac{64}{3}\,4^m .

よって Bmφ=∑jfjW(ej Uφ)B_m\varphi=\sum_jf_j\mathcal{W}(e_j\,\mathcal{U}\varphi) は ∥Bm∥≤(8/3)2m<5⋅2m\|B_m\|\le(8/\sqrt3)2^m\lt5\cdot2^m の有界作用素に拡張されます。

Hilbert–Schmidt の部分. ∥fj∥22≤4j+1aj\|f_j\|_2^2\le4^{j+1}a_j、また p>2p\gt2 から

∥hj∥22=∑k≥m−j+1∥gk∥22≤∑k≥m−j+14k+12−kp=4⋅2(m−j+1)(2−p)1−22−p\|h_j\|_2^2=\sum_{k\ge m-j+1}\|g_k\|_2^2\le\sum_{k\ge m-j+1}4^{k+1}2^{-kp}=\frac{4\cdot2^{(m-j+1)(2-p)}}{1-2^{2-p}}

は有限です。命題 6.4 から Xj=fj(x)hj(−i∇)X_j=f_j(x)h_j(-i\nabla) は Hilbert–Schmidt 作用素で、XjX_j の値域は L2(Fj)L^2(F_j) に含まれます。FjF_j は互いに素なので、任意の正規直交基底 (un)(u_n) について ∥∑j∈JXjun∥2=∑j∈J∥Xjun∥2\|\sum_{j\in J}X_ju_n\|^2=\sum_{j\in J}\|X_ju_n\|^2、したがって有限集合 JJ について ∥∑j∈JXj∥22=∑j∈J∥Xj∥22\|\sum_{j\in J}X_j\|_2^2=\sum_{j\in J}\|X_j\|_2^2 です。

∑j∥Xj∥22=(2π)−d∑j∥fj∥22∥hj∥22≤(2π)−d16⋅22−p1−22−p 2m(2−p)∑j22j2−j(2−p)aj≤(2π)−d642p−4 2m(2−p)\sum_j\|X_j\|_2^2=(2\pi)^{-d}\sum_j\|f_j\|_2^2\|h_j\|_2^2\le(2\pi)^{-d}\frac{16\cdot2^{2-p}}{1-2^{2-p}}\,2^{m(2-p)}\sum_j2^{2j}2^{-j(2-p)}a_j\le(2\pi)^{-d}\frac{64}{2^p-4}\,2^{m(2-p)}

(最後に ∑j2jpaj≤1\sum_j2^{jp}a_j\le1)なので、部分和は S2\mathcal{S}_2 のコーシー列で、定理 6.2 の完備性から Rm=∑jXj∈S2R_m=\sum_jX_j\in\mathcal{S}_2 が定まり、∥Rm∥22≤(2π)−d642p−42m(2−p)\|R_m\|_2^2\le(2\pi)^{-d}\frac{64}{2^p-4}2^{m(2-p)} です。D⊂L1∩L2\mathcal{D}\subset L^1\cap L^2 なので Rmφ=∑jfjW(hj Uφ)R_m\varphi=\sum_jf_j\mathcal{W}(h_j\,\mathcal{U}\varphi)(L2L^2 収束)です。

結論. 以上から D\mathcal{D} 上で f(x)g(−i∇)=Bm+Rmf(x)g(-i\nabla)=B_m+R_m で、右辺の有界作用素を AA とすると、D\mathcal{D} の稠密性から AA は mm によりません。∥A−Rm∥=∥Bm∥<5⋅2m→0\|A-R_m\|=\|B_m\|\lt5\cdot2^m\to0(m→−∞m\to-\infty)で RmR_m はコンパクトなので、命題 3.2[コンパクト作用素・Hilbert–Schmidt 作用素・トレースクラス] から AA もコンパクトです。補題 6.7 から、すべての mm で

sn(A)≤∥Bm∥+n−1/2∥Rm∥2≤5⋅2m+8(2π)−d/2(2p−4)−1/2 2m(1−p/2)n−1/2.s_n(A)\le\|B_m\|+n^{-1/2}\|R_m\|_2\le5\cdot2^m+8(2\pi)^{-d/2}(2^p-4)^{-1/2}\,2^{m(1-p/2)}n^{-1/2}.

2m≤n−1/p<2m+12^m\le n^{-1/p}\lt2^{m+1} となる mm を選ぶと 2−m<2n1/p2^{-m}\lt2n^{1/p} なので 2m(1−p/2)=2m(2−m)p/2<n−1/p2p/2n1/22^{m(1-p/2)}=2^m(2^{-m})^{p/2}\lt n^{-1/p}2^{p/2}n^{1/2} で、

sn(A)<(5+8(2π)−d/2(2p2p−4)1/2)n−1/p.s_n(A)\lt\Bigl(5+8(2\pi)^{-d/2}\Bigl(\frac{2^p}{2^p-4}\Bigr)^{1/2}\Bigr)n^{-1/p}.

証明の要は、二重の 2 進分解 f(x)g(ξ)=∑j,kfj(x)gk(ξ)f(x)g(\xi)=\sum_{j,k}f_j(x)g_k(\xi) を対角線 j+k=mj+k=m で切ることです。j+k≤mj+k\le m の側は大きさが 2j+k≲2m2^{j+k}\lesssim2^m の積の集まりで、xx 側の台が互いに素なので作用素ノルムが 2m2^m 程度に収まります。j+k>mj+k\gt m の側は L2L^2 ノルムが有限で、p>2p\gt2 のおかげで kk についての和 ∑2k(2−p)\sum2^{k(2-p)} が収束します。p=2p=2 では最後の和が発散し、評価は成り立ちません。

注意 6.9

2≤p<∞2\le p\lt\infty の Kato–Seiler–Simon の不等式は 定理 6.6、弱 LpL^p 版は 定理 6.8 で証明しました。p<2p\lt2 では同じ形の評価は成り立たず、Birman–Solomyak の定理(f,gf,g を格子 Zd\mathbb{Z}^d の各単位立方体上の L2L^2 ノルムを並べた ℓp(L2)\ell^p(L^2) 型のノルムで測る)がその代わりになります。こちらは主張のみとします。理由: 格子上の分解と、Sp\mathcal{S}_p(p<1p\lt1 を含む)の準ノルムでの和の評価を別に準備する必要があるためです。文献: B. Simon, Trace Ideals and Their Applications, 第 4 章。

7. 二重作用素積分と Daleckii–Krein 公式

Section titled “7. 二重作用素積分と Daleckii–Krein 公式”

この節では有限次元の Hermite 行列を扱い、定数が次元によらないことに注意を払います。A=∑iλiPiA=\sum_i\lambda_iP_i, B=∑jμjQjB=\sum_j\mu_jQ_j をスペクトル分解(PiP_i, QjQ_j は直交射影で ∑Pi=∑Qj=1\sum P_i=\sum Q_j=1)とし、関数 ϕ\phi に対して

TϕA,B(X)=∑i,jϕ(λi,μj) PiXQjT^{A,B}_\phi(X)=\sum_{i,j}\phi(\lambda_i,\mu_j)\,P_iXQ_j

と置きます(二重作用素積分の有限次元版)。f∈C1(R)f\in C^1(\mathbb{R}) の差分商を f[1](λ,μ)=∫01f′(μ+s(λ−μ)) dsf^{[1]}(\lambda,\mu)=\int_0^1f'(\mu+s(\lambda-\mu))\,ds と定めます(λ≠μ\lambda\ne\mu では (f(λ)−f(μ))/(λ−μ)(f(\lambda)-f(\mu))/(\lambda-\mu)、λ=μ\lambda=\mu では f′(λ)f'(\lambda)。R2\mathbb{R}^2 上連続)。

補題 7.1(二重作用素積分の基本性質)✓

  1. ∥TϕA,B(X)∥2≤max⁡i,j∣ϕ(λi,μj)∣ ∥X∥2\|T^{A,B}_\phi(X)\|_2\le\max_{i,j}|\phi(\lambda_i,\mu_j)|\,\|X\|_2(Hilbert–Schmidt ノルム)。
  2. f∈C1(R)f\in C^1(\mathbb{R}) なら f(A)−f(B)=Tf[1]A,B(A−B)f(A)-f(B)=T^{A,B}_{f^{[1]}}(A-B)。
証明(補題 7.1)

1: tr((PiXQj)∗PkXQl)=tr(QlQjX∗PiPkX)=δikδjl∥PiXQj∥22\mathrm{tr}\bigl((P_iXQ_j)^*P_kXQ_l\bigr)=\mathrm{tr}(Q_lQ_jX^*P_iP_kX)=\delta_{ik}\delta_{jl}\|P_iXQ_j\|_2^2 なので、ブロック PiXQjP_iXQ_j は Hilbert–Schmidt 内積で互いに直交し、∥Tϕ(X)∥22=∑∣ϕij∣2∥PiXQj∥22≤max⁡∣ϕij∣2∥X∥22\|T_\phi(X)\|_2^2=\sum|\phi_{ij}|^2\|P_iXQ_j\|_2^2\le\max|\phi_{ij}|^2\|X\|_2^2。2: Pi(A−B)Qj=(λi−μj)PiQjP_i(A-B)Q_j=(\lambda_i-\mu_j)P_iQ_j で、f[1](λi,μj)(λi−μj)=f(λi)−f(μj)f^{[1]}(\lambda_i,\mu_j)(\lambda_i-\mu_j)=f(\lambda_i)-f(\mu_j)(λi=μj\lambda_i=\mu_j でも両辺 00)なので

Tf[1](A−B)=∑i,j(f(λi)−f(μj))PiQj=f(A)∑jQj−∑iPif(B)=f(A)−f(B).T_{f^{[1]}}(A-B)=\sum_{i,j}\bigl(f(\lambda_i)-f(\mu_j)\bigr)P_iQ_j=f(A)\sum_jQ_j-\sum_iP_if(B)=f(A)-f(B).

定理 7.2(Hilbert–Schmidt ノルムでのリプシッツ性)✓

f∈C1(R)f\in C^1(\mathbb{R}) で f′f' が有界なら、任意の次元の Hermite 行列 A,BA,B について ∥f(A)−f(B)∥2≤∥f′∥∞∥A−B∥2\|f(A)-f(B)\|_2\le\|f'\|_\infty\|A-B\|_2。

証明(定理 7.2)

∣f[1]∣≤∥f′∥∞|f^{[1]}|\le\|f'\|_\infty なので 補題 7.1 の 1 と 2 からただちに出ます。

作用素ノルムでは同じ評価は成り立ちません。f(t)=∣t∣f(t)=|t| はリプシッツですが、∥∣A∣−∣B∣∥≤C∥A−B∥\||A|-|B|\|\le C\|A-B\| となる次元によらない CC は存在しないことが知られています(Kato, McIntosh, Davies)。Hilbert–Schmidt ノルムで評価が自明になるのは、補題 7.1 の 1 のように「Schur 積」がブロックの直交性と両立するからです。

定理 7.3(Daleckii–Krein 公式)✓

A,BA,B を Hermite 行列、A=∑iλiPiA=\sum_i\lambda_iP_i とし、f∈C1(R)f\in C^1(\mathbb{R}) とする。このとき

ddtf(A+tB)∣t=0=Tf[1]A,A(B)=∑i,jf[1](λi,λj)PiBPj.\frac{d}{dt}f(A+tB)\Big|_{t=0}=T^{A,A}_{f^{[1]}}(B)=\sum_{i,j}f^{[1]}(\lambda_i,\lambda_j)P_iBP_j .

とくに ddttrf(A+tB)∣t=0=tr(f′(A)B)\frac{d}{dt}\mathrm{tr}f(A+tB)|_{t=0}=\mathrm{tr}\bigl(f'(A)B\bigr)。

証明(定理 7.3)

補題 7.1 の 2 を A+tBA+tB と AA に使うと 1t(f(A+tB)−f(A))=Tf[1]A+tB,A(B)\frac1t\bigl(f(A+tB)-f(A)\bigr)=T^{A+tB,A}_{f^{[1]}}(B) です。t→0t\to0 でこれが Tf[1]A,A(B)T^{A,A}_{f^{[1]}}(B) に収束することを示します。ϕ(λ,μ)=p(λ)q(μ)\phi(\lambda,\mu)=p(\lambda)q(\mu) の形なら TϕA′,A(X)=p(A′)Xq(A)T^{A',A}_\phi(X)=p(A')Xq(A) は A′A' について連続です。連続関数 f[1]f^{[1]} は、スペクトルを含む正方形の上で、こうした積の有限和(2 変数多項式)で一様に近似できます(Stone–Weierstrass の定理 定理 4.3[Arzelà–Ascoli と Stone–Weierstrass])。補題 7.1 の 1 の一様な評価と合わせると(ε/3\varepsilon/3 論法)、A′↦Tf[1]A′,A(B)A'\mapsto T^{A',A}_{f^{[1]}}(B) は連続です。トレースについては、tr(PiBPj)=tr(PjPiB)=0\mathrm{tr}(P_iBP_j)=\mathrm{tr}(P_jP_iB)=0(i≠ji\ne j)と f[1](λi,λi)=f′(λi)f^{[1]}(\lambda_i,\lambda_i)=f'(\lambda_i) から tr Tf[1](B)=∑if′(λi)tr(PiB)=tr(f′(A)B)\mathrm{tr}\,T_{f^{[1]}}(B)=\sum_if'(\lambda_i)\mathrm{tr}(P_iB)=\mathrm{tr}(f'(A)B) です。

注意 7.4

無限次元では、2 つの可換な射影値測度 X↦EA(⋅)XX\mapsto E_A(\cdot)X と X↦XEB(⋅)X\mapsto XE_B(\cdot) の積測度で ∫ ⁣ ⁣∫ϕ(λ,μ) dEA(λ) X dEB(μ)\int\!\!\int\phi(\lambda,\mu)\,dE_A(\lambda)\,X\,dE_B(\mu) を Hilbert–Schmidt 作用素の空間上で定義できます(Birman–Solomyak の二重作用素積分)。定理 7.2 は A−B∈S2A-B\in\mathcal{S}_2 の自己共役作用素にそのまま拡張されます。S1\mathcal{S}_1 や作用素ノルムで同じことをするには ϕ\phi を積の積分で表す必要があり、ff が Besov 空間に属するといった条件が現れます(Peller の理論)。

定理 8.1(Klein の不等式)✓

A,BA,B を Hermite 行列、ff をスペクトルを含む区間 II 上の凸関数で、BB のスペクトル上で微分可能とする。このとき

tr[f(A)−f(B)−(A−B)f′(B)]≥0.\mathrm{tr}\bigl[f(A)-f(B)-(A-B)f'(B)\bigr]\ge0 .

ff が狭義凸なら、等号は A=BA=B のときに限る。

証明(定理 8.1)

A=∑iai∣ui⟩⟨ui∣A=\sum_ia_i|u_i\rangle\langle u_i|, B=∑jbj∣vj⟩⟨vj∣B=\sum_jb_j|v_j\rangle\langle v_j|(正規直交固有基底)とし、cij=∣⟨ui,vj⟩∣2c_{ij}=|\langle u_i,v_j\rangle|^2 と置くと ∑icij=∑jcij=1\sum_ic_{ij}=\sum_jc_{ij}=1 です。trf(A)=∑i,jcijf(ai)\mathrm{tr}f(A)=\sum_{i,j}c_{ij}f(a_i), trf(B)=∑i,jcijf(bj)\mathrm{tr}f(B)=\sum_{i,j}c_{ij}f(b_j), tr((A−B)f′(B))=∑i,jcij(ai−bj)f′(bj)\mathrm{tr}((A-B)f'(B))=\sum_{i,j}c_{ij}(a_i-b_j)f'(b_j) なので

tr[f(A)−f(B)−(A−B)f′(B)]=∑i,jcij[f(ai)−f(bj)−(ai−bj)f′(bj)]≥0\mathrm{tr}\bigl[f(A)-f(B)-(A-B)f'(B)\bigr]=\sum_{i,j}c_{ij}\bigl[f(a_i)-f(b_j)-(a_i-b_j)f'(b_j)\bigr]\ge0

(凸関数は接線の上にある)。狭義凸で等号なら、cij>0c_{ij}>0 となる組で ai=bja_i=b_j なので Avj=∑iai⟨ui,vj⟩ui=bjvjAv_j=\sum_ia_i\langle u_i,v_j\rangle u_i=b_jv_j、すなわち A=BA=B です。

定理 8.2(Gibbs の変分原理)✓

HH を Hermite 行列、Z=tr e−HZ=\mathrm{tr}\,e^{-H} とする。密度行列 ρ\rho(ρ≥0\rho\ge0, trρ=1\mathrm{tr}\rho=1)について

tr(ρH)+tr(ρlog⁡ρ)≥−log⁡Z\mathrm{tr}(\rho H)+\mathrm{tr}(\rho\log\rho)\ge-\log Z

(0log⁡0=00\log0=0)であり、等号は ρ=e−H/Z\rho=e^{-H}/Z のときに限る。

証明(定理 8.2)

σ=e−H/Z\sigma=e^{-H}/Z は正定値で log⁡σ=−H−log⁡Z\log\sigma=-H-\log Z です。f(t)=tlog⁡tf(t)=t\log t は [0,∞)[0,\infty) 上で狭義凸、(0,∞)(0,\infty) で f′(t)=log⁡t+1f'(t)=\log t+1 です。定理 8.1 を A=ρA=\rho, B=σB=\sigma に使うと

tr(ρlog⁡ρ)−tr(σlog⁡σ)−tr((ρ−σ)(log⁡σ+1))=tr(ρlog⁡ρ)−tr(ρlog⁡σ)≥0\mathrm{tr}(\rho\log\rho)-\mathrm{tr}(\sigma\log\sigma)-\mathrm{tr}\bigl((\rho-\sigma)(\log\sigma+1)\bigr)=\mathrm{tr}(\rho\log\rho)-\mathrm{tr}(\rho\log\sigma)\ge0

(trρ=trσ=1\mathrm{tr}\rho=\mathrm{tr}\sigma=1)。−tr(ρlog⁡σ)=tr(ρH)+log⁡Z-\mathrm{tr}(\rho\log\sigma)=\mathrm{tr}(\rho H)+\log Z を代入すれば主張で、等号条件も Klein の不等式からです。

左辺 tr(ρlog⁡ρ)−tr(ρlog⁡σ)\mathrm{tr}(\rho\log\rho)-\mathrm{tr}(\rho\log\sigma) は相対エントロピー S(ρ∥σ)S(\rho\|\sigma) で、その非負性が統計力学の自由エネルギー最小原理です。

系 8.3(Peierls–Bogoliubov の不等式)✓

Hermite 行列 A,BA,B について

tr eA+Btr eA≥exp⁡(tr(eAB)tr eA).\frac{\mathrm{tr}\,e^{A+B}}{\mathrm{tr}\,e^A}\ge\exp\Bigl(\frac{\mathrm{tr}(e^AB)}{\mathrm{tr}\,e^A}\Bigr).
証明(系 8.3)

定理 8.2 を H=−(A+B)H=-(A+B) と ρ=eA/tr eA\rho=e^A/\mathrm{tr}\,e^A に使います。tr(ρlog⁡ρ)=tr(ρA)−log⁡tr eA\mathrm{tr}(\rho\log\rho)=\mathrm{tr}(\rho A)-\log\mathrm{tr}\,e^A なので

−tr(ρA)−tr(ρB)+tr(ρA)−log⁡tr eA≥−log⁡tr eA+B,-\mathrm{tr}(\rho A)-\mathrm{tr}(\rho B)+\mathrm{tr}(\rho A)-\log\mathrm{tr}\,e^A\ge-\log\mathrm{tr}\,e^{A+B},

すなわち log⁡tr eA+B−log⁡tr eA≥tr(ρB)\log\mathrm{tr}\,e^{A+B}-\log\mathrm{tr}\,e^A\ge\mathrm{tr}(\rho B) です。

統計力学の言葉では、A=−βH0A=-\beta H_0, B=−βVB=-\beta V として「H0+VH_0+V の自由エネルギーは、H0H_0 の自由エネルギーに VV の H0H_0 の Gibbs 状態での期待値を足したもの以下」(Bogoliubov の不等式)で、平均場近似が自由エネルギーの上界を与える理由です。

Peierls–Bogoliubov の不等式は tr eA+B\mathrm{tr}\,e^{A+B} の下界でした。上界を与えるのが Golden–Thompson の不等式です。証明の道具は、非正規行列の固有値と特異値を比べる Weyl の不等式と、Lie–Trotter の積公式です。以下、T∈Mn(C)T\in M_n(\mathbb{C}) の特異値 s1(T)≥⋯≥sn(T)≥0s_1(T)\ge\cdots\ge s_n(T)\ge0 は T∗TT^*T の固有値の非負平方根とします。

補題 8.4(Weyl の優越不等式)✓

T∈Mn(C)T\in M_n(\mathbb{C}) の固有値を重複度込みで ∣λ1∣≥⋯≥∣λn∣|\lambda_1|\ge\cdots\ge|\lambda_n| と並べ、si=si(T)s_i=s_i(T) とする。このとき

  1. 1≤k≤n1\le k\le n について ∏i≤k∣λi∣≤∏i≤ksi\prod_{i\le k}|\lambda_i|\le\prod_{i\le k}s_i。
  2. 任意の p>0p>0 について ∑i∣λi∣p≤∑isip\sum_i|\lambda_i|^p\le\sum_is_i^p。
証明(補題 8.4)

1. Schur 分解(定理 3.1[行列の分解と標準形])で T=URU∗T=URU^*、RR は上三角で対角成分が順に λ1,…,λn\lambda_1,\ldots,\lambda_n とします。T∗T=UR∗RU∗T^*T=UR^*RU^* なので si(T)=si(R)s_i(T)=s_i(R) です。E:Ck→CnE:\mathbb{C}^k\to\mathbb{C}^n を最初の kk 座標への埋め込み、Rk=E∗RER_k=E^*RE(左上の k×kk\times k ブロック)とします。RR は上三角なので最初の kk 個の基底ベクトルの張る部分空間を保ち、RE=ERkRE=ER_k です。E∗E=1kE^*E=1_k と合わせて

Rk∗Rk=Rk∗E∗ERk=E∗R∗RE.R_k^*R_k=R_k^*E^*ER_k=E^*R^*RE .

Courant–Fischer の定理(定理 3.1[行列の固有値評価と不等式])から、i≤ki\le k について

si(Rk)2=max⁡dim⁡S=i min⁡x∈S,∥x∥=1⟨R∗REx,Ex⟩≤si(R)2s_i(R_k)^2=\max_{\dim S=i}\ \min_{x\in S,\|x\|=1}\langle R^*REx,Ex\rangle\le s_i(R)^2

です(S⊂CkS\subset\mathbb{C}^k が ii 次元なら ES⊂CnES\subset\mathbb{C}^n も ii 次元で、∥Ex∥=∥x∥\|Ex\|=\|x\|)。RkR_k は対角成分 λ1,…,λk\lambda_1,\ldots,\lambda_k の上三角行列なので

∏i≤k∣λi∣=∣det⁡Rk∣=det⁡(Rk∗Rk)1/2=∏i≤ksi(Rk)≤∏i≤ksi.\prod_{i\le k}|\lambda_i|=|\det R_k|=\det(R_k^*R_k)^{1/2}=\prod_{i\le k}s_i(R_k)\le\prod_{i\le k}s_i .

2. λi≠0\lambda_i\ne0 となる ii の個数を rr とします(r=0r=0 なら自明)。1 から i≤ri\le r で si>0s_i>0 です。xi=log⁡∣λi∣x_i=\log|\lambda_i|, yi=log⁡siy_i=\log s_i(i≤ri\le r)はどちらも非増加で、1 から Dk=∑i≤k(yi−xi)≥0D_k=\sum_{i\le k}(y_i-x_i)\ge0 です。g(u)=epug(u)=e^{pu} は凸で増加なので、弦の傾き ci=g(yi)−g(xi)yi−xic_i=\frac{g(y_i)-g(x_i)}{y_i-x_i}(xi=yix_i=y_i なら g′(xi)g'(x_i))は両端点のそれぞれについて非減少です(凸関数の三弦補題)。xi≥xi+1x_i\ge x_{i+1}, yi≥yi+1y_i\ge y_{i+1} なので c1≥c2≥⋯≥cr≥0c_1\ge c_2\ge\cdots\ge c_r\ge0。Abel の総和法(D0=0D_0=0)で

∑i≤r(g(yi)−g(xi))=∑i≤rci(Di−Di−1)=∑i<r(ci−ci+1)Di+crDr≥0.\sum_{i\le r}\bigl(g(y_i)-g(x_i)\bigr)=\sum_{i\le r}c_i(D_i-D_{i-1})=\sum_{i<r}(c_i-c_{i+1})D_i+c_rD_r\ge0 .

よって ∑i∣λi∣p=∑i≤r∣λi∣p≤∑i≤rsip≤∑isip\sum_i|\lambda_i|^p=\sum_{i\le r}|\lambda_i|^p\le\sum_{i\le r}s_i^p\le\sum_is_i^p です。

補題 8.5(冪のトレースの比較)✓

X,YX,Y を正定値行列、m≥1m\ge1 を整数とすると tr((XY)2m)≤tr((X2Y2)m)\mathrm{tr}\bigl((XY)^{2m}\bigr)\le\mathrm{tr}\bigl((X^2Y^2)^m\bigr)。

証明(補題 8.5)

XY=X1/2(X1/2YX1/2)X−1/2XY=X^{1/2}\bigl(X^{1/2}YX^{1/2}\bigr)X^{-1/2} は正定値行列 X1/2YX1/2X^{1/2}YX^{1/2} に相似なので、固有値 λi\lambda_i は正です。Schur 分解の上三角行列の冪を考えれば (XY)2m(XY)^{2m} の固有値は λi2m\lambda_i^{2m} で、補題 8.4 の 2 を p=2mp=2m で使うと

tr((XY)2m)=∑iλi2m≤∑isi(XY)2m=tr(((XY)∗XY)m)=tr((YX2Y)m).\mathrm{tr}\bigl((XY)^{2m}\bigr)=\sum_i\lambda_i^{2m}\le\sum_is_i(XY)^{2m}=\mathrm{tr}\bigl(((XY)^*XY)^m\bigr)=\mathrm{tr}\bigl((YX^2Y)^m\bigr).

(YX2Y)m=Y(X2Y2)m−1X2Y(YX^2Y)^m=Y(X^2Y^2)^{m-1}X^2Y とトレースの巡回性から、右辺は tr((X2Y2)m)\mathrm{tr}\bigl((X^2Y^2)^m\bigr) です。

定理 8.6(Golden–Thompson の不等式)✓

Hermite 行列 A,BA,B について tr eA+B≤tr(eAeB)\mathrm{tr}\,e^{A+B}\le\mathrm{tr}\bigl(e^Ae^B\bigr)。

証明(定理 8.6)

k≥1k\ge1 を固定し、X=eA/2kX=e^{A/2^k}, Y=eB/2kY=e^{B/2^k}(正定値)と置きます。補題 8.5 を正定値行列の組 (X2j,Y2j)(X^{2^j},Y^{2^j}) と m=2k−j−1m=2^{k-j-1} に j=0,1,…,k−1j=0,1,\ldots,k-1 の順で使うと

tr((XY)2k)≤tr((X2Y2)2k−1)≤tr((X4Y4)2k−2)≤⋯≤tr(X2kY2k)=tr(eAeB).\mathrm{tr}\bigl((XY)^{2^k}\bigr)\le\mathrm{tr}\bigl((X^2Y^2)^{2^{k-1}}\bigr)\le\mathrm{tr}\bigl((X^4Y^4)^{2^{k-2}}\bigr)\le\cdots\le\mathrm{tr}\bigl(X^{2^k}Y^{2^k}\bigr)=\mathrm{tr}\bigl(e^Ae^B\bigr).

左辺は tr((eA/2keB/2k)2k)\mathrm{tr}\bigl((e^{A/2^k}e^{B/2^k})^{2^k}\bigr) で、Lie–Trotter の積公式(命題 2.2[行列 Lie 群と指数写像])を部分列 2k2^k に使えば (eA/2keB/2k)2k→eA+B(e^{A/2^k}e^{B/2^k})^{2^k}\to e^{A+B} です。トレースは連続なので k→∞k\to\infty として主張を得ます。

A,BA,B が可換なら両辺は等しく、Golden–Thompson の不等式は「非可換性は tr eA+B\mathrm{tr}\,e^{A+B} を小さくする方向にしか働かない」と読めます。系 8.3 と合わせると、tr eA+B\mathrm{tr}\,e^{A+B} は可換な場合の公式で上下から挟まれます。3 つの行列では tr eA+B+C≤∣tr(eAeBeC)∣\mathrm{tr}\,e^{A+B+C}\le|\mathrm{tr}(e^Ae^Be^C)| は一般には成り立たず、Lieb はその代わりになる三重行列の不等式を与えました。

注意 8.7

補題 8.5 の反復で得た tr(XY)2k≤tr(X2kY2k)\mathrm{tr}(XY)^{2^k}\le\mathrm{tr}(X^{2^k}Y^{2^k}) は、実際には任意の正の整数冪で成り立ち、さらに tr((Y1/2XY1/2)rp)≤tr((Yr/2XrYr/2)p)\mathrm{tr}\bigl((Y^{1/2}XY^{1/2})^{rp}\bigr)\le\mathrm{tr}\bigl((Y^{r/2}X^rY^{r/2})^p\bigr)(r≥1r\ge1, p>0p>0)に一般化されます(Araki–Lieb–Thirring の不等式)。こちらは外積冪を使った特異値の対数的優越が必要で、この記事では扱いません。Golden–Thompson の不等式の証明には 2 の冪の場合で足ります。もう一つの柱である Lieb の凹性定理は、次節で作用素凸関数を用意して証明します。

9. 作用素凸関数と Lieb の凹性定理

Section titled “9. 作用素凸関数と Lieb の凹性定理”

この節では有限次元で考え、MnM_n は nn 次複素正方行列の全体、P,Q>0P,Q>0 は正定値を表します。Hermite 行列の間の不等式はすべて Löwner 順序です。

定義 9.1(作用素凸関数・同時凹写像)↗

区間 II 上の連続関数 ff が作用素凸であるとは、任意の nn と、スペクトルが II に含まれる Hermite 行列 A,B∈MnA,B\in M_n、0≤λ≤10\le\lambda\le1 について

f(λA+(1−λ)B)≤λf(A)+(1−λ)f(B)f\bigl(\lambda A+(1-\lambda)B\bigr)\le\lambda f(A)+(1-\lambda)f(B)

が成り立つことをいう。−f-f が作用素凸のとき ff を作用素凹という。正定値行列の組 (P,Q)(P,Q) から Hermite 行列への写像 Φ\Phi が同時凹であるとは、Φ(λP1+(1−λ)P2,λQ1+(1−λ)Q2)≥λΦ(P1,Q1)+(1−λ)Φ(P2,Q2)\Phi\bigl(\lambda P_1+(1-\lambda)P_2,\lambda Q_1+(1-\lambda)Q_2\bigr)\ge\lambda\Phi(P_1,Q_1)+(1-\lambda)\Phi(P_2,Q_2) が成り立つことをいう。

λA+(1−λ)B\lambda A+(1-\lambda)B の固有値は A,BA,B の最小固有値の凸結合以上、最大固有値の凸結合以下(Rayleigh 商)なので、左辺の ff は定義されています。行列の不等式 M≥NM\ge N は「全ベクトル xx で ⟨x,Mx⟩≥⟨x,Nx⟩\langle x,Mx\rangle\ge\langle x,Nx\rangle」なので、凹性は各 xx ごとの 2 次形式の凹性に帰着します。とくに、各 xx について ⟨x,Φ(⋅)x⟩\langle x,\Phi(\cdot)x\rangle が線形な関数族の下限や、凹関数の正の重みによる積分で書ければ、Φ\Phi は凹です。以下ではこの原理を繰り返し使います。

補題 9.2(逆行列の変分表示と作用素凸性)✓

P>0P>0 と x∈Cnx\in\mathbb{C}^n について

⟨x,P−1x⟩=max⁡y∈Cn(2 Re⟨x,y⟩−⟨y,Py⟩)\langle x,P^{-1}x\rangle=\max_{y\in\mathbb{C}^n}\bigl(2\,\mathrm{Re}\langle x,y\rangle-\langle y,Py\rangle\bigr)

であり、最大は y=P−1xy=P^{-1}x で達成される。したがって t↦t−1t\mapsto t^{-1} は (0,∞)(0,\infty) 上の作用素凸関数である。

証明(補題 9.2)

展開すると

⟨y−P−1x, P(y−P−1x)⟩=⟨y,Py⟩−2 Re⟨x,y⟩+⟨x,P−1x⟩≥0\bigl\langle y-P^{-1}x,\,P(y-P^{-1}x)\bigr\rangle=\langle y,Py\rangle-2\,\mathrm{Re}\langle x,y\rangle+\langle x,P^{-1}x\rangle\ge0

で、等号は y=P−1xy=P^{-1}x のときです。右辺の括弧の中は各 yy について PP のアフィン関数なので、その最大値は PP の凸関数です。

補題 9.3(並列和の同時凹性)✓

P,Q>0P,Q>0 について P:Q=(P−1+Q−1)−1P:Q=(P^{-1}+Q^{-1})^{-1} と置くと、

⟨x,(P:Q)x⟩=min⁡y∈Cn(⟨y,Py⟩+⟨x−y,Q(x−y)⟩).\langle x,(P:Q)x\rangle=\min_{y\in\mathbb{C}^n}\bigl(\langle y,Py\rangle+\langle x-y,Q(x-y)\rangle\bigr).

したがって (P,Q)↦P:Q(P,Q)\mapsto P:Q は同時凹である。

証明(補題 9.3)

M=P:QM=P:Q, y0=P−1Mxy_0=P^{-1}Mx と置きます。(P−1+Q−1)M=1(P^{-1}+Q^{-1})M=1 から x−y0=Q−1Mxx-y_0=Q^{-1}Mx なので Py0=Q(x−y0)=MxPy_0=Q(x-y_0)=Mx です。y=y0+hy=y_0+h と書くと、交差項 2 Re⟨h,Py0−Q(x−y0)⟩2\,\mathrm{Re}\langle h,Py_0-Q(x-y_0)\rangle が消えて

⟨y,Py⟩+⟨x−y,Q(x−y)⟩=⟨y0,Py0⟩+⟨x−y0,Q(x−y0)⟩+⟨h,(P+Q)h⟩.\langle y,Py\rangle+\langle x-y,Q(x-y)\rangle=\langle y_0,Py_0\rangle+\langle x-y_0,Q(x-y_0)\rangle+\langle h,(P+Q)h\rangle .

最後の項は非負で、h=0h=0 での値は ⟨P−1Mx,Mx⟩+⟨Q−1Mx,Mx⟩=⟨(P−1+Q−1)Mx,Mx⟩=⟨x,Mx⟩\langle P^{-1}Mx,Mx\rangle+\langle Q^{-1}Mx,Mx\rangle=\langle(P^{-1}+Q^{-1})Mx,Mx\rangle=\langle x,Mx\rangle です。右辺の括弧の中は各 yy について (P,Q)(P,Q) の線形関数なので、その最小値は同時凹です。

定理 9.4(冪関数の作用素凹性)✓

0≤α≤10\le\alpha\le1 なら t↦tαt\mapsto t^\alpha は [0,∞)[0,\infty) 上の作用素凹関数である。また log⁡\log は (0,∞)(0,\infty) 上の作用素凹関数である。

証明(定理 9.4)

α=0,1\alpha=0,1 は自明なので 0<α<10<\alpha<1 とします。s>0s>0 について A↦A(s+A)−1=1−s(s+A)−1A\mapsto A(s+A)^{-1}=1-s(s+A)^{-1} は、補題 9.2 を s+As+A に使えば作用素凹です。定理 3.3 の証明と同じく、補題 3.2 から

⟨x,Aαx⟩=sin⁡παπ∫0∞sα−1⟨x,A(s+A)−1x⟩ ds\langle x,A^\alpha x\rangle=\frac{\sin\pi\alpha}{\pi}\int_0^\infty s^{\alpha-1}\langle x,A(s+A)^{-1}x\rangle\,ds

で、被積分関数は各 ss で AA の凹関数、重みは正なので、左辺は AA の凹関数です。log⁡\log については、系 3.4 の証明と同じく (Aα−1)/α→log⁡A(A^\alpha-1)/\alpha\to\log A が α↓0\alpha\downarrow0 でノルム収束し、凹性の不等式は極限で保たれます。

一方 t2t^2 は作用素凸ですが(A2+B22−(A+B2)2=(A−B2)2≥0\frac{A^2+B^2}2-\bigl(\frac{A+B}2\bigr)^2=\bigl(\frac{A-B}2\bigr)^2\ge0)、例 3.5 のとおり作用素単調ではありません。t3t^3 は [0,∞)[0,\infty) 上で作用素凸ですらありません。作用素凸性も作用素単調性も、凸性・単調性よりずっと強い条件です。

次の対象は、正定値行列の重み付き幾何平均

P#tQ=P1/2(P−1/2QP−1/2)tP1/2(0≤t≤1)P\#_tQ=P^{1/2}\bigl(P^{-1/2}QP^{-1/2}\bigr)^tP^{1/2}\qquad(0\le t\le1)

です。P,QP,Q が可換なら P1−tQtP^{1-t}Q^t に一致します。

定理 9.5(重み付き幾何平均の同時凹性)✓

0<t<10<t<1 と P,Q>0P,Q>0 について

P#tQ=sin⁡πtπ∫0∞st−1(P−1+sQ−1)−1dsP\#_tQ=\frac{\sin\pi t}{\pi}\int_0^\infty s^{t-1}\bigl(P^{-1}+sQ^{-1}\bigr)^{-1}ds

が成り立つ。したがって (P,Q)↦P#tQ(P,Q)\mapsto P\#_tQ は同時凹である。

証明(定理 9.5)

Z=P−1/2QP−1/2>0Z=P^{-1/2}QP^{-1/2}>0 と置きます。ZZ の固有基底で見ると、各固有値 zz について 補題 3.2 が zt=sin⁡πtπ∫0∞st−1zs+zdsz^t=\frac{\sin\pi t}{\pi}\int_0^\infty s^{t-1}\frac{z}{s+z}ds を与えるので、行列の成分ごとの積分として Zt=sin⁡πtπ∫0∞st−1Z(s+Z)−1dsZ^t=\frac{\sin\pi t}{\pi}\int_0^\infty s^{t-1}Z(s+Z)^{-1}ds です。成分ごとの積分は定数行列による左右からの積と交換するので、

P#tQ=sin⁡πtπ∫0∞st−1P1/2Z(s+Z)−1P1/2 ds.P\#_tQ=\frac{\sin\pi t}{\pi}\int_0^\infty s^{t-1}P^{1/2}Z(s+Z)^{-1}P^{1/2}\,ds .

s>0s>0 では Z(s+Z)−1=(1+sZ−1)−1Z(s+Z)^{-1}=(1+sZ^{-1})^{-1} で、Z−1=P1/2Q−1P1/2Z^{-1}=P^{1/2}Q^{-1}P^{1/2} と (CDC)−1=C−1D−1C−1(CDC)^{-1}=C^{-1}D^{-1}C^{-1} から

P1/2(1+sZ−1)−1P1/2=(P−1/2(1+sZ−1)P−1/2)−1=(P−1+sQ−1)−1.P^{1/2}(1+sZ^{-1})^{-1}P^{1/2}=\bigl(P^{-1/2}(1+sZ^{-1})P^{-1/2}\bigr)^{-1}=\bigl(P^{-1}+sQ^{-1}\bigr)^{-1}.

これで積分表示が示されました。被積分関数は P:(s−1Q)P:(s^{-1}Q) で、(P,Q)↦(P,s−1Q)(P,Q)\mapsto(P,s^{-1}Q) は線形なので、補題 9.3 から各 ss で同時凹です。2 次形式 ⟨x,⋅ x⟩\langle x,\cdot\,x\rangle を取って正の重みで積分すれば、同時凹性が積分に移ります。

Lieb の凹性定理は、この同時凹性を「行列に左から掛ける作用素」と「右から掛ける作用素」に使うことで得られます(Ando の証明)。

定理 9.6(Lieb の凹性定理)✓

0≤t≤10\le t\le1, K∈MnK\in M_n とする。半正定値行列の組 (A,B)(A,B) の関数

(A,B)⟼tr(K∗AtKB1−t)(A,B)\longmapsto\mathrm{tr}\bigl(K^*A^tKB^{1-t}\bigr)

は同時凹である(A0=B0=1A^0=B^0=1 と約束する)。

証明(定理 9.6)

t=0,1t=0,1 では線形なので、0<t<10<t<1 とします。まず A,B>0A,B>0 とします。MnM_n に Hilbert–Schmidt 内積 ⟨X,Y⟩=tr(X∗Y)\langle X,Y\rangle=\mathrm{tr}(X^*Y) を入れ、MnM_n 上の線形作用素 LA(X)=AXL_A(X)=AX, RB(X)=XBR_B(X)=XB を考えます。A=∑iaiuiui∗A=\sum_ia_iu_iu_i^*, B=∑jbjwjwj∗B=\sum_jb_jw_jw_j^* をスペクトル分解とすると、Eij=uiwj∗E_{ij}=u_iw_j^* は tr(Eij∗Ekl)=δikδjl\mathrm{tr}(E_{ij}^*E_{kl})=\delta_{ik}\delta_{jl} から MnM_n の正規直交基底で、

LAEij=aiEij,RBEij=ui(Bwj)∗=bjEij.L_AE_{ij}=a_iE_{ij},\qquad R_BE_{ij}=u_i(Bw_j)^*=b_jE_{ij}.

よって LA,RBL_A,R_B は同じ正規直交基底で対角化される正定値作用素で、RB−1/2LARB−1/2Eij=(ai/bj)EijR_B^{-1/2}L_AR_B^{-1/2}E_{ij}=(a_i/b_j)E_{ij} から

(RB#tLA)Eij=bj1/2(aibj)tbj1/2Eij=aitbj1−tEij=LAtRB1−tEij.\bigl(R_B\#_tL_A\bigr)E_{ij}=b_j^{1/2}\Bigl(\frac{a_i}{b_j}\Bigr)^tb_j^{1/2}E_{ij}=a_i^tb_j^{1-t}E_{ij}=L_{A^t}R_{B^{1-t}}E_{ij}.

したがって RB#tLA=LAtRB1−tR_B\#_tL_A=L_{A^t}R_{B^{1-t}} で、

tr(K∗AtKB1−t)=⟨K, LAtRB1−tK⟩=⟨K, (RB#tLA)K⟩.\mathrm{tr}\bigl(K^*A^tKB^{1-t}\bigr)=\bigl\langle K,\,L_{A^t}R_{B^{1-t}}K\bigr\rangle=\bigl\langle K,\,(R_B\#_tL_A)K\bigr\rangle .

(A,B)↦(RB,LA)(A,B)\mapsto(R_B,L_A) は線形で、正定値行列の組を MnM_n 上の正定値作用素の組に写すので、定理 9.5(Hilbert 空間 MnM_n に適用)から右辺は (A,B)(A,B) について同時凹です。

半正定値の場合は A+εA+\varepsilon, B+εB+\varepsilon に適用し ε↓0\varepsilon\downarrow0 とします。(λ+ε)t→λt(\lambda+\varepsilon)^t\to\lambda^t は [0,∥A∥][0,\|A\|] 上一様なので、関数計算のノルムの等式から (A+ε)t→At(A+\varepsilon)^t\to A^t で、凹性の不等式は極限で保たれます。

系 9.7(相対エントロピーの同時凸性)✓

正定値行列 ρ,σ\rho,\sigma について S(ρ∥σ)=tr(ρlog⁡ρ−ρlog⁡σ)S(\rho\|\sigma)=\mathrm{tr}\bigl(\rho\log\rho-\rho\log\sigma\bigr) と置くと、(ρ,σ)↦S(ρ∥σ)(\rho,\sigma)\mapsto S(\rho\|\sigma) は同時凸である。とくに密度行列 ρi,σi>0\rho_i,\sigma_i>0 と確率 pip_i について S(∑ipiρi∥∑ipiσi)≤∑ipiS(ρi∥σi)S\bigl(\sum_ip_i\rho_i\big\|\sum_ip_i\sigma_i\bigr)\le\sum_ip_iS(\rho_i\|\sigma_i)。

証明(系 9.7)

ρ=∑iriuiui∗\rho=\sum_ir_iu_iu_i^*, σ=∑jqjwjwj∗\sigma=\sum_jq_jw_jw_j^* とし、cij=∣⟨ui,wj⟩∣2c_{ij}=|\langle u_i,w_j\rangle|^2(∑jcij=1\sum_jc_{ij}=1)と置くと

g(t)=tr(ρtσ1−t)=∑i,jcijritqj1−tg(t)=\mathrm{tr}\bigl(\rho^t\sigma^{1-t}\bigr)=\sum_{i,j}c_{ij}r_i^tq_j^{1-t}

は tt について微分可能で、g(1)=∑iri=trρg(1)=\sum_ir_i=\mathrm{tr}\rho、

g′(1)=∑i,jcijri(log⁡ri−log⁡qj)=tr(ρlog⁡ρ)−tr(ρlog⁡σ)=S(ρ∥σ)g'(1)=\sum_{i,j}c_{ij}r_i(\log r_i-\log q_j)=\mathrm{tr}(\rho\log\rho)-\mathrm{tr}(\rho\log\sigma)=S(\rho\|\sigma)

です(tr(ρlog⁡σ)=∑i,jrilog⁡qj cij\mathrm{tr}(\rho\log\sigma)=\sum_{i,j}r_i\log q_j\,c_{ij})。よって

S(ρ∥σ)=lim⁡t↑1trρ−tr(ρtσ1−t)1−t.S(\rho\|\sigma)=\lim_{t\uparrow1}\frac{\mathrm{tr}\rho-\mathrm{tr}\bigl(\rho^t\sigma^{1-t}\bigr)}{1-t}.

各 t∈(0,1)t\in(0,1) で右辺の分子は、線形関数 trρ\mathrm{tr}\rho から 定理 9.6(K=1K=1)の同時凹関数を引いたものなので同時凸で、正の数 1−t1-t で割っても同時凸です。同時凸関数の各点極限は同時凸です。

Klein の不等式(定理 8.1)が与えたのは S(ρ∥σ)≥0S(\rho\|\sigma)\ge0 でした。同時凸性はそれよりずっと強く、部分トレースのもとでの単調性 S(ρ12∥σ12)≥S(ρ1∥σ1)S(\rho_{12}\|\sigma_{12})\ge S(\rho_1\|\sigma_1) や、量子エントロピーの強劣加法性はここから導かれます(ユニタリ群上の平均で部分トレースを表し、同時凸性とユニタリ不変性を使います)。

注意 9.8

この記事で証明せずに残したのは、p<2p\lt2 の Birman–Solomyak の定理(注意 6.9)と Araki–Lieb–Thirring の不等式の一般形です。Kato–Seiler–Simon の不等式(定理 6.6)は、Schatten クラスの補間理論の全体ではなく、p=2p=2 と p=∞p=\infty をつなぐのに必要な形の三線補題(補題 6.5)だけで証明しました。前者はさらに格子上の分解を要するので主張だけにとどめています。

演習 10.1易

L2(0,∞)L^2(0,\infty) 上の Hardy 作用素 (Hf)(x)=1x∫0xf(y) dy(Hf)(x)=\frac1x\int_0^xf(y)\,dy について、定理 5.1 を p(x)=q(x)=x−1/2p(x)=q(x)=x^{-1/2} で使い、∥H∥≤2\|H\|\le2 を示してください。

解答

核は K(x,y)=x−11{y<x}K(x,y)=x^{-1}\mathbf{1}_{\{y<x\}} です。∫0xx−1y−1/2dy=2x−1/2\int_0^xx^{-1}y^{-1/2}dy=2x^{-1/2}、∫y∞x−1x−1/2dx=2y−1/2\int_y^\infty x^{-1}x^{-1/2}dx=2y^{-1/2} なので C1=C2=2C_1=C_2=2 で、∥H∥≤2⋅2=2\|H\|\le\sqrt{2\cdot2}=2。これが Hardy の不等式 ∫0∞∣1x∫0xf∣2dx≤4∫0∞∣f∣2\int_0^\infty\bigl|\frac1x\int_0^xf\bigr|^2dx\le4\int_0^\infty|f|^2 で、定数 44 は最良です。

演習 10.2標準

A,B≥0A,B\ge0 を有界作用素、X∈B(H)X\in\mathcal{B}(\mathcal{H})、0≤α≤10\le\alpha\le1 とします。定理 4.1 の証明と同じ方法で ∥AαXBα∥≤∥AXB∥α∥X∥1−α\|A^\alpha XB^\alpha\|\le\|AXB\|^\alpha\|X\|^{1-\alpha} を示し、X=1X=1 として 定理 4.2 を再び導いてください。

解答

A,BA,B を正定値として示せば、A+εA+\varepsilon, B+εB+\varepsilon で近似して一般の場合が出ます。F(z)=⟨u,AzXBzv⟩F(z)=\langle u,A^zXB^zv\rangle は帯 0≤Re z≤10\le\mathrm{Re}\,z\le1 で有界連続、内部で正則です。AiyA^{iy}, BiyB^{iy} はユニタリなので ∣F(iy)∣≤∥X∥∥u∥∥v∥|F(iy)|\le\|X\|\|u\|\|v\|、∣F(1+iy)∣=∣⟨u,Aiy(AXB)Biyv⟩∣≤∥AXB∥∥u∥∥v∥|F(1+iy)|=|\langle u,A^{iy}(AXB)B^{iy}v\rangle|\le\|AXB\|\|u\|\|v\|。三線定理から ∣F(α)∣≤∥X∥1−α∥AXB∥α∥u∥∥v∥|F(\alpha)|\le\|X\|^{1-\alpha}\|AXB\|^\alpha\|u\|\|v\| です。X=1X=1 なら ∥AαBα∥≤∥AB∥α\|A^\alpha B^\alpha\|\le\|AB\|^\alpha です。

演習 10.3標準

Hermite 行列 A,BA,B について φ(t)=log⁡tr eA+tB\varphi(t)=\log\mathrm{tr}\,e^{A+tB} が凸関数であることを、定理 7.3 と 系 8.3 を使って示してください。

解答

定理 7.3 から φ′(t)=tr(eA+tBB)/tr eA+tB\varphi'(t)=\mathrm{tr}(e^{A+tB}B)/\mathrm{tr}\,e^{A+tB} です。系 8.3 を A+tBA+tB と (s−t)B(s-t)B に使うと

φ(s)≥φ(t)+(s−t)tr(eA+tBB)tr eA+tB=φ(t)+(s−t)φ′(t)\varphi(s)\ge\varphi(t)+(s-t)\frac{\mathrm{tr}(e^{A+tB}B)}{\mathrm{tr}\,e^{A+tB}}=\varphi(t)+(s-t)\varphi'(t)

がすべての s,ts,t で成り立ちます。微分可能な関数が各点で接線の上にあれば凸なので、φ\varphi は凸です。統計力学では、自由エネルギーが結合定数について凹(−φ-\varphi)であることに対応します。

演習 10.4標準

Hermite 行列 A,BA,B について

tr eA⋅exp⁡(tr(eAB)tr eA)≤tr eA+B≤tr(eAeB)\mathrm{tr}\,e^A\cdot\exp\Bigl(\frac{\mathrm{tr}(e^AB)}{\mathrm{tr}\,e^A}\Bigr)\le\mathrm{tr}\,e^{A+B}\le\mathrm{tr}\bigl(e^Ae^B\bigr)

を示し、A,BA,B が可換で B=b⋅1B=b\cdot1(b∈Rb\in\mathbb{R})のとき両側が等号になることを確かめてください。また、2×22\times2 の例 A=(100−1)A=\begin{pmatrix}1&0\\0&-1\end{pmatrix}, B=(0110)B=\begin{pmatrix}0&1\\1&0\end{pmatrix} で 3 つの量を計算し、不等式が狭義であることを確かめてください。

解答

左の不等式は 系 8.3、右は 定理 8.6 です。B=b⋅1B=b\cdot1 なら 3 つの量はどれも eb tr eAe^b\,\mathrm{tr}\,e^A です。

例では (A+B)2=2⋅1(A+B)^2=2\cdot1 なので eA+B=cosh⁡2⋅1+sinh⁡22(A+B)e^{A+B}=\cosh\sqrt2\cdot1+\frac{\sinh\sqrt2}{\sqrt2}(A+B)、tr eA+B=2cosh⁡2≈4.356\mathrm{tr}\,e^{A+B}=2\cosh\sqrt2\approx4.356。eA=diag(e,e−1)e^A=\mathrm{diag}(e,e^{-1}) で tr(eAB)=0\mathrm{tr}(e^AB)=0 なので左辺は tr eA=2cosh⁡1≈3.086\mathrm{tr}\,e^A=2\cosh1\approx3.086。eB=cosh⁡1⋅1+sinh⁡1⋅Be^B=\cosh1\cdot1+\sinh1\cdot B なので tr(eAeB)=cosh⁡1⋅tr eA=2cosh⁡21≈4.762\mathrm{tr}(e^Ae^B)=\cosh1\cdot\mathrm{tr}\,e^A=2\cosh^21\approx4.762。よって 3.086<4.356<4.7623.086<4.356<4.762 で、どちらも狭義です。

演習 10.5標準

系 9.7 を使って、密度行列 ρ,σ>0\rho,\sigma>0 と、ユニタリ行列 U1,…,UNU_1,\ldots,U_N、確率 p1,…,pNp_1,\ldots,p_N について S(∑kpkUkρUk∗∥∑kpkUkσUk∗)≤S(ρ∥σ)S\bigl(\sum_kp_kU_k\rho U_k^*\big\|\sum_kp_kU_k\sigma U_k^*\bigr)\le S(\rho\|\sigma) を示してください。

解答

log⁡(UσU∗)=U(log⁡σ)U∗\log(U\sigma U^*)=U(\log\sigma)U^* と tr\mathrm{tr} の巡回性から S(UρU∗∥UσU∗)=S(ρ∥σ)S(U\rho U^*\|U\sigma U^*)=S(\rho\|\sigma) です。同時凸性から

S(∑kpkUkρUk∗∥∑kpkUkσUk∗)≤∑kpkS(UkρUk∗∥UkσUk∗)=S(ρ∥σ).S\Bigl(\sum_kp_kU_k\rho U_k^*\Big\|\sum_kp_kU_k\sigma U_k^*\Bigr)\le\sum_kp_kS(U_k\rho U_k^*\|U_k\sigma U_k^*)=S(\rho\|\sigma).

適切なユニタリの平均で部分トレース(と最大混合状態のテンソル積)を表せることと組み合わせると、部分トレースのもとでの単調性が得られます。

  • R. Bhatia, Matrix Analysis, Springer (GTM 169), 1997 — Chapter V(作用素単調関数、Löwner–Heinz)、Chapter IX(Golden–Thompson、トレース不等式)。
  • T. Kato, Perturbation Theory for Linear Operators, 2nd edition, Springer, 1976 — Chapter V §4(Heinz–Kato の不等式)。
  • B. Simon, Trace Ideals and Their Applications, 2nd edition, American Mathematical Society, 2005 — Chapter 2(Schatten クラスの変分表示と Hölder 不等式)、Chapter 4(Kato–Seiler–Simon、Birman–Solomyak)。
  • M. Cwikel, “Weak type estimates for singular values and the number of bound states of Schrödinger operators”, Annals of Mathematics 106 (1977) 93–100。
  • M. Sh. Birman, M. Solomyak, “Double operator integrals in a Hilbert space”, Integral Equations Operator Theory 47 (2003) 131–168。
  • T. Ando, “Concavity of certain maps on positive definite matrices and applications to Hadamard products”, Linear Algebra and its Applications 26 (1979) 203–241。
  • E. A. Carlen, “Trace inequalities and quantum entropy: an introductory course”, in Entropy and the Quantum, Contemporary Mathematics 529, American Mathematical Society, 2010, 73–140。
  • G. H. Hardy, J. E. Littlewood, G. Pólya, Inequalities, 2nd edition, Cambridge University Press, 1952 — Chapter IX(Hilbert と Hardy の不等式)。

このテーマの 100 本ノック(100 問)

執筆: Claude Code(Opus 5.5)、Claude Code(Sonnet 5.5)・監査 3 回(Claude Code)・最終 2026-10-12

この記事の誤りを報告する ・運営: 夢現技研合同会社 ・料金 ・利用条件 ・特定商取引法に基づく表記

© 2026 夢現技研合同会社 ・コード例は MIT ライセンスです。