研究・技術カタログ
九州大学 マス・フォア・

インダストリ研究所
Institute of Mathematics for Industry
Kyushu University

目 次【数学テクノロジー先端研究部門】 2011/4/1 開設
離散最適化とその応用.....................................................................................................................................3
神山 直之 教授
学位：博士（工学）（京都大学） 専門分野：離散最適化、グラフ理論、計算量理論

Modeling, Optimization, and Control of Energy and Other Complex Systems......................4
Hoa Dinh NGUYEN 准教授

学位：PhD（Information Science and Technology）（The University of Tokyo）
専門分野：‌Modeling, optimization and control toward low-carbon and autonomous energy, transportation and otherinterconnected, complex systems. Particular focuses are on renewables and distributed energy resources, smart grid,
intelligent transportation, multi-agent systems, graph theory, artificial intelligence, and decentralized optimization.

【応用理論研究部門】 2011/4/1 開設
離散微分幾何・可積分系....................................................................................................................................5
梶原 健司 教授
学位：博士（工学）（東京大学） 専門分野：離散微分幾何、可積分系、
パンルヴェ系、離散・超離散系

流体力学、特に渦運動とその応用....................................................................................................................6
福本 康秀 教授
学位：理学博士（東京大学） 専門分野：流体力学、電磁流体力学、
トポロジカル流体力学

理論計算機科学とその応用..............................................................................................................................7
溝口 佳寛 教授
学位：博士（理学）（九州大学） 専門分野：計算機科学

生物に学ぶ適応ネットワーク理論.....................................................................................................................8
手老 篤史 准教授
学位：博士（理学）（北海道大学） 専門分野：数理モデリング、適応ネットワーク理論

最適化を基本とした研究・教育活動を目指して...............................................................................................9
脇 隼人 准教授
学位：博士（理学）（東京工業大学） 専門分野：最適化理論、連続最適化

計算理論と整数論の間..................................................................................................................................10
浦本 武雄 助教
学位：博士（理学）（京都大学） 専門分野：理論計算機科学

双曲型偏微分方程式の逆問題......................................................................................................................11
髙瀬 裕志 助教
学位：博士（数理科学）（東京大学） 専門分野：偏微分方程式、逆問題、幾何解析【基礎理論研究部門】 2011/4/1 開設
代数解析：数学の落ち合うところ..................................................................................................................12
落合 啓之 教授
学位：博士（数理科学）（東京大学） 専門分野：代数解析

トポロジーとその応用....................................................................................................................................13
佐伯 修 教授
学位：博士（理学）（東京大学） 専門分野：位相幾何学、トポロジー、特異点論、微分位相幾何学、DNA結び目

確率論とその応用..........................................................................................................................................14
白井 朋之 教授
学位：博士（数理科学）（東京大学） 専門分野：確率論

数論的不変式論とそれにまつわる幾何........................................................................................................15
石塚 裕大 助教
学位：博士（理学）（京都大学） 専門分野：数論、数論的不変式論、
ディオファンタス幾何学

数理科学的観点からの血糖値管理...............................................................................................................16
小谷 久寿 助教
学位：博士（数理学）（九州大学） 専門分野：位相幾何学、整数論、数論的位相幾何学

4次元図形の設計図：視えない次元を視る....................................................................................................17
濵田 法行 助教
学位：博士（数理学）（九州大学） 専門分野：低次元トポロジー、曲面の写像類群、
シンプレクティックトポロジー【数理計算インテリジェント社会実装推進部門】 2019/4/1 開設
トポロジーとその諸科学への応用.................................................................................................................18
鍛冶 静雄 教授
学位：博士（理学）（京都大学） 専門分野：位相幾何学、
リー群、応用数学

1

Institute of Mathematics for Industry
Kyushu University

九州大学マス・フォア・インダストリ研究所

物質構造解析の数理......................................................................................................................................19
富安（大石） 亮子 教授
学位：博士（数理科学）（東京大学） 専門分野：応用代数・数論、数理結晶学、アルゴリズム

グラフ解析と最適化問題の高速計算及び実社会への応用.........................................................................20
藤澤 克樹 教授
学位：博士（理学）（東京工業大学） 専門分野：最適化問題、グラフ解析、高性能計算

圧縮性回転粘性流体の数学解析..................................................................................................................21
藤井 幹大 助教
学位：博士（数理学）（九州大学） 専門分野：流体力学に現れる偏微分方程式の数学解析【産業数理統計研究部門】 2022/4/1 開設
L1正則化に基づくスパース多変量解析.......................................................................................................22
廣瀨 慧 教授
学位：博士（機能数理学）（九州大学） 専門分野：スパース推定、L1正則化、多変量解析

Finite mixture models for statistical inference..............................................................................23
Hien Duy NGUYEN 教授

学位：PhD（University of Queensland, Australia）
専門分野：‌Mathematical Statistics, Statistical Computing, Statistical Learning, Bayesian Statistics, Signal
Processing, Stochastic Programming, Optimization Theory

位相的データ解析と層の理論.......................................................................................................................24
池 祐一 准教授
学位：博士（数理科学）（東京大学） 専門分野：位相的データ解析、超局所層理論

ベイズ推定による計測と数理モデリングの橋渡し.......................................................................................25
德田 悟 准教授
学位：博士（科学）（東京大学） 専門分野：ベイズ推定、モデリング、情報統計力学

統計的ダイバージェンスに基づく頑健なモデル選択規準...........................................................................26
倉田 澄人 助教
学位：博士（理学）（大阪大学） 専門分野：数理統計学、
モデル選択、
ロバストネス

小区分別統計的推測法の理論発展及びその活用.......................................................................................27
廣瀬 雅代 助教
学位：博士（工学）（大阪大学） 専門分野：統計科学、小地域推定、混合効果モデル【先進暗号数理デザイン室】 2015/4/1 開設
数学と暗号の交差点......................................................................................................................................28
縫田 光司 教授
学位：博士（数理科学）（東京大学） 専門分野：暗号数理、秘密計算、組合せ論的群論

耐量子暗号の安全性解析..............................................................................................................................29
池松 泰彦 助教
学位：博士（数理学）（九州大学） 専門分野：耐量子暗号、多変数多項式暗号【リエゾン戦略部門】 2023/4/1 開設「変な振る舞い」＝「素直な振る舞い」...........................................................................................................30
松江 要 教授
学位：博士（理学）（京都大学） 専門分野：力学系、数値解析、精度保証付き数値計算、特異摂動論、燃焼など

ソーシャル数理と動的最適化........................................................................................................................31
吉良 知文 准教授
学位：博士（機能数理学）（九州大学） ソーシャル数理、動的最適化、確率最適化、
マルコフ決定過程

数値シミュレーションによる現象の理解.......................................................................................................32
田上 大助 准教授
学位：博士（数理学）（九州大学） 専門分野：数値解析、計算力学【オーストラリア分室】 2015/3/1 開設
Modeling of Solid-to-Solid Phase-Transformations in Shape-Memory Alloys Homogenization andGamma-Convergence Problems for Nematic Elastomers........................................................................33
Pierluigi CESANA 准教授
学位：PhD（Applied Mathematics）（SISSA International School for Advanced Studies, Italy）
専門分野：Partial Differential Equations, Variational Problems

Algebraic specification...........................................................................................................................34
Daniel GAINA 准教授
学位：PhD（Japan Advanced Institute of Science and Technology）
専門分野：Universal Logic, Formal Methods, Category Theory2◯ 数学テクノロジー先端研究部門

離散最適化とその応用

神山直之

学位：博士（工学）（京都大学）
専門分野：離散最適化、グラフ理論、計算量理論
私の専門分野は離散最適化の理論的研究です。加えて、位相的にどのように接続しているかのみに注目したものです

離散最適化と関係の深いグラフ理論や計算量理論の研究（図２参照）。直観的には、道路等のネットワークを数学的にも行っています。これらの分野は数学と計算機科学の落ち表わしたものだと思っていただければよいでしょう。そして、合うところに位置し、互いに深く絡み合っています。これらのグラフ理論とは様々なグラフの持つ特性を一般的に明らか分野が含む問題は非常に多岐に渡るのですが、私は離散的にする学問です。例えば、どのようなグラフならば高い頑健な凸性である劣モジュラ性や、双対性を基本とする多面体

性を持っているかを考えたりします。このグラフ理論におい的手法などを通じて、統一的な理解を深めることを目標としては、私は古典的な問題である有向木の詰込問題におけるています。以下では、私が興味を持っている問題と併せて、最大最小定理や（図２参照）、
グラフ上の最適化問題である

これらの分野の簡単な説明をさせていただきます。

辺支配集合問題などの研究を行っています。

まず離散最適化ですが、そもそも最適化とは幾つかの解の候補から、目的関数を最大化もしくは最小化する解を求める問題です。離散最適化の研究においては、その中でも解の候補が離散的（つまりばらばら）な性質を持つ問題を扱います。これらの問題に対して、解の集合がどのような性質を満たすならば、全ての解の候補を調べずに効率的に最適解を見つけることができるかを明らかにし、実際にその解を求める方法（アルゴリズムと呼ばれます）を構築することが

図２：（左）有向グラフ。
（右）有向木。赤い点を根とし木のように広がる部分グ
ラフ。

大きな目標となります。この分野に対して、私は二つの集合
間の良い割り当てを求める安定マッチング問題や（図１参照）、
ものの流れを数理モデル化したネットワークフロー、そ

最後に計算量理論です。上記の二つの分野がどちらかとして離散最適化における抽象的な枠組みである劣モジュラ

いえば「できる」
ことを研究する学問だったのに対し、
この分関数やマトロイドに関係する最適化問題などの研究を行っ野は「できない」
ことを研究する分野です。具体的には、計算ています。

機における「計算」を数学的に定義しその能力の限界を研究するものです。この分野においては、
クレイ数学研究所の発表したミレニアム懸賞問題の一つである「Ｐ対ＮＰ問題」

1
2

2

2

が象徴的な問題だといえるでしょう。一見、最適化と計算量

1

1

く近年計算量理論の研究においても最適化の手法が使用

1
2

1
2

1

1

2

2

理論の分野は相反するものに見えますが、実はそうではな
されるようになってきています。私はこの計算量理論の研究

1

において、離散最適化の分野で用いられている手法、例えば
多面体的アプローチ等をどのように活用することができる

2

かを明らかにすることに興味があります。
また、これらの分野の研究は、現実社会の問題を数理モデ
ル化し解決するオペレーションズ・リサーチと呼ばれる分野と非常に相性が非常に良く、私自身も交通や都市計画、そし

図１：（左）安定マッチング問題の例。左と右の点間の良い割り当てを求める。
付記されている数字は相手に対する選好順序。
（右）割り当ての例。

てソーシャルネットワーク等の現実的な問題から生じる問題に対して得られた理論的成果を応用することに興味があり続いては、グラフ理論です。グラフとは点とそれらをつなます。

ぐ辺で構成されるもので、幾何学的な位置関係ではなく、

3

Institute of Mathematics for Industry
Kyushu University

九州大学マス・フォア・インダストリ研究所

Modeling, Optimization, and Control of Energy and Other Complex Systems
数値シミュレーションによる現象の理解

Hoa Dinh NGUYEN
学位：PhD (Information
田上
大助 Science and Technology)(The University of Tokyo)

専門分野：Modeling, optimization and control toward low-carbon and autonomous energy, transportation and other
学位：博士（数理学）（九州大学）
interconnected,
complex systems.
Particular focuses are on renewables and distributed energy resources, smart grid,
専門分野：数値解析、
計算力学
intelligent transportation,
multi-agent
systems, graph theory, artificial intelligence, and decentralized optimization.
た。
またこの成果を発展させ、ガラス溶融炉の最適設計を行う際のderivedmodels for optimizing andcontrolling such complex systems,

My research theme is Modeling, Optimization, and Control of Energy
and Other Complex Systems. An illustration of such complex systems is
depicted in Fig. 1.
The ultimate goals of this research are: (i) low-carbon, decentralized,
resilient, autonomous, and comfortable energy systems; and (ii)harmonization of interconnected and complex systems. To achieve that,the following studies will be conducted. First, mathematical models will
be built for representing the multi-scale, complicated dynamics andbehaviors in energy systems, as well as other systems in engineering andnature. Then mathematical frameworks will be developed based on thegiven specific system objectives.
指標の1つである炉全体の熱収支計算に対して整合流束法に基づ
A variety of branches in 熱収支計算に対する最適な誤差評価を導mathematics will be employed in this research,
く離散化手法を提案し、for example
dynamical systems
and their stability theories, linear andいた。
我々の誤差評価は、境界積分で定義される物理量を直接計算
nonlinear optimization, graph theory, control theory, multi-agent system,
するよりも整合流束法を用いて領域積分に置き換えて計算した方and artificial intelligence.
Particular applications include smart grid, smart cities, renewable and
ように数理的正当化した離散化手法に基づく数値シミュレーション
distributed energy systems, network systems, wireless power transfer,を用いて、製品品質の維持とエネルギー消費の抑制を両立したガ
intelligent and decarbonized transportation systems, bio systems for
ラス溶融炉の最適設計について検討を試みている。
clean energy production and carbon capture.
が高精度である、という事実の数理的正当化に対応している。このその他に我々が取り組んでいる研究の一つとして、変圧器内部の渦電流に代表される磁場の数値シミュレーションがある(図3参照)。
図1.数値シミュレーションの概念図。斜体文字は非圧縮粘性流れの場合の対応する例。

複雑な形状を持った領域や変化の著しい現象の数値シミュレーショ
ンを精度良く行うためには、自由度数が107(あるいはそれ以上)の水の流れや熱の伝達など、自然界や産業界で見られる様々な現近似方程式を効率良く扱うことがしばしば要求される。我々は静磁象を、計算機を用いた数値シミュレーションによって理解することに場問題に対して、既知量である電流密度の補正を考慮した混合型

興味を持っている。数値シミュレーションでは、まず現象を物理法則定式化による数理モデルを導入した。この数理モデルに対して反復に基づいて微分方程式で記述する 数理モデル化 を行い、次に微型領域分割法に基づく近似方程式を提案し、さらに混合型定式化の分方程式を計算機で扱うことのできる近似方程式に置き換える 離際に現れるLagrange乗数の性質を利用した簡略化を行った。提案

散化 を行い、最後に近似方程式の解法を計算機に実装し目的とすした手法は、必要となる行列ベクトル積演算を部分領域ごとに独立る現象を再現する 数値計算 を行う(図1参照)。我々は、提案するな静磁場問題を解くことに帰着できるため、計算効率向上のために離散化 手法が本来の現象に対する 精度の良い 離散化になって用いられる並列計算に適しているという利点がある。以上のことかいるか、実用に耐えうる 効率の良い 離散化になっているか、を数理ら、提案した手法に現れるある反復計算の収束特性が改善され、計的に検討する研究に取り組んでいる。また様々な問題に対して、提算効率が向上した。これにより、従来は解くことのできなかった大自案する 離散化 手法に基づいた 数値実験 を実際に行い、自然現象

由度モデルの解析に成功した。このように効率化した離散化手法にの理解や工業製品の設計へ応用することにも取り組んでいる。

基づく数値シミュレーションを用いて、磁場の振る舞いをより精度

我々が取り組んでいる研究の一つに、ガラス溶融炉内部における良く把握することを通して、変圧器の最適設計についての検討を試

ガラス原料の流動や熱伝達の数値シミュレーションがある(図2参みている。

照)。ガラス溶融炉内部にある高温(約1,500℃)のガラス原料は非

さらに我々は、粘弾性流れ問題、移動境界流れ問題、光波の干渉・

圧縮粘性流体と見なせることから、溶融炉内部における現象は温度散乱問題などの数値シミュレーションにも関心を広げて研究を続け

依存係数を持つ熱対流方程式によって数理モデル化できる。我々ている。これにより、数理的正当化がなされた数値シミュレーションは時間方向に後退Euler法を、空間方向に混合型適合有限要素法

がより多くの現象に対して可能になると期待できる。

Figure 1. Illustration for power systems decision and control over different time scales.

