グラジエント最適化の実務 ― 勾配を緩める操作は、608 時間かけてもカラムを倍にした 17 分に届かない【2026年版】
クロマトグラフィー

グラジエント最適化の実務 ― 勾配を緩める操作は、608 時間かけてもカラムを倍にした 17 分に届かない【2026年版】

分離が足りないとき、まず勾配を緩めるのは合理的に見える。しかし線形溶媒強度モデルを数値で解くと、勾配を緩める操作には天井がある。基準条件から勾配時間を 4,096 倍(608 時間)にしても分離度は 1.40 倍で止まり、同じ改善をカラム長 2 倍なら 17 分で超える。グラジエントは「溶質ごとにばらばらな保持係数を一つの値に揃える装置」であり、揃えたあとに残る手は限られている。

ピークが重なっているとき、多くの人がまず試すのは勾配を緩めることです。10 分の勾配を 20 分にする。 直感的で、装置に触らずにでき、たいてい少しは効きます。

問題は「少しは効く」で止まることです。この操作には天井があります。

線形溶媒強度モデルを数値で解いて確かめました。ある基準条件(100 × 2.1 mm、1.8 µm、分析時間 8.5 分)で 分離度 1.30 のピーク対があるとします。ここから勾配時間だけを延ばしていくと、

  • 2 倍(分析時間 14 分)で分離度 1.51
  • 4 倍(同 24 分)で 1.65
  • 4,096 倍(608 時間 = 25 日、分析時間 324 分)で 1.81

一方、カラムを 2 倍の長さにして勾配時間も 2 倍にすると、分析時間 16.75 分で分離度 1.83 に届きます。

つまり 勾配を緩める操作は、25 日かけても、カラムを倍にして 17 分で得られる分離に届きません。

なぜこうなるのか。そして、では何に時間を使うべきなのか。この記事は、その理由を計算で示します。

この記事の読み方

  • これから学ぶ方:1〜3 節を読んでください。「グラジエントとは何をしている装置なのか」が 1 節で、なぜ緩めても効かないのかが 3 節です。
  • すでに測っている方:4〜7 節が実務の核です。時間の使いどころ(4・5 節)、ドウェル容量が作る 死角(7 節)、日局の換算式が本当は何を保存しているか(8・9 節)。
  • 本記事の数値計算はすべて筆者のシミュレーションです。モデルと仮定は 2 節に明記し、 検証の結果は 10 節と「未確認・限界」に書きました。実測値ではありません。

1. グラジエントは「保持係数を揃える装置」である

グラジエント溶離の説明としてよく見るのは「弱い溶媒から始めて、だんだん強くしていく」というものです。 間違いではありませんが、何のためにそうするのかが抜けています。

逆相 LC の保持は、有機溶媒の体積分率 φ に対しておおむね直線的に効きます。Guillarme らは 線形溶媒強度(LSS: linear solvent strength)モデルの出発点をこう書いています。

log k = log k₀ − S × C

ここで k は保持係数、C は移動相中の有機溶媒の体積分率で 0 から 1 のあいだの値、k₀ は純水中(C = 0) における k の(外挿)値、S は溶媒強度パラメータで、与えられた化合物と固定した実験条件に対して一定である。

S は小分子でおよそ 4〜5 です。S = 5 なら、有機溶媒を 10 % 増やすごとに保持係数が 1/3.2 になる (10^0.5 ≈ 3.16)。20 % 増やせば 1/10 です。保持は組成に対して指数的に効きます。

だから、アイソクラティックでひとつの組成に固定すると、保持係数が 3 桁も 4 桁も違う溶質群が そのまま出てきます。早いものは溶媒ピークに埋もれ、遅いものは何時間も出てきません。

グラジエントがしているのは、この差を潰すことです。カラムの中で組成が上がっていくので、 保持の強い溶質ほど強い溶媒に出会ってから動き出す。結果として、どの溶質もほぼ同じ「実効的な保持係数」で 出てきます。

