平成21年度 東大院 物理学専攻 修士 数学 解答

出典と、この解答の位置づけ

本記事が扱う問題は、東京大学大学院理学系研究科物理学専攻の大学院入試問題からの引用です（引用元: https://www.phys.s.u-tokyo.ac.jp/）。

• 掲載している解答・解説は東京大学が公表したものではなく、夢現技研合同会社が独自に作成したものです。誤りが含まれる可能性があります。
• 本記事に問題文は転載していません。問題は上記の専攻公式サイトでご確認ください。
• 本記事の内容について、東京大学および同専攻へのお問い合わせはご遠慮ください。誤りのご指摘はお問い合わせまでお願いします。

大問2題で、2問すべてが必答です（試験時間は60分）。第1問は前半が $\theta$ に依存する2次実対称行列で、固有値が $\cos\theta$ と $\sin\theta$ という簡単な形になることを見抜けば、以降は Cayley–Hamilton の定理と対角化を使い回すだけで最後まで到達します。後半は固有値が非負な実対称行列に対するトレース不等式で、$\lambda_i \le \operatorname{Tr}X$ という一行の観察が証明と等号条件の両方を同時に片付けます。第2問は2準位系のシュレーディンガー方程式と同じ形の連立1階常微分方程式で、係数行列がエルミートであることによるノルム保存、対角成分を位相に繰り込む変換、そして最後にローレンツ型の積分を留数で評価する構成です。

第1問(iii) で行列多項式が $A-I$ に潰れることに気づくか、第2問設問4で係数行列が時間に依らない行列に比例するため時間順序積が単純な指数関数に落ちることに気づくか、この2点が分かれ目になります。

問題	分野	主題
第1問	線形代数	2次実対称行列の固有値と行列多項式、半正定値行列のトレース不等式
第2問	微分方程式・複素解析	2成分連立1階常微分方程式と遷移確率

第1問 2次実対称行列の固有値と行列多項式、トレース不等式

設定

設問1では、実変数 $\theta$（範囲は $0\le\theta\le\pi/2$）に依存する2次実対称行列

$$A=\frac{1}{4}\begin{pmatrix}\cos\theta+3\sin\theta & -\sqrt{3}\cos\theta+\sqrt{3}\sin\theta\\ -\sqrt{3}\cos\theta+\sqrt{3}\sin\theta & 3\cos\theta+\sin\theta\end{pmatrix}$$

を扱います。固有値を求め、$f(\theta)=(\operatorname{Tr}A)^3-\operatorname{Tr}(A^3)$ の最大最小を出し、行列多項式

$$B=A^4-A^2+A+(\cos^2\theta\sin^2\theta-1)I$$

が特異になる $\theta$ を決め、最後に $B^{-1}$ を $A$ の1次式として書きます。$I$ は2次単位行列です。

設問2では、$N$ 次実対称行列 $X$ の固有値 $\lambda_i$ $(i=1,\dots,N)$ がすべて非負であるとして、$(\operatorname{Tr}X)^n\ge\operatorname{Tr}(X^n)$ とその等号条件を調べます。

以下、第1問を通して $c=\cos\theta$、$s=\sin\theta$ と略記します。範囲 $0\le\theta\le\pi/2$ では $c\ge0$、$s\ge0$ です。

設問1

(i) 2次行列の固有値はトレースと行列式で決まるので、その2つを計算します。

$$\operatorname{Tr}A=\frac{1}{4}\bigl\{(c+3s)+(3c+s)\bigr\}=\frac{4c+4s}{4}=c+s .$$

行列式は、非対角成分が $\frac{\sqrt3}{4}(s-c)$ であることに注意して

$$\begin{aligned} 16\det A&=(c+3s)(3c+s)-3(s-c)^2\\ &=\bigl(3c^2+10cs+3s^2\bigr)-3\bigl(1-2cs\bigr)\\ &=(3+10cs)-(3-6cs)=16cs \end{aligned}$$

となり、$\det A=cs$ です。ここで $c^2+s^2=1$ を使いました。したがって固有多項式は

$$\lambda^2-(c+s)\lambda+cs=(\lambda-c)(\lambda-s)$$

と因数分解でき、答えは $\lambda=\cos\theta$ と $\lambda=\sin\theta$ の2つです。

