Rouse-Modell

Das Rouse-Modell ist eines der einfachsten, in der Polymerphysik verwendeten Modelle für die Dynamik von Polymeren. Es wird für Polymere verwendet, die keine Verschlaufungen aufweisen. Die Dynamik verschlaufter Polymere wird meist durch das Reptationsmodell beschrieben.

Das Rouse-Modell beschreibt das Polymer als ideale Kette von Massenpunkten (oft als englisch beads bezeichnet), die durch Federn verbunden sind. Die zeitliche Veränderung der Konformation wird dann durch brownsche Bewegung der einzelnen Massenpunkte realisiert, indem an jedem eine zufällige (thermische) Kraft angreift. Dabei werden keine Volumenausschlusseffekte beachtet, das Polymer kann sich also mit sich selbst überkreuzen. Das Modell wurde 1953 von Prince E. Rouse vorgeschlagen.^[1]

Das Rouse-Modell führt zu folgender stochastischer Differentialgleichung (Langevin-Gleichung) für die Position R → n {\displaystyle {\vec {R}}_{n}} des n {\displaystyle n} -ten Massenpunktes:

d R → n d t = k ζ ⋅ ( R → n − 1 − R → n + R → n + 1 − R → n ) ⏟ Wechselwirkung mit Nachbargliedern + f → n ( t ) ⏟ Zufallskraft {\displaystyle {\frac {\mathrm {d} {\vec {R}}_{n}}{\mathrm {d} t}}=\underbrace {{\frac {k}{\zeta }}\cdot \left({\vec {R}}_{n-1}-{\vec {R}}_{n}+{\vec {R}}_{n+1}-{\vec {R}}_{n}\right)} _{\text{Wechselwirkung mit Nachbargliedern}}+\underbrace {{\vec {f}}_{n}(t)} _{\text{Zufallskraft}}}

Dabei ist k {\displaystyle k} die Federkonstante der Federn im Modell, ζ {\displaystyle \zeta } der Reibungskoeffizient eines Massenpunkts im Lösungsmittel und N {\displaystyle N} die Anzahl der Kettenglieder. Der Term k ( R → n ± 1 − R → n ) {\displaystyle k({\vec {R}}_{n\pm 1}-{\vec {R}}_{n})} ist jeweils die lineare Rückstellkraft durch die Feder zum Vorgänger n − 1 {\displaystyle n-1} und zum Nachfolger n + 1 {\displaystyle n+1} (an den beiden freien Enden des Polymers entfällt jeweils einer der beiden Terme). Der Term f → n ( t ) {\displaystyle {\vec {f}}_{n}(t)} beschreibt eine thermische Zufallskraft (Brownsche Molekulardynamik), die keine Vorzugsrichtung besitzt und räumlich und zeitlich unkorreliert ist. Aus diesem Ansatz ergeben sich folgende Eigenschaften des Polymers:

• Diffusionskoeffizient des Massenschwerpunkts: D G = k B T N ζ {\displaystyle D_{G}={\frac {k_{\mathrm {B} }T}{N\zeta }}} ( T {\displaystyle T} : Temperatur, k B {\displaystyle k_{\mathrm {B} }} : Boltzmann-Konstante)
• Rotationsrelaxationszeit: τ R = ζ N 2 l 2 3 π 2 k B T {\displaystyle \tau _{R}={\frac {\zeta N^{2}l^{2}}{3\pi ^{2}k_{\mathrm {B} }T}}}
• mittlere quadratische Verschiebung^[2] eines einzelnen Segments:

⟨ R → n 2 ( τ ) ⟩ = ⟨ ( R → n ( t + τ ) − R → n ( t ) ) 2 ⟩ ≈ 2 N l 2 π 3 / 2 τ τ R {\displaystyle \left\langle {\vec {R}}_{n}^{2}(\tau )\right\rangle =\left\langle \left({\vec {R}}_{n}(t+\tau )-{\vec {R}}_{n}(t)\right)^{2}\right\rangle \approx {\frac {2Nl^{2}}{\pi ^{3/2}}}{\sqrt {\frac {\tau }{\tau _{R}}}}}

Eine wichtige Erweiterung wurde 1956 von Bruno Zimm veröffentlicht:^[3] Sein Modell (oft einfach als "Zimm Modell" bezeichnet) berücksichtigt auch hydrodynamische Wechselwirkungen zwischen den Massenpunkten der Kette. Diese sind Wechselwirkungen (Kräfte), die durch die Lösungsmittelmoleküle um das Polymer vermittelt werden: Die Massenpunkte ziehen die Lösungsmittelmoleküle bei ihrer Bewegung mit sich, was auch zu einer Kraft auf benachbarte Kettenglieder führt (siehe Abbildung rechts). Das Zimm-Modell führt zu einer besseren Beschreibung realer Polymere als das Rouse-Modell, die auch mit experimentellen Daten für verdünnte Lösungen gewisser Polymere übereinstimmt.

Die obige Langevin-Gleichung für das Rouse-Modell wird um einen Tensor (Matrix) H n m {\displaystyle \mathrm {H} _{nm}} erweitert, die die hydrodynamische Kraft zwischen dem n {\displaystyle n} -ten und m {\displaystyle m} -ten Segment darstellt:

d R → n d t = k ⋅ ∑ m H n m ( R → n − 1 − R → n + R → n + 1 − R → n ) ⏟ Wechselwirkung mit Nachbargliedern+Hydrodynamik + f → n ( t ) ⏟ Zufallskraft {\displaystyle {\frac {\mathrm {d} {\vec {R}}_{n}}{\mathrm {d} t}}=\underbrace {k\cdot \sum \limits _{m}\mathrm {H} _{nm}\left({\vec {R}}_{n-1}-{\vec {R}}_{n}+{\vec {R}}_{n+1}-{\vec {R}}_{n}\right)} _{\text{Wechselwirkung mit Nachbargliedern+Hydrodynamik}}+\underbrace {{\vec {f}}_{n}(t)} _{\text{Zufallskraft}}}

