Dans la cellule des organismes vivants, le solvant (l'eau) joue un rôle très important dans la stabilisation des structures tridimensionnelles des macromolécules biologiques et lors de leurs interactions. Les méthodes théoriques de simulations de modélisation moléculaire permettent de compléter les informations partielles sur l'hydratation des biomolécules obtenues par les méthodes expérimentales. Nous avons développé un nouveau modèle de solvatation semi-implicite pour représenter le solvant en modélisation moléculaire. Ce modèle décrit le solvant comme des particules microscopiques dont les propriétés diélectriques découlent des lois macroscopiques de l'électrostatique. Nous obtenons ainsi à l'équilibre électrostatique un fluide de particules de Lennard-Jones non polaires, polarisables par le champ électrique créé par le soluté. Ce modèle a l'intérêt de prendre en compte la structure moléculaire du solvant tout en calculant efficacement l'énergie libre électrostatique de solvatation du système. De plus, il est d'un faible coût numérique comparé aux méthodes explicites. Après avoir implémenté notre modèle dans un programme de dynamique moléculaire et l'avoir paramétré de façon simple, nous l'avons appliqué à plusieurs peptides, protéines et acides nucléiques (ADN et ARN de transfert). Les trajectoires de ces simulations sont stables sur une à deux nanosecondes, et les structures obtenues sont tout à fait en accord avec les méthodes expérimentales et les méthodes théoriques de solvatation explicites. Notre modèle permet également de retrouver les sites préférentiels d'hydratation des molécules étudiées identifiés expérimentalement ou théoriquement, malgré l'absence de liaisons hydrogène dans notre solvant. De plus, nous observons de bonnes corrélations entre les énergies libres électrostatiques de solvatation calculées avec notre modèle et celles calculées avec les méthodes de résolution de l'équation de Poisson-Boltzmann, et ces résultats paraissent très encourageants.

Après la découverte des rayonnements ionisants et leurs effets chimiques, il était important d'étudier et de comprendre les mécanismes de formations des radicaux libres et des produits moléculaires caractérisés par leurs courtes durées de vie. Ceci a encouragé les groupes de recherches à développer leurs outils pour qu'ils puissent réaliser ces études. De nos jours la radiolyse impulsionnelle se manifeste comme un outil fondamental permettant de sonder les effets chimiques ainsi que les mécanismes réactionnels dans le milieu étudié.Le laboratoire de Chimie Physique d'Orsay " LCP " est un laboratoire interdisciplinaire abritant la plateforme " ELYSE " qui est un centre de cinétiques rapides. Grâce au laser femtoseconde et à l'accélérateur d'électrons picoseconde, nous avons eu la possibilité, dans le domaine de la radiolyse, de remonter en temps, en étudiant les effets chimiques dans un milieu réactionnel, jusqu'à ~ 5 ps.Nous nous sommes intéressés par les réactions primaires induites par les rayonnements ionisants en solution et ELYSE représentait l'outil principal pour ces études. Les résultats obtenus concernent :- La détermination directe du rendement radiolytique du radical hydroxyle " HO* " en fonction du temps à l'échelle de la picoseconde ;- Etude de l'effet direct du rayonnement ionisant sur les solutions aqueuses concentrées ainsi que la vérification de la réaction de transfert d'électron ultrarapide entre le soluté et le trou positif " H2O*+ " issu lors de la radiolyse de l'eau ;- Etude à température ambiante de la réaction de transfert d'électron entre un électron solvaté (donneur d'électron) et un soluté organique (accepteur d'électron) en milieu visqueux ;- Etude à température ambiante de la solvatation de l'électron dans l'éthylène glycol et dans le propan-2-ol.

Différentes familles de liquides ioniques ont été sélectionnées pour leur capacité d'absorption de gaz et plus particulièrement de dioxyde de carbone. L'objectif de ces travaux est de définir les liquides ioniques les plus à même d'être utilisés en tant qu'absorbants dits alternatifs dans les procédés de captage de gaz, l'enjeu principal étant à terme de diminuer les coûts de production des procédés industriels. Les systèmes retenus sont des liquides ioniques purs issus de la combinaison de cations imidazolium, pyrrolidinium et phosphonium avec des anions de type carboxylate ou tris (pentafluoroethyl) trifluorophosphate ainsi que des mélanges binaires de liquides ioniques avec anion carboxylate + eau. L'absorption de différents gaz - dioxyde de carbone, azote, protoxyde d'azote et éthane - dans les liquides ioniques purs ainsi que dans les mélanges liquide ionique +eau a été mesurée dans une gamme de températures comprises entre 303.15 et 353.15 K et pour des pressions proches de l'atmosphérique. La sélectivité des liquides ioniques pour l'absorption de dioxyde de carbone par rapport aux autres gaz a pu être déterminée. La caractérisation de l'absorption de gaz est fonction des systèmes étudiés, les interactions entre les gaz et certains liquides ioniques (ou certains mélanges liquide ionique + eau) étant uniquement de type physique alors que dans d'autres systèmes, l'absorption est le résultat à la fois d'interactions physiques mais également de la présence de réaction chimique.

Face au problème des micropolluants émergents trouvés dans l’eau, l’utilisation de carbones activés est un moyen de réduire cette pollution à la source. Le but de ce travail est de mieux comprendre les mécanismes d’adsorption de certains micropolluants sur des tissus et feutres de carbones activés.Neuf molécules ont été étudiées dont des médicaments : la carbamazépine (CBZ), le diclofénac (DFN), l’ibuprofène (IBP) et l’ofloxacine (OFX), un produit anticorrosion : le benzotriazole (BZT), un perturbateur endocrinien : le bisphénol-A (BPA), deux herbicides : le mécoprop (MCP) et le pentachlorophénol (PCP) et une molécule utilisée comme indicateur de pollution des eaux usées : la caféine (CAF). Les adsorbants ultramicroporeux (tissu KIP1200 et feutre CSV4) et l’adsorbant mésoporeux (tissu BBV 800) (fournis par Dacarb, France) ont été caractérisés par adsorption d'azote à 77K et de CO2 à 273K, titrages acido-basiques (méthode de Boehm), mesure du pHpzc (point isoélectrique). Les cinétiques et isothermes d'adsorption ont été étudiées à 25°C à pH=7,5 dans un tampon phosphate NaHPO4/KH2PO4 (à 0,04M). La concentration résiduelle est analysée par HPLC.Les cinétiques d'adsorption ont été étudiées pour les 9 molécules à différentes concentrations initiales. Le temps pour atteindre l’équilibre d’adsorption dépend du volume des molécules ainsi que de leur affinité avec l’adsorbant. La quantité maximale adsorbable dépend du volume microporeux ainsi que de la surface spécifique de l’adsorbant, la quantité adsorbable est donc plus importante sur le tissu KIP 1200 que sur le feutre CSV 4. La vitesse de diffusion est la plus lente pour les adsorbants possédant un volume microporeux important, le tissu mésoporeux BBV 800 permet donc une adsorption rapide grâce à de plus larges pores qui permettent un accès plus rapide à la porosité.Des analyses en mélanges binaires et multi composés ont alors été réalisées pour connaître les paramètres clés gouvernant les cinétiques d’adsorption. Une compétition existe entre molécules dans certains cas (BZT et MCP par exemple) avec une première phase gouvernée par la cinétique d’adsorption liée à la diffusion dans les pores et la deuxième phase gouvernée par des phénomènes thermodynamiques entre le système soluté/solvant/carbone.Les isothermes d'adsorption ont été réalisés à 3 températures différentes et modélisées par des équations de Langmuir-Freundlich pour tous les micropolluants. Des paramètres thermodynamiques (enthalpie d’adsorption et enthalpie libre) ont alors été calculés et corrélés aux propriétés physico-chimiques des molécules. Une corrélation est mise en évidence entre l’enthalpie libre et la polarisabilité des molécules ainsi que les forces de Van der Waals déterminées avec le logiciel COSMO-RS mettant en évidence l’importance des forces non polaires dans le phénomène d’adsorption. Des mesures par calorimétrie d’adsorption à très faibles quantités adsorbées ont permis de mettre en évidence que l’entropie est le paramètre thermodynamique qui contrôle l’adsorption de molécules (BZT, PCP, CAF et OFX) sur le tissu KIP 1200. De fortes énergies d’interaction ont été mis en évidence entre les molécules (BZT, CAF et OFX) et les sites d’adsorption.Une étude d’adsorption-désorption de N2 et de CO2 sur des tissus KIP 1200 chargés en PCP, BZT, CAF et OFX a permis de mieux localiser le lieu de l’adsorption dans la porosité montrant une adsorption prioritairement dans les ultramicropores puis dans les supermicropores. Il a été montré également par cette méthode et par des mesures thermiques que l’eau est fortement adsorbée dans la porosité. / A lot of studies have revealed that some organic molecules such as pharmaceutical molecules, solvents, pesticides, etc.. are frequently found in water, at concentrations below µg/L, even after treatment at the exhaust of wastewater treatment plants. These molecules are highly toxic when accumulated in environment. One of the possibility for removing these micropollutants is the adsorption on activated carbons. Thus the aim of this work is to better understand the adsorption mechanism of some micropollutants onto activated carbon (ACs) in felt or fabric form.Nine micropollutants were studied, such as some pharmaceuticals: Carbamazepine (CBZ), Diclofenac (DFN), Ibuprofen (IBP), and Ofloxacin (OFX); one anti-corrosion compound : Benzotriazol (BZT); two herbicides : Mecoprop (MCP) and Pentachlorophenol (PCP) and an endocrine disruptor : Bisphenol A (BPA). Adsorption of Caffeine (CAF) which is an anthropic indicator of pollution in waste water, was also studied. The ACs (microporous KIP1200 fabric and CSV4 felt and mesoporous BBV 800 fabric, from Dacarb, France) were characterized by N2 adsorption-desorption at 77 K and CO2 adsorption at 273 K, pHpzc (point of zero charge) measurements and acido-basic titrations (Boehm method). The adsorption kinetics and isotherms were studied at pH 7.4 at 25°C in NaHPO4/KH2PO4 buffered solutions (about 0.04 M) using UV spectrometry and HPLC for the analysis of organic molecules in the remaining solution.Kinetics have been studied for 9 molecules at different initial concentrations. Time to reach adsorption equilibrium depends of the volume of the molecule and its affinity with the activated carbon. The maximum adsorbed quantity depends of the microporous volume and the specific area of the adsorbent, the adsorbed quantity is then bigger for KIP 1200 fabric than for CSV 4 felt. The speed of diffusion is slower for the adsorbent with high microporous volume, the mesoporous BBV 800 fabric gives place to a quick adsorption kinetics thanks to its large pores that gives an easy access to porosity.Binary and multi components kinetics have been done in order to understand key processes that drive kinetics adsorption. Competition between molecules have been shown (for BZT and MCP for example). Adsorption kinetics can be divided into two phases: the first one is driven by pore diffusion and the second one by thermodynamic phenomenon between the solute, the solvent and the AC.The adsorption isotherms of the molecules were studied at 13, 25 and 40°C; and the thermodynamic parameters (isoteric enthalpies and entropies Gibbs free energies) were determined. A correlation between Gibbs free energy and polarizability of molecules as well as Van der Waals energy calculated with Cosmotherm software shows the importance of non polar forces on adsorption phenomenon. Adsorption calorimetry experiments showed that entropy is the thermodynamic parameter that drives adsorption of molecules (BZT, PCP, CAF and OFX) onto KIP 1200 fabric.The pore size distributions of the carbons loaded with micropollutants were determined by DFT simulation from CO2 and N2 adsorption isotherms, to investigate the porosity accessible to the adsorbate. The accessible pore are firstly the ultramicropores and then supermicropores. With this technique and thermal experiments, it has been shown that water is highly bonded inside the porosity.

La théorie de la fonctionnelle de la densité classique est utilisée pour étudier la solvatation de solutés quelconques dans le solvant eau. Une forme approchée de la fonctionnelle d’excès pour l’eau est proposée. Cette fonctionnelle nécessite l’utilisation de fonctions de corrélation du solvant pur. Celles-ci peuvent être calculées par simulations numériques, dynamique moléculaire ou Monte Carlo ou obtenues expérimentalement. La minimisation de cette fonctionnelle donne accès à l’énergie libre de solvatation ainsi qu’à la densité d’équilibre du solvant. Différentes corrections de cette fonctionnelle approchée sont proposées. Une correction permet de renforcer l’ordre tétraédrique du solvant eau autour des solutés chargés, une autre permet de reproduire le comportement hydrophobe à longue distance de solutés apolaires. Pour réaliser la minimisation numérique de la fonctionnelle, la théorie a été implémentée sur une double grille tridimensionnelle pour les coordonnées angulaires et spatiales, dans un code de minimisation fonctionnelle écrit en Fortran moderne, mdft. Ce programme a été utilisé pour étudier la solvatation en milieu aqueux de petits solutés atomiques neutres et chargés et de petites molécules polaires et apolaires ainsi que de solutés plus complexes, une argile hydrophobe et une petite protéine. Dans chacun des cas la théorie de la fonctionnelle de la densité classique permet d’obtenir des résultats similaires à ceux théoriquement exacts obtenus par dynamique moléculaire, avec des temps de calculs inférieurs d’au moins trois ordres de grandeurs. / A classical density functional theory is applied to study solvation of solutes in water. An approx- imate form of the excess functional is proposed for water. This functional requires the knowledge of pure solvent direct correlation functions. Those functions can be computed by using molecular simulations such as molecular dynamic or Monte Carlo. It is also possible to use functions that have been determined experimentally. The functional minimization gives access to the solvation free energy and to the equilibrium solvent density. Some correction to the functional are also proposed to get the proper tetrahedral order of solvent molecules around a charged solute and to reproduce the correct long range hydrophobic behavior of big apolar solutes. To proceed the numerical minimization of the functional, the theory has been discretized on two tridimensional grids, one for the space coordinates, the other for the angular coordinates, in a functional mini- mization code written in modern Fortran, mdft. This program is used to study the solvation in water of small solutes of several kind, atomic and molecular, charged or neutral. More complex solutes, a neutral clay and a small protein have also been studied by functional minimization. In each case the classical density functional theory is able to reproduce the exact results predicted by MD. The computational cost is at least three order of magnitude less than in explicit methods.

Les réactions d’oxydation d’hydrocarbures en phase liquide (aussi appelées auto-oxydation) jouent un rôle essentiel dans un grand nombre de procédés de l’industrie pétrochimique car elles assurent la conversion du pétrole en composés chimiques organiques valorisables. Elles régissent également la stabilité à l’oxydation des carburants (vieillissement) et des produits chimiques dérivés du pétrole. Ces réactions d’oxydation en phase liquide relèvent de mécanismes radicalaires en chaîne impliquant des milliers d’espèces et de réactions élémentaires. La modélisation cinétique de tels systèmes reste actuellement un défi car elle nécessite de disposer de données thermodynamiques et cinétiques précises, qui sont rares dans la littérature. Le logiciel EXGAS, développé au LRGP, permet de générer automatiquement des modèles cinétiques détaillés pour des réactions d’oxydation d’hydrocarbures en phase gazeuse. Qu’il s’agisse d’une phase gazeuse ou liquide, les réactions élémentaires mises en jeu sont de même nature et la méthodologie de génération du mécanisme est la même. Pour passer d’un mécanisme en phase gaz à un mécanisme en phase liquide il convient d’adapter les valeurs des constantes d’équilibre et de vitesse (appelées constantes thermocinétiques) des réactions du mécanisme. L’objectif de cette thèse est de proposer une méthode pour corriger les constantes thermocinétiques de la phase gaz pour qu’elles deviennent applicables à la phase liquide. Cette correction fait intervenir une grandeur appelée énergie de GIBBS de solvatation molaire partielle. Une analyse de la définition précise de cette quantité nous a permis de montrer qu’elle s’exprime simplement en fonction d’un coefficient de fugacité et d’une densité molaire. Nous avons ensuite relié cette grandeur à des quantités thermodynamiques mesurables (coefficients d’activité, constantes de HENRY …) et nous nous sommes appuyés sur toutes les données qu’il nous a été possible de trouver dans la littérature pour créer la banque de données expérimentales d’énergies de GIBBS de solvatation molaires partielles la plus complète (intitulée CompSol). Cette banque de données a ensuite servi à valider l’utilisation de l’équation d’état UMR-PRU pour prédire ces énergies. Les bases d’une équation d’état de type SAFT, au paramétrage original, développé dans le cadre de cette thèse, ont été posées. Notre objectif était de simplifier l’estimation des paramètres corps purs de cette équation d’état en proposant une méthode de paramétrage ne nécessitant aucune procédure d’optimisation, claire et reproductible, à partir de données très facilement accessibles dans la littérature. Cette équation a été utilisée pour estimer les énergies de GIBBS de solvatation molaires des corps purs et les énergies de GIBBS de solvatation molaires partielles de systèmes {soluté+solvant}. Enfin, ces méthodes d’estimation des énergies de GIBBS de solvatation molaires partielles ont été combinées au logiciel EXGAS afin de modéliser l’oxydation du n-butane en phase liquide / Liquid phase oxidation of hydrocarbons (also called autoxidation) is central to a large number of processes in the petrochemical industry as it plays a key role in the conversion of petroleum feedstock into valuable organic chemicals. This phenomenon is also crucial in oxidation-stability studies of fuels and its derivatives (aging). These liquid-phase oxidation reactions entail radical mechanisms involving more than thousands of compounds and elementary reactions. Kinetic modelling of these kinds of reactions remains a significant challenge because it requires thermodynamic and kinetic parameters, which are not abundant in literature. The EXGAS software, developed at LRGP, is able to generate these kinds of models but only for oxidation reactions taking place in a gaseous phase. It is assumed that the nature of elementary reactions in the liquid and gaseous phases is the same. The unique need to transfer a kinetic mechanism from a gas phase to a liquid phase is to update kinetic rate constant values and equilibrium constant values (called thermokinetic constants) of mechanism reactions. Therefore, in the framework of this PhD thesis, a new method aimed at applying a correction term to thermokinetic constants of gaseous phases is proposed in order to obtain constants usable to describe liquid-phase mechanisms. This correction involves a quantity called partial molar solvation GIBBS energy. An analysis of the precise definition of this property led us to conclude that it can be simply expressed as a function of fugacity coefficients and liquid molar density. As a result, this property could also be expressed with respect to measurable thermodynamic quantities as activity coefficients or HENRY’s law constants. By combining all the experimental data related to these measurable properties that can be found in the literature, it was possible to develop a comprehensive databank of partial molar solvation GIBBS energies (called the CompSol database). This database was used to validate the use of the UMR-PRU equation of state to predict solvation quantities. Moreover, the bases of a new parameterization for SAFT-type equations of state were laid. It consists in estimating pure-component parameters of SAFT-like equation using a very simple, reproducible and transparent path for non-associating pure components. This equation was used to calculate partial molar GIBBS energy of solvation of pure and mixed solutes. Last, equations of state were combined with EXGAS software to model the oxidation of n-butane in the liquid phase

Diese Arbeit charakterisiert die Solvatation des organischen Farbstoffes Coumarin 153 in Gemischen aus jeweils einem Alkan und einem Alkohol. Dabei werden Methoden der statischen und zeitaufgeloesten optischen Spektroskopie sowie klassische molekulardynamische Simulationen fuer die Untersuchungen angewendet. Die experimentellen Ergebnisse zeigen, dass der Farbstoff im Gemisch selektiv durch den Alkohol solvatisiert ist. Die Staerke dieser Solvatation ist für den elektronischen Grund- und Anregungszustand des Farbstoffes unterschiedlich. Aus diesem Grund wird die Solvatationsdynamik in Alkohol/Alkan Gemischen durch einen Translationsdiffusionsprozess bestimmt. Die molekulardynamischen Simulationen veranschaulichen die selektive Solvatation des Coumarin 153 in einem Methanol/Hexan Gemisch. Die Solvathuelle enhaelt im Anregungszustand des Farbstoffes bis zu dreimal mehr Molekuele als im Grundzustand. Im Unterschied zur Solvatation in reinem Methanol spielen spezifische Bindungen wie Wasserstoffbrueckenbindungen in einem Methanol/Hexan Gemisch eine wesentliche Rolle.

Les propriétés de solvatation jouent un rôle important dans les problèmes chimiques et biochimiques. La théorie fonctionnelle de la densité moléculaire (MDFT) est l'une des méthodes frontières pour évaluer ces propriétés, dans laquelle une fonction d'énergie libre de solvatation est minimisée pour un soluté arbitraire dans une boîte de solvant cubique périodique. Dans cette thèse, nous travaillons sur l'évaluation du terme d'excès de la fonctionnelle d’énergie libre sous l’approximation du fluide de référence homogène (HRF), équivalent à l'approximation de la chaîne hypernettée (HNC) dans la théorie des équations intégrales. Deux algorithmes sont proposés: le premier est une extension d'un algorithme précédent, qui permet de traiter le cas d'un solvant moléculaire à trois dimensions (en fonction de trois angles d'Euler) au lieu d'un solvant linéaire (selon deux angles); L'autre est un nouvel algorithme qui intègre le traitement de la convolution angulaire de l'équation Ornstein-Zernike (OZ) moléculaire dans MDFT, et en fait développe la densité du solvant et le gradient fonctionnel en harmoniques sphériques généralisées (GSHs). On montre que le nouvel algorithme est beaucoup plus rapide que le précédent. Les deux algorithmes sont appropriés pour des solutés arbitraires tridimensionnel dans l'eau liquide, et pour prédire l'énergie libre et la structure de solvatation d'ions et de molécules. / Solvation properties play an important role in chemical and bio-chemical issues. The molecular density functional theory (MDFT) is one of the frontier numerical methods to evaluate these properties, in which the solvation free energy functional is minimized for an arbitrary solute in a periodic cubic solvent box. In this thesis, we work on the evaluation of the excess term of the free energy functional under the homogeneous reference fluid (HRF) approximation, which is equivalent to hypernetted-chain (HNC) approximation in integral equation theory. Two algorithms are proposed: the first one is an extension of a previously implemented algorithm, which makes it possible to handle full 3D molecular solvent (depending on three Euler angles) instead of linear solvent (depending on two angles); the other one is a new algorithm that integrates the molecular Ornstein-Zernike (OZ) equation treatment of angular convolution into MDFT, which in fact expands the solvent density and the functional gradient on generalized spherical harmonics (GSHs). It is shown that the new algorithm is much more rapid than the previous one. Both algorithms are suitable for arbitrary three-dimensional solute in liquid water, and are able to predict the solvation free energy and structure of ions and molecules.

