院試hub

東大 工学系研究科 物理工学専攻 物理学 2021年度 院試 過去問 解答例・解説

大学
東京大学
研究科・専攻
工学系研究科 物理工学専攻
科目
物理学
年度
2021年度

解説は全問、途中式と最終答は大問1問を無料で公開しています。PDFの収録範囲は商品ページで確認できます。問題本文は含みません。

東京大学 工学系研究科 物理工学専攻 物理学の解答PDF2012〜2026年度(15年度)

¥3,200

商品内容・見本・購入へ

第1問 量子ビットの時間発展とゲート生成

見通し

時間発展の本質は「エルミート演算子がユニタリーの生成子である」ことにある。[1][1] は形式的解、[2][2] はその帰結、[3][3] は確率保存(ノルム保存)を意味する。

指数の扱い

H^\hat H が時間に依存しないときのみ U^(Δt)=e−iH^Δt/ℏ\hat U(\Delta t)=e^{-i\hat H\Delta t/\hbar} と書ける。行列指数関数は対角化(固有値分解)で計算でき、U^2=I^\hat U^2=\hat I を満たす(エルミートな)ゲートは固有値が ±1\pm 1 なので、固有値 +1+1 を e0=1e^{0}=1、固有値 −1-1 を e−iπ=−1e^{-i\pi}=-1 に対応させる発想で H^\hat H を作れる。

典型ミス

位相ゲートやアダマールの生成では、全体位相(H^\hat H に ∝I^\propto \hat I を足すこと)を無視して良いかどうかを答案で曖昧にしないこと。設問が「U^(τ)=⋯\hat U(\tau)=\cdots」と等号で指定している場合は、全体位相込みで一致する H^\hat H をひとつ提示するのが安全である。

解答

[1] 時間発展演算子

ハミルトニアンが時間に依存しないとき、シュレーディンガー方程式 iℏddt∣ψ(t)⟩=H^∣ψ(t)⟩ i\hbar \frac{d}{dt}\ket{\psi(t)}=\hat H\ket{\psi(t)} の解は、初期時刻 t0t_0 に対して ∣ψ(t)⟩=U^(t−t0)∣ψ(t0)⟩,U^(t−t0)=exp⁡ ⁣(−iℏH^(t−t0)) \ket{\psi(t)}=\hat U(t-t_0)\ket{\psi(t_0)},\qquad \hat U(t-t_0)=\exp\!\left(-\frac{i}{\hbar}\hat H (t-t_0)\right) で与えられる(U^(0)=I^\hat U(0)=\hat I を満たす)。

[2] U^\hat U のユニタリ性

H^†=H^\hat H^\dagger=\hat H より U^(t−t0)†=exp⁡ ⁣(+iℏH^(t−t0)). \hat U(t-t_0)^\dagger =\exp\!\left(+\frac{i}{\hbar}\hat H (t-t_0)\right). H^\hat H は自分自身と可換なので指数の性質から U^†U^=exp⁡ ⁣(+iℏH^Δt)exp⁡ ⁣(−iℏH^Δt)=I^ \hat U^\dagger \hat U =\exp\!\left(+\frac{i}{\hbar}\hat H\Delta t\right) \exp\!\left(-\frac{i}{\hbar}\hat H\Delta t\right) =\hat I であり、U^\hat U はユニタリーである(同様に U^U^†=I^\hat U\hat U^\dagger=\hat I)。

[3] ノルム保存

[2][2] より ⟨ψ(t)∣ψ(t)⟩=⟨ψ(t0)∣U^†U^∣ψ(t0)⟩=⟨ψ(t0)∣ψ(t0)⟩ \braket{\psi(t)|\psi(t)} =\bra{\psi(t_0)}\hat U^\dagger \hat U\ket{\psi(t_0)} =\braket{\psi(t_0)|\psi(t_0)} で、内積(ノルム)は時間に依存しない。

以降は t0=0, t=τt_0=0,\ t=\tau とする

[4] H^=aσ^z\hat H=a\hat\sigma_z のとき

σ^z\hat\sigma_z の固有値は ±1\pm 1 なので U^(τ)=exp⁡ ⁣(−iℏaτσ^z)=(e−iaτ/ℏ00e+iaτ/ℏ). \hat U(\tau)=\exp\!\left(-\frac{i}{\hbar}a\tau \hat\sigma_z\right) =\begin{pmatrix} e^{-i a\tau/\hbar} & 0\\ 0 & e^{+i a\tau/\hbar} \end{pmatrix}.

[5] 位相シフト diag(1,eiφ)\mathrm{diag}(1,e^{i\varphi}) を与える H^\hat H