固有ベクトルが $\theta$ に依存しないことも確認しておきます。

$$P=\frac{1}{4}\begin{pmatrix}1&-\sqrt3\\-\sqrt3&3\end{pmatrix},\qquad Q=\frac{1}{4}\begin{pmatrix}3&\sqrt3\\\sqrt3&1\end{pmatrix}$$

と置くと $A=cP+sQ$ で、直接計算すると $P^2=P$、$Q^2=Q$、$PQ=QP=O$、$P+Q=I$ が成り立ちます。すなわち $P,Q$ は互いに直交する1次元射影で、$P$ は単位ベクトル $u=\frac12(1,-\sqrt3)^{T}$、$Q$ は $w=\frac12(\sqrt3,1)^{T}$ への直交射影です（$u\cdot w=0$、$\lvert u\rvert=\lvert w\rvert=1$）。よって

$$Au=\cos\theta\,u,\qquad Aw=\sin\theta\,w$$

であり、固有ベクトルは $\theta$ によらず $u,w$ に固定されています。これは以降の設問でそのまま使えます。

(ii) $A$ の固有値が $c,s$ なので $A^3$ の固有値は $c^3,s^3$ で、$\operatorname{Tr}(A^3)=c^3+s^3$ です。したがって

$$f(\theta)=(c+s)^3-(c^3+s^3)=3c^2s+3cs^2=3cs(c+s).$$

$v=c+s$ を独立変数にとります。$c^2+s^2=1$ より $2cs=v^2-1$ なので

$$f=\frac{3}{2}\,v\,(v^2-1)=\frac{3}{2}\bigl(v^3-v\bigr)\equiv g(v).$$

$v=c+s=\sqrt2\sin\!\left(\theta+\frac{\pi}{4}\right)$ で、$0\le\theta\le\pi/2$ のとき $\theta+\frac{\pi}{4}\in[\frac{\pi}{4},\frac{3\pi}{4}]$ ですから $\sin(\theta+\frac{\pi}{4})\in[\frac{1}{\sqrt2},1]$、すなわち $v$ は $[1,\sqrt2]$ を動きます。$v$ は $\theta$ の連続関数なので $[1,\sqrt2]$ の値をすべてとります。

この区間で $g'(v)=\frac32(3v^2-1)\ge\frac32\cdot 2>0$ なので $g$ は狭義単調増加です。よって最大値と最小値は端点で実現され、

$$\max f=g(\sqrt2)=\frac{3}{2}\bigl(2\sqrt2-\sqrt2\bigr)=\frac{3\sqrt2}{2},\qquad \min f=g(1)=0 .$$

$v=\sqrt2$ は $\theta=\pi/4$、$v=1$ は $\theta=0$ または $\theta=\pi/2$ に対応します。答えは、最大値が $\dfrac{3\sqrt2}{2}=\dfrac{3}{\sqrt2}\simeq2.121$（$\theta=\pi/4$ のとき）、最小値が $0$（$\theta=0,\ \pi/2$ のとき）です。

検算として $f=3cs(c+s)$ の形を直接見ると、$0\le\theta\le\pi/2$ では $c,s\ge0$ なので $f\ge0$ で、$f=0$ となるのは $cs=0$ すなわち $\theta=0,\pi/2$ に限られます。$\theta=\pi/4$ では $c=s=1/\sqrt2$ から $f=3\cdot\frac12\cdot\sqrt2=\frac{3\sqrt2}{2}$ で、上の結果と一致します。

(iii) $B$ は $A$ の多項式 $B=p(A)$ で、

$$p(x)=x^4-x^2+x+\bigl(c^2s^2-1\bigr)$$

です。$A$ は実対称なので直交行列 $R$ で $A=R\,\mathrm{diag}(c,s)\,R^{T}$ と対角化でき、多項式は対角成分ごとに作用して

$$B=p(A)=R\,\mathrm{diag}\bigl(p(c),p(s)\bigr)R^{T}$$

となります。ここで $s^2=1-c^2$ より $c^2s^2=c^2-c^4$ なので

$$p(c)=c^4-c^2+c+\bigl(c^2-c^4\bigr)-1=c-1,$$

同様に $c^2s^2=s^2-s^4$ を使って $p(s)=s-1$ です。したがって