Les ARN, seules molécules à la fois informatives et catalytiques, interviennent à tous les niveaux de la vie cellulaire. Pour certaines molécules d'ARN, comme les ARN messagers, la fonction réside dans la séquence. Pour d'autres, comme les ARN de transfert, la fonction réside dans l'adoption et le maintien d'une structure tridimensionnelle compacte. Ces derniers, dits ARN structurés, nécessitent des ions Mg2+ pour maintenir leur structure tertiaire et assurer leur fonction de catalyse. Dans l'équipe d'Eric WESTHOF, de nombreux modèles tridimensionnels de ces ARN structurés, élaborés à partir de contraintes fournies par des méthodes de cartographie chimique en solution et de comparaison de séquences, ont été construits. Ces modèles regroupent l'ensemble des données expérimentales existantes pour un système donné et permettent ainsi de proposer des mécanismes de réaction ou de reconnaissance pertinents. Toutefois, ces modèles ne prétendent pas atteindre la résolution atomique de structures cristallographiques. C'est pourquoi des méthodes de simulations de dynamique moléculaire (DM) sont développées afin d'affiner les structures modélisées et, ultérieurement, de mieux comprendre les relations structure-fonction de ces ARN.<br /> Les progrès récents dans les techniques de simulation de DM réalisés au laboratoire, nous permettent maintenant de simuler la dynamique de molécules de la taille d'un ARN de transfert ou de petits ribozymes avec un degré de précision jusqu'alors rarement atteint. Pour mieux comprendre les fonctions des ARN, qui sont par nature des processus dynamiques, nous utilisons les méthodes de DM pour compléter les informations structurales fragmentaires qui sont disponibles sur l'hydratation de ces ARN et le rôle d'agents structurants (ions monovalents, divalents et anions). En effet, cette connaissance préalable est indispensable pour la compréhension du repliement et de la stabilité de ces ARN. De plus, elle contribuera à la construction de modèles structuraux plus précis. Parmi les systèmes que nous nous proposons d'étudier se trouvent les ARN de transfert ainsi que des fragments d'ARN ribosomaux. Des complexes entre ARN et antibiotiques ainsi que des complexes ARN/protéines sont en cours d'étude.

