光散乱(DLS, SLS, ゼータ電位)
------------------------------------------------------------------------------ ■ 光散乱 光散乱現象を利用して、液体、溶液、縣濁液のさまざまな構造や物性を調べることができる。 低分子系にも適用でき、密度や濃度の揺らぎに関する情報(例えば、等温圧縮率や浸透圧縮率)が得られるが、溶液中のコロイドや高分子は低分子に比べて光を散乱する能力が高いので、光散乱法は高分子・コロイド溶液の研究に適した測定法である。また、光の波長が高分子の大きさ(10 - 数百 mm)を測定するのに適していることも、高分子・コロイド研究に光散乱法がよく用いられる理由である。 光散乱法には、静的光散乱と動的光散乱の2種類の方法がある。 ------------------------------------------------------------------------------ ◆ 静的光散乱(Static Light Scattering: SLS)
◇ 分子量 ・静的光散乱法は、高分子・コロイド試料の分子量測定の最も標準的な方法である。サイズ排除クロマトグラフィーなどのように分子量が既知の標準試料を必要とせず、粘度法のように粘度式が文献に与えられている必要もない絶対測定法であり、新規な高分子・コロイドの分子量(モル質量)決定が可能である。 ・質量分析やNMR法は、かなり高分子量の試料の絶対分子量測定が可能になってきたが、分子量分布のある試料では分子量によってイオン化される度合いが異なっていたり、末端基が明確に定まってない試料などの場合には正確な(平均)分子量が得られない。また、10 6 を超える高分子量試料の絶対分子量測定には静的光散乱法以外の適当な測定法がないのが現状である。 ・光散乱法によって分子量が測定できるのは、試料の試料濃度(単位:g/cm 3 )が一定の溶液では、分子量が高いほど散乱光強度が強くなる性質に起因している。 ・1個の高分子あるいはコロイド粒子から散乱される光の電場の振幅は1粒子を構成する電子数(すなわち分子量)に比例する。電場の2乗に比例する光の強度は分子量の2乗に比例する。散乱体積内の全粒子からの散乱光強度は粒子数に比例する。質量濃度が一定の溶液では粒子の数濃度は分子量に反比例して減少する。以上の結果、散乱体積内の全粒子からの散乱光強度は分子量に比例することになる。 ・以上の議論は同一あるいは異なる粒子に属する別の構成単位から散乱された光が干渉することを無視している。これらの粒子内および粒子間の干渉効果は、それぞれ散乱角あるいは粒子濃度をゼロに外挿することによって消すことができる。 ・溶液全体からの散乱光強度は粒子と溶媒からの散乱光強度の和として与えられるので、粒子成分から散乱された光に対するレイリー比は、溶液と溶媒のレイリー比の差(過剰レイリー比と呼ぶ)ΔRθによって与えられる。これより次式が成立する。 lim θ,c→0 ΔRθ/c = K*M ここで、cは粒子の質量濃度、Mは分子量、Kは光学定数と呼ばれる比例係数で、次式で定義される。 K ≡ ((4* π 2 *n 2 )/(NA*λ 4 ))*(∂n/∂c) 2 ここで、NAはアボガドロ定数、また∂n/∂cは溶液の屈折率の濃度増分で、構成単位の光散乱能と関係した量である。この∂n/∂cは、溶液と溶媒との屈折率差を精密に測定できる示差屈折率計を用いて実験的に求めることで、Kを実験的に決定し、分子量を求めることができる。ただし、試料と溶媒の屈折率が等しく、∂n/∂cがゼロの特別の場合には、試料に光散乱能がなく、光散乱実験は行えない。 ・試料に分子量分布がある場合には、Mは重量平均分子量で置き換えられる。複数種類の構成単位からなるコポリマー試料については。組成分布があると、一般には正しい分子量が求まらない(タンパク質のように組成分布のない場合にはその限りではない)。 ・構成単位の過剰屈折率が近い場合には、それほどの分子量の誤差は生じないが、厳密には屈折率の異なる最低3種類の溶媒を用いて光散乱(および示差屈折率)測定を行い、それぞれで得られた見かけの分子量から真の分子量を評価する方法がある。 ・混合溶媒を用いて光散乱測定をする場合には、選択的溶媒和と呼ばれる現象が起こり、粒子の近傍の溶媒組成は平均の組成とは異なっている可能性がある。その場合には、通常の方法で測定した∂n/∂cは構成単位の光散乱能を正確に表しておらず、計算式から正しい分子量を見積もれない。この困難は、透析と呼ばれる手法を用いて解消することができる。 ・透析とは溶媒成分のみが透過でき、粒子は透過できない半透膜を隔てて、溶液と大量の溶媒を接触させ、平衡化させる操作である。選択的溶媒和が起こっていれば、溶液中の粒子との親和性がより高い溶媒成分が流れ込み、溶液側の溶媒組成が変化する。この透析した溶液と溶媒を使って測定した∂n/∂cを計算式に代入して計算したKを使うと、正しい分子量が得られることが熱力学的に証明されている。 ・厳密には光散乱測定もこの透析した溶液について行うべきであるが、一般に、光散乱測定結果は透析前後でほとんど変化しないので、未透析の溶液についてのΔRθで代用しても有意の誤差は生じない。 ・光散乱で分子量を求める際に生じる測定誤差は、散乱角度と試料濃度の外挿操作に起因することが多い。外挿のしかたについてはいろいろな方法が提案されている。高分子電解質の水溶液系は、長距離の静電相互作用によって粒子間の干渉効果を消すことが難しい。通常は低分子塩を0.1 M程度添加して静電相互作用を遮蔽して分子量測定が行われる。
◇ 分子量測定の注意点 ・分子量測定は、散乱光量の絶対値を計測するため、分散状態は安定、分子量分布も狭い状態のサンプルが必要 特定の時間に大きなピークが集中:フィルター処理を実施する 大きくうねるようなバックグラウンドがある:温度平衡時間を設ける ・Debye plotは希釈系列から、濃度0のY軸切片を算出する方法、そのため、高濃度からの直線近似ではなく、できる限り希薄な状態での直線近似が好ましい
◇ 粒子の形状・コンホメーション ・粒子サイズが10 nm程度以上になると、同一粒子内の異なる部位から散乱された光が干渉を起こす。この効果を利用して、散乱粒子のサイズや形状を調べることができる。 ・位相のずれは、散乱光源である散乱要素間の相対位置によって決まるので、散乱体内の構成要素の分布状態で粒子内干渉効果は規程されている。 ・異方的な粒子の場合は、溶液中でいろいろな方向を向いているので、この干渉効果は粒子配向に関して平均化する必要がある。具体的な計算方法は省略するが、さまざまな形状(構成要素の分布)の粒子について粒子内干渉効果の計算が行われている。結果は式で定義される粒子散乱関数(分子内干渉因子あるいは形状因子とも呼ばれる)P(θ)によって記述される。 lim c→0 ΔR θ /(K*c) = M*P(θ) ここで、試料濃度をゼロに外挿する操作は、粒子間の干渉効果を消すためである。上式のMに[lim θ,c→0 ΔR θ /c = K*M]を代入すると、両辺のKが消失するので、P(θ)を求めるには、∂n/∂cの値は必要ない。 構成要素が均一に分布したさまざまな幾何学的形状を持つ粒子や、高分子のモデル鎖であるガウス鎖やみみず鎖に対するP(θ)の具体的な関数が知られている。これrなお関数型と光散乱実験により得られたP(θ)の実測値とを比較することにより、試料粒子の形状を推定できる。ただし、多分散性の試料の場合には、サイズ分布を考慮して理論と実験との比較をする必要がある。一般に、P(θ)はサイズ分布の影響を強く受けるので注意したい。 ・通常、P(θ)は単調な関数なので、粒子形状を一つの試料に対するP(θ)の曲線から一義的に決定することは困難な場合が多い。また、研究対象の粒子に対する適当なP(θ)の理論式がない場合も多い。そのようなときには、P(θ)そのものではなく、P(θ)対q 2 のプロットの初期勾配より次の式を使って計算される平均2乗回転半径と呼ばれる量を用いて粒子の形状を議論するのが有効である。 P(θ) = 1 - (1/3)**q 2 + Ο(q 4 ) このは粒子の形状とは無関係に次式で定義される量である。 ≡ (1/N0)*Σi=1 N0 Gi 2 > ただし、N0は1粒子中の構成要素数(重合度)、S Gi は粒子の重心からi番目の構成要素までの距離を表す。角かっこは粒子(例えば高分子鎖)の内部自由度に関する平均を表す(剛体の場合には平均操作は必要ない)。この定義により 1/2 は粒子の有効な半径とみなすことができる。 ・同種類で、分子量(すなわちN0)だけが異なっている高分子・コロイド試料について、とMを光散乱より測定し、両者の関係がわかると、試料粒子の形状・コンホメーションに関する知見が得られる。均一密度の球状粒子の場合には、半径RはMの1/3乗に比例するので、はMの2/3乗に比例する。 ・良溶媒中での屈曲性高分子は、分子内排除体積効果のためにガウス鎖としてふるまわずに、はMの1.2乗に比例することが知られている。したがって、対Mの両対数プロットを作り、その傾きを見積もると粒子形状(形態)を知ることができる。 ・ブロックコポリマーのP(θ)あるいはは、各ブロック鎖の(∂n/∂c) 2 の重みがついた平均量として与えられる。特別な場合として、片方のブロック鎖の(∂n/∂c)がゼロの場合には、散乱能のあるブロック鎖だけのP(θ)あるいはが求められる。 ・さまざまな形状の散乱体の粒子散乱関数と回転半径に関する文献 H. Yamakawa, Modern Theroy of Polymer Solutions, Chap. V, Haper & Row (1971) 斉藤信彦、高分子物理学、裳華房 (1967) 倉田道夫、近代工業化学18, 高分子工業化学?, 朝倉書店 (1975) M. B. Huglin (ed.), Light Scattering from Polymer Solutions, Academic Press (1972)
◇ 溶液中での粒子間相互作用 ・光散乱法では、粒子間の干渉効果を利用して、溶液中での粒子間に働く相互作用の知見を得ることができる。粒子間干渉効果を考慮に入れた2成分溶液の過剰レイリー比の式を次に示す。 lim θ→0(K*c)/ΔRθ = (1/(R*T))*(∂Π/∂c) = 1/M + 2*A2*c + 3*A3*c^2 + . ここで、R*Tは気体定数と絶対温度の積、Πは溶液の浸透圧(∂Π/∂cは浸透圧縮率の逆数と呼ばれる)、またA2とA3は第2、第3ビリアル係数である。後ろの等式は、浸透圧のビリアル展開より得られる。散乱角をゼロに外挿するのは、粒子内乾燥効果を消すためである。 ・分子量分布を有する試料の場合には、上式中のMは重量平均量、A2とA3についても分子量分布に関するある平均量となるが、∂Π/∂cとの関係式は厳密には成立しなくなる。 ・溶液中での粒子間には、一般に斥力と引力が働く。良溶媒中では引力が弱く、貧溶媒中では強い。第2ビリアル係数は斥力項と引力項とからなっており、引力の弱い良溶媒中ではA2の値が大きく、貧溶媒中では小さい。 ・高分子溶液の場合にはA2がちょうどゼロとなるときをシータ状態と呼び、一種の理想状態とみなす。見掛け上、高分子間の相互作用が消失した状態で、分子内排除体積効果も同時に消失して、屈曲性高分子鎖はガウス鎖としてふるまう。ただし、このシータ状態でも高次のビリアル係数項まで含んだ∂Π/∂cを濃度の関係として求め、適当な統計力学理論を利用すると、粒子間の斥力と引力を分離評価できる。 Y. Jinbo et al., Macromolecules, 27(1994)6080. R. Koyama, T. Sato, Macromolecules, 35(2002)2235. ・A 2 > 0: 「親溶媒性」、溶媒と溶質の馴染みが良い(分散状態良好) ・A 2 < 0: 「親溶質性」、溶媒と溶質の馴染みが悪い(凝集傾向)※ 「親溶媒性」と「親溶質性」は、媒と質の一字違いなので見間違えないように注意。
◇ 高次構造 ・低分子液体に対してX線散乱実験を行うと、結晶ほど鋭くはないが数本の回折ピークが観測される。これは、液体中で分子は主として斥力相互作用により近距離の規則的な配列をとっている証拠である。濃厚な縣濁液中の球状ラテックスでも、スケールが違うが、同様な「液体構造」をとっており、これは光散乱によって検出できる。 ・低イオン強度の水溶液中で、電荷を有する球状や棒状粒子は、長距離の強い静電相互作用によって、秩序性の高い規則的配置をとる。 ・有限濃度での剛体粒子系では、粒子間干渉効果のために次式で表される。 ΔRθ/(K*c) = M*P(θ)*S(q) ここで、S(q)は静的構造因子と呼ばれる系内の秩序構造を表す関数である。粒子形状が濃度に依存しないときには、M*P(θ)は濃度外挿より求められて、実際のΔRθ/(K*c)とM*P(θ)よりS(q)を求めることができる。 ・剛耐球や剛体棒の系に対するS(q)の計算方法は下記文献を参照。 E. E. Maier et al., Macromolecules, 25 (1992) 1125. J. P. Hansen, I. R. McDonald, Theory of Simple Liquids, Academic Press (1986). ・高分子容器では濃度が高くなると高分子鎖どうしは相互貫入できるので球状粒子のような規則構造はとらない。高分子の濃厚溶液に対する過剰レイリー比は通常、次式で解析される。 (K*c)/ΔRθ = (1/(R*T))*(∂Π/∂c)*[1 + (1/3)*ξ 2 *q 2 + Ο(q 4 )] ここで、ξは相関長と呼ばれる量で、高分子濃度溶液の性質を議論するときの基本量としてしばしば利用される。
◇ 粒子内・粒子間の配向相関(異方性散乱) ・検出器の前に水平偏向のみを透過するように偏光板(検光子)を挿入すると、異方性散乱成分ΔRθ,HVだけを検出できる。このΔRθ,HVは散乱要素間の配向相関と関係する量である。無限希釈状態では粒子内の要素間の配向相関のみが反映され、 limc,θ→0ΔRθ,HV/(K*c)より粒子の形態に関する知見が得られる。 ・例えば高分子では、limc,θ→0ΔRθ,HV/(K*c)は棒状形態では分子量に比例し、屈曲性鎖では分子量に依存しないので、その分子量依存性を調べれば高分子鎖の剛直性を評価できる。 ・液晶状態にある系や液晶相転移近傍の等方相系では、分子間に長距離の配向相関があるので、異方性光散乱によって、液晶の弾性率や配向揺らぎ、熱力学的安定性などについての知見が得られる。 ・スピノーダル点(配向揺らぎの発散点)を評価できる
◇ SLS測定の注意点 ・粒子径の結果からSLSが測定できるかどうかを判断する場合 多分散指数(PdI)が0.2以下が好ましい(小さいほどよい) ・測定時の濃度は、できるだけ低濃度の領域で測定する 高濃度領域では分子間相互作用によってDebye plotの傾きが変わるため ※ Debye plot: 角度固定で、c(濃度)に対して(K*c)/Rθをプロットしたもの ・分子量分布が広い(≒凝集体含む)のは苦手 大きな分子の散乱光が支配的になってしまうため ・解析には試料の∂n/∂c値が必要 ・測定する試料の濃度情報(mg/mL)も必要 ・粒子を混ぜたときのピーク比 散乱強度基準:体積比率が明確に分けられない 体積基準:散乱強度基準に比べ体積比率が確認しやすい
◇ SLS解析の注意点 ・4点以上で測定してみる 1点が直線から外れるデータであった場合、残りの3点で直線を引くことができるため ・通常の角セルでは問題ないが、特注の小容量セルを使用する場合には光路長の幅に注意する セルのエッジが散乱光を不安定にし、測定結果が安定に得られなくなる ・キュムラントフィットと粒子径分布フィットが共に一致しているかを確認する キュムラントフィットが一致していないとキュムラント径(Z-Average)は信頼性が低くなる ・粒子濃度を高くしたとき 粒子濃度でピーク形状が変わらない場合:粒度補正をするとピークが一致する。粒子間相互作用の影響はない。 粒子濃度でピーク形状が変わる場合:粒子間相互作用の影響がある ・粒子間拡散は濃度に応じてピークの形状が異なる ・制限拡散(粒子の拡散が遅くなる)は濃度が変化してもピークの形状は変化しない ------------------------------------------------------------------------------ ◆ 動的光散乱(Dynamic Light Scattering: DLC)
試料 直接測定する主な現象 得られる主な知見 微粒子分散液 並進拡散、回転拡散 流体力学的半径 形状異方性 高分子溶液(希薄) 重心運動、分子内運動 鎖の広がり方 鎖と溶媒の相互作用 高分子溶液(準希薄) 濃度揺らぎ 擬似網目のサイズ 絡み合いの緩和 高分子溶液(濃厚)、融液 分子内運動、分子間協同運動 鎖の運動性 ガラス転移過程 ゲル 網目の協同拡散、クラスターの運動 網目のサイズ ゾル-ゲル転移過程 ミセル 並進拡散、回転拡散、濃度揺らぎ、 孤立ミセルの形状 巨大ミセルの構造 液晶 分子配向秩序の揺らぎ、濃度ゆらぎ 弾性定数 相転移近傍での臨界緩和◇ 概要 ・動的光散乱は試料にレーザー光を入射し、試料から散乱された光を測定して、その時間変化を分析する。散乱光の時間変化は、試料内部で起こっているブラウン運動の速さに対応している。 ・溶媒中の粒子サイズが大きいとブラウン運動による拡散運動が穏やかなので、散乱光強度の揺らぐ速さが遅くなる。このとき動的光散乱では、散乱光の時間変化あら相関関数を算出する。相関関数から求められる緩和時間から、ブラウン運動の速さが見積もられる。 ・粒子が大きい:運動が遅いので散乱光強度の揺らぐ速さが遅く緩和時間が長くなり ・粒子が小さい:運動が速いので散乱光強度の揺らぐ速さが速く緩和時間が短くなる ・一般的な動的光散乱の装置で測定できる相関時間の範囲は10 -6 から10 2 秒程度である。 ・粒子の拡散係数を反映した散乱光のゆらぎを検出し、粒子径を測定する方法。拡散係数を求める方法は2つあり、光子相関法と周波数解析法に分類される。一般的に光子相関法を採用しているメーカーが多い。 ・光子相関法:相関計(コリレーター)を用いて自己相関関数を導き、この自己相関関数から、平均粒子径および多分散指数を得る方法 ・周波数解析法:散乱光の周波数変化をフーリエ変換し、そこから平均粒子径および多分散指数を得る方法
◇ 装置 ・可視レーザー光(出力は数10 mWから1 W程度)を試料に入射し、試料から散乱される散乱光を検出する ・入射光の進行方向と散乱光の進行方向の成す角が散乱角θである ・散乱装置の設定によって散乱ベクトルqの大きさが次のように定まる |q| = ((4* π *n)/λ)*sin(θ/2) ここで、nは試料の屈折率、λはレーザー光の波長である。 ・散乱光は微弱光の検出に優れた光電子増倍管やアバランシェフォトダイオードによって検出する。 ・微弱な光は、光子パルスとして測定されるため、プリアンプ/ディスクリミネーターによってパルスの増幅と整形/ノイズ除去を行なった後、デジタル通信としてコリレーター(相関計)に入力する ・コリレーターは光子パルス信号から、散乱光の自己相関関数を計算する ・計算手順の一例を下に示す 1. 入力された光パルスを一定のゲート時間ts間隔でカウントする 2. 散乱光強度の時間変化n(i)をデジタルデータとしてメモリに取り込む 3. デジタル計算によって光パルスの時間自己相関関数C(j*ts)を計算する C(j*ts) = = (1/N) * Σ ここでNはデータn(i)の総数であり、は時間平均を表す。 τ = j*tsとすると、相関関数は相関時間τの関数C(τ)となる。τ = 0のときには、相関関数は散乱光強度n(i)の2乗の平均C(0) = になり、τ = ∞のときには、相関関数は散乱光強度n(i)の平均の2乗 C(∞) = 2 になる。一般に、光子パルスの信号列n(i)はガウス統計に従うので、 = 2* 2 、すなわち、C(0) = 2*C(∞)という関係が成り立つ。この関係があるので、C(τ)は次式によって規格化される。 c(τ) = C(τ) / C(∞)
◇ 相関関数 動的光散乱で直接測定されるのは散乱光強度の揺らぎである。散乱光強度は散乱電場の2乗なので、散乱光強度の自己相関関数c(τ)は通常、散乱電場の2次の(2乗の)相関関数という意味で、g (2) (q,τ)と書かれる。すなわち、c(τ) = g (2) (q,τ)である。 一方、散乱電場の(1次の)自己相関関数g (1) (q,τ)はg (2) (q,τ)と次の関係がある。 g (1) (q,τ) = √ = √ 動的光散乱の理論によれば、試料内部の複素誘電率の変化量をδε(q,t)とすると g (1) (q,τ) ∝ g (2) (q,τ) ∝ という関係がある。誘電率の時間変化は試料内部の濃度(あるいは密度)の動的な揺らぎと関係があるので、相関関数を観測することで試料の内部に関する知見をえることができる。例えば式[g (1) (q,τ) ∝ ]は、g (1) (q,τ)が試料内部の濃度揺らぎに比例していることを示している。したがって、動的光散乱では直接測定される散乱光強度の自己相関関数c(τ)からg (1) (q,τ)を求め、g (1) (q,τ)を分析することで試料内部の動的な揺らぎについての知見を得ようとすることが多い。
◇ 動的光散乱と拡散係数 試料内のブラウン運動によって微粒子の(並進)拡散が起きているときには、相関関数 g(q,τ)に対して次式を用いて解析を行い、初期緩和速度Γを決定できる。 Γ = -lim τ→0 [∂ln( g(q,τ) )/∂τ] 具体的にはln( g(q,τ) )をτに対してプロットすると、τ=0 における接線の傾きが-Γとなる。このとき、緩和時間はτR = 1/Γ として求められ、τRから微粒子の(並進)拡散係数を決定できる。 入射光が照射された領域の内部には多数の微粒子が存在し、それぞれの粒子からの散乱光が合わさって検出器に届く。散乱光どうしは互いに干渉し合うので、微粒子の空間分布が時間変化すると検出器から出力される散乱光強度の大きさが時間変化する。動的光散乱はこのような原理によって試料中の微粒子の拡散過程を検出する。 次に試料内部の微粒子の拡散過程について考える。時刻tにおける微粒子の位置をr(t)と表すと、微粒子がブラウン運動によって拡散しているとき、その微粒子の移動距離の2乗の平均は経過時間τに比例するので、拡散の式は次式と表される。 =~ 6*D*τ ここでDは拡散係数である。一方、散乱現象において隣接する微粒子からのそれぞれの散乱光が干渉によって強め合うときの微粒子間距離は2* π /qである。そこで、緩和時間τRは微粒子が強め合う位置から弱め合う位置に距離 π /q だけ移動するのに要する時間であると考える。すなわち、 ~ π /qを用いると拡散の式から直ちに、次式の関係式が得られる(数係数の詳細な計算は割愛)。 D = 1/(q 2 *τR) このことから緩和時間τRから拡散係数Dを決定できることがわかる。ところで、微粒子が剛体球のときには、Einstein-Stokesの式を用いることができる。 D = k B *T / (6* π *η*α) ここで、k B はBoltzmann定数、Tは絶対温度、ηは溶媒の粘性率、αは剛体球の半径である。したがって、この式と[ D = 1/(q 2 *τ R ) ]を用いることによって緩和時間τRから微粒子の半径αが直ちに求まることが分かる。
◇ 試料作製法 動的光散乱を行なうために必要な試料の量は測定装置に依存するが、だいたい1から5 ml程度あれば測定が可能である。希薄な溶液の場合には濃度が1%以下でも有効な測定が十分可能なので、必要な溶質の量はだいたい0.1g以下である。 試料の調製にあたっては、塵の影響を考慮する必要がある。試料中に塵が含まれていると、そのブラウン運動を同時に測定することになってしまう。試料セルはあからじめ塵の含まれていない溶媒で洗浄し、試料セル内部に試料を入れる際にメンブランフィルター(孔径0.2 μm程度)を通す。試料の性質上、フィルターの使用が困難な場合には、それに代わる工夫をするとともに、試料内部に塵が含まれている可能性を常に念頭におく必要がある。溶媒のみを試料として測定し、どの程度塵が含まれているか確認する方法もある。
◇ 測定例 ・微粒子分散液 球形の微粒子の測定の場合は、緩和時間から微粒子の半径を決定することができる。粒径に分布がある場合には、簡便な方法としては、キュムラント解析法がある。これは次式のような式によって非線形最小二乗法によるフィッティングを行なう。 g(q,τ) = exp(-Γ*τ + (K2/2)*τ 2 + (K3/6)*τ 3 ) ここで、K2、K3は定数である。このとき緩和時間の分散はだいたいK2/Γ 2 で与えられ、粒径分布の目安となる。 溶液中の粒子が並進拡散しているときは、緩和速度Γが角度依存性を持つので複数の角度で測定を行い、D=1/(q 2 *τR)に基づいて、縦軸に緩和速度(Γ/s)を横軸に(散乱ベクトル) 2 q 2 /m 2 としてプロットすると、グラフの傾きから拡散係数を正確に決定することができる。 溶液中の粒子が棒状で光学的異方性を持つ場合には、偏光解消散乱(偏光子と検光子をクロスニコルの位置にして行なう光散乱)で測定をし、Γ = -lim τ→0 [∂ln( g(q,τ) )/∂τ] で解析を行なう。得られた緩和時間Γを散乱ベクトルの2乗に対してプロットするとグラフの傾きから並進拡散係数Dt, 切片から回転拡散係数Drが求まる。Kirkwood-Risemannの式より下記の関係がある。 Dt = [kB*T / (3* π *η*L)]*ln(L/d) Dr = [3*k B *T / ( π *η*L 3 )]*ln(L/d) ここでLは棒の長さ、dは直径である。この解析から粒子の異方性を表すアスペクト比 L/d を決定できる。 ・高分子溶液(希薄) 高分子溶液の測定は濃度によっていくつかの領域に分けることができる。高分子の系では、高分子の鎖どうしの溶液中での接触が無視できなくなる重なり合い濃度c*(cは高分子溶液の濃度)を境に高分子鎖の挙動が変化する。 重なり合い濃度より低い濃度(c < c*)の高分子溶液を希薄溶液と呼ぶ。希薄溶液では高分子鎖は孤立した鎖としてブラウン運動をする。動的光散乱で測定される散乱係数は次式となる。D = B*T / (6* π *η*RH)>*(1+k'*c) ここでRHは流体力学的半径、k'は濃度依存性の係数である。測定される拡散係数は濃度に依存するので、複数の濃度で測定し、c → 0 に外挿することで流体力学的半径RHを決定できる。 q > 1/RH となる条件で測定した場合には、高分子の分子内の緩和モードが測定され、拡散の緩和時間よりも速く緩和する。 ・高分子溶液(準希薄) 重なり合い濃度以上では、溶液中に過渡的な網目構造が形成される。しばしば二つの緩和過程が測定されるので、次式のような二つの緩和速度を持つ式を相関関数にフィッティングして解析が行なわれる。 g(q,τ) = A*exp(-Γ1*τ) + (1-A)*exp(-Γ2*τ) ここでAは係数(|A|=1とΓ2はそれぞれ速い緩和速度と遅い緩和速度である。重なり合い濃度以上では、鎖のブラウン運動は過渡的な網目構造で分断され、網目の協同拡散(ゲルモード)に起因する速い緩和Γc(=Γ1)が測定される。Γcは濃度依存性があり、良溶媒のときには図(省略)に示すような挙動を示す。過渡的な網目のサイズは高分子ゲルの場合と同様にして求めることができる。また同時に、過渡的網目構造の中を拡散する鎖の自己拡散に起因する遅い緩和ΓG(=Γ2)が測定される。粘性率の濃度依存性からもc*を決定できるが、動的光散乱で決定したc*よりも高濃度であることが多い。動的光散乱のほうが鎖どうしの重なり合いの効果が低濃度で現われるためであると考えられる。 ・高分子溶液(濃厚領域、融液) 濃厚な場合の相関関数は、緩和時間が広い分布を持つようになるので、Kohlraush-Williams-Watts(KWW)型の伸長指数関数にフィッティングをして解析が行なわれる。 g(q,τ) = exp ここでβ(= 0.7 Indicates a very broad distribution of particle sizes
◆ ゼータ電位(Zeta potential) ・電気二重層のすべり面の電位(表面電位ではない) ・液中での粒子同士の反発性の指標 ・ゼータ電位が0 mV(等電点)に近い場合には凝集や沈殿を起こす ・ゼータ電位が±30 mV以上の場合には分散安定性が良好
低い ← 電解質濃度 → 高い 長い ← 1/κ → 短い 低電解質濃度だと電場は遠くまで届くが、高電解質濃度だと電場は近くにしか届かない > 高電解質濃度領域下では電気泳動しづらい(電場が届かないから) ◇ 電気二重層の目安(z-z型) 濃度 1/κ 0.1 M 1 nm 0.01 M 3 nm 0.001 M 10 nm◆ 測定原理 ・動いている粒子に対して一定周波数の光(レーザー光)を照射する。散乱された光の周波数が、粒子が動く速度に応じて変化するドップラーシフトを利用 ・周波数のシフトΔf = 2*v*sin(θ/2)/λ ここでvは粒子の速度、θは散乱角、λはレーザーの波長である。実際には、レーザー光を2つに分けて、試料の入ったセルを通らない光の周波数F 1 と試料からの散乱光の周波数F 2 を比較して、周波数変化Δf (=F1 - F2)を求める(PALS法)。 ◇ ヘンリーの式 ゼータ電位は計測した電気泳動移動度より次のヘンリーの式を用いて求められる。 UE = (2*ε*z*F(Ka)/(3*η) ここで、UEは電気泳動移動度、zはゼータ電位、εは誘電率、ηは粘度、F(ka)はHenry係数(水系・非水系で切り分け)。 粒子径のパラメータはヘンリーの式に存在していない。電気泳動中は大きな粒子も小さな粒子も一定(移動度)の動きをする。
◆ 測定結果の注意点 理想的ではないPhase plotが得られた場合、以下の理由が考えられる
◇ ゼータ電位値が低い ・pHを前後に振って、変化するかを確認 ・pHを酸性側とアルカリ側に変化させた場合に、チャージが反転する場合は、最初のサンプルは等電点付近で、もともとチャージが弱かったたために、Phase plotの凹凸が小さかった可能性がある
◇ 粒子濃度が高く、電気泳動しづらい環境下にある ・周波数グラフを確認。必要であれば希釈してみる ・電気泳動している粒子は一定に動く(≒周波数変化も一定)。一定の周波数でピークが変化しているか確認する。確認できれば粒子が同じ速度で動いており、電気泳動している粒子を計測できている ・一定の周波数でピークが変化していることが確認できない場合、粒子が様々な速度で動いており、電気泳動している粒子を計測できない ・サンプルの濃度が高い状態で電気泳動させると、電位泳動、粒子同士の衝突、ブラウン運動など様々な粒子の動きが発生する。そのため周波数の変化が一定でなくなる。
◇ セル電極と装置電極の接続が悪く、電場が安定的に形成されていない ・電流/電圧結果のグラフを確認する ・電圧が反転したときに電流もそれに応じて反転しているかを確認する。電極同士が接触いしていない場合、電流が0付近で変化してないことが確認できる。必要であればセルを変える
◆ 試料濃度の影響 ・高濃度でも周波数変化が確認されるが、ゼータ電位が小さくなる場合がある この場合には、粒子間相互作用の影響が示唆される。粒子間相互作用の影響により粒子が動けず電気泳動しづらい場合には、希釈することで粒子間相互作用が弱まり電気泳動し易くなる。 原液などでは、粒子がほとんど動いてはいないが、周波数変化が一定の変化をしていることがある。希釈すると、動き出してゼータ電位値が大きくなる。
◆ 測定範囲(ゼータサイザーの場合) ◇ 粒子径 ・粒子径:3.8 nm - 100 μm ・粒子の光学的特性(屈折率など)に影響される ・5 nm以下の粒子径の場合、濃度等の制約が大きくなる ・測定できる粒子径の上限は、沈降するかどうか ・単一組成の粒子の場合、粒子径が異なっても、ゼータ電位は同じであると考えられる ◇ サンプル濃度 ・サンプル濃度は、光がセルを通り抜けること(13°の前方散乱を捉えるため) ・サンプルによって測定可能な濃度が異なる 粒子径、粒度分布の広さ、光学的特性(屈折率など)、濁度 ・測定に必要な散乱光量が得られること 10 kcps以上の散乱光量が望ましい 粒子径が大きくなれば、より薄い濃度でも測定できる 濃すぎても光がセルを通り抜けられなくなるので測定不可 粒子径が大きくなると、より薄い濃度でも、光が通り抜けられなくなる ・粒子の電気泳動をお互いに阻害しないこと ◇ 粘度 ・粘度が高くなりすぎると電気泳動されなくなる ・サンプル濃度が高すぎるとサンプル粘度が高くなることがある ・濃度・粘度等に問題がある場合は、希釈によって測定が可能になることがある。このときpH, 電気伝導度(イオン濃度)、添加剤濃度に注意が必要 ◇ 電気伝導度 ・200 mS/cm 以下 ・ジュール熱、電気分解が起きないこと、電極へのダメージ、サンプルへのダメージといった問題が防げること ・電気伝導度と設定電圧 < 5 mS/cm: 150 V5 to 30 mS/cm: 50 V > 30 mS/cm: 10 V ◇ 誘電率 ・10以下の場合、もしくは、ポリカーボネイトを溶解させる溶媒を使用するときには、ユニバーサルディップセル(セル材質: PEEK、非水系メイン、水系も使用可能)が必要 ・キャプラリーセル(セル材質:PC、水系メイン、低級アルコールは使用可能{メタノール、エタノール、イソプロピルアルコール})が使用できるのは、誘電率が10以上の場合で、ポリカーボネートが溶解しないとき ・比誘電率11以上は水系として考えても大丈夫
◆ 環境によるゼータ電位の変化 ・pHが変われば、ゼータ電位は変わる アルミナの場合は下記のようになる。 □ 酸性溶液中では、表面のプロトン化によってプラス Al-H2O + \ O / Al-H 2 O + □ 塩基性溶液中では、表面からプロトンが乖離してマイナス Al-O - \ O / Al-O - ・イオン濃度が増加すると、同じ距離での反発エネルギーは減少し、ゼータ電位は低くなる ・イオンの種類によっては吸着する場合もある 「等電点が変化している > 表面が変化している > イオンが吸着している」となる ・ゼータ電位の減少 > ・リポソームが電気泳動しづらくなっている ポリマーによる立体障害は、電荷密度によって異なり、電荷密度が低いと架橋、高いと粒子の電荷を中和する方向に進む ※ 分散媒が異なるとゼータ電位値が異なる。符号が反転する場合もある。 ※ 添加剤は入れすぎてもよくない。
◆ 電気泳動移動度の求め方 ・移動速度 = distance (電極間距離)/ time = μm/s ・電場強度 = volts (印加電圧) / distance = V/cm ・電気泳動移動度の単位 = 移動速度 / 電場強度 = m 2 V -1 s -1 * 10-8 or μmcmV -1 s -1 ------------------------------------------------------------------------------ ■ 参考文献 [1] 野瀬卓平ら『若手研究者のための有機・高分子測定ラボガイド』講談社サイエンティフィック [2] ゼータサイザーナノ講習会資料、Malvern・Panalytical [3] 微粒子工学大系 第1巻 基本技術 -----------------------------------------------------------------------
- 第一原理計算(概論)
- 電子顕微鏡を扱う方向け(SEM, TEM, STEM, EDS)
- 集束イオンビーム(FIB)
- SEM & TEM用試料準備
- 電子線マイクロプローブを扱う方向け(EPMA)
- EPMA(要点)
- EPMA用の試料準備方法
- 電子回折を扱う方向け
- MEM解析を扱う方向け
- STMを扱う方向け
- 分子軌道法
- バンド計算法
- 分子動力学法
- 入力ファイルのQ&A
- 入力ファイルの変換
- cif2cell
- プリミティブセル&空間群の探し方
- 各種の第一原理計算コードの特徴
- IR&Raman
- XAS&EELS
- Akai-KKR
- Akai-KKR: input file (making)
- Akai-KKR: action on error
- Akai-KKR (Igor)
- Akai-KKR (Igor, up-down)
- Akai-KKR (modify)
- Akai-KKR (ref)
- Akai-KKR (sequential search)
- Akai-KKR (igor 4.06, macro)
- Akai-KKR (dos plot)
- Akai-KKR (dos)
- SPR-KKR
- SPR-KKR (GGA+U)
- SPR-KKR (XAS and XMCD)
- SPR-KKR (ARPES)
- ABCAP
- Elk
- Elk (gallery)
- Elk (XAS Mn2CoGa)
- Elk (RDMFT) (Linux package)
- Elk_3Dplot2VESTA
- exciting
- WIEN2kをはじめる前に
- WIEN2k
- WIEN2k tips
- WIEN2k memo
- WIEN2k (GGA+U)
- WIEN2k (XMCD)
- WIEN2k (Mac, MacBook Air)
- WIEN2k for Mac (XAS, XMCD)
- WIEN2k(ver.16)+phonopy(ver.1.11.2)
- WIEN2k(ver.16)+alamode(ver.1.4.1)
- fold2Bloch (WIEN2k)
- DVXa
- DVXa + DVME
- cluster calculation
- クラスターの作り方(XcrySDen+VESTA+displat)
- Igor macro (Gaussian broadening)
- qvasp
- VASP
- VASP (gallery)
- VASP (INCAR format)
- vasp.5.3.3 compile
- vasp.5.4.1 compile
- VASP.5.4.1.05Feb16
- VASP (VCA: 仮想結晶近似)
- VASP+phonopy (Fe2VAl-dfpt)
- VASP+phono3py (Fe2VAl, kappa lattice)
- VASP+phono3py (Fe2VTaAl, kappa lattice)
- VASP (relaxation)
- VASP+phonopy (Fe2VAl-dfpt, 2x2x2 primitive cell)
- phonopy & phono3py
- VASP (sequential calculation)
- VASP (sequential calculation, simple version)
- vasp_phonopy_run
- vasp_phonopy_qha_run
- vasp_phono3py_run
- vasp_phono3py_run_kdeplot
- vasp_boltztrap_run
- vasp_autoset_run
- vasp_autoset_run_band
- vasp_autoset_run_boltztrap
- vasp_autoset_run_boltztrap_phonopy
- vasp_autoset_run_boltztrap_phonopy_band
- vasp_autoset_run_boltztrap_phonopy_phono
- ABINIT
- ABINIT (gallery)
- ABINIT (構造最適化)
- ABINIT (fold2Bloch)
- ABINIT (DFT+U + boltztrap)
- ABINIT (HSE06, PBE0, B3LYP)
- ABINIT (1-shot GW)
- ABINIT (DMFT) (Linux package)
- ABINIT 翻訳
- ABINIT (The Lesson 1-4)
- ABINIT (spin, DFT+U, GW, BSE, RF)
- ABINIT(光触媒研究用)
- ABINIT (input files)
- Ueff data for LDA(GGA)+U
- unfolding (BandUP)
- ABINIT&PWscf+Yambo
- PWscf
- PWscf (Pseudo-potential)
- PWscf (gallery PAW)
- PWscf (gallery USPP)
- PWscf (atomic PP)
- PWscf (XANES)
- PWscf (XANES LiCoO2)
- PWscf (Phonon by DFPT)
- PWscf (Free energy, QHA)
- PWscf v.6.3 (organic semiconductor)
- PWscf v.6.3 (QHA)
- PWscf v.6.3 (Cp)
- PWscf v.6.3 (TDDFT, EELS)
- PWscf v.6.3 (PWcond)
- PWscf v.6.3 (GIPAW)
- PWscf v.6.3 (NEB, ESM, RISM)
- PWscf v.6.3 + phonopy-1.11.2
- PWscf v.6.4 (Raman)
- PWscf v.6.4 (Car-Parrinello)
- PWscf v.6.4 (TDDFPT, GW, BSE, epsilon)
- PWscf v.6.4.1 (UPF,VCA) (ubuntu package)
- PWscf v.6.4.1 (VCA+vacancy)
- PWscf v.7.2 (Wannier90)
- PWscf (EPW)
- PWscf (thermo_py)
- PWscf+BerkeleyGW
- CASINO
- CASINO (density map)
- CASINO (density map OPIUM)
- Phase(2012年6月以降無更新)
- STATE-senri
- ecalj
- GPAW
- BigDFT
- CPMD v.3.17
- CPMD v.4.3
- CP2k v.5.1 (ubuntu package)
- CP2k v.6.1 (Linux package)
- CP2k v.6.1.0
- CP2k v.7.1.0
- NWChem v.6.8.1 (debian package)
- NWChem v.7.0.0 (ubuntu package)
- NWChem v.6.8.1
- NWChem v.7.0.0
- OpenMX (Linux package)
- OpenMX v.3.7
- Siesta v.3.2.5
- Siesta v.4.X
- DFTB+ v.17.1 (win)
- DFTB+ v.18.2 (win)
- DFTB+ v.19.1
- DFTB+ v.20.2.1
- DFTB+ v.21.1
- DFTB+ v.22.1
- DFTB+ v.23.1
- DFTB+ v.24.1
- DFTBaby
- hotbit
- hotcent & tango
- skprogs
- skpar
- adpt_core
- DFTB+ Parameterization Tool set
- xtb v.6.2.3
- sTDA v.1.6.1
- MOPAC v.7
- MOPAC2016
- MOPAC v.22.1.X
- QM/MM
- QM/MM (GAMESS)
- QM/MM (MOPAC + NAMD)
- QM/MM (ORCA + NAMD)
- QM/MM (CPMD + Gromacs)
- QM/MM (ORCA + Gromacs)
- QMMMW (PWscf + Lammps)
- GULP
- GULP (parameter fitting -1-)
- GULP (parameter fitting -2-)
- GULP (converts inputs)
- LAMMPS (Linux)
- LAMMPS (Mac)
- Lammps (phonon)
- Lammps (Examples)
- Lammps (Tutorial&plot)
- Lammps (GCMC, Mg(OH)2)
- Lammps (ZTC)
- Lammps (Zeolite)
- Lammps (GCMC, Zeolite)
- Lammps (hybrid)
- Lammps (GCMC) (for education)
- Lammps (ZTC) (for education)
- Lammps (LJ) (for education)
- Buckingham database
- Lammps (Metal) (for education)
- Lammps (MoS2) (for education)
- Lammps (LiCoO2) (for education)
- Lammps (ReaxFF) (for education)
- Lammps (flow) (for education)
- Lammps (couple)
- Lammps (GNN)
- OCTA 8.4
- Dock 6.11
- OVITO
- Gromacs
- Gromacs (tutorial -1-)
- Gromacs (tutorial -2-)
- Gromacs (tutorial,Jpn)
- Gromos
- NAMD
- Charmm
- Tinker
- AMBER, GROMACS, NAMD
- 固溶体の計算
- supercell (SQS)
- 熱力学(Legendre変換)
- 計算状態図(CALPHAD)
- pycalphad v.0.10.3
- OpenCalphad v.5.0
- OpenCalphad v.6.0
- フェーズフィールド法
- PFDB
- PRISMS-PF v.2.2
- OpenPhase v.0.9.3
- MOOSE
- MMSP
- SymPhas
- mef90/vDef
- 第一原理フェーズフィールド法
- code test
- BoltzTraP
- Phonopy-Spectroscopy
- DynaPhoPy
- Alamode v.1.4.1
- Alamode v.1.2.0
- Alamode v.1.1.0
- ShengBTE
- PhonTS
- Linuxのインストール
- andLinux
- Ubuntu 11.04 (32bit)
- Ubuntu 11.04 (64bit)
- Ubuntu 11.10
- Ubuntu 12.04 LTS (64bit)
- Ubuntu 14.04 LTS (64bit)
- Ubuntu 16.04 LTS (64bit)
- Ubuntu 16.04 LTS (64bit) beginners' course
- CentOS 5 x86_64 (64 bit)
- CentOS 6 i686 (32bit)
- CentOS 6.4 x86_64 (64bit)
- Fedora18 x86_64 (64bit)
- Scientific Linux 32bit
- MacBook Air (Mac OS X 10.8)
- MacBook Pro (Mac OS X 10.7.5)
- コンパイラの設定
- ライブラリの設定
- PC省電力化
- 並列化
- 高速化
- チーム開発方法
- GitHub(共同でのコード開発)
- web検索システム作成方法
- プログラミングの方法
- 数値計算法
- テトラヘドロン法
- 自作コマンド(Linux)
- 分散表示変換プログラム
- 遠隔操作(VNC, putty, ssh)
- サーバー
- XPS, ESCA, UPS, ARPES
- XPS(放射光)
- XPS装置の管理&利用
- XPS(解析)
- AES
- XAFS
- XAFS tips
- Athena&Artemis
- FPMS (XANES)
- FDMNES (XAFS)
- FDMNES (XES)
- XES
- エスカルゴ(写真1)
- エスカルゴ(写真2)
- X線反射率測定(XRR)
- XRR software
- 小角散乱(SAXS)
- SunBIM
- 全散乱 (Total scattering)
- PDFgui
- DISCUS
- EPSR
- ePDF software
- 逆モンテカルロ法(RMC)
- RMCA
- RMC++
- RMC++ (Linux)
- RMCProfile
- HRMC (Hybrid Reverse Monte Carlo)
- fullrmc
- 粉末X線回折(XRD)
- XRD (有機物)
- QualX2 + POW_COD
- EXPO2014
- 中性子粉末回折
- X線及び中性子回折用の群論
- VESTA
- リートベルト法
- GSAS (EXPGUI)
- JANA&FullProf
- MEM解析
- MEM (Alchemy)
- MEM > DOS (test)
- EDA
- ReciPro
- QSTEM
- QSTEM (CUI -1-)
- QSTEM (CUI -2-)
- QSTEM (CUI -3-)
- QSTEM (image fitting)
- QSTEM + Lammps
- QSTEM + MOPAC2016
- Dr.Probe
- Dr.Probe (CUI -1-)
- Prismatic
- Prismatic (CUI -1-)
- EELSModel
- JEM-2100 操作方法(TEM, STEM and EDS)
- ImageJ
- TEM (gallery)
- SEM (gallery)
- STM
- AFM
- グロー放電発光表面分析
- 二次イオン質量分析計(SIMS)
- 固体酸化物型燃料電池(SOFC)
- 水素吸蔵合金
- Zeolite Templated Carbon (ZTC)
- 形成エネルギーの計算法
- sxdefectalign
- CoFFEE
- vickers hardness
- フォノン 理解への道
- フォノン分散の実験的解析
- 比熱
- 伝熱
- 溶解度パラメーター(SP)
- 実験を始める前の心構え
- ゾルゲル法
- 試料作製法
- 真空アーク溶解装置
- 試料作製での理論
- 粉末作製の理論
- メカニカルアロイング
- 粉末冶金・焼結での理論
- 有機合成法
- 生化学
- 無機化合物の特性相関
- 実験における測定法の選択
- 実験装置の分類
- 化学実験での注意事項
- クロマトグラフィー 概要
- クロマトグラフィー GC
- クロマトグラフィー HPLC
- クロマトグラフィー LC/MS
- 核磁気共鳴(NMR)-概要
- 核磁気共鳴(NMR)-理論
- 核磁気共鳴(NMR)-実験
- 核磁気共鳴(NMR)-二次元
- 赤外分光(IR)
- ラマン(Raman)
- ラマン(無機材料)
- 質量分析(MS)
- 蛍光X線分光法(XRF)
- 紫外可視分光(UV-Vis)
- 円二色性(CD)・旋光分散(ORD)
- 光化学
- 時間分解赤外分光法
- 光散乱(DLS, SLS, ゼータ電位)
- 熱分析(DSC, TG-DTA)
- 電子スピン共鳴(ESR)
- 電気化学測定(特にCV)
- CHN分析
- ICP発光分析
- 単結晶X線構造解析
- 固体力学・レオロジー測定
- 内部摩擦
- 電気物性および光・電子機能測定
- ガス吸着法
- pyGAPS
- pyGAPS (DFT fit)
- pyGAPS (Usage, DFT fit)
- pyGAPS (win)
- poreblazer (Surface area calculation)
- Zeo++
- RASPA
- RASPA2
- RASPA3
- MUSIC
- MUSIC (GCMC, MOF-5, H2 gas)
- 研究生活
- 研究生活(実態)
- 研究室の情報(次を探している方向け)
- 労働研究
- 研究開発(R&D)
- 組織学
- 利用申請書の書き方
- 科研費の書き方
- ビジネスマナー(電話対応&メール)
- 英文メール
- 教育指導論
- 英語講演のポイント
- 英語論文執筆のポイント
- 論文執筆
- TeX (LaTeX)
- 英語
- 記憶術
- 大学に行けない人向け
- windows7 home
- Mac
- Software (plot and graph)
- ghostscript(&gnuplot)
- beamer(学会発表用資料の作成)& LaTeX2e
- xmgrace
- Igor
- Igor (ZEM-3 > Excel > Igor)
- Kumi macro
- Vine toolkit ABINIT
- Psi4 (gcolab)
- GAMESS
- Gaussian
- Gaussian (励起状態計算)
- GaussSum
- フロンティア軌道論
- 基礎無機化学 理解への道
- 基礎有機化学 理解への道
- 相図 理解への道
- 真空技術
- 炭素材料と燃焼
- 生体材料
- 超伝導材料
- 光触媒材料
- 熱電変換材料
- 熱伝導率
- 熱伝導率の測定
- 磁気熱量効果
- 有機半導体
- タンパク質計算科学
- 理論計算 理解への道1
- 理論計算 理解への道2
- 理論計算 理解への道3
- 計算機科学 理解への道
- 計算機購入 理解への道
- 空間群 理解への道
- 密度汎関数法 理解への道
- 擬ポテンシャル 理解への道
- atomic 翻訳
- 構造最適化のための理論
- 分子動力学 理解への道
- 分子動力学(MD)でのポテンシャル作成 -1-
- 分子動力学(MD)でのポテンシャル作成 -2-
- eam_database_fit
- adp_database_fit
- potfit
- potfit (table, VASP)
- potfit (table, PWscf)
- potfit (table, xTB)
- MEAMfit v1.02
- atomicrex
- MPC
- 第一原理計算 理論の学習の仕方
- 企業が求める人材
- 放射光
- 原子分解能ホログラフィー
- 3D-AIR-IMAGE
- X線イメージング
- 統計
- 信頼性解析
- データ駆動科学
- 数理的アプローチ
- 機械学習
- 機械学習とpython
- MI
- Google Colaboratory
- XenonPy
- Structure prediction
- CrySPY
- Database
- 拡散係数
- 表面拡散係数
- データベース(物理系)
- データベース(化学系)
- データベース(その他)
- 建築データベース
- 建築
- 化学用フリーソフト
- マジックナンバー
- 研究助成
- 公募
- equipment
- レーザーフラッシュ出力処理コード
- PARATEC
- 新規の研究Q&A
- 結果の解釈方法(まとめ)
- Youtube
- リンク
- アルバム
- ブログ
- 掲示板
- お問合せ
- 小説
- 著名な先生方にtwitterの返信を書くとしたら
- Cygwin
- LINPACK (Benchmark)
- DFT (isotherm)
- 量子コンピューター
- WSL2 (linux on win11)
- NFE-HR
- 名言
- Copilotに書かせてみた1
- Copilotに書かせてみた2
- Copilotに書かせてみた3
- Copilotに書かせてみた4
- Copilotに書かせてみた5
プロフィール
ブログ
アクセス数 ページビュー数- 第一原理計算(概論)
- 電子顕微鏡を扱う方向け(SEM, TEM, STEM, EDS)
- 集束イオンビーム(FIB)
- SEM & TEM用試料準備
- 電子線マイクロプローブを扱う方向け(EPMA)
- EPMA(要点)
- EPMA用の試料準備方法
- 電子回折を扱う方向け
- MEM解析を扱う方向け
- STMを扱う方向け
- 分子軌道法
- バンド計算法
- 分子動力学法
- 入力ファイルのQ&A
- 入力ファイルの変換
- cif2cell
- プリミティブセル&空間群の探し方
- 各種の第一原理計算コードの特徴
- IR&Raman
- XAS&EELS
- Akai-KKR
- Akai-KKR: input file (making)
- Akai-KKR: action on error
- Akai-KKR (Igor)
- Akai-KKR (Igor, up-down)
- Akai-KKR (modify)
- Akai-KKR (ref)
- Akai-KKR (sequential search)
- Akai-KKR (igor 4.06, macro)
- Akai-KKR (dos plot)
- Akai-KKR (dos)
- SPR-KKR
- SPR-KKR (GGA+U)
- SPR-KKR (XAS and XMCD)
- SPR-KKR (ARPES)
- ABCAP
- Elk
- Elk (gallery)
- Elk (XAS Mn2CoGa)
- Elk (RDMFT) (Linux package)
- Elk_3Dplot2VESTA
- exciting
- WIEN2kをはじめる前に
- WIEN2k
- WIEN2k tips
- WIEN2k memo
- WIEN2k (GGA+U)
- WIEN2k (XMCD)
- WIEN2k (Mac, MacBook Air)
- WIEN2k for Mac (XAS, XMCD)
- WIEN2k(ver.16)+phonopy(ver.1.11.2)
- WIEN2k(ver.16)+alamode(ver.1.4.1)
- fold2Bloch (WIEN2k)
- DVXa
- DVXa + DVME
- cluster calculation
- クラスターの作り方(XcrySDen+VESTA+displat)
- Igor macro (Gaussian broadening)
- qvasp
- VASP
- VASP (gallery)
- VASP (INCAR format)
- vasp.5.3.3 compile
- vasp.5.4.1 compile
- VASP.5.4.1.05Feb16
- VASP (VCA: 仮想結晶近似)
- VASP+phonopy (Fe2VAl-dfpt)
- VASP+phono3py (Fe2VAl, kappa lattice)
- VASP+phono3py (Fe2VTaAl, kappa lattice)
- VASP (relaxation)
- VASP+phonopy (Fe2VAl-dfpt, 2x2x2 primitive cell)
- phonopy & phono3py
- VASP (sequential calculation)
- VASP (sequential calculation, simple version)
- vasp_phonopy_run
- vasp_phonopy_qha_run
- vasp_phono3py_run
- vasp_phono3py_run_kdeplot
- vasp_boltztrap_run
- vasp_autoset_run
- vasp_autoset_run_band
- vasp_autoset_run_boltztrap
- vasp_autoset_run_boltztrap_phonopy
- vasp_autoset_run_boltztrap_phonopy_band
- vasp_autoset_run_boltztrap_phonopy_phono
- ABINIT
- ABINIT (gallery)
- ABINIT (構造最適化)
- ABINIT (fold2Bloch)
- ABINIT (DFT+U + boltztrap)
- ABINIT (HSE06, PBE0, B3LYP)
- ABINIT (1-shot GW)
- ABINIT (DMFT) (Linux package)
- ABINIT 翻訳
- ABINIT (The Lesson 1-4)
- ABINIT (spin, DFT+U, GW, BSE, RF)
- ABINIT(光触媒研究用)
- ABINIT (input files)
- Ueff data for LDA(GGA)+U
- unfolding (BandUP)
- ABINIT&PWscf+Yambo
- PWscf
- PWscf (Pseudo-potential)
- PWscf (gallery PAW)
- PWscf (gallery USPP)
- PWscf (atomic PP)
- PWscf (XANES)
- PWscf (XANES LiCoO2)
- PWscf (Phonon by DFPT)
- PWscf (Free energy, QHA)
- PWscf v.6.3 (organic semiconductor)
- PWscf v.6.3 (QHA)
- PWscf v.6.3 (Cp)
- PWscf v.6.3 (TDDFT, EELS)
- PWscf v.6.3 (PWcond)
- PWscf v.6.3 (GIPAW)
- PWscf v.6.3 (NEB, ESM, RISM)
- PWscf v.6.3 + phonopy-1.11.2
- PWscf v.6.4 (Raman)
- PWscf v.6.4 (Car-Parrinello)
- PWscf v.6.4 (TDDFPT, GW, BSE, epsilon)
- PWscf v.6.4.1 (UPF,VCA) (ubuntu package)
- PWscf v.6.4.1 (VCA+vacancy)
- PWscf v.7.2 (Wannier90)
- PWscf (EPW)
- PWscf (thermo_py)
- PWscf+BerkeleyGW
- CASINO
- CASINO (density map)
- CASINO (density map OPIUM)
- Phase(2012年6月以降無更新)
- STATE-senri
- ecalj
- GPAW
- BigDFT
- CPMD v.3.17
- CPMD v.4.3
- CP2k v.5.1 (ubuntu package)
- CP2k v.6.1 (Linux package)
- CP2k v.6.1.0
- CP2k v.7.1.0
- NWChem v.6.8.1 (debian package)
- NWChem v.7.0.0 (ubuntu package)
- NWChem v.6.8.1
- NWChem v.7.0.0
- OpenMX (Linux package)
- OpenMX v.3.7
- Siesta v.3.2.5
- Siesta v.4.X
- DFTB+ v.17.1 (win)
- DFTB+ v.18.2 (win)
- DFTB+ v.19.1
- DFTB+ v.20.2.1
- DFTB+ v.21.1
- DFTB+ v.22.1
- DFTB+ v.23.1
- DFTB+ v.24.1
- DFTBaby
- hotbit
- hotcent & tango
- skprogs
- skpar
- adpt_core
- DFTB+ Parameterization Tool set
- xtb v.6.2.3
- sTDA v.1.6.1
- MOPAC v.7
- MOPAC2016
- MOPAC v.22.1.X
- QM/MM
- QM/MM (GAMESS)
- QM/MM (MOPAC + NAMD)
- QM/MM (ORCA + NAMD)
- QM/MM (CPMD + Gromacs)
- QM/MM (ORCA + Gromacs)
- QMMMW (PWscf + Lammps)
- GULP
- GULP (parameter fitting -1-)
- GULP (parameter fitting -2-)
- GULP (converts inputs)
- LAMMPS (Linux)
- LAMMPS (Mac)
- Lammps (phonon)
- Lammps (Examples)
- Lammps (Tutorial&plot)
- Lammps (GCMC, Mg(OH)2)
- Lammps (ZTC)
- Lammps (Zeolite)
- Lammps (GCMC, Zeolite)
- Lammps (hybrid)
- Lammps (GCMC) (for education)
- Lammps (ZTC) (for education)
- Lammps (LJ) (for education)
- Buckingham database
- Lammps (Metal) (for education)
- Lammps (MoS2) (for education)
- Lammps (LiCoO2) (for education)
- Lammps (ReaxFF) (for education)
- Lammps (flow) (for education)
- Lammps (couple)
- Lammps (GNN)
- OCTA 8.4
- Dock 6.11
- OVITO
- Gromacs
- Gromacs (tutorial -1-)
- Gromacs (tutorial -2-)
- Gromacs (tutorial,Jpn)
- Gromos
- NAMD
- Charmm
- Tinker
- AMBER, GROMACS, NAMD
- 固溶体の計算
- supercell (SQS)
- 熱力学(Legendre変換)
- 計算状態図(CALPHAD)
- pycalphad v.0.10.3
- OpenCalphad v.5.0
- OpenCalphad v.6.0
- フェーズフィールド法
- PFDB
- PRISMS-PF v.2.2
- OpenPhase v.0.9.3
- MOOSE
- MMSP
- SymPhas
- mef90/vDef
- 第一原理フェーズフィールド法
- code test
- BoltzTraP
- Phonopy-Spectroscopy
- DynaPhoPy
- Alamode v.1.4.1
- Alamode v.1.2.0
- Alamode v.1.1.0
- ShengBTE
- PhonTS
- Linuxのインストール
- andLinux
- Ubuntu 11.04 (32bit)
- Ubuntu 11.04 (64bit)
- Ubuntu 11.10
- Ubuntu 12.04 LTS (64bit)
- Ubuntu 14.04 LTS (64bit)
- Ubuntu 16.04 LTS (64bit)
- Ubuntu 16.04 LTS (64bit) beginners' course
- CentOS 5 x86_64 (64 bit)
- CentOS 6 i686 (32bit)
- CentOS 6.4 x86_64 (64bit)
- Fedora18 x86_64 (64bit)
- Scientific Linux 32bit
- MacBook Air (Mac OS X 10.8)
- MacBook Pro (Mac OS X 10.7.5)
- コンパイラの設定
- ライブラリの設定
- PC省電力化
- 並列化
- 高速化
- チーム開発方法
- GitHub(共同でのコード開発)
- web検索システム作成方法
- プログラミングの方法
- 数値計算法
- テトラヘドロン法
- 自作コマンド(Linux)
- 分散表示変換プログラム
- 遠隔操作(VNC, putty, ssh)
- サーバー
- XPS, ESCA, UPS, ARPES
- XPS(放射光)
- XPS装置の管理&利用
- XPS(解析)
- AES
- XAFS
- XAFS tips
- Athena&Artemis
- FPMS (XANES)
- FDMNES (XAFS)
- FDMNES (XES)
- XES
- エスカルゴ(写真1)
- エスカルゴ(写真2)
- X線反射率測定(XRR)
- XRR software
- 小角散乱(SAXS)
- SunBIM
- 全散乱 (Total scattering)
- PDFgui
- DISCUS
- EPSR
- ePDF software
- 逆モンテカルロ法(RMC)
- RMCA
- RMC++
- RMC++ (Linux)
- RMCProfile
- HRMC (Hybrid Reverse Monte Carlo)
- fullrmc
- 粉末X線回折(XRD)
- XRD (有機物)
- QualX2 + POW_COD
- EXPO2014
- 中性子粉末回折
- X線及び中性子回折用の群論
- VESTA
- リートベルト法
- GSAS (EXPGUI)
- JANA&FullProf
- MEM解析
- MEM (Alchemy)
- MEM > DOS (test)
- EDA
- ReciPro
- QSTEM
- QSTEM (CUI -1-)
- QSTEM (CUI -2-)
- QSTEM (CUI -3-)
- QSTEM (image fitting)
- QSTEM + Lammps
- QSTEM + MOPAC2016
- Dr.Probe
- Dr.Probe (CUI -1-)
- Prismatic
- Prismatic (CUI -1-)
- EELSModel
- JEM-2100 操作方法(TEM, STEM and EDS)
- ImageJ
- TEM (gallery)
- SEM (gallery)
- STM
- AFM
- グロー放電発光表面分析
- 二次イオン質量分析計(SIMS)
- 固体酸化物型燃料電池(SOFC)
- 水素吸蔵合金
- Zeolite Templated Carbon (ZTC)
- 形成エネルギーの計算法
- sxdefectalign
- CoFFEE
- vickers hardness
- フォノン 理解への道
- フォノン分散の実験的解析
- 比熱
- 伝熱
- 溶解度パラメーター(SP)
- 実験を始める前の心構え
- ゾルゲル法
- 試料作製法
- 真空アーク溶解装置
- 試料作製での理論
- 粉末作製の理論
- メカニカルアロイング
- 粉末冶金・焼結での理論
- 有機合成法
- 生化学
- 無機化合物の特性相関
- 実験における測定法の選択
- 実験装置の分類
- 化学実験での注意事項
- クロマトグラフィー 概要
- クロマトグラフィー GC
- クロマトグラフィー HPLC
- クロマトグラフィー LC/MS
- 核磁気共鳴(NMR)-概要
- 核磁気共鳴(NMR)-理論
- 核磁気共鳴(NMR)-実験
- 核磁気共鳴(NMR)-二次元
- 赤外分光(IR)
- ラマン(Raman)
- ラマン(無機材料)
- 質量分析(MS)
- 蛍光X線分光法(XRF)
- 紫外可視分光(UV-Vis)
- 円二色性(CD)・旋光分散(ORD)
- 光化学
- 時間分解赤外分光法
- 光散乱(DLS, SLS, ゼータ電位)
- 熱分析(DSC, TG-DTA)
- 電子スピン共鳴(ESR)
- 電気化学測定(特にCV)
- CHN分析
- ICP発光分析
- 単結晶X線構造解析
- 固体力学・レオロジー測定
- 内部摩擦
- 電気物性および光・電子機能測定
- ガス吸着法
- pyGAPS
- pyGAPS (DFT fit)
- pyGAPS (Usage, DFT fit)
- pyGAPS (win)
- poreblazer (Surface area calculation)
- Zeo++
- RASPA
- RASPA2
- RASPA3
- MUSIC
- MUSIC (GCMC, MOF-5, H2 gas)
- 研究生活
- 研究生活(実態)
- 研究室の情報(次を探している方向け)
- 労働研究
- 研究開発(R&D)
- 組織学
- 利用申請書の書き方
- 科研費の書き方
- ビジネスマナー(電話対応&メール)
- 英文メール
- 教育指導論
- 英語講演のポイント
- 英語論文執筆のポイント
- 論文執筆
- TeX (LaTeX)
- 英語
- 記憶術
- 大学に行けない人向け
- windows7 home
- Mac
- Software (plot and graph)
- ghostscript(&gnuplot)
- beamer(学会発表用資料の作成)& LaTeX2e
- xmgrace
- Igor
- Igor (ZEM-3 > Excel > Igor)
- Kumi macro
- Vine toolkit ABINIT
- Psi4 (gcolab)
- GAMESS
- Gaussian
- Gaussian (励起状態計算)
- GaussSum
- フロンティア軌道論
- 基礎無機化学 理解への道
- 基礎有機化学 理解への道
- 相図 理解への道
- 真空技術
- 炭素材料と燃焼
- 生体材料
- 超伝導材料
- 光触媒材料
- 熱電変換材料
- 熱伝導率
- 熱伝導率の測定
- 磁気熱量効果
- 有機半導体
- タンパク質計算科学
- 理論計算 理解への道1
- 理論計算 理解への道2
- 理論計算 理解への道3
- 計算機科学 理解への道
- 計算機購入 理解への道
- 空間群 理解への道
- 密度汎関数法 理解への道
- 擬ポテンシャル 理解への道
- atomic 翻訳
- 構造最適化のための理論
- 分子動力学 理解への道
- 分子動力学(MD)でのポテンシャル作成 -1-
- 分子動力学(MD)でのポテンシャル作成 -2-
- eam_database_fit
- adp_database_fit
- potfit
- potfit (table, VASP)
- potfit (table, PWscf)
- potfit (table, xTB)
- MEAMfit v1.02
- atomicrex
- MPC
- 第一原理計算 理論の学習の仕方
- 企業が求める人材
- 放射光
- 原子分解能ホログラフィー
- 3D-AIR-IMAGE
- X線イメージング
- 統計
- 信頼性解析
- データ駆動科学
- 数理的アプローチ
- 機械学習
- 機械学習とpython
- MI
- Google Colaboratory
- XenonPy
- Structure prediction
- CrySPY
- Database
- 拡散係数
- 表面拡散係数
- データベース(物理系)
- データベース(化学系)
- データベース(その他)
- 建築データベース
- 建築
- 化学用フリーソフト
- マジックナンバー
- 研究助成
- 公募
- equipment
- レーザーフラッシュ出力処理コード
- PARATEC
- 新規の研究Q&A
- 結果の解釈方法(まとめ)
- Youtube
- リンク
- アルバム
- ブログ
- 掲示板
- お問合せ
- 小説
- 著名な先生方にtwitterの返信を書くとしたら
- Cygwin
- LINPACK (Benchmark)
- DFT (isotherm)
- 量子コンピューター
- WSL2 (linux on win11)
- NFE-HR
- 名言
- Copilotに書かせてみた1
- Copilotに書かせてみた2
- Copilotに書かせてみた3
- Copilotに書かせてみた4
- Copilotに書かせてみた5