Progrep2015 の履歴(No.4)

履歴一覧
ソースを表示
Progrep2015 は削除されています。
- 1 (2015-04-02 (木) 12:09:11)
- 2 (2015-04-02 (木) 12:48:45)
- 3 (2015-04-02 (木) 20:52:45)
- 4 (2015-04-03 (金) 18:40:51)

2015年4月3日

2015年4月3日†

↑

前回の宿題†

(p,p'x)および(p,2p)における有効質量の取り扱いを整理する。
g(k',k,E,ρ)を入力とするよう、pikoTDX を改良(k' と k は線形補間。ρはsuphod. k_F<0.05 は k_F=0.05 とする)。

↑

報告†

↑

有効質量について†

ノート&ref(): File not found: "EffMassP2P.pdf" at page "Progrep2015";
- 結論: (p,pN)では不要
- 参考資料&ref(): File not found: "effmass(20010810).pdf" at page "Progrep2015"; &ref(): File not found: "RCNPinformal140130.pdf" at page "Progrep2015";

↑

コード改良†

改良のメモ(ノート)&ref(): File not found: "inclusiveP2Pmod.pdf" at page "Progrep2015";
- 完了(σ~のtable待ち)
  - 従来の計算はオプションにて可能。
  - σ~の形式(案)は下記。今のところ、k_Fを昇り順に与えるよう制限。
  - 今のところ、κ', κ, θは内挿せず、一番近い分点の値を用いる。外挿が必要な局面が発生すると計算中止。
  - k_Fについては4点補間。外挿が必要な場合は、端点の値で代用。
- 2核子散乱の性質について調査&ref(): File not found: "report150401(NNprof_p2p).pdf" at page "Progrep2015";
  - off-shell性は問題にならない模様(朗報)。

↑

その他†

川瀬氏の14O(p,2p)解析にポテンシャルを提供。計算が進行中。上坂さんにコード提供を予告。
早川氏の科研費(直線偏光ガンマ線で生成された中性子非対称性測定による新しいプローブ)が内定。基盤(B)で3年間。1600万円程度。
渡邉君の学振書類の議論、お疲れ様でした。

↑

議論†

↑

コードのデザイン†

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 (済)

↑

σ~の読み込み形式(案)†

      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

↑

報告ページ(二位の会?)の必要性について†

↑

三位の会の日程について†

明さんが、事実上金曜日しか参加不可能。

↑

次回までの宿題†

wiki(反応会)の整備
論文の執筆を進める。
次回予定: 4月9日13時