• Search Research Projects
  • Search Researchers
  • How to Use
  1. Back to project page

2016 Fiscal Year Annual Research Report

2ループ輻射補正計算システムの構築

Research Project

Project/Area Number 15H03668
Research InstitutionHigh Energy Accelerator Research Organization

Principal Investigator

石川 正  大学共同利用機関法人高エネルギー加速器研究機構, 計算科学センター, 准教授 (90184481)

Co-Investigator(Kenkyū-buntansha) 加藤 潔  工学院大学, 基礎・教養教育部門(公私立大学の部局等), 教授 (50152707)
安井 良彰  東京経営短期大学, その他部局等, 准教授 (50389839)
台坂 博  一橋大学, 大学院商学研究科, 准教授 (80399295)
Project Period (FY) 2015-04-01 – 2019-03-31
Keywords素粒子論 / ファインマン・ループ積分 / 数理解析 / 多倍長計算 / 精密計算
Outline of Annual Research Achievements

素粒子物理学においてはヒッグス粒子の発見により、標準模型の精密な検証および標準模型を超える粒子を探索するという時代に突入した。将来の電子・陽電子衝突型のリニアコライダー実験により精密実験が切望される。標準模型における2ループの電弱相互作用まで含んだ理論計算が必要であり、標準模型の整合性について深く議論しなければならない。しかしながら、1ループの輻射補正に比べ、ダイアグラムの数が増えしかもファイマン・ループ積分も複雑になる。電弱相互作用の2ループ補正に関しては、近似式を用いる方法しかない。ミュオン異常磁気能率の2ループ補正電弱相互作用を例として、我々は理論的定式化のプログラム作成とファインマン・ループ積分を実行するための専用システムの研究開発を行った。専用システムについては計算能力増強を行った。
ミュオン異常磁気能率の2ループ・ファインマン・ダイアグラムは1780個あり、これらは14種類のトポロジーに分類される。トポロジー毎の次元正則化を数式処理により紫外部分と有限部分に次元正則化パラメータで展開する方法と、次元正則パラメータをそのまま変数として数値解析により数値結果を求める方法を開発して検証を行っている。繰り込みに関してはon-shell繰り込みを用いて対応するカウンター項の1ループ・ファインマン・ダイアグラムは70個あり、2ループ積分とカウンター項の紫外発散は数値的に相殺されることを確認し、梯子タイプのフェルミオン・ループがあるものについては数値結果を公表するに至った。
電弱相互作用計算では、ボゾンの質量/ミュオンの質量の二乗の桁、すなわち5,6桁ほどが異なるトポロジーで数値的な相殺が発生する。我々はこのため、多倍長積分計算可能な専用システムを構築すると共に用いている2重指数変数変換方式などの数値積分アルゴリズムについてもプログラム開発を行った。

Current Status of Research Progress
Current Status of Research Progress

2: Research has progressed on the whole more than it was originally planned.

Reason

ミュオン異常磁気能率の2ループ電弱相互作用の補正計算に関しては、2ループバーテックス型のトポロジー14種類に対して、次元正則化のパラメータでローラン展開して数値計算する方法と次元正則化のパラメータを変数とおいて数値解析して求める方法で検証を行う方法を確立した。バーテックス型梯子タイプの2ループにおいては、紫外発散部分及び有限部分について、第一の方法と比較して検証を行うことができた。
QEDダイアグラムのような赤外発散のあるダイアグラムでは、次元正則化のパラメータを負にすることにより数値計算で求めることができることを確認した。QEDの2ループ・ダイアグラムに関しては、1957年の計算は6桁の一致した。
三角形のフェルミオン・ループ・ダイアグラムでZ粒子の結合がある場合には、アノマリーが発生する。1992年の近似計算と比較することができて、第一世代、第二世代のフェルミオン・ループ・ダイアグラムのそれぞれの和はほぼ一致することがわかった。第三世代ではトップ・クォーク粒子のダイアグラムの導出不明な近似式とは一致していない。これらについては一部、高速で多倍長計算可能な専用システムにおいて検証も行った。これらについて国際会議で発表し論文を投稿した。
多倍長計算を高速で処理できる専用シウテムにはアクセラレータとして用いるFPGAの16ボードを追加して全体として64ボードとして運用を開始した。
数値積分においては二重指数変数変換を用いた台形公式で数値計算を行っており数値積分の誤差についても経験的評価式に沿っていることを確認した。次元正則化のパラメータを変数とおいて数値解析して求める方法について、研鑽を積むために、紫外発散を含むファインマン・ループ積分を行い、これらについて研究成果として論文を発表した。

Strategy for Future Research Activity

