Unified microscopic theory of equilibrium thermodynamics and ion association in aqueous and non-aqueous electrolytes with explicit hard-core size

Résumé Technique : Théorie Microscopique Unifiée de la Thermodynamique d'Équilibre et de l'Association Ionique

Énoncé du Problème
La caractérisation de la thermodynamique d'équilibre dans les solutions chargées repose historiquement sur la théorie de Debye-Hückel (DH), une approche de couplage faible (WC) qui néglige la taille de cœur dur (HC) des ions. Bien que la théorie DH soit précise pour les électrolytes monovalents dilués, elle échoue à décrire les solutions condensées et les systèmes fortement chargés où le compactage ionique et les corrélations à courte portée dominent. À l'inverse, bien que les simulations de Monte-Carlo (MC) et les théories numériquement intensives (par exemple, la théorie des fonctionnelles de densité et les théories d'équations intégrales) puissent traiter ces régimes, elles manquent souvent de tractabilité analytique. Les approches antérieures par intégrale de chemin, incluant la théorie CCDH (cumulant-corrected DH) développée par les auteurs, ont intégré avec succès les interactions HC mais ont conservé une clôture gaussienne pour les identités de Schwinger-Dyson (SD). Cette clôture nécessitait un traitement WC des interactions à courte portée (électrostatiques et HC) fortement couplées, limitant la précision des prédictions pour l'énergie excédentaire, les longueurs de criblage et l'association ionique dans les régimes concentrés ou fortement couplés.

Méthodologie
Les auteurs formulent une théorie unifiée basée sur un formalisme d'intégrale de chemin qui traite la charge ionique et les interactions HC sur un pied d'égalité. Le cœur de la méthodologie est le formalisme de l'Auto-Cohérence de Debye-Hückel (SCDH), une version améliorée de l'approche CCDH. Les composantes méthodologiques clés incluent :

1. Schéma de Décomposition : Le potentiel d'interaction de Coulomb est divisé en une composante à courte portée ($v_s$) et une composante à longue portée ($v_l$) en utilisant une longueur de séparation variationnelle $\sigma$. Cela permet un traitement asymétrique des échelles d'interaction : les interactions à longue portée sont gérées via une approche gaussienne WC, tandis que les interactions à courte portée sont traitées via une approximation virielle.
2. Détermination Variationnelle : La longueur de séparation $\sigma$ est déterminée par voie variationnelle en imposant l'invariance du potentiel grand canonique sous les variations de $\sigma$, conduisant à une identité formellement exacte impliquant les fonctions de corrélation à deux points (2PCF).
3. Amélioration Analytique : Une mise à niveau computationnelle significative implique l'évaluation analytique de la fonction de corrélation de niveau gaussien ($G_{l,0}$). Cela réduit le nombre de transformées de Fourier numériques requises, améliorant considérablement la vitesse d'exécution et la convergence numérique.
4. Modèle d'Association Ionique : La théorie incorpore un modèle de mise en paire d'ions de Bjerrum amélioré. La constante d'association est dérivée des fonctions de distribution de paires du SCDH plutôt que du potentiel de Coulomb pur, permettant au modèle de rendre compte de l'atmosphère ionique entourant la paire. La fraction de paires est calculée à l'aide d'une intégrale de configuration avec une distance de coupure déterminée par le minimum de la fonction de distribution radiale.
5. Relations Thermodynamiques : Le formalisme dérive des identités exactes reliant la 2PCF à la fonction de corrélation totale et aux quantités thermodynamiques (énergie excédentaire, pression, coefficient osmotique) via des conditions de moments et les équations de Schwinger-Dyson.

Résultats Clés
La théorie a été validée par rapport à un large éventail de données de simulations MC et de mesures expérimentales :

• Distributions de Paires : Pour les ions monovalents avec des tailles d'hydratation typiques ($d = 3,0$ Å et $4,0$ Å), la théorie SCDH prédit des distributions de paires avec une précision quantitative jusqu'à des concentrations molaires ($n_i \approx 2,0$ M). Cela représente une amélioration significative par rapport à la théorie DH, qui échoue même dans les régimes dilués lorsque la taille des ions est non négligeable.
• Fonctions Thermodynamiques : Dans le régime de fraction de compactage ionique $\eta \lesssim 0,1$, le formalisme capture avec précision la dépendance de la taille des ions pour l'énergie excédentaire et la pression pour des tailles d'ions allant de $d = 14,3$ Å jusqu'à $d = 1,6$ Å. La théorie reproduit l'effet de "sur-criblage" ($\kappa > \kappa_{DH}$) observé dans les simulations.
• Limitations : La théorie montre une précision réduite pour les tailles d'ions extrêmement petites ($d \approx 1,03$ Å) ou équivalemment pour les basses températures ($T \lesssim 102$ K), où l'approximation des cumulants surestime le couplage des charges de même signe. Les auteurs notent que l'accès au point critique liquide-vapeur dans ce régime nécessiterait d'étendre le calcul à des ordres de cumulants plus élevés.
• Coefficients Osmotiques : En ajustant le rayon d'hydratation, la théorie reproduit la dépendance non monotone de la concentration des coefficients osmotiques pour divers électrolytes aqueux et non aqueux. Elle identifie correctement la transition d'un régime dominé par la pression électrostatique attractive (submolaire) vers un régime dominé par la pression HC répulsive (molaire).
• Association Ionique dans les Solvants Non Aqueux : La théorie élucide le mécanisme d'association ionique dans les solvants de faible polarité. Elle prédit que l'association substantielle des ions est confinée au régime submolaire ($n_i \lesssim 50$ mM). À des concentrations plus élevées, la répulsion HC entrave l'association, menant à une quasi-saturation de la fraction de paires. Ce comportement est piloté par la compétition entre l'attraction forte des charges opposées (due à la faible constante diélectrique) et les contraintes de volume exclu.

Signification et Revendications
L'article prétend fournir une théorie microscopique unifiée qui comble le fossé entre la tractabilité analytique et la précision requise pour les électrolytes concentrés. En évitant le traitement de couplage faible des interactions à courte portée grâce à une technique de séparation variationnelle, le formalisme SCDH offre une alternative computationnellement efficace aux simulations MC pour prédire les propriétés thermodynamiques et les distributions de paires dans des systèmes avec des tailles d'ions explicites.

Les auteurs soulignent que leur approche explique avec succès le "sous-criblage" des interactions macromoléculaires dans les électrolytes non aqueux en le liant au régime de concentration spécifique où l'association ionique se produit. La théorie reproduit quantitativement les coefficients osmotiques expérimentaux à travers diverses polarités de solvants et tailles d'ions, validant l'image physique d'une compétition abrupte entre l'attraction électrostatique et la répulsion de volume exclu. Ce travail établit un cadre où l'association ionique n'est pas seulement un phénomène de basse concentration, mais est strictement limitée par les contraintes stériques à des densités plus élevées, une conclusion cohérente avec les expériences de microscopie à force atomique (AFM).