12月13日
13:05-13:45

ヤヌス粒子の自己組織化（九大院理）岩下靖孝、木村康之【はじめに】「パッチ粒子」と呼ばれるコロイド粒子は表面に物性の異なる領域(パッチ)を持ち、相互作用が異方的になる[1]。そのため分散系においてパッチ構造に応じた多様な構造を自己組織的に形成すると予想されており、実際に等方的なコロイド粒子にはない新規な結晶構造が形成された[2]。また生体高分子を含む多様な異方的粒子分散系の基礎的な性質、特に相挙動などを解明するためのモデル系としても研究されている[3]。
しかしパッチ粒子分散系の相挙動に関しては、理論・シミュレーションによる研究が盛んに行われている一方で、実験による研究例はほとんど無い。その主な原因は、実験で異方的な相互作用を制御することの困難さにある。そこで我々は室温付近で液-液相分離する溶液を分散媒に選び、その臨界点近傍で誘起される粒子間相互作用を利用した。粒子としては、最も単純なパッチ粒子である「ヤヌス粒子」(半球パッチを持つ粒子)を用いた。上記相互作用は粒子の表面物性に依存するため、ヤヌス粒子間の相互作用は異方的となり、かつ相互作用は臨界点に近づくほど強くなるため、温度により可逆的かつ連続的に制御できる。
【結果と考察】分散媒として水 - 2,6-Lutidine 臨界溶液を用いた。ヤヌス粒子は、粒径 1.5mm のシリカ粒子の半球面を金薄膜で被覆し、金面をイオン性チオールで保護したものである。溶媒中に沈殿したヤヌス粒子の 2 次元分散系を顕微鏡観察し、臨界点に近づくにつれ分散状態がどのように変化するか調べた。その結果、始めは分散していた粒子が 4 粒子からなる tetramer などの
小さなクラスターを形成し (図(a))、さらに臨界点に近づくと tetramer が凝集し一次元的なクラスターを形成する「階層構造」が現れることが分かった (図(b))。これらの結果をシミュレーションとも比較し、形成機構を解明した。
また六方格子におけるパターン形成に関しても興味深い方向秩序を発見した(図(c))。以上の研究により、ヤヌス粒子、あるいはパッチ粒子系特有の構造形成機構の一端を解明した。

a)

b)

c)

10mm

10mm

20mm
図(a)(b) ヤヌス粒子の凝集構造。黒い方が金面。(c) 六方格子中のストライプパターン。
【参考文献】[1] F. Sciortino and E. Zaccarelli, Curr. Opin. Solid State Mat. Sci., 15, 246 (2011); E. Bianchi, R.
Blaak, and C. N. Likos, Phys. Chem. Chem. Phys., 13, 6397 (2011)[2] Q. Chen, S. C. Bae, and S. Granick, Nature, 469, 381 (2011)[3] S. C. Weber and C. P. Brangwynne, Cell, 149, 1188 (2012); P. Li et al., Nature, 483, 336 (2012)

Copyright (C) 2013 ソフトマター研究会 All Rights Reserved.

12月13日
13:45-14:00

球状錯体超分子 M12L24 の自発形成シミュレーション(1 産総研,2 東北大,3 東大）○米谷慎 1,山口智彦 1,佐藤宗太 2,藤田誠 3【はじめに】
近年、金属イオンと単純なリガンド分子から、100 成分に迫る球状錯体超分子の自発形成が報告され、生体系における高度な自己組織化の理解の足がかりとして期待されている。本
研究では、分子シミュレーションを用いて、配位子・金属錯体系における自発形成過程の解明を目指しており、最も単純な球状錯体として知られる、Pd(II)と三座リガンドから極性溶
媒中で自己集合形成される M6L8 錯体系に適用した結果を報告している 1)。
今回は、より構成要素数の大きな、二座リガンド分子による M12L24 錯体系 2)に適用した結果を報告する。
【結果と考察】リガンド分子モデルは、AMBER 力場による united-atom モデルを用い、電荷分布は
B3LYP/6-31G(d)レベルによる RESP 電荷を用いた。Pd(II) とリガンド分子間の配位結合はCationic Dummy Atom Method 3)を用いてモデル化した。この方法では、配位結合を、その立体対称性（Pd(II)の場合は平面四配位）に対応して配置した dummy atom 上に分配した点電荷によるクーロン相互作用でモデル化する。配位結合による錯体自発形成プロセスを加速するための手法として、Langevin Dynamics (LD)計算、generalized reaction field、effective
ポテンシャルの 3 つを組み合わせた粗視化溶媒モデルを用いた。generalized reaction field
における遠距離比誘電率は DMSO を想定して 47.0 とした。
effective ポテンシャルとしては、
Lennard-Jones(LJ) 型 を truncated and shifted し て soft repulsive と し た 、Weeks-Chandler-Andersen(WCA)型ポテンシャルを用いた。LD 計算は、上記の WCA ポテンシ
ャルで Lennard-Jones ルーチンを書き変えた汎用プログラム GROMACS-4.0.7 を用いて行った。
上記のモデルを用い、Pd(II) を 48 分子、二座リガンド分子を 96 分子(M12L24 4 超分子分)
ランダムに配置した構造を初期構造として LD 計算を設定温度 343K (70℃) で行ったところ、
下図の様に球状錯体 M12L24 の自発形成が確認できた。形成過程を Pd(II)の配位数の変化から解析すると、実験的に示唆されている 4)、自己集合、構造進化(下図)、構造捕捉の、３つのステージから成ることが示された。

【参考文献】(1) M. Yoneya, T. Yamaguchi, S. Sato, M. Fujita, J. Am. Chem. Soc., 134, 14401 (2012).
(2) M. Tominaga, K. Suzuki, M. Kawano, T. Kusukawa, T. Ozeki, S. Sakamoto, K.
Yamaguchi, M. Fujita, Angew. Chem. Int. Ed., 43, 5739 (2004).
(3) Y. Pang, J. of Mol. Modeling, 3, 196 (1999).
(4) S. Sato, Y. Ishido, M. Fujita, J. Am. Chem. Soc., 131, 6064 (2009).

Copyright (C) 2013 ソフトマター研究会 All Rights Reserved.

12月13日
14:00-14:15

マイクロリンクルの溝に閉じ込められた液晶における欠陥構造（産総研ナノシステム）大園拓哉、福田順一【はじめに】
本研究で用いるネマチック液晶の内部配向構造は特に、フランク弾性エネルギーを最小化することが知られている。ここで配向分布は、液晶が有しているバルクの弾性定数と、境界条件（幾何構造＋界面アンカリング）に強く影響される。本研究では、自己組織化マイクロリンクル[1]という微細凹凸構造を用いて、これまでに無い特殊な境界条件を実現し、その結果発現する興味深い内部配向構造とそれに伴った欠陥構造について調査する。
【結果と考察】本研究で用いるネマチック液晶の内部配向構造は特に、フランク弾性エネルギーを最小化することが知られている。ここで配向分布は、液晶が有しているバルクの弾性定数と、境界条件（幾何構造＋界面アンカリング）に強く影響される。本研究では、自己組織化マイクロリンクル[1]という微細凹凸構造を用いて、これまでに無い特殊な境界条件を実現し、その結果発現する興味深い内部配向構
造とそれに伴った欠陥構造について調査する。

溝に閉じ込められた液晶
マイクロリンクル

ネマチックの場合、+1/2 の欠陥線がジグザク状に周期的に折れ曲がる[Figure 1]。さらにこの各直線要素では、下部の容易軸方向からの配向方“T”の釘マークは各面に向の僅かな回転が確認され、 投影したその付近の液晶の配向方向を示し、釘これはツイスト変形に起因す の頭が各面から飛び出す方向を示す周期的液晶配向の偏光線欠陥る。そのツイスト変形では右
顕微鏡像巻きと左巻きが交互に出現し Figure 1. ネマチック液晶に誘起されるジグザグ欠陥ている。補足的実験と理論的 の概略考察により、この構造形成には液晶弾性異方性が重要であることが説明できる[2]。ス
メクチック A の場合、フォーカルコニックドメインが溝の左右に互い違いに並ぶ複雑な構造が確認された[3]。
これらのように、マイクロリンクルという特殊な境界条件を用いることで、液晶配向における自発的な周期構造の発生する例を発見した。これらの欠陥構造は、コロイド配列などへの応用が考えられ、その動的性質[4]や制御[5]を含め、現在検討中である。
【参考文献】[1] Ohzono, et al. Soft Matter, 5, 4585, 2009. [2] Ohzono, Fukuda, Nature Commun. 3, 701,
2012. [3] Ohzono, Takenaka, Fukuda, Soft Matter, 8, 2395, 2012.[4] Ohzono, , Fukuda,
Soft Matter, 8, 11552, 2012. [5] Ohzono, Fukuda, Suzuki, Yamaguchi, Phys. Rev. E, 86,
030701(R), 2012.

Copyright (C) 2013 ソフトマター研究会 All Rights Reserved.

12月13日
14:15-14:30

レーザー照射が引き起こす液滴運動のモード分岐（京大院理 A、同志社大生命医 B）

高畠

芙弥 A、吉川研一 B、市川

正敏 A【はじめに】
化学的あるいは熱的な非平衡状態において、マランゴニ効果と呼ばれる界面の不安定化現象に起因して物体の自己推進運動が引き起こされることが知られている。これらの自己推進
運動は、等方的な条件ではアクティブブラウン運動と呼ばれるような不規則な運動を示すが、外部境界条件や運動物体の形状に異方性を導入すると、その非対称性に応じた規則運動を取り出すことができる[1, 2]。さらに、最近の研究により、与えられた非対称性だけでなく自発的な対称性の破れにも伴って運動モード間の分岐が起こることが明らかになりつつある[3]。
本研究では液滴の自己推進運動において、上記のような自発的な対称性の破れが生じる実験系を見出し、そのメカニズムを調べた。
【結果と考察】
液体基板上に浮かべた液滴にレーザー光を照射して局所的に加熱することで、液滴周囲に温度勾配を構築することができる。この時、熱マランゴニ対流によってレーザー焦点を中心とする液滴の振動運動が誘起された。さらに、レーザー出力を変化させると、直線的な往復運動から回転運動に転化する様子が観測された。この運動モード分岐は、温度勾配が十分に大きくなり液滴に働く駆動力が大きくなった際に、運動の左右対称性が自発的に破れることで生じることが分かった。また、マランゴニ対流に起因する駆動力の時間遅れを考慮に入れた力学モデルにより、実験結果を定性的に再現することができたので、報告する。

図：(a) 液滴重心の軌跡。油滴の運動は直線的な往復運動から回転運動へと転化する。油滴の体積は 20μl。(b) レーザー出力に対する並進運動エネルギーの変化。ある領域では角度成分が動径成分よりも大きくなり、回転運動が支配的になっている。

【参考文献】[1] Y. Sumino et al., Phys. Rev. Lett. 94, 068301 (2005).
[2] S. Nakata et al., Langmuir 13, 4454-4458 (1997).
[3] F. Takabatake et al., J. Chem. Phys. 134, 114704 (2011).

Copyright (C) 2013 ソフトマター研究会 All Rights Reserved.

12月13日
14:45-15:00