を用いて得られる近似方程式の解に対する最適な誤差評価を導い

図2.あるガラス溶融炉内部における熱対流問題の数値シミュレーションによって得られた、ある時刻におけるガラス原料の温度分布。上図は炉に電極を挿入した場合、下図は炉に電極を挿入しない場合である。電極を挿入し電流を流すことで発生するJoule熱が対流の駆動力となるため、温度分布に違いが生じていることが分かる。

図3.ある変圧器内部における時間調和な渦電流問題の数値シミュレーショ
ンによって得られた、渦電流の流線。2つの図は、変圧器に対して互いに異なる側にある視点から鳥瞰した場合である。変圧器の壁面内部に設置された磁気シールド(水色や黄緑色の流線で囲まれた青い長方形部分)を囲むように渦電流が生じていることが分かる。

4

◯ 応用理論研究部門

離散微分幾何・可積分系

梶原健司

学位：博士（工学）（東京大学）
専門分野：離散微分幾何、可積分系、パンルヴェ系、離散・超離散系
ソリトンと呼ばれる、粒子性をあわせ持つ非線形波動の研究に端で表示される「楕円超幾何函数」が現れます。この函数は超幾何型を発する「可積分系」の研究を軸に研究活動を行っています。数理

特殊函数の頂上に位置するものと考えられています。離散微分幾

的には、ソリトンを記述する基礎方程式（ソリトン方程式）は非線形

何における応用の一つとして、ある種の複素正則函数の離散化をパ偏微分方程式であるにもかかわらず厳密に解けるという、奇跡的なンルヴェ・離散パンルヴェ方程式が記述します。図２は冪函数Z１/２の性質を持っています。その奇跡の背後には「無限自由度の対称性を離散化で、グリッドはサークルパターンで特徴付けられます。

もつ無限次元の空間」の数理があり、その数理を共有する函数方程式の族を「可積分系」と呼びます。背後の数理を深く理解することによって、可積分系は多くの分野への応用が可能です。以下、私が関わっている３つの例を挙げます。
1. 離散化と超離散化：ソリトン方程式の可積分性を保存したまま
独立変数を離散化して差分方程式を構成したり、従属変数まで離散
図３：‌竜巻に見られる渦糸（左）、渦糸のループと渦輪の離散モデルによる計算結果（右）

化してセルオートマトンを構成する方法（「超離散化」）が構築されています。図１上は典型的なソリトンの相互作用を表し、左からやってきた大きく速いソリトンが小さく遅いソリトンを追い抜いていま

3. 離散可積分微分幾何と幾何学的形状生成：空間中の曲線・曲面

す。一方、図１下はソリトンを記述するセルオートマトンです。
１次元やその変形を記述する基礎方程式としてさまざまな可積分系が現の箱の列と玉があるとし、各時刻で左から玉を右の最も近い空箱にれることが知られています。さらに最近、離散可積分系と整合する

移し、全ての玉を移動したら時刻を１進めます。この単純なモデル

曲線・曲面論の構築が進み、私は幾何学的形状生成、特に何らかのがソリトンを記述し、超離散化の手続きで偏微分方程式と直接の対意味で「よい」曲線・曲面の生成の研究を工業意匠設計や建築の専応がつきます。離散化・超離散化の理論は数値解析や交通流など広門家と協力して進めています。図３は天然の渦糸である竜巻と、そ範な分野に応用され、大きな成功を収めており、私が展開する他分の離散モデルの計算結果。背後の構造を保存しているので、コスト

野、特に幾何に関係する分野との連携活動のバックボーンを成しての低い計算で高品質な結果が得られます。図５（左）は工業意匠設います（離散微分幾何）。最近では超離散化の適用範囲が可積分で計において、美的要素をもつ形状要素として日本で提唱された「対はない反応拡散系にも拡大され、性質のよい反応拡散セルオートマ

数型美的曲線」と呼ばれる平面曲線の族で、日本刀や蝶の羽など、トンモデルが構成できるようになっています。さらに、超離散系の背私たちが「美しい」と思う形状から抽出されました（図４）。最近、私後の構造が最近急速に発達している「トロピカル幾何学」でよく記たちは可積分幾何の立場から全く新しい理論的枠組みを提唱し、そ述されることが判明し、純粋数学の発達も促しています。

れを用いて構造を保存した、高速生成可能で高品質な離散化や（図６）、空間曲線への一般化を構築しました（図５（右））。この理論を曲面に拡張し、得られた美的要素をもつ曲線・曲面の族を幾何学的形
状要素として実装し、意匠設計での標準化を目指して研究を進めます。これらの形状は「美しさ」という人間の感覚の要素を取り込んだ数理モデルとしても重要な例です：よい方程式は、よい形状を生成する。

図１：‌ソリトンとセルオートマト
ンモデル

図２：離散冪函数とサークルパターン

2. 離散パンルヴェ方程式と楕円曲線：離散パンルヴェ方程式と呼ばれる可積分な２階常差分方程式の族は、数理物理学や特殊函数論で重要な役割を果たし、背後に極めて豊富な構造があります。そ図４：‌日本刀や蝶に見られる「美しい曲線」 図５：対数型美的曲線とその空間曲線版 図６：対数型美的曲線とその離散化の解はよく知られたベッセル函数や超幾何函数などの特殊函数の一般化と見なすことができます。これらの頂上にE型の対称性を（１）８持つ「楕円パンルヴェ方程式」があり、特殊解として楕円テータ函数

5Institute of Mathematics for Industry
Kyushu University

九州大学マス・フォア・インダストリ研究所

流体力学、特に渦運動とその応用

福本

康秀

学位：理学博士（東京大学）
専門分野：流体力学、電磁流体力学、
トポロジカル流体力学
流体運動を象るのは渦と波である。飛行機の背後に長くのびる

平均流の３次元相互作用を計算する新しい数学的枠組みの構築、２本の筋状の雲、浴槽の排水口の上に立つ気泡の柱などのパター

その回転・成層効果の取り込み、電磁流体力学への拡張を進めていンは渦構造の典型例である。渦は無限自由度力学系で、大小様々なる。最近では、燃焼火炎面のダイナミックスなど圧縮性流体におけスケールのモードが非線形的に相互作用しながら自律的に特定のる渦運動に取り組んでいる。2013年３月には、福岡で、議長として配位に発達、あるいはカオス・乱流状態に至る。この過程で、渦は、運渦運動に関するIUTAMシンポジウムを開催した。2015年１月よ動量や物質を「輸送する」、「混ぜる」
といった機能を獲得する。空気り、国際学術誌Fluid Dynamics Research誌（IOPP）のEditor-砲が特定の空気塊を遠くまで搬送することはお馴染であるが、心臓

in-Chiefを務めている。

ポンプによる血流の駆動をはじめ生命体は様々な形で渦輪を利用

文部科学省グローバルCOEプログラム「マス・フォア・インダス

している。さらに、翼端渦の形成と安定性が、航空機の安全航行、風トリ教育研究拠点」（ 九州大学数理学府、拠点リーダー若山正人、

力発電の効率向上から飛翔ロボットの設計にまでかかわることに代

2008～2012年度）拠点サブリーダーを務めた。2008年度か表されるように、渦の非線形ダイナミックスはエネルギーや地球環ら始まったマス・フォア・インダストリ・フォーラム（FMfI）
や2010年

境問題から先端テクノロジー分野にわたる今日的課題解決の鍵を度開始のスタディ・グループ（SGW：企業の未解決問題に数学研究握る。

者・大学院生が短期間集中的に取り組むが合宿）の運営に参画し、マ
渦運動の数理解析が主た

ス・フォア・インダストリ研究所（IMI、2011年度～）の立ち上げ、IMIるテーマである。特に、３次元の共同利用・共同研究拠点認定（2013年度～）に関わった。また、

渦運動の理論において世界

IMIが幹事拠点として受託した文部科学省委託事業「数学アドバンに先行する結果をもつ。平成

ストイノベーションプラットフォーム」（AIMaP、2017～2021年

６年には、渦糸の３次元運動度）の運営に携わり、全国12の有力数学・数理科学研究機関の協力に関する論文に対して、若手を得て、数学と諸科学分野・産業技術との連携を推進する様々な取を対 象とする日本 流 体 力 学

り組みを進めている。目下、社会科学や人文科学への産業数学フロ

会竜門賞の第１回受賞者に選

ンティアの拡大に努めている。

ばれた。最近、
レイノルズ数の博士後期課程を含む大学院生の指導に力を注ぎ、外国人留学生

大・小両極限で、実験とよく合の受け入れも積極的に行っている。
５名以上の指導学生を国内企業

う「渦輪」の進行速度の公式に長期研究インターンシップ（３カ月以上）に派遣した。研究者レベを導出することに成功した。まルでの国際交流も盛んに行っている。日本学術振興会特定国派遣

た、「 曲率不安定性」という渦研究者（長期）に採用され、1996年に10か月間、ケンブリッジ大学

輪の新しい不安定機構を発見を訪問し、Moffatt教授と渦輪の運動に関する共同研究を行った。

した。以降、無限自由度ハミルトン力学系の視点から、渦の点・連続

帰国後、有力外国人研究者の訪問が絶えることなく、最前線の情報

スペクトル、渦の非線形運動を計算する新しいラグランジュ的アプを直接交換している。2001年以降、日本学術振興会招へい研究

ローチの開発を進めている。

者（短期）として外国人を10名以上受け入れた。2016年には、代流体運動の偏微分方程式による解析を創始したのは18世紀の表として、ニュージーランドの産業数学グループとSGWの相互乗

オイラーであるが、渦運動の概念はヘルムホルツの論文（1858）

り入れを行い、
また、FMfIのアジア太平洋産業数学コンソーシアム

まで１世紀待つ。ヘルムホルツは、粘性がないとき、「渦線が流体に（APCMfI）加盟機関による海外開催（2016年～）を進めている。

凍結して運ばれる」
ことを示した。このことはとりもなおさず渦線のこうした活動の積み重ねの結果、創始者Moffatt教授や第一人者

絡み目型・結び目型が時間的に変化しないことを意味する。ヘルム

Ricca教授（Milan 大）を中心とする「トポロジカル流体力学」、そホルツの法則のトポロジー的意味を汲み出しその応用をはかるとして、
ラ・
トローブ大学などAPCMfIや欧米の指導的研究者と「産業いうのが、Arnold（‘66）に始まる20世紀後半の流体力学の１つ数学」の世界的ネットワークを築き上げるに至った。このネットワーの流れである。
しかし、現状は、２次元流に留まる。流体粒子の変位

クを若手研究者の交流や育成に役立てていきたい。

を基本変数とするラグランジュ的記述によってはじめて、トポロジカ
ル不変量を厳密に保ちながら渦運動を扱うことができ、分子運動から、固体、流体、弾性体、プラズマの運動までを貫く共通の土壌を与えてくれる。ラグランジュ的記述のもつ高い拡張性を活かして、波と6◯ 応用理論研究部門

理論計算機科学とその応用

溝口佳寛

学位：博士（理学）（九州大学）
専門分野：計算機科学
グラフ構造は、複雑な状況を直観的に説明するための基本的な道具として計算機が生まれる前からその構造についての研究や応用が行なわれていました。計算機誕生後は、計算機の内部表現(データ構造)としてプログラム中でグラフが利用されるばかりか、
グ
ラフそのものを変換して計算を行なう新しい計算システムなども考案され、従来の静的なグラフ構造の研究(グラフ理論)だけではな
く、グラフを操作することに対する性質の変化を究明する研究(グラフ変換理論)が特に重要になっています。また、グラフ変換操作を既存の電子計算機上で実現することを前提にするのではなく、ある種の基本的グラフ変換操作がハードウェア的に実現された計算機を想定してのアルゴリズム構築の視点も必要です。特に近年は、分子計算や量子計算等の自然計算の枠組みでの計算機の立案・設計に

図2. 変換規則の適用場所による計算の分岐と合流。

ついても研究が進められており、グラフ変換による計算の基礎理論の必要性が高まっています。

直線状に並んだ頂点をつないだグラフの変形による計算は非線計算機による数式処理は数式を木構造で表現しパターンに適合

形の物理現象のモデルとして応用されているセルオートマトンによした部分木を書き換えることにより式の計算や変形を行います。数る計算と深く関わりがあります(図3)。セルオートマトンの計算機構

式中の同一の項を共有するグラフ構造で式を表現しグラフ変換をを直接実現する素子の可能性として、遺伝子等のDNA鎖を用いた行うことで計算効率を上げることが可能になります(図1)。そこで分子計算による実現が考えられています。グラフ変換系やセルオー

は、グラフ構造のパターン照合アルゴリズムや式の意味を変えない

トマトンの計算可能性、計算万能性、計算等価性の判定などの理論

グラフ表現方法が重要になります。

は、実現される計算機の能力の有効性を評価します。

図3. セルオートマトンによるフラクタル図形生成例
図1．
項を共有し、
グラフ変換で効率良く計算を行う。

グラフ変換の同時並列計算の定式化や計算の合流性の証明に
グラフ変換は変換計算領域に重複がない限り同時並列計算が可は、離散数学、二項関係による関係計算理論、そして、圏論などの代能です。交通網等における最短経路の計算、情報通信ネットワーク

数的手法が必要となります。また、計算能力の評価については、
オー

におけるネットワーク信頼度の計算、電気回路におけるインピーダ

トマトンや言語理論、計算量の理論、そして、数理論理学が必要となンスの計算などには、グラフ変換を利用した計算アルゴリズムが開ります。理論計算機科学は、現在の電子計算機を効率良く利用する

発されています。変換を利用した計算は変換規則の適用場所や適ための、データ構造やアルゴリズムの理論だけでなく、新しい計算用順序により計算過程が分岐するため、いつも同じ計算結果が得ら機構の立案・設計を支援するための理論でもあります。

れるとは限りません(図2)。その計算の分岐が必ず合流することを

現在の情報化社会を支える電子計算機の利用技術革新や次世

理論的に証明すること、あるいは、同等の計算を行う必ず合流する

代の新しい計算機構開発を支えるための数学、理論計算機科学の変換規則群を導くことが重要になります。

研究を行います。

7

Institute of Mathematics for Industry
Kyushu University

九州大学マス・フォア・インダストリ研究所

生物に学ぶ適応ネットワーク理論

手老篤史

学位：博士（理学）（北海道大学）
専門分野：数理モデリング、適応ネットワーク理論
生物は長く過酷な生存競争を生き延びてきただけあって、生命現究者の中垣俊之教授（はこだて未来大学）
により行われました（図２

象の内には巧みで洗練された技術が多くあります。私の研究ではこa-c）。
しかし脳が無く自分の周りの情報しか持たない粘菌がどのよのような生命現象を数式で書き表すことによって生命の持つ技術をうにして大域的な情報が必要なネットワーク問題を解くことができ

理解・抽出し、工学・工業的な応用へと繋げる事を目標としています。

たのでしょうか。

鉄道網や蟻の列、血管網に葉脈など、輸送に関する様々なネット

私はこの現象を数式で表すことによりこの現象を再現しました

ワークがあります。これらは全て使われている経路が発達し、使われ（図２d-f）。その結果、
この最短経路解が求まるパラメータを境界にて無い経路は縮退していく性質があります。このようなものを適応してネットワークの形が大きく変わることがわかりました。適応ネッ

ネットワークと呼びます。これらの適応ネットワークの形状は状況に

トワークの成長率が太い経路ほど強かった場合、図３aにあるように応じて最終的な形状が（毛細血管や大動脈のように）大きく異なり

経路は１本だけ残ります。これは初期状態で太い管が成長しやすく、ます。このような適応ネットワークの形成を理解することが本テーマ

その後、より一層成長しやすくなっていくからである。反対に太い経の目的です。

路ほど維持コストが高ければ太い管は細い管に血流を奪われてしまい、全ての経路が残ります(図３c)。このように粘菌が迷路を解くという研究の私の数理モデルはこの２つのネットワーク形成の境界

図1を与えているのです。このためネットワークトポロジーが変化する場
合には必ず図３bの状況を通過するので、適応ネットワークの組み換えを知るには図３(a-c)のパラメータに注目すれば良いのです。
図3
（a）