$$B=R\,\mathrm{diag}(c-1,\;s-1)\,R^{T}=A-I$$

であり、$B$ は $\cos^2\theta\sin^2\theta$ という項のおかげで4次式が1次式に潰れています。固有値は $\cos\theta-1$ と $\sin\theta-1$ で、

$$\det B=(\cos\theta-1)(\sin\theta-1).$$

$B$ が逆行列を持たないのは $\det B=0$、すなわち $\cos\theta=1$ または $\sin\theta=1$ のときです。$0\le\theta\le\pi/2$ の範囲では前者は $\theta=0$、後者は $\theta=\pi/2$ のみです。答えは $\theta=0$ と $\theta=\dfrac{\pi}{2}$ です。

念のため $\theta=0$ を直接見ると $A=P$ で、$P^2=P$ と $c^2s^2=0$ から $B=P-P+P-I=P-I$ です。$P$ の固有値は $1,0$ なので $P-I$ の固有値は $0,-1$ となり、確かに特異です。

(iv) 前問より $\theta\ne0,\pi/2$ のとき $B=A-I$ は可逆です。Cayley–Hamilton の定理から

$$A^2-(c+s)A+cs\,I=O,\qquad \text{すなわち}\quad A^2=(c+s)A-cs\,I$$

が使えます。$X=a_1A+a_0I$ と置いて $(A-I)X=I$ を要求すると

$$\begin{aligned} (A-I)(a_1A+a_0I)&=a_1A^2+(a_0-a_1)A-a_0I\\ &=a_1\bigl\{(c+s)A-cs\,I\bigr\}+(a_0-a_1)A-a_0I\\ &=\bigl\{a_1(c+s-1)+a_0\bigr\}A-\bigl(a_1cs+a_0\bigr)I . \end{aligned}$$

これが $I$ に等しくなるためには

$$\begin{cases} a_1(c+s-1)+a_0=0,\\ a_1cs+a_0=-1 \end{cases}$$

が成り立てば十分です（$A$ と $I$ が1次独立かどうかを問わず、この2式から行列としての等式が従います）。上式から下式を引くと

$$a_1\bigl(c+s-1-cs\bigr)=1 .$$

ここで

$$c+s-1-cs=-\bigl(1-c-s+cs\bigr)=-(1-c)(1-s)$$

であり、$\theta\ne0,\pi/2$ では $c\ne1$ かつ $s\ne1$ なのでこれは $0$ になりません。よって

$$a_1(\theta)=-\frac{1}{(1-\cos\theta)(1-\sin\theta)},\qquad a_0(\theta)=-a_1(c+s-1)=\frac{\cos\theta+\sin\theta-1}{(1-\cos\theta)(1-\sin\theta)} .$$

まとめると答えは

$$B^{-1}=\frac{-A+(\cos\theta+\sin\theta-1)I}{(1-\cos\theta)(1-\sin\theta)} .$$

固有値で検算します。$A$ の固有ベクトル $u,w$ 上で

$$a_1c+a_0=\frac{-c+(c+s-1)}{(1-c)(1-s)}=\frac{-(1-s)}{(1-c)(1-s)}=\frac{1}{c-1},$$

同様に $a_1s+a_0=\dfrac{1}{s-1}$ となり、これは $B=A-I$ の固有値 $c-1,\ s-1$ の逆数です。$B$ と $a_1A+a_0I$ は同じ直交基底で対角化されるので、両者は互いに逆行列です。

なお $\theta=\pi/4$ では $A=\frac{1}{\sqrt2}I$ となって $I$ と $A$ が1次従属になるため、$(a_1,a_0)$ の組は一意ではありません。上の式はそのうちの1つを与えており、実際に代入すると

$$a_1\cdot\frac{1}{\sqrt2}+a_0=\frac{-\frac{1}{\sqrt2}+\sqrt2-1}{\left(1-\frac{1}{\sqrt2}\right)^2} =\frac{-\left(1-\frac{1}{\sqrt2}\right)}{\left(1-\frac{1}{\sqrt2}\right)^2}=\frac{1}{\frac{1}{\sqrt2}-1}$$

で、$B^{-1}=\left(\frac{1}{\sqrt2}-1\right)^{-1}I$ に一致します。

設問2