Dabei ist zu beachten, dass der Tensor H n m {\displaystyle \mathrm {H} _{nm}} von den Positionen R → 0 , . . . , R → N − 1 {\displaystyle {\vec {R}}_{0},...,{\vec {R}}_{N-1}} aller Segmente abhängt. Dadurch ist die obige Langevin-Gleichung nichtlinear und kann nicht mehr einfach gelöst werden. Bruno Zimm ersetzte H n m ( R → 0 , . . . , R → N − 1 ) {\displaystyle \mathrm {H} _{nm}({\vec {R}}_{0},...,{\vec {R}}_{N-1})} daher durch seinen Gleichgewichtsmittelwert ⟨ H n m ⟩ e q {\displaystyle \langle \mathrm {H} _{nm}\rangle _{eq}} , der berechnet werden kann. Daraus leiten sich folgende Eigenschaften eines Zimm-Polymers her:

• Diffusionskoeffizient des Massenschwerpunkts: D G = 8 k B T 3 6 π 3 η s N ⋅ l {\displaystyle D_{G}={\frac {8k_{\mathrm {B} }T}{3{\sqrt {6\pi ^{3}}}\eta _{s}{\sqrt {N}}\cdot l}}} ( η s {\displaystyle \eta _{s}} : Viskosität des Lösungsmittels)
• Rotationsrelaxationszeit: τ R = η S ( N l ) 3 3 π k B T {\displaystyle \tau _{R}={\frac {\eta _{S}({\sqrt {N}}l)^{3}}{{\sqrt {3\pi }}k_{\mathrm {B} }T}}}
• mittlere quadratische Verschiebung eines einzelnen Segments: ⟨ R → n 2 ( τ ) ⟩ = 2 Γ ( 1 / 3 ) N l 2 π 2 ( τ τ R ) 2 / 3 {\displaystyle \left\langle {\vec {R}}_{n}^{2}(\tau )\right\rangle ={\frac {2\Gamma (1/3)Nl^{2}}{\pi ^{2}}}\left({\frac {\tau }{\tau _{R}}}\right)^{2/3}}

In diesem Abschnitt soll auf einige reale Polymere hingewiesen werden, die sich gut mit einem der oberen Modelle beschreiben lassen:

• Einzelsträngige DNA ist ein relativ biegsames Polymer und zeigt auf kurzen Zeitskalen in verdünnter Lösung Zimm-artige Dynamik für einzelne Segmente, wie mit Fluoreszenz-Korrelations-Spektroskopie gezeigt werden konnte.^[4]
• Doppelsträngige DNA ist deutlich steifer als einzelsträngige DNA, daher spielen die hydrodynamischen Wechselwirkungen eine wesentlich geringere Rolle, sodass ihre Monomerdynamik in verdünnter Lösung mit dem Rouse-Modell gut beschrieben werden kann.^[4]

• Michael Rubinstein: Polymer Physics Oxford University Press, USA, 2003, ISBN 978-0-19-852059-7
• Masao Doi: Introduction to polymer physics Oxford University Press, USA, 1996, ISBN 0-19-851772-6
• Masao Doi, S. F. Edwards: The theory of polymer dynamics Oxford University Press, USA, 1988, ISBN 978-0-19-852033-7

1. ↑ Prince E. Rouse, A Theory of the Linear Viscoelastic Properties of Dilute Solutions of Coiling Polymers, J. Chem. Phys. 21, 1272 (1953), cited over 1000 times by 2010.
2. ↑ Die mittlere quadratische Auslenkung ⟨ r ( τ ) ⟩ {\displaystyle \left\langle r(\tau )\right\rangle } (englisch mean squared displacement oder kurz MSD) ist eine gute Größe zur Charakterisierung der Dynamik eines sich zufällig bewegenden Teilchens. Sie misst den mittleren Abstand, den ein Teilchen in einer Zeitspanne τ {\displaystyle \tau } von seinem Ursprung zurücklegt. Für ein frei diffundierendes Teilchen (in 3 Dimensionen) ergibt sich die Beziehung ⟨ r 2 ( τ ) ⟩ = 6 D ⋅ τ {\displaystyle \left\langle r^{2}(\tau )\right\rangle =6D\cdot \tau } mit dem Diffusionskoeffizienten D {\displaystyle D} . Abweichungen von diesem linearen Verhalten (normale Diffusion) werden als anomale Diffusion bezeichnet. Dann gilt oft (wie auch im vorliegenden Fall) die allgemeinere Form ⟨ r 2 ( τ ) ⟩ = Γ ⋅ τ α {\displaystyle \left\langle r^{2}(\tau )\right\rangle =\Gamma \cdot \tau ^{\alpha }} .
3. ↑ Bruno H. Zimm, Dynamics of Polymer Molecules in Dilute Solution: Viscoelasticity, Flow Birefringence and Dielectric Loss, J. Chem. Phys. 24, 269 (1956).
4. ↑ ^{a b} Roman Shusterman, Sergey Alon, Tatyana Gavrinyov, Oleg Krichevsky: Monomer Dynamics in Double- and Single-Stranded DNA Polymers. In: Physical Review Letters. Band 92, Nr. 4, Januar 2004, ISSN 0031-9007, doi:10.1103/PhysRevLett.92.048303 (ens.fr [PDF; 158 kB]).