Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
橋本 慎太郎; 永井 泰樹*
Journal of the Physical Society of Japan, 92(12), p.124202_1 - 124202_11, 2023/12
Moの娘核種であるテクネチウム99m(Tc)は6時間の半減期で線を出す放射性核種(RI)であり、SPECTなどの診断用RIとして世界中で用いられている。親核種のMoは原子炉で生成されているが、原子炉施設の老朽化等で安定供給が危ぶまれており、医療用加速器で濃縮Mo試料に陽子を照射し、Mo()Tc反応を通じて生成する方法が提案されている。しかし、試料に微量に含まれるMo以外のMo同位体から、Tc以外のTcが生成される問題がある。そこで本研究では、高純度のTcを生成できるMo組成を調べるために、天然組成のMoの照射実験データに基づいて、全てのTc同位体生成の反応率を決定し、任意の組成をもつMo試料のTc同位体生成量をPHITSを用いて評価する手法を開発した。本手法の有効性を、濃縮Moで生成された各Tcの実測値との比較によって検証するとともに、高純度のTc生成にはMoの各不純物含有量の低減が必須であることを明らかにした。本成果は、照射時間等の施設固有の制約を考慮しながら、高純度のTcを生成する際に必要なMo原料の組成決定に重要な役割を果たす。
久峩 尚也*; 椎葉 拓郎*; 佐藤 達彦; 橋本 慎太郎; 黒岩 靖淳*
Journal of Nuclear Science and Technology, 10 Pages, 2023/00
本論文では、放射線挙動解析コードPHITSの高エネルギー光子線治療に対する精度検証結果を報告する。具体的には、Clinac 21EX電子ライナックにより発生させた光子線治療場における線量深さ分布,ビームプロファイル,アウトプットファクターを測定し、PHITS及び同様のモンテカルロ計算コードであるEGSnrcの結果と比較した。その結果、PHITS計算と測定は医学物理分野における許容レベル以内で、PHITS計算とEGSnrc計算は約1%の範囲内で一致することが分かった。これらの結果から、PHITSが高エネルギー光子線治療の医学物理計算に応用可能であることが確認できた。
佐藤 達彦; 岩元 洋介; 橋本 慎太郎; 小川 達彦; 古田 琢哉; 安部 晋一郎; 甲斐 健師; 松谷 悠佑; 松田 規宏; 平田 悠歩; et al.
Journal of Nuclear Science and Technology, 9 Pages, 2023/00
放射線挙動解析コードPHITSは、モンテカルロ法に基づいてほぼ全ての放射線の挙動を解析することができる。その最新版であるPHITS version 3.31を開発し公開した。最新版では、高エネルギー核データに対する親和性や飛跡構造解析アルゴリズムなどが改良されている。また、PHIG-3DやRT-PHITSなど、パッケージに組み込まれた外部ソフトウェアも充実している。本論文では、2017年にリリースされたPHITS3.02以降に導入された新しい機能について説明する。
古田 琢哉; 古場 裕介*; 橋本 慎太郎; Chang, W.*; 米内 俊祐*; 松本 真之介*; 石川 諒尚*; 佐藤 達彦
Physics in Medicine & Biology, 67(14), p.145002_1 - 145002_15, 2022/07
被引用回数:0 パーセンタイル:49.94(Engineering, Biomedical)炭素線治療は従来の放射線治療よりも腫瘍部への線量集中性に関する優位性を持つが、二次的ながんの発生原因となり得る正常組織への照射を完全に無くすことは困難である。そのため、発がんリスクを照射炭素ビームの核反応によって生成される二次粒子による線量まで含めて評価するには、計算シミュレーション解析が有効となる。本研究では、PHITSコードを中核とした炭素線治療の線量再構築システムを開発した。このシステムでは、治療計画を記録したDICOMデータから自動でPHITSの入力ファイルを作成し、PHITSシミュレーションの実行によって腫瘍および周辺正常組織での線量分布を計算する。PHITSの様々な機能を利用することで、粒子毎の線量寄与や二次粒子の発生場所の特定など、詳細な解析が実施可能である。開発したシステムの妥当性は、水中での線量分布の実験結果や人体等価ファントムへの治療計画との比較により確認した。今後、本システムは量子科学技術研究開発機構において、過去の治療データを用いたシミュレーション解析による遡及的研究に利用される予定である。
Walter, H.*; Colonna, M.*; Cozma, D.*; Danielewicz, P.*; Ko, C. M.*; Kumar, R.*; 小野 章*; Tsang, M. Y. B*; Xu, J.*; Zhang, Y.-X.*; et al.
Progress in Particle and Nuclear Physics, 125, p.103962_1 - 103962_90, 2022/07
被引用回数:32 パーセンタイル:96.94(Physics, Nuclear)原子核-原子核衝突や原子核の状態方程式の研究において、反応計算モデルは重要なツールとなり、世界中で開発が進んでいる。本論文は、原子力機構のJQMD-2.0を含め、現在開発中の複数のコード開発者の協力により、これらコードを同じ条件で比較することで共通点や差異を明らかにしたプロジェクトTransport Model Evaluation Project (TMEP)を総括したものである。参加したコードはBoltzmann-Uehling-Uhlenbeck(BUU)法に基づく13のコードと、Quantum Molecular Dynamics (QMD)法に基づく12のコードであった。プロジェクトでは、Au原子核同士を衝突させてその終状態を観測する現実的な計算や、一辺が640nmの箱に核子を詰めて時間発展させる仮想的な計算を行った。その結果、BUU法コードとQMD法コードは計算原理が異なるため、計算の設定に関係なく系統的な差異が生じることが明らかになった。その一方で、同じ方法を採用するコード間の比較では、時間発展を細かく計算することでコード間の差は埋まっていき、一定の収束値を持つことが示された。この結果は今後開発される同分野のコードのベンチマークデータとして有用なものであるだけでなく、原子核基礎物理学の実験や理論研究の標準的な指針としても役に立つことが期待される。
岩元 洋介; 橋本 慎太郎; 佐藤 達彦; 松田 規宏; 国枝 賢; elik, Y.*; 古立 直也*; 仁井田 浩二*
Journal of Nuclear Science and Technology, 59(5), p.665 - 675, 2022/05
被引用回数:3 パーセンタイル:73.26(Nuclear Science & Technology)粒子・重イオン輸送計算コードPHITSは、加速器施設の放射線遮蔽設計に広く利用されている。本研究では、NEA/NSCの遮蔽積分実験ベンチマーク・アーカイブデータベース(SINBAD)の中から、陽子加速器施設の中性子遮蔽実験に関するデータを用いて、最新のPHITSバージョン3.24の精度検証を行った。検証では、中性子源に用いる陽子入射核反応が(1)43及び68MeV陽子-リチウム、(2)52MeV陽子-炭素、(3)590MeV陽子-鉛、(4)500MeV陽子-タングステン、及び(5)800MeV陽子-タンタルの条件で、遮蔽体における中性子エネルギースペクトル及び反応率を計算した。その結果、PHITSの標準設定である核反応モデルINCL4.6を用いた計算結果は(1)及び(2)のケースの100MeV以下で実験値を再現しなかったが、高エネルギー核データライブラリJENDL-4.0/HEを用いた場合は全ケースで実験値を良く再現した。本検証により、陽子加速器施設における中性子遮蔽設計において、JENDL-4.0/HEを用いたPHITSの計算が有効であることが明らかになった。
永井 泰樹*; 川端 方子*; 橋本 慎太郎; 塚田 和明; 橋本 和幸*; 本石 章司*; 佐伯 秀也*; 本村 新*; 湊 太志; 伊藤 正俊*
Journal of the Physical Society of Japan, 91(4), p.044201_1 - 044201_10, 2022/04
被引用回数:0 パーセンタイル:48.58(Physics, Multidisciplinary)近年、神経内分泌腫瘍を治療するための医療用RIとしてLuが注目されており、加速器施設で重陽子を濃縮Yb試料に照射し、高純度のLuを製造する方法が検討されている。ただし、Yb試料には様々な同位体が微量に含まれており、Lu以外のLu同位体が不純物として生成される。医療用として利用するLuは一定の純度が求められるため、定量的にLuの不純物を評価する手法が求められていた。本研究では、実験値を基に関与する全てのYb()Lu反応の断面積を決定し、粒子輸送計算コードPHITSと組み合わせることで、Yb試料が任意の組成比をもつ場合の各Lu同位体の生成量を推定する新しい手法を開発した。他に、天然組成のYb試料へ25MeV重陽子を照射した実験も行い、本手法の有効性を検証した。また、市販の濃縮Yb試料に照射した条件で計算を行い、重陽子のエネルギーを15MeVとすることで、純度99%以上のLu生成を達成できることを示した。開発した手法は、加速器を用いて高純度の医療用Luを製造する際、必要な濃縮Yb試料の同位体組成を議論する上で重要な役割を果たすものである。
川端 方子*; 本石 章司*; 太田 朗生*; 本村 新*; 佐伯 秀也*; 塚田 和明; 橋本 慎太郎; 岩本 信之; 永井 泰樹*; 橋本 和幸*
Journal of Radioanalytical and Nuclear Chemistry, 330(3), p.913 - 922, 2021/12
被引用回数:4 パーセンタイル:73.26(Chemistry, Analytical)加速器中性子を利用し、将来有望な医療RIであるCu及びCuの大量合成に関して、実験に基づく生成量評価と臨床に必要な生成量のPHITSコード使用による評価、そして大量のターゲット物質からのCu同位体の分離精製手法開発を行った。実験では55.4gという大量の、天然亜鉛試料を利用し、加速器中性子の照射による生成率について定量的な評価を行った。また、同時に臨床的に必要な量のCu及びCuの製造の可能性について、それぞれ100gの濃縮ZnおよびZn試料に対して、40MeV、2mAの重陽子照射により発生する中性子を照射する数値シミュレーションを行い、その生成量を推定した。また、大量の天然亜鉛試料からCu及びCuを分離するために、熱分離と樹脂分離を組み合わせた分離方法を開発し、73%の分離効率と97%の亜鉛回収率が得られることを実証した。これらの結果、このような加速器中性子による生成と本分離手法の組み合わせで、臨床応用のために必要なCu及びCuの大規模生産を提供することができることがわかった。
佐藤 達彦; 橋本 慎太郎; 稲庭 拓*; 高田 健太*; 熊田 博明*
International Journal of Radiation Biology, 97(10), p.1450 - 1460, 2021/10
被引用回数:6 パーセンタイル:78.4(Biology)確率論的マイクロドジメトリ運動学(SMK)モデルは、粒子線治療(CRT)やホウ素中性子捕捉療法(BNCT)の生物学的線量を最も精度よく推定できるモデルの一つである。しかし、SMKモデルは、計算時間が掛かるため治療計画システムに導入することは困難であった。そこで、テイラー展開や高速フーリエ変換を用いて単純化したSMKモデルを開発し、それらを原子力機構が中心となって開発している汎用放射線挙動解析コードPHITSに組み込んだ。これにより、SMKモデルに基づく生物学的線量を簡単に短時間で導出可能となり、PHITSを用いたより高度な治療計画システムの開発が可能となった。
Chang, W.*; 古場 裕介*; 古田 琢哉; 米内 俊祐*; 橋本 慎太郎; 松本 真之介*; 佐藤 達彦
Journal of Radiation Research (Internet), 62(5), p.846 - 855, 2021/09
被引用回数:2 パーセンタイル:27.7(Biology)炭素線治療の治療計画をモンテカルロシミュレーションによって再評価するためのツール開発の一環として二つの重要な手法を考案した。一つは治療計画装置に含まれる校正済みのCT-水阻止能表を反映しつつ、患者CT画像から物質識別を行う手法である。もう一つの手法は、粒子およびエネルギー毎に生体物質と水との阻止能比を考慮し、水等価線量を導出することを目的としたものである。これらの手法の有効性を確認するため、生体物質および水で構成される均質および不均質ファントムに対し、SOBPサイズ8cmの炭素線を400MeV/uで照射するシミュレーションをPHITSで計算し、得られた線量深さ分布を従来の治療計画装置による結果と比較した。その結果、従来の治療計画装置で採用されている生体物質を全て水に置換する手法では、異なる物質において二次粒子の発生率が一次粒子の阻止能に単純に比例するため、二次粒子による線量寄与を評価する上で不適切であることが分かった。一方、各物質の二次粒子の発生率を適切に考慮した新しい手法は炭素線治療の再評価において、重要な役割を果たすことが期待できることが分かった。
橋本 慎太郎; 佐藤 達彦
EPJ Web of Conferences, 239, p.03015_1 - 03015_4, 2020/09
被引用回数:0 パーセンタイル:0.1モンテカルロ法に基づいた粒子輸送計算コードは、加速器施設の遮蔽計算等の目的で盛んに利用されている。計算結果の信頼性を評価するためには、計算で用いた試行回数によって決まる統計的不確かさだけでなく、計算で使用される全断面積モデル等の不確かさに起因する系統的不確かさも評価する必要がある。我々は、全断面積の実験値を参照しながら最小自乗法に基づいたKALMANコードを利用することで、内部パラメータに不確かさの幅をもつ全断面積モデルを開発した。新規モデルで計算した全断面積の不確かさは実験値の誤差と同程度であり、不確かさの幅の範囲で内部パラメータを変動させながら輸送計算を実行することで、計算結果における系統的不確かさを適切に評価できることが分かった。
岩元 洋介; 明午 伸一郎; 橋本 慎太郎
Journal of Nuclear Materials, 538, p.152261_1 - 152261_9, 2020/09
被引用回数:13 パーセンタイル:90.11(Materials Science, Multidisciplinary)様々な放射線環境下の材料の照射損傷の指標として、はじき出し原子数(dpa)が用いられている。これまでの照射損傷モデルではNorgertt-Robinson-Torrens(NRT)モデルが標準とされていたが、最新の非熱的再結合補正(athermal-recombination-correction, arc)を考慮したモデルは照射損傷過程をより正確に再現できるため、arcモデルを用いた材料の照射損傷の評価が今後標準となることが期待されている。本研究では、様々な材料に対するarcモデルを粒子・重イオン輸送計算コードPHITSに組み込み、核分裂、核融合及び加速器施設を想定した幅広いエネルギー範囲における再スケーリング係数(NRT-dpa/arc-dpa)を計算した。計算の結果、核施設の照射環境を想定したエネルギー分布を持つ中性子照射、及びエネルギー600MeV-50GeVを持つ陽子照射の場合、各材料の再スケーリング係数は、照射条件によらず、各材料でほぼ同じ値であることがわかった。本成果は、幅広いエネルギー範囲における様々な核施設の従来の放射線損傷評価において、NRT-dpaを再スケーリングするための大変有益な情報となる。
塚田 和明; 永井 泰樹*; 橋本 慎太郎; 湊 太志; 川端 方子*; 初川 雄一*; 橋本 和幸*; 渡辺 智*; 佐伯 秀也*; 本石 章司*
Journal of the Physical Society of Japan, 89(3), p.034201_1 - 034201_7, 2020/03
被引用回数:2 パーセンタイル:22.58(Physics, Multidisciplinary)ポリエチレン遮へい中のZnOに、50MeV重陽子とBeによる()反応で生成した中性子を照射することで、Ga, Ga, Zn、並びにCuの特異な生成を実験的に確認した。特に、ポリエチレン遮へい内で得られた収率は、遮へいなしの実験と比較して、約20倍の収量を示した。一方、鉛遮へい内の金属Zn試料の照射におけるGa, Ga, Zn、並びにCuの収量と、ZnO及び金属Zn試料の照射におけるCu, Ni及びZnの収量は、遮へいによる影響はほとんど受けていない。この実験結果は、遮へい条件を調整することで、中性子反応に限らず陽子反応を含む多様で大量の放射性同位元素を、一度の照射で同時に合成できるという加速器中性子の注目すべき特性を示すものである。また、PHITSコードを利用した生成量予測を試み、本実験結果と比較することで、本生成量の特異性について評価した。
中山 梓介; 岩本 修; 岩本 信之; 橋本 慎太郎
核データニュース(インターネット), (123), p.53 - 59, 2019/06
2019年5月19日から24日まで中国ナショナルコンベンションセンターにて、2019年科学と技術のための核データ国際会議(ND2019)が開催された。核データ国際会議は3年に一度開催される核データ研究分野の中で規模の最も大きな国際会議である。本稿ではND2019の会議報告の一部として、本会議で行われた核データ評価及び理論に関する発表の概要を記した。
橋本 慎太郎; 佐藤 達彦
Journal of Nuclear Science and Technology, 56(4), p.345 - 354, 2019/04
被引用回数:4 パーセンタイル:49.54(Nuclear Science & Technology)モンテカルロ法に基づいた粒子輸送シミュレーションは、加速器施設の遮蔽計算等の様々な目的で利用されている。モンテカルロ法による計算結果の信頼性を定量的に評価するためには、試行回数によって決まる統計的不確かさに加えて、計算で使われる反応断面積や計算の入力情報となる遮蔽材の密度がもつ誤差の影響を系統的不確かさとして求める必要がある。本研究の評価方法は分散分析に基づいており、これを鉄やコンクリートを遮蔽材とする中性子の遮蔽計算の解析に適用し、統計的及び系統的不確かさの両方を評価可能であることを示した。その際、入力情報の値を誤差の範囲で変動させるランダム条件法の他、三条件(中央値,上限値及び下限値)のみ変動させる三条件法を提案した。ランダム条件法は計算体系に関わらず適切に系統的不確かさが評価できるものの長い計算時間を必要とするのに対し、入力情報の誤差の影響が複雑な場合を除いて、三条件法は計算時間を抑えてランダム条件法と同じ評価結果を与えることがわかった。さらに、試行回数を増加させた場合の収束状況を判断できる新しい基準値を示し、必要最小限の計算時間で収束した評価結果が得られることを明らかにした。
小川 達彦; 佐藤 達彦; 橋本 慎太郎; 仁井田 浩二*
Physical Review C, 98(2), p.024611_1 - 024611_15, 2018/08
被引用回数:1 パーセンタイル:11.79(Physics, Nuclear)GeVからTeV級の原子核-原子核衝突による粒子生成は、重イオン加速器の安全設計や宇宙線被ばく評価などで非常に重要で、様々な反応モデルが考案されてきた。しかし、従来のモデルは核子間に働く相互作用が座標系に依存しており、運動する入射核と静止しているターゲット核が同じ座標系に移行すると本来は安定な核が壊れる問題があった。ここでは、人為的なバイアスを計算アルゴリズムに導入することで安定性を補っていたが、その副作用として残留核の質量や粒子の放出量が過小評価される問題があった。そこで、核子間の相互作用を座標系に依存しないよう記述した反応モデルJAMQMDを構築した。このモデルでは、原子核-原子核衝突反応を再現した際に核は座標系によらず安定を保ち、重陽子を含む多様な残留核の形成と二次粒子の生成も実験値をよく再現することも確認できた。さらに、3年前に開発された原子核-原子核衝突反応モデルJQMD Ver.2は3AGeV以下の入射エネルギーでしか適用できなかったが、JAMQMDでは入射エネルギーの制限をなくすことができた。そのため、JAMQMDは、放射線防護研究のみならず基礎研究の解析や計画などにも有益なモデルといえる。
佐藤 達彦; 岩元 洋介; 橋本 慎太郎; 小川 達彦; 古田 琢哉; 安部 晋一郎; 甲斐 健師; Tsai, P.-E.; 松田 規宏; 岩瀬 広*; et al.
Journal of Nuclear Science and Technology, 55(6), p.684 - 690, 2018/06
被引用回数:652 パーセンタイル:100(Nuclear Science & Technology)粒子・重イオン輸送計算コードPHITS Version 3.02を開発して公開した。核反応モデルや原子反応モデルを改良することにより、適応エネルギー範囲や精度を向上させた。また、計算効率を向上させる新しいタリーや、放射性同位元素を線源とする機能、医学物理分野にPHITSを応用する際に有用となるソフトウェアなど、様々なユーザーサポート機能を開発して、その利便性を大幅に向上させた。これらの改良により、PHITSは、加速器設計、放射線遮へい及び防護、医学物理、宇宙線研究など、幅広い分野で3000名以上のユーザーに利用されている。本論文では、2013年に公開したPHITS2.52以降に追加された新機能を中心に紹介する。
橋本 慎太郎; 佐藤 達彦; 岩元 洋介; 小川 達彦; 古田 琢哉; 安部 晋一郎; 仁井田 浩二*
核データニュース(インターネット), (120), p.26 - 34, 2018/06
放射線の挙動を模擬できる粒子・重イオン輸送計算コードPHITSは、加速器の遮へい設計や宇宙線・地球科学、放射線防護研究等の広範な分野で利用されている。PHITSでは、放射線が引き起こす物理現象を輸送過程と衝突過程の2つに分けて記述しており、さらに核反応が含まれる衝突過程を「核反応の発生」、「前平衡過程」、「複合核過程」の3段階的に分けて計算することで模擬している。これらの計算では断面積モデルKurotamaや核反応モデルINCL4.6やJQMDといった数多くの物理モデルが使用されており、PHITSの信頼性を高めるために、我々は各モデルの高度化を行ってきた。本稿は、日本原子力学会2018年春の年会の企画セッション「我が国における核データ計算コード開発の現状と将来ビジョン」における報告を踏まえ、PHITSにおける核反応モデルの役割を解説するとともに、近年我々が行ってきた様々なモデルの高度化や今後の展開についてまとめたものである。
小川 達彦; 橋本 慎太郎; 佐藤 達彦; 仁井田 浩二*; 釜江 常好*
Proceedings of 3rd International Conference on Particle Physics and Astrophysics (ICPPA 2017) (Internet), p.391 - 398, 2018/04
被引用回数:0 パーセンタイル:0.11宇宙線重イオンに起因する高エネルギー線生成の反応において、重イオンは高速で地表に向かって進み、その破砕片から生じる即発線はドップラー効果により数100MeV程度の高いエネルギーを持ち、地表でも観測できる。そのため、近年、宇宙放射線物理においてこの反応への高い関心が寄せられている。しかし、この高エネルギー線の生成を計算で再現するためには、破砕反応を個々のイベントに分けて計算し、核破砕片の励起エネルギーや角運動量を計算し、さらに原子核構造データをもとにした脱励起計算をする必要がある。粒子重イオン輸送計算コードPHITSは、他のコードにはないこのような過程を模擬できるモデル等を含んでいる。例えば、近年PHITSは破砕片生成を計算するJAERI Quantum Molecular Dynamics (JQMD)モデルを改良し、破砕片の電荷や質量分布をより正確に再現することが可能となった。また、近年実装されたモデルEBITEMは、励起エネルギーや角運動量の情報に基づき、蒸発過程後の脱励起を再現できる。以上のように、最新のPHITSは、宇宙線重イオンに起因する線生成に対して優れた再現性を発揮できる。本研究は、PHITSを用いた宇宙線重イオンに起因する反応の統合的なシミュレーションにより、重イオンに起因する高エネルギー線生成の再現について貴重なアプローチを図ったものである。
湊 太志; 塚田 和明; 佐藤 望*; 渡辺 智*; 佐伯 秀也*; 川端 方子*; 橋本 慎太郎; 永井 泰樹*
Journal of the Physical Society of Japan, 86(11), p.114803_1 - 114803_6, 2017/11
被引用回数:8 パーセンタイル:57.21(Physics, Multidisciplinary)核診断に使われるTcの親核であるMoの生成量の測定と計算を行った。MoはMo(,)Mo反応によって作られ、中性子源はC(,)反応を用いた。中性子のエネルギー範囲は熱領域から約40MeVまでである。測定されたMoの生成量は、最新のC(,)反応の実験データおよびJENDL評価済み断面積を用いて予測された数値計算結果と一致していることが分かった。次に、経済的に適したMoの生成法を模索するため、Mo生成の条件を変えた新しい系統的な数値計算を実施した。考慮した条件は、MoOサンプルの質量、炭素標的とサンプル間の距離、重陽子ビームの半径、照射時間である。得られたMoの生成量より、一つの加速器でMoの日本の需要量の約30%を担うことができることが分かった。また、1日2回MoからTcを溶出することで、Moの需要量の約50%まで引き上げることが可能であることが分かった。