（b）

（c）

私がこの研究に用いた生物は真正粘菌変形体という生き物です（図１）。粘菌は単細胞生物ですが、体内に栄養等を輸送する適応

初期値依存で1本残る

ネットワークを持っています。また、内部に核をたくさん持った集団

最短経路が選択

全ての経路で残る（毛細血管型）的な性質も持った生き物です。例えばナイフで幾つかに切断すればそのそれぞれが別個体として生存できるし、反対にくっつけてしま

このようにして適応ネットワークの共通法則が得られたため、それえば１つの個体として生活することができます。粘菌はこのようにを一般の鉄道網に応用しました（図４a-c）。これからはさまざまな適

切ったり貼ったりが自由にできることから適応ネットワークを理解す

応ネットワークに対してこの共通法則を応用し、工学・実用的な結果るのに優れた素材になっています。この粘菌が作る輸送ネットワー

へと繋げていくことを目標とします。

クは迷路を解き、最適なネットワークを求めるという実験が共同研

また、他にも生物の多様なリズムを用いた行動制御や生物の自発的な形態形成についての研究も行っております。

図2

図4

8

◯ 応用理論研究部門

最適化を基本とした研究・教育活動を目指して脇隼人

学位：博士（理学）（東京工業大学）
専門分野：最適化理論、連続最適化
自然界のみならず、社会活動や日常生活においても、「最学生指導においては「分野の縮小再生産はしない」と考え適化」が現れるので学ぶ価値のある研究分野といえよう。最ていて、できるだけ自分の守備範囲の外の最適化について

適化の中でも、連続量を扱う最適化を連続最適化と呼び、私研究してもらっている。例えば、整数計画問題の応用およびは主に連続最適化を対象として研究や教育を行っている。

そのためのアルゴリズムの開発などがそれにあたる。幸い、

例えば、大学一年生の微分積分学で習う「ラグランジュの未学生が優秀なためこちらがいろいろと教わって論文を一緒

定乗数法」は連続最適化で議論される手法の一つである。

に書いている。

以下では、私が携わった研究・教育および産学連携について

理学部数学科では最適化を教える講義がないが、大学院

簡単に紹介する。

では用意されている。こちらもできるだけ幅広くテーマを選んでいて、最適化理論の基本を重点的に教えている。

これまで凸最適化、特に、行列変数を持つ半正定値計画問題と呼ばれる最適化問題を中心に研究をしてきたが、IMIに

これまで、IMIのシステムを利用して、自動車業界、製造来るまでは、半正定値計画問題を使って別の最適化問題を業、行政などと幾つかのテーマで研究費をもらって共同研

解くアルゴリズムを提案・実装する研究に取り組んでいた。

究を行った。いわゆる産学連携である。私の方針としては、

現在は、もう少し半正定値計画問題に関して理論寄りの研学生に補助はしてもらうができるだけ一緒に手を動かすよう究にシフトしている。具体的には、制御理論で現れる最適化にしている。学生を安価な労働力と捉えられないようにする問題に対して、得られる半正定値計画問題をテーマに研究ためだ。

を行っている。この半正定値計画問題はしばしば悪条件にな

このように取り組んでいると、産学連携についてはいつもることがあるのだが、それがシステム論の言葉で記述できる

難しさを感じている。その一因は、「論文を書いて自分たちことに気づき、それ以来このテーマに取り組んできる。調べの獲得した知を後世に残す」という共通認識を持てないとこてみると双対問題が興味深い数学的構造を有しており、幾ろにある。我々にとっては当たり前の思想であるが、相手側

つか論文を書いている。
「最適化理論の深化よりも、他の研は必ずしもそうではない。産学連携を通じて産業界・教員お究分野でどうやってうまく最適化理論を使うか」、ということよび学生が論文を書き、その学生が学位を取り、学術界・産を考えることが性に合っているようだ。

業界で活躍してもらうのが産学連携の理想だと考えている。
この理想に同意してもらえる相手と優先的に共同研究を実施しようと思っている。

9

Institute of Mathematics for Industry
Kyushu University

九州大学マス・フォア・インダストリ研究所

計算理論と整数論の間

浦本武雄

学位：博士（理学）（京都大学）
専門分野：理論計算機科学
私の主要な研究領域は、広義には理論計算機科学です。

興味を持たれていました。Schützenbergerはこの特徴付

その関連で最近は特に代数的言語理論と古典類体論の融けを、有限半群論の手法を用いて与え、それがEilenberg理合領域の研究に取り組んでいます。

論の出発点となります。Schützenberger以降も正規言語に関する組合せ的・数理論理学的決定問題が、有限半群論（1）代数的言語理論の問題に帰着させて解決できる場合のあることが発見され代数的言語理論は形式言語理論の一領域で、歴史的にはました。その方法を体系化したものがEilenberg理論です。

言語（＝有限文字列集合）の計算階層分類の関連で1960

現在では以下の図のように正規言語の様々な階層は、有限年代から発展してきたものです。例えば我々は素朴に言えば半群（モノイド）の階層と自然に対応することがわかっていま「２進数1101110110101が素数であるか」を判定す

す
（双対性）
：

るのが難しいと感じ、一方で「1101110110101の中に含まれる0の個数は3の倍数か」を判定するのは容易である

regular languages

finite monoidsと感じますが、それは何故でしょうか？ 形式言語理論（或はpiecewise testable

J-trivial

FO definableaperiodic

計算複雑性理論）ではこのような「有限文字列の性質を判定する問題の難しさ」に正確な（Turing機械の実行ステップ数などによる）定義・指標を与え、その難しさの階層を分類することを一つの主要な関心としています。

図2：正規言語と有限モノイドの階層の対応（双対性）（2）Galois理論・古典類体論との遭遇
Schützenbergerの定理は、技術的には次のように主張
します：「正規言語LがFOで定義できるための必要十分条
件は、Lから構成される有限統語的モノイドM（L）が非自明な部分群を含まないこと（= aperiodic）である。」簡単にいうと、正規言語がFOの論理式という特定の表現形式で記述

図1：バイナリ列に含まれる0の個数が３の倍数かを判定する
有限オートマトンできるための必要十分条件を、それに対応する有限モノイド
M（L）の純半群論的性質で特徴づけられるということです。

代数的言語理論はこのような言語の研究の中でも、代この現象は代数学においてGaloisが示した次の主張と似数的な手法を援用することを特徴としています。特に有限ていないでしょうか：「多項式f（z）の根が四則演算とベキ根

オートマトンによって受理される言語である正規言語に限で表されるための必要十分条件は、（z）fの最小分解体Kの定して言えば、
有限半群論の手法を援用した美しい分類理論

ガロア群Gal（K/Q）が可解であることである。」実際、
この（Eilenberg理論）が知られており、私の研究もこの正規言類似はただの類似ではなく、圏論と呼ばれる抽象的な言葉

語の分類理論に関わります。

を使って正当化することができ、その結果としてEilenberg

Eilenberg理論の内容を端的に言い表す良い例は、

理論とGalois理論とは正確な意味で統一することができま

正 規 言 語 の 論 理 的 表 現 に 関 する問 題 で 、1 9 6 5 年 にす。そしてさらにこの統一は、代数体のアーベル拡大の分類

Schützenbergerが証明した定理です。正規言語は一般

や、
アーベル拡大での素数・素イデアルの分岐を記述する類にBüchiのMonadic Second-order Logic（MSO）
とい体論に対しても、興味深い見方を与えることがわかってきまう論理体系の論理式で表現できることが知られていました

した。現在はこの観察を進め、古典的な代数的言語理論のが、MSOのうち一階の量化子のみを使う部分体系（FO）でアイデアと、明示的類体論や非可換類体論との関連づけを表現できる正規表現を特徴づける問題が1960年あたりに模索しています。

10◯ 応用理論研究部門

双曲型偏微分方程式の逆問題髙瀬裕志

学位：博士（数理科学）（東京大学）
専門分野：偏微分方程式、逆問題、幾何解析
双曲型偏微分方程式は解の有限伝播性やホイヘンスの原理等重要な数学的特徴を持つ。数理科学の分野での重要性はもちろんのこと、力学や波動をはじめとする様々な自然現象のモデル化が可能であるため、これまで物理学や工学の分野においても重要な役割を担ってきた。中でも、存在す
る解の情報から支配される方程式系を決定する逆問題は、数理科学の分野では非適切性問題として極めて重要であ

図２：コンパクト多様体と共形コンパクト多様体

る。加えて物理現象の背景を詳細に理解できるツールとしても強く注目されている。

結果の原因となる未知量を特定する逆問題は、当然なが偏微分方程式の逆問題における代表的な研究課題は、未ら必ずしも肯定的に解決されるわけではない。否定的な例

知の波源を特定する波源項決定逆問題と、方程式が記述すとして、映画「ハリー・ポッター」に出てくる透明マントのようる未知の物理的性質を特定する係数決定逆問題である。考な光のクローキング技術がある。これは、ある対象の周りのえている領域の一部で波動を観測した観測者が、これらの有界な領域において媒質を変化させることで光の進行を曲

未知量を観測領域の近く或いは考えている領域全体で特定

げ、外部からの光が伝播しない見えないゾーンを作り出すこ

できるかどうかを調べることが目的となる（図１）。

とでその対象を不可視化する技術である（図３）。偏微分方
程式論において媒質の変化はしばしば低階項であるポテン
シャルとして定式化される。見えないゾーンにおいてはそのポテンシャルの影響により外部から光が伝播しない状況を、
ポテンシャル項付き波動方程式のコーシー問題の観点から理解することもできる。このように否定的な状況を改めて考図１：逆問題察することが数学的に面白い構造を発見する上で助けにな

様々な偏微分方程式のタイプがある中で、特に双曲型偏ることもあるのかもしれない。

微分方程式は座標に依存しない幾何学的な記述が可能である場合が多い。
したがって例えば不定値な計量を導入した擬リーマン多様体上における解析との相性が非常に良い。近年では逆問題の様相を幾何学的視点から捉えるべく、
コンパクトなローレンツ多様体上の幾何解析を用いた研究が進められている。さらに、物理学からの要請として多様体の境界で計量が発散する共形コンパクトな多様体上での逆問題解析も注目を集めている（図２）。共形コンパクト多様体上では方程式の主要部が退化するため解析が容易ではな
く、退化型方程式に対する新たな逆問題理論の構築が必要

図３：クローキング

である。

11

Institute of Mathematics for Industry◯ 基礎理論研究部門

Kyushu University

九州大学マス・フォア・インダストリ研究所

代数解析：数学の落ち合うところ落合啓之

学位：博士（数理科学）（東京大学）
専門分野：代数解析
私の研究分野は代数解で再び離任することと析学で、特にD加群を用いなった の は 残 念 です

た研究を主としています。

が、「転んでも起きな

初期の頃、修士論文ではい」という不屈の精神

ある非ホロノミー系を満で新地へ向かいたい

たす佐藤超関数を調べ、と思います。私の今ま学位論文ではあるホロノでの活動のうち、社会ミー 系 を 満 た す 半 単 純との接点と言えるのは、数学が伝統的に社会と向き合ってきた諸分

リー群の大域指標の定め野のうちのいくつかですが、数学の社会への情報発信や連携は機る輪体を調べました。もう会をいただいてつとめています。例えば、高校教育（のうち、おもに少し先を説明するために、

入試問題説明会や教科書の編集に関すること）、保険数理（アクチュ

私の経歴を簡単にまとめてみます。1998年、それまで９年間お世アリ会の研究会員）、執筆（数学セミナーや数学小景などでの数学話になった立教大学を離れ、九州大学に赴任しました。
３年後に東工の平易な紹介）、講演（公開講座、市民講座）などです。昨年度から安大へ移り、さらに2002年の秋に名大へ異動、そこで７年を過ごし、生健一氏を代表者とするクレスト（科学技術振興機構）の『デジタル2009年の秋からまた九大に戻ってきました。この間、概均質ベクト

映像数学の構築と表現技術の革新』の九大班の班員もつとめていル空間、１変数ならびに多変数の超幾何関数を始めとする特殊関ます。

数、そしてゼータ関数と出会います。私自身は佐藤幹夫氏の直接の近年は幾何学的表現論の枠組みで積分が果たす役割を中心に研指導を受けたことはないのですが、無意識のうちにその影響下にい究を進めています。日本数学会では函数解析分科会の「表現論と調ると思っています。ゼータ関数も研究の対象にしていますが、強い和解析グループ」に属していますが、必ずしもその範囲に限らないて言えば、私自身の興味は、数に対するものよりも関数に対するも

対象や手法にも興味を持って活動して来ました。今後も柔軟に数学のといえるでしょう。今まで論文の対象としてきたものも、関数で表をして行きたいと考えています。

わせるものやそれらの間の関係式として表わせるものが多いと思います。そして、関数そのものとそれを決定・記述する様式の双方か
らアプローチしています。これが代数解析的な手法の特徴のひとつです。また、図形全般に対する興味はともかくとして、不変式論、
ならびにその双対と言える軌道分解には、研究の初期の頃からかなりの思い入れがあります。座標を使わない内在的な叙述と座標を使った
明示的な記述を、あえて片方だ
けで徹底したり行ったり来たりしたり、あるいは組み合わせ的技法を絡めたりという部分に面
白 み を 感 じ て い ま す 。今 やKazhdan-Lusztig

予想も常識のようになりつつあり、見渡せる景色は前世紀よりも確実に広くなっていると思います。
このパンフレットはIMI に関するものですから、私とマス・フォア・イ
ンダストリとの関わりについて説明します。2010年４月に、GCOEの事業推進担当者のひとりである岩崎克則氏が北大へ転任したため、岩崎氏のつとめていた「機能数理の基礎」ユニットの事業推進担
当者を引き継いだことで、マス・フォア・インダストリとのつながりが始まりました。このたび、IMIの発足に伴い、数理学研究院を１年半

12◯ 基礎理論研究部門

トポロジーとその応用

佐伯修

学位：博士（理学）（東京大学）
専門分野：位相幾何学、トポロジー、特異点論、微分位相幾何学、DNA結び目
位相幾何学（トポロジー）とは、図形をゴムのようなものでできてませんでした。1980年代終盤、数学者のC.ErnstとD.W.Sumnersいると考え、
グニャグニャと連続的に変形できるものは同じと思う、

は、数学、特にトポロジーにおける結び目理論の最新の結果を駆使

とても柔らかい幾何学のことです。言い換えると、幾何学的対象のしてその酵素の仕組みを解明しました。もともと結び目理論は19性質のうち、それを連続的に変形しても変わらないものを研究す

世紀のGaussによる電磁気学やKelvin卿による渦原子仮説に端る、純粋数学の一分野です。たとえばコーヒーカップとドーナツは位を発すると言われていますが、その後化学者・物理学者は結び目を

相幾何学的には同じものと考えられます。それらには「穴」がちょう

忘れ、数学者だけが興味を持って研究してきた歴史があります。

ど１つずつありますが、この「穴」の個数は、図形の連続変形で変わ

我々はその流れを受け、現代数学の中でも最近もっとも盛んに研らない量の１つの典型的な例となっています。

究されている分野の１つである結び目理論を駆使することにより、
トポイソメラーゼと呼ばれる酵素によるDNA組み換えの解析を、数学的側面から研究しており、こうした解析の産業技術への応用も目指しています。
なお、こうした研究と同時に、可微分写像の特異点論も活発に研究しています。特に、滑らかな物体間の写像の特異点が、そうした物体の位相幾何学的性

ドーナツとコーヒーカップは位相幾何学的には同じ質を深く反映する具体的

こうした位相幾何学は、図形が柔軟性を持つ場合に威力を発揮し事実などを数多く発見してます。その最たるものが、紐を結んでできる結び目や絡み目です。

きています。特に、特異ファ

紐を結ぶことは日常生活でも重要な行為で、結び目が原始時代からイバーと呼ばれる1点の逆

人間の生活と深く結びついていることは間違いありません。実際、像に着目する研究では世

