研究の初期の段階においては,多重ゼータ値の間に成立する種々の関係式と1変数多重対数関数のみたす接続関係をMellin変換と逆Mellin変換により対応付けすることを試みた.これにより多重ゼータ値の線型関係式のクラスとしては最も大きいと目される「多重ゼータ値に関する大野関係式」と「1変数多重対数関数の接続関係式(z→1-zとしたときの接続関係式)」の対応が明らかになった.このような具体的なレベルにおける1多重ゼータ値の関係式と多重対数関数の対称性との関係が明らかにされたことは大変意義深い,この結果は、論文「Relations for Multiple Zeta Values and Mellin Transforms of Multiple Polylogarithms, Pub1. RIMS, Kyoto Univ. 40(2004),537-564」として出版された.
この研究を契機として,1変数多重対数関数のみたす1変数KZ方程式の対称性(z=0,1,∞の特異点において正規化された基本解の間の変換理論=3次の対称群に同型)を形式的な水準(これは方程式の係数行列を単なる非可換変数と見倣して方程式を解析することを意味する)で記述することを試みた.各特異点近傍における正規化された基本解は,いわゆる,DrinfeldAssociatorにより接続され,基本解同士の接続関係からDrinfeldAsscoiatorのみたす双対関係式と6角形関係式が導かれる.この結果は「The Sum Formula of Multiple Zeta Values and Connection Problem of the Formal Knizhinik-Zamolodchikov Equation, Zeta Functions, Topology and Quantum Physics ed. by T. Aoki et al. Developments in Mathematics 14, Springer(2005)145-170」において公表された.
この論文の発表の後,研究の中心は多変数多重対数関数の接続関係を多変数KZ方程式の対称性(正規化された基本解同士の間の変換理論)から導くことを試みている.現在まで,2変数KZ方程式の対称性を一般的な水準で記述することに半ば成功している(基本解についての分解定理,解析性定理,また,Drinfeld Associatorが5角形関係式をみたすことなどを示すことができている.)一つの予想は,「基本解の分解定理から多重対数関数の調和積の関係式が導かれる」である.これらの結果は,2006年9月の日本数学会の秋季総合分科会(於大坂市立大学)における企画特別講演の予稿集で公表されている.

平成12年度から14年度において、申請者は『多重ゼータ値(MZV)』に関する研究を行なった。これは奥田順一(早稲田大学大学院博士課程在学)との共同研究であり、その結果は昨年度と今年度に発表した二編の論文、"New Approach to Ohno Relation for MultipleZeta Values"(by Jun-ichi Okuda and Kimio Ueno, arXiv : math.NT/0106148)、"Relations for Multiple Zeta Values and Mellin Transforms of Multiple Polylogarithms"(by Jun-ichi Okuda and Kimio Ueno, arXiv : math.NT/0301277)に纏めた。後者は前者の拡大版であり、これを論文雑誌に投稿した。これらの論文の続編は現在準備中である。研究の内容を一言で云えば、MZVのみたす関係式がメリン変換-逆メリン変換を通じて、多重高次対数関数(multiple polylogarithm、MPLと略す)のモノドロミー問題として解釈できることを示すことにある。発表した論文では、MZVの大野関係式とMPLのLanden接続公式の関係を書いたが、MZVとMPLの関係はそれにとどまるものではない。研究の経過については、これまで、日本数学会、京都大学数理解析研究所における短期共同研究会、あるいは、香川大学教育学部における講演会などにおいて発表してきたが、平成15年3月3日から6日まで近畿大学で開催された国除会議"Zeta Functions, Topology Quantum Physics"において共同研究者の奥田順一が理論の全容を初めて発表し(3月4日、奥田順一"Multiple Zeta Values and Mellin Transforms of Multiple Polylogarithms")、国内外の研究者の注目を集めた。研究分担者の米田元は、平成12年度から14年度において重力場方程式の解析に取り組んで論文を多数生産した。また、福島延久は可積分な格子模型の研究に取り組んだ。さらに、平成12年度と13年度、西澤道知は(平成13年度から東京大学に学振特別研究員として移籍)多重ガンマ関数とその楕円的類似の研究に新境地を切り拓いた