Force acting on a moving colloidal particle in a polymer solution
Takashi Taniguchi, Yu-ichi Arai, Remco Tuinier; and Tai-Hsi Fan
Kyoto Univ., DSM ChemTech, Utrecht Univ., and Univ. of Connecticut
Introduction: Colloid and polymeric systems has been used in various industrial products, e.g., paints,
cosmetics, displays and so forth. In a colloid and polymer dispersion system, it has been well-known
that a depletion zone of polymer segment in the vicinity of the surface of a colloidal particle is formed,
which yields an attractive interaction among colloidal particles. So far, the equilibrium properties suchas the strength of the attractive interaction as a function of distance and phase diagrams of
colloid-polymer mixtures has been intensively studied. On the other hand, there are a limited numberof studies on the dynamic properties of such systems. Therefore, here we focus on the flow behavior
and the change of polymer concentration profile around a moving colloidal particle with a constantvelocity. The aim of the present work is to clarify the effect of polymer concentration distortion by theflow on the frictional force acting on the particle.
Results and Discussions: When a single colloidal particle is moving with a constant velocity in adilute polymer solution, the flow can change the depletion layer and therefore the frictional force
acting on the particle and the flow patterns are tightly related each other. We investigated the changeof polymer depletion distribution around a particle moving with a constant velocity to a quiescentpolymer solution at bulk and the resultant frictional force acting on it. To understand the effect of thespatially dependent viscosity on the frictional force, firstly we investigate the constant viscosity case. I
found that the spherically formed depletion in quiescent state is distorted by the effect ofhydrodynamic effect. In front of the moving particle, polymer concentration becomes higher than thatat bulk. In the rear of it, on the other hand, the depletion is enhanced and becomes wider by the flow
effect. The frictional force along the particle moving direction is contributed from two forces, viscous
force Fv and osmotic force FS. One of the control parameters of the present system is the Peclet
number Pe.As Pe becomes larger, Fv dominates, while the osmotic one becomes small. At a very highPe, FS still gives a finite contribution, but Fv becomes very large (҃Pe). We also investigated the casethat the viscosity of polymer solution depends on the local concentration of polymer and therefore
varies spatially. We found that the flow behavior around the colloidal particle is quite different fromthe case where the viscosity is spatially constant. To understand the effect of spatially dependent
viscosity on the frictional force, we compare our numerical results with the results obtained by twolayer model [1] proposed by Fan et al. As a result, it is found that Fv under taking into account the
distortion of concentration profile around the particle is larger than the two layer model and the casewithout any change of volume fraction profile from the equilibrium state. It is found that the origin of
the larger frictional force in the case of the distorted concentration profile by flow comes from a larger
shear deformation in the vicinity of the colloidal surface (i.e., the apparent slippage effect [2]).
[1] T. H. Fan et al., Phys. Rev. E, 75, 011803 (2003).
[2] Remco Tuinier and Takashi Taniguchi, J. Phys: Cond. Matt., 17, L9 (2005).

Copyright (C) 2013 ソフトマター研究会 All Rights Reserved.

12月13日
15:00-15:15

高分子潤滑のマルチスケールシミュレーション（兵庫県大・シミュレーション）安田（京大・工）山本修悟
量一【はじめに】
複雑液体の非線形レオロジー特性を計算機シミュレーションによって調べる場合には，通常，非平衡分子動力学法（Non-Equilibrium Molecular Dynamics, NEMD）が用いられる．この方法では，流体のバルクの状態を模擬した微小な MD セル（境界壁は周期境界条件が用いられる）において一様一定な速度勾配𝛾̇ と温度分布 𝑇 をもった理想的な環境における分子集団の運動状態を再現し，そこでの応力 𝜎𝑥𝑦 を計算することで応力と流動・温度の関係
𝜎𝑥𝑦 = 𝜎𝑥𝑦 (𝛾̇ , 𝑇 )を明らかにする．物質の本質的な物性を明らかにする上では有効な手法であり，様々な速度勾配，温度，密度などを用いた多数のパラメータに対して NEMD を実行する
ことで対象とする物質の構成方程式の基礎となるデータベースを作成することも可能である．
しかし，複雑液体を材料流体とする多くのシステムでは，その流動状態や温度分布は境界条
件とも結合し局所的に大きく変化するため，単純なシステムにおいてもシステム全体の挙動が材料固有の物性によって予測されるそれとは大きくかけ離れることがしばしば起こる．このような問題を解決するためには，境界条件も含めたシステム全体での熱流動を適切な分子モデルを基に解析することができるシミュレーション手法の開発が必要とされている．
一つの解決策としては，システム全体を分子レベルで構成してその振舞いをフル MD シミュレーションによって解析する方法が考えられる．この方法は原理的には可能であるが，MD
シミュレーションは一般に計算量が膨大であるため，実際の製品・機器の時空間スケールでのフル MD シミュレーションは現実的には不可能である．そこで期待されるのが，マクロな挙動を計算する計算流体力学法（Computational Fluid Dynamics, CFD）と物質固有の物性を求めることができるミクロな MD シミュレーションをうまく相補的に接続することで適切な分子モデルを基にマクロな熱流動までを計算することができるマルチスケールシミュレーショ
ン法の開発である．
本研究では，著者らがこれまでに開発した運動量輸送に対するマルチスケール法を拡張し，巨視的なエネルギー輸送をカップリングした高分子液体の熱流動を計算することができるマ
ルチスケール法を開発する．そして，その手法を用いて，高分子液体の熱流体潤滑問題を解析し，レオロジー特性と内部構造と粘性発熱による温度上昇の関係を詳細に調べる．
【結果】
これまでのところ以下の 3 点が主な結果として明らかになった．
(1) 高分子液体の発熱は，平板の速度が小さい場合にはあまり平板の速度には影響を受けない
が，平板の速度がある閾値をこえると急激に温度上昇が起こる．
(2) また，温度上昇と同時に粘性の低下（シア・シニング）や高分子鎖の形態変化や運動状態の変化も起こる．
(3) 平板速度が大きい場合，高分子液体の温度，応力，高分子鎖の形態変化は複雑な挙動を示
すにもかかわらず，そのような非線形領域においても線形の応力光学則が良く満たされる．
【参考文献】(1) S. Yasuda and R. Yamamoto, Phys. Fluids 20, 113101 (2008).
(2) S. Yasuda and R. Yamamoto, PRE 81, 036308 (2010).
(3) T. Murashima, S. Yasuda, T. Taniguchi, and R. Yamamoto, JPSJ 82, 012001 (2013).

Copyright (C) 2013 ソフトマター研究会 All Rights Reserved.

12月13日
15:15-15:30

シア下での水和水粘弾性応答の分子層分解計測（関西医大物理学教室）影島賢巳【はじめに】
親水性表面での水和水については、その生命科学上などの重要性とも相まって、理論・実験の両面から精力的に研究が行われている。局所的な物性であるゆえに実験的研究には困難がつきまとったが、動的計測に適した原子間力顕微鏡(AFM)で新たな可能性が開かれつつある。しかし、却って解釈困難な実験報告も相次ぎ、事態はむしろ混迷の度を深めている。注
目すべき報告として、水和層のガラス状態の可能性を指摘した理論研究[1]や、AFM を用いた
shear thinning の実験的計測[2]などがある。ただし、文献 2 では、分子層ごとに分解した形での計測はなされていない。そこで、シア応力を加えた水和水の粘弾性を、AFM で層分解して計測し、物性解明の糸口を見出すことを目指す。
【結果と考察】純水中のマイカ基板を置き、図 1 のように、磁石つきの AFM 力センサーに一定振幅のトルクを与えて、SiN 製の探針を基板に平行に 30kHz で振動させ、水和水にシアを加えながら基板に接近させた。探針の振動振幅は約 0.1nm 程度と見積もられる。力センサーの振動をロックイン検出して、図２のように層ごとに分解された形で粘弾性が計測された。探針―基板間のギャップ中にあると推定される水層の数をグラフ中に記入してある。1 層だけの水を挟んだギャップの間隔が通常の層間隔より狭いが、2 層目が加わった際のギャップ増加はむし
ろ大きいことなどが読み取れるが、計測ノイズなどの影響もあるのでまだ検証が必要である。
粘弾性データの解釈についても当日触れる予定である。

図１

測定の模式図図２

マイカ上の水和水で計測したシア粘弾性の距離依存性【参考文献】(1) P. Gallo, M. Rovere, and E. Spohr, J. Chem. Phys. 113 (2000) 11324.
(2) T.-D. Li and E. Riedo, Phys. Rev. Lett. 100 (2008) 106102.

Copyright (C) 2013 ソフトマター研究会 All Rights Reserved.

12月13日
15:30-15:45

相溶高分子ブレンドにおける絡み合い：PPA による解析（ 山形大院理工）滝本淳一、Sathish K. Sukumaran【はじめに】
	  高分子液体のダイナミクス・レオロジーが分子鎖間の絡み合いに支配されていることは広く認識されており、管模型とその拡張は多くの成功を収めてきた。しかし、絡み合いや管に関してはその存在が最初から仮定されており、そもそも絡み合いとはどういうものかについては、実は良く理解されていなかった。これに対し近年 PPA（Primitive Path Analysis）と呼ばれる手法が開発され 1）、分子動力学シミュレーションによって得られた分子配置から分子鎖間の絡み合いを抽出・可視化し、定量的に調べることが可能になった。その結果、たとえば剛直な（特性比の大きい）高分子ほど絡み合い易く、管が細いことなども明らかになっている。そこで本研究では、剛直な高分子と柔軟な高分子の相溶ブレンドのシミュレーショ
ンによって得られた配置に PPA を適用して、ブレンド中での分子鎖の絡み合いを議論する。
【結果と考察】	  シミュレーションには、よく知られた Kremer–Grest のビース・スプリングモデルに角度
ポテンシャル U(θ ) = k(cosθ − cosθ 0 )2 を加えたものを用いる。ここでθは結合の外角で、最安定な外角 θ 0 =90˚（柔軟鎖:A）および 0˚（剛直鎖:B）の 2 種類を用い、k は 2kBT に固定する。
重合度は A,B 共に N=200 とし、剛直鎖の分率 φB を変化させた 5 種類のブレンド（鎖数合計50 本）を用い、十分熱平衡化して得られた分子配置から PPA により絡み合い点間重合度 Ne
管直径と a を求める（PPA には M. Kröger による Z-code2）を用いた）。
	  まず、ブレンド中での各成分鎖の末端間距離の 2 乗平均を調べると、それぞれの単一成分
系に於ける値とほとんど変化していないことがわかった。つまり、排除体積効果の遮蔽はブレンド中でも単一成分系中と同じである。次に各成分の Ne と a を求めると、図に示すように、
ブレンド中では２つの成分の管直径 a がほぼ等

9.5

しくなっていることがわかった。また、図中の LN

9Colby-Rubinstein の単一成分系に対する仮説をブ

8.5

レンドに拡張して a を予想した結果であり（両者はほとんど同じ）、シミュレーション結果と良く一致する。一方 CW は Chen-Watanabe らによる実験から提案されたブレンド則 3)である。我々の結tube diameter a

お よ び CR は 、 そ れ ぞ れ Lin-Noolandi お よ び

7.5

a(CR)
a(LN)

7

果との不一致の原因は検討中であるが、実験では

6.5

2 成分のモノマーの大きさが大きく異なる（シミ6ュレーションでは同一）ことなどが考えられる。

a(CW)

8

0

【参考文献】(1) R. Everaers, S. K. Sukumaran, et al., Science, 303, 823 (2004)(2) M. Köger, Comput. Phys. Comm., 168, 209 (2005)(3) H. Watanabe et al. Macromolecules, 44, 1570 (2011)

0.2

0.4

φ

0.6
B

0.8

1

12月13日
15:45-16:00

多様な高分子の混み合い
多様な高分子の混み合い効果が細胞内粘性へ及ぼす
混み合い効果が細胞内粘性へ及ぼす影響
効果が細胞内粘性へ及ぼす影響
九大院理 A, 東北大院工 B
西澤賢治 A, 藤原慶 B, 柳澤実穂 A, 水野大介 A【はじめに】
細胞内の力学的性質は、従来アクチン等の細胞骨格によって決定されると考えられてきた。
一方で細胞内は、様々なタンパク質等からなる高濃度の高分子溶液で満たされており、近年
これら高分子の混み合いも細胞内の力学特性へ強く影響することが示唆されてきた 。そこで
本研究では、単成分の高分子溶液と多成分高分子からなる３種類の細胞抽出液を用い、高分子濃度を変化させながら光褪色後蛍光回復法法（FRAP）とマイクロレオロジー計測を行い、高分子濃度と粘性率の相関を評価した。さらに、生細胞内において同様の測定を行い、上記のモデル細胞質系と比較することで高分子の混み合い効果が細胞内粘性へ及ぼす影響へ迫る。
【結果と考察】単成分高分子溶液として球状タンパク質である BSA 溶液を、多成分高分子溶液として３種類の細胞抽出液（原核細胞の抽出液、真核細胞であるカエル卵と HeLa 細胞の抽出液）を用い、高分子の質量濃度（c）を変化させながら粘性測定を行った。その結果を図１に示す。全ての試料において、高分子濃度増加に伴う粘性率の上昇が観察された。また、その濃度依存性はガラス転移点近傍の振る舞いを表す経験的な式：Vogel-Fulcher-Tamman 式に従った（１）。
そのフィッティングから粘性が発散する高分子濃度（c’）を見積もると、単一系（BSA）では
乱雑最密充填時の濃度と一致した。多成分系で粘性率が発散する濃度は細胞種によらず一定となり、またその濃度は単成分系の約半分であった。さらに生細胞内の粘性率は、多成分系から推定される粘性率よりも低い値を示した。以上から、多成分高分子溶液では高分子の形
状や大きさ、表面電荷などの不均一性によって粘性が高くなること、細胞内では例えば力生成などによって粘性が低く抑えられていること、が示されたので報告したい。