(i) $X$ は実対称なので直交行列で対角化でき、固有値 $\lambda_i$ は実数です。$X=R\,\mathrm{diag}(\lambda_1,\dots,\lambda_N)R^{T}$ とすると $X^n=R\,\mathrm{diag}(\lambda_1^n,\dots,\lambda_N^n)R^{T}$ なので

$$\operatorname{Tr}X=\sum_{i=1}^{N}\lambda_i\equiv S,\qquad \operatorname{Tr}(X^n)=\sum_{i=1}^{N}\lambda_i^{\,n}$$

です。仮定 $\lambda_i\ge0$ より、各 $i$ について

$$0\le\lambda_i\le\sum_{j=1}^{N}\lambda_j=S$$

が成り立ちます（他の固有値をすべて足しても減らない）。とくに $S\ge0$ です。

$n=1$ のときは両辺とも $S$ で等号が成立し、主張は成り立ちます。$n\ge2$ のときは、$x\mapsto x^{n-1}$ が $[0,\infty)$ 上で単調増加なので $0\le\lambda_i\le S$ から

$$\lambda_i^{\,n-1}\le S^{\,n-1}.$$

両辺に $\lambda_i\ge0$ をかけて（不等号の向きは保たれます）

$$\lambda_i^{\,n}\le\lambda_i\,S^{\,n-1}\qquad(i=1,\dots,N)$$

を得ます。$i$ について総和をとると

$$\operatorname{Tr}(X^n)=\sum_{i=1}^{N}\lambda_i^{\,n}\le S^{\,n-1}\sum_{i=1}^{N}\lambda_i=S^{\,n}=(\operatorname{Tr}X)^n$$

となり、任意の自然数 $n$ について $(\operatorname{Tr}X)^n\ge\operatorname{Tr}(X^n)$ が示されました。証明に使った仮定は、$X$ が実対称で対角化可能であること（$\operatorname{Tr}(X^n)=\sum\lambda_i^n$ が言えること）と、全固有値が非負であることの2つだけです。

別の見方として、多項定理で

$$S^{\,n}=\Bigl(\sum_j\lambda_j\Bigr)^{\!n}=\sum_{k_1+\cdots+k_N=n}\frac{n!}{k_1!\cdots k_N!}\,\lambda_1^{k_1}\cdots\lambda_N^{k_N}$$

と展開すると、$\lambda_i\ge0$ より全項が非負で、そのうち $k_i=n$（他は $0$）の $N$ 項がちょうど $\sum_i\lambda_i^n$ を与えます。残りの項が非負であることが不等式の中身です。

(ii) $n\ge2$ とします。(i) の証明では各 $i$ ごとの不等式 $\lambda_i^{\,n}\le\lambda_iS^{\,n-1}$ を足し合わせただけなので、和について等号が成立することは、すべての $i$ で個別に等号が成立すること

$$\lambda_i\bigl(\lambda_i^{\,n-1}-S^{\,n-1}\bigr)=0\qquad(i=1,\dots,N)$$

と同値です。$\lambda_i\ge0$、$S\ge0$ で、$n\ge2$ のとき $x\mapsto x^{n-1}$ は $[0,\infty)$ 上で単射なので、この条件は「各 $i$ について $\lambda_i=0$ または $\lambda_i=S$」と書けます。

ここで $\lambda_k>0$ となる添字が2つ以上あったとします。それらを $k\ne l$ とすると $\lambda_k=\lambda_l=S>0$ であり、

$$S=\sum_i\lambda_i\ge\lambda_k+\lambda_l=2S$$

から $S\le0$、すなわち $S=0$ となって $\lambda_k>0$ に矛盾します。よって正の固有値は高々1個です。

逆に、固有値のうち高々1個だけが正で残りが $0$ のとき、その値を $\lambda\ (\ge0)$ とすれば $(\operatorname{Tr}X)^n=\lambda^n=\operatorname{Tr}(X^n)$ で等号が成立します（すべてが $0$、つまり $X=O$ の場合は両辺 $0$ で成立）。

したがって答えは、$n\ge2$ での等号成立は「$\lambda_i$ のうち非零のものが高々1個、すなわち少なくとも $N-1$ 個の固有値が $0$ である場合」に限られます。行列の言葉では $\operatorname{rank}X\le1$、つまり単位ベクトル $v$ と $\lambda\ge0$ を用いて $X=\lambda\,vv^{T}$ と書ける場合（$X=O$ を含む）です。$n=1$ を除外する理由もこれで明らかで、$n=1$ なら固有値によらず常に等号が成立してしまいます。