4次元多様体の鍵を握る特異ファイバー

野生のゴリラも結び目を作ることができることが知られています。

界的な第一人者を自負しており、それを最初に定式化した著書も出

ごく最近になって、こうした結び目が、DNA（デオキシリボ核酸）の研版しました。最近はこうした理論を、多値関数データのための視覚究に深く関わることが明らかになってきています。

的解析に応用すべく、研究を行っています。

遺伝情報の担い手であるDNAは、生物の細胞内で捩れた紐状のまた、可微分写像の特異点論以外にも、位相的埋め込みの第一障形をしていて、輪になっていることもあり、
さらにそれが結ばれたり

害類、余次元1写像の分離性質、複素超曲面の孤立特異点の位相絡んだりしていることもあります。こうしたDNAの結び目・絡み目

幾何、
ファイバー結び目、4次元多様体、余次元1の埋め込み、空間は、酵素の働きによって作られることは知られていましたが、その仕曲線の微分幾何学的不変量、結び目解消数など、位相幾何学の組みの詳細については実験技術の限界のために明らかにされてい様々な分野を幅広く研究しています。こうした研究は、種々の分野で応用できる可能性があり、たとえば物質・材料の性質をミクロなレベルから考察する際には威力を発揮することが期待されます。なお、一般化されたフィボナッチ数列の漸近挙動についても精力的に研究を行っており、解析関数の零点に関する研究もあります。
こうした幅広い研究を行っていることは、教育に対する効果もかなり大きく、
これまでに指導した学生の産業技術への貢献実績もあります。また、異分野連携で活躍する数学博士人材を育成するための卓越大学院プログラムのコーディネーターも務めております。
今後も、トポロジーという純粋数学を基に、数学と産業技術の関わりを深めてゆければ誠に幸いです。

酵素によるDNA組み換えの例

13

Institute of Mathematics for Industry
Kyushu University

九州大学マス・フォア・インダストリ研究所

確率論とその応用
確率論とその応用

白井
白井 朋之
朋之
学位
学位
：
：博士
博士
（数理科学）（数理科学）（（東京大学）
東京大学）
専門分野
専門分野
：
：確率論
確率論
ルーレッ
ルーレッ
ト、
ト、
双六、
双六、
宝く
宝く
じ、
じ、
株価の値動きなど、
株価の値動きなど、
確率的現象は日常生
確率的現象は日常生

ランダ
ランダム
ムな
な要
要素
素は
は

活においても身近に観察できるものです。
活においても身近に観察できるものです。
これらの現象(対象)にはこれらの現象(対象)には

まっ
まった
たく
くあ
ありま
りませ
せ

自然に確率的な要素が備わっていますが、自然に確率的な要素が備わっていますが、
一見ランダムな要素が一見ランダムな要素が

ん。
ん。
しかし、
しかし、
多くの理
多くの理

まったくないように見える問題にもランダムな側面があったり、まったくないように見える問題にもランダムな側面があったり、
意外
意外

論
論的
的考
考察
察や
やシミュ
シミュ

な形で確率論的手法が使えることがあります。
な形で確率論的手法が使えることがあります。
ここではそのような
ここではそのようなレーションなどによレーションなどによ例を3つばかり紹介します。
例を3つばかり紹介します。

り、り、
その零点達があ
その零点達があ(1)日本では(1)日本では「角谷予想」「角谷予想」
、
、
海外では
海外では
「コラッツ予想、「コラッツ予想、
3n+1問題」
3n+1問題」
なな

るランダム行列の固るランダム行列の固どと呼ばれる初等整数論の有名な未解決問題をご存知でしょどと呼ばれる初等整数論の有名な未解決問題をご存知でしょう
う

有値と同じ振る舞い
有値と同じ振る舞いか?
か? 任意に選んだ自然数に対して、任意に選んだ自然数に対して、「「自然数が偶数ならば2で割自然数が偶数ならば2で割をすることがわかっをすることがわかっ

り、り、
奇数ならば3をかけて1足す」
奇数ならば3をかけて1足す」という操作を繰返すと必ず最という操作を繰返すと必ず最てきています。
てきています。

後は1になるという予想です。
後は1になるという予想です。
例えば、
例えば、
初めに7を選んだとする
初めに7を選んだとする

私の興味の一つは、
私の興味の一つは、
上
上

と、
と、
7、
7、
22、
22、
11、
11、
34、
34、
17、
17、
52、
52、
26、
26、
13、
13、
40、
40、
20、
20、
10、
10、
5、
5、
16、
16、

に挙げたような一見ラに挙げたような一見ラ

以下の自然
8、
8、
4、
4、
2、
2、
1となります。
1となります。
2009年の時点で20×2
2009年の時点で20×2 以下の自然

ンダムには見えない現象の裏に潜むランダムな現象を見つけて、ンダムには見えない現象の裏に潜むランダムな現象を見つけて、
そそ数に対しては、数に対しては、コンピュータによって予想が正しいことが確かめ
コンピュータによって予想が正しいことが確かめ

れを確率論的手法を用いて研究することです。
れを確率論的手法を用いて研究することです。

58
58

リーマンのゼータ関数の値の挙動:
リーマンのゼータ関数の値の挙動:

られています。
られています。
この問題自身にはどこにも確率的要素はありま
この問題自身にはどこにも確率的要素はありま
せんが、
せんが、
実はランダムな構造が潜んでいます。
実はランダムな構造が潜んでいます。
(2)以下の文字列は、(2)以下の文字列は、
単純換字暗号とよばれる古典的な方式で、単純換字暗号とよばれる古典的な方式で、あある文章を暗号化したものです。
る文章を暗号化したものです。
つまり、
つまり、
aをf
aをfに、
に、
bをtにというよう
bをtにというように、に、aからzの順番を入れ変えて一対一に対応させることによっaからzの順番を入れ変えて一対一に対応させることによって暗号化しています。
て暗号化しています。
fyeaxjqfzjoxceddedcjbujcxbjhxgpjbegxrjukjzebbed
fyeaxjqfzjoxceddedcjbujcxbjhxgpjbegxrjukjzebbed
cjopjsxgjzezbxgjudjbsxjofdljfdrjukjsfhedcj
cjopjsxgjzezbxgjudjbsxjofdljfdrjukjsfhedcj
暗号化の方法が上の方式であるとわかっているとすると、暗号化の方法が上の方式であるとわかっているとすると、
マルコ
マルコ

ランダム行列(GUE)の固有値の2点相関関数
ランダム行列(GUE)の固有値の2点相関関数

フ連鎖という確率論の基本的な道具で(ある程度)解読すること
フ連鎖という確率論の基本的な道具で(ある程度)解読することができます。
ができます。
以下はその実行例です。
以下はその実行例です。

(1)の問題は1960年代に数学者の角谷静夫先生が興味を持たれ(1)の問題は1960年代に数学者の角谷静夫先生が興味を持たれて、て、
イェール大学の同僚に伝えたところ、イェール大学の同僚に伝えたところ、イェールの数学教室ではイェールの数学教室では

0:0:fyeaxjqfzjoxceddedcjbujcxbjhxgpjbegxrjukjzebbedcjopjsxgjzezbxgjudjbsxjofdljfdrjukjsfhedcj
fyeaxjqfzjoxceddedcjbujcxbjhxgpjbegxrjukjzebbedcjopjsxgjzezbxgjudjbsxjofdljfdrjukjsfhedcj
1:1:ugado
ugadocut
cutbohassash
bohassashieiehoi
hoiponmponmianol
ianolek
ektaiiash
taiiashbm
bmfon
fontationtationesesifo
ifobusy
busyusl
uslek
ekfupashfupash

一ヶ月くらいの間皆がこの問題に嵌ってしまったそうです。
一ヶ月くらいの間皆がこの問題に嵌ってしまったそうです。
他の大学
他の大学

2:2:upico
upicofut
futbodissisd
bodissisdreredor
dorhong
hongrinol
rinolekektirrisd
tirrisdbg
bgmon
montitron
titronesesrmo
rmobusy
busyusluslek
ekmuhisdmuhisdでも同じような現象が起こり、でも同じような現象が起こり、一時期流行ったゲーム
一時期流行ったゲーム
「テ
「テ
トリス」
トリス」
と同
と同

3:3:avico
avicowat
watbohissishbohissishrerehor
horpong
pongrinol
rinolededtirrish
tirrishbgbgmon
montitron
titronesesrmo
rmobasybasyasl
asled
edmapishmapish4:4:aviceavicewat
watmepirrirp
mepirrirpsosopes
pesleny
lenysined
sinedob
obtissirp
tissirpmymyhen
hentitsen
titsenororsheshemarkmarkardardob
obhalirp
halirp

じように、
じように、数学の研究を鈍らせるための某国の陰謀だったのではな
数学の研究を鈍らせるための某国の陰謀だったのではな

5:5:apice
apicewaswasbeginning
beginningtotoget
getderyderytirel
tirelofofsitting
sittingby
byher
hersistersisterononthe
thebank
bankanl
anlofofhading
hadingいかという笑い話もあったそうです。
いかという笑い話もあったそうです。

6:6:alice
alicewaswasbeginning
beginningtotoget
getveryverytired
tiredofofsitting
sittingby
byher
hersistersisterononthe
thebank
bankandandofofhaving
having
7:7:alice
alicewaswasbeginning
beginningtotoget
getveryverytired
tiredofofsitting
sittingby
byher
hersistersisterononthe
thebank
bankandandofofhaving
having(2)であげた問題は、(2)であげた問題は、
アメリカのある州立刑務所の囚人の書いた壁
アメリカのある州立刑務所の囚人の書いた壁

8:8:amice
amicewaswasbeginning
beginningtotoget
getveryverytired
tiredofofsitting
sittingby
byher
hersistersisterononthe
thebanlbanlandandofofhaving
havingの落書きが何を意味するのかを調べるために、の落書きが何を意味するのかを調べるために、刑務所の心理学者
刑務所の心理学者

9:9:alice
alicewaswasbeginning
beginningtotoget
getveryverytired
tiredofofsitting
sittingby
byher
hersistersisterononthe
thebank
bankandandofofhaving
having

マルコフ連鎖の実行例
マルコフ連鎖の実行例（1000回毎の様子）（1000回毎の様子）

によってスタンフォー
によってスタンフォード大学統計学教室の技術相談窓口に持ちこまド大学統計学教室の技術相談窓口に持ちこま

Googleなどの検索サイ
Googleなどの検索サイ
トの人気ページのランク付けにもマルトの人気ページのランク付けにもマルれたことがきっかけだっただそうです。
れたことがきっかけだっただそうです。

コフ連鎖が使われていることはよく知られています。
コフ連鎖が使われていることはよく知られています。

(3)の研究の方向性は、(3)の研究の方向性は、
プリンス
プリンス
トン高等研究所のコモンルームでトン高等研究所のコモンルームで(3)2000年にアメリカのクレイ数学研究所が、(3)2000年にアメリカのクレイ数学研究所が、ミレニアム問題とミレニアム問題とのティータイムに、のティータイムに、物理学者であるダイソンから解析数論の研究者
物理学者であるダイソンから解析数論の研究者

して100万ドルの懸賞金をかけて7つの問題を提出しました。
して100万ドルの懸賞金をかけて7つの問題を提出しました。

であるモンゴメリーへ述べられた一つのコメンであるモンゴメリーへ述べられた一つのコメン
トから始まりました。
トから始まりました。

そのうちの一つであるポアンカレ予想は最終的にはペレルマン
そのうちの一つであるポアンカレ予想は最終的にはペレルマン

数学者レンニは、数学者レンニは、「「数学者とはコーヒーを定理にかえる機械である」
数学者とはコーヒーを定理にかえる機械である」

によって解かれましたが、によって解かれましたが、
いわゆるリーマン予想は未解決の難いわゆるリーマン予想は未解決の難と言ったそうですが、と言ったそうですが、一杯のコーヒー(紅茶?)が重要な研究の方向
一杯のコーヒー(紅茶?)が重要な研究の方向

問です。
問です。
リーマン予想とは、
リーマン予想とは、付けをした場面でした。
付けをした場面でした。

で定義されるリーマンで定義されるリーマンのゼータ関数とよばれる複素関数ののゼータ関数とよばれる複素関数の（非自明な）（非自明な）
零点はすべて零点はすべて

IMIにおいても、
IMIにおいても、自由な発想により実りある異文化間の交流がで自由な発想により実りある異文化間の交流がでという線上にあるだろうという予想です。
という線上にあるだろうという予想です。
この問題にもこの問題にもきることを心から願っています。
きることを心から願っています。

14◯ 基礎理論研究部門

数論的不変式論とそれにまつわる幾何石塚裕大

学位：博士（理学）（京都大学）
専門分野：数論、数論的不変式論、ディオファンタス幾何学
私は主に数論的不変式論と呼ばれる分野を中心として研究活動をしています。大枠としては代数学の中の数論に属する分野です。その動機と手法から紹介します。
（図1）
図2：線形行列式表示の例

関連した別の研究として、ある数論的な性質を満たす曲線などの具体例を構成したり、あるいは具体的に与えられた曲線の性質を計算する研究も行っています。
その一例が、平面四次曲線の双接線についての研究です。平面四次曲線の双接線は、その曲線と二点で接する直線、あるいは一点で四重に接する
直線を指します（図3）。平面四次曲線が有理数係数で定義されていても、双接線が有理数係数で定義できるとは限りません：たとえば虚数が係数に必要かもしれません。伊藤哲史
氏、大下達也氏、谷口隆氏、内田幸寛氏との共同研究では、
より強く双接線が局所大域性を 図3：ある平面四次曲線（青）の双接線（緑）
満たさない有理数体上の平面四次曲線を構成しました。実数やp進数の範囲では双接線が見つかるものの、有理数係数では見つからない曲線です。

図1：数論的不変式論の議論の流れ

まず興味の対象は、代数体やそのイデアル類群、（有理数体上の）楕円曲線のモーデル・ヴェイユ群など、数論において重要だが扱いの難しい対象です。代数体や楕円曲線は無限に存在し、それぞれ例を作れと言われると比較的簡単に
構成できます。
しかし個々の例に対してイデアル類群やモーデル・ヴェイユ群を計算することは、
しばしばかなり難しくな
ります。そしてパラメータを少し変更するとそれぞれの量は大きく変わります。
そこで、個々の例ではなく、代数体や楕円曲線全体を考えたときにイデアル類群などがどう振る舞うか、という問題を考えます。たとえばイデアル類群の『平均』位数や、モーデル・ヴェイユ群の階数が大きい楕円曲線の『割合』などを問題にします。こうした問題意識を持つ研究は総じて数論統計と呼ばれ、数論幾何、表現論、ランダム行列の理論などと結びつきながら発展しています。数論的不変式論の主な結果である M. Bhargava と A. Shankar の結果も、
モーデ
ル・ヴェイユ群の（高さに関する）平均階数の上界を与えるものです。
次にその手法です。まず数論的対象を、より取り扱いやす
い代数群の線形表現の軌道として解釈します。この解釈を踏まえると、数論的対象の数え上げの問題が、ある基本領域
中の格子点の数え上げに帰着される場合があります。そうした場合には、解析数論の技術と結びつけることで、もとの数論的対象についても平均・割合などが議論できるのです。私もこの流れを踏まえて、有理数体上の平面三次曲線のなかで、線形行列式表示と呼ばれる表示を持つものの割合を研究しました（例：図2）。実際には平面三次曲線で、ヤコビ多様体が正の階数を持つものの割合も調べています。

こうした方向性では、平面三次曲線の変曲点についても同様の考察を行っているほか、フェルマーの四次曲線についてその法を四としたガロア表現の決定も行っています。
手法としては近年の結果を利用した数論幾何的な考察とともに、代数的な計算に数式処理ソフトを活用することで結果を得ています。
これらに通底するのは、代数幾何や不変式論の古典的な問題意識が、数論的な設定で新しい意味を帯びる点です。た
とえば数論的不変式論で用いられた数論的対象と表現の軌
道の対応の一部は、複素数体だと非常に単純な対応になってしまい、構造が却って見にくくなります。私はそういった数論的な設定で現れる構造について、数式処理の助けを借りながら研究を続けています。