数値で確かめました。カラム 100 × 2.1 mm・1.8 µm・0.4 mL/min(カラム空隙容量 0.2251 mL、t₀ = 0.5628 分)、 ドウェル容量 0.10 mL、5 → 60 %B を 8.90 分、S = 5 という条件です。

log k₀ 初期組成での k′₀ 溶出時間(分) 溶出時の %B カラム中点での k
1.8 35.5 4.62 28.5 4.476
2.0 56.2 5.24 32.4 4.654
2.2 89.1 5.87 36.3 4.774
2.4 141.3 6.50 40.2 4.853
2.6 223.9 7.14 44.1 4.904
2.8 354.8 7.79 48.1 4.937
3.0 562.3 8.43 52.1 4.958
3.2 891.3 9.08 56.1 4.971

初期組成での保持係数は 35.5 から 891.3 まで 25 倍ちがうのに、カラム中点での保持係数は 4.476 〜 4.971、 わずか 1.11 倍の幅に収まっています。

これがグラジエントの正体です。溶質ごとにばらばらな保持係数を、一つの値に揃えてしまう装置。 だから、グラジエントではどのピークもほぼ同じ幅で出てきます。アイソクラティックのように 「後ろのピークほど太くて低い」ということが起きません。

そして、揃えられたこの値には名前があります。k*(実効保持係数)です。

2. 勾配の急さは、たった一つの数に集約される

勾配の急さを表すのが勾配傾斜 b です。Beyaz らはこう定義しています。

b ≡ S Δφ V_m / (F t_G)

ここで Δφ は初期(φ_o)と最終(φ_f)の溶離強度の差、V_m はカラム空隙容量(mL)、F は流量(mL/min)、 t_G は勾配時間(分)である。

V_m/F はカラム空隙時間 t₀ ですから、b = S・Δφ・t₀ / t_G と書いても同じです。

ここが実務上いちばん重要な点です。b に入っているのは「溶質の性質(S)」「勾配の幅(Δφ)」 「カラム容積(V_m)」「流量(F)」「勾配時間(t_G)」の 5 つだけで、カラム長も内径も粒子径も、 V_m を通してしか効きません。そして k* は b の逆数に比例します。

k* ≈ 1 / (1.15 b)

先ほどの条件では b = 0.1739、1/(1.15b) = 5.000。シミュレーションで実際に測ったカラム中点の保持係数は 4.476〜4.971 でした。理論値とよく合っています。

もうひとつ独立な検証をしておきます。Guillarme らは溶出時点での保持係数 k_e について、 k_e = 1/(2.3 b + 1/k_i) という式を挙げています(k_i は初期組成での保持係数)。 log k₀ = 3.0 の溶質でこの式に入れると 2.4889。筆者のシミュレーションが返した溶出時の k_e は 2.488 でした。有効数字 4 桁で一致します。

本記事のモデルと仮定

計算はすべて、次のモデルの数値解です。外部から借りているのは上の保持式だけで、 実効保持係数もピーク幅も分離度も、すべてシミュレーションから出しています。

  • カラム長を 1 に規格化し、溶質位置 x の常微分方程式 dx/dt = 1/(t₀(1+k)) を陽解法で解く
  • 位置 x・時刻 t の組成は、入口組成を x·t₀ だけ遅らせたもの(溶媒はプラグで進むと仮定)
  • ピーク幅は σ = t₀(1+k_e)/√N。バンド圧縮(G 因子)は 1 と置いた(分離度をやや過小評価する側)
  • 理論段数は N = L/(2.5 dp)、カラム空隙率は 0.65 と仮定
  • 分離度は Rs = (tR₂ − tR₁) / (2(σ₁ + σ₂))
  • カラム外拡散、試料溶媒の影響、二次的保持機構は考慮していない

3. 勾配を緩めても、分離度には天井がある

では本題です。分離が足りないとき、勾配を緩めるとどこまで行けるのか。

条件は 2 節と同じ、選択性 α = 1.05 のピーク対(同じ S で log k₀ が log 1.05 だけ違う 2 成分)です。 基準は k* ≈ 5 になる勾配時間 8.90 分。

