Chaîne idéale

La chaîne idéale est le modèle le plus simple utilisé pour décrire les polymères^[1], comme les acides nucléiques et les protéines. Le polymère est considéré comme une marche aléatoire^[2], et toute interaction entre les monomères est négligée. Bien que ce modèle soit élémentaire, il donne un aperçu de la physique des polymères.

Dans ce modèle, les monomères sont des tiges rigides de longueur fixée l, et leur orientation est complètement indépendante des positions et orientations de leurs voisins, dans la mesure où deux polymères peuvent coexister à la même place. Dans certains cas, le monomère peut être interprété physiquement comme un amino-acide dans un polypeptide. Dans d'autres cas, un monomère est simplement un segment de polymère, représenté comme une unité libre, discrète et donc, l est la longueur de Kuhn. Par exemple, la chromatine est modélisée comme un polymère, pour lequel chaque segment mesure environ 14 à 46 kilobases de long^[3].

Modèle

Soit un polymère constitué de N unités monomères, dont la longueur totale une fois déplié est :

${\displaystyle L=Nl}$.

Dans cette approche simple, où l'on considère que les monomères ne sont pas en interaction, l'énergie du polymère est choisie de telle sorte qu'elle soit indépendante de sa forme. Ceci qui signifie qu'à l'équilibre thermodynamique, toutes les configurations du polymère ont autant de chances les unes que les autres de s'établir lors du mouvement du polymère, relativement à la distribution de Maxwell–Boltzmann.

Soit ${\displaystyle {\vec {R}}}$ le vecteur d'un bout à l'autre de la chaîne idéale, et ${\displaystyle {\vec {r_{1}}},...,{\vec {r_{N}}}}$ les vecteurs correspondant à chaque monomère. Ces vecteurs aléatoires ont des composantes dans les trois directions de l'espace. La plupart des expressions de cet article suivent l'hypothèse que le nombre de monomères ${\displaystyle N}$ est grand, c'est pourquoi le théorème central limite s'applique. La figure ci-dessous représente le schéma d'une chaîne idéale (courte).

Les deux extrémités de la chaîne ne coïncident pas, mais fluctuent l'une autour de l'autre, d'où :

${\displaystyle \langle {\vec {R}}\rangle =\Sigma _{i=1}^{N}\langle {\vec {r_{i}}}\rangle ={\vec {0}}}$.

Tout au long de cet article, les symboles ${\displaystyle \langle \rangle }$ désigneront la valeur moyenne (par rapport au temps) de variables ou vecteurs aléatoires, comme ci-dessus.

Comme ${\displaystyle {\vec {r_{1}}},...,{\vec {r_{N}}}}$ sont indépendants, il vient du théorème central limite que ${\displaystyle {\vec {R}}}$ suit une distribution normale^[4] (ou distribution de Gauss) : précisément, en 3D, ${\displaystyle R_{x},R_{y}}$ et ${\displaystyle R_{z}}$ suivent une distribution normale de moyenne ${\displaystyle 0}$ et de variance :

${\displaystyle \sigma ^{2}=\langle R_{x}^{2}\rangle -\langle R_{x}\rangle ^{2}=\langle R_{x}^{2}\rangle -0}$.

${\displaystyle \langle R_{x}^{2}\rangle =\langle R_{y}^{2}\rangle =\langle R_{z}^{2}\rangle ={\frac {Nl^{2}}{3}}}$.

Et ainsi : ${\displaystyle \langle {\vec {R}}^{2}\rangle =Nl^{2}=Ll}$. Le vecteur d'un bout à l'autre de la chaîne suit la fonction de densité de probabilité suivante :

${\displaystyle P({\vec {R}})=\left({\frac {3}{2\pi Nl^{2}}}\right)^{3/2}e^{-{\frac {3{\vec {R}}^{2}}{2Nl^{2}}}}}$.

La distance moyenne d'une extrémité à l'autre du polymère est donnée par :

${\displaystyle {\sqrt {\langle {\vec {R}}^{2}\rangle }}={\sqrt {N}}l={\sqrt {Ll}}}$.

Le rayon de giration est une quantité fréquemment utilisée en physique des polymères, et peut s'exprimer comme suit :