ファインマン・ループ積分は次元正則のパラメータでてローラン展開できることを用いている。数値を得るには先に展開した式を求めてから数値化するのが一般的であり、数式処理を用いて式を導いた。もう一つの方法として次元正則のパラメータを変数とおいて複数数値積分した後に得た後に数値解析する方法である。ミュオン異常磁気能率の2ループの電弱相互作用補正計算にこの2つの方法を適用して検証を行った。後者の数値的次元正則法については他のファインマン・ループ積分についてとりわけ紫外発散及び赤外発散について理論的考察を深めるため多倍長計算を高速に計算可能な開発した専用システムを用いて数値実験を行うことが必要であると考えている。
またミュオン異常磁気能率の2ループの電弱相互作用の補正完全計算については一部の数値結果を研究成果として公表したが完全な数値結果を求める予定である。すでに5つの非線形ゲージ項を導入しており、これらの5つをパラメータのゲージ不変性により検証を行うが、これらの計算には長いCPU時間が必要となる。ゲージ相殺が6,7桁に及ぶような多倍長計算が必要と考えられるダイアグラムについては多倍長計算専用システムを用いる。専用システム自体の高速化についてもハード・ソフト面で考察も引き続き行う。
平成29年度は2体から2体についての散乱過程における2ループ補正計算の自動化に向けての詳細設計を行う。ファインマン・ループ積分の積分空間内に分母が0となる超平面が一般的には存在する。ボックス型では6次元積分空間に5次元の超平面がある。3次元積分の場合までは数値的に求めることができるようになるが適応型アルゴリズムなどの研究開発が急がれる。
また国際リニアコライダーでは偏極の取扱うことが可能な電子・陽電子衝突実験であり、フェルミオンの偏極を取扱うことができるよう理論定式化の考察を行う。

  • Research Products

    (10 results)

All 2017 2016

All Journal Article (3 results) (of which Int'l Joint Research: 2 results,  Peer Reviewed: 2 results,  Acknowledgement Compliant: 2 results) Presentation (7 results) (of which Int'l Joint Research: 1 results)

  • [Journal Article] 多倍長精度積分計算を加速させる専用システム GRAPE9-MPX の開発とその応用2017

    • Author(s)
      台坂 博・中里直人・石川 正・湯浅富久子・似鳥啓吾
    • Journal Title

      電子情報通信学会 信学技報 IEICE Technical Report

      Volume: 116 Pages: 13-18

    • Acknowledgement Compliant
  • [Journal Article] High Performance and Increased Precision Techniques for Feynman Loop Integrals2016

    • Author(s)
      K Kato, E de Doncker, T Ishikawa, J Kapenga, O Olagbemi and F Yuasa
    • Journal Title

      Journal of Physics:Conference Series

      Volume: 762 Pages: 012070

    • DOI

      10.1088/1742-6596/762/1/012070

    • Peer Reviewed / Int'l Joint Research / Acknowledgement Compliant
  • [Journal Article] Adaptive Integration and Singular Boundary Transformations2016

    • Author(s)
      Elise de Doncker, Fukuko Yuasa, Tadashi Ishikawa, John Kapenga and Fola Olagbemi
    • Journal Title

      Procedia Computer Science

      Volume: 80 Pages: 1428-1438

    • DOI

      doi:10.1016/j.procs.2016.05.462

    • Peer Reviewed / Int'l Joint Research
  • [Presentation] 素粒子物理学における多倍長計算2017

    • Author(s)
      石川 正
    • Organizer
      第5回 大規模並列数値計算技術に関する研究集会
    • Place of Presentation
      計算科学研究機構(兵庫県・神戸市)
    • Year and Date
      2017-03-27 – 2017-03-27
  • [Presentation] ファインマン積分の数値計算法:特異性と変数変換の最適化2017

    • Author(s)
      台坂 博
    • Organizer
      日本物理学会 第72回年次大会
    • Place of Presentation
      大阪大学 豊中キャンパス(大阪府・豊中市)
    • Year and Date
      2017-03-17 – 2017-03-17
  • [Presentation] ファインマン積分の数値計算法:1-loop 4点テンソル積分の自動計算2017

    • Author(s)
      湯浅富久子
    • Organizer
      日本物理学会 第72回年次大会
    • Place of Presentation
      大阪大学 豊中キャンパス(大阪府・豊中市)
    • Year and Date
      2017-03-17 – 2017-03-17
  • [Presentation] Automatic calculation of two-loop ELWK corrections to Muon (g-2)2016

    • Author(s)
      N. Nakazawa
    • Organizer
      4th Computational Particle Physics Workshop
    • Place of Presentation
      IPC生産性国際交流センター(神奈川県・三浦郡葉山町)
    • Year and Date
      2016-10-10 – 2016-10-10
    • Int'l Joint Research
  • [Presentation] 素粒子反応の自動計算にむけて:多倍長精度計算システムによる2-loop 4点積分の計算2016

    • Author(s)
      台坂 博
    • Organizer
      日本物理学会 2016年秋季大会
    • Place of Presentation
      宮崎大学木花キャンパス(宮崎県・宮崎市)
    • Year and Date
      2016-09-22 – 2016-09-22
  • [Presentation] 素粒子反応の自動計算にむけて:1-loop 4点テンソル積分の完全数値的な取り扱い2016

    • Author(s)
      湯浅富久子
    • Organizer
      日本物理学会 2016年秋季大会
    • Place of Presentation
      宮崎大学木花キャンパス(宮崎県・宮崎市)
    • Year and Date
      2016-09-22 – 2016-09-22
  • [Presentation] 多倍長精度積分計算を加速させる専用システムの開発とその応用II2016

    • Author(s)
      台坂 博
    • Organizer
      日本応用数理学会 2016年度 年会
    • Place of Presentation
      北九州国際会議場(福岡県・北九州市)
    • Year and Date
      2016-09-12 – 2016-09-12

URL: 

Published: 2018-01-16  

Information User Guide FAQ News Terms of Use Attribution of KAKENHI

Powered by NII kakenhi