対角行列 H^=diag(E0,E1)\hat H=\mathrm{diag}(E_0,E_1) なら U^(τ)=diag(e−iE0τ/ℏ,e−iE1τ/ℏ)\hat U(\tau)=\mathrm{diag}(e^{-iE_0\tau/\hbar},e^{-iE_1\tau/\hbar}) である。 U^(τ)=diag(1,eiφ)\hat U(\tau)=\mathrm{diag}(1,e^{i\varphi}) を満たす一例として H^=(000−ℏφ/τ) \hat H= \begin{pmatrix} 0 & 0\\ 0 & -\hbar\varphi/\tau \end{pmatrix} を取ればよい。

[6] NOT 変換 U^NOT=σ^x\hat U_{\mathrm{NOT}}=\hat\sigma_x を与える H^\hat H

σ^x\hat\sigma_x はユニタリーかつエルミートで σ^x2=I^\hat\sigma_x^2=\hat I を満たす。 したがって U^(τ)=exp⁡ ⁣(−iℏH^τ)=σ^x \hat U(\tau) =\exp\!\left(-\frac{i}{\hbar}\hat H\tau\right) =\hat\sigma_x となる H^\hat H の一例は H^=πℏ2τ (I^−σ^x)=πℏ2τ(1−1−11) \hat H=\frac{\pi\hbar}{2\tau}\,(\hat I-\hat\sigma_x) =\frac{\pi\hbar}{2\tau} \begin{pmatrix} 1 & -1\\ -1 & 1 \end{pmatrix} である(固有値が 0, πℏ/τ0,\ \pi\hbar/\tau となり、指数を取ると +1,−1+1,-1 を与える)。

[7] H^=a(n⋅σ^)+bI^\hat H=a(\boldsymbol n\cdot\hat{\boldsymbol\sigma})+b\hat I のとき

n=(sin⁡θcos⁡ϕ,sin⁡θsin⁡ϕ,cos⁡θ)\boldsymbol n=(\sin\theta\cos\phi,\sin\theta\sin\phi,\cos\theta) とおくと n⋅σ^=nxσ^x+nyσ^y+nzσ^z\boldsymbol n\cdot\hat{\boldsymbol\sigma}=n_x\hat\sigma_x+n_y\hat\sigma_y+n_z\hat\sigma_z は (n⋅σ^)2=I^(\boldsymbol n\cdot\hat{\boldsymbol\sigma})^2=\hat I を満たす(∣n∣=1|\boldsymbol n|=1 を用いる)。 さらに I^\hat I は全てと可換なので U^(τ)=e−ibτ/ℏexp⁡ ⁣(−iℏaτ(n⋅σ^)). \hat U(\tau) =e^{-ib\tau/\hbar}\exp\!\left(-\frac{i}{\hbar}a\tau(\boldsymbol n\cdot\hat{\boldsymbol\sigma})\right). 指数関数をテイラー展開し、偶数次と奇数次に分けると exp⁡ ⁣(−iγ(n⋅σ^))=I^cos⁡γ−i(n⋅σ^)sin⁡γ,γ=aτℏ. \exp\!\left(-i\gamma(\boldsymbol n\cdot\hat{\boldsymbol\sigma})\right) =\hat I\cos\gamma-i(\boldsymbol n\cdot\hat{\boldsymbol\sigma})\sin\gamma, \qquad \gamma=\frac{a\tau}{\hbar}. よって U^(τ)=e−ibτ/ℏ(I^cos⁡γ−i(n⋅σ^)sin⁡γ). \hat U(\tau) =e^{-ib\tau/\hbar}\Bigl(\hat I\cos\gamma -i(\boldsymbol n\cdot\hat{\boldsymbol\sigma})\sin\gamma\Bigr).

[8] アダマール変換 U^H\hat U_H を与える H^\hat H

U^H\hat U_H はエルミートかつユニタリーで U^H2=I^\hat U_H^2=\hat I を満たす。 よって [6][6] と同様に H^=πℏ2τ (I^−U^H)=πℏ2τ((1001)−12(111−1)) \hat H=\frac{\pi\hbar}{2\tau}\,(\hat I-\hat U_H) =\frac{\pi\hbar}{2\tau} \left( \begin{pmatrix}1&0\\0&1\end{pmatrix} -\frac{1}{\sqrt2}\begin{pmatrix}1&1\\1&-1\end{pmatrix} \right) とすれば exp⁡ ⁣(−iℏH^τ)=U^H\exp\!\left(-\frac{i}{\hbar}\hat H\tau\right)=\hat U_H が成り立つ。

[9] 制御位相シフト U^C−φ\hat U_{C-\varphi} を与える H^\hat H

U^C−φ=diag(1,1,1,eiφ)\hat U_{C-\varphi}=\mathrm{diag}(1,1,1,e^{i\varphi}) を与える一例は H^=diag ⁣(0,0,0,−ℏφ/τ) \hat H=\mathrm{diag}\!\left(0,0,0,-\hbar\varphi/\tau\right) である(U^(τ)=exp⁡(−iH^τ/ℏ)=U^C−φ\hat U(\tau)=\exp(-i\hat H\tau/\hbar)=\hat U_{C-\varphi})。