15

Institute of Mathematics for Industry
Kyushu University

九州大学マス・フォア・インダストリ研究所

数理科学的観点からの血糖値管理小谷久寿

学位：博士（数理学）（九州大学）
専門分野：位相幾何学、整数論、数論的位相幾何学
ここのところ、私は他大学の医療分野の研究者や応用数

そこで、この共同研究では、数理的なアプローチと医療現学者の方たちと血糖値管理に関する分野の枠を超えた共同場での知見とを組み合わせ、実際の医療現場で役立つよう研究に取り組んでいます。もう少し具体的に書くと、
これはな血糖値管理をアシストする標準的手法の開発や熟練看護

手術後集中治療室（ICU）患者のデータを用いて、現実の医師の投薬判断アルゴリズムの推定を行うことを一つの目標

療現場に役立つような血糖値を管理するアルゴリズムの開としています。

発を目指した研究です。この研究は上に掲げている私の専門分野のお話とは異なるものですが、ここでは産業社会への数学・数理科学の応用を目指したIMI関連の研究として、

血糖値・インスリン
データ

それについてお話しします。

血糖値管理支援システム
投薬判断アルゴリズムの推定

手術直後ICU患者は手術侵襲によるストレスや強心剤等による影響で血糖値の急上昇が見られます。高血糖は多臓
器不全・昏睡・予後悪化等を引き起こす原因となるため、
イン
スリン投与により血糖値を適正値内に維持することが重要

このような問題に対する数理的なアプローチとして、これであると考えられています。
しかし、
インスリンによる血糖値までにも血糖値管理に関連する血糖値・インスリンの数理モ

管理はインスリン作用の遅延やインスリン感度の時間変化

デル等が考案されてきましたが、生体における動態を表現

等の様々な要因によりコントロールが難しく、さらにインスする数理モデルには観測できない変数が多く含まれていた
り、患者ごとに固有のパラメータを推定する必要がある等、難しさがあったり、それぞれできることや目的が異なっています。また、近年データ科学分野で様々な新しいデータ駆動型の研究手法が開発され続けています。この共同研究では、従来の手法と新しい手法を取り入れつつ、今回の研究の目的に即したより良いアルゴリズムの開発に貢献すべく努めています。
血糖値管理のような医学的知見が重要である共同研究においては、数学・数理科学的な立場だけからではなく、外科
医師、ICU専門看護師を交えて協働的に研究を進めることが重要です。医療従事者の方々とのコミュニケーションを大事にしながら着実に共同研究を進めて、数学・数理科学の医リン投与により逆に低血糖に陥ってしまった場合には深刻

療分野・産業社会への応用へと繋げていきたいと思います。

な予後不良を引き起こしてしまうという問題があります。現在、一定の条件のもとICUでの血糖管理は看護師によって行われていますが、標準的な管理アルゴリズムは未だに確立されておらず、血糖値のコントロールは看護師の経験に依る部分が大きいのが現状となっています。ただ、そうした
血糖値管理は看護師の負担も大きく、手法の標準化が望まれています。

16◯ 基礎理論研究部門

4次元図形の設計図：視えない次元を視る濵田

法行

学位：博士（数理学）（九州大学）
専門分野：低次元トポロジー、曲面の写像類群、シンプレクティックトポロジー
私は、低次元トポロジーとよばれる純粋数学の一分野を研究しています。より詳しくは、４次元の図形を曲面の上に“描く”
理論があり、
この４次元と２次元が相互に影響し合う関係に強い魅力を感じ、研究しています。
４次元を描く方法の説明の前に、次元を落とした簡単なモデルを見てみましょう。次の図のように細長い長方形の紐を用意して、両端を張り合わせることを考えます。張り合わせ方は二通りあり、一方からは普通の輪っか（アニュラスといいます）ができ、もう一方からは捩じった輪っか（いわゆるメビウスの帯）ができます。このとき図のように、アニュラス、
メ
ビウスの帯の下に仮想的な円周をそれぞれ置くと、円周の各点の上に短い線分が並んだような構造があることが見てとれます。円（１次元）の上に線分（１次元）が並んで、全体としては１＋１＝２次元の図形となっているわけです。

曲面の上に（特異）曲面が並んだ空間

こうして４次元多様体と曲面の写像類群が関係づき、
４次元多様体の制約から写像類群の性質が得られたり、逆に写像類群の議論から４次元多様体の性質が得られたりと相互補完的な研究があります。
私がもっとも得意とするのは、上で説明した写像類群における「設計図」を作ることです。これまでに、基礎的な例をみつけたり、再利用が可能な使い勝手の良い設計図を作った
り、またそれらをもとに複雑な設計図を作ることで存在が知られていなかった４次元多様体を構成したりしてきました。

アニュラスとメビウスの帯

想像力を膨らませてこの考え方を発展させると、曲面（２次元）の上に曲面（２次元）
が並んだ空間というものが考えられ、この空間の次元は２＋２＝４次元になります。さらにもう少し条件を緩くして、次の図のようにときどき特異点をもった曲面があらわれることも許すことにすると、非常に豊かなクラスの４次元図形（専門用語でシンプレクティック４次元
多様体といいます）
を扱えるようになります。
上に並んでいる曲面をファイバーといいますが、ファイバーの上に特異点に対応する閉曲線を描けば、考えている４次元多様体の情報をすべてもった「設計図」が得られます。もう少し正確には、
この設計図はファイバーの曲面に付随する写像類群という群の言葉で記述されるものです。

4次元多様体の設計図の実際の例

ところで４次元多様体の分類の仕方として、２つの多様体が“連続的に”
移り合うとき同じとみなす同相の立場と、“滑らかに”移り合うとき同じとみなす微分同相の立場があります。この２通りの区別の違いは非常に繊細であり、４次元トポロジーでは中心的な話題です。近年の私の（共同）研究では、ある標準的な４次元多様体に同相であるが微分同相ではない４次元多様体というものを多く発見しており、４次元における同相と微分同相のギャップについての新たな知見を得ています。

17

Institute of Mathematics for Industry◯ 数理計算インテリジェント社会実装推進部門

Kyushu University

九州大学マス・フォア・インダストリ研究所

トポロジーとその諸科学への応用

鍛冶静雄

学位：博士（理学）（京都大学）
専門分野：位相幾何学、
リー群、応用数学
私は図形をおおざっぱに調べるトポロジー（位相幾何学）

上にあげた特徴を眺めてみると、細かな差異に敏感な反という分野の研究をしています。どれくらいおおざっぱかと面、全体像を捉える直感が働かないコンピューターと対極

いうと、ドーナッツとマグカップの形が見分けられないくらい的に見えてきます。実際にトポロジーは、人間のもつ俯瞰的です（佐伯修教授の研究・技術カタログを参照）。数学に限らなものの見方を定式化し、コンピューターにその能力を授ず多くの学問は、対象を分類することを一つの目標としていける可能性を秘めています。

ます。どんな性質に着目するかによって分類の基準が変わっ

私は、空間のトポロジーを代数的に調べることに興味があてきて、猫と犬というように大雑把に分けることもあれば、ります。主にリー群や等質空間といった対称性の高い空間

秋田犬と柴犬と言うように詳しく分類することもあります。
トの、位相不変量を定義・計算することを専門としています。まポロジーでは図形、
より数学的にいうと空間を対象とします

た、代数はコンピューターの得意とする対象であり、応用とが、おなじみの合同や相似よりずっと荒く、連続的に変形でして図形をコンピューターで処理するアルゴリズムの研究

きるものは区別しません。そんな"いいかげん"なことで役ににも携わっています。一例として、下図にペンギン算をお見立つのかと思われるかもしれませんが、実はさまざまな場面せします。形を代数化することで、異なった形状を足し合わで活躍します。

せたり平均をとったりすることが可能になります。また様々な形の"地図"ができ（座標が導入され）
ることで、
どの形はど

まず、人間の視覚や空間把握は非常に大雑把ですが、多くの形に似ているといった解析が可能になります。この考え方の面でまだ機械を凌駕しています。細かな変化に気がつくを応用して、
コンピューターで形状をデザインする方法の開のは苦手な反面、手書きの歪んだ文字や、笑顔でも泣き顔

発や、医用画像・センサーデータの解析も行なっています。

でも同一人物を認識することができます。また全く未知の図形、いわゆるビッグデータなどの高次元や時に無限次元の空間を扱う時、
まずは大局的な情報が重要になります。新発見の島を調べるのに、いきなりルーペを持ち出してもあまり得るところはありません。それよりも高台に登っておお
まかな形を知りたいでしょう。
もう一つ別の観点からは、柔らかい性質は適用範囲が広い、というのも利点に挙げられます。硬い性質を扱う数学的事実、例えば、球の体積の公式は、真に球なものにしか適用できませんが、世の中に存在するものはイデアからはかけ離れていて、常に誤差が生まれます。一方で、代表的な位相的性質である
“穴の数”（ 正確にはベッチ数）は、歪んだ図形でも問題なく扱うことができます。
3羽の黄色いペンギンをもとに、足し合わせによってペンギンのバリエーションを生み出す（アルマジロのモデルは The Stanford 3D Scanning Repository より）

18◯ 数理計算インテリジェント社会実装推進部門

物質構造解析の数理

富安（大石） 亮子
学位：博士（数理科学）（東京大学）
専門分野：応用代数・数論、数理結晶学、アルゴリズム
物質構造解析研究は、調和解析・信号処理・最適化・統計学など
様々な手法を用いたアルゴリズム開発としての側面がありますが、私の場合、新しい知見を得るのに数論・代数学といった純粋数学のアイデア・技術を組み合わせることをよく行っています。
以下では、学術論文に加え特許出願に至った（したがって割と非数学者向けの）3つのプロジェクトを紹介します。
（１）‌格子決定（Ab-initio indexing）に関わるCONOGRAPH の方法“Ab-initio”とは専門用語で、物質構造（この場合は結晶格子）

図1: SDRを用いた構造解析法の概要

に関する事前情報を用いない解析のことです。以下の手法開発を行い結晶学分野でそれぞれ論文を書いた後、粉末回折用［1］・電子（3）‌黄金角の方法の一般化・高次化

線後方回折用プログラム［2］
を開発、web上で配布しています。開葉序・ひまわり頭部に見られるパターンのモデリングに使用され発には当時所属していたKEKの研究室、共同研究を実施していた

る黄金角の方法の一般化は様々な文献で試みられ、取り上げられ（株）
日本製鉄からの支援を受けました。

た未解決問題でしたが、マルコフ理論の高次元版であるproduct・‌大きな観測誤差下での格子対称性（ブラベー格子）の決定…格of linear formsと呼ばれる数の幾何・数論の問題に帰着させるこ子基底簡約の応用とで一般曲面・一般次元に適用可能となりました。

・ピークサーチ・‌実験データによくフィットする解の検出に用いるfigure ofmerit・消滅則の一般的性質の導出（トポグラフを用いて記述）
・‌解析の解の一意性（Ambiguity）の高速判定法…数論の2次形式論の応用開発したソフトウェアCONOGRAPHは解析の成功率を向上させました。得られた数学の結果は様々な回折データからの格子決定に対応できます。
（2）半正定値計画緩和法（SDR）に基づく結晶・磁気構造解析の方法
非線形最適化におけるローカルミニマムの問題はよく知られていまず。SDRは双対定理で保証された状況下において2次計画問
題の大域的最適解を求める方法です。一般に、結晶構造のフーリエ

図2: 適用結果（3Dは断面図、着色は各点近傍のパッキング密度による）変換すなわち「構造因子」の偏角を求める位相回復の問題（絶対値は観測値から得られる）は2次計画問題として表現できます。

モデリング・パターン生成・メッシュ生成において今回開発した方もともと未知構造解析の解の一意性を調べるために始めた研究

法・理論の応用を展開する予定です。

ですが、偶然、実験家との共同研究で役に立つ出口（磁気構造解析［1］https://z-code.kek.jp/zrg/［3］）も見つかりました。SDRの優れた性能は物質構造科学分野［2］J. Appl. Cryst. (2021) 54 (2), 624-635

で代数計算・離散数学の応用を進めるためにも有用と考えていま［3］Scientific Reports (2018) 8:16228.

す。

19

Institute of Mathematics for Industry
Kyushu University

九州大学マス・フォア・インダストリ研究所

グラフ解析と最適化問題の高速計算及び実社会への応用藤澤克樹

学位：博士（理学）（東京工業大学）
専門分野：最適化問題、グラフ解析、高性能計算
新しいスーパーコンピュータの応用として大規模なグラフ解析が注目を集

従来は計算中心であった高性能計算(ハイパフォーマンスコンピューティン

めている。グラフは点集合と枝集合から構成される。具体的には図1のようにグ)の分野においても大規模なデータ処理を中心に扱うアプリケーション

