数理解析研究所講究録
第 849 巻 1993 年 85-9085

結び目不変量の物理学への一つの応用*
Tetsuo Deguchi
哲生

出口東京大学理学部物理学教室

東京都文京区本郷 7 丁目 3 番地 1 号
環状高分子とは長い高分子鎖が閉 じているもののことであるが、環状になったことによって統計力学的にどのような影響があるのであろうか。 環状高分子の具体例としては、 環状 DNA 分子が考えられるが、 その他にも様々な環状の高分子がある。 分子が幾っか集まて一つの単位を構成し、 それが大体$10^{5}$っ個以上鎖状にっながっているものを想定 している。
環状の DNA は生体内で リラ形 (スーパーコイリング)

ックスした形と縮み上がったようなの二っの形状をとることが知られている。

そしてこのスーパーコイリングの形状をうまく選ぶことによって、$D$*数理解析研究所短期共同研究 代数的組み合わせ論と低次元トポロジー、 1992 年 12 月 14 日–16 日

86NA の複製や蛋白質へ情報を伝達したりする生化学反応をコントロールしているらしい。 そして、 トポアイソメラーゼという酵素が D$N$

A 二重螺旋を切ったりっなげた りしてそのリンキングナ ンバーを変えて、 リラックス形とスーパー コイリング形の間の変換をやっているらしい。 ここでもしランダムにっなぎ変えが起こったらどうなるか、ということを調べ、 そしてこれを手始めにしてトポロジーの効果を組織的に調べて DNA の形状変化制御の機構を解析 しよう、 というのが DNA 分子に対して環状高分子の応用研究をすることの動機である。 このような問題意識は生物物理学と生化学の境界領域では196$0$

年代以来ずっと引き継がれて今日に至っている。 スーパーコイリングをうまく制御出来れば、 生命の起源を理解することも可能となるであろう。

分子数が小さい場合の通常の統計力学的取り扱いとは異なり、 高分子の統計力学は一般に独特な取り扱いをする必要がある。 例えば鎖が長いために長距離の相関があるとみることもできる。高分子ネ ットワークのからみあいから、 ゾルゲル相転移などの臨界現象が説明されている。 (現在この現象は紙おむつに応用されている。) 環状高分子を研究するには閉じることによるトポロジー的制限をうまく取り扱う必要があり、 ここが理論的には非常に興味深い。 以下で議論されるように不変量をうまく用いてこの効果が議論出来るようになれば、 これは多体問題と しても非自明な取り扱いになるであろう。

以上のような背景のもとで、 統計物理学あるいは物性物理学の基87
礎論の一っの研究テーマとして、 環状高分子の統計力学を考察する。

一番簡単な問題として、 閉じたランダムウオークを発生してその結び目型を判定するということを考える。 基本的には次のようなやり方である。 ランダムウオークは 3 次元の座標データの列として与えられる。 これを 2 次元平面に射影して、 何等かの不変量を計算し、 その結果から結び目型を当てる。 そしてランダムウオー

クを多数発生させて調べて、 統計的法則性を見出そう、 というわけである。 ラン
ダムウオークのステ ッ プの数はま ず $100$ 以上、 大体して試行回数は1000 ぐらい、 そ10000 回以上を目安としている。

この問題の要素は次の 3 点である。

1. 閉じたランダムウオークの発生
2. 結び目型の判定

3. 統計法則の発見
我々 は 2 に関して議論する。 1 に関してもガオッシアンウオーク、2 量化法、 モンテカルロ法、 分子動力学的方法など、 様々な方法が行われていることを注意しておきたい。 (3 に関しては現在研究が進行中。)

19 7$0$年代以来、 幾っかのグループがランダム結び目の数値実験を実行しているが、 結び目の判定に用いられる方法はいっもアレクサンダー多項式の特殊値$K$$\triangle_{K}(-1)$

であった。 ($\triangle_{K}(t)$とは結び目のアレクサンダー多項式のことで、 これは通常ある変数$t$のロー88ラン多項式である。)