La théorie de la fonctionnelle de la densité classique est utilisée pour étudier la solvatation de solutés quelconques dans le solvant eau. Une forme approchée de la fonctionnelle d'excès pour l'eau est proposée. Cette fonctionnelle nécessite l'utilisation de fonctions de corrélation du solvant pur. Celles-ci peuvent être calculées par simulations numériques, dynamique moléculaire ou Monte Carlo ou obtenues expérimentalement. La minimisation de cette fonctionnelle donne accès à l'énergie libre de solvatation ainsi qu'à la densité d'équilibre du solvant. Différentes corrections de cette fonctionnelle approchée sont proposées. Une correction permet de renforcer l'ordre tétraédrique du solvant eau autour des solutés chargés, une autre permet de reproduire le comportement hydrophobe à longue distance de solutés apolaires. Pour réaliser la minimisation numérique de la fonctionnelle, la théorie a été implémentée sur une double grille tridimensionnelle pour les coordonnées angulaires et spatiales, dans un code de minimisation fonctionnelle écrit en Fortran moderne, mdft. Ce programme a été utilisé pour étudier la solvatation en milieu aqueux de petits solutés atomiques neutres et chargés et de petites molécules polaires et apolaires ainsi que de solutés plus complexes, une argile hydrophobe et une petite protéine. Dans chacun des cas la théorie de la fonctionnelle de la densité classique permet d'obtenir des résultats similaires à ceux théoriquement exacts obtenus par dynamique moléculaire, avec des temps de calculs inférieurs d'au moins trois ordres de grandeurs.