院試hub

東京科学大(旧東工大) 理学院 地球惑星科学系 専門科目 2023年度 院試 過去問 解答例・解説

大学
東京科学大学
研究科・専攻
理学院 地球惑星科学系
科目
地球惑星科学系 専門科目
年度
2023年度

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

東京科学大学 理学院 地球惑星科学系 専門科目の解答PDF2021〜2025年度(5年度)

¥2,880

商品内容・見本・購入へ

第1問 地球内部・海洋・大気

方針

本問は暗記項目の羅列ではなく、観測・分配・相転移・物質循環を「なぜそう推定できるか」まで結びつけて書くと安定する。特に 1-2 は Hf と W の親和性の違い、1-3 はクラペイロン勾配の符号、1-5 は元素ごとに供給源と沈殿・変質反応を分けることが採点上の要点である。

検算と注意点

スラブ滞留では 410 km 相転移の正のクラペイロン勾配だけを述べても不十分で、滞留の障壁として重要な 660 km 付近の負のクラペイロン勾配に触れる必要がある。中央海嶺では「高温だから融ける」ではなく、上昇に伴う減圧でソリダスを横切ると書く。海水中の S は硫化物ではなく、酸化的な現在の海洋では硫酸イオンが主成分である点を間違えやすい。

試験で書くべきポイント

短答では、各設問に一つずつ因果の鎖を入れる。例えば「S 波が通らないので外核は液体」「Hf はマントル、W は核へ分配」「水蒸気凝結の潜熱で湿潤断熱減率になる」のように、観測または過程と結論を一文内で結ぶと、行数が限られていても失点しにくい。

第2問 化学分析法

方針

前半は電磁波と顕微鏡の基本式、後半は鉱物分析の原理を問う構成である。数値問題では E=hν=hc/λE=h\nu=hc/\lambda を最初に書く。鉱物式は酸化物重量を mol に直し、酸素数または陽イオン比で規格化する。

検算

波長 10−5 m10^{-5}\,\mathrm{m} は赤外線なので、エネルギーが可視光の数 eV より小さい 0.1 eV0.1\,\mathrm{eV} 程度になるのは妥当である。試料 1 は (Mg+Fe)/Si≃1(\mathrm{Mg+Fe})/\mathrm{Si}\simeq1 で輝石型、試料 2 は Mg/Si≃2\mathrm{Mg}/\mathrm{Si}\simeq2 でかんらん石型となり、酸素数もそれぞれ O3\mathrm{O_3}, O4\mathrm{O_4} に対応する。

典型ミス

電子の波長を非相対論式だけで答えると、加速電子の運動量増加を落とす。X 線分析では連続 X 線は背景であり、元素同定に直接使うのは特性 X 線である。黒体放射では、温度上昇でピークが長波長側へ行くと逆に書かないこと。

試験で書くべきポイント

Bragg 条件の導出では、角度 2θ2\theta は入射線と回折線の開きであり、式に入る θ\theta は結晶面との角であることを明記するとよい。分析法の説明は「何を測るか」と「地球惑星科学で何に使うか」を一組で書く。

第3問 線形方程式と熱方程式

方針

前半は固有ベクトル展開で各成分を 1 次元の指数減衰に分ける問題である。差分化後は Euler 法の増幅因子 σi=1−hλi\sigma_i=1-h\lambda_i を見れば安定性が判定できる。後半は偶関数の Fourier 余弦変換を使って熱方程式の基本解を導く問題である。

主要式

熱方程式は ∂tu=a∂x2u\partial_t u=a\partial_x^2u なので、波数 kk のモードは exp⁡(−ak2t)\exp(-ak^2t) で減衰する。初期分布を余弦変換し、逆変換で戻すと、ガウス核 G(x−ξ,t)=12πatexp⁡{−(x−ξ)24at} G(x-\xi,t)=\frac{1}{2\sqrt{\pi at}} \exp\left\{-\frac{(x-\xi)^2}{4at}\right\} が現れる。係数の 2πat2\sqrt{\pi at} は、全実軸で積分したとき 11 になる正規化定数である。

検算

安定条件で h→0h\to0 とすれば σi→1\sigma_i\to1 だが ∣σi∣<1|\sigma_i|<1 を満たし、時間発展はゆっくり減衰する。逆に hλi>2h\lambda_i>2 では σi<−1\sigma_i<-1 となり、符号を反転しながら発散する。熱核は t↓0t\downarrow0 でデルタ関数に近づき、u(x,0)=φ(x)u(x,0)=\varphi(x) を回復する。

