Lanthanides et Actinides : Propriétés et Applications
📝 Situation du Projet : L'Enjeu Stratégique des Aimants Permanents
Vous avez été recruté comme Ingénieur Procédés Senior au sein du pôle Innovation d'un géant de la chimie fine, basé à La Rochelle, leader mondial de la séparation des Terres Rares. Le contexte géopolitique et énergétique mondial exerce une pression sans précédent sur la production de Néodyme (\(\text{Nd}\)). Cet élément est le composant critique des aimants permanents \(\text{NdFeB}\) (Néodyme-Fer-Bore) qui équipent les rotors des éoliennes offshore de nouvelle génération et les moteurs synchrones des véhicules électriques. Une éolienne de 3 MW nécessite à elle seule plus de 600 kg de Néodyme. La pureté exigée par les fabricants d'aimants est drastique : 99.9% minimum, car toute impureté de Praséodyme (\(\text{Pr}\)) dégrade la coercivité magnétique à haute température.
Le défi technique est colossal. Le minerai source, une monazite importée, a subi une attaque sodique puis une dissolution nitrique. La solution obtenue contient un mélange complexe de Lanthanides légers (\(\text{La}, \text{Ce}, \text{Pr}, \text{Nd}, \text{Sm}\)). Les étapes précédentes ont permis d'isoler le couple \(\text{Nd}/\text{Pr}\) des autres éléments. Vous vous trouvez donc face à une solution aqueuse concentrée contenant uniquement ces deux "frères jumeaux" chimiques. Le problème réside dans leur séparation : le Néodyme (\(Z=60\)) et le Praséodyme (\(Z=59\)) ont des rayons ioniques quasi-identiques (à 1 pm près !) et des propriétés chimiques si semblables que les méthodes de précipitation classiques sont inopérantes. Votre mission est de concevoir et dimensionner la batterie d'extraction liquide-liquide (SX) qui réalisera cet exploit, en exploitant la très faible différence de leurs coefficients de partage dans un solvant organophosphoré.
Vous devez dimensionner l'unité de fractionnement final. Cela implique de :
1. Justifier atomistiquement la difficulté de séparation (structure électronique 4f).
2. Calculer les paramètres thermodynamiques (coefficients de partage et facteur de séparation).
3. Déterminer le nombre d'étages théoriques nécessaires pour atteindre la pureté de 99.9%, en utilisant le modèle de Kremser à contre-courant.
"Attention, l'agent extractant utilisé est de type organophosphoré (TBP ou HDEHP) dilué dans du kérosène. Ce sont des composés inflammables et irritants. De plus, la solution aqueuse contient de l'acide nitrique concentré (\(\text{HNO}_3\)). Assurez-vous que vos calculs de débits évitent tout engorgement (flooding) des décanteurs."
Pour mener à bien votre étude de dimensionnement, vous disposez des données thermodynamiques et cinétiques issues des essais en laboratoire. Ces valeurs ont été validées pour l'isotherme d'extraction à 25°C et constituent la base non négociable de vos calculs.
📚 Référentiel Normatif & Scientifique
Les standards suivants régissent la nomenclature et la sécurité des produits manipulés :
IUPAC Red Book (Nomenclature Inorganique) ISO 11014 (Fiche de Données de Sécurité)Le tableau ci-dessous résume les caractéristiques fondamentales des deux éléments. Notez la proximité des numéros atomiques (Z) et l'ordre de grandeur des coefficients de distribution (\(D\)). Le coefficient \(D\) représente le rapport de concentration à l'équilibre : \(D = [C]_{\text{org}} / [C]_{\text{aq}}\). Une valeur supérieure à 1 indique une affinité pour la phase organique.
| DONNÉES ATOMIQUES | |
| Praséodyme (\(\text{Pr}\)) | Z = 59 (Voisin "Léger") |
| Néodyme (\(\text{Nd}\)) | Z = 60 (Voisin "Lourd") |
| Configuration Gaz Rare de Référence | [\(\text{Xe}\)] (Xénon, Z=54) |
| PARAMÈTRES D'EXTRACTION (à l'équilibre) | |
| Coeff. de distribution \(D_{\text{Nd}}\) | 2.54 (Adimensionnel - Affinité Moyenne) |
| Coeff. de distribution \(D_{\text{Pr}}\) | 1.82 (Adimensionnel - Affinité Faible) |
| Pureté cible (\(P_{\text{cible}}\)) | 99.9 % (Qualité "Aimant") |
📐 Paramètres de Dimensionnement du Procédé
Ces paramètres hydrodynamiques définissent les conditions opératoires du pilote. Le rapport de phase (\(\text{O}/\text{A}\)) est fixé par les débits pour optimiser le transfert de masse sans créer d'émulsion stable.
- Débit Phase Aqueuse (\(Q_{\text{aq}}\)): 1000 L/h (Flux entrant contenant les terres rares)
- Débit Phase Organique (\(Q_{\text{org}}\)): 1500 L/h (Solvant extractant frais ou recyclé)
- Concentration Initiale \(\text{Nd}\) (\(x_{F, \text{Nd}}\)): 12 g/L
- Concentration Initiale \(\text{Pr}\) (\(x_{F, \text{Pr}}\)): 4 g/L
⚖️ Contraintes de Fonctionnement
Pour garantir la rentabilité économique du procédé, les pertes en Néodyme dans le raffinat (rejet aqueux) doivent être minimes.
E. Protocole de Résolution
Pour dimensionner cette unité de séparation complexe, nous suivrons une démarche rigoureuse allant de l'atomistique fondamentale jusqu'au génie chimique appliqué.
Analyse Électronique & Contraction
Étude des configurations électroniques des ions \(\text{Nd}^{3+}\) et \(\text{Pr}^{3+}\) pour justifier la difficulté de séparation (similitude chimique) et l'origine de la contraction des lanthanides.
Thermodynamique de l'Extraction
Calcul du facteur de séparation intrinsèque (\(\alpha\)) à partir des coefficients de distribution expérimentaux (\(D\)). C'est la clé de la "sélectivité" du solvant.
Dimensionnement du Procédé (Kremser)
Calcul du Facteur d'Extraction (\(E\)) prenant en compte les débits. Application de l'équation de Kremser pour déterminer le nombre d'étages théoriques (\(N\)) nécessaires en contre-courant.
Validation de la Pureté
Vérification finale des ratios de concentration en sortie pour s'assurer que la pureté cible de 99.9% est techniquement atteignable avec le dimensionnement proposé.
Lanthanides et Actinides : Propriétés et Applications
🎯 Objectif
L'objectif premier de cette étude est d'établir la justification fondamentale de la difficulté de séparation entre le Néodyme et le Praséodyme. En tant qu'ingénieur, nous ne pouvons pas nous contenter de constater que les coefficients de distribution sont proches ; nous devons comprendre pourquoi au niveau atomique. Cela passe par l'écriture rigoureuse des configurations électroniques pour identifier le remplissage de la sous-couche 4f, responsable de l'effet d'écran imparfait et de la "contraction des lanthanides". C'est cette contraction qui induit la minuscule différence de rayon ionique exploitée par notre procédé d'extraction.
📚 Référentiel
Pour déterminer l'état fondamental des atomes et des ions, nous nous basons sur les principes quantiques suivants :
- Principe de stabilité (Aufbau) : Remplissage des orbitales par énergie croissante (règle de Klechkowski).
- Règle de Hund : Maximisation du spin total dans une sous-couche dégénérée.
- Principe d'exclusion de Pauli : Deux électrons ne peuvent avoir les mêmes 4 nombres quantiques.
Les électrons de valence qui régissent la chimie classique (liaisons covalentes, ioniques) sont ceux des couches externes. Or, chez les Lanthanides, le remplissage se fait sur la sous-couche 4f, qui est une couche "interne", enfouie sous les couches pleines 5s et 5p et la couche de valence 6s. Ces électrons 4f sont donc masqués de l'environnement extérieur. C'est pourquoi tous les lanthanides ont un comportement chimique quasi-identique. Notre seule "prise" pour les séparer est la variation de taille : plus le noyau est chargé (Z augmente), plus il attire ces électrons internes, réduisant le rayon. Nous devons vérifier que Nd et Pr sont bien consécutifs et stables à l'état trivalent.
Contrairement aux métaux de transition (bloc d), les électrons 4f ont une forme très diffuse et n'écrantent pas efficacement la charge nucléaire pour les électrons externes. Ainsi, lorsque l'on passe du \(\text{La}\) (\(Z=57\)) au \(\text{Lu}\) (\(Z=71\)), chaque proton ajouté au noyau attire plus fortement le nuage électronique global, causant une diminution progressive et régulière du rayon ionique (de 103 pm à 86 pm). Cette contraction augmente la densité de charge de l'ion, renforçant légèrement ses interactions électrostatiques avec les ligands extractants.
📋 Données d'Entrée
| Donnée | Symbole | Valeur |
|---|---|---|
| Numéro Atomique Praséodyme | \(Z_{\text{Pr}}\) | 59 |
| Numéro Atomique Néodyme | \(Z_{\text{Nd}}\) | 60 |
| Numéro Atomique Xénon (Référence) | \(Z_{\text{Xe}}\) | 54 |
Pour trouver rapidement le nombre d'électrons 'f' d'un lanthanide trivalent, regardez simplement le chiffre des unités de son numéro atomique, puis ajustez. Mais la méthode rigoureuse est la suivante :
Exemple pour Nd(60) :
📝 Calculs Détaillés des Configurations
1. Configuration du Praséodyme (Pr, Z=59)
Détail de la manipulation :
1. Identification du Gaz Rare : On repère le gaz rare précédant le Praséodyme. Il s'agit du Xénon (Z=54).
2. Décompte des Électrons de Valence : On soustrait le numéro atomique du gaz rare à celui de l'élément :
3. Remplissage (Klechkowski) : On remplit d'abord la couche \(6\text{s}\) (2 électrons), puis on place les 3 restants dans la \(4\text{f}\).
4. Ionisation Trivalente : Pour former l'ion \(+3\), on arrache **toujours** les électrons de la couche la plus externe (\(n\) le plus grand) en premier, donc les 2 électrons \(6\text{s}\), puis le dernier électron nécessaire dans la \(4\text{f}\).
Interprétation Post-Calcul : L'ion \(\text{Pr}^{3+}\) possède 2 électrons célibataires dans ses orbitales f.
2. Configuration du Néodyme (Nd, Z=60)
Le Néodyme suit le Praséodyme. Il a un proton et un électron de plus.
Ionisation identique.
Interprétation Post-Calcul : L'ion \(\text{Nd}^{3+}\) possède 3 électrons f.
Les calculs confirment que \(\text{Pr}\) et \(\text{Nd}\) sont des voisins immédiats dans le tableau périodique. La structure électronique externe (celle qui réagit) est identique (\(6\text{s}^0\) pour les ions). La seule différence est un électron profond supplémentaire pour le Nd. Cela valide l'hypothèse de la contraction des lanthanides : le rayon de \(\text{Nd}^{3+}\) est donc très légèrement inférieur à celui de \(\text{Pr}^{3+}\). Cette densité de charge supérieure explique pourquoi \(D_{\text{Nd}} > D_{\text{Pr}}\) (le Nd se lie un peu mieux aux extractants).
Ne confondez pas la configuration de l'atome neutre (utilisée en métallurgie) et celle de l'ion trivalent (chimie en solution). L'erreur classique est d'oublier que les électrons 's' partent en premier lors de l'ionisation !
🎯 Objectif
Après avoir validé l'origine atomique de la différence de comportement, nous devons maintenant la quantifier pour le procédé industriel. Avoir des coefficients de distribution (\(D\)) différents ne suffit pas ; il faut mesurer l'écart relatif entre eux pour savoir si la séparation sera facile ou pénible. Nous allons calculer le Facteur de Séparation (\(\alpha\)), qui est en extraction liquide-liquide ce que la volatilité relative est en distillation : l'indicateur clé de performance thermodynamique du système solvant/soluté.
📚 Référentiel
- Loi de Distribution de Nernst : À l'équilibre et à T constante, le rapport des activités d'un soluté entre deux phases non miscibles est constant.
- Définition de la Sélectivité : Capacité d'un solvant à extraire préférentiellement une espèce par rapport à une autre.
Le coefficient de distribution \(D\) nous dit "où va l'élément".
Ici, \(D_{\text{Nd}}=2.54\) et \(D_{\text{Pr}}=1.82\). Les deux sont supérieurs à 1, ce qui signifie que les deux éléments préfèrent aller dans l'huile (phase organique) que rester dans l'eau. C'est un problème ! Pour les séparer, nous voudrions idéalement que l'un aille dans l'huile et l'autre reste dans l'eau. Comme ce n'est pas le cas, nous devons jouer sur la vitesse relative de transfert ou la préférence relative. Le facteur \(\alpha\) mesure cette préférence. Si \(\alpha\) est proche de 1 (par exemple 1.05), la séparation est quasi-impossible (mélange azéotropique). Si \(\alpha > 2\), c'est facile. Ici, nous nous attendons à une valeur faible.
Le facteur de séparation intrinsèque, noté \(\alpha\) (ou parfois \(\beta\)), est défini comme le rapport des coefficients de distribution de deux solutés. Il représente l'enrichissement de la phase organique en l'élément le plus extractible par rapport à l'autre. C'est une grandeur thermodynamique indépendante des débits (tant qu'on est à l'équilibre et hors saturation).
Par convention, pour obtenir une valeur supérieure à 1 (plus lisible), on place toujours le coefficient le plus élevé au numérateur :
Où \(D_{\text{Nd}}\) est le coeff. de distribution du Néodyme (élément clé) et \(D_{\text{Pr}}\) celui du Praséodyme (impureté).
📋 Données d'Entrée
| Paramètre | Symbole | Valeur |
|---|---|---|
| Coeff. Distribution Néodyme | \(D_{\text{Nd}}\) | 2.54 |
| Coeff. Distribution Praséodyme | \(D_{\text{Pr}}\) | 1.82 |
Un \(\alpha\) de 1.5 signifie que le Nd est 50% plus "soluble" dans l'organique que le Pr. En une seule étape, l'enrichissement est médiocre. C'est ce \(\alpha\) qui dicte exponentiellement le nombre d'étages nécessaires.
📝 Calcul Détaillé
1. Application Numérique
Détail de la manipulation : On divise simplement le coefficient de distribution de l'espèce la plus extraite (Nd, D=2.54) par celui de l'espèce la moins extraite (Pr, D=1.82). Cette division simple nous donne le ratio de sélectivité du solvant.
Interprétation Post-Calcul : Nous obtenons un facteur de séparation de **1.40**. Dans le monde de l'hydrométallurgie standard (Cuivre, Cobalt, Zinc), ce serait considéré comme très mauvais (difficile). Mais pour des Lanthanides voisins, c'est une valeur "standard" et exploitable, à condition d'utiliser un grand nombre d'étages.
La valeur de 1.40 confirme la faisabilité technique mais souligne la complexité économique : il faudra beaucoup de mélangeurs-décanteurs pour amplifier cette petite différence thermodynamique jusqu'à obtenir une pureté de 99.9%.
Ce \(\alpha\) est valable à 25°C. Attention aux variations de température dans l'usine ! Une baisse de T pourrait rapprocher les coefficients D et faire tendre \(\alpha\) vers 1, annulant toute séparation.
🎯 Objectif
Nous entrons maintenant dans le dimensionnement génie des procédés pur. Sachant que la séparation est difficile (\(\alpha = 1.40\)), combien d'étages théoriques (de mélangeurs-décanteurs en série) devons-nous installer pour atteindre l'objectif de récupération de 99.9% du Néodyme ? Nous allons utiliser le modèle analytique de Kremser pour l'extraction à contre-courant, outil puissant pour dimensionner rapidement des cascades d'extraction en régime dilué.
📚 Référentiel
- Génie Chimique : Opérations Unitaires de Transfert de Masse.
- Équation de Kremser-Brown-Souders : Modèle algébrique pour le calcul du nombre d'étages d'équilibre théoriques (NET).
Le nombre d'étages dépend de la lutte entre deux forces : la thermodynamique (le \(D\) qui tire le soluté vers l'organique) et l'hydrodynamique (le débit de solvant qu'on oppose). Le produit des deux donne le Facteur d'Extraction (\(E\)).
- Si \(E > 1\) : La capacité d'extraction du solvant dépasse la charge aqueuse. Le soluté est "aspiré" vers la sortie organique.
- Si \(E < 1\) : Le solvant est saturé ou insuffisant, le soluté "fuit" dans la sortie aqueuse.
Ici, nous voulons extraire le Nd (\(E_{\text{Nd}}\) doit être > 1). Nous allons calculer ce facteur E, puis l'utiliser dans la formule logarithmique de Kremser pour trouver N.
L'équation de Kremser suppose que les coefficients de partage sont constants (isotherme linéaire) et que les phases sont non miscibles. Elle relie la performance de séparation (fraction non extraite) au nombre d'étages \(N\) et au facteur d'extraction \(E\). Elle remplace avantageusement la construction graphique de McCabe-Thiele pour les solutions diluées.
1. Rapport de Débit (R)
2. Facteur d'Extraction (E)
3. Équation de Kremser (Nombre d'Étages N)
Pour une entrée de solvant pur (\(y_{\text{in}}=0\)) :
📋 Données d'Entrée
| Donnée | Valeur |
|---|---|
| Débit Organique \(Q_{\text{org}}\) | 1500 L/h |
| Débit Aqueux \(Q_{\text{aq}}\) | 1000 L/h |
| Coeff \(D_{\text{Nd}}\) | 2.54 |
| Performance Cible (Récupération) | 99.9 % |
Si vous voulez récupérer 99.9% du produit, cela signifie qu'il en reste 0.1% dans la phase aqueuse de sortie. Le rapport de réduction de concentration \(x_{\text{in}}/x_{\text{out}}\) est donc :
C'est ce chiffre "1000" que l'on injecte dans la formule.
📝 Calculs Détaillés
1. Calcul du Rapport de Débit (R)
Détail de la manipulation : On calcule le ratio hydrodynamique en divisant le débit de solvant par le débit d'alimentation aqueuse.
Interprétation : On utilise 1.5 fois plus de solvant que d'eau.
2. Calcul du Facteur d'Extraction du Néodyme (\(E_{\text{Nd}}\))
Détail de la manipulation : On multiplie le coefficient de distribution thermodynamique (\(D\)) par le rapport de débit (\(R\)) calculé précédemment.
Interprétation : \(E_{\text{Nd}} = 3.81 \gg 1\). L'extraction sera très efficace et favorable.
3. Calcul du Nombre d'Étages (\(N\))
Détail de la manipulation :
1. Définition du Ratio de Réduction : On veut 99.9% de récupération. Cela veut dire que la concentration en sortie \(x_{\text{out}}\) est 0.1% de \(x_{\text{in}}\). Donc le rapport \(x_{\text{in}}/x_{\text{out}}\) vaut 1000.
2. Structure de la formule : Le terme principal est \(\frac{x_{\text{in}}}{x_{\text{out}}}(1 - 1/E)\). Il représente la force motrice ajustée.
3. Calcul Logarithmique : On calcule d'abord la valeur numérique à l'intérieur des crochets, on prend son logarithme népérien, puis on divise par \(\ln(E)\).
Interprétation Post-Calcul : Le calcul montre qu'il suffit de seulement 5 étages pour récupérer la quasi-totalité du Néodyme. C'est très peu ! Cela s'explique par la forte affinité (\(D=2.54\)) et le débit généreux de solvant (\(R=1.5\)).
Le résultat de 5 étages semble excellent pour le rendement. MAIS attention ! Qu'en est-il de la pureté ? Le Praséodyme a aussi un \(D > 1\) (\(D_{\text{Pr}}=1.82\)). Avec le même débit de solvant, son facteur d'extraction sera :
Comme \(E_{\text{Pr}} > 1\), le Pr sera LUI AUSSI extrait presque totalement en 5 étages. Conclusion : On récupère tout le Nd, mais il est totalement pollué par le Pr.
Ne jamais dimensionner une séparation sur le seul critère du rendement de l'espèce d'intérêt. Il faut toujours vérifier le comportement de l'impureté (voir Question 4).
🎯 Objectif
L'étape précédente a révélé un piège potentiel : 5 étages suffisent pour le rendement (quantité), mais nous suspectons que la sélectivité (qualité) est désastreuse. L'objectif final est de valider ou d'invalider ce dimensionnement en calculant précisément la contamination du produit final organique par le Praséodyme. Si la pureté cible de 99.9% n'est pas atteinte, nous devrons proposer une modification architecturale du procédé (Lavage/Scrubbing).
📚 Référentiel
- Bilan Matière Global : Conservation de la masse sur l'ensemble de la batterie.
- Formule de Fraction Non-Extraite (\(\phi\)) : Dérivée de Kremser pour calculer la "fuite" d'un soluté dans le raffinat.
Nous avons vu que \(E_{\text{Pr}} = 2.73\). Comme cette valeur est bien supérieure à 1, la mathématique des étages théoriques prédit que le Praséodyme va monter dans la phase organique avec le Néodyme. Nous allons calculer exactement quel pourcentage de Pr se retrouve dans le solvant avec 5 étages. Si ce pourcentage est élevé (proche de 100%), cela prouvera que notre "séparation" n'est en fait qu'un transfert global sans purification.
Pour calculer la fraction d'un soluté qui reste dans l'eau (raffinat) après N étages, on utilise la forme simplifiée de Kremser pour le raffinat :
Cette formule est valide quand \(E \neq 1\). Le taux d'extraction (pourcentage passant dans l'organique) est donc :
📋 Données d'Entrée
| Donnée | Valeur |
|---|---|
| Facteur d'extraction Pr (\(E_{\text{Pr}}\)) | 2.73 |
| Nombre d'Étages (\(N\)) | 5 |
Quand \(E \gg 1\), le terme \(E^{N+1}\) au dénominateur devient énorme, rendant \(\phi\) très petit. C'est mathématique : tout passe dans l'organique.
📝 Calculs Détaillés
1. Facteur d'Extraction du Praséodyme (\(E_{\text{Pr}}\))
Détail de la manipulation : On recalcule le facteur d'extraction spécifique au Praséodyme en utilisant son propre coefficient de distribution \(D_{\text{Pr}}\) et le même ratio de débit \(R\) que précédemment.
2. Fraction de Pr restant dans l'eau (\(\phi_{\text{Pr}}\)) avec N=5
Détail de la manipulation :
1. On injecte \(E=2.73\) et \(N=5\) dans la formule.
2. On calcule la puissance :
3. On effectue la division finale.
Interprétation Post-Calcul : Seulement 0.42% du Pr reste dans l'eau. Cela signifie que 99.58% du Praséodyme est extrait avec le Néodyme !
Avec 5 étages en extraction simple, nous récupérons 99.9% du Nd mais aussi 99.6% du Pr. La pureté du produit final sera quasiment identique à celle de l'alimentation. La séparation a échoué.
Solution Industrielle : Pour séparer des éléments avec un \(\alpha\) aussi faible (1.4), l'extraction simple est impossible. Il faut mettre en œuvre un procédé à REFLUX avec une section de LAVAGE (Scrubbing).
1. On extrait TOUT (Nd+Pr).
2. On "lave" la phase organique chargée avec une solution acide qui redissout préférentiellement le Pr (moins affine).
3. Cela nécessite beaucoup plus d'étages (typiquement 40 à 60 au total).
Ne jamais conclure un dimensionnement sur une simple validation de rendement. La pureté est souvent le critère limitant.
📄 Livrable Final (Note de Calculs EXE)
| Ind. | Date | Objet de la modification | Rédacteur |
|---|---|---|---|
| A | 10/10/24 | Calcul préliminaire (Extraction simple) | Dr. Chimie |
| B | 24/10/24 | Ajout section Lavage (Scrubbing) | Ing. Procédés |
- Contraction des lanthanides confirmée (Rayon Nd < Rayon Pr).
- Facteur de séparation expérimental \(\alpha = 1.40\).
- Isothermes linéaires supposées (domaine dilué).
| Débit Aqueux | 1000 L/h |
| Débit Organique | 1500 L/h |
| Charge Totale (Ln) | 16 g/L |
Dimensionnement pour séparation Nd (Pur) / Pr (Rejet).
Ing. Expert Chimie
Directeur R&D
Laisser un commentaire