図１：粘性率 η の高分子質量濃度 c への依存性【参考文献】(1) J. Mattson, et al., Nature, 462: 83 (2009).

Copyright (C) 2013 ソフトマター研究会 All Rights Reserved.

12月13日
16:15-16:30

ベシクルの自己生産の三次元画像解析 	 	
 （ 東北大学理学研究科）神保	
 岳大、佐久間	
 由香、今井	 正幸	
 （山口大学物理・情報科学科）浦上	 直人	
 （Jozef	
 Stefan	 Inst.,	
 Univ.	
 of	
 Ljubljana）Primoz	
 Ziherl	 	 【 は じ め に 】 	
 生体膜のモデルであるベシクルは、顕微鏡下での実験や膜弾性理論によってその物理的性質を広く研究されてきている。特に、ベシクルが余剰面積を得ることによってみせる自己生産
過程は、生命にとって必須の機能を単純な分子の集合体で実現できるという点で注目を浴びてきた[1]。また、脂質二重膜の三次元的な形状変化を実験から定量的に測定することはこれまで困難であったが、最近になって高速共焦点レーザー顕微鏡を用いて得られたベシクルの三次元画像を解析することで、そのダイナミクスを定量化することが可能になった[2]。	
 ベシクルの自己生産過程がどのような形状変化を経るのかを知ることは、自己生産の理解のために必要不可欠である。今回我々は三次元画像解析の手法を DPPC と DLPE の二成分系からなる自己生産ベシクル[1]に用いることで、自己生産過程の定量化に成功し、それを理論的に説明することを試みた。	
 	
 【 結 果 と 考 察 】 	
 DPPC と DLPE からなる二成分ベシクルは温度コントロールによる相転移で余剰面積を与えられ、自己生産過程が実現される[1]。この過程を共焦点レーザー顕微鏡によって三次元解析することで、reduced	
 volume	
 (v)	
 と	
 reduced	 area	 difference	 (∆a)	
 がどのような軌跡を描いて変化するのかを明らかにすることが出来た。実験によって得られた	
 (v,	
 ∆a)	
 の軌跡は Area	
 Difference	 Elasticity(ADE)モデルから計算されるエネルギー最小の形状とよく一致し、異なる形状へ遷移する際には ADE モデルに特徴的な階段状の軌跡を描いていることが明らかになった。これらの実験結果から、観察された変形過程においてベシクルは、形状から測定される	
 ∆a	
 とは異なる最適な膜面積差	
 ∆a0	
 を温度変化とともに増加させている事が示された。	
 
また、	
 (v,	
 ∆a0)	
 の値とベシクルの形状を対応させる相図から、実験で観察される形状の転移点における	 ∆a0	
 の値を求め、ベシクルの変形過程における	
 ∆a0	
 の初期条件と	
 ∆a0	
 の増加を見積もり、自己生産が行われる変形過程とそうでない変形過程を比較することで、自己生産過程を実現するのに必要な	
 ∆a0	
 の初期条件と	 ∆a0	
 の増加を明らかにした。これらの結果を実験に用いた脂質の幾何学的な形状を基にして議論し、自己生産に必要な構造パラメーターを抽出する。	
 	 【 参 考 文 献 】 	
 (1) A. Sakashita, N. Urakami, P. Ziherl and M. Imai: Soft Matter 8 (2012) 8569.

(2) Y. Sakuma and M. Imai: Phys. Rev. Lett. 107 (2011) 198101.

Copyright (C) 2013 ソフトマター研究会 All Rights Reserved.

12月13日
16:30-16:45

脂質膜上の分子流動性とリン酸化分子反応ネットワーク（独）物質・材料研究機構	
  貝塚芳久【はじめに】
脂質膜上の分子自己組織化は、異なる分子流動性を持つ複数種の相を形成することが広 く 知 ら れ る 。こ れ ら 各 相 の 分 子 構 造 が 生 命 機 能 に 果 す 役 割 を 明 ら か に す る た め に は 、それぞれの構造で起こる分子の反応についても考える必要がある。本研究では、分子リン酸化を行う酵素タンパク分子が脂質二重膜上で流動性の異なる分子クラスター構
造内で異なる反応を行う、とするモデルを実験的に示すことを目的とする。

【結果と考察】
脂質二重膜上で異なる分子流動性を持つリン酸化酵素分子の反応を蛍光顕微鏡によりイ メ ー ジ ン グ・計 測 し た 。こ の 酵 素 分 子 は ア ミ ノ 酸 チ ロ シ ン 基 を リ ン 酸 化 す る が 、「 分子間反応」と「分子内反応」を起こす。実験結果は、分子流動性の低い固相・ゲル相に近い状況の酵素分子は「分子内反応」を起こし、分子流動性の高い相の酵素分子では「分子間反応」が優位であることを示した。この酵素分子のリン酸化作用は分子間反応が本来優位である（リン酸化速度が速い、分子三次元構造による熱力学的要因）
が、分子流動性が下がれば分子内反応も起こる、として説明できる。この結果と関連する研究結果を含めて、以下のように考察している：（１）細胞膜で酵素分子の拡散が制限されるクラスターが形成されるが、その内部でも我々が人工的なシステムで観察したものと同様のリン酸化反応が起こっていることがわかった。
（２）「分子内反応」と「分子間反応」により、酵素分子の活性が変化することが、生化学的解析により明らかにされている。他の「脱リン酸化酵素」による酵素分子の脱リン酸化は、酵素分子の活性を上げる効果があるとする実験結果を得たが、酵素の反応速度論による計算で実験結果を再現した。
（３）この実験・計算結果は、即ち「脱リン酸化酵素」を系に含めても、リン酸化酵
素の活性が上がることで、系全体のリン酸化が上昇する、ことを示すが、これは遺伝学的な結果と一致している。本研究では、この現象を生化学的に初めて示すことができた。

Copyright (C) 2013 ソフトマター研究会 All Rights Reserved.

12月13日
16:45-17:00

リン脂質ベシクルの化学駆動（東北大・院理・物理）兒玉 篤治、佐久間 由香、今井 正幸、川勝 年洋【はじめに】
近年、化学的に非対称な環境に置かれた油滴やコロイド粒子の自己駆動現象が相次いで報告され注目されている。そのような状況の中、我々はリン脂質の一つである
1,2-dioleoyl-sn-glycero-3-phosphocholine(DOPC)より調製したジャイアントベシクルが特定のイオン種を含む溶液のマイクロインジェクションによりその chemical source に向かって集まる chemophoresis 現象を示すことを見出した。また、ベシクルがある程度の余剰面積を有する場合、ベシクルの移動とカップルして劇的に変形する様子も観察された。今回、我々は特に水酸化ナトリウム(NaOH)水溶液をインジェクションした場合について、このベシクルの移動、変形に注目して実験、解析を行った。
【結果と考察】NaOH 水溶液をベシクル溶液にインジェクションするとベシクルは次第に速度を増しながら
ほぼ直線的にマイクロピペット先端に向かって移動した。NaOH 水溶液の濃度が 10 mM、直径 10 μm 程度のベシクルについてその速度を求めたところ、数 μm/s から数十 μm/s であり、特にマイクロピペット近傍では 100 μm/s 近くに達するものもあった。我々は実験から得られたベシクルの速度プロファイルとインジェクションを考慮した拡散方程式から求めた OH-イ
オンの濃度プロファイルを比較した結果、相関があると判断し、ベシクルの界面エネルギーが系の OH-イオン濃度に依存し、その界面エネルギーが小さくなる方向にベシクルは移動するというモデルを立てて検証を行った。その結果、実験で得られたベシクルの速度とモデルから得られたベシクルの速度は比較的良く一致することが分かった。このモデルではベシク
ルの先端に大きな力が働くことが示唆されるが、このような状況ではベシクルからチューブが引き出されることが期待される(1)。我々はある程度大きな余剰面積をもつベシクルに対して NaOH 水溶液インジェクションしたところ、チューブ形成とマイクロピペット近傍でそのチューブが伸びる様子が観察され、実際にベシクル先端に大きな力がかかっていることが示唆された。また、ベシクルがマイクロピペットから比較的離れている場所では Prolate 型のベシクルなどが数秒のうちにパーリングを示すような様子も観察された。このような状況ではベシクルの周りの OH-の濃度勾配は小さいため、ベシクル先端にかかる力も小さいと考えられる。この場合、ベシクル膜の外葉のみがインジェクション溶液に曝されて膜外葉の面積が増加する。その結果、膜の内側と外側の面積差が増加することでベシクルのパーリングが誘起されたと考えた。発表では NaOH 溶液以外の結果についても触れる。
【参考文献】(1) Derényi, I.; Jülicher F.; Prost J. Phys. Rev. Lett. 2002, 88, 238101.

Copyright (C) 2013 ソフトマター研究会 All Rights Reserved.

12月13日
17:00-17:15

球状ベシクルに内包されたベシクルの形状決定機構の解明
お 茶 大 院 理 A, 東 大 物 性 研 B,東 北 大 理 C
坂 下 あ い A, B, , 今 井 正 幸 C, 野 口 博 司 B	 	  近 年 生 体 膜 研 究 は よ り 複 雑 な 系 へ と 展 開 さ れ て い る 中 、我 々 は ミ
ト コ ン ド リ ア に 注 目 し た 。ミ ト コ ン ド リ ア は ほ ぼ 全 て の 細 胞 に 含 ま
れ る 細 胞 小 器 官 で あ り 、外 膜 と 内 膜 の 二 重 の 生 体 膜 か ら 成 る 。主 な機 能 は 生 体 活 動 に 必 要 な エ ネ ル ギ ー の 生 産 で 、内 膜 に 特 殊 な 陥 入 構
造をとることで生産効率を高めている。
	  従 来 こ の 陥 入 構 造 は タ ン パ ク 質 に よ っ て 制 御 さ れ て い る も の と
考 え ら れ て き た が 、我 々 は 外 膜 の 存 在 こ そ が 形 状 決 定 の 鍵 を 握 っ ていると考える。
	  そ こ で 本 研 究 で は 、球 状 ベ シ ク ル に 他 の ベ シ ク ル が 内 包 さ れ て い
る マ ル チ ラ メ ラ ベ シ ク ル の 変 形 過 程 を 実 験 的 に 撮 影 し 、そ れ を 元 に
シ ミ ュ レ ー シ ョ ン を 用 い て 形 状 決 定 機 構 の 解 明 を 目 指 し た 。我 々 が用 い た 手 法 は 動 的 三 角 格 子 モ デ ル [1]で あ り 、 現 在 実 験 と 最 も 良 く
一 致 す る ADE model [2,3]の 効 果 を 取 り 入 れ て い る 。
	  そ の 結 果 、類 似 研 究 [ 4 ] に 報 告 さ れ て い る 二 重 の s t o m a t o c y t e 状の 陥 入 構 造 に 加 え て 、新 た に ミ ト コ ン ド リ ア で 見 ら れ る よ う な ひ だ
状 の 構 造 や 、細 胞 分 裂 の 初 期 に 見 ら れ る 胚 の よ う な 区 画 構 造 の 再 現に 成 功 し た 。ま た 実 験 結 果 を 元 に 再 現 し た シ ミ ュ レ ー シ ョ ン 結 果 は 、ADE model を ベ ー ス に 独 自 に 開 発 し た 理 論 モ デ ル と も 非 常 に 良 い一致を示したため、合わせて報告する。
文献[1] H. Noguchi and G. Gompper, Phys. Rev. E 7 2 , 011901 (2005).
[2] S. Svetina and B. Žekš, Eur. Biophys. J 1 7 , 101 (1989).
[3] A. Sakashita, N. Urakami, P. Ziherl, and M. Imai, Soft Matter 8 , 8569
(2012).
[4] O. Kahraman, N. Stoop, M. M. Müller, Nes J. Phys. 1 4 095021 (2012).

Copyright (C) 2013 ソフトマター研究会 All Rights Reserved.

12月13日
17:15-17:30