典型ミス

差分法の条件を 1−hλi<11-h\lambda_i<1 だけで終えると、負側の発散 σi<−1\sigma_i<-1 を見落とす。熱方程式では kk 積分を 00 からではなく全実軸で扱っているため、係数 1/(2π)1/(2\pi) と積の公式の 2 倍を丁寧に追う必要がある。

解答

3-1-1. v(t)=∑i=1Nfi(t)wi v(t)=\sum_{i=1}^N f_i(t)w_i を代入すると、 ∑i=1Nfi′(t)wi=−∑i=1Nfi(t)Awi=−∑i=1Nλifi(t)wi. \sum_{i=1}^N f_i'(t)w_i =-\sum_{i=1}^N f_i(t)Aw_i =-\sum_{i=1}^N \lambda_i f_i(t)w_i . 固有値が異なり固有ベクトルが一次独立なので、 fi′(t)=−λifi(t) f_i'(t)=-\lambda_i f_i(t) であり、 fi(t)=Ciexp⁡(−λit). f_i(t)=C_i\exp(-\lambda_i t). したがって一般解は v(t)=∑i=1NCie−λitwi. v(t)=\sum_{i=1}^N C_i e^{-\lambda_i t}w_i .

3-1-2. 差分式は vn+1=vn−hAvn=(I−hA)vn v_{n+1}=v_n-hAv_n=(I-hA)v_n である。よって B=I−hA,vn=Bnv0. B=I-hA,\qquad v_n=B^n v_0 .

3-1-3. Bwi=(I−hA)wi=(1−hλi)wi. Bw_i=(I-hA)w_i=(1-h\lambda_i)w_i . したがって BB の規格化固有ベクトルは zi=wiz_i=w_i と選べ、対応する固有値は σi=1−hλi \sigma_i=1-h\lambda_i である。

3-1-4. P=(z1 z2 ⋯ zN) P=(z_1\ z_2\ \cdots\ z_N) とおくと B=P(σ10⋱0σN)P−1. B=P \begin{pmatrix} \sigma_1&&0\\ &\ddots&\\ 0&&\sigma_N \end{pmatrix} P^{-1}. したがって Bn=P(σ1n0⋱0σNn)P−1. B^n=P \begin{pmatrix} \sigma_1^n&&0\\ &\ddots&\\ 0&&\sigma_N^n \end{pmatrix} P^{-1}.

3-1-5. 任意の初期ベクトルに対して vn→0v_n\to0 となるには ∣σi∣<1(i=1,…,N) |\sigma_i|<1\qquad (i=1,\ldots,N) が必要である。ここで σi=1−hλi\sigma_i=1-h\lambda_i だから、λi>0\lambda_i>0 のもとでは −1<1−hλi<1⟺0<hλi<2. -1<1-h\lambda_i<1 \quad\Longleftrightarrow\quad 0<h\lambda_i<2 . 全ての ii について 0<h<2max⁡iλi 0<h<\frac{2}{\max_i\lambda_i} を満たす必要がある。