様々な応用分野おいて解析の対象とする事象の関係(Relationships)を点(データ・インテンシブアプリケーション)が増 加している。G r a p h 5 0 0と枝で表現していく。例えば道路交通ネットワークでは点は交差点、枝は交差(http://www.graph500.org)は並行探索、最短路探索をはじめとする最点間の道路に該当する。またTwitterなどのソーシャルネットワークの解析で適化、極大独立集合などのグラフ解析、などの複数のグラフ処理カーネルかは、点はユーザ、枝はユーザ間のフォロー関係(あるいはメッセージ送信)など

らなるベンチマークにより計算機の性能を評価しランキングを行う。
グラフ解に関連させることが多い。さらに各枝を連結させてグラフを構成して(Step析はサイバーセキュリティ、創薬、
データマイニング、
ネットワーク解析などの

1)、
目的に応じて最短路検索などのグラフ解析を行う(Step 2)。またグラフ

分野において必要とされる重要な計算カーネルとして位置づけられている。

解析の結果は元の応用問題の分析や理解のために使用される(Step 3)。実際にカーナビゲーションシステムでは道路ネットワークがグラフデータとして内蔵されていて、ユーザの指示に応じて出発地点と目的地点間の最短路検索を行っている。このように社会における実データをグラフデータに変換して、計算機で高速処理する需要が非常に高まっている。

図３：Graph500 の技術を応用した Twitter ネットワーク解析
このGraph500に対するBFSの高速実装は図3のようにTwitterネット
ワークの解析等に用いることができる。図3ではTwitterのユーザとフォロー
図１：グラフ解析の利用方法と応用分野関係を表したFellowship network 2009(点数:4100万、枝数24億)を用

グラフ解析に関しては数年以内に数千万規模の並列性を備えたポストペタ

構築している。これによって、あるユーザから何ホップ以内に何人のユーザがスパコン上での高性能な超大規模グラフ処理技術が確立され、ストレージの存在しているかについて高速に探索することが可能になる。この場合では24

階層性が深化したポストペタスパコン上での高性能な超大規模グラフ処理技

億枝のグラフに対してわずか0.069秒でBFS木の構築に成功している。

いて特定のユーザからBFSを行い、そのユーザを根(root)とするBFS木を術が開発されていくと予想されている。これによって防災計画の策定、災害時

また最適化問題の高速計算と実社会への応用にも取り組んでおり、例えばの避難と誘導及び情報収集と解析、
スマートグリッドによる高度かつ安定な電半正定値計画問題(SDP)は組合せ最適化、
システムと制御、データ科学、金力供給など、安全安心な社会基盤実現に貢献することが可能となる。以下の融工学、量子化学など非常に幅広い応用を持ち、現在最適化の研究分野で最分野に対してグラフ解析を適用することを想定している(図2)。

も注目されている最適化問題の一つとなっている。SDPに対しては高速かつ

1.交通データに対する経路探索：動的に変化する交通量等から高速な重要度

安定した反復解法である内点法アルゴリズムが存在しているが、巨大な線形方程式系の計算(行列要素の計算と行列のCholesky分解)が大きなボトル判定を行うことによって、交通管制等に活用する。
2.ソーシャルネットワークデータ(マイクロブログやSNSなど)やウェブデータ

ネックとなっている。最近の結果では多数GPUの活用や計算と通信のオーに対する動的な重要度、影響度の判定。各点の周辺、及び広域内におけるバーラップ技術を応用することによって、主要なボトルネックの1つである線形方程式系のCholesky分解の高速化と世界最大規模のSDPを高速に解く

影響（情報の伝播力）を推定する。

ことに成功した(最大で1.713PFlopsの性能を達成 : 図４)。

3.その他：疫病の拡散、人口の増減、経済動向等の分析。ライフライン等の基盤計画（電力、水、食料）。生命科学系（創薬、遺伝子）。ビジネス系（金融、
データマイニング）。安全保障分野（組織構成の解明、事件事故の事前予測）。

図４：半正定値計画問題に対する大規模並列計算

図２：大規模グラフ解析とその応用

20

◯ 数理計算インテリジェント社会実装推進部門

圧縮性回転粘性流体の数学解析藤井幹大

学位：博士（数理学）（九州大学）
専門分野：流体力学に現れる偏微分方程式の数学解析
大気や海洋などの大規模な地球流体の運動では、地球の自転の影響によるコリオリ力が流体の運動に作用することが特徴的である。実際、流体の基礎方程式であるNavierStokes方程式を回転座標系の問題に書き換えた際に流体の加速度を記述する項にコリオリ力を表す異方的な線形項が現れる（図１）。流体が密度変化しない「非圧縮性回転流」
は物理的にも自然な状況であり、この場合では数学的にも多くの研究が行われてきた。特にコリオリ力は線形解に分散性の効果をもたらすため、流体のエネルギーには一切影響しないにも関わらず、大きい初期値に対して回転速度を十分速く設定することで時間大域解を構成できることが知ら
れている。

図2

今後の研究では、木星の大赤斑（図２）の解析などにも用いられ応用上重要な回転浅水波方程式の解析に取り組む。
上記で述べた圧縮性回転Navier-Stokes方程式の研究は３次元流体であったが、回転浅水波方程式は２次元流であ

図1るため低周波での非線形評価がさらに難しくなる。さらに、
一方で、流体の密度変化を考慮した「圧縮性回転流」の場

回転速度を無限大とする特異極限を考えると、３次元の場合合の数学的な研究はあまり行われていないようである。地とは異なった振る舞いをすることも期待され多くの興味深

球流体力学では非圧縮性流体で十分であるため圧縮性性い問題がある。

流体の解析は軽視されてきたという背景が想像されるが、圧縮性回転流には回転浅水波方程式などの物理的にも重要な応用を持つ方程式があるため、その数学解析は重要な課題である。最近の私の研究により、圧縮性回転流は非圧
縮性回転流と異なる性質を持ち特徴的な挙動を示すことが明らかとなった。大きな違いは密度変化を考慮したことでコリオリ力が流体のエネルギーに影響を与える点である。特に、コリオリ力は０階の線形項であるため線形解の時間減衰率を悪くさせる効果を持つため、非線形解の時間大域的可解性の構成では困難点が多いが、技術的な工夫を施すことで解決することが出来た。

21

Institute of Mathematics for Industry◯ 産業数理統計研究部門

Kyushu University

九州大学マス・フォア・インダストリ研究所

L1正則化に基づくスパース多変量解析

廣瀨慧

学位：博士（機能数理学）（九州大学）
専門分野：スパース推定、L1正則化、多変量解析
近年、ビッグデータ解析が重要視されていますが、データの一般化であり、更に因子回転よりもスパースな解が得ら量が多くなった反面、データの冗長性も増してしまい、必要れるということを理論的に示すことができました。私はこののない情報を取り除いて有効な情報のみをうまく抽出する

関係性を見出しただけでなく、パラメータをある程度高速にことが必要とされています。それを実現する、極めて有効な推定できるアルゴリズムを考え、さらにソフトウェアパッケー

方法の一つが、L１正則化法をはじめとする、スパース推定

ジfanc（https://cran.r-project.org/web/packages/です。スパース推定とは、数万・数億にものぼる数のパラメー

fanc/index.html）を作りました。実際にこのパッケージをタが存在するときに、ほとんどを「ゼロ」と推定することがで使った論文もいくつかあります。

きる方法であり、非ゼロ要素に対応する変数のみが有効となります。この方法の良い所は、
たとえ数万オーダーの次元（２）‌データが大量に欠損する場合の因子分析モデルの最尤推定のデータであっても、わずか数分で計算が完了してしまうく
らいに高速なところです。計算効率が良い上に、統計的にもNTTと共同研究している際、データに大量の欠損がある

良い性質がたくさんあるため、L１正則化は多くの統計学者

場合の因子分析をしなければならないということがありましを魅了しているのではないかと思います。

た（図参照）。具体的には、
アンケートを行う際、すべてのアン
ケート項目に答えるのでなく、その中のごく一部を選び、その選んだ項目のみ回答するというデータの取り方をしまし私は、上記のL１正則化法を使った統計解析、
とくに、多変た。

量解析に興味を持っています。多変量解析とは、大量の変数があった時、似た変数をうまくまとめたりすることにより、変データが欠損する場合、パラメータはEMアルゴリズムに

数間の関係性を見出す方法です。この方法は、古くから現在

よって推定できますが、欠損数が多い場合、通常のEMアルまでずっと使われ続けている基本的な手法です。私は、多変ゴリズムだと計算時間がかかってしまいます。そこで私は、

量解析の中でも、因子分析と呼ばれる手法に興味を持って大量欠損時の因子分析モデルにおけるEMアルゴリズムをいます。因子分析は、
もともと心理学者が作った方法なので提案しました。このアルゴリズムは、従来のEMアルゴリズム

すが、近年は生命科学でも使われています。また、因子分析

よりも、数百倍、場合によっては数千倍もスピードが早い事を拡張したテンソル分解は、機械学習の分野で使われ始めがわかりました。

ています。因子モデルの研究は今後ますます発展していくと考えられます。以下、私の最近行った２つの研究を紹介します。
（１）因子分析のスパース推定
因子分析の面白い点は、モデルそのものに識別性がないというところです。このような解が一意に存在しないモデル
は、統計学者からすると「奇異なモデル」となります。実際にどのようにパラメータを推定するかというと、因子回転と呼ばれる、因子分析独自の最適化問題を解きます。このような問題設定自体、他のモデルにはない独特なものですが、それが50年以上も使われてきたスタンダーダードな方法です。
私は、この因子分析モデルにL１正則化を適用すると何が起こるかを考えてみました。すると、正則化法が、因子回転22◯ 産業数理統計研究部門

Finite mixture models for statistical inference

Hien Duy NGUYEN
学位：PhD（University of Queensland, Australia）
専門分野：Mathematical‌
Statistics, Statistical Computing, Statistical Learning, Bayesian
Statistics, Signal Processing, Stochastic Programming, Optimization Theory
My research in mixture model computation, estimation, and inference
has found widespread application in real-world scenarios. For example,
I have collaborated with neuroscientists and cell biologists to analyze
heterogeneous biological phenomena, worked with quantum physicists
to characterize switching behaviors of quantum circuitry, assisted
economists in characterizing subpopulations of experimental outcomes,
partnered with civil engineers to study regional differences in traffic
behavior, supported fisheries scientists in characterizing growth stages ofaquatic species, and collaborated with image scientists to segment andcharacterize imaging data, among other practical applications.

Many real-world datasets are heterogeneous and multipopulational
phenomena. In such contexts, it is insufficient to capture the overallvariation among the data using a single statistical model. Therefore,a cohesive approach to modeling the multiple subpopulations within
the superpopulation is necessary. In such scenarios, a useful approach
involves modeling each subpopulation and their contributions to thesuperpopulation through a weighted averaging construction, known as afinite mixture model. These models are highly flexible and interpretable,
enabling them to capture and provide inference for known heterogeneities
in the data while also identifying new heterogeneous phenomena thatwere previously concealed.
The class of finite mixture models is extensive, and choosing
between different mixture models can be challenging. In my work, I havestudied model selection procedures required to make mathematicallyprincipled choices among competing finite mixture models. I have madeprogress in two key directions to address this problem. Firstly, I employ
sequences of hypothesis tests to determine the number of components or
subpopulations required in each mixture model. This approach relies ona new hypothesis testing method called universal inference, which offersa straightforward and assumption-light mechanism for deciding whether
a model accurately represents the observed data. Using these universal
inference tests, I have developed a way to construct confidence intervals
for the number of underlying subpopulations in the data, providing
insight into the complexity of the overall superpopulation.

Figure 1: Traffic crash rate clustering ofdifferent regions in Victoria, Australia.

Secondly, by leveraging modern stochastic programming techniques
for optimizing random objects, I have developed new penalization
methods for selecting between different finite mixture models within
broader model selection and decision problems. My novel information
criterion, known as PanIC, offers a more assumption-light alternative
to existing methods like the Bayesian information criterion or Akaike
information criterion. PanIC provides a single-number summary forchoosing between competing models, guaranteed to asymptotically selectthe correct model as the dataset size increases.
Figure 2: Quantization of mandrill
photograph using different mixture models.

Beyond their utility for modeling heterogeneous processes, finite
mixtures and their regression variants, the mixture of experts (MoEs)
also serve as excellent functional approximations of probability density
functions (PDFs) and conditional PDFs that characterize statistical
relationships. My colleagues and I have contributed to understanding
the approximation theoretic properties of mixture models and MoEs
for various classes of PDFs. We have provided sufficient conditions for
ensuring that PDFs, conditional PDFs, or mean functions of conditional
PDFs can be effectively approximated using a sufficiently large numberof components in a finite mixture model construction. These results,
often referred to as universal approximation theorems, are valuable fordetermining whether a class of functions serves as an adequate basis formodeling an underlying mathematical phenomenon.

Figure 3: Mixture-based false discovery
rate control of p-values for a mouse brain
morphometry experiment.

23

Institute of Mathematics for Industry
Kyushu University

九州大学マス・フォア・インダストリ研究所

位相的データ解析と層の理論

池祐一

学位：博士（数理科学）（東京大学）
専門分野：位相的データ解析、超局所層理論
私は、位相的データ解析と層理論について研究していまします。すると、穴がrに関してどのくらい継続（persist）
しす。これらはそれぞれ実応用と純粋な数学理論として全く異ているかの情報を見ることで、大きい穴とゴミのような穴をなるように思えるかもしれませんが、最近では互いに影響し区別することができます。この手法はパーシステントホモロ

合いながら研究が進んでいます。

ジーと呼ばれています。
このようにデータからパーシステントホモロジーを用いて（１）位相的データ解析（パーシステントホモロジー）

取り出した「大まかな形」をあらわす特徴量は、機械学習の図１にあるような２つの点群を見たときに、われわれ人間

入力などに使うことができ、様々なタスクに役立てることがは穴の有無で２つを区別することができます。このようなできます。私は、ニューラルネットワークにパーシステントホ

データの「大まかな形」を取り出して、その情報をコンピューモロジーを使うことでネットワークの状態を調べたり、ネットタに扱わせようとする手法が位相的データ解析です。数学

ワークの学習にパーシステントホモロジーを用いたりするこにおいては、連続な空間の「大まかな形」はトポロジーというとに興味を持っています。

分野で研究されていて、穴の数はホモロジーという道具で調べられます。この枠組みを使って、点群の形を調べる最も（２）層理論とパーシステントホモロジー

単純なアイデアは、各点を中心とするある半径rの円盤たち

上のような研究の他に、層理論とパーシステントホモロの和集合を考えて、そのホモロジーを見ることです（図２）。

ジーのつながりについても研究しています。層は代数幾何やしかし、この方法では円盤の半径によって、ゴミのような小トポロジーにおいて役立つ数学的対象ですが、近年は層をさな穴しか見えなかったり（図２（a））、穴が全然見えなかっ使ってパーシステントホモロジーを理解しようという試みが

たり
（図２
（c））
します。

いろいろと現れています。
私は、
パーシステントホモロジーの間の距離を層理論的に解釈して、その理解を位相的データ解析や他の数学分野に役立てようとしています。最近は、層の間の距離とジグザグ
パーシステンス間の距離との関係を調べたり、層の間の距
離とシンプレクティック幾何学におけるエネルギーとを関係づけたりしました。
図１：２つの点群

図2：異なる半径の円盤の和集合

このやり方には、見える形が半径rの大きさに大きく依存し
図3：パーシステントホモロジー・層理論・シンプレクティック幾何学の関係てしまう・適切な半径rをデータから決めることは一般には難しいという問題があります。そこで、１つの半径rを設定することは諦めて、rを動かして点群の「形」、特に穴の数をあら

応用数学と純粋数学の狭間での研究によって、それらのわすホモロジーがどのように変わっていくかを考えることに垣根を越えた交流が起これば良いと考えています。

24◯ 産業数理統計研究部門

ベイズ推定による計測と数理モデリングの橋渡し德田悟

学位：博士（科学）（東京大学）
専門分野：ベイズ推定、モデリング、情報統計力学
古くは17世紀に発見されたケプラーの法則が象徴するといった凝縮系物理学のモデルの選択において有用性を示

ように、単純な数式を用いて計測データを表す数理モデリ

してきました。

ングは様々な物理現象に対する理解を深めてきました。
しかし、扱う現象が複雑になり、高度な計測技術を駆使する現代（2）観測ノイズの存在
データの計測には必ず観測ノイズを伴い、観測ノイズが科学では理解が困難な計測データが顕在化しています。私は計測データに根ざしたモデリング原理を確立し、あらゆる大きいほど、データから推定するパラメータの値も不確実現象を曖昧さなく理解する指針を打ち出す研究に挑んでいになります。我々はこのような誤差伝搬の性質がモデルのます。特にベイズ推定という統計的手法の数理を探求する

評価にも影響することに着目し、ノイズの強度と妥当なモデとともに、凝縮系物理学を始め、広く自然科学の研究者とのルを同時推定する方法論を開発しました。また、実証研究を

異分野協働による実証研究を推進しています。これまでの通じてその有用性を示しました。パラメータの値やモデル研究を通じて、私は理解が困難なデータの要因となり得るの推定はノイズの強度（データの質）とデータ量に依存しま

次の三課題に着目しています。

す。我々はベイズ推定と統計力学の数学的な対応（情報統計力学）に立脚し、理論解析を進めることで、ベイズ推定に対する計測データの量や質のスケーリング則を解明しまし（1）モデルの不定性モデルは現象の本質を表しますが、何が本質かは常に自た。

明ではありません。多くの場合、その判断は研究者の洞察に委ねられ、時に研究者間で見解の相違が生じています。例え（3）モデル不一致の存在ば図1の計測データが示す振動現象は摩擦が無視できれば理想と現実、
つまりモデルと計測データの間には必ずズレ

単振動、そうでなければ減衰振動を表す関数によってモデリがあります。まず、モデルは真実の近似表現です。また、真実

ングできますが、どちらが妥当かは場合によります。

とデータの間には観測ノイズや系統誤差が生じます。私は

我々は計測データに対するモデルの妥当性を確率としてこれらの内、
ランダムなノイズ以外をまとめてモデル不一致定量化するベイズ推定を用いて、こうした不定性を解消すると総称しています。モデル不一致の起源を帰属することは実証研究を進めています。実際、速度分布関数やバンド構造

容易でなく、それらを具体的な数式で記述することはさらに困難です。加えて、ベイズ推定に正当性を与える従来の漸近理論はモデル不一致がない理想的な状況を仮定するものです。我々は実証研究を通じ、モデル不一致を系統的に対処する方法論を開発するとともに、その正当性を裏付ける
新たなベイズ推定の漸近理論の構築を進めています。

図１：計測データを構成する数理モデル、
観測ノイズとモデル不一致の具体例

25

Institute of Mathematics for Industry
Kyushu University

九州大学マス・フォア・インダストリ研究所

統計的ダイバージェンスに基づく頑健なモデル選択規準

倉田澄人

学位：博士（理学）（大阪大学）
専門分野：数理統計学、モデル選択、ロバストネス
例えば突出した能力、例えば災害級の現象、例えば観測機いた規準の頑健性が分かってきました（図2）。あらゆる現象

器の故障、例えば人的なミス…等々、現実のデータには様々に対してモデルが作成出来る以上、良きモデル選択規準をな由来を持った、他の観測値から見て大きく外れた値を取検討することは、あらゆる分野に対して意義を持ちます。頑るデータが付き物です。これを「外れ値」と呼びます。この外健性をはじめ、選択に資する様々な性質を詳しく調べることれ値ですが、何を以て外れ値とするか、どこまでが外れ値でで、文理を問わず幅広い分野を支えられる研究を目指してなくどこからが外れ値であるかという明確な定義、線引きはいます。

与え難いものです。また、災害やミスを完全に零にするということが非常に難しいため、外れ値の発生を防ぐことは事実上不可能と言わざるを得ません。その為、外れ値が混ざっていてもその影響を小さく抑えられる「頑健（ロバスト）」な分
析手法が、分析において重要な意味を持つと考えられます。
私は、モデル選択問題を中心に、頑健な手法について研究を行っています。モデル選択規準は、自然現象や人間の行動を表現する数理モデルの候補の中から最も好ましいモデ
ルを択ぶ尺度です。私はこれまで特に、確率分布間の遠さを（図1）一つは黒い線で画かれた「真の分布」、
もう一つは赤測る「統計的ダイバージェンス」を用いて、現象や行動の根い線の外れ値を発生させる分布、これらの異なる確率分布底に在ると考えられる「真の分布」とモデルとの近さを測るを考え、観測値はこの二つの分布の混合から発生しているという観点から導出されるモデル選択規準を検討してきまと想定します。もし、黒の分布のみからデータが出てきていした。モデル選択問題を考えるにあたっても、頑健性は重要るときと、混合分布から出てきているときで分析結果が大な性質と言えるでしょう。
しかし、従来広く用いられている、

きく異なれば、その手法は外れ値に対して敏感（頑健ではな

Kullback-Leiblerダイバージェンスに基づいたモデル選択い）と考えることが出来ます。

規準は、データに外れ値が混ざっていると選択精度を落としてしまう傾向にあることが指摘されています。そこで本研究
では、統計的推測において頑健性に優れることが知られている、「BHHJダイバージェンス」と呼ばれるダイバージェンスに基づいて、従来のAICやBICを拡張したモデル選択規準を導出し、検討してきました。
外れ値は頻繁に、パラメータ推定やモデル選択へ悪影響（図2）一般化線形モデルの数値シミュレーションにおける、を及ぼします。モデル選択における頑健性を議論する為に

幾つかのダイバージェンス（Divergence A-D）に基づいたは、モデル選択規準の値の、外れ値によって生じる変動をモデル選択規準のaccuracy rate（正しいモデルを選択評価する必要があると考えられます。本分野ではしばしば、した率）です。三つのグラフは異なる外れ値の混合率（0%、

データを発生させている分布に外れ値が含まれている場合

5%、20%）に対応しており、それぞれの横軸はデータ数をといない場合との差を評価することで、外れ値の混合に対指します。Divergence AとBに基づく規準は外れ値がある

する敏感さを評価します。観測値の大部分は「真の分布」かと精度が大きく下がっており、混合による悪影響を敏感に受ら発生していると想定し、そして外れ値は「真の分布」とはけていることが分かります。一方でDivergence Dの規準

別の分布から出現しているものだと解釈します（図1）。私はは、全ての混合率に対してほぼ同じaccuracy rateになっ

これまで、多くのダイバージェンスに基づいたモデル選択規ており、強い頑健性が見て取れます。

準の頑健性を比較検討してきました。結果として、BHHJダ
イバージェンスを筆頭に、幾つかのダイバージェンスに基づ26◯ 産業数理統計研究部門

小区分別統計的推測法の理論発展及びその活用

廣瀬雅代

学位：博士（工学）（大阪大学）
専門分野：統計科学、小地域推定、混合効果モデル
データサイエンス全盛のこの時代に、Evidence Basedが多くの場面で活用される。
しかし、従来用いられる方法はPolicy Making（EBPM）の重要性が認識されつつある。

実用面からの問題を複数抱えている。私は、その経験的最

その中で小区分ごとの実態把握に基づくより丁寧な計画立

良線形不偏予測量とその予測誤差の問題を考慮に入れた、案への期待も大きくなっており、調査データによって小区分

統計的精度の損失がほとんどない、改善可能性のある手法

ごとの特性値を推定する際、区分ごとの情報のみならず、統の開発を目指している。また、リサンプリング法と同様の統計的モデルを介して他の区分からの情報を有効に活用する

計的精度を保ちながら、計算負荷を軽減できる、より実用的ことで、統計的精度向上が期待できる場合がある。

な経験的ベイズ信頼区間法の提案も行っており、さらなるそのため、
このような統計的推定手法を活用したエビデン

効率的な統計的手法開発とその理論保証に取り組んでいス資料作成・提供は、有用な計画立案の礎になり、社会的課る。

題解決への期待を高めることができると考えている
（図１参
照）。

（２）実データ分析への応用

統計的モデルを用いたモデルに基づくアプローチの研究開発した統計的手法を実データに適用し、質の良いエビ

は、海外の官庁統計分野でも盛んに行われ発展を遂げていデンス資料作成に取り組んでいる。これまでにも、過去に開る。官庁分野のみならず、他の様々な応用分野で重要視さ発した統計的推定手法を、わが国のある都市住民を対象とれている研究でもある。ここで、モデルに基づくアプローチ

した意識調査データに適用し、小区分別に防災意識レベルに関する私の研究内容を２つ簡単に紹介する。

の実態把握を図っている。さらに、国勢調査小地域集計との共通項目データを活用して、国勢調査結果との誤差を調べている。その結果、（国勢調査小地域集計を真とした場合,）（１）小区分別統計的推測法の理論
小区分ごとに特性値を推定する場合、統計的モデルを活

慣習的に用いられてきた推定手法と比べたときの適用手用したモデルに基づくアプローチでは、仮定モデル下の線法の一種の有用性を示すことに成功している。詳細は廣瀬ら

形不偏なクラスの中で平均二乗（予測）誤差を最小化する（2018、
日本統計学会誌）を参照されたい。

最良線形不偏予測量に対する経験的最良線形不偏予測量

図１ EBPMのための資料作成における統計的研究の立場の一例

27

Institute of Mathematics for Industry◯ 先進暗号数理デザイン室

Kyushu University

九州大学マス・フォア・インダストリ研究所

数学と暗号の交差点

縫田光司

学位：博士（数理科学）（東京大学）
専門分野：暗号数理、秘密計算、組合せ論的群論
インターネット経由での安全な買い物など、現代の便利な 



生活の基盤となっている情報技術のさらに基盤となっている技術の一つが暗号技術です。さらにその暗号技術の基盤として数学が役立っています。1970年代後半に考案され

たRSA暗号では初等整数論が、1980年代半ばに考案された楕円曲線暗号では名称の通り楕円曲線の性質が活用さ





れています。さらに、来たるべき量子コンピュータ時代に備える次世代暗号（耐量子計算機暗号）には、線型代数やグレブナー基底や高次元整数格子などまた別のさまざまな数学が現れます。
数学者の目線で暗号分野を眺めたとき、「より新しく多彩

用いられていなかった群論的な手法を用いてこの複雑な操な数学を実戦投入する」楽しさはもちろんのこと、それ以外

作を回避することを目指しています。この研究には、数学にに「素朴な数学的道具でも、使い方次第で大きく花開く」とおける私の専門である組合せ論的群論の知見も活用できる

いう楽しさや、「 安全性」など暗号分野の様々な概念に対しと見込んでおり、その意味でもやりがいのある研究です。

て「より直感に沿った、かつ理論的にもより取り扱いやすい

他にも秘密計算に関する最近の研究として、従来この分「良い」定義を探究する」
といった楽しさも感じられます。

野で安全と考えられてきた「直感的にもっともらしいある種の構成法」が実は安全でない場合があることを示す「病的な暗号分野で私が近年主に取り組んでいる研究テーマは反例」を発見しています（国際会議PKC 2021で発表）。こ「秘密計算」と呼ばれる技術です。秘密計算は、複数人がれもある意味で「いかにも数学者的」な研究でしょう。

データを持ち寄って行う情報処理において、「データを互い

これらの事例に限らず、暗号分野には（他の応用分野でもに見せない」のに「欲しい計算結果はちゃんと得られる」といそうでしょうが）本当に多種多様な数学が関わっていて、私う手品のような性質を実現する暗号技術です。秘密計算に自身の研究でも「まさかこんな数学が役に立つとは」と驚い

よく用いられる数学的・暗号的な道具のうち、「 完全準同型たことが何度もありました。そうした数学と暗号との豊かな暗号」という技術は、暗号化したままの状態で暗号文の中身

関連性についての情報発信もより一層進めていきたいとのデータに対する計算を行える特殊な暗号化技術です。私

思っています。

たちの論文（国際会議EUROCRYPT 2015で発表）
で構成した完全準同型暗号では「有限体上の関数はどれも多項式で表示できる」という性質が要となりました。この性質自体は代数学の初歩的な内容なのですが、それを「上手いタ47

イミングで使う」
ことで暗号学的に有意義な成果が得られた/=

という実例となっています。

"3

この完全準同型暗号については、既存の構成手法におい

&#>




化の妨げとなっています。私の最近の研究では、従来研究で06:5


;#'%







1-!
,2><>
*>
9+




28

&#
>

)$(

てどうしても外すことのできない複雑な操作が存在し、効率8.

◯ 先進暗号数理デザイン室

耐量子暗号の安全性解析

池松泰彦

学位：博士（数理学）（九州大学）
専門分野：耐量子暗号、多変数多項式暗号
現代の情報セキュリティを支えるRSA暗号・楕円曲線暗号

MPKCは、その他の耐量子暗号（格子暗号や同種写像暗は、素因数分解問題・離散対数問題などの計算困難数学問

号など）と比べても高速な処理性能を有し、さらに署名方式

題を安全性の根拠として構成されます。
しかし、大規模な量においては署名長が最も短くなることから、IoTで取り扱う

子コンピュータが将来的に開発されると、Shorの量子アル機器やICカードのような低資源デバイスでの実装に向いてゴリズムにより、素因数分解問題・離散対数問題などは多項いると考えられています。また最近では、署名方式の一つで式時間で解読でき、結果として、RSA暗号などは安全でなくあるRainbowから仮想通貨を作るという面白い試みも行なることが知られています。

われています。

量子コンピュータでも解読困難とされる暗号（耐量子暗私は、主にMPKCの安全性解析について研究を行なって号）は、上記のような背景から、現在世界中で活発に研究開います。それには、
グレブナー基底や代数幾何などによる数発が行われています。この研究には、素因数分解問題・離散学的議論が欠かせません。
しかし、グレブナー基底アルゴリ

対数問題とは異なる数学問題が扱われ、広範囲な数学理論

ズム、特にF４/F５アルゴリズムによる攻撃計算量の解析は、が必要となります。

理論的にはまだまだ未解明であり、実験による解析に基づいています。私が現在研究している暗号方式EFCにおいても、Hybrid攻撃を用いることで、安全性が想定されていたものより低下することを実験的に示しましたが（下図）、それらを理論的に説明することが未解決であり、今後の研究課題となっています。

私は、耐量子暗号の中でも、特に多変数多項式理論を用いた、多変数多項式暗号（MPKC）に興味を持ち、研究を行なっています。MPKCは、有限体上の多変数二次連立方程式の求解問題（MQ問題）の困難性を安全性の根拠として構成されます。このMQ問題はNP完全問題であることが証明されており、そのことから、量子コンピュータでも解読困難な暗号が構成できると期待されています。
またMPKC以外にも格子暗号や同種写像暗号についても研究を行なっており、これらを組み合わせた新しい耐量子暗号の開発にも関心があります。

29

Institute of Mathematics for Industry◯ リエゾン戦略部門

Kyushu University

九州大学マス・フォア・インダストリ研究所「変な振る舞い」＝「素直な振る舞い」

松江要

学位：博士（理学）（京都大学）
専門分野：力学系、数値解析、精度保証付き数値計算、特異摂動論、燃焼など
物体の動き、空気の流れ、物体や空間の温度など、多くの自然現象は微分

図２：無限への発散＝地平線への収束

方程式の解で記述されます。典型的に、解は徐々に（定数や時間周期的な振動する関数などで記述される）一定の状態に近づく、あるいは指数関数的に発散しますが、時折この枠組みに当てはまらない「特異」な動きをするものがは時間微分）
について、時刻t=0におけるuの値をあります。例えば u’= u２（’
正の値に取ると、解はある有限の値Tに対して、t→Tの時に無限大に発散し、それ以降方程式は「解けなく」なります。
これは典型的な物理現象では起こり得ない「変」な振る舞いで、解が「有限時間爆発」すると言います。有限時間爆発は熱源着火に付随する熱暴走、連鎖的化学反応による反応暴走、エネルギー集中に伴う極短時間における梁の急激・爆発的な振動など、特異な現象の数学的対象物としてしばしば確認されます。一方、上の方程式が示すように「系＝現象の支配法則」自体は特異性がない場合も多く、与えられた系に対して「爆発は起こるのか」、起こるとすれば(i)：平面の放物面への埋め込み。焦点と平面上の点Mを通る線分を引き、Mを放物面の交点P(M)（ の平面上への射影）と対応させる。Mが無限遠方に行く時、その方向を保ってP(M)は放物面の境界：地平線に近づく。(ii):解が地平線へ近づく様子。近づき方を記述する集合：「局所安定多様体」の記述をもって、発散する解、特に爆発解の記述が実現される。

「いつ、
どこで、
どのように」起こるのかは非自明な問いであり、このような振る舞いを制御するには現象そのものの理解を深める必要があります。他方、爆発現象の場合は「無限大」あるいは「無限遠」を直接相手にするため、数学的にも数値的にも捉えることが難しく、この現象そのものが長年の数学的な研究対象とされています（図１参照）。

このアプローチの利点の１つとして、数値計算・数値解析など、他の理論や技術との親和性が非常に高いことが挙げられます。
（１）：爆発の有無を含め

図１：解は爆発する？しない？

た解の振る舞いを網羅する、（２）：初期値の摂動で壊れてしまう爆発解を、数学的厳密性をもって数値計算する、（３）：「複雑」な爆発の様子も地平線の様子で正確にわかる、などです。これらは精度保証付き数値計算、特異点論などを組み合わせることで実現できます。
図3：爆発解を『視る』

ある2次元常微分方程式の解。初期値をわずかに変えて初期値問題を解き、ある時刻で止めた。この後、これらの解はどう振る舞うだろう？
正解は、(i)：時間を進めると(u,v)はともに0に収束する。(ii)：(u,v)はともに有限の時刻で発散する、すなわち「有限時間爆発」する。方程式を構成する右辺（ベクトル場）
は滑らかである。ここからこの描像が想像できるだろうか？
(i):上に示した系の解全てを放物面に埋め込んだ時の様子。赤い領域を通る解は全ての時間において有界で、原点、あるいは有界な点に収束する。それ以外の領域を通る解は爆発解に対応する。爆発の様子は、地平
線上の4点のどこかに収束することをもって記述される。(ii):別の系における「初期値の摂動で壊れる」爆発解。力学系理論と精度保証付き数値計算を組み合わせることで、振る舞い方だけでなく、各点における爆発時
刻なども正確に求められる。

ここで私は、「無限遠」など変・特異な振る舞いを誘発する因子の捉え方を変えることで、物体の動きや温度と同じように「素直」に爆発現象などの「変な振る舞い」を捉え、統一的記述を実現することを研究テーマの１つとしています。具体的な１つのアイデアは「空間全体を半球面やカットされた放物面の中に埋め込み、その境界として無限遠を表現する」というものです。複素関数論
におけるリーマン球面、位相空間論におけるコンパクト化に類するものですが、ここでは系の持つスケール性にも着目し、適切な（境界つき）曲面への全解の爆発だけでなく、時間に付随して起こる様々な変な振る舞い：「有限空間の埋め込みを構成し、その境界：「地平線」
として無限遠を表現します。こ

時間特異性」は捉え所が難しいものです。
しかし、上記のように見方を変えるの考え方自体は代数幾何学などで見られる特異点や無限遠点の捉え方に端ことで、実は別な角度から見た素直な振る舞いに過ぎないことが徐々に明らを発するものですが、全ての解のとり得る（定性的な）振る舞いを包括的に記かになっています。この考え方は、古くはマルチスケールダイナミクス、最近

述する力学系理論と組み合わせることで、「無限大に発散する解＝地平線上

では転換点（Tipping point）なども、力学系理論や幾何学の知識で「素直なの集合に収束する解」という対応を得られます。さらに、地平線への収束の仕解釈」が可能となっています。これらを包括的に記述する思想は、様々な点で方から、変換する前の解の爆発の様子を正確に記述することも可能になりま「変なもの」＝「素直なもの」という見方を得ることに繋がるでしょう。

す。これは力学系理論の幾何学的側面、代数幾何学由来の系の適切な変換、漸近解析を組み合わせて実現されるものです。

30◯ リエゾン戦略部門

ソーシャル数理と動的最適化

吉良知文

学位：博士（機能数理学）（九州大学）
専門分野：ソーシャル数理、動的最適化、確率最適化、マルコフ決定過程
私の専門分野は数理最適化です。特に、意思決定が繰り

実際、各自治体では試行錯誤に多くの人手と時間を要し返し必要な問題（多段意思決定過程）
や不確実性を含む問題ており、自治体によってはきょうだいが別々の保育所となる（確率モデル）に関心があります。私は、この種の問題を効

ケースが増えるなどの問題が生じていました。そこで、IMI

率良く解決する手法である動的計画法の理論とその応用のの富士通ソーシャル数理共同研究部門（2014年９⽉〜

研究に従事してきました。また、社会的課題に対して、数理技2017年８⽉）では、富士通研究所とともにこの課題に取り術（最適化やゲーム理論など）を用いて、公平で納得度の高

組み、公平な利用調整を実現する新しい方法（展開形ゲー

い制度や施策を設計する研究（ソーシャル数理）に従事してムの理論を用いる）
を提案しました。提案した方法は富士通

います。これまで社会的課題の現場と協働して技術開発を株式会社で製品化されており、実際に多くの自治体で利用おこなってきました。以下では、２つの事例を紹介します。

されています（2020年６⽉時点で35の自治体が導入済

保育所の利用調整（選考）では、申請者の優先順位（保み）。

育の必要性を点数として算出）だけでなく、きょうだいの同「物流クライシス」が叫ばれる昨今、物流業界は人手不足

一保育所への入所希望も考慮します。
“ Matching withが深刻化しており、より少ないトラックでより多くの荷物を

Couples”という問題であり、学術的にも難しいことが知ら運ぶことができる共同輸送の必要性が更に高まっていまれています。図１は申請者が不満に思うパターンを示していす。そこで、多数の輸送ルートが登録されたデータベースのます。この種の不満が生じない優れた選考結果を「安定マッ中から協力効果が高い共同輸送の組合せを瞬時に列挙してチング」といいますが、その存在すら保証されておらず、納提案する共同輸送マッチング技術を開発しました。例えば、得がいく調整は容易ではありません。

東京→金沢の輸送レーンをもつ企業からのマッチング依頼に対して、他の２企業と協力した三角輸送（東京→金沢、金不満 (a) 第１希望に自分より点数が低い子がいる

沢→大阪、大阪→東京など）などを提案します（図２）。その際、実車率（走行距離に占める貨物積載区間の割合）が高くなる組合せを瞬時に列挙します。単純な総当りでは多数の依頼を処理するサービスは困難です。距離の公理などを用いて実車率の上界を上手く計算することにより、正確さを損なうことなく探索範囲を絞り込んでいます。本技術は日本パ
レットレンタル株式会社が提供する共同輸送マッチングサー
ビスTranOpt（トランオプト）にコアエンジンとして搭載され不満 (b) 第３希望に兄妹そろって入園できたはずており、既に約180社の企業が利用しています（2023年
10⽉時点）。

図１：申請者が不満に思うパターン

図２ 三角輸送のマッチング依頼に対する処理

31

Institute of Mathematics for Industry
Kyushu University

九州大学マス・フォア・インダストリ研究所

数値シミュレーションによる現象の理解

田上大助

学位：博士（数理学）（九州大学）
専門分野：数値解析、計算力学
た。またこの成果を発展させ、ガラス溶融炉の最適設計を行う際の指標の1つである炉全体の熱収支計算に対して整合流束法に基づく離散化手法を提案し、熱収支計算に対する最適な誤差評価を導いた。我々の誤差評価は、境界積分で定義される物理量を直接計算するよりも整合流束法を用いて領域積分に置き換えて計算した方が高精度である、という事実の数理的正当化に対応している。このように数理的正当化した離散化手法に基づく数値シミュレーション
を用いて、製品品質の維持とエネルギー消費の抑制を両立したガ
ラス溶融炉の最適設計について検討を試みている。
その他に我々が取り組んでいる研究の一つとして、変圧器内部の渦電流に代表される磁場の数値シミュレーションがある(図3参照)。
図1.数値シミュレーションの概念図。斜体文字は非圧縮粘性流れの場合の対応する例。

複雑な形状を持った領域や変化の著しい現象の数値シミュレーショ
ンを精度良く行うためには、自由度数が107(あるいはそれ以上)の水の流れや熱の伝達など、自然界や産業界で見られる様々な現近似方程式を効率良く扱うことがしばしば要求される。我々は静磁象を、計算機を用いた数値シミュレーションによって理解することに場問題に対して、既知量である電流密度の補正を考慮した混合型

興味を持っている。数値シミュレーションでは、まず現象を物理法則定式化による数理モデルを導入した。この数理モデルに対して反復に基づいて微分方程式で記述する 数理モデル化 を行い、次に微型領域分割法に基づく近似方程式を提案し、さらに混合型定式化の分方程式を計算機で扱うことのできる近似方程式に置き換える 離際に現れるLagrange乗数の性質を利用した簡略化を行った。提案

散化 を行い、最後に近似方程式の解法を計算機に実装し目的とすした手法は、必要となる行列ベクトル積演算を部分領域ごとに独立る現象を再現する 数値計算 を行う(図1参照)。我々は、提案するな静磁場問題を解くことに帰着できるため、計算効率向上のために離散化 手法が本来の現象に対する 精度の良い 離散化になって用いられる並列計算に適しているという利点がある。以上のことかいるか、実用に耐えうる 効率の良い 離散化になっているか、を数理ら、提案した手法に現れるある反復計算の収束特性が改善され、計的に検討する研究に取り組んでいる。また様々な問題に対して、提算効率が向上した。これにより、従来は解くことのできなかった大自案する 離散化 手法に基づいた 数値実験 を実際に行い、自然現象

由度モデルの解析に成功した。このように効率化した離散化手法にの理解や工業製品の設計へ応用することにも取り組んでいる。

基づく数値シミュレーションを用いて、磁場の振る舞いをより精度

我々が取り組んでいる研究の一つに、ガラス溶融炉内部における良く把握することを通して、変圧器の最適設計についての検討を試

ガラス原料の流動や熱伝達の数値シミュレーションがある(図2参みている。

照)。ガラス溶融炉内部にある高温(約1,500℃)のガラス原料は非

さらに我々は、粘弾性流れ問題、移動境界流れ問題、光波の干渉・

圧縮粘性流体と見なせることから、溶融炉内部における現象は温度散乱問題などの数値シミュレーションにも関心を広げて研究を続け

依存係数を持つ熱対流方程式によって数理モデル化できる。我々ている。これにより、数理的正当化がなされた数値シミュレーションは時間方向に後退Euler法を、空間方向に混合型適合有限要素法

がより多くの現象に対して可能になると期待できる。

を用いて得られる近似方程式の解に対する最適な誤差評価を導い

図2.あるガラス溶融炉内部における熱対流問題の数値シミュレーションによって得られた、ある時刻におけるガラス原料の温度分布。上図は炉に電極を挿入した場合、下図は炉に電極を挿入しない場合である。電極を挿入し電流を流すことで発生するJoule熱が対流の駆動力となるため、温度分布に違いが生じていることが分かる。

図3.ある変圧器内部における時間調和な渦電流問題の数値シミュレーショ
ンによって得られた、渦電流の流線。2つの図は、変圧器に対して互いに異なる側にある視点から鳥瞰した場合である。変圧器の壁面内部に設置された磁気シールド(水色や黄緑色の流線で囲まれた青い長方形部分)を囲むように渦電流が生じていることが分かる。

32◯ オーストラリア分室

Modeling of Solid-to-Solid Phase-Transformations in Shape-Memory Alloys
Homogenization and Gamma-Convergence Problems for Nematic Elastomers

Pierluigi CESANA
学位：PhD (Applied Mathematics)(SISSA International School for Advanced Studies, Italy)
専門分野：Partial Differential Equations, Variational Problems
Fig. 1a) Self-similar triple-star patterns observed in Pb 3(VO4)2 and MgCd
resulting from the rotation of twin boundaries and matching of different
crystal variants via kinematic compatibility.