直流電場に誘起されるマイクロサイズ液滴の周期運動（ 京大院理, 東工大総理工 A ,同志社大生命医 B ）栗村朋, 瀧ノ上正浩 A ,吉川研一 B ,市川正敏【はじめに】
ナノ〜マイクロメートル程度のスケールにおいては粘性が支配的であり、熱ゆらぎの効果も大きい。それ故、objective lensof microscope

マクロな機械を単にサイズダウンしたものは散逸が大mineral oilwater droplet

きくエネルギー効率が悪いため実用に適さない。一方で、
V

バクテリアの回転べん毛モーターなどはこの様なスケ

Lglass plate  
Schematicールでもゆらぎに抗して自律的な周期運動をしている。

Fig.

このようなスケールにおいて効率的に駆動する機械を

representation

つくるためには、従来の動作原理からまったく異なる運experimental setup.

１

	
 

of

the

動機構で構成する必要があるだろう。	
 	
  最近、我々により sub-mm スケールにおいて油相中に界面活性剤で安定化した水滴が DC
電圧下で周期運動することが報告されている(1)(2)。今回この系をサイズダウンすることにより、低電圧で動作させ、その運動機構を解明した。	
 【結果と考察】Fig. 1 のようにマイクロメートルスケールの水滴を油中に作成し、直流電場の中での振る舞いを、スケール依存性に特に着目して調べた。電極間距離 L を固定して電圧を上昇させると、静止していた液滴は動き始め、電極間を往復運動する。このとき、静止状態から往復運動に移る閾
値電圧が、電極間距離 L の 3/2 乗に比例することが明らかになった。往復運動の数理モデルを考案し、実験と対比しながら理論的な考察を行い、リミットサイクルへのホップ分岐で実験を説明出来る事が分かった。さらに、電圧にノイズを加えた場合に、振動状態の拡大と安定化が起こることがわかった。(3)	
 	  以上のように、ノイズに強く、利用も出来るこの液滴運動の性質は、マイクロメータースケールの物体を利用した駆動機械の新しい動作原理として有望であることがわかった。この系を更にサイズダウンして、その性質を明らかにする事が、今後の課題である。	
 	 L = 141µm
t=0 [s]0

V = 13.0V

V = 13.7V
1

2

3【参考文献】(1)	
 Hase, et al., Phys. Rev. E. 74, 046301(2006).
(2) Takinoue, et al., Appl. Phys. Lett. 96, 104105 (2010).
(3)	
 Kurimura, et al. Phys. Rev. E. 88, 042918 (2013).

100μm

Fig.2 Spatio-temporal diagram

Copyright (C) 2013 ソフトマター研究会 All Rights Reserved.

12月14日
09:00-09:40

温度勾配下での巨大紐状分子の輸送と秩序形成（ 京大白眉センター、科学技術振興機構さきがけ）前多	  裕介【 は じ め に 】 温度勾配の下で分子が高温側から低温側へと輸送される現象を Soret 効果（Ludwig-Soret effect）という。Soret 効果は古くから知られているが、近年になって荷電コロイドや荷電高分子を用いた実験と理論の両面から精力的に研究が進められている(1)。本講演では、PEG 高分子溶液中における DNA のソーレ効果を中心に、輸送のメカニズム、パターン形成、分子操作技術への展開について述べる。
【結果と考察】	  ポリエチレングリコール(PEG)高分子ポリマーを体積分率 1~5%で添加した溶液中の DNAの Soret 効果を解析した。1480 nm 赤外線レーザーを集光し、最大温度差５℃、典型的な温度勾配 ∇T =2.5 K/µm の温度分布 T (r) を形成し、蛍光色素でラベルされた DNA の濃度分布を蛍光顕微鏡で計測した（図 A）。定常的な温度勾配の下で DNA は低温側に輸送され、c(r) = c0 exp "#−ST (T (r) − T0 )$% に従う濃度分布を形成する。ここで PEG も同様に Soret 効果
で輸送され濃度勾配をうみだすため、DNA は Soret 効果のみならず PEG 濃度勾配に起因す
る拡散泳動という効果を受ける(2)。２つの輸送が拮抗し、溶液中の PEG 濃度が上昇するにつれて、DNA はリング局在(2%~4%)または高温側へ濃縮される(4%~5%)こと（図 B）、大きさの異なる分子をふるい分けることを見いだした（図 C）。DNA の熱泳動と拡散泳動を考慮すると、定常状態の濃度分布は c(r) = c0 exp "#−ST (T (r) − T0 ) + c0PEG − c PEG (r) V $% と表すことが()

でき、実験で見られた振る舞いを説明することが出来る（3, 4）。
	  また、高分子溶液中のソーレ効果は、
対象物の電磁気学的な性質によらずに生体高分子や細胞などを捕捉する特性をもつ。これまでに赤外線レーザーをス
キャニングすることで、２次元面上でDNA、タンパク質、細菌などを様々なパターンへと捕捉・制御すること、時間的にも変調させることに成功している(5)。講演の最後に本技術について紹介すると共に、ソフトマターのソーレ効果
について多角的な議論を行いたい。
図. 高分子溶液中のソーレ効果
A. 実 験 の 模 式 図 , B. DN A の リ ン グ 局 在 と 濃 縮 , C. RN A を Soret 効 果 で 分 離【参考文献】(1) D Braun and A Libchaber. Phys. Rev. Lett. 89, 188103 (2002)(2) HR Jiang, H Wada, N Yoshinaga, M Sano. Phys. Rev. Lett. 102, 208301 (2009)(3) YT Maeda, A Buguin, A Libchaber. Phys. Rev. Lett. 107, 038301 (2011)(4) YT Maeda, T Tlusty, A Libchaber. Proc. Natl. Acad. Sci. USA 109, 17972 (2012)(5) YT Maeda. Applied Physics Letters, in press (2013)

Copyright (C) 2013 ソフトマター研究会 All Rights Reserved.

12月14日
09:40-10:20

濃度勾配のある高分子溶液中で熱泳動により駆動されるコロイドの輸送（ 産総研、ナノシステム研究部門）関	  和彦【 は じ め に 】最近、熱勾配下で温度の高い方から低い方へ熱泳動するコロイドが、ポリマー溶液中では温度の高い方へ集まる現象が観測されている。[1]さらに別の条件では、コロイド
は温度の最も高いスポットではなくその周辺にリング状に集まることも観測されている。[2]
これらの実験ではレーザーが用いられており、局所的に温度勾配が作られ、温度やコロイド
およびポリマー濃度がマイクロメートルスケールで測定されている。[1-3]実験結果は、高分子濃度とコロイドのサイズを軸とした２次元の相図として整理されており、コロイドの濃度
分布のパターンを統一的に示している。本講演では、輸送方程式を導出し、得られた方程式
がリング形成を含めたコロイドの空間分布を再現することを示す。さらに、相図の理論的な解析から、熱泳動の強さの指標である Soret 係数のコロイドサイズ依存性についての知見が得られることを示す。	
 【 結 果 と 考 察 】 熱泳動により駆動されるコロイドに対する高分子の濃度勾配の効果として、排除体積相互作用を考慮した。コロイドと濃度勾配を伴って分布している粒子との排除体積相互作用の効果を理論的に考慮するために、コロイドと高分子のサイズが同じであると仮定した格子モデルを用いてコロイドに対する輸送方程式の導出を行なった。その結果コロイドの定常流 J は、コロイドの濃度 c(r, t) と高分子の濃度 p(r, t) を用いて以下の式で表すことができた、	
 J = #$ D {1− p(r)}{∇c(r) + c(r)STc ∇T (r)} + c(r)∇p(r)%& 。	
  	
  	
  	
  	
  	
  	
  (1)	 ここで、 D はコロイドの拡散係数、 STc は Soret 係数を表す。この式は、高分子による排除体
積相互作用のために、コロイドの拡散係数は因子 1− p(r) 減少し高分子の濃度勾配 ∇p(r) によるドリフトが引き起こされることを示している。実験で得られているコロイドの空間分布は、リング状の空間分布も含めてこの式から得ることができる。さらに、コロイドのサイズと高分子のサイズが異なっている場合に拡張するために、(1)式をエントロピーと化学ポテンシャ
ルを用いた輸送方程式から再導出した。この方法を用いて、コロイドのサイズと高分子のサイズが異なる場合の定常流の拡張式を導出し、その結果を相図に表し実験結果と比較する事
により、Soret 係数のコロイドサイズ依存性についての考察を行なった。	
 【参考文献】(1) H. R. Jiang, H. Wada, N. Yoshinaga and M. Sano, Phys. Rev. Lett. 102, 208301 (2009).

(2) Y. T. Maeda, A. Buguin and A. Libchaber, Phys. Rev. Lett. 107, 038301, (2011).
(3) Y. T. Maeda, T. Tlusty, and A. Libchaber, Proc. Natl. Acad. Sci. USA 109,
17972 (2012).
(4) K. Odagiri, K. Seki and K. Kudo, Soft Matter, 8, 2775 (2012).

12月14日
10:40-11:20

高分子溶液の温度勾配下における濃度勾配形成現象（Ludwig-Soreteffect）の測定と解析（東海大理）喜多 理王【はじめに】
混合流体に外場として安定な温度勾配を与えると, 温度勾配をドライビングフォースとす
る成分の拡散が生じ安定な濃度度勾配が形成される。これはルードヴィッヒ・ソレー効果とも呼ばれ, 2 成分系では成分 1 の流束 J1 は, 濃度勾配∇c1 と温度勾配∇T を用いて現象論的方
程式

J1=－D∇c1－c1DT∇T によって記述される[1]。高分子やコロイドなどソフトマターのルードヴィッヒ・ソレー効果は, 測定手法の進歩とともに現象を特徴づける輸送係数が得られるようになってきたが, 複雑な分子間相互作用を有する水溶性高分子やミセル系などでは測定が難しいことが多く また結果の解釈が困難な場合もある。本報告では PNIPAM や PEO
また生体高分子において得られたルードヴィッヒ・ソレー効果の最近の研究成果を説明する。
【結果と考察】
温度勾配による単一の粒子の輸送速度 v1 は, 熱（物質）拡散係数 DT を用いて v1=－DT∇Tであるが, 高分子鎖の v1 の直接計測は難しい。そこで高分子溶液では濃度勾配の大きさを計測することが多い。また濃度勾配の形成過程の解析から自己拡散係数 D が得られる。定常状
態ではみかけ上流束が消失し(J1=0), DT と D の比が DT / D =‐c1)(dc1/dT)となり, またソレー
係数 ST が ST≡DT / D と定義される。したがって ST と D の測定結果から DT が得られる。気体
子運動論的解釈では, 成分の運動量や運動エネルギーから, 重い成分が温度勾配の低温側へ偏ることが導かれる。溶液では成分の質量やサイズに加えて, 分子間相互作用やエネルギーの散逸に関与する粘性率などが濃度勾配の大きさと向きに決定的な役割を果たすことが分かってきた。最近では電荷や流体力学的相互作用などの種々の効果が調べられている。我々はPNIPAM や PEO 溶液また界面活性剤や多糖類において, レーザー干渉計（TDFRS 法[2]）を用いた測定により ST , DT , D の濃度, 温度, 分子量依存性などを調べてきた[3-5]。その結果,(1) PNIPAM 水溶液のΘ温度近傍における濃度勾配の亢進現象, (2) オリゴ糖や多糖類が水溶
液中で温度勾配の高温側へ偏るという異常な（負の）ルードヴィッヒ・ソレー効果, (3) 尿素など分極の大きな分子を溶液に加えたときの水構造の変化とそれに対応した濃度勾配形成方向の反転現象, (4) 合成高分子とは異なる分子量依存性など, 水素結合が関与していると考えられるソフトマテリアルに特有なさまざまな現象が見出されてきた。発表ではこれらの中から興味深い結果を報告する。
【参考文献】(1) S. R. de Groot and P. Mazur, Non-equilibrium Thermodynamics. (Dover, New York, 1984).

(2) W. Koehler and R.

Schaefer, New Developments in Polymer Analytics II; vol 151, 1-59 (Springer: Berlin, 2000).

(3) B. J. de Gans, R. Kita, S.

Wiegand and J. Luettmer-Strathmann, Phys. Rev. Lett. 91, 245501 (2003).
4554 (2005).

(4) R. Kita and S. Wiegand, Macromolecules 38,(5) Y. Kishikawa, H. Shinohara, K. Maeda, Y. Nakamura, S. Wiegand and R. Kita, Phys. Chem. Chem. Phys.

14, 10147 (2012).

Copyright (C) 2013 ソフトマター研究会 All Rights Reserved.