勾配時間(分) 基準からの倍率 k*(中点) 最終ピーク(分) 分離度 Rs 基準比
4.45 0.5 7.63 5.14 1.009 0.78
8.90 1 13.54 8.50 1.297 1.00
17.80 2 23.63 14.25 1.513 1.17
35.60 4 40.43 23.84 1.650 1.27
71.20 8 67.36 39.25 1.729 1.33
142.40 16 108.37 62.83 1.771 1.37
284.79 32 166.48 96.44 1.793 1.38
569.59 64 240.94 139.86 1.804 1.39
36,454 4,096 324 1.815 1.40

時間を 2 倍にして得られるのは 17 %。4 倍にしても 27 %。そこから先は、いくら時間を積んでも 40 % で止まります。

最後の行は極端に見えますが、意味があります。勾配時間 608 時間――3 週間以上ポンプを回し続けても、 分離度は 1.815 で頭打ちです。この操作の上限が、ここにあります。

理由は 1 節の性質にあります。勾配を緩めると k* が上がりますが、分離度に効くのは k*/(1+k*) という形です。k* が 5 のとき 0.833、10 で 0.909、50 で 0.980、100 で 0.990。1 に近づくだけで、 1 を超えません。しかも k* が上がるとピーク幅そのものも太るので、正味の取り分はさらに小さくなります。

勾配を緩める操作は、すでに 8 割方使い切っている資源を、残り 2 割ぶん絞り出す操作です。

4. 同じ時間をカラムに使うと、倍以上返ってくる

では、同じ時間を別のことに使ったらどうなるか。分析時間を揃えて比べます。

内容 k* 理論段数 N 最終ピーク Rs 基準比
基準 100 × 2.1 mm / 1.8 µm、tG = 8.90 分 13.5 22,222 8.50 分 1.297 1.00
A 勾配だけ緩める(tG = 22.10 分) 12.2 22,222 16.75 分 1.564 1.21
B カラム長 2 倍(200 mm)+ tG も 2 倍 13.5 44,444 16.75 分 1.834 1.41
C 選択性を α 1.05 → 1.10 に変える 13.5 22,222 8.57 分 2.534 1.95

まったく同じ 16.75 分を使って、勾配を緩めれば 21 % 増、カラムを倍にすれば 41 % 増。差は 2 倍です。

案 B の 1.414 倍という数字は、理論段数が 2 倍になったこと(Rs ∝ √N)とぴったり一致します。 シミュレーションが正しく動いていることの確認にもなっています。

そして 3 節の天井(1.40 倍)を思い出してください。勾配を緩める操作は、時間を無限に使っても 案 B の 1.41 倍に届きません。17 分で済むことに 608 時間かけても負ける、というのはそういう意味です。

なぜカラムのほうが強いのか。分離度は理論段数の平方根に比例し、理論段数はカラム長に比例します。 一方、勾配を緩めて増えるのは k*/(1+k*) だけで、これは 1 で頭打ちになります。 片方は上限のない量、もう片方は上限 1 の量です。

代償は圧力です。同じ流量でカラム長を 2 倍にすれば、背圧も 2 倍になります。 すでに耐圧の上限近くで回しているなら、この手は使えません。そのときは流量を下げることになり、 分析時間はさらに延びます。装置の耐圧が「打てる手の数」を決めている、という話は LC装置本体の選び方 ― カタログの耐圧より、ドウェル容量がメソッドの移管可否を決める【2026年版】 で扱いました。

5. 選択性だけが、時間を使わずに効く

表の案 C は、時間をまったく増やさずに分離度を 1.95 倍にしています。選択性 α を 1.05 から 1.10 に 変えただけです。

分離度への効き方が (α−1)/α という形をしているためで、α = 1.05 では 0.0476、α = 1.10 では 0.0909。 ちょうど 1.91 倍です(シミュレーションは 1.95 倍。保持がわずかに動くぶんの差)。