[10] 制御 NOT U^CNOT\hat U_{\mathrm{CNOT}} を与える H^\hat H

U^CNOT\hat U_{\mathrm{CNOT}} は U^CNOT2=I^4\hat U_{\mathrm{CNOT}}^2=\hat I_4 を満たす(エルミートかつユニタリー)。 したがって [6][6] と同様に H^=πℏ2τ (I^4−U^CNOT)=πℏ2τ(00000000001−100−11) \hat H=\frac{\pi\hbar}{2\tau}\,(\hat I_4-\hat U_{\mathrm{CNOT}}) =\frac{\pi\hbar}{2\tau} \begin{pmatrix} 0&0&0&0\\ 0&0&0&0\\ 0&0&1&-1\\ 0&0&-1&1 \end{pmatrix} とすれば exp⁡ ⁣(−iℏH^τ)=U^CNOT\exp\!\left(-\frac{i}{\hbar}\hat H\tau\right)=\hat U_{\mathrm{CNOT}} を与える。

最終答

時間に依存しない H^\hat H では U^(Δt)=exp⁡ ⁣(−iℏH^Δt)\hat U(\Delta t)=\exp\!\left(-\frac{i}{\hbar}\hat H\Delta t\right)。H^†=H^⇒U^†U^=I^\hat H^\dagger=\hat H\Rightarrow \hat U^\dagger\hat U=\hat I で、⟨ψ(t)∣ψ(t)⟩\braket{\psi(t)|\psi(t)} は保存される。H^=aσ^z⇒U^(τ)=diag(e−iaτ/ℏ,eiaτ/ℏ)\hat H=a\hat\sigma_z\Rightarrow \hat U(\tau)=\mathrm{diag}(e^{-ia\tau/\hbar},e^{ia\tau/\hbar})。diag(1,eiφ)\mathrm{diag}(1,e^{i\varphi}) は H^=diag(0,−ℏφ/τ)\hat H=\mathrm{diag}(0,-\hbar\varphi/\tau) で生成できる。σ^x,U^H,U^CNOT\hat\sigma_x,\hat U_H,\hat U_{\mathrm{CNOT}} のように U2=I^U^2=\hat I を満たすゲートは H^=πℏ2τ(I^−U)\hat H=\frac{\pi\hbar}{2\tau}(\hat I-U) の一例で生成できる。

第2問 等方弾性体の応力と弾性波反射

押さえるべき対応

p(k)\boldsymbol p^{(k)} は法線が +xk+x_k の面に作用する牽引(トラクション)で、成分は σjk=pj(k)\sigma_{jk}=p_j^{(k)} と同一視できる。従って p(2)=0\boldsymbol p^{(2)}=\boldsymbol 0 は σ12=σ22=0\sigma_{12}=\sigma_{22}=0(この設定では σ32=0\sigma_{32}=0 も自動的)という自由表面条件である。

位相整合の考え方

境界条件は「境界上の任意の x1x_1 と任意の時刻 tt」で成立する必要がある。したがって、境界上の指数因子 ei(k1x1−ωt)e^{i(k_1x_1-\omega t)} の形が一致しない項どうしは打ち消し合えない。この観察だけで ω\omega の保存、k1k_1 の保存(スネル則の本質)が一気に出る。

なぜ αt=0,π/4\alpha_t=0,\pi/4 が特別か

横波は ∇⋅u=0\nabla\cdot\boldsymbol u=0 なので、自由表面条件は実質的に σ12=σ22=0\sigma_{12}=\sigma_{22}=0 の2本を同時に満たす必要がある。入射・反射の2つの振幅だけで2条件を満たせるのは、係数行列が特異になる(sin⁡2αt\sin 2\alpha_t か cos⁡2αt\cos 2\alpha_t がゼロになる)角度に限られ、原題の 0≤αt<π/20\le\alpha_t<\pi/2 では αt=0,π/4\alpha_t=0,\pi/4 がそれに当たる。垂直入射の端点 αt=0\alpha_t=0 も除外しない。

第3問 理想ボース気体とボース凝縮

計算の核

凝縮の有無は「μ=0\mu=0 としても励起状態に入れられる粒子数 Nex(T)N_{\mathrm{ex}}(T) が有限かどうか」で決まる。低エネルギーでの状態密度 D(ε)D(\varepsilon) のべきが臨界挙動を支配し、3次元自由気体では D(ε)∝ε1/2D(\varepsilon)\propto \varepsilon^{1/2} のため Nex(T)∝T3/2N_{\mathrm{ex}}(T)\propto T^{3/2} が有限になる一方、2次元自由気体では D(ε)D(\varepsilon) が定数となり赤外発散が出る。