(1) 1997年-1998年、文部省科学研究費補助金(基盤C)「ゲージ項をもつDirac作用素の境界条件と不変量」により、ゲージ項をもったDirac作用素の指数公式を、グラスマン型の境界条件の下で調べた。とくに4次元半球面上における指数をAtiyah-Singerを用いずに計算した。
(2) さらに、物理でNinomiya-Tanの定理(Chiral anomalyの公式)と呼ばれる結果をゲージ付きDirac作用素の指数の不変性として精密化した。すなわち、ゲージ項をもつDirac作用素のS^4上の指数が、半球面上の幾何的Dirac作用素の(ゲージ項により定まる)グラスマン型の境界条件を与えた時の指数と等しいことを示した。これはゲージ項の効果が境界条件に吸収されることを言っている。この結果は早稲田大学理工学研究所Technical Reportに報告し、さらにRoskild大学のB,B,Booss教授を通してAMS Series"Contemporary Mathematics""Geometric Aspects of Partial Differetial Equationsに掲載される予定である。
(3) 以上の結果で、変数分離法によるS^3のスピノールのS^4へのゼロモード延長問題を、固有関数展開の形をさらに整理し、Bergmann核関数の類似により述べようと試みた。その結果コーシー核、それに続いて積分定理が得られた。そこからスピノール関数論と呼ぶべきものの構成を始め、特異点を持つゼロモードスピノールに対しローラン展開の類似を得た。C^2上の理論はほぼ完成している。これは次年度の萌芽研究に申請中である。

従来の数値相対論の方法であるADM形式と比較して,Ashtekar形式の特長の一つに変数の逆数を含まないという利点がある。この利点により縮退点通過の計算が可能になるのではないかと考えた。解析的な研究の結果,縮退点を直接通過しようとする方法では,変数またはその微分の一部が必ず発散することが分かった。そこで,変数を一時的に複素に拡張することで,縮退点を回避しながら計算し,また実領域に戻ってくる方法を考案した。これらについて縮退点通過の判定条件を設定し,実際に数値計算を試みた。その結果,縮退点を直接通過しようとする方法ではやはり計算は適切に機能しないことが分かった。そして,複素拡張による縮退点回避の方法では,その回避の仕方により適切に機能する場合としない場合が混在して存在することが確かめられた。この計算により,Ashtekar形式と複素拡張による縮退点回避のテクニックを組み合わせた方法により,縮退点通過が可能だと結論づけられる。解が離散的になるのは,本来片側にしか課さない境界条件を両側に課したことが原因であろう。どのような回避の仕方の時に計算が適切に機能するかは現在調査検討中である。また同様の試みをADM形式を使って行ってみたところ,解の分布に若干の違いが認められるが,Ashtekar形式による場合と概ね同様に通過が成功するものが発見された。これらの結果より,ここでの縮退点通過のポイントは複素多様体に拡張利用した縮退点回避だと結論づけられる。このような結果をまとめる論文を投稿準備中であるが,離散的に存在する解について,固有値問題特有の保存チャージを見つけることが出来ればさらに価値あるものになると思うので,研究中である

当研究の目的はEinstein方程式を数値的に解くときに適した離散化の方法を見つけることにある。数値的に解こうとするとき、Einstein方程式の大きな特徴は２つあり、１つは非線形性が強いこと、１つは1st classの束縛条件があることである。前者は数値的な安定性に関係し、後者は束縛条件の破れに関係する。特に後者について、離散化の良し悪しを評価する方法としてConstraint Accuracy Order(CAO)という方法を提案した。学会発表する予定であったが、学会が中止となってしまったので、2020年度以降に成果発表したい。このCAOを用いて、既存の様々な離散化手法に適用・評価し、数値実験でその実用性を確かめることに成功した。

