研究概要 |
今年度は、研究代表者が開発したシミュレーションコードを用い、新規購入したハイパフォーマンスコンピュータ上で、ヘキサティック・スメクティックB(HexB)液晶相の等温・等圧シミュレーションを系統的に行ってきた。本研究で用いているソフトマターのためのシンプレクティック解法は、ハミルトニアンの構造から熱の出入りが詳細に計算できるため、相対エントロピー・定圧比熱などの熱力学量の計算が、従来の方法に比して簡便に、しかも精度良く求められる。複数の状態のエントロピーを比較するには、同じ初期配置から出発して異なる状態へ持っていき、その際の熱の出入りを測ることにより行う。複数の熱力学量を求めることによって、動力学的な特徴に加え、定量的に各状態を区別することが可能であり、熱平衡状態に加え、今のところ、少なくとも10以上の HexB 準安定状態が見つかっている。 HexB 液晶相は、6回対称ボンド配向秩序を持つことを特徴とするが、実験によれば、6, 12, 18,...6n回対称ボンド配向秩序パラメーター間には指数的なスケーリングが成り立っている。今回系統的に行っているシミュレーションでも、6回, 12回, 18回対称ボンド配向秩序パラメータ C6, C12, C18 をすべて計算しており、実験と同様のスケーリング則が成り立つことがシミュレーションで初めて確認された。とくに今回のモデル分子では、スケーリングの指数が二次元から三次元へのクロスオーバーを見せており、二次元系の融解理論との比較が興味深い。また、ゆらぎの大きい相転移近傍の温度にするとあまり分子数の大きくない系では、複数の HexB 準安定状態を遷移する。このような事が観察できるのは、ゆらぎがどんなに大きくても安定に計算できるソフトマターのためのシンプレクティック解法のみである。これらのデータを使えば、位相空間のトポロジーも解明できる。
|