調和トラップの直感

箱では「体積 L3L^3」が状態数のスケールを与えるが、調和トラップではエネルギー殻の縮退度が増えるため D(ε)∝εD(\varepsilon)\propto \varepsilon となり、2次元でも Nex(T)∝T2N_{\mathrm{ex}}(T)\propto T^2 が有限になって凝縮が可能になる。

第4問 誘起双極子による電磁波散乱

どこが放射(1/r1/r)か

双極子の時間変化が作る波(放射場)は、遠方で E,B∝1/rE,B\propto 1/r となる成分である。スカラーポテンシャルの 1/r21/r^2 成分は近接場で、遠方では寄与が落ちる。

横波になる理由

−∇ϕ-\nabla\phi と −∂tA-\partial_t\boldsymbol A の er\boldsymbol e_r 成分が遠方で相殺され、E\boldsymbol E が eθ\boldsymbol e_\theta のみになる。結果として B\boldsymbol B は eφ\boldsymbol e_\varphi 方向で、S\boldsymbol S は er\boldsymbol e_r 方向(放射は外向き)になる。

周波数混合

α(t)\alpha(t) が振動すると p(t)p(t) は入射周波数 ω0\omega_0 の振幅変調になり、和周波・差周波 ω0±ωv\omega_0\pm\omega_v が現れる。散乱電場は p¨\ddot p に比例するため、強度は概ね各成分の ω4\omega^4 スケールで重み付けされる。

物理学の解答・解説PDF

2012〜2026年度(15年度) · ¥3,200

解答PDFの商品ページへ

年度別の解答例

  • 2026年度

    剛体回転と磁気トルク / 調和振動子と突然近似 / ゴム弾性の統計力学

  • 2025年度

    減衰振動と揺動応答 / 一次元障壁散乱 / 平均場イジング模型

  • 2024年度

    剛体振り子と磁気モーメント測定 / 1次元光共振器 / 固体の比熱

  • 2023年度

    量子井戸と対称性 / 支点加振振り子 / 古典気体の回転平衡

  • 2022年度

    回転磁場によるスピン遷移と Ramsey 干渉 / 一次元格子振動 / 二元合金の相分離

  • 2021年度(このページ)

    量子ビットの時間発展とゲート生成 / 等方弾性体の応力と弾性波反射 / 理想ボース気体とボース凝縮

  • 2020年度

    固定端ばね鎖の固有振動 / 一軸性結晶中の電磁波 / 理想気体と吸着等温式

  • 2019年度

    午前・軌道力学と重力勾配 / 午前・誘導加熱と円環電流 / 午後・デルタ障壁と二粒子干渉

  • 2018年度

    午前・半円筒上の滑りと転がり / 午前・球形コンデンサと漏れ電流 / 午後・デルタ井戸と二粒子交換

  • 2017年度

    午前・摩擦つき連結物体と安定性 / 午前・電気双極子と鏡像法 / 午後・井戸中の二フェルミ粒子

  • 2016年度最終答まで公開

    午前・振り子と二重振り子 / 午前・導体リングと電磁誘導 / 午後・調和振動子のコヒーレント状態

  • 2015年度最終答まで公開

    午前・跳ね返る剛体球 / 午前・線電荷と円筒導体の鏡像法 / 午後・二重井戸と二粒子トンネル

  • 2014年度最終答まで公開

    午前・中心力とラプラス--ルンゲ--レンツベクトル / 午前・電磁場とLC振動 / 午後・スピン軌道相互作用

  • 2013年度最終答まで公開

    粗い円柱のつり合いと滑りから転がりへの移行 / 平行平板間の空間電荷制限電流 / 調和振動子の遷移と非調和摂動

  • 2012年度最終答まで公開

    円板に巻き付く糸につながれた質点 / 円柱導体のジュール熱とポインティングベクトル / スピン1の磁気共鳴と四極子摂動

2021年度の出題内容

大問分野主題解説の小見出し最終答
第1問量子力学量子ビットの時間発展とゲート生成見通し / 指数の扱いあり
第2問—等方弾性体の応力と弾性波反射押さえるべき対応 / 位相整合の考え方あり
第3問—理想ボース気体とボース凝縮計算の核 / 調和トラップの直感あり
第4問—誘起双極子による電磁波散乱どこが放射(1/r)か / 横波になる理由あり

同じテーマを出した他大学・他年度

前年度(2020年度)との違い

大問数
2020年度 4問 → 2021年度 4問
2021年度で新しく出たテーマ
テイラー展開

2020年度の解答例

東京大学の他の解答PDF

院試対策ガイド

東京大学の公式過去問

公開年度範囲
2022〜2026
公式の解答例
解答例の掲載なし

東京大学の公式過去問PDF