The main focus of my research work is on rigorous mathematicalmodeling and analysis of multiscale and multiphysics systems inmaterials science. My investigations explore ways information anddisorder emerge and evolve generating complexity and patterns insmart materials such as martensite, nematic elastomers and liquid
crystals. Understanding the microscale features and mechanisms of
multifunctional materials and predicting their interactions on the overallmacroscopic properties is of strategic importance in the design ofmaterials for engineering applications. Two specific lines from my pastand current research are summarized below.

b) Numerical realization of a microstructure as a 2D fragmentation
process.

1) ‌Solid-to-solid phase-transformations in shape-memory
alloysAustenite-to-Martensite phase transformation is observed in various
metals, ceramics and biological systems. It is the activation mechanismof the Shape-Memory effect. Despite the vast potential to the shapememory effect, practical implementations have been slow, and to-date,mostly limited to NiTi. It is strategically important to improve and stabilize
the shape-memory effect in known materials and develop new modelingstrategies. A problem I have considered is the analysis and modelingof disclinations (topological defects at the lattice level characterized byrotational mismatch). Disclinations as in Fig. 1-a are characterized bya self-similar triple-star pattern resulting in intense rotational stretches.
Mathematically, I have shown that such configurations arise as a solutionof differential inclusion problems with special rotational symmetry
and rigidity. By identifying the basic algebraic structure underlying thedifferential inclusion I have computed exact solutions both in linearized
and finite elasticity models shedding light on the mechanism that drives
formation of triple-stars. Moreover, I have investigated onset of criticality
and self-organization in the evolution of martensite via sequential
avalanching. Here the modeling strategy describes the nucleation ofmartensitic variants as a branching random walk process (see Fig. 1-b).
The question that I addressed is the behavior of certain features of theself-similar structure thus formed and the computation of power laws
for the length interfaces in a martensitic transformation. This project
involved collaboration with the groups of J. Ball and B. Hambly (Oxford,
stochastic modeling of martensite); E. Vives and A. Planes (Barcelona)
and T. Inamura (TiTech) on experiments on avalanches and disclinations
in martensite. Work is in progress on the investigation of the activation
mechanisms that drive avalanches in metals with the ultimate goalof mechanically characterizing the dynamics of solid-to-solid phasetransformations.