${\displaystyle R_{G}={\frac {{\sqrt {N}}l}{\sqrt {6}}}}$.

Il est important de noter que la distance moyenne d'une extrémité à l'autre du polymère, donnée ci-dessus, qui dans le cas de notre modèle simplifié représente typiquement l'amplitude des fluctuations du système, devient négligeable comparée à la longueur totale du polymère déplié à la limite thermodynamique. Ce résultat est une propriété générale des systèmes statistiques.

Remarque mathématique : la démonstration rigoureuse de l'expression de la densité de probabilité ${\displaystyle P({\vec {R}})}$ n'est pas si évidente qu'elle l'est apparue précédemment : de l'application du théorème central limite, il apparaît que ${\displaystyle R_{x},R_{y}}$ et ${\displaystyle R_{z}}$ suivent une loi de distribution normale centrée de la variance ${\displaystyle Nl^{2}/3}$. L'expression donnée ci-dessus pour ${\displaystyle P({\vec {R}})}$ n'est pas la seule qui soit compatible avec une telle distribution pour ${\displaystyle R_{x},R_{y}}$ et ${\displaystyle R_{z}}$. Cependant, comme les composantes des vecteurs ${\displaystyle {\vec {r_{1}}},...,{\vec {r_{N}}}}$ sont décorrélées pour le trajet aléatoire que nous considérons, il vient que ${\displaystyle R_{x},R_{y}}$ et ${\displaystyle R_{z}}$ sont également décorrélés. Cette condition additionnelle peut uniquement être remplie si ${\displaystyle {\vec {R}}}$ suit la distribution de ${\displaystyle P({\vec {R}})}$. Autrement, ce résultat peut être démontré par l'application d'une généralisation multidimensionnelle du théorème central limite, ou par des arguments de symétrie.

Généralités du modèle