Einstein方程式の数値シミュレーションの精度向上を目指して研究した．今年度は，主に次の2つの方法について調査研究を行った．(1) Einstein方程式を差分方程式にする際に適した差分化の方法を探す．(2) 発展方程式に拘束条件の破れを加えることでより安定的な数値計算を可能にする．(1)について成果は，Einstein方程式の２次摂動まで考慮した数値計算を行う際の，差分化の工夫を行い，精度向上することを確かめた．(2)については，拘束伝搬方程式の係数行列に非平坦な時空を代入してから，数値計算により固有値を求め，その値が数値計算の安定性と密接に関連していることを確かめた．

Einstein方程式の数値シミュレーションの精度向上を目指して研究した．今年度は，主に次の2つの方法について調査研究を行った．(1) Einstein方程式を差分方程式にする際に適した差分化の方法を探す．(2) 発展方程式に拘束条件の破れを加えることでより安定的な数値計算を可能にする．(1)について成果は，Einstein方程式の２次摂動まで考慮した数値計算を行う際の，差分化の工夫を行い，精度向上することを確かめた．(2)については，拘束伝搬方程式の係数行列に非平坦な時空を代入してから，数値計算により固有値を求め，その値が数値計算の安定性と密接に関連していることを確かめた．

Einstein方程式の数値シミュレーションを行う時の，適切な方程式および離散方程式の構築を目的として研究した．束縛条件の破れを抑えるために，時間発展方程式に束縛条件を付加する方法について，従来より汎用的な方法を提案した．これにより，背景時空が非平坦な場合や，近い背景時空が無い場合についても，束縛条件を付加したときの効果が事前に予測できるようになった．これらの成果を挙げて，発表を行った．また束縛条件を離散方程式においても保存する，いわゆる構造保存形にスキームについても研究し，Maxwell方程式を土台として構築し成果を挙げて発表した．

Einstein方程式の数値シミュレーションを行う時の，適切な方程式および離散方程式の構築を目的として研究した．束縛条件の破れを抑えるために，時間発展方程式に束縛条件を付加する方法について，従来より汎用的な方法を提案した．これにより，背景時空が非平坦な場合や，近い背景時空が無い場合についても，束縛条件を付加したときの効果が事前に予測できるようになった．これらの成果を挙げて，発表を行った．また束縛条件を離散方程式においても保存する，いわゆる構造保存形にスキームについても研究し，Maxwell方程式を土台として構築し成果を挙げて発表した．

Einstein方程式の数値シミュレーションを行う時の，適切な離散式を構築を目的として研究した．同じ束縛条件つき時間発展問題，という意味で同様の構造を持つMaxwell方程式について適切な離散式の構築に成功し，成果発表をした．さらにEinstein方程式について，特に時間発展のみに着目した限定的な意味での適切な離散式の構築に成功し，成果発表を行った．また，Einstein方程式の時間発展式に拘束条件を加える手法についても，従来の研究を推し進め，非平坦背景での拘束伝搬方程式の固有値を求め，数値安定性との一致についての結果を得て，成果発表を行った．

Einstein方程式の数値シミュレーションを行う時の，適切な離散式を構築を目的として研究した．同じ束縛条件つき時間発展問題，という意味で同様の構造を持つMaxwell方程式について適切な離散式の構築に成功し，成果発表をした．さらにEinstein方程式について，特に時間発展のみに着目した限定的な意味での適切な離散式の構築に成功し，成果発表を行った．また，Einstein方程式の時間発展式に拘束条件を加える手法についても，従来の研究を推し進め，非平坦背景での拘束伝搬方程式の固有値を求め，数値安定性との一致についての結果を得て，成果発表を行った．