c) analytical construction of a rank-1 4-phase microstructure in NLCEs.

2) ‌Homogenization and Gamma-Convergence problems fornematic elastomers
Nematic Liquid Crystal Elastomers (NLCEs) are a class of soft ShapeMemory Alloys that combine the entropic elasticity of a network of crosslinked polymeric chains with the peculiar optical properties of nematic
liquid crystals. A thorough understanding of the manipulation of opticalbirefringence in thin-films of NLCEs by mechanical, electric and thermalmeans is a tremendous mathematical task which has strategical potential
applications in materials design and fundamental sciences. Focusing
on the strain-order coupling in NLCEs, I have investigated mechanismsthat rule the low energy states in mechanically and geometricallyconstrained systems such as artificial muscles, sensors and actuators.
The mathematical language required to tackle NLCEs problems is thatof calculus of variations, Gamma-convergence and relaxation. These
are sophisticated techniques at the intersection of the analysis of PDEs,
functional analysis and measure theory based on energy minimization
approach and which are particularly suitable for the study of singularly
perturbed variational problems. Collaborations are in progress with theexperimental Lab of K. Urayama (KyotoTech).

33

Institute of Mathematics for Industry
Kyushu University

九州大学マス・フォア・インダストリ研究所

Algebraic specification

Daniel GAINA
学位：PhD（Japan Advanced Institute of Science and Technology）
専門分野：Universal Logic, Formal Methods, Category Theory
flexibly react to client demands by allocating, for example, new server
units to meet higher rates of service requests. The model implemented
over the cloud is pay-per-usage, which means that the users will pay
only for using the services. Therefore, the cloud service providers haveto maintain a certain level of quality of service to keep up the reputation.
Generally speaking, reconfigurable systems are safety- and securitycritical systems with strong qualitative requirements, and consequently,
formal verification is needed.

Universal logic is a general study of logical structures with nocommitment to any particular logical system in the same way that
universal algebra is a general study of algebraic structures. The term
“universal” refers to the collection of global concepts that allow one tounify the treatment of the logical systems and avoid repetition of similar
results. One major approach to universal logic, in terms of both numberof research contributions and significance of the results, is institution
theory. This relies upon a category-based definition of the informal
notion of logical system, called institution, which includes both syntax
and semantics as well as the satisfaction relation between them. As
opposed to the bottom-up methodology of conventional logic tradition,
the institution theory approach is top-down: the concepts describe thefeatures that a logic may have and they are defined at the most appropriate
level of abstraction; the hypothesis are kept as general as possible andthey are introduced only on by-need basis. This has the advantage ofproving uniformly results for a multitude of logical systems. It leads toa deeper understanding of the logic ideas since the irrelevant details ofparticular logics are removed and the results are structurally obtained byclean causality. My research interests cover, roughly, institution theory
and its applications to computing science.

Table 2. Lifecycle of reconfigurable software(3) System development
I am currently maintaining the Constructor-based Theorem Prover(CITP), a proof management tool built on top of an algebraic specification
language for verifying safety properties of transition systems. The
methodology supported by the tool is not intended for formalizing
mathematics, but for the application to the development of softwaresystems. In order to achieve the targeted goal, the following importantresearch directions are pursued:(a) proposing more expressive logical systems to allow engineers tospecify easily and accurately the software systems,(b) develop decision procedures that can reason efficiently about thesemore sophisticated logics, and(c) improvements of the proof assistant interface to help the user
understand the current state of the proof and interact with the tool in amore natural way.
The interest is in the design of software systems as one can see in thetable below.

(1) Foundation of system specification and verification
There are many contributions of institution theory to computing
science, the most visible one is providing mathematical foundations
for the formal methods techniques, i.e. specification, development
and verification of systems. In algebraic specification, one of the mostimportant classes of formal methods, it is a standard and mandatory
practice to have an institution to underlie each language basic featureand construct; institution theory sets a standard style for developing an
algebraic specification language that initially requires to define a logicalsystem formalized as an institution and then develop all the language
constructs as mathematical entities in the framework provided by theunderlying institution.

Table 1. Algebraic specification language development(2) Reconfigurable software systems
The main direction of my research consists of developing logicalstructures supporting the efficient development of correct reconfigurable
software systems, i.e. systems with reconfigurable mechanisms
managing the dynamic evolution of their configurations in response toexternal stimuli or internal performance measures. A typical example
of reconfigurable system is given by the cloud-based applications that

Table 3. Deductive verification process

The system will be specified at the most appropriate level ofabstraction depending on the requirements for its behavior. The result of
the verification performed with the tool will determine if improvements ofthe design are required or not.

34

九州大学マス・フォア・インダストリ研究所〒819-0395 福岡市西区元岡744
TEL: 092-802-4402 FAX: 092-802-4405
URL http://www.imi.kyushu-u.ac.jp/
UpdateApril 1st, 2024