多項定理による見方でも同じ結論が出ます。等号は混合項がすべて消えること、とくに $k_i=1,\ k_j=n-1$ の項 $n\lambda_i\lambda_j^{\,n-1}=0$ が任意の $i\ne j$ で成り立つことを要求し、これは $\lambda_i\lambda_j=0$ $(i\ne j)$、すなわち非零固有値が高々1個であることと同じです。

第2問 2成分連立1階常微分方程式と遷移確率

設定

実変数 $t$ の複素関数 $f_1(t), f_2(t)$ が

$$i\frac{d}{dt}\begin{pmatrix}f_1(t)\\ f_2(t)\end{pmatrix} =\begin{pmatrix}a(t)&b(t)\\ b(t)&c(t)\end{pmatrix} \begin{pmatrix}f_1(t)\\ f_2(t)\end{pmatrix}$$

を満たします。$a(t),b(t),c(t)$ は $t$ の実関数（以下では解が微分可能であるように連続とします）、$i=\sqrt{-1}$ です。係数行列は実対称、したがってエルミートなので、これは2準位系のシュレーディンガー方程式と同じ構造をしています。

以下では

$$f(t)=\begin{pmatrix}f_1(t)\\ f_2(t)\end{pmatrix},\qquad H(t)=\begin{pmatrix}a(t)&b(t)\\ b(t)&c(t)\end{pmatrix}$$

と書き、上の方程式を $i\,df/dt=H(t)f$ と表します。設問3以降は $a=c=0$、すなわち $H(t)=b(t)\sigma_x$（$\sigma_x$ は非対角成分が $1$ の実対称行列）の場合です。

設問1

成分で書くと

$$\frac{df_1}{dt}=-i\bigl(a f_1+b f_2\bigr),\qquad \frac{df_2}{dt}=-i\bigl(b f_1+c f_2\bigr)$$

です。$N(t)=\lvert f_1\rvert^2+\lvert f_2\rvert^2=f_1^{*}f_1+f_2^{*}f_2$ を微分すると

$$\frac{dN}{dt}=\Bigl(\dot f_1^{*}f_1+f_1^{*}\dot f_1\Bigr)+\Bigl(\dot f_2^{*}f_2+f_2^{*}\dot f_2\Bigr) =2\,\mathrm{Re}\Bigl(f_1^{*}\dot f_1+f_2^{*}\dot f_2\Bigr)$$

となります。上の方程式を代入すると

$$\begin{aligned} f_1^{*}\dot f_1+f_2^{*}\dot f_2 &=-i\Bigl(a\lvert f_1\rvert^2+b f_1^{*}f_2+b f_2^{*}f_1+c\lvert f_2\rvert^2\Bigr)\\ &=-i\Bigl(a\lvert f_1\rvert^2+c\lvert f_2\rvert^2+2b\,\mathrm{Re}\bigl(f_1^{*}f_2\bigr)\Bigr) \end{aligned}$$

です。2行目では $f_1^{*}f_2+f_2^{*}f_1=f_1^{*}f_2+(f_1^{*}f_2)^{*}=2\,\mathrm{Re}(f_1^{*}f_2)$ を使いました。丸括弧の中身は、$a,b,c$ が実関数であることから実数です。よって $f_1^{*}\dot f_1+f_2^{*}\dot f_2$ は純虚数で、その実部は $0$、すなわち

$$\frac{dN}{dt}=0\qquad(\text{すべての }t)$$

です。$N$ は微分が恒等的に $0$ の実関数なので定数であり、$\lvert f_1(t)\rvert^2+\lvert f_2(t)\rvert^2$ は $t$ に依存しません。

証明で本質的に使ったのは、$a,b,c$ が実で係数行列の非対角成分が等しいこと、つまり $H(t)^{\dagger}=H(t)$ であることです。行列形式で書けば $\dot f=-iHf$、$\dot f^{\dagger}=i f^{\dagger}H^{\dagger}=if^{\dagger}H$ から

