Mthodes quantiques Utilisation des mthodes quantiques Dtermination de

  • Slides: 86
Download presentation
Méthodes quantiques

Méthodes quantiques

Utilisation des méthodes quantiques • Détermination de la structure électronique (orbitales, charges atomiques, potentiel

Utilisation des méthodes quantiques • Détermination de la structure électronique (orbitales, charges atomiques, potentiel électrostatique, aspects stéréoélectroniques) • Détermination des énergies • Détermination des géométries • Recherche des états de transition • Effets de solvant

La fonction d’onde: fonction des coordonnées des noyaux et des électrons Ytot( x 1,

La fonction d’onde: fonction des coordonnées des noyaux et des électrons Ytot( x 1, y 1, z 1, ……. . xi, yi, zi, ……… ) Solution de l’équation de Schrödinger Htot. Ytot = EYtot

Passage aux unités atomiques (u. a. ): 1 u. a. = 627. 5 kcal

Passage aux unités atomiques (u. a. ): 1 u. a. = 627. 5 kcal mol-1

1ère approximation: Born-Oppenheimer: Les noyaux sont supposés fixes par rapport aux électrons Ytot =

1ère approximation: Born-Oppenheimer: Les noyaux sont supposés fixes par rapport aux électrons Ytot = Ynucl(coord. noyaux)Yel(coord. elec) Potentiel nucléaire Constante pour un jeu de coordonnées nucléaires Equation de Schrödinger électronique Hel. Yel = Eel. Yel

soit où et est une somme d’opérateurs monoélectroniques est une somme de termes se

soit où et est une somme d’opérateurs monoélectroniques est une somme de termes se rapportant à des paires d’électrons (termes biélectroniques)

Cette équation ne peut être résolue de façon exacte que pour l’atome à un

Cette équation ne peut être résolue de façon exacte que pour l’atome à un électron. Solution: n: nombre quantique principal l ≤ n-1: nombre quantique azimutal -l ≤ m ≤+l: nombre quantique magnétique l =0, 1, 2, 3, 4, ……. sous-couches s, p, d, f, …. .

l est associé au moment cinétique orbital m est associé à Introduction du spin:

l est associé au moment cinétique orbital m est associé à Introduction du spin: L’interprétation de la structure fine du spectre de l’hydrogène nécessite d’introduire un moment cinétique intrinsèque à l’électron, le spin , dont la norme est définie par le nombre quantique ½: La composante Sz de ce moment est msħ ms peut prendre deux valeurs, ± 1/2.

La fonction d’onde de l’électron se met sous la forme d’un produit d’une fonction

La fonction d’onde de l’électron se met sous la forme d’un produit d’une fonction des coordonnées spatiales Fn, l, m(x, y, z) et d’une fonction des « coordonnées » de spin h(ms): est appelée spinorbitale

Système polyélectronique: On ne peut plus obtenir de solution analytique exacte de la fonction

Système polyélectronique: On ne peut plus obtenir de solution analytique exacte de la fonction d’onde 2ème Approximation: Approximation orbitalaire Fonction d’onde totale exprimée comme un produit de fonctions monoélectroniques, c. a. d. dépendant chacune des coordonnées d’un seul électron. Ces fonctions monoélectroniques sont les spinorbitales définies précédemment. Mais ce produit n’est pas un simple produit……….

(i) Les électrons sont indiscernables (ii) Les électrons sont des fermions: la fonction d’onde

(i) Les électrons sont indiscernables (ii) Les électrons sont des fermions: la fonction d’onde Ψ doit être antisymétrique par rapport à l’échange de deux électrons => Principe de Pauli (deux électrons ne peuvent pas avoir tous leurs nombres quantiques égaux) La fonction d’onde est prise comme un produit antisymmétrisé des spinorbitales, le déterminant de Slater. Les spinorbitales sont les orbitales moléculaires du système.

orbitale électron Energie associée au déterminant de Slater:

orbitale électron Energie associée au déterminant de Slater:

l’intégration porte sur un électron quelconque, ici l’électron 1 l’intégration porte sur deux électrons

l’intégration porte sur un électron quelconque, ici l’électron 1 l’intégration porte sur deux électrons quelconques, ici les électrons 1 et 2 Jij = intégrale biélectronique coulombienne Kij = intégrale biélectronique d’échange

On applique alors une méthode de résolution variationnelle: le meilleur choix de y est

On applique alors une méthode de résolution variationnelle: le meilleur choix de y est obtenu en minimisant l’énergie E par rapport à l’ensemble des paramètres définissant les orbitales moléculaires, avec la contrainte que les spinorbitales sont orthonormales, c. a. d. Ce qui amène aux équations de Hartree-Fock (équation aux valeurs propres de l’opérateur de Fock F):

Opérateur de Coulomb Opérateur d’échange

Opérateur de Coulomb Opérateur d’échange

Donc F s’exprime en fonction des OM qu’on recherche D’où la nécessité d’une procédure

Donc F s’exprime en fonction des OM qu’on recherche D’où la nécessité d’une procédure itérative Dans la pratique, on part d’un ensemble d’orbitales moléculaires i d’essai, on construit le déterminant de Slater, on calcule l’énergie correspondante. Cela implique le calcul des intégrales monoélectroniques hii et des intégrales biélectroniques Jij et Kij. On diagonalise F, ce qui fournit un nouvel ensemble d’orbitales i qui sont fonctions propres de F, et on répète le processus jusqu’à convergence.

Résolution des équations de Hartree-Fock: - soit numérique (possible uniquement avec les atomes) -

Résolution des équations de Hartree-Fock: - soit numérique (possible uniquement avec les atomes) - soit analytique Résolution analytique en utilisant une 3ème Approximation: Approximation LCAO-MO: Chaque orbitale moléculaire s’exprime comme une combinaison linéaire de fonctions de base cp Ce sont les coefficients cpi qui sont les paramètres variationnels du calcul HF

1ère approximation: Born-Oppenheimer Equation de Schrödinger électronique Hel. Yel = Eel. Yel 2ème Approximation:

1ère approximation: Born-Oppenheimer Equation de Schrödinger électronique Hel. Yel = Eel. Yel 2ème Approximation: Approximation orbitalaire = déterminant de Slater Méthode de résolution variationnelle Equations de Hartree-Fock 3ème Approximation: Approximation LCAO-MO

Approximation LCAO-MO: Chaque orbitale moléculaire s’exprime comme une combinaison linéaire de fonctions de base

Approximation LCAO-MO: Chaque orbitale moléculaire s’exprime comme une combinaison linéaire de fonctions de base cp Ce sont les coefficients cpi qui sont les paramètres variationnels du calcul HF Dans les programmes standards les fonctions de base cp sont des fonctions centrées sur les atomes.

Choix des fonctions de bases 3 considérations à prendre en compte: (i) La nature

Choix des fonctions de bases 3 considérations à prendre en compte: (i) La nature des fonctions doit avoir un sens d’un point de vue physique; (ii) Il faut choisir une forme de fonctions qui permet un calcul efficace des intégrales qui apparaissent dans les équations HF; (iii) Il faut garder un nombre de fonctions minimum (à cause de la progression de la longueur du calcul en N 4);

(i) Choix naturel: fonctions de Slater du type (ii) Cependant, les fonctions de Slater

(i) Choix naturel: fonctions de Slater du type (ii) Cependant, les fonctions de Slater ne permettent pas un calcul aisé des intégrales, d’où l’idée de les remplacer par des fonctions gaussiennes de type en coordonnées polaires ou en coordonnées cartésiennes

(iii) Le choix de fonctions gaussiennes conduirait, pour avoir des résultats pertinents, à un

(iii) Le choix de fonctions gaussiennes conduirait, pour avoir des résultats pertinents, à un nombre total de fonctions très élevé, et par conséquent à un nombre d’intégrales également très élevé, d’où un coût de calcul rédhibitoire. On a donc été conduit à regrouper ces fonctions de base dites « primitives » en groupes de fonctions dites « contractées » , chaque fonction contractée étant une combinaison linéaire de fonctions primitives. Les cp intervenant dans les orbitales moléculaires sont donc les fonctions contractées

On aboutit donc à la hiérarchisation suivante: • gaussienne « primitive » : où

On aboutit donc à la hiérarchisation suivante: • gaussienne « primitive » : où lx + ly +lz détermine le type de fonction: lx + ly +lz = 0: fonction de type s lx + ly +lz = 1: fonction de type p lx + ly +lz = 2: fonction de type d etc… • fonction de base « contractée » : • orbitale moléculaire construite sur la base des fonctions contractées

La flexibilité de la base peut être augmentée par des fonctions de polarisation et

La flexibilité de la base peut être augmentée par des fonctions de polarisation et des fonctions diffuses. Fonctions de polarisation: fonctions ayant un nombre quantique azimutal supérieur à celui des orbitales de valence, c. a. d. des orbitales p pour H, d pour C, N, O, f pour les métaux de transition. Fonctions diffuses: fonctions ayant un exposant très petit. Servent à décrire les régions loin du noyau, nécessaires pour décrire des anions.

1ère approximation: Born-Oppenheimer Equation de Schrödinger électronique Hel. Yel = Eel. Yel 2ème Approximation:

1ère approximation: Born-Oppenheimer Equation de Schrödinger électronique Hel. Yel = Eel. Yel 2ème Approximation: Approximation orbitalaire = déterminant de Slater Méthode de résolution variationnelle Equations de Hartree-Fock 3ème Approximation: Approximation LCAO-MO

Choix de la base Choix de la géométrie Choix des OM d’essai → Résolution

Choix de la base Choix de la géométrie Choix des OM d’essai → Résolution des équations de Hartree-Fock (n) E , suffisamment semblable à E(n-1) – E(n) ≤ seuil oui Fin du calcul, convergence non

L’analyse de population selon Mulliken

L’analyse de population selon Mulliken

La densité électronique (c. a. d. la probabilité de trouver un électron) d’une spinorbitale

La densité électronique (c. a. d. la probabilité de trouver un électron) d’une spinorbitale moléculaire à une position est donnée par le carré de cette spinorbitale or

Si on somme sur toutes les spinorbitales moléculaires occupées et on intègre sur tout

Si on somme sur toutes les spinorbitales moléculaires occupées et on intègre sur tout l’espace, on obtient le nombre total N d’électrons On généralise en introduisant un nombre d’occupation ni pour chaque spinorbitale moléculaire et en sommant sur toutes les spinorbitales

La matrice D dont les éléments sont Dpq est appelée matrice densité Dans l’analyse

La matrice D dont les éléments sont Dpq est appelée matrice densité Dans l’analyse de population selon Mulliken on utilise la matrice produit (D. S) pour déterminer les contributions atomiques des électrons. On a des contributions spécifiques à chaque fonction de base atomique, et des contributions partagées entre deux fonctions de base.

(i) Elément diagonal: (ii) Elément non diagonal: C’est la moitié du nombre d’électrons partagés

(i) Elément diagonal: (ii) Elément non diagonal: C’est la moitié du nombre d’électrons partagés entre cp et cq

On somme ensuite pour chaque atome l’ensemble des contributions à cet atome. Cela exige

On somme ensuite pour chaque atome l’ensemble des contributions à cet atome. Cela exige d’imposer un schéma de partition des contributions extradiagonales impliquant des fonctions de base localisées sur des atomes différents. Mulliken a proposé une équipartition entre les deux atomes impliqués. D’où: charge électronique de l’atome A: charge de l’atome:

Population de recouvrement entre deux atomes A et B: L’analyse de population selon Mulliken

Population de recouvrement entre deux atomes A et B: L’analyse de population selon Mulliken a des limites dues à la partition arbitraire entre les contributions à deux atomes différents

Inconvénients de l’analyse de population suivant Mulliken • Éléments diagonaux peuvent être supérieurs à

Inconvénients de l’analyse de population suivant Mulliken • Éléments diagonaux peuvent être supérieurs à 2! • Éléments non diagonaux peuvent devenir <0 ! • Pas de justification pour le partage égal des contributions non diagonales entre deux atomes. Les plus électronégatifs devraient recevoir une contribution plus importante. • Lorsqu’on a de très petits exposants sur des fonctions centrées sur un atome, ces fonctions vont décrire la fonction d’onde loin de cet atome. La partition de Mulliken ré-attribue la densité correspondante à l’atome, ce qui peut conduire à des effets bizarres. N’utiliser cette analyse de population que pour des comparaisons.

Le potentiel électrostatique

Le potentiel électrostatique

q 0 • • q + Energie potentielle électrostatique à de l’origine et à

q 0 • • q + Energie potentielle électrostatique à de l’origine et à de la charge q: En u. a. : Peut s’écrire: est le potentiel créé en Où • + • q 1 • q 2 par la charge q située en

Cas d’une distribution continue de charges: Cas d’une molécule: - contribution des noyaux: -

Cas d’une distribution continue de charges: Cas d’une molécule: - contribution des noyaux: - contribution de la densité électronique:

Analyse de population basée sur le potentiel électrostatique On calcule le potentiel électrostatique de

Analyse de population basée sur le potentiel électrostatique On calcule le potentiel électrostatique de la molécule à partir de sa fonction d’onde, et on détermine les charges par ajustement entre ce potentiel électrostatique et un potentiel électrostatique obtenu à partir des charges, avec la contrainte que la somme des charges est égale à la charge totale de la molécule. Des contraintes supplémentaires peuvent être incluses, par exemple la reproduction du moment dipolaire.

Approximation de Born-Oppenheimer Méthodes de type fonction d’onde E=F[Y(r)] Equation de Schrödinger électronique Hel.

Approximation de Born-Oppenheimer Méthodes de type fonction d’onde E=F[Y(r)] Equation de Schrödinger électronique Hel. Yel = Eel. Yel Approximation orbitalaire: = déterminant de Slater Méthode de résolution variationnelle Equations de Hartree-Fock Approximation LCAO-MO Approximations supplémentaires Méthodes semi-empiriques (EHT, MNDO, AM 1, PM 3) Addition de déterminants supplémentaires Méthodes CI, MC-SCF, MPn, CC Vers la solution exacte

Au-delà de Hartree Fock

Au-delà de Hartree Fock

La corrélation électronique Dans la méthode HF, où la fonction d’onde est un déterminant

La corrélation électronique Dans la méthode HF, où la fonction d’onde est un déterminant de Slater, chaque électron voit tous les autres comme un champ approché. La méthode HF ne rend pas compte du mouvement instantané des électrons => HF ne rend pas compte de la corrélation électronique, notamment pour les électrons de spins opposés Question: comment modifier la fonction d’onde HF pour obtenir une énergie plus basse? Les méthodes utilisées reposent sur le principe de la superposition de configurations

a 0 >>> a 1, a 2, a 3, …. : corrélation dynamique a

a 0 >>> a 1, a 2, a 3, …. : corrélation dynamique a 0 ≈ a 1, a 2, a 3, ……: corrélation non dynamique

Les méthodes de résolution • méthodes CI • méthodes MC-SCF • méthodes perturbationnelles: MPn

Les méthodes de résolution • méthodes CI • méthodes MC-SCF • méthodes perturbationnelles: MPn • méthodes « Coupled Cluster »

Les méthodes CI et MC-SCF sont des méthodes variationnelles: on est sûr d’obtenir une

Les méthodes CI et MC-SCF sont des méthodes variationnelles: on est sûr d’obtenir une limite supérieure de l’énergie exacte. Par contre avec les méthodes MPn, rien ne garantit que l’on obtienne une borne supérieure de l’énergie de corrélation. De plus la convergence, en termes d’ordre de perturbation, est aléatoire. Néanmoins la méthode MP 2 et plus généralement les méthodes MPn restent une méthode de choix, notamment lorsqu’on doit traiter des interactions faibles. Pour la description de systèmes standard elle a cependant été « détrônée » par les méthodes DFT

Méthodes « Coupled Cluster » où et agissant sur la fonction d’onde HF génère

Méthodes « Coupled Cluster » où et agissant sur la fonction d’onde HF génère tous les déterminants ayant i excitations

On peut en théorie inclure toutes les excitations d’un type donné (i=1, 2, 3,

On peut en théorie inclure toutes les excitations d’un type donné (i=1, 2, 3, …) jusqu’à un ordre infini du développement de Taylor. En fait, on doit se limiter et la troncature s’effectue au niveau de l’opérateur T. Exemples: CCSD: CCSDT: L’effet de T 3 peut être évalué par une méthode de perturbation: CCSD(T)

Approximation de Born-Oppenheimer Méthodes de type fonction d’onde E=F[Y(r)] Equation de Schrödinger électronique Hel.

Approximation de Born-Oppenheimer Méthodes de type fonction d’onde E=F[Y(r)] Equation de Schrödinger électronique Hel. Yel = Eel. Yel Approximation orbitalaire: = déterminant de Slater Méthode de résolution variationnelle Equations de Hartree-Fock Approximation LCAO-MO Approximations supplémentaires Méthodes semi-empiriques (EHT, MNDO, AM 1, PM 3) Addition de déterminants supplémentaires Méthodes CI, MC-SCF, MPn, CC Vers la solution exacte

Méthodes de la fonctionnelle de la densité E=F[r(r)]

Méthodes de la fonctionnelle de la densité E=F[r(r)]

E=F[r(r)] Idée sous-jacente: partir d’une observable physique pour déterminer l’énergie Théorème d’existence de Hohenberg

E=F[r(r)] Idée sous-jacente: partir d’une observable physique pour déterminer l’énergie Théorème d’existence de Hohenberg et Kohn (1964): l’énergie d’un système, dans son état fondamental est complètement déterminée par la densité électronique r(r) Cette densité électronique obéit aussi au principe variationnel Mais……………deux problèmes (i) Comment choisir de façon rationnelle des densités toujours meilleures? (ii) On ne connaît pas la forme analytique reliant l’énergie à la densité.

Décomposition de l’énergie: E[r] = T[r] + Ene[r] +Eee[r] T: énergie cinétique des électrons

Décomposition de l’énergie: E[r] = T[r] + Ene[r] +Eee[r] T: énergie cinétique des électrons Ene: énergie d’attraction électrons-noyaux Eee: énergie de répulsion électrons-électrons.

E[r] = Tni[r] + Ene[r] + J[r] + Exc[r] DJ[r], K[r] et DT[r] sont

E[r] = Tni[r] + Ene[r] + J[r] + Exc[r] DJ[r], K[r] et DT[r] sont inclus dans Exc[r] L’artifice du système fictif d’électrons non-interagissants, mais ayant la même densité que le système réel permet de réintroduire un déterminant de Slater composé de spinorbitales i appelées orbitales de Kohn-Sham. Si on sait définir Exc, on a simplement un problème variationnel à résoudre: comme dans la méthode HF on minimise E via des équations de type HF. Comme dans la méthode HF on développe les orbitales de Kohn-Sham sur une base de fonctions.

Les méthodes semi-empiriques

Les méthodes semi-empiriques

Approximations communes • Seuls les électrons de valence sont considérés • Base minimale, décrite

Approximations communes • Seuls les électrons de valence sont considérés • Base minimale, décrite par des fonctions de Slater • Approximation ZDO (Zero Differential Overlap): (On néglige tous les produits de fonctions de base dépendant de la même coordonnée électronique si ces fonctions sont localisées sur des atomes différents A et B)

Du fait de cette approximation, certaines intégrales mono et bi électroniques, en particulier les

Du fait de cette approximation, certaines intégrales mono et bi électroniques, en particulier les intégrales biélectroniques à 3 et 4 centres sont nulles Nécessité d’introduire des paramètres pour «compenser» ces approximations et obtenir des résultats «raisonables» Les acronymes correspondent à des méthodes différant par la façon dont la paramétrisation est effectuée

MNDO: • Répulsion entre atomes distants de 2 à 3 Å en général trop

MNDO: • Répulsion entre atomes distants de 2 à 3 Å en général trop grande: ØMolécules stériquement encombrées trop instables; ØEnergies d’activation pour rupture ou formation de liaisons trop élevées ØMolécules hypervalentes trop instables ØStructures non classiques (e. g. cation alkyl) trop instables par rapport aux structures classiques ØInteractions faibles pas fiables; on ne peut pas prédire les liaisons hydrogène • Et bien d’autres défauts…. .

AM 1: • Bien meilleure: ØEnergies d’activation pour rupture ou formation de liaisons bien

AM 1: • Bien meilleure: ØEnergies d’activation pour rupture ou formation de liaisons bien améliorées, restent cependant généralement trop élevées, et les géométries sont souvent trop «serrées» ØMolécules hypervalentes améliorées, mais il reste des erreurs plus grandes que sur les autres types de systèmes ØGroupes alkyles en général trop stables ØForce des liaisons hydrogènes prédite à peu près correctement, mais géométrie souvent fausse • restent des défauts, par exemple: ØLiaisons peroxide trop courtes, géométries des composés avec des phosphores, problème conformationnel pour l’éthanol

PM 3: • S’appuie sur un très grand nombre de données expérimentales ØDonne en

PM 3: • S’appuie sur un très grand nombre de données expérimentales ØDonne en général de meilleures énergies et de meilleures géométries, mais restent des problèmes d’analyse conformationnelle (e. g. éthanol, hydrazine) ØLiaisons hydrogènes encore trop courtes d’environ 0. 1Å ØLiaisons de van der Waals mal décrites ØLiaisons Si-X (X=Cl, Br, I) sous-estimées ØDifficultés patentes à décrire l’azote, en général obtenu trop pyramidalisé ØRotation peptidique (rotation autour de la liaison C-N dans les amides) sous-estimée ( parfois de 15 kcal/mol)

Prudence quand on utilise ces méthodes, les erreurs sur les énergies étant en général

Prudence quand on utilise ces méthodes, les erreurs sur les énergies étant en général très aléatoires Méthodes semi-empiriques ont les mêmes avantages et inconvénients que les méthodes à champ de force: marchent relativement bien pour des structures proches des structures connues, mais elles sont souvent incapables de prédire des types de structures inconnues, car elles n’ont pas été paramétrisées pour cela. Cependant on peut les utiliser (mais sans garantie!) pour des études de réactivité Aucune possibilité d’amélioration systématique

Optimisation des géométries

Optimisation des géométries

Au voisinage du point d’équilibre: 1 -dimension E 0 x = constante de force

Au voisinage du point d’équilibre: 1 -dimension E 0 x = constante de force • dérivée première de l’énergie est nulle: point stationnaire • énergie minimum => k>0

1 -dimension aux points stationnaire k<0 maximum de l’énergie (fréquence imaginaire, E c. a.

1 -dimension aux points stationnaire k<0 maximum de l’énergie (fréquence imaginaire, E c. a. d. nombre complexe imaginaire) k>0 minimum de l’énergie x = coordonnée de réaction A + B-C [A…. . B…. . C]‡ A-B + C

2 -dimensions équivalent à:

2 -dimensions équivalent à:

N-dimensions On peut développer l’énergie en une série de Taylor: : déplacement cartésien à

N-dimensions On peut développer l’énergie en une série de Taylor: : déplacement cartésien à partir du point de référence Pour un point stationnaire ; le gradient de l’énergie est nul = Matrice des constantes de forces = matrice hessienne

On peut réécrire E sous la forme: soit On a alors soit La matrice

On peut réécrire E sous la forme: soit On a alors soit La matrice C est la matrice des vecteurs propres de H ayant les l comme valeurs propres : modes normaux de vibration

N-dimensions i = 1, 2, ……… 3 N-6 Points stationnaires: i = 1, 2,

N-dimensions i = 1, 2, ……… 3 N-6 Points stationnaires: i = 1, 2, ……… 3 N-6, c. a. d. grad E = 0 diagonalisati on Qi li > 0 pour tous les i => minimum d’énergie lq < 0 , li > 0 pour tous les i différents de q => état de transition Qk : chemin de réaction à l’état de transition

Thermochimie

Thermochimie

Les effets de solvant

Les effets de solvant

Plusieurs modèles peuvent être utilisés: (i) Modèle discret Molécule de solvant Molécule de soluté

Plusieurs modèles peuvent être utilisés: (i) Modèle discret Molécule de solvant Molécule de soluté Inconvénients: • approche rapidement prohibitive en termes de durée de calcul; • incertitude quant à la pertinence de la géométrie optimisée

(ii) Modèle du continuum Le soluté est considéré comme une distribution de charges cavité

(ii) Modèle du continuum Le soluté est considéré comme une distribution de charges cavité électriques placée dans une cavité immergée dans un milieu continu soluté de constante diélectrique e Cette distribution polarise le milieu, => induit un champ électrique R dit de réaction, qui en retour interagit avec la distribution de charges du soluté, ce qui conduit à une stabilisation. R dépend de e, de la forme et de la dimension cavité, et de la distribution de charges (ce qui conduit à un processus itératif sur cette distribution de charges)

La distribution de charge du soluté peut être décrite par un dipôle moléculaire, ou

La distribution de charge du soluté peut être décrite par un dipôle moléculaire, ou par une distribution multipolaire. Dans le cas du modèle PCM (Polarized Continuum Model), la polarisation du diélectrique par le soluté est représentée par la création d’une distribution de charges apparentes sur la surface de la cavité. Inconvénients des méthodes du continuum: • l’énergie de solvatation dépend des paramètres de la cavité (taille, forme) • problèmes possibles d’adéquation de la géométrie du soluté à la géométrie de la cavité => utilisation d’une cavité obtenue à partir des sphères de Van der Waals des atomes du soluté se généralise

(iii) Modèle du semi-continuum (ou mixte) Molécule de la première couche de solvatation continuum

(iii) Modèle du semi-continuum (ou mixte) Molécule de la première couche de solvatation continuum cavité soluté Inconvénients: taille du calcul, optimisation de géométrie peut être délicate