Progrep2015 の変更点
Top/Progrep2015
- 追加された行はこの色です。
- 削除された行はこの色です。
- Progrep2015 は削除されています。
- Progrep2015 の差分を削除
''目次'' [[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); ***コード改良 [#ed3b4b65] -改良のメモ(ノート)&ref(inclusiveP2Pmod.pdf); --完了(σ~の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 を与えることも可能にする。 ※粒子1がpかnかは指定可能とする(directかexchangeかを区別する必要あり)。 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] -wiki(反応会)の整備 -PMD論文の執筆を進める。 -次回予定: 4月9日13時