α をわずかに動かすだけで、時間を 2 倍にしても届かない改善が手に入ります。だからメソッド開発の順序は、 勾配の傾きより先に「選択性を変える操作」を尽くすべきです。具体的には、

  • 有機溶媒の種類を変える(アセトニトリル ⇄ メタノール)。溶出強度だけでなく選択性が変わります
  • pH を変える。解離性化合物では最も大きく効きます
  • 固定相の化学を変える(C18 → フェニル、ペンタフルオロフェニル、極性埋め込み型など)
  • カラム温度を変える

いずれも分析時間を延ばしません。カラム選択の考え方は HPLC/UHPLCカラムの選び方 ― 分離モード・固定相・粒子・細孔径で迷わない実践ガイド【2026年版】 に、 移動相と緩衝液の実務は 移動相・緩衝液調製の実務 ― その pH は、どこで測った値か にまとめてあります。

ただし、選択性の変更には日局の制約がかかります。〈2.00〉は固定相について 「置換基の変更は認められない(例えば,C18 が C8 に変更されるなど)」と明記しています。 公定法の枠内で動かせるのは、pH ±0.2 単位、塩濃度 ±10 %、温度 ±5 ℃(グラジエントの場合)までです。 選択性がいちばん効くのに、規制下ではいちばん動かせない。この非対称は JP19 クロマトグラフィー総論〈2.00〉 ― SN比の改正と、公定法で「条件を変えてよい」範囲【2026年版】 で条文ごとに整理しています。

6. ピークキャパシティも、勾配時間に比例しない

分離度は「特定の 2 本」の話です。多成分試料ではピークキャパシティ(その分析時間に何本のピークを 並べられるか)で考えます。こちらも計算しました。

勾配時間(分) 時間倍率 ピーク幅(分) ピークキャパシティ 基準比
2.22 0.25 0.0245 91.7 0.54
4.45 0.5 0.0339 132.3 0.78
8.90 1 0.0525 170.6 1.00
17.80 2 0.0890 200.9 1.18
35.60 4 0.1598 223.8 1.31
71.20 8 0.2926 244.3 1.43
142.40 16 0.5280 270.7 1.59
284.79 32 0.9058 315.4 1.85

時間を 32 倍にして、ピークキャパシティは 1.85 倍。勾配時間を延ばすとピーク幅も一緒に太るので、 「窓が広がったぶん」がそのまま本数にはなりません。

逆向きに読むと使い道があります。時間を半分にしても、ピークキャパシティは 78 % 残ります。 スクリーニングのように「全部分けきる必要はないが、数をこなしたい」場面では、 勾配を急にして本数を捨てるほうが割に合うことがあります。

本当に本数が要るなら、次元を増やすほうが桁で効きます。二次元 LC は 2D-LC徹底解説 ― ピークキャパシティの積が「絵に描いた餅」になる理由と、その回避策【2026年版】 で扱いました。

7. ドウェル容量のあいだに出る溶質には、勾配が一度も触れていない

ここまでは「勾配が効いている溶質」の話でした。効いていない溶質がいます。

勾配はポンプの混合点で作られ、配管を通ってカラムに届きます。その容量がドウェル容量で、 時間に直すと t_D = V_D/F。この t_D のあいだ、カラムには初期組成が流れ続けます。

Beyaz らは自分たちの保持時間予測式について、こう注記しています。

式(2) は t_o + t_D と t_o + t_D + t_G のあいだに溶出するピークについてのみ有効であることに注意が重要である。

この境界は導出できます。勾配の先端は移動相の速度で進むので、時刻 t における位置は (t − t_D)/t₀。 一方、溶質は初期保持係数 k′₀ で t/(t₀(1+k′₀)) 付近を進みます。追いつく位置を解くと x = t_D / (k′₀ t₀)。これが 1 以上、つまり

k′₀ ≤ t_D / t₀

なら、溶質はカラムを出るまで勾配に追いつかれません。アイソクラティックで出てくるわけです。

数値でも確かめました(t₀ = 0.5628 分、t_D = 0.25 分、境界 0.4442)。

k′₀ 溶出時間(分) 溶出時の %B
0.20 0.675 5.00(初期組成のまま)
0.40 0.788 5.00(初期組成のまま)
0.50 0.844 5.19
0.80 1.000 6.16
2.00 1.494 9.21

では、実機の境界はどこにあるか。