3-2-1. I(ξ)=∫−∞∞e−b2cos⁡(bξ) db I(\xi)=\int_{-\infty}^{\infty}e^{-b^2}\cos(b\xi)\,db を微分すると I′(ξ)=−∫−∞∞be−b2sin⁡(bξ) db. I'(\xi)=-\int_{-\infty}^{\infty}b e^{-b^2}\sin(b\xi)\,db . be−b2=−(1/2)(e−b2)′b e^{-b^2}=-(1/2)(e^{-b^2})' を用いて部分積分すれば、端点項は消えて I′(ξ)=12∫−∞∞(e−b2)′sin⁡(bξ) db=−ξ2I(ξ). I'(\xi) =\frac12\int_{-\infty}^{\infty}(e^{-b^2})'\sin(b\xi)\,db =-\frac{\xi}{2}I(\xi). よって I′I=−ξ2⟹I(ξ)=Cexp⁡(−ξ24). \frac{I'}{I}=-\frac{\xi}{2} \quad\Longrightarrow\quad I(\xi)=C\exp\left(-\frac{\xi^2}{4}\right). I(0)=πI(0)=\sqrt{\pi} より I(ξ)=πexp⁡(−ξ24). I(\xi)=\sqrt{\pi}\exp\left(-\frac{\xi^2}{4}\right).

3-2-2-1. u(x,t)=W(x)T(t)u(x,t)=W(x)T(t) とおくと 1aT′T=W′′W. \frac{1}{a}\frac{T'}{T}=\frac{W''}{W}. 時間と空間に分かれた量が常に等しいので、これを −k2-k^2 とおく。時間方向は T(t)=exp⁡(−ak2t) T(t)=\exp(-ak^2t) であり、時間とともに減衰する条件を満たす。空間方向は W′′+k2W=0 W''+k^2W=0 で、偶関数条件から sin⁡(kx)\sin(kx) 成分は消える。したがって U(x,t)=A(k)cos⁡(kx)exp⁡(−ak2t). U(x,t)=A(k)\cos(kx)\exp(-ak^2t).

3-2-2-2. φ(ξ)\varphi(\xi) は偶関数であるから、 cos⁡{k(ξ±x)}=cos⁡(kξ)cos⁡(kx)∓sin⁡(kξ)sin⁡(kx) \cos\{k(\xi\pm x)\} =\cos(k\xi)\cos(kx)\mp \sin(k\xi)\sin(kx) を用いると、φ(ξ)sin⁡(kξ)\varphi(\xi)\sin(k\xi) は奇関数になる。全実軸積分では奇関数部分が消えるので ∫−∞∞φ(ξ)cos⁡{k(ξ+x)} dξ=∫−∞∞φ(ξ)cos⁡{k(ξ−x)} dξ. \int_{-\infty}^{\infty}\varphi(\xi)\cos\{k(\xi+x)\}\,d\xi = \int_{-\infty}^{\infty}\varphi(\xi)\cos\{k(\xi-x)\}\,d\xi .

3-2-2-3. 与えられた表示に A(k)A(k) を代入し、積分順序を入れ替えると u(x,t)=12π∫−∞∞φ(ξ)[∫−∞∞cos⁡(kξ)cos⁡(kx)e−ak2t dk]dξ. u(x,t)=\frac{1}{2\pi}\int_{-\infty}^{\infty}\varphi(\xi) \left[\int_{-\infty}^{\infty} \cos(k\xi)\cos(kx)e^{-ak^2t}\,dk \right]d\xi . 積の公式と前問の等式により、内側の積分は ∫−∞∞e−ak2tcos⁡{k(x−ξ)} dk \int_{-\infty}^{\infty}e^{-ak^2t}\cos\{k(x-\xi)\}\,dk だけを考えればよい。3-2-1 の結果を b=at k,η=x−ξat b=\sqrt{at}\,k,\qquad \eta=\frac{x-\xi}{\sqrt{at}} として用いると ∫−∞∞e−ak2tcos⁡{k(x−ξ)} dk=πatexp⁡{−(x−ξ)24at}. \int_{-\infty}^{\infty}e^{-ak^2t}\cos\{k(x-\xi)\}\,dk = \sqrt{\frac{\pi}{at}} \exp\left\{-\frac{(x-\xi)^2}{4at}\right\}. したがって u(x,t)=∫−∞∞φ(ξ) 12πatexp⁡{−(x−ξ)24at} dξ. u(x,t)= \int_{-\infty}^{\infty} \varphi(\xi)\, \frac{1}{2\sqrt{\pi at}} \exp\left\{-\frac{(x-\xi)^2}{4at}\right\}\,d\xi .

最終答

v(t)=∑iCie−λitwiv(t)=\sum_i C_i e^{-\lambda_i t}w_i、B=I−hAB=I-hA、σi=1−hλi\sigma_i=1-h\lambda_i、安定条件は ∣σi∣<1|\sigma_i|<1、すなわち 0<h<2/max⁡iλi0<h<2/\max_i\lambda_i。熱方程式の解はガウス核との畳み込みで u(x,t)=∫φ(ξ){2πat}−1exp⁡[−(x−ξ)2/(4at)] dξu(x,t)=\int \varphi(\xi)\{2\sqrt{\pi at}\}^{-1}\exp[-(x-\xi)^2/(4at)]\,d\xi。

第4問 剛体回転と衝撃波

方針

剛体部分は、拘束条件 x=hθx=h\theta を先に立てれば 1 自由度のラグランジアンになる。衝撃波部分は、保存量を「質量流束、運動量流束、エネルギー流束」の順に積分形で扱うと式変形が見通しやすい。

主要式

円板の降下加速度 x¨=g1+I/(Mh2) \ddot{x}=\frac{g}{1+I/(Mh^2)} は、並進運動に加えて回転運動へも重力エネルギーが分配されることを表す。慣性モーメントが大きいほど同じ重力で角加速度が小さくなり、降下加速度も小さくなる。

衝撃波では s2−s1=cvln⁡(βαγ) s_2-s_1=c_v\ln\left(\frac{\beta}{\alpha^\gamma}\right) を最後に評価する。圧力増加だけでなく密度増加も同時に起こるため、β>1\beta>1 だけではエントロピー増加の証明にならない。

検算と注意点

I=0I=0 とすれば円板は回転にエネルギーを使わず、x¨=g\ddot{x}=g となる。薄い輪の極限 I=Mh2I=Mh^2 では x¨=g/2\ddot{x}=g/2 で、よく知られたヨーヨー型の結果と一致する。衝撃波では M1=1M_1=1 で α=β=1\alpha=\beta=1 となり、衝撃波が消える極限と合う。

典型ミス

ポテンシャルエネルギーは下向き座標を正に取っているため U=−MgxU=-Mgx である。符号を逆にすると運動方程式の加速度が上向きになる。衝撃波では u2=u1/αu_2=u_1/\alpha を使わずに式を進めると、密度比と速度比を取り違えやすい。

第5問 化学結合

方針

混成軌道の問題では、電子数だけでなく「どの軌道がどの方向を向くか」を述べる。酸塩基計算では、5-9 は弱酸単独、5-10 は共役酸塩基対の緩衝液として分ける。

主要式

5-9 の計算は Ka=x2C−x K_a=\frac{x^2}{C-x} が出発点である。5-10 は Henderson--Hasselbalch 型の pH=pKa+log⁡[F−][HF] \mathrm{pH}=\mathrm{p}K_a+\log\frac{[\mathrm{F^-}]}{[\mathrm{HF}]} を使う。pH 3.5 と pKa 3.20 の差が 0.30 なので、比はほぼ 2 である。

検算と注意点

5-9 で [H+]≃0.014 M[\mathrm{H^+}]\simeq0.014\,\mathrm{M} に対し、求めた全濃度は 0.330 M0.330\,\mathrm{M} であり、弱酸として解離度は数 % 程度で自然である。5-10 の CsF は 0.66 mol 程度なので、分子量 152 を掛けて約 100 g となる。桁として 10 g ではなく 100 g である点に注意する。

典型ミス

NH4+\mathrm{NH_4^+} の結合を「1 本だけ配位結合なので特別」としたままにしない。生成後は 4 本の N--H 結合が等価に扱える。HF が弱酸である理由を「F が電気陰性だから H が取れやすい」とだけ書くと逆方向で、強い H--F 結合を明記する必要がある。

第6問 火山岩と相平衡

方針

火山岩分類は観察記載、全岩化学組成、相図をつなげて答える。薄片記載では「干渉色・消光・双晶・不透明性」を鉱物名に結びつけ、相図では最初の液、共存固相、液組成の移動方向を分けて書く。

主要式

液相濃集元素が固相に入らない場合、源岩濃度 C0C_0、液濃度 ClC_l、融解度 FF は Cl=C0F C_l=\frac{C_0}{F} で結ばれる。結晶分化後の残液濃度は、成分質量を残液質量で割ればよい。かんらん石を 40 wt% 除く計算では、残液質量が 60 g になる点が最も重要である。

検算と注意点

1250 ∘C^{\circ}\mathrm{C} の液は Na2O+K2O\mathrm{Na_2O+K_2O} が高いので、融解度は小さく見積もられる。1400 ∘C^{\circ}\mathrm{C} の液はより多量の融解で希釈されるため、アルカリ濃度が低い。これは計算結果 6.6%6.6\% と 19%19\% の大小関係と整合する。

典型ミス

粗面岩質という語だけからホットスポット的なアルカリ岩と決めつけない。ここでは島弧という造構場を保ったまま、岩石系列と結晶分化で説明する。相図問題では、無水条件と含水条件を混同しない。水は単に温度を少し下げるだけでなく、低温で SiO2\mathrm{SiO_2} に富む液を生じやすくする。

第7問 電磁波吸収と量子統計

方針

電磁波部分は、入射エネルギーを Poynting ベクトルで、吸収を減衰振動子への仕事率で評価する。統計力学部分は、エネルギー固有値の体積依存性 ϵi∝V−2/3\epsilon_i\propto V^{-2/3} だけを使えば、Bose 統計の詳細に立ち入らずに状態方程式が出る。

主要式

仕事率の時間平均で吸収が残るのは、変位の電場に対する位相遅れを表す CC の項である。抵抗がなければ C=0C=0 となり、1 周期平均の吸収は消える。

検算

ω≪Ω\omega\ll\Omega では振動子は外力にほぼ追随し、位相遅れは小さい。そのため吸収が低周波で ω2\omega^2 に落ちるのは自然である。自由粒子の関係 pV=2E/3pV=2E/3 は古典理想気体と同じ形で、非相互作用粒子の運動エネルギーが二次分散を持つことの帰結である。

典型ミス

電磁波のエネルギー密度を電場分だけ ϵ0E2/2\epsilon_0E^2/2 としてしまうと、磁場分を落として係数が半分になる。平面波では電場と磁場のエネルギーが等しく、全エネルギー密度は ϵ0E2\epsilon_0E^2 である。

第8問 ラプラス変換と Chebyshev 多項式

方針

8-1 は、減衰を持つ 2 階方程式をラプラス変換で畳み込み応答に直す問題である。8-2 は、与えられた母関数で定義される多項式が Chebyshev 多項式であることを示し、直交性を使って arccos⁡x\arccos x を展開する。

主要式

確率過程の分散では、白色雑音の相関 E[f(r)f(r′)]=σ2δ(r−r′) E[f(r)f(r')]=\sigma^2\delta(r-r') により二重積分が一重積分へ潰れる。ここを書けば、最後の式は通常の指数関数の積分だけになる。

Chebyshev 展開では、重み (1−x2)−1/2(1-x^2)^{-1/2} は x=cos⁡θx=\cos\theta によって単なる dθd\theta に変わる。この変換が直交性の本体である。

検算

t↓0t\downarrow0 では分散式の括弧内は 0 から始まるため、初期条件 x(0)=0x(0)=0 と合う。長時間では指数項が消え、 V[x(t)]∼σ2k2t V[x(t)]\sim \frac{\sigma^2}{k^2}t となり、ランダム加速度を積分した粒子位置が拡散的に広がることを表す。Chebyshev 係数では arccos⁡x\arccos x が θ=π/2\theta=\pi/2 を中心に奇対称な余弦級数になるため、偶数次係数が 0 になる。

典型ミス

ラプラス変換で初期条件項を落とすのは、初期条件が 0 だからであり、一般には sx(0)sx(0) と x′(0)x'(0) が残る。8-2-5 の a0a_0 は通常の Fourier 余弦級数の定義に従い、他の係数と同じ 2/π2/\pi を掛けない点に注意する。

地球惑星科学系 専門科目の解答・解説PDF

2021〜2025年度(5年度) · ¥2,880

解答PDFの商品ページへ

年度別の解答例

  • 2025年度

    ルジャンドル多項式と多次元ガウス積分 / 中心力による捕獲断面積と抵抗を受ける斜方投射 / 自己重力ガス球とカノニカル分布

  • 2024年度

    確率・ルジャンドル多項式・中心力 / 解析力学と電磁場 / 地球内部熱史と放射平衡

  • 2023年度(このページ)

    地球内部・海洋・大気 / 化学分析法 / 線形方程式と熱方程式

  • 2022年度

    化学熱力学と同位体分別 / 地球表層史と付加体 / 微分方程式と特殊関数

  • 2021年度

    隕石と母天体の熱史 / プレートテクトニクスと火成岩 / 線形代数・複素積分・確率

2023年度の出題内容

大問分野主題解説の小見出し最終答
第1問固体地球科学地球内部・海洋・大気検算と注意点あり
第2問—化学分析法—あり
第3問—線形方程式と熱方程式主要式あり
第4問力学剛体回転と衝撃波主要式 / 検算と注意点あり
第5問—化学結合主要式 / 検算と注意点あり
第6問地質火山岩と相平衡主要式 / 検算と注意点あり
第7問量子力学電磁波吸収と量子統計主要式あり
第8問—ラプラス変換と Chebyshev 多項式主要式あり

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

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

東京科学大学の他の解答PDF

院試対策ガイド

東京科学大学の公式過去問

公開年度範囲
主に2021〜2025年実施(数学系は2020〜2024年実施)
公式の解答例
解答例の掲載なし
英語の扱い
2027年4月・2026年9月入学の理工学系修士は、数学系を除きTOEIC L&R・TOEFL iBT・Home Editionの外部スコアを使用。数学系は独自英語筆答が必須で、外部スコアによる免除はありません。提出期限は系別です。

東京科学大学の公式過去問PDF