12月14日
11:20-12:00

温度勾配下での荷電コロイドの結晶化
名古屋市立大学大学院 薬学研究科

山中淳平, 豊玉彰子, 奥薗透【はじめに】
表面に電荷を持つ「荷電コロイド粒子」は, 粒子間の静電相互作用が十分大きいとき,
液体媒体中で規則正しく配列して「結晶」構造を形成する.コロイド系には, (1)構成単位である一粒子を, 光学顕微鏡によりその場・実時間観察できること, (2)粒子間相互作用の大きさが広範囲に調節できること, (3)特性時間が長いため,動的過程の観察が容易であるなどの利点があり, 相転移のモデル系として注目を集めてきた.また近年では,光学材料などへの応用も盛んである. 原子・分子系の単結晶育成 (ブリッジマン法など) では, 温度勾配下での一方向凝固が用いられる. 一方, 荷電コロイドの結晶化の駆動力である静電相互作用の大きさは, 主に粒子の体積分率φ, 粒子の電荷数 Z, 媒体の塩濃度 C などで決まり, 温度 T は弱い変数である.しかし C, φ, Zのいずれか温度依存性を利用すれば, 温度により結晶化を制御できる 1. 本講演では Z の温度依存性を利用した荷電コロイド系の温度誘起結晶化および温度勾配下での一方向結晶成長 2,3 を報告する.
【結果と考察】シリカ（SiO2）微粒子の電荷数は高 pH ほど大きい.また, 弱塩基であるピリジン Pyの解離度は高温ほど大きい. したがって Py を添加したシリカコロイドでは, 昇温により Z が増加し, 適切な条件では結晶化する（熱運動の寄与を静電相互作用の増加が上回る）.Fig.1 に一方向結晶成長の例を示す 3．試料を石英セル(0.114.5cm)に入れ, 左端をヒーターブロックに接触させて温度勾配を与えている（a：φ＝0.035,[Py] ＝ 35μM, b: φ＝0.05,[Py] ＝ 27.5μM.ヒーター温度＝60℃）.結晶成
長曲線は熱伝導により説明できた.荷電コロイ
ドのゾーンメルト , 不純物排除などについても報告する.

Fig.1 Unidirectional crystallization ofcolloidal silica +Py dispersions undertemperature gradients.

【文献】
1) J. Yamanaka, T. Koga, H. Yoshida, N. Ise, T. Hashimoto, Slow Dynamics inComplex Systems; Tokuyama, M. and Oppenheim, I. eds., Woodbury, NewYork (1998), p. 144.
2) A.Toyotama, J.Yamanaka, M.Yonese, T.Sawada, F.Uchida, J.Am.Chem.Soc. 129(2007) 3044.
3) A. Toyotama, J. Yamanaka, M. Shinohara, S. Onda, T. Sawada, M. Yonese, F. Uchida, Langmuir, 25(2009) 589.

Copyright (C) 2013 ソフトマター研究会 All Rights Reserved.

12月14日
12:00-12:15

温度勾配下における二分子膜ラメラ相のダイナミクス（首都大学東京）栗田玲（東京大学生産技術研究所）田中 肇【はじめに】
非イオン界面活性剤・水系における二分子膜は、お互いにぶつかり合うことによって安定化している。そのため、ラメラ相のエネルギーは温度に敏感である。この二分子膜に温度勾配を与えると、大きなエネルギー勾配が形成される。そのエネルギー勾配に対して、どのような挙動を示
すかを実験的に調べた。

【結果と考察】結果と考察】

ラメラ相のエネルギーを最適化するためには、膜間距離を変更する必要があるが、壁にホメオトロピックに配向したラメラ相は膜という幾何学的な拘束を受けている
ため、簡単には膜間距離を変更することができない。そのため、膜が折れ曲がっている線欠陥を移動させることによって、膜の枚数の空間分布を変化させていることを見つけた。さらに、線欠陥の向きと温度勾配に相関を調べるため、くさび形セルを用いた。くさび形セルを用いることにより、一方向の線欠陥だけを選択的に形成することが可能である。その結果、線欠陥と温度勾配に特異的な関係を見つけることに成功した。これらの結果について、くわしく発表する。

図：左側が３５度、右側が２５度の温度勾配を与えたときの線欠陥の移動。温度勾配を与えてからの時間は、左図から０分、３３０分、６００分。

Copyright (C) 2013 ソフトマター研究会 All Rights Reserved.

12月14日
16:00-16:15

光で設計した
設計した共連続構造
した共連続構造の共連続構造の傾斜を
傾斜を有するポリマーブレンド
するポリマーブレンドにポリマーブレンドに
関する研究
する研究
京工繊大・院工宮田 貴章

自然界は熱力学的な開放系であり、常に平衡から遠く離れているか、非平衡のまま凍結される場合が多い。そのため、自然界で見られている構造がほとんどの場合、不均一性を示している。具体的に言えば、傾斜構造や階層構造のような空間的分布を有し、これらの特異的な構造から均一系に見られない機能が発現すると考えられる。高分子、液晶やミセルなどのソフトマターの立場から考え、上述した外界の不均一的な刺激に対する高分子系の応答性を明らかにし、これにより新たな機能を示す高分子材料の設計を目的とする。私たちは以前から非反応系の相分離に対する温度勾配の影響[11]や反応系ポリマーブレンドの相分離に及ぼす光強度の勾配の役割について検討してきた[22]。熱分解を避けるために、光強度の勾配を制御パラメータとして用いてきた。ここで、光強度の勾配で誘起した相分離過程と、生成したモルフォロジーについて紹介する。さらに、傾斜構造の応用として非対称性の表面を有する高分子の有効利用についても述べる。
Projector
実験方法と
験方法と結果本研究の基本的なアイデアは高分子の重合反応を光で誘起し、光強度を変化させることにより、高分子/モノマーの混合物
Computer 1
は相溶（一相）領域から二相領域に突入させ、相分離散を誘起する。その際、弱い光の強度で相分離を引き起こすと、混合系
Lens
Computer 2が浅くクエンチされ、準安定領域に突入し、核生成・成長が起こる。一方、強い光で相分離を引き起こす場合、混合系は深くクエンチされ、スピノーダル分解が起こり、重合の収率（高分子の生成量）によって、相分離が途中で停止され、共連続構造
Microscope
CCD cameraが得られる[22]。従って、光強度と Lambert-Beer 則を組み合わせ
Fig.1 Computer-assisted
ると、照射光強度の勾配に沿って、クエンチ深さの傾斜を変化
irradiation(CAI) apparatus
させることができ、スピノーダル構造や生成した核や dropletsの特性長に空間的な傾斜が形成される。さらに、Z 方向に沿って光を試料に入照射させる場
合、光の進行方向に共連続構造や droplet 構造の傾斜が形成される[33]。光強度を変えることによりこの構造の空間傾斜を変化させることができる。さらに、Z 方向の傾斜が無視でき、
XY 平面内のみの傾斜構造を作製するために、当研究室で開発したコンピューター支援光照4]を利用することにより、Z 方向の光強度の傾
射 (Computer-Assisted Irradiation, CAI 法) 法[4斜を無視できるくらい小さく、XY 平面内に光強度の傾斜を高くする条件下で相分離を引き起こすと XY 平面内のみモルフォロジーの傾斜を発現させることができた。その一例を図図２に示す。すなわち、このようにして、二方向（Y と Z）の傾斜共連続構造を有する高分子を設計できた。
このように CAI 法を重合誘起相分離と用い合わせると、ミクロン域において予めプログラムした任意の傾斜構造が設計できるだけでなく、室温で実験を行なわれるので、温度勾配を利用した場合と比べて、高温側の熱分解などの問題を避けることができた。現在、CAI 法で設計した様々な相分離構造を有効に利用するため、空間的傾斜や異方性の共連続構造を有する高分子の表面的特性やナノフィラーの選
択分散の特性について実験的に検討している。詳細は当日に発表する。

References. [1] J. Okinaka, Q. Tran-Cong Physica D 84, Fig. 2. PEA/PMMA mixture withbi-directionally graded co- continuous
23-30 (1995); [2] H. Nishioka et al. Macromolecules 33,
morphology.
4301 (2000); [2] N. Kimura et al. Soft Matter 9, 8428
(2013); [3] H. Nakanishi et al. Soft Matter 2, 149 (2006); [4] Ishino et al. Macromol. Rapid
Commun., 27, 758 (2006).

Copyright (C) 2013 ソフトマター研究会 All Rights Reserved.

12月14日
16:15-16:30

拘束や温度勾配がある系の multiplicative なブラウン運動（筑波大学物理・名古屋大学理学研究科）黒岩健

宮崎州正○

本講演では、manifold 中に拘束された系や、温度勾配などの非平衡条件下におかれた系のブ
ラウン運動に対するランジュバン方程式について議論を行う。
ランジュバン方程式は、揺らぎがある系のダイナミクスを記述するための代表的な手法であり、あらゆる科学の分野で広く使われている。ランジュバン方程式は、物理量の時間発展を、遅い変数で駆動される部分とランダム力の和として表すだけの単純なものであるが、その導出は必ずしも自明ではない。微視的な運動方程式から粗視化による導出や、現象論的な推論により書き下すなどの方法があるが、いずれにしても遅い部分とランダム力をどのように分離するかは自明ではない。特にランダム力自身が遅い変数にも依存する場合（これをmutliplicative ノイズと呼ぶ）には、確率過程における積分の定義にもかかわる数学的な解釈も絡み合って、さらに難しい問題となる。この問題は、代表的な確率積分名にちなんで、
Ito-Stratonovich ジレンマと呼ばれている。
この multiplicative ノイズをもつ系は、ソフトマターの世界に満ち溢れている。たとえば
高分子鎖中のモノマーや膜中のたんぱく質のブラウン運動は、運動に制約があるという点で拘束条件付きのブラウン運動であるが、これはノイズに強い変数依存性を課することになるため multiplicative となる。また、温度勾配や濃度勾配がある系などはノイズが場所に依存するから当然 multiplicative である。いずれも単純な例題であり多くの研究が報告されているが、それらを記述するための「正しい」ランジュバン方程式の導出についての議論は収束していない。
本公演では我々は、拘束条件や非平衡条件下における「正しい」ランジュバン方程式を導出する手法について報告すると同時に、その適用限界や問題点について皆さんから意見をいただきたい。
【参考文献】(1) T. Kuroiwa, K. Miyazaki, arXiv:1309.4189,

Copyright (C) 2013 ソフトマター研究会 All Rights Reserved.

12月14日
16:30-16:45

温度勾配下のコレステリック液晶滴における配向／重心回転（早大理工）吉岡【はじめに】

潤、鈴木

雄仁、伊藤

文哉、多辺

由佳

コレステリック(Ch)液晶に対して温度勾配を印加すると配向方向が一方向に回転する現象は Lehmann 回転と呼ばれており、その機構は Leslie の理論によって現象論的に説明されている[1]。しかしながら実際に Lehmann 回転を観測した実験の再現性は極めて低く、未だ十分な実験的検証がなされているとは言えない。近年、この Lehmann 回転の存
在を示唆するものとして、等方(I)相‐Ch 相共存領域(I+Ch)において系に温度勾配を印加することによって、等方液体媒質中に分散している Ch 液晶が形成する液滴（Ch 液晶滴）上の縞状組織が回転する現象が報告された[2]。しかしながらこの現象の実態に関しては未だ不明な点が多く、実際に Leslie の理論を用いてその回転運動の機構を説明できるかどうか判別がついていない。本研究では上記の回転運動の機構を解明するため、偏光顕微鏡観察および蛍光褪色回復法による流動場測定によって、回転運動のダイナミクスを解析した。
【結果と考察】 本研究では、Ch 液晶試料としてネマチック液晶混合物（5CB-No270032 混合系、LCC 社）にカイラル添加剤 R811 を 0.5~2.0wt.%程度添加したものを用いた。試料の相系列
は R811 の濃度にあまり依存せず、I-58℃-I+Ch-54℃-Ch である。上記の試料を用いて、I 相からの降温によって I-Ch 共存領域において Ch 液晶滴を作成したところ、図 1 に示すように縞状、同心円状の組織を有した 2 種類の液晶滴が発現し（縞型、同心円型液晶滴）、温度勾配を印加する
ことでどちらの液晶滴においても回転運動が観測された。上記の 2 つの液晶滴における回転運

動の特性を解析するために、我々は系のカイラリティの変調に伴う回転速度の変化を解析した。図 2 に印加した温度勾配T から回転運動（角速度）への変換効率T（回転効率）の
カイラリティ依存性を示す。同心円型液晶滴においては、
Lehmann 回転と同様に、系における R811 の濃度が増加
しカイラリティが強まるにつれて回転効率が増加する。一方、縞型液晶滴においてはカイラリティが強まるにつれて回転効率が減少していく Lehmann 回転とは正反対の挙動が観測された。さらに、蛍光褪色回復法を用いて系の流動場を解析したところ、縞型液晶滴においては縞状組織の縞型
同心円型
図 1：Ch 液晶滴の偏光顕微鏡像回転に伴った回転流が発生していることが判明した。以上の結果より、我々は同心円型液晶滴における回転運動は
Lehmann 回転と同様の配向回転、縞型液晶滴における回転運動は Lehmann 回転では説明できない液晶滴自体の回転である、と結論付けた。詳細は発表で述べる。
【参考文献】(1) S. Chandrasekhar, “Liquid Crystals”,
Cambridge University Press 2nd edition (1992) (2) P.
Oswald and A. Dequidt, Phys. Rev. Lett. 100, 217802(2008)

図 2：回転効率のカイラリティ
依存性（○印：同心円型液晶滴）（□印：縞型液晶滴）

Copyright (C) 2013 ソフトマター研究会 All Rights Reserved.

12月14日
16:45-17:00

薄いキラル液晶セルが示す秩序構造の光学的性質（ 産総研ナノシステム）福田順一（リュブリャナ大学）Slobodan Žumer
【はじめに】
	  薄い平行平板セルに強いキラリティを有する液晶を閉じ込めると，様々な特異な秩序構造を示すことを，我々は連続体理論に基づく数値計算によって明らかにしてきた[1]．本研究では，それらの秩序構造がどのような光学的性質を示すか，より具体的には，垂直な単色入射
光がどのような強度分布で反射されるかを，平面波展開に基づく数値計算によって調べた[2]．
【結果と考察】	  図に示しているのは，特異な秩序構造の一例である，液晶配向の線欠陥の 2 重らせん構造と，そこに垂直に無偏光の単色光が入射した際の反射光の強度プロファイルの計算例である．
反射光の強度プロファイルは入射光の波長に依存するが，ここで示したケースにおいては，欠陥の２重らせんのらせん軸方向の強度変調は急激に減衰し，セルから離れたところでは強度プロファイルがほぼ１次元的なものになっていることがわかる．
	  他の構造に関する計算結果については，講演で紹介する．これらの構造の格子定数は可視光の波長程度（数百 nm）であり，それらの実空間観察は決して容易ではないが，スロベニ
アの実験グループが実験を既に試み始めており，これらの計算結果がそのような実験検証を促進し，実験結果に対する理論的基盤を与えることを期待している．
	  本研究は Slovenian Research Agency, COE“NAMASTE”(Slovenia)，および科研費基盤 C (25400437)の支援を受けている．

図: 左は薄い液晶セル中の秩序構造の一例である２重らせん配向欠陥（黒）．右はそれぞれ，液晶セルの光が入射した側の面から 0, 1, 3 だけ離れた面内における反射光の強度プロファイル．ここで，長さはねじれ配向のピッチが 4π になるように規格化している．セルの厚さは 12, 入射光の波長は 8.84 である．

【参考文献】[1] J. Fukuda and S. Žumer, Phys. Rev. Lett. 104, 017801 (2010); ibid. 106, 097801(2011); Nature Commun. 2, 246 (2011).
[2] J. Fukuda and S. Žumer, submitted to Mol. Cryst. Liq. Cryst.

Copyright (C) 2013 ソフトマター研究会 All Rights Reserved.

12月15日
09:00-09:15±³NS®

¡Usã ~}qtÄ¨µt´¸µä¯Ð ´HÕ