$$\frac{d}{dt}\bigl(f^{\dagger}f\bigr)=i f^{\dagger}Hf-i f^{\dagger}Hf=0$$

となり、同じ結論が1行で出ます。$H$ がエルミートでなければ成り立ちません。

設問2

$$\alpha(t)=\int_0^{t}a(\tau)\,d\tau,\qquad \gamma(t)=\int_0^{t}c(\tau)\,d\tau$$

と置きます。$a,c$ は実で連続なので $\alpha,\gamma$ は実の $C^1$ 関数で、$\dot\alpha=a$、$\dot\gamma=c$、$\alpha(0)=\gamma(0)=0$ です。定義式は $f_1=e^{-i\alpha}\tilde f_1$、$f_2=e^{-i\gamma}\tilde f_2$ ですが、$\lvert e^{-i\alpha}\rvert=1$ なので $\tilde f_1=e^{i\alpha}f_1$、$\tilde f_2=e^{i\gamma}f_2$ と逆に解けます。すなわちこれは可逆な変数変換で、位相を付け替えているだけです。

第1成分の式に代入します。積の微分から

$$i\dot f_1=i\Bigl(-i\dot\alpha\,e^{-i\alpha}\tilde f_1+e^{-i\alpha}\dot{\tilde f}_1\Bigr) =a\,e^{-i\alpha}\tilde f_1+i\,e^{-i\alpha}\dot{\tilde f}_1 =a f_1+i\,e^{-i\alpha}\dot{\tilde f}_1$$

であり、これがもとの方程式の右辺 $af_1+bf_2=af_1+b\,e^{-i\gamma}\tilde f_2$ に等しいので、$af_1$ が両辺で消えて

$$i\,e^{-i\alpha}\dot{\tilde f}_1=b\,e^{-i\gamma}\tilde f_2, \qquad\text{すなわち}\qquad i\dot{\tilde f}_1=b\,e^{\,i(\alpha-\gamma)}\,\tilde f_2 .$$

第2成分も同様で、$i\dot f_2=c f_2+i e^{-i\gamma}\dot{\tilde f}_2=b f_1+c f_2$ から

$$i\dot{\tilde f}_2=b\,e^{\,i(\gamma-\alpha)}\,\tilde f_1 .$$

そこで

$$\tilde b(t)=b(t)\,e^{\,i\{\alpha(t)-\gamma(t)\}} =b(t)\exp\left[\,i\int_0^{t}\bigl\{a(\tau)-c(\tau)\bigr\}d\tau\right]$$

と定義すると、$b$ と $\alpha-\gamma$ が実であることから $b\,e^{i(\gamma-\alpha)}=\tilde b^{*}$ となり、2式はまとめて

$$i\frac{d}{dt}\begin{pmatrix}\tilde f_1\\ \tilde f_2\end{pmatrix} =\begin{pmatrix}0&\tilde b(t)\\ \tilde b(t)^{*}&0\end{pmatrix} \begin{pmatrix}\tilde f_1\\ \tilde f_2\end{pmatrix}$$

という主張の形になります。求める表式は上の $\tilde b(t)$ です。

対角成分が消えたのは、$e^{-i\alpha}$、$e^{-i\gamma}$ が対角成分 $a$、$c$ の作る位相をちょうど打ち消すように選ばれているからです。新しい係数行列もエルミートで、$\lvert\tilde b\rvert=\lvert b\rvert$、$\tilde b(0)=b(0)$、また $\lvert\tilde f_j\rvert=\lvert f_j\rvert$ なので設問1のノルム保存もそのまま引き継がれます。$a=c$ の場合は $\tilde b=b$ で、位相以外は何も変わりません。

設問3

$a=c=0$、$b(t)=b_0$（実定数）のとき方程式は

$$i\dot f_1=b_0 f_2,\qquad i\dot f_2=b_0 f_1$$

です。和と差 $g_{\pm}=f_1\pm f_2$ をとると2式は分離して

$$i\dot g_{\pm}=\pm b_0\,g_{\pm}, \qquad\text{よって}\qquad g_{\pm}(t)=g_{\pm}(0)\,e^{\mp i b_0 t}$$

となります。$g_{\pm}(0)=f_1(0)\pm f_2(0)$ を使って戻すと