一般相対性理論の基礎方程式であるアインシュタイン方程式は非線形偏微分方程式なので，特殊な対称性を課さない限り厳密解を得るのは難しい．特に時空のダイナミクスを追うのは，数値シミュレーションに頼るところが多い．この成果は重力波の観測成功などで大きな成果を挙げている．時空の場の方程式であるアインシュタイン方程式を時間と空間に分解したあと，離散化して数値スキームを作成する．そのときに用いられる方法は，他の方程式に特化して作成されたスキームを借用することが多い．そのため，構造保存型になっておらず，時間発展と共に拘束条件の破れが発生してしまう．そこでアインシュタイン方程式を数値シミュレーションをするための，構造保存型スキームの作成を目的として研究した．ハミルトニアンの正準形式を先に離散化し，それから時間発展方程式を求めていくという手法をとる離散変分法が知られている．この離散変分法を用いて，アインシュタイン方程式のスキーム作成に取り組んでいるが，まだ完全なものは得られていない．拘束条件つきの時間発展方程式の類題として，電磁気学の基礎方程式であるマクスウェルの方程式について，離散変分法を用いて構造保存型のスキームの作成に成功した．これは，拘束条件つきの時間発展方程式としては，初めて離散変分法を適用した例である．しかし，単に離散変分法を用いるだけでなく，様々な工夫も必要だということも分かった．これらの成果について，学会発表，論文投稿を行った．これをステップアップし，アインシュタイン方程式のスキーム作成にとりかかっている．だがマクスウェル方程式とアインシュタイン方程式では，ハミルトニアンの段階で，やや異なった構造を持っており，マクスウェル方程式で成功した離散変分法をそのまま適用するだけでは，部分的には成功するものの，完全な構造保存のスキームは出来ないことも分かった．現在その改良に取り組んでいるところである．

一般相対性理論の基礎方程式であるアインシュタイン方程式は非線形偏微分方程式なので，特殊な対称性を課さない限り厳密解を得るのは難しい．特に時空のダイナミクスを追うのは，数値シミュレーションに頼るところが多い．この成果は重力波の観測成功などで大きな成果を挙げている．時空の場の方程式であるアインシュタイン方程式を時間と空間に分解したあと，離散化して数値スキームを作成する．そのときに用いられる方法は，他の方程式に特化して作成されたスキームを借用することが多い．そのため，構造保存型になっておらず，時間発展と共に拘束条件の破れが発生してしまう．そこでアインシュタイン方程式を数値シミュレーションをするための，構造保存型スキームの作成を目的として研究した．ハミルトニアンの正準形式を先に離散化し，それから時間発展方程式を求めていくという手法をとる離散変分法が知られている．この離散変分法を用いて，アインシュタイン方程式のスキーム作成に取り組んでいるが，まだ完全なものは得られていない．拘束条件つきの時間発展方程式の類題として，電磁気学の基礎方程式であるマクスウェルの方程式について，離散変分法を用いて構造保存型のスキームの作成に成功した．これは，拘束条件つきの時間発展方程式としては，初めて離散変分法を適用した例である．しかし，単に離散変分法を用いるだけでなく，様々な工夫も必要だということも分かった．これらの成果について，学会発表，論文投稿を行った．これをステップアップし，アインシュタイン方程式のスキーム作成にとりかかっている．だがマクスウェル方程式とアインシュタイン方程式では，ハミルトニアンの段階で，やや異なった構造を持っており，マクスウェル方程式で成功した離散変分法をそのまま適用するだけでは，部分的には成功するものの，完全な構造保存のスキームは出来ないことも分かった．現在その改良に取り組んでいるところである．

アインシュタイン方程式は時空分解すれば，拘束条件つき時間発展方程式であり，数値シミュレーションする際は，スキームが必要となる．これを構造保存型にすることで数値シミュレーションの安定化を図るのが目的である．離散変分法を用いて，アインシュタイン方程式の構造保存型スキームを得ることに成功した．論文は投稿したが，実効性の検証が不十分だという理由で，差し戻されて改良，実証実験を実行中である．また，離散変分法を用いる際，２階微分の扱いに自由度があり，これがスキームの一意性を阻害していることが分かった．この点を考察，改良しつつ，論文発表，学会発表を計画している．