±³NS®0vÚbåçR³ÝÝ·\¡Us0++¾*È¼-0å

CW±³NS®
-JÖÃ°¦

0¢-]

Ú

ÝW0ÎÍ-±`!-¶skÝWå­Râ\ruv

ÙÅ^0.-«å±ªÑ
%!±³NS®
Ë.-oe¡(%-

*#.+

£Ç

ßeR

£
z

¡(QJ-+å7E:4Ñp)ßeRå±ÖPÜ-!¥bÀF36[m²±\rÜY±³NS®0wá²vÚb-!x!

*

Æ¦
+

|

*ÉV£L¡Usng-(1)%åÔ»ç;B1>7E\rÜY`³NS®
-'åÔ»ç;B1CW0+å(c¿±\r->7E0Tå±³NS®

¡Us0-0Ê& Ô»ç;B1CW B !å§ÞØÞ
_<9E5?C w(r)çR³ÝÝ h(r)%!l

+

B=
*K

-£ÑÜ

1
2

h(r)d =1
2

e w(r)/kT

1 d+.-åd!{R·ÅÀå·\!Z¹Ü/-R

1CWfåB < 0 ."_åB > 0 ."_-\ra_t5=@DG5AEÅ!å»Hå>7E£Òà rc %â¾('åI

·\$

rc %

»ç;B

×²-æOh+ h(r)0-Þ

yK0È¼å»Mår > rc

½i!

c¿Rº-¤Ó

h(r) =0®, ÛÒà

yK0È¼-å(*,å>7E½i

1
Ae r/ cos(krr-Ô»ç;B1CW B 0 1 jâ¾È¼-æ|åâ¡!qUs!R>7E»ç;B1CWz©0/Ï(¡ 238 K

ÁåQ¡

I(cÌ)X0-
å

! B > 0 ,å¡+q+ 373 Kæ

¡

!f2G8G-ædÂ¬
(1) K. Koga, Osmotic Second Virial Coefficient of Methane in Water. J. Phys. Chem. B, 117, 12619(2013).

Copyright (C) 2013 ソフトマター研究会 All Rights Reserved.

12月15日
09:15-09:30

ミセル水溶液における内部構造変化の複屈折測定（長岡技科大）伊藤 雅利，吉武 裕美子，高橋 勉【はじめに】
ひも状のミセルを形成する界面活性剤水溶液は，低いせん断速度では単一緩和 Maxwell モデル
で近似される線形挙動を示すが，臨界せん断速度以上において溶液内部のミセル構造が変化し非線形な挙動を呈するようになる．この現象は一般に SIS（Shear-Induced Structure）と呼ばれている．特に，CTAB/NaSal 水溶液では，流動初期に粘度が定常値の 100 倍程度まで増加するせん断粘ちょう化現象（Shear-hardening）や，流路内でせん断速度が不連続になるシアバンド形成（Shear-banding）
が報告されており，近年，精力的に研究が成されている(1),(2)．本研究では，回転型レオメータに設置した同心二重円筒型流路で，CTAB/NaSal 水溶液に一定せん断速度の流動を与えた場合における内部構造変化をクロスニコル観察法により観察した．流路内の観察と溶液のせん断応力挙動とを比較することで内部構造の変化が流動特性に及ぼす影響を検討する．
【結果と考察】本溶液にせん断速度 ̇ = 5s-1 を与えた際の応力測定結果を Fig.1 左図に示す．ひずみ量に応じて，①shear-thickening 域，②shear-thinning 域，③応力平衡域へと遷移していく様子がわかる．その際の流路隙間 2.0mm における複屈折分布の観察画像を Fig.1(a)~(f)に示す．各ひずみ量での複屈折（位相差）が色相として現れている．また画像の下部が一定速度で回転する内円筒壁面である．
流動開始とともに画像の色相が複屈折の強い色へと変化していくが，ひずみ量が 10 を超えると急激に複屈折が弱くなる．Fig.1(c)から，この複屈折の低下は高い配向度の内部構造が空間的に崩壊することに起因していると考えることができ，これは急激な応力低下の挙動とも一致している．その後，Fig.1(d)のように流路内の複屈折分布は非常に不安定かつ不均一な様相を示し，応力もゆるやかに減少していく．さらに，ひずみが大きくなるにつれて，高い配向度の領域が回転壁面側へと移動していくことから，Fig.1(e)(f)に示されるようにシアバンドを形成することで応力が平衡状態へと落ち着くと考えられる．

Fig.1

Stress behavior and visualized birefringence images in CTAB/NaSal solution.

【参考文献】
1) Fardin M.A., Lerouge S., The European Physical Journal E, 35, 9, 91 (2012).
2) Bueno V.L., Kohlbrecher J., Fischer P., Rheologica Acta, 52, 4, pp.297-312 (2013).

Copyright (C) 2013 ソフトマター研究会 All Rights Reserved.

12月15日
09:30-09:45

高分子網目における結合率と弾性率の関係性（ 東大物性研）西 健吾、野口 博司、柴山 充弘 (東大院工)酒井 崇匡【 は じ め に 】 	
 パーコレーション理論はゲル形成過程における網目構造のモデルとして用いられ、こ
こから予測されるゲル化過程に関する種々の物性値が現実とよく一致することが知ら
れ て い る 。し か し な が ら 、こ の 理 論 で は 網 目 構 造 を 解 析 的 に 扱 う 手 段 が な い た め 、種 々の物性の絶対値を求めることができず、コンピューターによる計算が欠点であった。
そこで本研究では格子上のランダムな網目における弾性率を近似的に算出する一般的な手続きを示す。
【結果と考察】格子において格子点間をランダムに結合率 p で結合させていくことを考える。正方格子の例を図 1 に示す。この結合率 が 低 い 場 合 系 を 覆 う 無 限 大 ク ラ ス タ ー は 存 在 し な い (図 1(a))が 、 結 合 率 が 高 く な る と 無 限 大 ク ラ ス タ ー が 出 現 し (図1(b))、 最 後 に は 全 て の 格 子 点 間 が 結 合 す る (図 1(d))。 こ のような一連の過程をボンドパーコレーションと呼び、ゲル化過程における網目構造のモデルとして用いられる。
本 発 表 で は 、 こ の 網 目 の 弾 性 率 を 求 め る た め に 、 Kirkpatr
ick1 や Feng2、 Sahimi3 ら の 伝 導 網 に 関 す る 議 論 を 高 分 子 網目へと適用する。その詳細については当日発表にて行う。

図 1. 各結合率における正方格子【シミュレーション方法】
	  今 回 は 結 合 確 率 p で 格 子 点 間 を WLC(Worm-Like-Chain)
で 結 ん だ 三 角 格 子 と 立 方 格 子 を 作 成 し た 。 格 子 状 に WLCを並べた初期網目では力の分布が偏っている。そこで全体のエネルギーが小さくなるように、網目の構造を最もエ ネ ル ギ ー 的 に 安 定 な 状 態 へ と 構 造 最 適 化 さ せ た (図 2)。
変形に際しても歪みに伴って安定状態が変化するので構造最適化を行いながら変形を進めた。
図 2 三角格子の構造最適化。左が最適化前、右が最適化後。
【結果と考察】1.0
	  図 3 に ダ イ ヤ モ ン ド 格 子 の シ ミ ュ レ ーション結果と本理論の結果を示す。図 3 よ0.8
り全結合率領域において本理論とシミュレ
Prediction : slope ~ 1.7-1.9ーション結果が一致していることがわかる。
This theory: slope ~ 1.58
0.6
Simulation: slope ~ 1.60
ま た inset に ゲ ル 化 点 付 近 の 弾 性 率 の 挙 動 を
1.58
図 示 し た が 、 本 理 論 で は G ~ |p - pc| 、 シ ミ
p-p [-]
Simulation
0.4
ュ レ ー シ ョ ン で は G ~ |p - pc|1.6 で 振 る 舞 っ てThis theory
お り 、 こ れ は de Gennes ら が 予 想 し て い る
0.2
Exact value:
pc=0.2488
G ~ |p - pc|1.7~1.9 に 十 分 近 い 。 こ こ か ら シ ミ ュThis theory:
pc=0.267
レーションと本理論が妥当であることがわ
0.0
か る 。 ま た 本 理 論 は 2D の 正 方 格 子 に 対 し
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1.0
ても全く同様に適用できることが確認でき
p [-]ているが、その詳細については当日の発表にて言及する。
	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  図 3 立方格子の初期弾性率【参考文献】
1

8
6

G/G0 [-]

4
2

0.1

8
6

G/G0 [-]

4
2

0.01

7

8 9

0.1

2

3

4

5

6

7

c

(1) Kirkpatrick, S., Rev. Mod. Phys., 1973, 45, (4), 574-588.
(2) Feng, S.; Thorpe, M. F.; Garboczi, E., Phys. Rev. B, 1985, 31, (1), 276-280.
(3) Sahimi, M.; Hughes, B. D.; Scriven, L. E.; Davis, H. T., Phys. Rev. B, 1983, 28, (1), 307-311.

Copyright (C) 2013 ソフトマター研究会 All Rights Reserved.

12月15日
09:45-10:00

均一な第一網目構造を持つ超高強度 DN ゲルの一軸伸長変形挙動解析
123（北大院生命科学 , 北大院先端生命 , 東大院工 ）松田昂大 1, 中島祐 2, 黒川孝幸 2, 野々山貴行 2, 酒井崇匡 3, 鄭雄一 3, グン剣萍 2【はじめに】
DN ゲルは、大きく膨潤した第一網目鎖と延性の大きい第二網目鎖からなる相互侵入網目構造を持つ極めて強靭なハイドロゲルで、一軸引張変形において降伏挙動やヒステリシスを示す。強靭化の原因であるエネルギー散逸は第一網目構造の内部破壊に起因するが、従来の DN
ゲルの第一網目は極めて不均一で、内部破壊過程を評価することは困難であった。
そこで筆者らは、第一網目に均一な Tetra-PEG ゲル網目を用いて DN ゲルを作成した。
Tetra-PEG ゲルは膨潤度が小さく DN ゲルの第一網目に適さないため、ゲル内部に合成した直鎖電解質高分子による浸透圧効果を利用する分子ステント法を用いて Tetra-PEG ゲルを膨潤させた。この内部に第二網目を合成することで、均一な第一網目構造を持ち強靭で降伏挙動を示す DN ゲルを作成した(St-Tetra-PEG DN ゲル) 1)。この DN ゲルの変形挙動や内部破壊過程を一軸引張試験によって評価した。
【結果と考察】Figure 1 に、St-Tetra-PEG DN ゲルの応力-変形率曲線を示す。分子ステントに用いた AMPSのモノマー濃度 CSt が大きいほど Tetra-PEG ゲルの膨潤率 α[mm/mm]は大きく、CSt=0.35 M 以
上では従来の DN ゲルと同様に降伏挙動を示した。この降伏点と膨潤率の関係を調べたとこ
ろ、降伏応力 σy と降伏変形率 λy はそれぞれ 1/α2, 1/α に比例した。前者は降伏応力が引張方向に垂直方向の高分子鎖の面密度に比例することを示す。また後者は、膨潤前から降伏までの引張方向の高分子鎖の伸長率が膨潤度によらず一定であることを示す。この伸長率は 11 で、調製条件より期待される高分子鎖の両末端間距離と伸び切り鎖長の比 7.1 に近い値を示した。
これは、降伏点までに第一網目の高分子鎖が十分伸びていることを示唆する結果である。

1

Cst=0.20M (α=1.05)
0.35M (α=1.81)
0.50M (α=2.32)
0.70M (α=3.00)
1.00M (α=4.10)
2.00M (α=5.07)

サイクル引張試験では、従来の DN ゲルとにおいて従来の DN ゲルに比べ初期弾性率の低下やヒステリシスが共に小さかった。これにより、第一網目の均一性が高いほど内部破壊が少ないことが示唆された。当日はこれら降伏応力や降伏変形率、ヒステリシスエネ

0.8Nominal Stress σ[MPa]同様にヒステリシス挙動を示したが、降伏前

0.6
0.4
0.2

ルギーなどの物理的意味を定量的に議論する。
【参考文献】

0
1

3
5
Deformation Ratio λ [mm/mm](1) T.Nakajma, Y. Fukuda, T. Kurokawa, T.

Figure 1 The uniaxial tensile stress -deformation

Sakai, UI. Chung, J. P. Gong, ACS macro lett.,

ratio curves of the St-Tetra-PEG/PAAm DN gels

2013, 2, 518

for various stent monomer concentration CSt.

Copyright (C) 2013 ソフトマター研究会 All Rights Reserved.

12月15日
10:00-10:15

ソフトな分子とハードな量子ドットのナノ複合体が実現する
超効率的電子移動（京大院理）金賢得【はじめに】
光励起によって半導体量子ドット内部に生成したキャリア（電子または正孔）の非平衡ダイ
ナミクスはその興味深い物理的・化学的特性から大きな注目をあびている。その一方で、量子ドット内部に生成したキャリアをどのように効率的に量子ドット外部に取り出すかは、ほとんど手つかずの未解決問題であり、そこにまた多くの新たな非平衡現象が存在する可能性がある。
我々は、量子ドット単体だけでなく、ハードマターである量子ドットにソフトマターである
蛍光色素を修飾したソフト―ハードマター複合ナノ構造体において、光励起後の量子ドットから蛍光分子への電子移動を実時間で第一原理的に数値実験できる新しい手法を開発した。
この手法によって、量子ドット内部の Auger 相互作用（電子－正孔の強い相互作用）を考慮
した蛍光分子へ電子移動を実時間で追うことが可能となった。
【結果と考察】我々は、開発した手法を用いて、右図のようなハードマターであるCdSe量子ドットにソフトマターであるメチレンブルー分子が付着した「ソフト―ハードマター複合体」において量子ドットからメチレンブルーへの電子移動を調べた。
計算の結果、マーカス理論における電子移動速度の反転（減少）領域が存在しない新しいタイプの超効率的電子移動が発生することを発見した。
また、量子ドット中に電子とともに生成した正孔が電子移動に伴った余剰エネルギーを受け取り励起することを数値実験的に実証した。我々は、このような電子移動に伴う電子と正孔のエネルギー授受は、量子ドット内に閉じ込められた電子と正孔の強いAuger相互作用によって発生すると結論付けた。
この電子移動速度に減少領域のない電子移動は、従来のマーカス理論に従わない新しいタ
イプの電子移動という意味で物理的に重要であるとともに、ナノ素子界面・表面における効率的な電子移動を目指す実用面からも価値ある結果である。
【参考文献】(1) Haiming Zhu, Ye Yang, Kim Hyeon-Deuk, Marco Califano, Nianhui Song, Youwei
Wang, Wengqing Zhang, Oleg V. Prezhdo, Tianquan Lian, Nano Letters, in press (2014).
(2) Kim Hyeon-Deuk and Oleg V. Prezhdo, submitted.

Copyright (C) 2013 ソフトマター研究会 All Rights Reserved.

12月15日
10:15-10:30

ガラス転移における平均場シナリオの理論的検証(筑波大数理, 名古屋大理 A )池田 晴國, 宮崎 州正 A

【はじめに】
結晶化の機会を逃したまま融点以下まで冷却された液体を過冷却液体と呼ぶ。過冷却液体をさらに冷却していくと緩和時間が劇的に増大し、ついには観測時間を越えてしまう。この現象はガラス転移と呼ばれ、分子性液体だけではなく金属液体や高分子溶液等、多くの物質を急冷した際に見られる普遍的な現象である[1]。ガラス自体は金属についで我々に身近な物質であるにも関わらず、ガラス転移の微視的起源は未だ明らかにされていない。
ガラス転移は通常の相転移とは異なった現象である。結晶化や強磁性転移等の、熱力学的相転移における緩和時間は、静的な相関長と共に増大して行く。この為、多数の粒子が強調的に動くことが遅い緩和を引き起こす原因であると信じられている。これに対してガラス転移では、分子の配置は液体同様乱雑なままに、緩和時間だけが増大して行く。ガラス転移の背後には如何なる機構が隠されているのだろうか?
我々はガラス転移の背後に隠れた熱力学的機構があるという立場に立って研究を進めている。この立場に立った際に参考になるのはスピングラスの平均場模型である。スピングラスでは過冷却液体と同様に配位が乱れたまま動力学が遅くなる。そして、少なくとも平均場近似のレベルでは遅い
緩和とそれを引き起こす熱力学的機構の関係が完全に明らかにされている。
ガラス転移がある種のスピングラス模型が示す転移と同等のユニバーサルクラスに属するという仮定のもとに、ガラス転移の平均場シナリオが提唱された[2]。このシナリオでは過冷却液体の温度を減少させていくと、動的転移点と呼ばれる温度で自由エネルギーに指数関数的な量の極小が表われ初めることを主張する。極小の数が指数関数的なので、全ての極小を足し合わせれば一つの極小にある状態よりもエントロピーが高い状態を作ることが出来る。この為熱力学的には全ての極小について重ね合わせた乱雑な状態がより安定となる。一方で、動力学はこの極小の影響を受け、極小の間を熱活性化過程によって遷移していくことによって遅い緩和が実現されると考えられる。
上のシナリオは、運動論と液体論という独立に発展して来た二つの理論が、ガラス転移の研究においては互いに整合する可能性を示唆するという点で興味深い。ガラス転移の研究を通して二つの理論の利点を相互に参照し、より精度の高い汎用的な理論を構築できることが期待されるからである。
【結果と考察】上で説明したシナリオは、ある種のスピングラスの平均場模型においては厳密に成立している。しかし当然ながら、同じシナリオが実際のガラス転移において成立する保証は無い。そこでこれを理論的に検証することが本研究の目的である。
我々は、自由エネルギーに多数の極小が表われ初める温度と、動力学が遅くなり初める温度を場の理論を用いた摂動計算によって算出した[3]。その結果、近似の最低次の範囲では、二つの結果が整合していることが明らかになった。この結果はガラス転移における平均場描像の成立を示
唆する強力な状況証拠になると考えられる。
【参考文献】[1][2][3]

P. G. Dobenedetti, et al, Nature materials, 410, 259(2001)
T. R. Kirkpatrick, et al, Phys.Rev.A, 40, 1045(1989)
G. F. Mazenko, Phys.Rev.E 83, 041125(2011)

Copyright (C) 2013 ソフトマター研究会 All Rights Reserved.

12月15日
10:45-11:00

弱電解質拡散場中の高分子電解質ドメインのダイナミクス
名古屋市立大学

奥薗

透、豊玉彰子、山中淳平【はじめに】 電解質濃度の勾配がある系におけるコロイド粒子の運動については、既に多くの研究があるが、それらのほとんどは強電解質の系が仮定されている。また、貧溶媒系の高分子電解質のダイナミクスに関する研究も多くはない。このような電解質の系は複雑であり、古くから研究されているにもかかわらず、理論的にも、実験的にも不明な点が多い。
本研究では、弱電解質の濃度勾配下における高分子電解質の相分離ドメインのダイナミクスを記述するモデルを構築し、数値シミュレーションを行う。モデルは、高分子濃度、電解質濃度、および正、負のイオン濃度に関する時間発展方程式によって構成される。自由エネルギーには、フローリー・ハギンス型の自由エネルギー密度に加え、高分子と電解質の相互作用を取り入れる。これは、小貫ら（1,2）によって提案されている自由エネルギーとほぼ等価である。ただし、弱電解質の効果は時間発展方程式の反応項として取り入れる。さらに、静電相互作用および流体力学的相互作用も考慮する。
【結果と考察】 相分離ドメインが存在するような状況において、電解質のリザーバーに接触した系を設定し、2 次元の系で数値シミュレーションを行った結果、電解質の一方向拡散に伴う相分離
ドメインの移動が認められた。移動の方向は、相分離ダイナミクスによる拡散的な効果と流体力学的な効果の両方が関与しており、それらの効果のバランスによって決まっていることがわかった。

【参考文献】(1) A. Onuki and R. Okamoto, J. Phys. Chem. B 113, 3988 (2009).
(2) T. Araki and A. Onuki, J. Phys. Condens. Matter 21, 424116 (2009).

Copyright (C) 2013 ソフトマター研究会 All Rights Reserved.

12月15日
11:00-11:15

クーロン系ソフトマターシミュレータによる電気浸透流の解析（豊田中研・京大触媒電池）吉田広顕，金城友之，鷲津仁志【はじめに】
マイクロ・ナノスケールの固液界面における高効率のエネルギー輸送・変換のため，長距離クーロ
ン力が活用される事例は多い．電気浸透流 (EO) ポンプは代表的な応用例の一つであり，表面に電荷を有するマイクロチャネルを電解質溶液で満たすと固液界面に電気二重層が生じ，電場の印
可によりイオンを動かすと同時に溶媒を引きずり流動を生じさせるという機構である．昨年，谷口らにより提案された EO ポンプ(1)は，ナノビーズ修飾したポリカーボネートのトラックエッチド (TE) 膜に交流電場を印加するというソフトマテリアルの特色を活かしたデバイスであり，複雑な MEMS 的工作を必要とせず，直流高圧電場による水の電気分解を抑えられる，このような EO 現象の解析のためには，サブミクロン孔径のチャネルにおける電解質溶液の電場下における流動を扱わなければならない．また，カルボン酸基により負に帯電していると思われるポリカーボネート表面の状態予測や，界面すべりなどの分子レベルでの流動性の議論のためには，分子〜メソスケールにおけるシミュレーシ
ョン手法の開発が必要である．本報告では，筆者らによるマルチスケールシミュレータによる EO ポ
ンプの解析事例を紹介する．
【シミュレーション】
まず，イオン輸送，溶媒の流れ，外部電場を記述する
方程式が互いに連成した連続体モデルを解析するための計算手法を作成した(2)．通常，イオン分布に Boltzmann
分布を仮定した簡略化したモデル方程式が用いられるが，非平衡過程や，バルク状態を定義できないマイクロある
いはナノチャネル系においてはこの仮定が適切でないケースがある．本研究では，イオン輸送過程を
Nernst-Planck 方程式に基づいて直接数値解析し，マイ
クロチャネル内の電気浸透流を再現した (図 1)．
さらに，分子レベルにおいて電気二重層の生成およびイオンによる溶媒の引きずり効果を解析するため，Brown
粒子を Langevin 動力学により扱い，溶媒の流れを格子
Boltzmann 法により解析する手法を提案した (3) ．さらに，
イオン周囲の溶媒和を散逸粒子動力学で記述する手法(4)など，各スケールにおける解析手法を提案している．

【参考文献】(1) 谷口学，奥村泰志，樋口博紀，菊池裕嗣，ソフト
マター研究会予稿集, 2012, P-27.

図 1 イオン流動シミュレータによる電気浸透流の解析例(a)静電ポテンシャル(b) 正味の電荷(c)圧力(d) 流れ場．壁面
での表面電荷密度一定、左右に外部電位差を付与．

(2) Yoshida, H.; Kinjo, T.; Washizu, H., in Proc. 3rd Euro. Conf. Microuidics, 231 (2012).
(3) Yoshida, H.; Kinjo, T.; Washizu, H., in Proc. 3rd Intl. Conf. Mol. Sim. (2013).
(4) Kinjo, T.; Yoshida, H.; Washizu, H., J. Phys. Soc. Jpn. suppl. (in print).

Copyright (C) 2013 ソフトマター研究会 All Rights Reserved.

12月15日
11:15-11:30

細胞膜陥入構造“カベオラ”の力学（理化学研究所）立川正志【はじめに】
カベオラは特定の細胞の細胞膜に多数存在する 100nm 程度の陥入構造で，Caveolin1 タンパ
ク質により裏打ちされていること，ラフトと同様な脂質組成を持つこと，多数のレセプターが集積していること等知られているが，カベオラそのものの機能的意義はこれまで明らかにされてこなかった．近年，低浸透圧下においてカベオラが解消(平坦化)することが発見され，この効果が細胞膜面積ひいては細胞体積を一時的に増加させる役割，つまり細胞膜のバッフ
ァの役割を果たすという説が提出された(1)．その後，末次志郎准教授(東大)らの実験により，低浸透圧下では解消するカベオラも存在する一方，無視できない数のカベオラがエンドサイ
トーシスにより細胞内に取り込まれることが見出された．
【結果と考察】本研究では，この同じ刺激に対して，なぜカベオラは２種類の全く相反する応答を示すかを考えるために，低浸透圧刺激 onset でカベオラにかかる力を Time dependent Ginzburg Landau
方程式を用いて計算した．その結果，カベオラに収縮力がかかるか拡張力がかかるかはカベ
オラの形態(陥入構造の深さと断面半径の比)に依存することを見出した．このことは，カベオラの刺激応答の２値性は，その形態に依存して決まる力が生み出していることを示唆している．実際，さらなる実験結果から，低浸透圧下で収縮したカベオラと拡張され広がったカ
ベオラが観察されており，この力の２値性の存在を裏付けている．
また，Gompper の方法を用いたポリゴン膜シミュレーション(2)に，膜を曲げる作用を持
つ足場タンパク質を導入してカベオラ様陥入構造を再現し，浸透圧変化に対する形態変化の応答を数値的に調べた．結果，集積している足場タンパク質の量に依存して応答が異なり，量が少ない場合は浸透圧低下とともに構造が解消するが，多い場合は陥入構造のくびれが進行し，小胞形成へと進む事が見出された．陥入構造の深さは足場タンパク質の量と正の相関を示すため，この結果は上の計算とも整合している．
これまで，カベオラの理論モデルは脂質組成に注目したドメインモデルのみが議論されており(3)，低浸透圧に対してはカベオラ解消のみが示されていたが，本研究で，足場タンパ
ク質による形態形成に注目することにより，カベオラの多面的な形態変化を記述することに成功した．

【参考文献】(1) B. Sinha, et al., Cell 144, 402 (2011).
(2) G. Gompper and D.M. Kroll, J. Phys. Condens. Matter 9, 8795 (1997).
(3) P. Sens and M. S. Turner, Phys. Rev. E, 031918 (2006).

Copyright (C) 2013 ソフトマター研究会 All Rights Reserved.

12月15日
11:30-11:45

アクティブソフトマターの剪断流中でのダイナミクス（ 1 京 都 大 理 、 2 D ü s s e l d o r f 	
 	
 UU n i v e r s i t ä t 、 3 東 京 大 物 性 研 、 4 東 京 大 理 、 	 5
U n i v e r s i t y 	 	
 oo f 	 	 EE d i n b u r g h 、 6 お 茶 大 ソ フ ト マ タ ー セ ン タ ー ） 	
 多 羅 間 充 輔 1 2 3 4 , 	 	
 AA . M . 	
 	
 MM e n z e l 2 1 , 	 	
 BB . 	 	 tt e n 	 	
 HH a g e n 2 , 	
 	
 	
 	
 
R . 	
 	
 WW i t t k o w s k i 5 , 	 	 太太 田 隆 夫 1 4 6 , 	 	
 HH . 	 	 LL ö w e n 2 	 【はじめに】

近 年 、ア ク テ ィ ブ マ タ ー と よ ば れ る 自 発 的 に 運 動 す る 粒 子 の 運 動 が 盛 ん に 研 究
さ れ て い る 。そ の 例 は ソ フ ト マ タ ー の 中 に 数 多 く み ら れ 、細 胞 や 微 生 物 な ど の
生 物 か ら 、自 発 運 動 す る 液 滴 や ヤ ヌ ス 粒 子 と 呼 ば れ る コ ロ イ ド な ど の 非 生 物 まで 多 岐 に 渡 る 。そ れ ら の 静 的 な 環 境 下 で の ダ イ ナ ミ ク ス に つ い て は 、こ れ ま でに も 理 論 研 究 が 進 め ら れ て き て お り 、自 発 的 な 並 進 運 動 と 変 形 と の 相 互 作 用 によ り 円 運 動 が 自 発 的 に 起 こ る な ど １ 粒 子 で あ っ て も 非 常 に 興 味 深 い 。し か し 現
実 に は 多 く の 場 合 、こ れ ら の 粒 子 は 化 学 走 性 や 光 走 性 、重 力 な ど の 外 場 の あ る
環境下で運動している。とくに重要な問題の一つにマイクロチューブ中の
Couette 流や Poiseuille 流などの流れの中での運動がある。そのような流れの 中 で の ダ イ ナ ミ ク ス は 、自 発 運 動 す る 剛 体 粒 子 に 対 し て は す で に 研 究 さ れ て
き た が 、変 形 を 伴 う 柔 ら か い ア ク テ ィ ブ マ タ ー に つ い て は 、わ れ わ れ が 知 る 限
り い ま ま で 調 べ ら れ て こ な か っ た 。わ れ わ れ は 、流 れ の 中 の ア ク テ ィ ブ ソ フ ト
マ タ ー の ダ イ ナ ミ ク ス を 理 論 的 に 理 解 す る た め 、そ の 運 動 を 記 述 す る モ デ ル を考案し、その解析を行った。

【結果と考察】ア ク テ ィ ブ ソ フ ト マ タ ー の 例 は 多 岐 に 渡 り 、そ の 運 動 機 構 な ど も 様 々 で あ る た
め 、出 来 る 限 り 系 の 詳 細 に は 依 ら な い 一 般 的 な 理 解 を 試 み た 。そ の た め 変 数 を
ベ ク ト ル お よ び テ ン ソ ル で 表 現 し 、対 称 性 の 議 論 の み か ら 意 味 の あ る 運 動 方 程式 を 導 出 し た 。自 発 的 な 並 進 運 動 の ほ か に 、粒 子 の 自 発 的 な 回 転 運 動 も 考 慮 した 。線 形 剪 断 流 を 一 つ の 例 と し て 、異 な る 自 発 速 度 お よ び 剪 断 の 強 さ に 対 し て方程式を数値的に解くことにより、様々な運動を得た。

【参考文献】(1) M.	
 Tarama,	
 A.M.	
 Menzel,	
 B.	
 ten	
 Hagen,	
 R.	
 Wittkowski,	
 T.	
 Ohta,	
 and	 H.	
 Löwen,	
 J.	
 Chem.	
 Phys.	
 139,	
 104906	
 (2013).	
 Copyright (C) 2013 ソフトマター研究会 All Rights Reserved.

12月15日
11:45-12:00

対称的な性質を持つ壁の間に挟まれた 2 元混合系の相分離：バルク
臨界点付近の Capillary transition（ 京都大学理学研究科）藪中	  俊介、岡本	  隆一、小貫	  明【はじめに】狭い空間に閉じ込められた流体の相挙動は、様々な系において調べられていて、その相挙動は、壁の形状、壁と分子の間の相互作用に大きく影響される事が分かっている。実際、壁と分子の相互作用により、多くの場合、気液系の液相や、二元混合系のうちの一成分が壁に選択的に吸着される。粒子溜めと平衡を保った状況では、壁に挟まれた流体では、外部の粒子溜めの化学ポテンシャルを変化させる（粒子
溜めの中の壁の好まない成分の組成を増やす）事で、壁の好まない成分が壁の間に侵入 す る Capillary transition と い う 一 次 相 転 移 が 起 こ る 事 が 知 ら れ て い る 。本 研 究 では 、 Capillary transition 線 付 近 で の 臨 界 流 体 の 相 分 離 を 研 究 し た 。
【結果と考察】本研究では、臨界揺らぎによる繰り込みの効果を岡本、小貫によって最 近 提 案 さ れ た Local Functional theory を 用 い て 取 り 込 ん で 計 算 を 行 っ た 。 ま ず 、
平均組成を固定した上で 1 次元での平衡解を計算し、平均組成と化学ポテンシャルの関 係 を 見 る と 、 Capillary transition の 臨 界 温 度 Tc ca 以 下 で は Van der W aals loopが 現 れ た 。 次 に 三 次 元 計 算 で 、 Van der W aals loop の 傾 き が 負 で あ る 、 ス ピ ノ ダ ル
線の内側の領域に系を温度クエンチした。すると、初期には、壁に垂直な方向（横方€向）に関して系はほぼ均一であるが、ある程度時間が経つと、横方向の揺らぎに関して、初期状態が不安定化し、相分離が起こり、壁の好まない成分の扁平なドメインが壁で挟まれた領域の中央に現れた。さらに、この相分離が流体効果によって加速されている事も確認した。

図１：1 次元解から得られた、化学ポテンシャルと縦方向の組成の平均の間の関係。点線部分は、横方向の揺らぎに対し不安定である

図 ２：壁 に 挟 ま れ た 空 間 で の 相分離過程。色は組成場を表し、矢印が流れ場の大きさを示す。

図 ３：相 分 離 状 態 の 各 ド メ イ ン での 、中 央 に お け る 組 成 、及 び 縦 方
向 の 平 均 の 組 成 。一 次 元 の 平 衡 解
か ら 計 算 さ れ た 共 存 線（ 赤 、青 の実線）の上に乗っている。

【参考文献】(1) R. Okamoto and A. Onuki, J. Chem. Phys. 136, 114704 (2012)(2) S. Yabunaka, R. Okamoto, and A. Onuki, Phys. Rev. E 87 032405 (2013).

Copyright (C) 2013 ソフトマター研究会 All Rights Reserved.