$$\begin{aligned} f_1(t)&=\frac{g_+(t)+g_-(t)}{2} =\frac{f_1(0)+f_2(0)}{2}e^{-ib_0t}+\frac{f_1(0)-f_2(0)}{2}e^{ib_0t}\\ &=f_1(0)\cos(b_0t)-i f_2(0)\sin(b_0t) . \end{aligned}$$

答えは $f_1(t)=f_1(0)\cos(b_0t)-i\,f_2(0)\sin(b_0t)$ です。同じ計算で

$$f_2(t)=-i f_1(0)\sin(b_0t)+f_2(0)\cos(b_0t)$$

も得られ、行列でまとめると

$$f(t)=U(t)f(0),\qquad U(t)=\cos(b_0t)\,I-i\sin(b_0t)\,\sigma_x =\begin{pmatrix}\cos b_0t & -i\sin b_0t\\ -i\sin b_0t & \cos b_0t\end{pmatrix}$$

です。検算として $t=0$ で $f_1=f_1(0)$ となること、また $\sigma_x^2=I$ と $b_0,\ t$ が実であることから

$$U^{\dagger}U=\bigl(\cos b_0t\,I+i\sin b_0t\,\sigma_x\bigr)\bigl(\cos b_0t\,I-i\sin b_0t\,\sigma_x\bigr) =\cos^2 b_0t\,I+\sin^2 b_0t\,\sigma_x^2=I$$

すなわち $U$ がユニタリであることを確認できます。これは設問1のノルム保存と整合します。$b_0=0$ なら $f_1(t)=f_1(0)$ で定数、これも正しい極限です。

設問4

$a=c=0$ で $b(t)$ が時間に依存する場合ですが、係数行列は $H(t)=b(t)\sigma_x$ と、時間に依らない行列 $\sigma_x$ のスカラー倍です。異なる時刻の $H$ が互いに可換なので、設問3と同じ分離がそのまま通ります。$g_{\pm}=f_1\pm f_2$ に対して

$$i\dot g_{\pm}=\pm b(t)\,g_{\pm} \qquad\Longrightarrow\qquad g_{\pm}(t)=g_{\pm}(t_1)\exp\left[\mp i\int_{t_1}^{t}b(\tau)\,d\tau\right] .$$

$b(t)$ は $t\to-\infty$ で十分はやく減衰するので $\int_{-\infty}^{t}\lvert b(\tau)\rvert d\tau<\infty$ であり、

$$\Phi(t)=\int_{-\infty}^{t}b(\tau)\,d\tau$$

が収束します。そこで $t_1\to-\infty$ の極限をとると（この極限で $g_\pm(t_1)\to g_\pm(-\infty)=f_1(-\infty)\pm f_2(-\infty)$ が存在します）

$$g_{\pm}(t)=\bigl\{f_1(-\infty)\pm f_2(-\infty)\bigr\}e^{\mp i\Phi(t)}$$

となり、$f_1=(g_++g_-)/2$ から

$$f_1(t)=f_1(-\infty)\cos\Phi(t)-i\,f_2(-\infty)\sin\Phi(t), \qquad \Phi(t)=\int_{-\infty}^{t}b(\tau)\,d\tau$$

が答えです。同様に $f_2=(g_+-g_-)/2$ から

$$f_2(t)=-i f_1(-\infty)\sin\Phi(t)+f_2(-\infty)\cos\Phi(t)$$

です。$b(t)=b_0$ を形式的に代入すると $\Phi$ は発散してしまいますが、$b$ が有限の区間の外で $0$ になる場合には、その区間より前で $f$ が定数、区間を通過した後は $\Phi$ が全積分値で止まる、という正しい振る舞いになります。$t\to-\infty$ で $\Phi\to0$ なので $f_1(t)\to f_1(-\infty)$ となり、定義とも整合します。$\Phi$ が実であることから $\lvert f_1\rvert^2+\lvert f_2\rvert^2$ が保存することも直接確かめられます。

設問5

$$b(t)=\frac{\beta\cos\omega t}{t^2+t_0^2},\qquad \beta,\omega,t_0>0$$

の場合です。$\lvert b(t)\rvert\le\beta/(t^2+t_0^2)$ で右辺は $\mathbb{R}$ 上可積分なので、設問4の仮定（$t\to-\infty$ での十分はやい減衰）が満たされ、さらに $t\to+\infty$ でも積分が収束します。したがって