（Ｎ＋１）次元のＥｉｎｓｔｅｉｎ方程式をＮ次元超平面の時間発展として解析していくためには，数値解析による手法が不可欠である．数値解析をするためには，連続の方程式であるＥｉｎｓｔｅｉｎ方程式を差分方程式に直す必要がある．これには従来ではＣｒａｎｋ－Ｎｉｃｏｌｓｏｎ法やＲｕｎｇｅ－Ｋｕｔｔａ法がよく用いられている．しかしこれらは，他の方程式に対して安定的に数値解析可能なように開発された方法であり，Ｅｉｎｓｔｅｉｎ方程式のために開発されたものではない．ある近似の範囲でＥｉｎｓｔｅｉｎ方程式が有する性質を仮定して有効性が認められるのみである．具体的にはｃｏｎｓｔｒａｉｎｔが保存するという意味で構造保存の性質を有していない．Ｅｉｎｓｔｅｉｎ方程式の数値解析が必要とされる一般相対論的が支配的な状況では，強い重力場や長時間の時間発展などを扱う．その場合に上記の近似が無意味になり，数値解析をする時の最大の問題となるのはｃｏｎｓｔｒａｉｎｔの破れである．よって，構造保存の差分方程式を開発することは重要である．そこでＤｉｓｃｒｅｔｅ　ｖａｒｉａｔｉｏｎａｌ　Ｄｅｒｉｖａｔｉｖｅ　Ｍｅｔｈｏｄ（ＤＶＤＭ）という方法を使って，アインシュタイン方程式の差分化を行った．最初にＨａｍｉｌｔｏｎｉａｎ形式である必要があるので，Ｏｒｉｇｉｎａｌ　Ａｒｎｏｗｉｔｔ－Ｄｅｓｅｒ－Ｍｉｓｎｅｒ形式をベースにした．ＤＶＤＭを用いる例としては，２階微分を含む非線形の方程式では例がないので，それらの扱いを可能にする工夫を要した．その結果，構造保存の性質を持った差分化に成功した．数値解析によりその有効性を実験的に確かめた．また，同様なことをマクスウェルの方程式についても行い，従来の差分方程式と今回作成した差分方程式の違いを明瞭にした．これらの結果を応用数理学会（ＪＳＩＡＭ），日本数学会（ＪＭＳ），一般素体論と重力研究会（ＪＧＲＧ）において発表した．（発表は連携研究者の土屋拓也氏．）また結果をまとめたものを論文投稿中である．

Ｅｉｎｓｔｅｉｎ方程式の数値シミュレーションを行う際には，様々な工夫が必要である．現在の世界の数値相対論では，状況に応じて，様々な工夫を施しながら得たい精度まで達するように試行錯誤するという作戦が主なものである．しかし，なるべく汎用的に，長時間安定的かつ高精度な数値結果を得るためにどのような方法が可能か調査研究と検証を行った．その結果，他の方法に比べ，比較的汎用的に効果のあるｃｏｎｓｔｒａｉｎｔを使用した補正方法に注目した．従来の研究では，定数に固定されていた補正係数を，時間とともに変化するように制御するための，様々な作戦を検討し，実証を行った．ｃｏｎｓｔｒａｉｎｔの破れ量に対応して，柔軟に補正の大きさをコントロールする方法が数値的安定性に効果の高いことが分かった．具体的にはｃｏｎｓｔｒａｉｎｔの破れが小さいときは補正の大きさをゼロとし，破れが歩いていど大きくになるまでは，指数的に補正も大きくするが，ある程度以上に破れが大きいときは別の補正方法に切り替える作戦である．効果があることは大体わかったが，これを下支えする理論，あるいは系統的にこの作戦を行う処方箋などについては，いまだ未完成であり，これが完成した時点で学会発表，論文投稿したいと考えている．また，ｃｏｎｓｔｒａｉｎｔの破れ量に応じた制御ではなく，時間発展スキームとも関連した制御方法も検討して，研究実証を行っているところである．本来は双曲型微分方程式である時間発展方程式が，状況によってはそうでなくなり，放物型などに変化することがある．そうなると，Ｃｏｕｒａｎｔ条件などが変わってきてしまうので，数値計算の破綻につながるという可能性がある．それを避けるために，動的に補正係数を変化させたり，時間と空間の刻み幅を変化させたりする方法も検討した．まだ試行錯誤の段階で，発表できるような成果には至っていないが，検討的に研究調査すれば必ず成果が得られそうだという感触を得るまでにいたっている．

