Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
小川 達彦; 平田 悠歩; 松谷 悠佑; 甲斐 健師
Computer Physics Communications, 316, p.109758_1 - 109758_15, 2025/11
被引用回数：0 パーセンタイル：0.00任意の物質・任意のイオン照射において、電離や励起等を個々にシミュレートしながら放射線挙動を解析できる飛跡構造解析モデルITSART Ver.2を開発した。一般のイオン輸送アルゴリズムは、数千の原子反応をまとめて1イベントとして計算するため、計算が高速ではあるが電離など原子レベルの反応を識別することはできない。また従来の飛跡構造解析モデルは水や生体分子のみにしか使えず、入射粒子も陽子線や炭素線などの医療目的の粒子線に限られていたことから、本モデルの汎用性は画期的である。本モデルを汎用的飛跡構造解析モデルとして開発するにあたり、多くの工夫が必要であることを本研究では突き止め、それらを実現した。具体的には、電離で発生する二次電子のエネルギー・角度分布、標的核への運動量移行、原子脱励起モデルとの接続であり、これらを実装したことで多くのベンチマークデータを再現することができた。特に重い元素標的に照射した場合のビーム軸外縁における線量や、二次電子スペクトルに見られる原子固有のオージェ電子ピークが再現できたことは、標的物質が水だけであった従来研究に比べて特に顕著な進歩である。こうして開発されたITSART Ver.2はPHITSの次期リリースから実装される予定であり、炉材料や半導体などの放射線影響を解析するための強力なツールとして期待される。
松永 佳大*; 相澤 広記*; 浅森 浩一; 小川 大輝*; 宇津木 充*; 吉村 令慧*; 山崎 健一*; 内田 和也*; 松島 健*; 井上 智裕*; et al.
Tectonophysics, 915, p.230880_1 - 230880_11, 2025/11
We estimated the resistivity structure in the region of the 1997 Northwestern Kagoshima earthquake doublet (March 26, Mw6.1; May 13, Mw6.0) by dense magnetotelluric (MT) observations to investigate the influence of low-resistivity zones on the initiation and termination of large inland earthquake ruptures. In addition to MT data from previous studies, we acquired new broadband MT and telluric data at 42 sites. With a total of 86 sites in approximately 50 km 
50 km area, we constructed a high-resolution resistivity structure through 3-D inversion. A comparison between the slip distribution and resistivity structure showed that the rupture of the mainshock in March initiated at the western edge of a low-resistivity zone, propagated westward (with a maximum slip of 0.8 m), and was ultimately arrested by another low-resistivity zone. Based on geological structures and 
He/
He ratios in hot springs, we suggest that magma-derived fluids, located beneath a shallow, low-permeability smectite-rich zone and granodiorite, contributed to the initiation of the mainshock rupture. Another smectite-rich zone, along with underlying high pore pressure fluids, played a role in arresting the rupture. The rupture of the mainshock in May initiated at the edge of another low-resistivity zone and propagated in two directions (with a maximum slip of 0.4 m): one along the edge of the low-resistivity zone and the other toward the low-resistivity zone, before finally stopping within the low-resistivity zone. Our findings are consistent with those from the 2016 Kumamoto earthquake (M7.3), supporting the idea that resistivity structures can help assess the spatial potential of large inland earthquakes.
松谷 悠佑; 吉井 勇治*; 楠本 多聞*; 小川 達彦; 大西 世紀*; 平田 悠歩; 佐藤 達彦; 甲斐 健師
Physical Chemistry Chemical Physics, 27(14), p.6887 - 6898, 2025/04
被引用回数：2 パーセンタイル：83.51(Chemistry, Physical)水の放射線分解により生成されるラジカルは、DNA損傷誘発、染色体異常、発がんなど、放射線による生物影響の評価において重要な役割を果たす。粒子および重イオン輸送コードシステム(PHITS)では、あらゆる荷電粒子について水中の原子相互作用を推定できる飛跡構造解析モードと、ラジカルをシミュレート可能な電子線専用の化学コード(PHITS-Chem)が先行研究にて開発された。本研究では、あらゆるイオンビームに適用可能なPHITS-Chemコードを開発すると同時に、化学種間の反応をより効率的に検出する空間分割法や化学種の4次元可視化機能を整備した。更新されたPHITS-Chemコードは、文献にて報告される陽子線、
粒子線、炭素イオン線のG値と比較することにより検証され、PHITSオリジナル3次元描画ソフトPHIG-3Dによりラジカルの拡散動態を直感的に評価することに成功した。また、空間分割法の導入により、計算精度を維持しながら計算時間を大幅に短縮(約28倍高速化)することにも成功した。開発したPHITS-Chemコードは、粒子線治療においてラジカルにより誘発される生物効果の正確かつ直感的な理解に貢献することが期待される。
小川 達彦; 平田 悠歩; 松谷 悠佑; 甲斐 健師; 佐藤 達彦; 岩元 洋介; 橋本 慎太郎; 古田 琢哉; 安部 晋一郎; 松田 規宏; et al.
EPJ Nuclear Sciences & Technologies (Internet), 10, p.13_1 - 13_8, 2024/11
放射線挙動解析コードPHITSは、モンテカルロ法に基づいてほぼ全ての放射線の挙動を解析することができる放射線挙動解析計算コードである。その最新版であるPHITS version 3.34の、飛跡構造解析機能に焦点を置いて説明する。飛跡構造解析とは、荷電粒子が物質中を運動する挙動を計算する手法の一つで、個々の原子反応を識別することにより原子スケールでの追跡を可能にするものである。従来の飛跡構造解析モデルは生体を模擬する水だけにしか適用できず、遺伝子への放射線損傷を解析するツールとして使われてきた飛跡構造解析であるが、PHITSにおいてはPHITS-ETS、PHITS-ETS for Si、PHITS-KURBUC、ETSART、ITSARという飛跡構造解析モデルが補い合うことにより、生体の放射線影響だけでなく、半導体や材料物質など任意物質に対する適用が可能になっている。実際にこれらのモデルを使って、放射線によるDNA損傷予測、半導体のキャリア生成エネルギー計算、DPAの空間配置予測など、新しい解析研究も発表されており、飛跡構造解析を基礎とするボトムアップ型の放射線影響研究の推進に重要な役割を果たすことが期待できる。
菊地 晋; 近藤 俊樹; 土井 大輔; 清野 裕; 小川 賢吾*; 中川 武志*
Proceedings of 14th International Topical Meeting on Nuclear Reactor Thermal-Hydraulics, Operation, and Safety (NUTHOS-14) (Internet), 12 Pages, 2024/08
In this study, thermal behavior of NaI in different atmospheric conditions was investigated using thermal analysis technique to understand the formation behavior of gaseous iodine. Based on the results, it was revealed that the thermal behavior of NaI in inert atmosphere is phase change and approximately 10 mass% of I2 could be formed from NaI in 20 vol% oxygen atmosphere. The chemical equilibrium simulation was also performed to understand the mechanism of I2 formation. The estimated value of I2 formation in 20 vol% oxygen atmosphere was comparable to the result of thermal analysis.
平田 悠歩; 甲斐 健師; 小川 達彦; 松谷 悠佑; 佐藤 達彦
Nuclear Instruments and Methods in Physics Research B, 547, p.165183_1 - 165183_7, 2024/02
被引用回数：0 パーセンタイル：34.66(Instruments & Instrumentation)蛍光体の粒子線に対する発光効率は消光効果により低下することが知られている。蛍光体検出器を用いて正確な線量分布を得るためには、消光効果のメカニズムを理解することが不可欠である。本研究では、PHITSに実装された任意の物質に対する飛跡構造解析モードに基づいて蛍光体の発光強度を推定する新しいモデルを開発した。開発したモデルにより、BaFBr検出器の消光効果のシミュレーションが可能となり、その結果を対応する測定データと比較することにより検証した。このモデルは、様々なの蛍光体検出器の開発に貢献することが期待される。
佐藤 達彦; 岩元 洋介; 橋本 慎太郎; 小川 達彦; 古田 琢哉; 安部 晋一郎; 甲斐 健師; 松谷 悠佑; 松田 規宏; 平田 悠歩; et al.
Journal of Nuclear Science and Technology, 61(1), p.127 - 135, 2024/01
被引用回数：183 パーセンタイル：99.97(Nuclear Science & Technology)放射線挙動解析コードPHITSは、モンテカルロ法に基づいてほぼ全ての放射線の挙動を解析することができる。その最新版であるPHITS version 3.31を開発し公開した。最新版では、高エネルギー核データに対する親和性や飛跡構造解析アルゴリズムなどが改良されている。また、PHIG-3DやRT-PHITSなど、パッケージに組み込まれた外部ソフトウェアも充実している。本論文では、2017年にリリースされたPHITS3.02以降に導入された新しい機能について説明する。
平田 悠歩; 甲斐 健師; 小川 達彦; 松谷 悠佑*; 佐藤 達彦
Japanese Journal of Applied Physics, 62(10), p.106001_1 - 106001_6, 2023/10
被引用回数：5 パーセンタイル：44.36(Physics, Applied)半導体検出器の設計を最適化するには、半導体物質内において放射線がキャリア(励起電子)に変換されるまでの過程を理論的に解析する必要がある。本研究では、任意の半導体物質に対し、放射線により生じる二次電子の挙動を極低エネルギー(数eV)まで追跡し、励起電子が生成される過程を模擬できる機能(ETSART)を開発し、PHITSに実装した。具体的には、ETSARTを用いて計算した電子の飛程はICRU37で推奨されたデータ別の計算結果と一致することを確認した。さらに、半導体検出器の特性を表す重要な指標である、一つの励起電子の生成に必要な平均エネルギー(
値)について検討し、これまで値とバンドギャップエネルギーの関係は単純な直線モデルで考えられていたが、その関係は非線形関数であることを明らかにした。ETSARTは半導体検出器の最適化設計や応答解析に留まらず、新しい半導体物質の特性評価への応用も期待できる。
鶴見 美和*; 榎戸 輝掲*; 一方井 祐子*; Wu, T.*; Wang, D.*; 篠田 太郎*; 中澤 知洋*; 辻 直樹*; Diniz, G.*; 片岡 淳*; et al.
Geophysical Research Letters, 50(13), p.e2023GL103612_1 - e2023GL103612_9, 2023/07
被引用回数：9 パーセンタイル：75.81(Geosciences, Multidisciplinary)Gamma-ray glows are observational evidence of relativistic electron acceleration due to the electric field in thunderclouds. However, it is yet to be understood whether such relativistic electrons contribute to the initiation of lightning discharges. To tackle this question, we started the citizen science "Thundercloud Project," where we map radiation measurements of glows from winter thunderclouds along Japan's sea coast area. We developed and deployed 58 compact gamma-ray monitors at the end of 2021. On 30 December 2021, five monitors simultaneously detected a glow with its radiation distribution horizontally extending for 2 km. The glow terminated coinciding with a lightning flash at 04:08:34 JST, which was recorded by the two radio-band lightning mapping systems, FALMA and DALMA. The initial discharges during the preliminary breakdown started above the glow, that is, in vicinity of the electron acceleration site. This result provides one example of possible connections between electron acceleration and lightning initiation.
佐藤 達彦; 松谷 悠佑*; 小川 達彦; 甲斐 健師; 平田 悠歩; 津田 修一; Parisi, A.*
Physics in Medicine & Biology, 68(15), p.155005_1 - 155005_15, 2023/07
被引用回数：11 パーセンタイル：86.95(Engineering, Biomedical)PHITSには、マイクロドジメトリ機能と呼ばれるモンテカルロ法と解析関数を組み合わせて巨視的な空間内における微視的な領域の線量分布を計算する機能が備わっている。本論文では、そのマイクロドジメトリ機能を同じくPHITSに組み込まれた最新の飛跡構造解析モードを使って改良した結果について報告する。
小川 達彦; 平田 悠歩; 松谷 悠佑; 甲斐 健師
Isotope News, (784), p.13 - 16, 2022/12
入射荷電粒子が二次電子を生じる過程を原子サイズで明示的に計算する飛跡構造解析計算は、放射線生物影響,材料照射効果,放射線検出などの研究にとって重要な技術であり、近年主著者らの研究で新しい飛跡構造解析計算コードが開発された。従来の飛跡構造解析計算は標的物質の誘電関数を基に断面積を計算するため、誘電関数が良く測定されている水以外に、適用できるモデルは限られていた。本研究では誘電関数を使うことなく、二次電子エネルギー分布の系統式と阻止能を基に飛跡構造解析計算を行う手法により、誘電関数の測定値の有無にかかわらず、任意の物質で飛跡構造解析計算を実行することを可能とした。こうして開発したモデルで、陽子による水中の動径線量分布や二次電子生成量を計算したところ、従来のコードや実験値とよく一致した。このモデルは原子力機構の放射線輸送計算コードであるPHITS Ver3.25以降に実装され、任意物質に適用できる世界初の汎用飛跡構造解析コードとしてユーザーに提供されている。
平田 悠歩; 甲斐 健師; 小川 達彦; 松谷 悠佑; 佐藤 達彦
Japanese Journal of Applied Physics, 61(10), p.106004_1 - 106004_6, 2022/10
被引用回数：7 パーセンタイル：46.16(Physics, Applied)検出器や半導体メモリなどのSiデバイスにおいて、パルス波高欠損やソフトエラーなどの放射線影響が問題となっている。このような放射線影響のメカニズムを解明するためには、放射線による精密なエネルギー付与情報が必要である。そこで、Siにおける電子線のエネルギー付与をナノスケールで計算できる電子線飛跡構造解析機能を開発しPHITSに実装した。開発した機能の検証として電子の飛程や付与エネルギー分布を計算したところ、既報のモデルと一致することを確認した。また、一つのキャリア生成に必要なエネルギー(値)について、実験値を再現する二次電子生成のエネルギー閾値は2.75eVであることを見出すとともに、このエネルギー閾値は解析的に計算された結果および実験値と一致することがわかった。本研究で開発した電子線飛跡構造解析機能はSiデバイスに対する放射線影響の調査に応用することが期待される。
村瀬 清華*; 片岡 龍峰*; 西山 尚典*; 西村 耕司*; 橋本 大志*; 田中 良昌*; 門倉 昭*; 冨川 喜弘*; 堤 雅基*; 小川 泰信*; et al.
Journal of Space Weather and Space Climate (Internet), 12, p.18_1 - 18_16, 2022/06
被引用回数：3 パーセンタイル：31.12(Astronomy & Astrophysics)巨大な太陽フレアによってもたらされる太陽風により磁気圏内にエネルギーが溜まり、そのエネルギーが一気に解放されるサブストームが発生する。そのサブストームが発生する際、高エネルギー電子が大量に中間圏まで降り注ぐ事象(EEP)がしばしば観測されるが、その詳細な発生メカニズムは解明されていない。本研究では、あらせ衛星により観測された2つのEEPに対して、3次元グローバル電磁流体力学的(MHD)シミュレーションや放射線挙動解析コードPHITSを使った解析によりその発生メカニズムを検討した。その結果、カレントシート散乱とwave-particle散乱がEEPの初期及びサブストーム発生後に重要な役割を果たしていることが示唆された。
松谷 悠佑; 甲斐 健師; 佐藤 達彦; 小川 達彦; 平田 悠歩; 吉井 勇治*; Parisi, A.*; Liamsuwan, T.*
International Journal of Radiation Biology, 98(2), p.148 - 157, 2022/02
被引用回数：24 パーセンタイル：75.95(Biology)放射線による生物学的効果の研究を計算手法で進める場合、細胞内およびDNAスケールにおいて生体と等価物質である液体水中の各原子相互作用を明示的に考慮して、飛跡構造を解析することが重要となる。粒子・重イオン輸送計算コードPHITSは、独自の電子線飛跡構造解析モード(etsmode)ならびに、陽子・炭素イオンを模擬可能な世界的に有名なKURBUCアルゴリズムを使用することで、飛跡構造を詳細に計算できる。本研究では、電子線,陽子線,炭素線に関するPHITSの飛跡構造解析モードについて、物理的な特性(飛程・動径線量・微視的エネルギー付与)を評価し、文献の実験データや他のシミュレーションの文献値と一致することを検証した。また、電子線を対象に、早期の生物影響であるDNA一本鎖切断やDNA二本鎖切断、さらには複雑なDNA損傷であるクラスター損傷の発生数を推定し、電子線の入射エネルギー依存性を評価した。その結果、DNA損傷タイプはナノメートルスケールの電離励起の空間パターンに深く関与し、約500eVの電子線で複雑なDNA損傷の収量が大きくなることを明らかにした。PHITS飛跡構造解析モードの開発や検証並びに放射線生物研究へ応用した本成果は、放射線による生物学的効果の発生メカニズムの解明へ向けた研究を発展させるものである。
小川 達彦; 平田 悠歩; 松谷 悠佑; 甲斐 健師
Scientific Reports (Internet), 11(1), p.24401_1 - 24401_10, 2021/12
被引用回数：22 パーセンタイル：75.32(Multidisciplinary Sciences)入射荷電粒子が二次電子を生じる過程を明示的に計算する飛跡構造解析計算は、放射線生物影響,材料照射効果,放射線検出などの研究にとって重要な技術であり、これまで多くの研究に応用されてきた。しかし、従来の飛跡構造解析計算は標的物質の誘電関数を基に行われており、水以外の誘電関数はあまりよく知られていない。そこで本研究では誘電関数を使うことなく、二次電子エネルギー分布の系統式と阻止能を基に飛跡構造解析計算を行う手法を考案した。これにより、誘電関数の測定値の有無にかかわらず、任意の物質で飛跡構造解析計算を実行することを可能とした。まずこの手法で、標準的な検証データである陽子の水中の飛程や、エネルギー付与の動径方向分布を計算したところ、従来のコードや実験値とよく一致した。そこで、従来のコードが扱えなかった生体等価ガスを例にとり、同ガス中でのlineal energyを本研究で開発した手法により計算したところ、これも実験値の分布をよく再現し、水として近似した場合とは明確な差がみられた。このモデルは原子力機構の放射線輸送計算コードであるPHITS Ver3.25以降に実装され、任意物質に適用できる世界初の汎用飛跡構造解析コードとしてユーザーに提供される予定である。
松谷 悠佑; 甲斐 健師; 小川 達彦; 平田 悠歩; 佐藤 達彦
放射線化学(インターネット), (112), p.15 - 20, 2021/11
Particle and Heavy Ion Transport code System (PHITS)は、放射線挙動を模擬する汎用モンテカルロコードであり、原子力分野のみならず工学,医学,理学などの多様な分野で広く利用されている。PHITSは2010年に公開されて以降、機能拡張や利便性向上のために改良が進められてきた。今日までに、電子線,陽電子線,陽子線,炭素線の4種類の荷電粒子を対象として、液相水中における個々の原子との反応を模擬できる飛跡構造解析モードの開発を進めてきた。本モードの開発により、DNAスケールまで分解した微視的な線量付与の計算が可能となった。加えて、飛跡構造解析モードの高精度化へ向けて、任意物質中において多様な粒子タイプに適用可能な汎用的飛跡構造解析モードの開発も進められている。本稿で解説するPHITS飛跡構造解析モードに関するこれまでの開発経緯と将来展望により、PHITSコードの原子物理学,放射線化学,量子生命科学分野への応用がより一層期待される。
和田 有希*; 榎戸 輝揚*; 中村 佳敬*; 森本 健志*; 佐藤 光輝*; 牛尾 知雄*; 中澤 知洋*; 湯浅 孝行*; 米徳 大輔*; 澤野 達也*; et al.
Journal of Geophysical Research; Atmospheres, 125(4), p.e2019JD031730_1 - e2019JD031730_11, 2020/02
被引用回数：30 パーセンタイル：79.11(Meteorology & Atmospheric Sciences)During 2017-2018 winter operation of the Gamma-Ray Observation of Winter Thunderclouds experiment in Japan, two downward terrestrial gamma-ray flashes (TGFs) that triggered atmospheric photonuclear reactions were detected. They took place during winter thunderstorms on 5 December 2017 and 9 January 2018 at Kanazawa, Ishikawa Prefecture, Japan. Each event coincided with an intracloud/intercloud discharge, which had a negative-polarity peak current higher than 150 kA. Their radio waveforms in the low-frequency band are categorized as a distinct lightning type called energetic in-cloud pulse (EIP). Negative-polarity EIPs have been previously suggested to be highly associated with downward TGFs, and the present observations provide evidence of the correlation between them for the first time. Furthermore, both of the downward TGFs followed gamma-ray glows, minute-lasting high-energy emissions from thunderclouds. It is suggested that the negative EIPs took place with downward propagating negative leaders or upward positive ones developed in highly electrified regions responsible for the gamma-ray glows.
小川 達彦; 石川 法人; 甲斐 健師
Nuclear Instruments and Methods in Physics Research B, 461, p.272 - 275, 2019/12
被引用回数：7 パーセンタイル：50.12(Instruments & Instrumentation)物質中に局所的にエネルギーを付与する重イオン照射は、新しい材料改質方法の開発や放射線損傷原理の解明などを目的として広く用いられている。ここで、重イオンによって弾かれた電子(
線)と、それがさらに散乱する二次電子のエネルギー付与について、イオンの進行と垂直方向に対する分布は研究されてきた。一方、重イオンの進行方向に対するエネルギー付与分布は一定と仮定されてきたが、電子は重イオンの進行方向に押されるため、表面付近はエネルギー付与が少ないと推定される。そこで、本研究では、飛跡構造解析コードRITRACKSを用いて、重イオン照射を受けた物質内での二次電子のエネルギー付与の空間分布を詳細に求めた。その結果、水中では表面2nmまでの領域はエネルギー付与が少なく、それ以降は一定にエネルギーが付与されることを明らかにした。この結果を電子密度でスケーリングすることで水以外にも適用可能と期待され、核分裂生成物により損傷を受ける原子炉の燃料被覆管や、重イオンビーム照射で改質を受ける材料物質において、表面1-2nmでは放射線の影響が深部より少なくなり得ることを示唆している。
和田 有希*; 榎戸 輝揚*; 中村 佳敬*; 古田 禄大; 湯浅 孝行*; 中澤 知洋*; 森本 健志*; 佐藤 光輝*; 松元 崇弘*; 米徳 大輔*; et al.
Communications Physics (Internet), 2(1), p.67_1 - 67_9, 2019/06
被引用回数：62 パーセンタイル：92.79(Physics, Multidisciplinary)Two types of high-energy events have been detected from thunderstorms. One is "terrestrial gamma-ray flashes" (TGFs), sub-millisecond emissions coinciding with lightning discharges. The other is minute-lasting "gamma-ray glows". Although both phenomena are thought to originate from relativistic runaway electron avalanches in strong electric fields, the connection between them is not well understood. Here we report unequivocal simultaneous detection of a gamma-ray glow termination and a downward TGF, observed from the ground. During a winter thunderstorm in Japan on 9 January 2018, our detectors caught a gamma-ray glow, which moved for 
100 s with ambient wind, and then abruptly ceased with a lightning discharge. Simultaneously, the detectors observed photonuclear reactions triggered by a downward TGF, whose radio pulse was located within 
1 km from where the glow ceased. It is suggested that the highly-electrified region producing the glow was related to the initiation of the downward TGF.
和田 有希*; Bowers, G. S.*; 榎戸 輝揚*; 鴨川 仁*; 中村 佳敬*; 森本 健志*; Smith, D.*; 古田 禄大*; 中澤 知洋*; 湯浅 孝行*; et al.
Geophysical Research Letters, 45(11), p.5700 - 5707, 2018/06
被引用回数：32 パーセンタイル：77.95(Geosciences, Multidisciplinary)An on-ground observation program for high energy atmospheric phenomena in winter thunderstorms along Japan Sea has been performed via lightning measurements of 
-ray radiation, atmospheric electric field and low-frequency radio band. On February 11, 2017, the radiation detectors recorded -ray emission lasting for 75 sec. The -ray spectrum extended up to 20 MeV and was reproduced by a cutoff power-law model with a photon index of 1.36
, being consistent with a Bremsstrahlung radiation from a thundercloud (as known as a -ray glow). Then the -ray glow was abruptly terminated with a nearby lightning discharge. The low-frequency radio monitors, installed 50 km away from Noto School, recorded intra/inter-cloud discharges spreading over 60km area with a 300 ms duration. The timing of the -ray termination coincided with the moment when a sequence of intra/inter-cloud discharges passed 0.7 km horizontally away from the radiation monitors. The atmospheric electric-field measurement presented that negative charge was located in the cloud base and not neutralized by the lightning discharge. This indicates that the -ray source was located at an higher region than the cloud base.