Alors que le modèle élémentaire décrit ci-dessus est totalement inapproprié pour la description des polymères réels à l'échelle microscopique, il peut cependant faire preuve de pertinence à l'échelle macroscopique dans le cas d'un polymère en solution, dont les monomères forment un mélange idéal avec le solvant (car dans ce cas, les interactions entre monomère et monomère, entre molécule de solvant et molécule de solvant, et entre monomère et molécule de solvant sont identiques, de telle sorte que l'énergie du système peut être considérée constante, validant les hypothèses du modèle). La pertinence de ce modèle est, malgré cela, limitée à l'échelle macroscopique, par le fait qu'il néglige les effets stériques (c'est-à-dire qu'il ne prend pas en compte la place qu'occupe chaque atome dans une molécule).

Les autres modèles de polymères fluctuants ne prenant pas en compte les interactions entre monomères ni les effets stériques, comme le worm like chain model, convergent tous asymptotiquement vers le modèle de chaîne idéale une fois la limite thermodynamique atteinte. Pour cette analogie, on introduit un segment de Kuhn, correspondant à un monomère de longueur équivalente, à considérer dans la chaîne idéale analogue. Le nombre de segments de Kuhn à considérer dans la chaîne idéale analogue est égal au quotient de la longueur totale du polymère déplié, par la longueur d'un segment de Kuhn.

Élasticité entropique d'une chaîne idéale

Si les deux extrémités libres d'une chaîne idéale sont attachées par un dispositif de micro-manipulation, alors le dispositif subit une force exercée par le polymère. L'énergie de la chaîne idéale est constante, et donc sa moyenne temporelle, l'énergie interne, est aussi constante, ce qui signifie que cette force provient nécessairement d'un effet purement entropique.

La force entropique est très similaire à la pression subie par les parois d'une boîte contenant un gaz idéal. L'énergie interne d'un gaz idéal dépend uniquement de sa température, et pas du volume de la boîte qui le contient. Ce n'est donc pas un effet énergétique qui tend à augmenter le volume de la boîte comme le fait la pression du gaz. Ceci implique que la pression d'un gaz idéal a une origine purement entropique^[5].

Quelle est donc l'origine microscopique d'une telle force ou d'une telle pression entropiques ? La réponse la plus commune concerne l'effet des fluctuations thermiques, qui tend à amener un système thermodynamique vers un état macroscopique, qui correspond à un maximum dans le nombre d'états microscopiques (ou micro-états) qui sont compatibles avec cet état macroscopique. En d'autres termes, les fluctuations thermiques tendent à amener un système vers son état macroscopique d'entropie maximale.

Que cela signifie-t-il dans le cas de la chaîne idéale ? Tout d'abord, pour la chaîne idéale, un état microscopique est caractérisé par une superposition des états individuels de chaque monomère ${\displaystyle {\vec {r}}_{i}}$(avec i variant de 1 à N). Dans son solvant, la chaîne est constamment sujette à des chocs dus aux mouvements des molécules de solvant, et chacun de ces chocs envoie le système de son état microscopique actuel jusqu'à un autre état microscopique, très similaire. Pour un polymère idéal, comme cela a été montré précédemment, il y a plus d'états microscopiques compatibles dans le cas d'une distance d'un bout à l'autre du polymère faible, que dans le cas d'une longue distance. Ainsi, pour une chaîne idéale, maximiser l'entropie signifie réduire la distance entre ses deux extrémités libres. De ce fait, une force qui essaie de plier la chaîne est exercée par la chaîne idéale entre ses deux extrémités^[6].

Dans cette section, la moyenne de cette force sera dérivée. L'expression générale obtenue à la limite thermodynamique sera ensuite discutée.

Chaîne idéale sous contrainte de longueur : ensemble de Helmholtz

Dans cette sous-section, nous considérerons le cas d'une chaîne idéale dont les deux extrémités sont attachées à des points fixes. Le vecteur ${\displaystyle {\vec {R}}}$ joignant ces deux points caractérise l'état macroscopique (ou macro-état) de la chaîne idéale. Chaque macro-état correspond à un certain nombre de micro-états, que nous appellerons ${\displaystyle \Omega ({\vec {R}})}$ (les micro-états sont définis dans l'introduction de cette section). Comme l'énergie de la chaîne idéale est constante, chacun de ces micro-états a autant de chances qu'un autre de s'établir. L'entropie associée à un macro-état est ainsi égale à :

${\displaystyle S({\vec {R}})=k_{B}log(\Omega ({\vec {R}}))}$, avec ${\displaystyle k_{B}}$ la constante de Boltzmann.

L'expression ci-dessus donne l'entropie absolue (quantique) du système. Une détermination précise de ${\displaystyle \Omega ({\vec {R}})}$ requerrait un modèle quantique pour la chaîne idéale, ce qui est au-delà de la portée de cet article. Cependant, nous avons déjà calculé la densité de probabilité ${\displaystyle P({\vec {R}})}$ associée au vecteur de la totalité de la chaîne non contrainte, ci-dessus. Comme tous les micro-états de la chaîne idéale ont la même probabilité de s'établir, ${\displaystyle P({\vec {R}})}$ est proportionnel à ${\displaystyle \Omega ({\vec {R}})}$. Ceci nous conduit à l'équation suivante, pour l'expression de l'entropie classique (relative) de la chaîne idéale :

${\displaystyle S({\vec {R}})=k_{B}log(P({\vec {R}}))+cste}$, où ${\displaystyle cste}$ est une constante fixée.

Soit ${\displaystyle {\vec {F}}}$ la force exercée par la chaîne sur le point où son extrémité est attachée. De l'expression précédente de l'entropie, on peut déduire une expression de cette force. Maintenant, supposons qu'au lieu d'être fixée, les positions des deux extrémités de la chaîne idéale soient contrôlées par un opérateur. Les opérateurs contrôlent l'évolution du vecteur entre les extrémités ${\displaystyle {\vec {R}}}$. Si l'opérateur fait varier ${\displaystyle {\vec {R}}}$ d'une petite quantité ${\displaystyle \mathrm {d} {\vec {R}}}$, alors la variation d'énergie interne est nulle, puisque l'énergie de la chaîne est constante. Cette condition peut être écrite comme suit :

${\displaystyle 0=\mathrm {d} U=\delta W+\delta Q}$.

${\displaystyle \delta W}$ est définie comme la quantité élémentaire de travail mécanique transférée par l'opérateur à la chaîne idéale, et ${\displaystyle \delta Q}$ est définie comme la quantité élémentaire de chaleur transférée par le solvant à la chaîne idéale. Maintenant, si l'on suppose que la transformation imposée par l'opérateur est quasi-statique (c'est-à-dire infiniment lente), alors la transformation du système sera réversible par rapport au temps, et l'on peut supposer que durant le passage d'un macro-état ${\displaystyle {\vec {R}}}$ à un macro-état ${\displaystyle {\vec {R}}+\mathrm {d} {\vec {R}}}$, le système passe par une série de macro-états à l'équilibre thermodynamique. Ceci a deux conséquences :

• la quantité de chaleur reçue par le système pendant la transformation peut être reliée à la variation de son entropie :

${\displaystyle \delta Q=T\mathrm {d} S}$, où ${\displaystyle T}$ est la température de la chaîne ;

• afin que la transformation reste infiniment lente, la moyenne de la force exercée par l'opérateur sur les extrémités de la chaîne doit équilibrer la moyenne de la force exercée par la force sur ses propres extrémités. Soit ${\displaystyle {\vec {f}}_{op}}$ la force exercée par l'opérateur et ${\displaystyle {\vec {f}}}$ la force exercée par la chaîne, alors on obtient :

${\displaystyle \delta W=\langle {\vec {f}}_{op}\rangle .\mathrm {d} {\vec {R}}=-\langle {\vec {f}}\rangle .\mathrm {d} {\vec {R}}}$.

Il en vient que :

${\displaystyle \langle {\vec {f}}\rangle =T{\frac {\mathrm {d} S}{\mathrm {d} {\vec {R}}}}={\frac {k_{B}T}{P({\vec {R}})}}{\frac {\mathrm {d} P{\vec {R}}}{\mathrm {d} {\vec {R}}}}}$

${\displaystyle \langle {\vec {f}}\rangle =-k_{B}T{\frac {3{\vec {R}}}{Nl^{2}}}}$.

L'équation ci-dessus est l'équation d'état de la chaîne idéale. Comme l'expression dépend du théorème central limite, il est exact uniquement dans le cas limite d'un polymère contenant un nombre élevé d'unités monomères (donc, à la limite thermodynamique). Il est également seulement valide dans le cas d'une distance entre les extrémités de la chaîne faible par rapport à la longueur totale du contour du polymère, où le comportement est semblable à celui d'un ressort hookéen. Le comportement pour de plus grandes plages de force peut être modélisé en utilisant traitement d'ensemble canonique identique à celui de la magnétisation des spins paramagnétiques. Pour les forces arbitraires, la relation force-extension est donnée par la fonction de Langevin ${\displaystyle {\mathcal {L}}}$ :

${\displaystyle {\frac {R}{Nl}}=\coth \left({\frac {fl}{k_{B}T}}\right)-{\frac {k_{B}T}{fl}}={\mathcal {L}}\left({\frac {fl}{k_{B}T}}\right)}$

où l'extension est ${\displaystyle R=\left\vert {\vec {R}}\right\vert }$.

Pour des extensions arbitraires, la relation force-extension peut être approximée comme suit :

${\displaystyle {\frac {fl}{k_{B}T}}={\mathcal {L}}^{-1}\left({\frac {R}{Nl}}\right)\approx 3{\frac {R}{Nl}}+{\frac {1}{5}}\left({\frac {R}{Nl}}\right)^{2}+\sin \left({\frac {7R}{2Nl}}\right)+{\frac {\left({\frac {R}{Nl}}\right)^{3}}{1-{\frac {R}{Nl}}}}}$,

où ${\displaystyle {\mathcal {L}}^{-1}}$ est l'inverse de la fonction de Langevin, ${\displaystyle N}$^[7] est le nombre de liaisons dans la molécule (donc si la molécule a ${\displaystyle N}$ liaisons, elle est composée de ${\displaystyle N+1}$ unités monomères).

Finalement, le modèle peut être étendu pour des plages de force encore plus larges, en incluant un module d'élasticité le long du contour du polymère, c'est-à-dire en permettant à chaque unité monomère de la chaîne de répondre élastiquement à la force appliquée^[8].

Chaîne idéale sous contrainte de force : ensemble de Gibbs

À travers cette sous-section, comme dans la précédente, les deux extrémités du polymère sont attachées par un dispositif de micro-manipulation. Cette fois, cependant, le dispositif ne maintient pas les deux extrémités de la chaîne dans une position fixe, mais il les maintient plutôt grâce à une force de traction ${\displaystyle {\vec {f}}_{op}}$ constante sur la chaîne idéale. Dans ce cas, les positions des deux extrémités de la chaîne fluctuent autour d'une position moyenne ${\displaystyle \langle {\vec {R}}\rangle }$. La chaîne idéale réagit par une force constante opposée ${\displaystyle {\vec {f}}={\vec {f}}_{op}}$.

Pour une chaîne idéale dans l'ensemble de Gibbs, sous contrainte de force, un macro-état du système est caractérisé par le vecteur ${\displaystyle {\vec {f}}}$.

La différence entre une chaîne idéale sous contrainte de longueur et une chaîne idéale sous contrainte de force est très proche de la différence entre les ensemble microcanonique et ensemble canonique. Le changement consiste à passer d'un état où une valeur fixe est imposée sur un certain paramètre, à un état où le système est laissé libre d'échanger sur ce paramètre avec l'extérieur. Le paramètre en question est l'énergie, pour les descriptions micro-canonique et canonique, alors qu'il s'agit de la longueur de la chaîne dans le cas d'une chaîne idéale.

Comme dans les ensembles micro-canonique et canonique, les deux descriptions de la chaîne idéale diffèrent seulement par la façon dont sont traitées les fluctuations du système. Ainsi, ils sont équivalents à la limite thermodynamique. L'équation d'état de la chaîne idéale reste la même, excepté qu'elle est désormais sujette à des fluctuations.

Références

1. James E. Mark, Physical properties of polymer handbook, Springer, 2006 (ISBN 978-0-387-31235-4, 0387312358 et 9780387690025, OCLC 619279219, lire en ligne)
2. Meyer B. Jackson, Molecular and cellular biophysics, Cambridge University Press, 2006, 512 p. (ISBN 978-0-521-62441-1, 9780521624701 et 052162441X, OCLC 61757063, lire en ligne)
3. Karsten Rippe, « Making contacts on a nucleic acid polymer », Trends in Biochemical Sciences, vol. 26, n^o 12,‎ décembre 2001, p. 733–740 (ISSN 0968-0004, DOI 10.1016/s0968-0004(01)01978-8, lire en ligne, consulté le 16 mai 2018)
4. Martial Mazars, « Statistical physics of the freely jointed chain », Physical Review E, vol. 53, n^o 6,‎ 1^er juin 1996, p. 6297–6319 (DOI 10.1103/PhysRevE.53.6297, lire en ligne, consulté le 16 mai 2018)
5. (en) Jerome Harris Weiner, Statistical mechanics of elasticity, New York/Chichester/Brisbane etc., Wiley, 1983, 439 p. (ISBN 0-471-09773-X et 9780471097730, OCLC 8845964, lire en ligne)
6. (en) Roland G. Winkler, « Deformation of semi-flexible chains », Journal of Chemical Physics,‎ 8 février 2003, p. 2919-2928 (lire en ligne)
7. (en) Michael Rubinstein, Polymer physics, Oxford/New York, Oxford University Press, 2003, 440 p. (ISBN 0-19-852059-X, 9780198520597 et 0198520603, OCLC 50339757, lire en ligne)
8. S. B. Smith, L. Finzi et C. Bustamante, « Direct mechanical measurements of the elasticity of single DNA molecules by using magnetic beads », Science, vol. 258, n^o 5085,‎ 13 novembre 1992, p. 1122–1126 (ISSN 0036-8075, PMID 1439819, lire en ligne, consulté le 15 mai 2018)

• (en) Cet article est partiellement ou en totalité issu de l’article de Wikipédia en anglais intitulé « Ideal chain » (voir la liste des auteurs).

Voir aussi

• modèle du ver élastique (modèle plus sophistiqué que la chaîne idéale)

• Portail de la chimie