また、 ジョーンズ多項式やその他の最近の新しい不変量を計算するということは全くなされていない。
これにはっぎの二っの理由がある。 (1) 数値的発散あるいはメモリーの制限。 (2) 計算時間が指数時間であるものは事実上計算不可能である。

アレクサンダー多項式の数値が使われていて多項式値が用いられ)ていないのは、 もし多項式値の計算(変数実行するとメモリーが足りなくなてしまうためである。 さてもしっ$t$をそのまま計算機で絶対値が 1 でない数値をこの変数に代入すると、 1$0$

$0$$0$個もの交差点がある場合 $t^{\pm 1000}$ というような項が計算の途中で出現してしまうので、 数値的に非常に大ぎな数と非常に小さな数の両方を取り扱う必要が生ずる。 そこで絶対値が 1 で非自明な数、 っまり $-1$ が代入されるのである。 アレクサンダー多項式は行列式を用いて計算出来るので$N$の3乗の計算時間で求められる。

$N$は交差点の数である。

っまり多項式時間であり、 実際に計算可能である。
ジョーンズ多項式についてはこれを計算するには指数時間がかかってしまうことが知られている。 例えば交差点の数が一っ増えると計算時間は 2 倍になるというわけである。 この困難のため、 ほとんど誰もジョーンズ多項式を数値計算的研究に応用 していなかったと考えられる。 (実はこの研究会の席上でJaeger 氏か ら、 ジ項式の一部分の多項式 (truncated Jones polynomial)間で計算する方法をョー

ンズ多を多項式時

Przytycka and Przytycki が示 しているという

89ことを教えて頂いた。 但し彼らの方式は多項式値の計算で、 数値的発散やメモリーの問題が解決できないようである。)

この問題に対して実は少し異なったやり方で解決法がある。 多項式を計算するの多項式を直接用いずに、 その係数を多項式時間で求めることができるのである。 (Deguchi

年4and TSuruSaki

1 9 9 2

月の日本物理学会での講演;preprint 1992, tO appear in Phys.

Lett. A. )トリ

ックは量子群の準古典展開を考え、 ジ ョー ンズ多項式 (あるいは量子群の表現から導かれる様々 な不変量)
開の係数 (整数値)$r$、次の係数は$N$の$r$の準古典展を計算するということである。 すると、 展開の乗の計算時間で求められるのである。 ここでは交差点の数そして交差点が 1$0$$N$個程度以下であれば実質的にはこの展開の係数をっかって結び目を分類することができる。 この方法
により、 数値的発散の困難あるいはメモリーの問題、 そして計算時間の問題の二っの問題が同時に解決できる。
こうして、 大きな結び目を計算機で判定するという問題は実質的には解決できた。 現在このアルゴリズムを用いてランダム結び目の統計的性質をいろいろ調べている状況である。 絡み且に対しても全く同様なアプローチができるので、 例えば環状高分子間の トポロジカル相互作用も調べることができる。 これはいわゆるゴム弾性の説明に用いられるようにエントロピーの力を求めていることに相当する。
さて、 応用研究の側面の他にも意外なおもしろい結び付きがあった。 じっはこの準古典展開の係数はバシリエフ不変量を与えている90のである。 というわけで、 バシリエフ不変量は多項式不変量の多項式計算可能な部分を抽出した不変量であると言える。 このことは、 昨年の夏にジョーンズ氏にプレプリントを送った際にす ぐに手紙の返事が来たが、 その中で実質的に指摘されたといってよい。 (バルナータンの論文を参考にするように、 との指示があった。) 多項式値と整数値の違いは、 結び目理論の不変量の一般理論(バシリエフ不変量は結び目不変量の統一理論といってもよいのかもしれない)においても、 実質的応用研究においても重要な違いとなっている。

最後に、 研究会に参加させて頂きましたことを、 河野先生に感謝いたします。