- 追加された行はこの色です。
- 削除された行はこの色です。
''目次''
#contents
*2015年4月3日 [#a9f1eab6]
[[2015年4月3日>#m1605301]]
*2015年4月3日 [#m1605301]
**前回の宿題 [#b9a42568]
-(p,p'x)および(p,2p)における有効質量の取り扱いを整理する。
-g(k',k,E,ρ)を入力とするよう、pikoTDX を改良(k' と k は線形補間。ρはsuphod. k_F<0.05 は k_F=0.05 とする)。
**報告 [#tdda34ca]
***有効質量について [#pb5b442d]
-ノート&ref(EffMassP2P.pdf);
--結論: (p,pN)では不要
--参考資料&ref(effmass(20010810).pdf);&ref(RCNPinformal140130.pdf);
--参考資料&ref(effmass(20010810).pdf); &ref(RCNPinformal140130.pdf);
***コード改良 [#ed3b4b65]
-改良のメモ(ノート)&ref(inclusiveP2Pmod.pdf);
--完了(gのtable待ち)
--完了(σ~のtable待ち)
---従来の計算はオプションにて可能。
---σ~の形式(案)は下記。今のところ、k_Fを昇り順に与えるよう制限。
---今のところ、κ', κ, θは内挿せず、一番近い分点の値を用いる。外挿が必要な局面が発生すると計算中止。
---k_Fについては4点補間。外挿が必要な場合は、端点の値で代用。
--2核子散乱の性質について調査&ref(report150401(NNprof_p2p).pdf);
---off-shell性は問題にならない模様(朗報)。
***その他 [#g93cb74f]
-川瀬氏の14O(p,2p)解析にポテンシャルを提供。計算が進行中。上坂さんにコード提供を予告。
-早川氏の科研費(直線偏光ガンマ線で生成された中性子非対称性測定による新しいプローブ)が内定。基盤(B)で3年間。1600万円程度。
-渡邉君の学振書類の議論、お疲れ様でした。
**議論 [#f0261cdb]
***コードのデザイン [#cd91197d]
0. 運動量保存を満たすかどうかを毎回チェックする。
1. 順運動学と逆運動学の選択(*)
2. E1, θ1, φ1, θ2, φ2 を与える。
2a.範囲指定も可とする。
2b.プロファイルの読み込みも可能にする。
2c.θ2 の代わりに k_B を与えることも可能にする。
3. ポテンシャルは KD がデフォルト。
3a. 当面、粒子0 と、粒子1 (順運動学でEを動かさない場合のみ)は機番
から読み込みも可能とする。
3b. 将来的には、エネルギーも振ったテーブルを読み込めるようにする。*
4. NN 散乱については、以下のオプションを想定*
4a. 等方近似
4b. in-med 等方近似
4c. in-med NN 散乱の角分布[密度も読み込む]
4d. free NN 散乱の角分布
4e. g 行列(m の処理も正しく行う) [密度も読み込む]
5. Perey factor*
6. 1粒子波動関数の外部読み込み*
*はinclusive codeにも影響するもの。
確認しておきたいもの:
・核内核子の運動量の積分
・局所半古典近似の導入
・クーロンの影響
・純量子力学的散乱波(分業)
以下はpikoTDX0aで実装
0. 運動量保存を満たすかどうかを毎回チェックする。(済)
2. E1, θ1, φ1, θ2, φ2 を与える。
→ φ1 とφ2 は固定。θ2 の代わりに k_B を振る(今のまま)。
3. ポテンシャルは KD がデフォルト。(今のまま)
4. NN 散乱については、以下のオプションを想定 (済※)
4a. 等方近似
4b. in-med 等方近似
4c. in-med NN 散乱の角分布[密度も読み込む] ※挙動は未確認
5. Perey factor (済)
***σ~の読み込み形式(案) [#j7eeafd1]
READ(KIBNN,*) FKAPFMIN,FKAPFMAX,DKAPF
READ(KIBNN,*) FKAPIMIN,FKAPIMAX,DKAPI
READ(KIBNN,*) TH2NMIN,TH2NMAX,DTH2N
C
NKAPFMAX=(FKAPFMAX-FKAPFMIN)/DKAPF+1.01
NKAPIMAX=(FKAPIMAX-FKAPIMIN)/DKAPI+1.01
NTH2NMAX=(TH2NMAX-TH2NMIN)/DTH2N+1.01
C
IFKF=0
FKFOLD=-999.0D0
100 CONTINUE
READ(KIBNN,*, END=200) FKF
IF(FKF.LT.FKFOLD) STOP 3 ! condition on the input
FKFOLD=FKF
IFKF=IFKF+1
FKF2NMS(IFKF)=FKF
DO IKAPF=1,NKAPFMAX
DO IKAPI=1,NKAPIMAX
DO ITH2N=1,NTH2NMAX
READ(KIBNN,*) FKAPF,FKAPI,TH2N
1 ,DSIG2NIM(IKAPF,IKAPI,ITH2N,IFKF)
END DO
END DO
END DO
GO TO 100
200 CONTINUE
***報告ページ(二位の会?)の必要性について [#t36e4edd]
***三位の会の日程について [#i398d1d2]
-明さんが、事実上金曜日しか参加不可能。
**次回までの宿題 [#e322b6a2]
-次回予定: 4月*日*時