$$\Theta\equiv\Phi(+\infty)=\int_{-\infty}^{\infty}b(t)\,dt =\beta\int_{-\infty}^{\infty}\frac{\cos\omega t}{t^2+t_0^2}\,dt$$

が有限に定まります。

この積分を留数で評価します。$\sin\omega t/(t^2+t_0^2)$ は奇関数で積分が $0$ なので

$$\int_{-\infty}^{\infty}\frac{\cos\omega t}{t^2+t_0^2}\,dt =\int_{-\infty}^{\infty}\frac{e^{\,i\omega t}}{t^2+t_0^2}\,dt$$

と書けます。積分路として実軸の $[-R,R]$ と上半平面の半円 $C_R$（半径 $R>t_0$）を合わせた閉曲線をとります。$\omega>0$ かつ上半平面では $\mathrm{Im}\,z\ge0$ なので $\lvert e^{i\omega z}\rvert=e^{-\omega\,\mathrm{Im}\,z}\le1$ で、また $\lvert z^2+t_0^2\rvert\ge R^2-t_0^2$ ですから

$$\left\lvert\int_{C_R}\frac{e^{\,i\omega z}}{z^2+t_0^2}dz\right\rvert \le\frac{\pi R}{R^2-t_0^2}\xrightarrow[R\to\infty]{}0$$

と半円の寄与は消えます。閉曲線の内部にある極は $z=it_0$（1位）だけで、留数は

$$\lim_{z\to it_0}\frac{e^{\,i\omega z}}{z+it_0}=\frac{e^{-\omega t_0}}{2it_0}$$

です。留数定理から

$$\int_{-\infty}^{\infty}\frac{\cos\omega t}{t^2+t_0^2}\,dt =2\pi i\cdot\frac{e^{-\omega t_0}}{2it_0}=\frac{\pi}{t_0}e^{-\omega t_0}, \qquad\text{よって}\qquad \Theta=\frac{\pi\beta}{t_0}e^{-\omega t_0}.$$

$f_1(-\infty)=1$、$f_2(-\infty)=0$ を設問4の結果に入れると $f_1(t)=\cos\Phi(t)$、$f_2(t)=-i\sin\Phi(t)$ なので、$t\to+\infty$ で $\Phi\to\Theta$ より

$$\lvert f_1(+\infty)\rvert^2=\cos^2\Theta=\cos^2\!\left(\frac{\pi\beta}{t_0}e^{-\omega t_0}\right), \qquad \lvert f_2(+\infty)\rvert^2=\sin^2\Theta=\sin^2\!\left(\frac{\pi\beta}{t_0}e^{-\omega t_0}\right)$$

が答えです。

確認をいくつかしておきます。方程式 $i\dot f=bf$ の形から $b$ は時間の逆数の次元をもち、$b=\beta\cos\omega t/(t^2+t_0^2)$ より $\beta$ は時間の次元です。よって $\beta/t_0$ と $\omega t_0$ はどちらも無次元で、三角関数と指数関数の引数として整合しています。また2つの答えの和は $\cos^2\Theta+\sin^2\Theta=1$ で、$f_1(-\infty)=1,\ f_2(-\infty)=0$ での $\lvert f_1\rvert^2+\lvert f_2\rvert^2=1$ と一致し、設問1のノルム保存を満たします。

極限も素直です。$\omega t_0\gg1$、すなわち振動 $\cos\omega t$ の周期がローレンツ型の包絡線の幅 $t_0$ よりずっと短いときは $\Theta\simeq0$ で、$\lvert f_2(+\infty)\rvert^2\simeq\Theta^2$ が指数関数的に小さくなります。速い振動が打ち消し合って遷移が起きない、という描像です。逆に $\omega\to0$ では $\Theta\to\pi\beta/t_0$ となり、パルス面積 $\int b\,dt$ そのものが遷移を決めます。$\Theta$ が $\pi$ の整数倍のときは $\lvert f_2(+\infty)\rvert^2=0$ で、系は始状態に完全に戻ります（ラビ振動と同じ構造です）。

出典: 東京大学大学院理学系研究科 物理学専攻 平成21年度 修士課程 入学試験問題 数学。問題文は要約して引用しています。