⚗️ Dissolution cristalline

La première loi de Fick énonce que le flux diffusif J est proportionnel au gradient de concentration : J = -D (dC/dx). Dans le modèle du noyau rétrécissant utilisé ici, le gradient à travers la couche limite d'épaisseur delta est approximé comme (Cs - Cbulk) / delta, donnant un flux J = D * (Cs - Cbulk) / delta. Le taux molaire total de dissolution est alors J multiplié par la surface du cristal (2 pi r pour le modèle circulaire 2D montré). À mesure que le rayon r rétrécit, la surface diminue et le taux chute même si tous les autres paramètres restent constants.

Les comprimés pharmaceutiques doivent se dissoudre à un taux contrôlé pour libérer les principes actifs dans la circulation sanguine au bon moment et à la bonne dose — ceci est entièrement gouverné par la diffusion de Fick à travers le milieu de dissolution. Les géologues utilisent la cinétique de dissolution pour modéliser comment les paysages karstiques calcaires se forment sur des milliers d'années. En construction, les particules de clinker de ciment se dissolvent dans l'eau pour former du silicate de calcium hydraté, la phase liante qui donne au béton sa résistance. La production de sucre et de sel dépend toutes deux de cycles optimisés de dissolution et de cristallisation.

L'agitation n'accélère la dissolution que lorsque la diffusion à travers la couche limite est l'étape limitante, ce qui est le cas le plus courant pour les composés peu solubles. Quand la dissolution est limitée par la réaction de surface plutôt que par la diffusion, l'agitation a peu d'effet car le goulot d'étranglement est le détachement des ions du réseau cristallin, pas leur transport loin de la surface. La simulation modélise le régime contrôlé par la diffusion, où l'agitation comprime l'épaisseur de la couche limite et augmente directement le flux de diffusion.

Le modèle du noyau rétrécissant a été formalisé en génie chimique par Octave Levenspiel dans son manuel emblématique de 1962 « Chemical Reaction Engineering ». Le modèle décrit une particule réactive ou en dissolution qui rétrécit de l'extérieur vers l'intérieur à mesure que le noyau non réagi ou non dissous devient plus petit. Des travaux théoriques antérieurs avaient été posés par Adolf Fick, qui publia ses lois de diffusion en 1855 à partir d'expériences sur le sel diffusant dans l'eau, établissant une analogie avec la loi de conduction thermique de Fourier.

La dissolution cristalline est l'inverse de la cristallisation, et les deux sont gouvernées par la même force motrice thermodynamique : la différence entre la concentration réelle et la concentration d'équilibre. Les phénomènes apparentés incluent le mûrissement d'Ostwald, où de petits cristaux se dissolvent et se redéposent sur des plus grands pour réduire l'énergie de surface totale ; les réactions de précipitation, où des ions dissous se recombinent pour former une nouvelle phase solide ; et la lixiviation en hydrométallurgie, où des métaux précieux sont sélectivement dissous à partir de minerais. La simulation est aussi conceptuellement liée à la corrosion et à la dissolution électrochimique des métaux.

Les formulateurs pharmaceutiques conçoivent l'habitus cristallin, la distribution granulométrique, et les matériaux de revêtement pour obtenir des profils de dissolution ciblés — par exemple des enrobages gastro-résistants qui résistent à la dissolution dans l'acide gastrique mais se dissolvent rapidement dans l'intestin. Dans la fabrication de semi-conducteurs, la gravure humide sélective dissout des plans cristallins spécifiques de silicium ou d'arséniure de gallium pour structurer des caractéristiques à l'échelle nanométrique. Les cuves de lixiviation hydrométallurgique utilisent l'agitation et le contrôle de la température pour maximiser le taux de dissolution des minéraux de cuivre ou d'or à partir du minerai broyé, appliquant exactement les mêmes principes de diffusion de Fick montrés dans cette simulation.

Les chercheurs utilisent la microscopie à force atomique in situ (AFM) et des techniques de rayons X synchrotron pour observer des faces cristallines individuelles se dissoudre à l'échelle nanométrique en temps réel, révélant une dynamique de recul de marches et de fossettes de gravure que les modèles globaux comme le noyau rétrécissant ne peuvent pas capturer. Les simulations de dynamique moléculaire résolvent désormais l'énergie de détachement des ions individuels des sites de coude sur la surface du cristal. En science pharmaceutique, des modèles d'apprentissage automatique sont entraînés pour prédire les taux de dissolution intrinsèques à partir de la seule structure cristalline, dans le but d'accélérer le criblage précoce des candidats médicaments sans expériences physiques.