装置のドウェル容量 カラム 流量 t₀(分) t_D(分) 境界 k′₀
0.10 mL 100 × 2.1 mm / 1.8 µm 0.4 0.563 0.250 0.444
0.60 mL 100 × 2.1 mm / 1.8 µm 0.4 0.563 1.500 2.665
2.00 mL 100 × 2.1 mm / 1.8 µm 0.4 0.563 5.000 8.884
2.00 mL 150 × 4.6 mm / 5 µm 1.0 1.620 2.000 1.234
0.10 mL 150 × 4.6 mm / 5 µm 1.0 1.620 0.100 0.062

同じ 2.0 mL のドウェル容量が、従来カラムでは k′₀ = 1.2 までしか食わないのに、 細いカラムでは 8.9 まで食います。7 倍の差です。

これは UHPLC 化の落とし穴そのものです。カラムを小さくすると、装置は何も変えていないのに、 勾配に触れない溶質の範囲が広がります。初期組成での保持が弱い成分――極性代謝物、無機イオン対、 親水性不純物――が、そっくり「勾配前」に落ちます。

対処は 3 つです。

  1. 初期 %B を下げる(k′₀ を上げて境界の外に出す)
  2. ドウェル容量を減らす(配管・ミキサーの見直し)
  3. 境界の内側にいる成分は、グラジエントで最適化しても無駄だと理解する(そこはアイソクラティック領域)

ドウェル容量そのものの測り方と実機の値は LC装置本体の選び方 ― カタログの耐圧より、ドウェル容量がメソッドの移管可否を決める【2026年版】 にまとめました。

8. 日局の 3 段階換算は、実は「k* を固定せよ」と言っている

第十九改正日本薬局方 一般試験法〈2.00〉4.2 は、グラジエント法のカラムを変えるときの手順を 3 段階で定めています。条文はこう書いています。

グラジエントシステムにおける試験条件の変更はイソクラティックシステムの場合より慎重さが求められる.

カラムの大きさを変えること,すなわちカラム容量の変更は,選択性をコントロールするグラジエント容量に 影響する.カラム容量に比例してグラジエント容量を変え,グラジエント条件をカラム容量に合わせて調整する.

換算式は 2 本です。

F₂ = F₁ × [(dc₂² × dp₁) / (dc₁² × dp₂)]

tG₂ = tG₁ × (F₁/F₂) × [(L₂ × dc₂²) / (L₁ × dc₁²)]

条文はこれを「グラジエント容量とカラム容量の比を一定に保つ」ためだと説明しています。 では、その操作は何を保存しているのか。確かめました。

条文の例(150 × 4.6 mm / 5 µm、F = 2.0 mL/min)から(100 × 2.1 mm / 3 µm、F = 0.6947 mL/min)へ 移し、勾配時間を換算式どおり 0.4 倍にすると、

変更前 変更後
勾配傾斜 b 0.22279881 0.22279881 1.0000000000
理論段数 N 12,000 13,333 1.111

b が小数点以下 10 桁まで完全に一致します。b が同じなら k* も同じ。つまり

日局の 3 段階換算は、「k* を厳密に固定したまま、理論段数だけを動かす」操作です。

条文にこの言葉は一度も出てきません。しかしそう読むと、〈2.00〉のもう一方の規定―― L/dp の比を −25 % から +50 % の範囲に保つ――の役割がはっきりします。 L/dp は理論段数の目安ですから、片方の式で k* を凍結し、もう片方の許容幅で理論段数の変動を管理している。 2 つの規定は別々の量を担当しているわけです。

条文の例では L/dp が 30.0 → 33.3(+11.1 %)で、許容幅の内側に収まっています。

9. 換算係数から内径は消える

もうひとつ、条文の 2 本の式を代入して整理すると出てくることがあります。

tG₂/tG₁ に F₂ の式を代入すると、dc² が分子と分母で打ち消し合って、

tG₂ / tG₁ = (L₂ × dp₂) / (L₁ × dp₁)

