量子力学的振動子の分配関数と自由エネルギーについての計算

このQ&Aのポイント
  • N個の量子力学的振動子からなる系の分配関数Z(T,V,N)は、Z(T,V,N)=[2sinh(hω/2kT)]^(-N)で表される。
  • 問題1では、ヘルムホルツの自由エネルギーF(T,V,N)を求める。問題2では、(a)エントロピーS=-dF/dTと(b)内部エネルギーU=F+TSを求める。問題3では、問題2で求めた内部エネルギーUの温度Tに対するグラフの概形を書く。問題4では、量子系ではT→0の極限でUが0にならない理由を考える。
  • 問題1のヘルムホルツの自由エネルギーF(T,V,N)はF=NkTln2sinh(hω/2kT)となる。問題2(a)のエントロピーの計算結果は-NKln2sinh(hω/2kT)-(Nhω/2T){cosh(hω/2kT)/sinh(hω/2kT)}であり、この結果の正しさについて不安がある。問題3以降の解法の手順も教えてほしい。
回答を見る
  • ベストアンサー

分配関数の計算に関して

初歩的な内容かもしれませんが、計算が合っているかどうか不安でしたのでわかる方よろしくお願いします。 ・N個の量子力学的振動子からなる系をカノニカルアンサンブルのもとで考えると、この系の分配関数Z(T,V,N)は    Z(T,V,N)=[2sinh(hω/2kT)]^(-N) として表されている。 ここでの問いとして 1:ヘルムホルツの自由エネルギーF(T,V,N)を求めよ。 2:(a)エントロピーS=-dF/dTおよび(b)内部エネルギーU=F+TSを求めよ。 3:2で求めた(b)で求めたU(内部エネルギー)のT(温度)に対するグラフの概形を書け。 4:古典力学系ではT→0の極限で、Uは0になるが、3で得られた量子系ではUは0にならない。それはなぜか。 1はF=NkTln2sinh(hω/2kT)とすぐわかったのですが、 2(a)のエントロピーを求める計算結果が -NKln2sinh(hω/2kT)-(Nhω/2T){cosh(hω/2kT)/sinh(hω/2kT)} となり、この結果で合っているかどうかが不安です。 よろしければこの結果であっているかどうか、また違っていた場合のその箇所考え方等をご教授願います。 もしよろしければ3以降の問題を解法の手順等もお願いします。 ※hはエイチバーとして判断して下さい。

質問者が選んだベストアンサー

  • ベストアンサー
  • eatern27
  • ベストアンサー率55% (635/1135)
回答No.1

>-NKln2sinh(hω/2kT)-(Nhω/2T){cosh(hω/2kT)/sinh(hω/2kT)} ん~、第二項が微妙に違いますね。 1/TをTで微分した時に負符号が出てくる事を見落としたとか?計算し直してみてください。 >もしよろしければ3以降の問題を解法の手順等もお願いします。 3は、T→0での振る舞いと,T→∞での振る舞いが分かるように図示すればいいでしょう。 4は、零点振動がどうのこうのって話ですね。

minmin529
質問者

お礼

うやむやがなくなりました、ありがとうございます。