Ｅｉｎｓｔｅｉｎ方程式の数値シミュレーションにおける安定性について，全般的に考察，研究した．まず束縛条件の２乗和積分を定義し，それを利用して発展方程式への拘束条件の付加をする方法について，理論的かつ実証的に研究し，有効性が確認できたので，学会発表と論文発表を行った．また有限要素法によるシミュレーションへの応用も検討した．

Ｅｉｎｓｔｅｉｎ方程式の数値シミュレーションを安定させるための方法の1つとして，発展方程式への拘束条件の付加について研究した．その結果，つぎのような利点を持つ方法を開発した．従来の方法では事前評価時に背景を仮定する必要があったが，新しい方法では，背景に依らずに事前評価が可能になった．この方法をＥｉｎｓｔｅｉｎ方程式の数値シミュレーションへ適用する方法を提示し，その効果を実証した．

安定な数値シミュレーションを行うためのEinstein方程式の定式化について，束縛条件発展方程式の固有値解析が有効であることを理論，数値計算両面から確かめた．メキシコで行われた数値相対論の国際会議において成果を講演した．また，高次元数値相対論についての真貝寿昭との共著を執筆し掲載決定した．その概要は以下の通りである．従来の４次元時空モデルより高い次元の時空モデルは，様々な解釈を提供し，かつ異なるダイナミックな面を持っている．そのような最近の興味による高次元時空の数値相対論の研究のために，束縛条件の時間発展における振る舞いの次元の依存を調べる．N+1 Arnowitt-Deser-Misner進化方程式にはNに依存する物質項があります、しかし、束縛条件と束縛条件発展方程式は空間次元に不変の形式を持つ．これは我々が，4次元の場合のそうだったように，N+1次元数値相対論において，安定性，正確さについての問題があることを示唆するものである．しかし同時に，従来の束縛条件発展方程式の解析が，高次元においても有用であることも分かる．

　従来の数値相対論の方法であるＡＤＭ形式と比較して、Ａｓｈｔｅｋａｒ形式の特長の一つに変数の逆数を含まないという利点がある。この利点により縮退点通過の計算が可能になるのではないかと考えた。解析的な研究の結果、縮退点を直接通過しようとする方法では、変数またはその微分の一部が必ず発散することが分かった。そこで、変数を一時的に複素に拡張することで、縮退点を回避しながら計算し、また実領域に戻ってくる方法を考案した。これらについて縮退点通過の判定条件を設定し、実際に数値計算を試みた。その結果、縮退点を直接通過しようとする方法ではやはり計算は適切に機能しないことが分かった。そして、複素拡張による縮退点回避の方法では、その回避の仕方により適切に機能する場合としない場合が混在して存在することが、確かめられた。この計算により、Ａｓｈｔｅｋａｒ形式と複素拡張による縮退点回避のテクニックを組み合わせた方法により、縮退点通過が可能だと結論づけられる。解が離散的になるのは、本来片側にしか課さない境界条件を両側に課したことが原因であろう。どのような回避の仕方の時に計算が適切に機能するかは現在調査検討中である。また、同様の試みをＡＤＭ形式を使って行ってみたところ、解の分布に若干の違いが認められるが、Ａｓｈｔｅｋａｒ形式による場合と概ね同様に通過が成功するものが発見された。これらの結果より、ここでの縮退点通過のポイントは複素拡張による縮退点回避だと結論づけられる。このような結果をまとめる論文を投稿準備中であるが、離散的に存在する解について、固有値問題特有の保存チャージを見つけることが出来れば、さらに価値あるものになると思うので、研究中である。