次のHPにつきまして、下記をご教示頂きましたら幸いです。
http://www.rs.kagu.tus.ac.jp/beachey/NP/NP_2006_ …
P53の
a = R*(1 + ε); 長径
b = R*(1 -ε/2); 短径
の場合、楕円体の表面積、表面エネルギー、クーロンエネルギーは
どのようにしたら導かれるのでしょうか?
P54のε=0.3の場合、ΔE=6MeVは
どのようにしたら導かれるのでしょうか?
A 回答 (3件)
- 最新から表示
- 回答順に表示
No.3
- 回答日時:
>やはり、Legendre関数での変形展開をしないと、正確な表面エネルギーは
>求めることができないことがわかりました。
>Legendre関数で変形展開をして、正確な表面エネルギーを求める資料はないでしょうか?
何処からLegendre関数が出てきて、Legendre関数を使えば正確な表面エネルギーを求める事ができると思うのか、また、なんでこの楕円体の模型が「正確でない」とするのか全く分かりませんが、
Legendre関数に関する基本的な事はそこら辺の物理数学の本にも載っているんじゃないですかね。そういうのを見ながら式をこねくりまわしてみては。
No.2
- 回答日時:
ん、すいません。
お礼・補足がある事に気付いてませんでした。>そうなのですか。実際に計算をされたのですね。計算する際の公式を
>ご教示頂きましたら幸いです。
あぁ、それは申し訳ないです。
この質問とは全く関係ない理由で、図書館で原子核関連の本を読み漁っていたのですが、偶然、表面エネルギーとクーロンエネルギーの解析的な計算結果が載っているのを見つけたので、基本的にはそれをもとに計算しました。(その時、たくさんの本を読んでいたので、何て本だったかは・・・)
少なくとも、表面積の方は、解析的な計算結果の導出からεの2次くらいまでの計算までちゃんとやったはずです。この導出に関しては、例えば、
http://next1.cc.it-hiroshima.ac.jp/MULTIMEDIA/ca …
が参考になるでしょう。例題8.4において、
r=(bsinucosv,bsinusinv,acosu)
としたのが、今考えている楕円体ですね。
クーロンエネルギーの方については、解析的な解の導出はやっていないのは確かなのですが、εの2次までの計算をやったかどうかはもう忘れました。#1を見る限り、どうも「確かめた」という口ぶりなので、多分、確かめたのだろうとは思います。
確かめていませんが、クーロンエネルギーは
ρ^2 ∫dVdV'/|x-x'|
を計算すればいいはずです。ρは原子核の電荷密度(=3Ze/4πR^3)、積分範囲はxについてもx'についても(変形した)原子核の内部です。
積分変数をx,y,zととるのは大変ですので、上手く置換してください。まぁ、解析解も求まるらしいですが、必要なのはεの低次の項だけですし、適当なところでテイラー展開するのが楽でしょう
なお、#1へのお礼のプログラムに関しては、
>m = a^2*(b^2^c^2/(b^2*(a^2 - c^2)));
これを、
m = a^2*((b^2-c^2)/(b^2*(a^2 - c^2)));
あるいは、もっと単純に
m = 0;
などと書き換えればいいのではないでしょうかね。
この回答への補足
お返事ありがとうございます。
>少なくとも、表面積の方は、解析的な計算結果の導出からεの2次くらいまでの計算>までちゃんとやったはずです。
計算できました。ありがとうございます。
おっしゃる通り、2/5εまでは計算できましたが、-4/105εは得られませんでした。
やはり、Legendre関数での変形展開をしないと、正確な表面エネルギーは
求めることができないことがわかりました。
Legendre関数で変形展開をして、正確な表面エネルギーを求める資料はないでしょうか?
No.1
- 回答日時:
その式は微妙に嘘っぽいです。
a=R(1+ε)
b=R/√(1+ε)
とした時に、表面積やクーロンエネルギーが与えられているようになりますので、これで計算してみてください。
※bについてはεの1次までなら与えられているものと一緒ですが、εの2次の項で違いが出てきます。
※上のようにするのは、原子核の体積が変化しないようにするためです。
>P54のε=0.3の場合、ΔE=6MeVは
>どのようにしたら導かれるのでしょうか?
一応、a_cは実験的に知られている値ですので、その値を使ったんじゃないですかね。(手元の本にはa_c≒0.699MeVとあります)
この回答への補足
お返事ありがとうございます。
>※bについてはεの1次までなら与えられている
>ものと一緒ですが、εの2次の項で違いが出てきます。
そうなのですか。実際に計算をされたのですね。計算する際の公式を
ご教示頂きましたら幸いです。
下記HPは、おっしゃる通り下記式を使用してます。
a=R(1+ε)
b=R/√(1+ε)
http://yuri.phys.tohoku.ac.jp/~hagino/lectures/n …
追伸
お返事が大変遅れましたことをお詫び申し上げます。
自分で下記HPを参考にしまして、mathematicaで計算してみました。
しかし、答え4πR^2(1+2/5ε^2-4/105ε^3+、、、、)が導けません。
後の処理をご教示頂きましたら幸いです。
http://documents.wolfram.com/v5/Demos/FormulaGal …
http://www.asahi-net.or.jp/~jb2y-bk/NaturalSci/m …
a = R*(1 + e);
b = R/Sqrt[1 + e];
c = b;
m = a^2*(b^2^c^2/(b^2*(a^2 - c^2)));
q = ArcSin[Sqrt[1 - c^2/a^2]];
S = 2*Pi*(c^2 + ((b*c^2)/Sqrt[a^2 - c^2])*EllipticE[q, m] + b*Sqrt[a^2 - c^2]*EllipticF[q, m]);
Print[FullSimplify[S]];
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 物理学 中性子の結合エネルギーの意味を教えてください。 2 2022/08/19 10:31
- 物理学 中心を同じに点に持つ半径aの導体球(導体1)、内半径b、外半径Cの導体球殻(導体2)があるとして、導 1 2023/08/12 23:36
- 工学 至急お願いします。 誘電体と接する導体表面に面密度のσ正の電荷を一様に与えると、境界面には応力が発生 1 2022/07/31 02:27
- 地球科学 地球内核最中心部における核融合反応の可能性について質問です。 現在地球内核の熱源は地球始原時のエネル 2 2022/05/08 15:57
- 物理学 電磁気の問題で質問です。 2 2022/07/16 17:41
- Excel(エクセル) Excelで睡眠時間をもとに判定したい 6 2022/08/19 13:47
- 統計学 統計学 二項分布の正規近似について 2 2023/02/10 11:58
- 物理学 共鳴領域での中性子の核分裂スペクトル 2 2023/03/29 11:39
- 物理学 電磁気学の問題です。 真空中に置かれた半径 a[m ]の導体球の周りを 、 内半径b[m ] 、 外 2 2023/08/10 19:02
- 物理学 万有引力と重力の位置エネルギーについて 例えば、地球の表面から真上に質量mの球を初速v₀で投げた時の 7 2022/04/18 23:15
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
振動数vをもつ一つの振動子につ...
-
位置エネルギー U
-
普段の生活の中での位置エネルギー
-
「U = mgh」の「U」は何の略な...
-
一分子の基底状態と励起状態の...
-
縮退をわかりやすくお願いします
-
ラマン分光法・ラマンシフトに...
-
超高温 vs 超低温(絶対零度)...
-
フェルミディラック分布関数の...
-
エネルギースペクトルとは?
-
波数(k)を用いた空間座標表示を...
-
永久機関が作れない理由
-
人間のジャンプ時の衝撃値は?
-
RLC回路
-
水路が分岐た場合の水圧について
-
磁石のエネルギー
-
出力とエネルギーの違いとは
-
消音(ノイズキャンセリング)す...
-
磁石の仕事
-
回路のエネルギー収支について ...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
E=mc²を誰にもわかるように説...
-
位置エネルギー U
-
車に詳しい方
-
「U = mgh」の「U」は何の略な...
-
高校物理の力学の質問
-
力学的エネルギーの保存でレー...
-
水路が分岐た場合の水圧について
-
再生可能エネルギーは、再生で...
-
人間のジャンプ時の衝撃値は?
-
スイッチSを閉じて十分に時間が...
-
何で暇だとエロいことを考えて...
-
どんなに能力高い人でも、多忙...
-
縮退をわかりやすくお願いします
-
共振が起こるとなぜエネルギー...
-
干渉して打ち消しあった光の波...
-
泡が壁面にくっつくのは…
-
もう、何も頑張りたくないです...
-
消音(ノイズキャンセリング)す...
-
物理のエッセンスで「失われた...
-
人体からの発熱量の計算方法
おすすめ情報