カラム内径が完全に消えます。条文の例で確かめると (100 × 3)/(150 × 5) = 300/750 = 0.400000。 条文の表に載っている調整因子 0.4 と、丸めなしで一致します。

つまり 3 段階の手順は、最終的に「長さ × 粒子径の比」というひとつの掛け算に潰れます。 内径の変更は流量の式が吸収しているので、勾配時間には影響しません。 実際、内径だけを 4.6 → 2.1 mm に変えた場合の調整因子は 1.000000、つまり勾配表はそのままです。

ただし、これは流量も条文どおり換算した場合に限ります。ポンプの下限や検出器の都合で流量を据え置くと、 話が変わります。

変更 流量も換算 流量は据え置き
150 × 4.6 / 5 µm → 100 × 2.1 / 3 µm 0.400 0.139 2.88 倍
150 × 4.6 / 5 µm → 100 × 4.6 / 3 µm 0.400 0.667 0.60 倍
150 × 4.6 / 5 µm → 150 × 2.1 / 5 µm 1.000 0.208 4.80 倍

流量を換算しなかっただけで、勾配時間が 4.8 倍ずれます。3 段階のうち 1 段を飛ばすと、 k* の固定という肝心の性質が壊れます。

10. 文献の式をそのまま使う前に、対数の底を確認する

最後に、この記事を書きながら実際に踏んだ落とし穴を共有します。同じ記号 S が、文献によって 違う対数の底で定義されています。

  • Guillarme らは log k = log k₀ − S × C(常用対数)と定義しています。式(4) に係数 2.3 が現れるのも 常用対数基準だからです。
  • Beyaz らの論文は保持時間の式を 自然対数(ln)で書き、本文でも「S と ln k′_w」を求めると 書いています。

S の値は基準によって 2.303 倍ちがいます。気づかずに混ぜるとどうなるか、数値で示します。 常用対数基準で決めた S を、そのまま自然対数の式に入れた場合です。

log k₀ 数値シミュレーション 底をそろえた式 底を混ぜた場合 ずれ
2.0 7.9214 7.9214 12.1664 1.54 倍
2.5 11.3452 11.3451 19.4642 1.72 倍
3.0 14.9112 14.9112 27.4613 1.84 倍
4.0 22.5468 22.1541 44.0449 1.95 倍

保持時間が最大 2 倍近くずれます。しかも「なんとなく大きい値」なので、桁違いのようには見えません。 検算しないと通ってしまう種類の誤りです。

底をそろえれば、筆者のシミュレーションと文献の解析式は 0.0001 分まで一致します (log k₀ = 4.0 の行だけずれているのは、この溶質が勾配の窓の外――勾配終了後の保持段階――で 溶出しているためで、7 節に引用した「式(2) の有効範囲」の注記どおりの挙動です)。

教訓は単純です。文献から式を持ってくるときは、記号の定義と対数の底を本文で確認する。 そして、可能なら独立な方法で検算する。

11. では、どういう順番で決めるか

計算結果を実務の手順に落とすと、こうなります。

  1. まず走査グラジエントを 1 本走らせる(5 → 95 %B など広く)。溶出位置から各成分の おおよその挙動が分かります。
  2. 境界 k′₀ = t_D/t₀ を計算する。この値より弱い保持の成分は、勾配で最適化しても動きません。 初期 %B を下げるか、そこはアイソクラティックとして扱います(7 節)。
  3. k* が 5 前後になる勾配時間を選ぶ。b = S·Δφ·t₀/t_G から逆算できます。ここが出発点で、 最適点ではありません。
  4. 足りない分離は、まず選択性で取りにいく。溶媒種・pH・固定相・温度。時間を使いません(5 節)。
  5. それでも足りなければ、理論段数を買う。カラムを長くするか、粒子径を下げます。耐圧と相談。
  6. 勾配を緩めるのは最後。しかも 2 倍にして 17 %、4 倍にして 27 %、そこから先は 40 % が天井 だと知ったうえで(3 節)。
  7. 公定法から動かすときは〈2.00〉4.2 の 3 段階を、3 段とも実行する。流量の換算を省くと k* が動きます(9 節)。