関連するQ&A

  • 分配関数(状態和)がわかりません。

    統計力学とかで出てくる分配関数(状態和)がありますが、物理的な意味がよくわかってません。 Σexp(-β・ei)とありますがどういう意味なんでしょうか? またある問題でエネルギー準位ε=(n+1/2)hνのN個の独立な調和振動系子の系があり この調和振動子一個に対する状態和が Z=1/{2sinh(hν/2kB・T)} となることを示せという問題があるんですが問題の意味すらよくわかりません。 一個に対する状態和?という感じです。 どうかお願いします。

  • ヘルムホルツの自由エネルギー

    温度T一定、体積V一定の系は、ヘルムホルツの自由エネルギー   F = E-TS を最小化しようとする、 つまり、力学的エネルギーEを小さくすることと エントロピーSを増大させようとする効果がせめぎあって 最小のFになろうとすると理解しています。 しかし、T一定のときは、Eの値は変わりようがないのでは?といういの疑問で行き詰っています。 根拠は以下の通りです。 系がギブス分布をとる場合は、Eは熱平衡での平均エネルギーなので、   E =(1/z) Σ Hi exp(-Hi/kT) です。 Hiは系が状態毎に持つエネルギーで、これを全状態分数え上げていますし、 その確率も、(1/z) exp(-Hi/kT)で決まっています。 よって、温度Tが決まれば、Eの値は変わらないと思うのです。 また、状態の出現確率が、(1/z) exp(-Hi/kT)で決まってしまうので、 無秩序さであるエントロピーSも決まってしまうようにも思います。 状態の出現確率が、(1/z) exp(-Hi/kT)から変化するのか?などと いろいろな状況を考えてみては混乱しています。 どうか、Fを極小化するために、EやSがどうなっているのか教えてください。 よろしくお願いします。

  • 熱力学関数です。

    熱力学関数の問題です。 熱力学関数の証明問題です。 ・ヘルムホルツの自由エネルギーから内部エネルギーを求める式:U=-T^2[∂(F/T)/∂T]vを導け。 ・ギブスの自由エネルギーからエンタルピーを求める式:H=T^2[∂(G/T)/∂T]pを導け。 どのようなプロセスを踏んで導けがよいのかわかりません。 わかる方は回答よろしくお願いします。

  • 熱統計力学 分配関数について

     E1,E2,E3という3つのエネルギー状態を持てる粒子数N個の系を考える。系の温度はτとする。Nが十分に大きいとき、系の全体のエネルギーがU0となる場合の粒子数Nを求めよ。 上記のような問題が出されました。私の解答としては、 分配関数 Z(τ)=exp(-E1/τ)+exp(-E2/τ)+exp(-E3/τ) よって平均エネルギーは U={E1exp(-E1/τ)+E2exp(-E2/τ)+E3exp(-E3/τ)}/Z(τ) したがって NU=U0;よりN=U0/U  という解答を出しましたが、平均エネルギーがUを単純にN倍したのが系全体のエネルギーU0になるとして良いのか不安です。また、Nを十分大きくとる理由がわかりません。熱力学はまだ習ったばかりで正直この解法も正しいのか疑問です・・・・。  どんなわずかなことでも良いので回答を下さい。

  • エンタルピーの計算

    エンタルピーの計算問題があるのですが、いまいちよく分かりません。 問題 体積一定のボンベ熱量計を用いて炭化タングステンWCを燃焼させたら1 molのWC燃焼に伴う内部エネルギー変化ΔU(300 K)は-1192 kJであった。300 KにおけるΔHはいくらか? このような問題なのですが・・・。 自分が考えた過程を書きます。 H=U+PV の式の両辺にΔを作用させて ΔH = ΔU + VΔP + PΔV ただし体積一定より ΔH = ΔU + VΔP ・・・(1) ここで状態方程式 PV = nRT の両辺にΔを作用させて VΔP + PΔV = RTΔn + nRΔT ただし体積、温度一定より VΔP = RTΔn よってこの式を(1)に代入して ΔH = ΔU + RTΔn ・・・(2) それで、ここからΔnを求めたいのですがよく分かりません。 自分としては U = 3/2nkT (k:ボルツマン定数)の両辺にΔを作用させてΔU = 3/2kTΔn これに値を代入して 1192 = 3/2*1.381*10^-23*300*Δn よってΔn=1.918*10^23になるのか、とも考えましたが、さすがにΔnが大きすぎます・・・。 どなたかお分かりになる方、ご教授ください。お願いいたします。

  • 熱力学に関する問題が解けません

    温度T,圧力p,体積v,内部エネルギーu,エントロピーs,エンタルピーhとし、体積変化だけが仕事に関係する。また系は準静的である 熱力学第一法則をもとに、独立な状態量をp,Tとおくことによって        Tds={(∂h/∂p)T -v}dp+(∂h/∂T)p dT となることを示せ。またこれをもとに        (∂h/∂p)T =v-T(∂v/∂T)p となることを示せ。 という問題です。自分で途中まで解いてみたんですが上式を導けません。誰か教えてください。お願いします。 下に途中まで解いた結果を書きます。 du=(∂u/∂P)T dP+(∂u/∂T)p dT ・・・(1) Tds=du+pdv であり(1)を代入 Tds=(∂u/∂p)T dp+(∂u/∂T)p dT+pdv ・・・(2) u=h-pv より(2)に代入 するところから間違っていると思います。 ちなみに( )のすぐ隣の文字は下文字です。( )T←下文字

  • 熱力学の変数

    気体の入った容器(体積V )の熱力学的性質が気体分子の総運動エネルギーUと体積Vで記述できないことを示したいのですが、どうすればよいでしょうか。 たとえば、一番単純な単原子分子理想気体である場合でも, エントロピーはU,V,Nで表されるので, 明らかにN(物質量)が足りないのは分かるのですが、 その証明がうまくいきません。 そもそも、単原子分子理想気体ではNが足りないのは良いけれども、 気体によっては物質量に無関係にエントロピーが定まるものがあるかもしれないということを排除するにはどうすれば良いのでしょうか。 また、気体の内部エネルギーがNによらないで、ある値U1(>0)が定まるとすると、 どんどん気体を抜いていきN=0の状態でもU1>0となり矛盾する ということも考えたのですが、なんか腑に落ちません。 どうすればよいでしょうか。教えてください。

  • 熱力学―エントロピーの導出について

    熱力学―エントロピーの導出について 熱力学が苦手でつまずいているのですが・・・ この問題の詳しい答えを教えてください。 「nモルの気体からなる閉じた系が圧力P,体積V,温度T,の状態にある。  この系の準静的変化を考えて,系のエントロピーを導出したい。  閉じた系内の気体の状態方程式は次式で表される。  {P+(an^2/v^2)}(V-nb)=nRT また,この気体の内部エネルギーUは次式で表される。  U=CT-a(n^2/V)  Rは気体定数,a,bおよびCはいずれも定数である。  温度T0,体積V0の基準状態でのエントロピーをS0とすれば,  温度T,体積V,のときのエントロピーSはどのように表されるか。」 ・・・面倒かもしれませんが詳しい導出の手順をお願いします。

  • E=mc2 のEは何エネルギー

    相対論でいわれるEは、エントロピー項を含んでいるのでしょうか。 Eは、次の熱力学的エネルギーのどれに相当するのでしょうか。 ・内部エネルギー U ・エンタルピー H ・ヘルムホルツ自由エネルギー F ・ギブズ自由エネルギー G 宜しくご教示お願い申し上げます。

  • 熱力学の問題で。

    温度Tの熱浴と熱平衡にある系のエネルギーUが、U=Vu(T)で与えられるときに準静的過程で体積変化させたときの熱の出入りとエントロピーの式の求め方がわかりません。Vは体積、uは単位体積あたりのエネルギーで温度の関数で系の圧力はp=u(T)/3です。 また、定積過程の場合も知りたいです。 計算したらRlogVになってしまいました。うーん・・・