「勾配を緩める」は、悪い手ではありません。安く、速く、副作用が少ない。ただし取り分が小さく、 上限がある操作です。それを知ったうえで最後に使うのと、最初に使って粘るのとでは、 かかる時間がまるで違います。

分析者の視点

  1. グラジエントは「保持係数を揃える装置」だと理解する。だからピーク幅がほぼ一定になる。 アイソクラティックの直感をそのまま持ち込まない。
  2. 勾配の急さは b ひとつに集約される。b = S·Δφ·t₀/t_G。カラム長も内径も粒子径も、 t₀ を通してしか効かない。
  3. k* ≈ 1/(1.15b)。まずここを 5 前後に合わせる。合わせたあとで最適化に入る。
  4. 時間を使う先を間違えない。同じ時間なら、勾配を緩めるより理論段数を買うほうが 2 倍返ってくる。
  5. 選択性だけが時間を使わずに効く。ただし公定法では動かせる幅が最も狭い。
  6. 境界 k′₀ = t_D/t₀ を必ず計算する。この内側の成分に勾配は触れていない。
  7. カラムを小さくすると、装置を変えていなくてもこの境界が広がる。UHPLC 化のときに必ず確認する。
  8. 日局の換算は 3 段とも実行する。流量の換算を飛ばすと勾配時間が最大 4.8 倍ずれる。
  9. 文献の式は対数の底を確認してから使う。S の基準違いで保持時間が 2 倍近くずれる。
  10. モデルはモデルである。本記事の数値は自分の試料の予測値ではない。方向と大きさの感覚を持つためのもの。

まとめ

  • グラジエントは、溶質ごとにばらばらな保持係数を一つの値 k* に揃える装置である。 初期保持が 25 倍ちがう溶質群でも、カラム中点での保持係数は 1.11 倍の幅に収まる。
  • 勾配の急さは b = S·Δφ·t₀/t_G ひとつに集約され、k* ≈ 1/(1.15b)
  • 勾配を緩める操作には天井がある。基準から時間 2 倍で分離度 +17 %、4 倍で +27 %、 4,096 倍(608 時間)でも +40 % で止まる。
  • 同じ分析時間 16.75 分で、勾配を緩めれば +21 %、カラム長を 2 倍にすれば +41 %。 後者は前者の天井すら超える。
  • 選択性は時間を使わずに効く。α を 1.05 → 1.10 にするだけで分離度は約 2 倍。 ただし公定法では最も動かせない部分である。
  • ピークキャパシティも比例しない。時間 32 倍で 1.85 倍。逆に時間を半分にしても 78 % 残る。
  • k′₀ ≤ t_D/t₀ の溶質には、勾配が一度も触れていない。旧型装置+細いカラムでは境界が k′₀ = 8.9 まで広がる(従来カラムなら 1.2)。
  • 日局〈2.00〉4.2 の 3 段階換算は、b(=k*)を小数点以下 10 桁まで保存する操作であり、 換算係数は (L₂·dp₂)/(L₁·dp₁) に潰れて内径が消える。ただし流量を換算しなければ壊れる。
  • 文献の S は対数の底が統一されていない。混ぜると保持時間が最大 1.95 倍ずれる。

用語集

  • アイソクラティック溶離:移動相組成を一定に保つ溶離。
  • グラジエント溶離:分析中に移動相組成を変えていく溶離。
  • 保持係数 k:溶質がカラム内で固定相にとどまっていた時間と、移動相とともに動いていた時間の比。
  • k*(実効保持係数):グラジエントで溶質がカラム中点を通過するときの保持係数。溶質によらずほぼ一定になる。
  • 勾配傾斜 b:勾配の急さを表す無次元量。b = S·Δφ·t₀/t_G。
  • S(溶媒強度パラメータ):組成変化に対する保持の感度。小分子でおよそ 4〜5。対数の底の定義に注意。
  • LSS(線形溶媒強度)モデル:log k が有機溶媒の体積分率に対して直線になるとする保持モデル。
  • t₀(カラム空隙時間):保持されない成分がカラムを通り抜ける時間。カラム空隙容量 ÷ 流量。
  • ドウェル容量(グラジエント遅延容量):移動相の混合点からカラム入口までの容量。時間に直すと t_D。
  • ピークキャパシティ:ある分析時間のなかに、分離度 1 で並べられるピークの本数。
  • 選択性 α:隣り合う 2 成分の保持係数の比。
  • 理論段数 N:カラムの分離効率を表す量。分離度は √N に比例する。

出典

  1. 第十九改正日本薬局方 一般試験法〈2.00〉クロマトグラフィー総論(厚生労働省告示第193号)。 4.2「液体クロマトグラフィー:グラジエント溶離」の条文、流量・グラジエント時間・注入量の換算式、 条文中の換算例、グラジエント遅延容量の補正式。 https://www.mhlw.go.jp/content/11120000/001689414.pdf
  2. D. Guillarme, T. Bouvarel, F. Rouvière, S. Heinisch, "A simple mathematical treatment for predicting linear solvent strength behavior in gradient elution: Application to biomolecules", Journal of Separation Science 2022, 45(17), 3276–3285. doi:10.1002/jssc.202200161 (LSS の保持式、勾配傾斜 b の定義、溶出時保持係数 k_e の式。オープンアクセス) https://pmc.ncbi.nlm.nih.gov/articles/PMC9543774/
  3. A. Beyaz, W. Fan, P. W. Carr, A. P. Schellinger, "Instrument parameters controlling retention precision in gradient elution reversed-phase liquid chromatography", Journal of Chromatography A 2014, 1371, 90–105. doi:10.1016/j.chroma.2014.09.085(保持時間予測式とその有効範囲、勾配傾斜の定義) https://pmc.ncbi.nlm.nih.gov/articles/PMC4388777/

数値計算は筆者による。モデルと仮定は 2 節に、限界は「未確認・限界」に記載した。

関連記事

未確認・限界

  • 本記事の数値はすべてモデル計算であり、実測値ではない。特定の試料・カラムの予測値として 使ってはならない。方向と大きさの感覚を得るためのものである。
  • S = 5、α = 1.05、カラム空隙率 0.65、換算段高 2.5 はいずれも仮定値である。実際の値は 化合物・固定相・移動相によって変わる。ただし結論(勾配を緩める操作に天井がある)は これらの値の選び方に依存しない。
  • バンド圧縮因子 G を 1 と置いた。実際にはグラジエントでピークはわずかに圧縮されるので、 本記事の分離度・ピークキャパシティはやや過小評価の側にある。
  • カラム外拡散、試料溶媒の影響、二次的保持機構、カラム内での勾配の分散は考慮していない。
  • LSS モデルは直線性を前提にしている。Guillarme らも「保持モデルが十分に直線的であること」を 条件に挙げ、「これらの条件は小分子やペプチドでは常に満たされるわけではない」と明記している。
  • 案 B(カラム長 2 倍)の圧力の影響は定性的にしか扱っていない。実機では耐圧と摩擦発熱が 効いてくるが、それらは本記事のモデルに入っていない。
  • メーカー公式のグラジエント設計資料(Agilent の技術資料など)と LCGC の Dolan による 連載記事は、いずれも HTTP 403 で取得できなかった。したがって「k* を 5 前後にする」という 実務上の目安について、本記事は一次資料を引用していない。b と k* の関係式は上記論文 2 件から取り、 目安の値そのものは筆者の計算(k* = 5 が表の中でどこに位置するか)として提示している。
  • USP〈621〉および Ph. Eur. 2.2.46 の原文は有償のため未入手。本記事の薬局方の記述は 日本薬局方〈2.00〉のみに基づく。三極間で許容幅が一致しているかは確認していない。
  • 温度勾配・流量勾配・多段勾配は扱っていない。本記事は単一の直線勾配に限る。

制作メモ:ディープリサーチで収集→編集側で原著検証のうえ執筆。数値計算は筆者が実装し、 独立な 3 つの方法(文献の解析式との突き合わせ、理論値 1/(1.15b) との比較、 カラム長 2 倍で分離度が厳密に √2 倍になることの確認)で検証した。