Gilles Fevotte

Nesrine Gherras

Jacques Moutte

Batch cooling solution crystallization of ammonium oxalate in the presence of impurities: study of solubility, supersaturation and steady-state inhibition

Batch cooling solution crystallization of ammonium oxalate in the presence of impurities: Study of solubility, supersaturation, and steady-state inhibition

Gilles Févotte, Nesrine Gherras, Jacques Moutte

pollution of the pure solvent by impurities. It is finally proposed to consider the impuritydependent solubility as the appropriate reference for defining supersaturation.

Introduction.

This paper reports partial results of a broader study devoted to the investigation of the effects of impurities during batch cooling crystallization processes. For this purpose, the crystallization of ammonium oxalate in water with nickel sulfate as impurity was selected as a model system. In a previous paper, the nucleation and growth kinetics of the pure solute/solvent system were first characterized [START_REF] Gherras | On the use of process analytical technologies and population balance equations for the estimation of crystallization kinetics. A case study[END_REF] . Nickel sulfate was previously reported to inhibit the growth of ammonium oxalate (AO) crystals. Before analyzing and modeling the batch cooling crystallization of AO in the presence of NiSO 4 , it was essential to characterize the equilibrium features of the impure system, that is to say the solubility data and the way in which the final steady-state of the solid/liquid system is reached through crystallization. In fact, it is well known that the equilibrium solute concentration might differ from the solubility in the presence of impurities. Such possible difference between the steady-state concentrations of pure and impure systems is sometimes referred to as supersaturation threshold in the literature.

Considering that the "pure solubility" denoted as . Surprisingly, such variations of solubility are rarely taken into account in studies dealing with process engineering and/or kinetics of crystallization in real media. This lack of interest for the effects of impurities is all the more surprising that, notably in the industrial context, the presence of impurities is quite inevitable while both increases of solubility and growth inhibition effects are often observed.

The latter raise obvious concerns regarding the quality and the productivity of batch crystallizers.

Inhibition of crystal growth induced by the presence of impurities is widely reported in the literature, but this problem is usually investigated through the experimental study of single crystals performed in specific experimental devices [START_REF] Sangwal | Additives and Crystallization Processes: From Fundamentals to Applications[END_REF][START_REF] Kubota | Effect of Impurities on the Growth Kinetics of Crystals[END_REF] . Growth inhibition is of tremendous importance, particularly in the case of complex multicomponent systems. This is a critical issue, for example, in biological media [START_REF] Heijna | Growth inhibition of protein crystals: A study of lysozyme polymorphs[END_REF][START_REF] Bhamidi | The influence of a homologous protein impurity on lysozyme crystal growth[END_REF][START_REF] Halberstein | Studies in Natural Products Chemistry (Atta-ur-Rahman)[END_REF] , geological systems [START_REF] Kalin | Passive mine water treatment: the correct approach?[END_REF][START_REF] Weijnen | The adsorption of additives at the gypsum crystal surface: A theoretical approach: II. Determination of the surface coverage required for growth inhibition[END_REF] etc., and in the case of industrial crystallization operations [START_REF] Bauer | Purity as an issue in pharmaceutical research and development[END_REF] Indeed, due to the complexity of the many upstream chemical and separation unit operations, impurities are almost unavoidable in industrial processes, at least as secondary products of the many chemical reactions implemented upstream. Growth inhibitions can result in the cessation of the growth of crystals even when supersaturation is not fully consumed. The yield of solid is thus reduced. The effect of impurities is also undesirable because unwanted growth rate reductions require increasing the duration of the crystallization operation. When supersaturation thresholds are observed at the end of solution crystallization operations performed in the presence of impurities, it is therefore not clear whether they should be attributed to growth inhibition effects and/or to solubility drifts. As far as process understanding and modeling is concerned, this is clearly an important issue. Through experimental and theoretical investigations, this work is intended to shed some light on the notion of supersaturation threshold.

In the presence of solubility variations it seems rather logical to compute supersaturation with respect to the actual solubility in impure solvent. Neglecting the presence of impurities in crystallizing media might lead to overestimate (or underestimate) supersaturation and, consequently, to underestimate (or overestimate) the effective basic crystallization rates.

Despite its obvious importance in crystallization kinetics and process dynamics studies, this particular question is almost never addressed in papers dealing with the evaluation of crystallization kinetic effects induced by impurities. Considering the problem of the definition of appropriate reference solubility in impure inhibiting media, two different standpoints can actually be adopted to interpret the observation of steady-state supersaturation thresholds. The first approach assumes that supersaturation thresholds actually result from erroneous interpretation of solubility increases. The second approach consists in confronting experimental and theoretical aspects of the notion of supersaturation thresholds on the ground of crystal growth inhibition models (e.g. Kubota Cooling ammonium oxalate (AO) crystallization operations were performed in pure water and in water containing dissolved NiSO 4 impurities. The choice of dissolved nickel sulfate was essentially motivated by the availability of the many published data concerning the effect of Ni 2+ ions on the growth rate of AO in water [START_REF] Sangwal | Effect of impurities on metastable zone width for the growth of ammonium oxalate monohydrate crystals from aqueous solutions[END_REF][START_REF] Sangwal | Effect of cationic impurities on solubility and crystal growth processes of ammonium oxalate monohydrate: Role of formation of metaloxalate complexes[END_REF] . Supersaturation measurements were performed using in situ ATR FTIR spectroscopy. The operating mode and the calibration procedure required by the spectroscopic technique were presented elsewhere 1 , and will not be recalled here. The lab-scale instrumented crystallizer used for the experiments is presented in Fig. 1. More details about the experimental setup are given in the preceding reference.

AO solubility measurements were done in pure water and in water containing dissolved nickel sulfate, as presented below. Batch cooling unseeded crystallization experiments were then carried out with constant cooling rates. The kinetics of crystallization of AO in pure water was identified and the modeling results were published [START_REF] Gherras | Etude expérimentale et modélisation par bilans de populations des cinétiques de nucléation et de croissance d'opérations discontinues de cristallisation par refroidissement en absence et en présence d'impuretés[END_REF] . Two different crystallization regimes were highlighted and the corresponding nucleation and growth kinetics were characterized, according to the rate of cooling (i.e. the rate of generating supersaturation).

To investigate the variations of solubility in impure water, the following set of impurity concentrations was explored: C imp =[6 10-3; 2 10-2; 4 10-2] mol.L -1 .

Solubility in pure water, experimental results.

Two different operating procedures were used to measure the solubility of AO in pure water and in the presence of dissolved nickel sulfate; both experimental strategies taking advantage of the availability of "continuous" in situ ATR FTIR concentration measurements through the application of successive isothermal/poly-thermal dissolution and crystallization steps [START_REF] Févotte | New perspectives for the on-line monitoring of pharmaceutical crystallization processes using in situ infrared spectroscopy[END_REF] . As shown in Fig. 2, a first method consists in cooling AO from a temperature of about 62°C. After reaching the metastable zone, the onset of primary nucleation, at a temperature of about 50°C, can be observed in Fig. 2. The decreasing solute concentration was then quantitatively monitored using ATR spectroscopy. For rapidly growing crystals, provided that the cooling rate is sufficiently low, it can be assumed that the concentration trajectory matches the solubility (i.e., supersaturation is negligible). The measured concentration trajectory thus provides a continuous estimate of the solubility. The same approach can also be reproduced through heating, assuming that the level of undersaturation is also negligible during the dissolution process. In the case of pure water, the fact of applying both heating and cooling procedures ensures that the solubility equilibrium is truly reached. Indeed, should the crystallization and/or dissolution dynamics be underestimated the cooling concentration trajectory will be located above the heating trajectory. The coincidence of both curves therefore confirms that the estimation of the solubility curve is correct. This overall continuous solubility measurement method is referred below to as Dynamic Heating/Cooling Solubility Approach (DHCSA)

Figure 2 shows the results obtained after applying the DHCSA: the 2 curves measured through ATR FTIR spectroscopy are almost identical to data published by [START_REF] Menczel | The molar enthalpies of solution and solubilities of ammonium, sodium and potassium oxalates in water[END_REF] . The relative uncertainty between the measured solubility data is only of the order of 0.12%. It should however be noted that data previously reported [START_REF] Seidell | Solubilities of Inorganic and Metal Organic Compounds[END_REF] are unexpectedly about 5.5% lower than our experimental results. The following polynomial expression was finally fit to the spectroscopic measurements: Figure 2: Continuous measurements (DHCSA) of the solubility of AO in pure water, using a slow cooling crystallization procedure (Method 1) followed by a slow heating dissolution process (Method 2), and comparison with two literature data sets. The average relative uncertainty between cooling (method 1) and heating (method 2) and data after Menczel et al. [START_REF] Menczel | The molar enthalpies of solution and solubilities of ammonium, sodium and potassium oxalates in water[END_REF] is about 0.12%.

) ( * 0 T C = 2.
Figure 2 shows that both heating and cooling plots differ only imperceptibly, and therefore allow computing an accurate average solubility curve. A second approach for the determination of solubility consists in implementing successive positive or negative temperature increments. At every steady-state temperature stage, the time variations of Mid-IR spectra are monitored until a constant spectrum can be recorded. Using calibration models, the recorded spectra allows computing the equilibrium solute concentration. This second approach is referred below to as Static Heating/Cooling Solubility Approach (SHCSA).

As in the case of DHSCA, it is clear that both steady state crystallization and dissolution measurements are necessary to confirm the achievement of the "true thermodynamic equilibrium".

Solubility with dissolved NiSO 4 , experimental results

Both SHCSA and DHSCA were applied to the evaluation of the solubility of ammonium oxalate in the presence of varying concentrations of NiSO 4 . A saturated AO solution was prepared at 60°C and slowly cooled with a rate of 2°C/h until a final temperature of 10°C.

Successive temperature steps were then performed from 10°C to 60°C. During each step the temperature was kept constant during a 30 h period. 240 ATR spectra were recorded after 30 h stabilization with a sampling period of 30s (i.e. for a period of 2 hours.) As displayed in Brzoska 11 , is also presented in Fig. 3 in insert, the reported results are in good agreement with our experimental results.

In the presence of different impurity concentrations, the parameters of the following polynomial expression were estimated, where

) ( * T C imp
is the solubility in kg AO monohydrate/kg de solution:

) , ( * imp imp C T C = aT 2 + b T +c (2)
The estimated parameters a,b,c are displayed in Table 1.

Figure 3 shows that below a temperature located around 50-55°C, the difference between the AO solubilities, with and without NiSO 4 impurities, increases for increasing impurity concentrations and decreasing temperatures. It also appears that, within the limits of the arbitrary fit of the second order polynomial eq.( 2), the solubility of AO with 6 10 -3 mol.L -1 dissolved NiSO 4 matches the solubility of the pure system. Above 55°C, as presented in more details in the sequel, the solubility drift observed with the two "high" impurity concentrations is difficult to estimate and, as a first approximation, can be considered negligible (i.e. within the limits of the concentration measurements)

In a study dedicated to the effect of cationic impurities on the solubility of AO, the following empirical relationship was proposed [START_REF] Sangwal | Effect of cationic impurities on solubility and crystal growth processes of ammonium oxalate monohydrate: Role of formation of metaloxalate complexes[END_REF] where A S is a constant depending the cationic impurity specie:

imp S imp C A C C + = * 0 * (3)
Unfortunately, Eq. ( 3) was evaluated only at 30°C. For Ni 2+ , the parameter value A S =1.676 ±0.06 was proposed by the authors with * 0 C = 0.496 mol/L (the solubility of AO in pure water.) Considering the results presented above, Eq.( 3) predicts an increase in solubility proportional to the concentration of impurity, which is not observed in Fig. 3.

Denoting * 0 * * C C C imp imp - = ∆
, the difference between the solubility in impure solvent and the solubility in pure water * 0 C , Fig. 4 presents experimental results expressed as a function of temperature and impurity concentration. These results were computed using Eq. ( 2) and the estimated parameters are given in Table 1. As published by Sangwal and Mielniczek-Brzoska 11 , the solubility increases proportionally to

the impurity concentrations at low temperatures. Fig. 4 shows that the slope of the variations agrees quite well with the value of A s (see Eq.( 3)) proposed by Sangwal and Mielniczek-Brzoska 11 . However, above 50°C, the A s seems to decrease and tends to zero for increasing temperatures. Moreover, due to the measurements uncertainties, the general trend of the solubility drift at low impurity concentration remains unclear. Anyway, the results suggest that above 50°C the solubility in the presence of impurities might be lower than the solubility in pure water.

The variations of solubility with the concentration of impurities can have several explanations. First, the composition of the solution could vary due to the formation of new chemical species resulting from solute/solvent/impurity interactions (e.g. formation of new complexes, chemical reactions, formation of new solvates, etc.) Second, the impurity could incorporate into the crystal lattice which, consequently, could behave as a new solid phase.

Also, in the presence of impurities, both solid and liquid phase variations might also take place in the neighborhood of the growing crystal surface. This is the reason why, to evaluate the possible incorporation of impurities in the crystal, analyses of the solid were performed.

Final suspension samples were withdrawn from the reactor, filtered and washed quickly with cold water in order to remove residual mother liquor. For an impurity concentration of 0.02 mol.L -1 NiSO 4 , X-ray diffraction did not allow concluding to any alteration of the AO crystal structure. XPS showed traces of nickel and almost undetectable traces of sulfur at the crystal surface, while Raman spectroscopy revealed very slight differences that could be related to small changes of the crystal lattice (i.e., appearance of tiny bands and slight shifts). The characterization of the solids therefore suggests that small amounts of Ni and sulfate species are present at the crystal surface, without being incorporated into the crystal lattice.

Speciation of the AO/NiSO 4 /H 2 O system

According to the previous observations, one can attribute the variations of solubility to the effect of the dissolution of NiSO 4 on the speciation of the AO solution and, perhaps, to the alteration of the solid surface by adsorbed impurities. Actually, as far as the chemistry of ionic impure solutions is concerned, ammonium oxalate dissociates in several species which are likely to change with temperature and through interactions with dissolved impurities. A theoretical study of the speciation of the solution was thus performed, firstly to evaluate the impact of dissolved nickel sulfate on the chemical composition of the crystallizing suspensions at the solubility equilibrium, secondly to check the effect of temperature on the speciation of the solution.

The dissociation of AO in pure water involves the following chemical equilibria 11 : system are those available in the latest update (slop07.dat) of the SUPCRT database for calculation of temperature dependent properties of aqueous species according to the Helgeson-Kirkham-Flowers (HKF) equation of state [START_REF] Tanger | Calculation of the Thermodynamic and Transport-Properties of Aqueous Species at High-Pressures and Temperatures -Revised Equations of State for the Standard Partial Molal Properties of Ions and Electrolytes[END_REF] . Concerning oxalate-bearing species, this database integrates the datasets from refs [START_REF] Shock | Organic-Acids in Hydrothermal Solutions -Standard Molal Thermodynamic Properties of Carboxylic-Acids and Estimates of Dissociation-Constants[END_REF][START_REF] Prapaipong | Metal-organic complexes in geochemical processes: Temperature dependence of the standard thermodynamic properties of aqueous complexes between metal cations and dicarboxylate ligands[END_REF] . Table 2 gives the equilibrium constants computed with ARXIM at different temperatures according to the HKF equations. The values are given for the reactions of formation from the following primary species:

[(NH 4 ) 2 C 2 O 4 , H 2 O] solid → [(NH 4 ) C 2 O 4 ] - Soluble complex + NH 4 + -log k 1 (25°C) =0.92 [(NH 4 ) C 2 O 4 ] - → NH 4+ + C 2 O 4 -2 -log k 2 (25°C) = 1.43 C 2 O 4 -2 + H 2 O → HC 2 O 4 -+ OH - -log k 3 (25°C) = 10.2 HC 2 O 4 -+ H 2 O → H 2 C 2 O 4 + OH - -log k 4 (25°C) = 8.8 NH 4 + + H 2 O → NH 4 OH + H 3 O + -log k 5 (25°C)= 4.75
H 2 O, H + , C 2 O 4 -2 , NH 4 + , Ni +2 , SO 4 -2 .
Besides aqueous species, the equilibrium constants for the dissolution of ammonium oxalate salt, not found in the database, were calculated as a function of temperature, using the experimental solubility data.

Considering the high ionic strength of AO-saturated solutions, activity coefficients should be computed using specific interaction models (e.g. Pitzer or SIT approach), but, as such interaction coefficients are currently not available for all species in consideration, activity coefficients have been computed according to a more general truncated Davies model [START_REF] Samson | Modeling chemical activity effects in strong ionic solutions[END_REF] , which gives at least a semi-quantitative understanding of the role of oxalate complexes. The equation proposed for the activity coefficient of the charged solute species i is as follows:

(4)

where A and B are temperature dependent solvent parameters, I is ionic strength of solution, and z i is the charge of the solute species considered. This adaptation is considered to expand the validity of the Davies equation to ionic strength up to 1.5 mol.L -1 .

Figure 5 shows the computed repartition of dissolved species, in the case of the pure AOwater system, as a function of pH. As batch cooling experiments are performed between 50 and 10°C, pH evolves between 6 and 7.

Estimation of the dissolution enthalpy from solubility data

The following Eq. ( 5) describes the temperature-dependence of the solubility in ideal solutions (i.e., the interactions between solute and solvent molecules is assumed to be equal to interactions the solute and the solvent molecules themselves 20 ):

        - ∆ = T T RT H T C f f f 1 ) ( ln * (5)
where T f is the solute fusion temperature and ∆H f is the enthalpy of fusion.

For non-ideal solutions, it was proposed to approximate Eq. ( 5) as a function of the solution enthalpy sol H ∆ . For exemple, [START_REF] Menczel | The molar enthalpies of solution and solubilities of ammonium, sodium and potassium oxalates in water[END_REF] expressed the solubility of salt hydrates as follows:

      - ∆ = T T RT H T C m m Sol 1 ) ( ln * υ (6)
The enthalpy of solution is expected to depend on the ratio between the solid and the solvent. This is why the heat of solution at infinite dilution, ∞ ∆ sol H , is defined as the heat exchanged when one mole of solute is dissolved in an infinite volume of solvent. Parameter υ is expected [START_REF] Menczel | The molar enthalpies of solution and solubilities of ammonium, sodium and potassium oxalates in water[END_REF] to account for the incomplete dissolution of ammonium oxalate into oxalate and ammonium ions, and also for possible slight variations of the solubility product. For concentrated solutions, υ is usually set between 1 and 2. υ=1.5 was assumed in the sequel.

The plot of the experimental results ln(C*) vs. 1/T is displayed in Fig. 7. The plot fits satisfactorily the model Eq. ( 6) and leads to the estimates of ∆H sol presented in Table 3. The average value for ∆H sol in pure water is thus found to be very close to values reported by 14 : 33.0 ± 4.9 kJ/mol. Table 3 shows that the solution enthalpy decreases for increasing impurity concentrations.

As explained in the sequel, DSC analyses which were tried to evaluate the fusion enthalpy turned out to be unreliable. To the best of our knowledge, no thermodynamic data is available that would allow predictive model-based computation of the solubility from the knowledge of the concentration of impurity. Anyway the main trends of our experimental observations are rather well explained by the first approach of solution thermodynamics developed above. In particular, the fact that increasing impurity concentrations increase the solubility -the latter effect being amplified for decreasing temperature -is quite consistent with adsorption theories developed in the context of crystal growth 2,21 .

4.

Notions of inhibiting adsorption of impurities and residual supersaturation.

Experimental observation of supersaturation thresholds and inhibition effects

After determining the corresponding solubility curves, supersaturation profiles were computed from the ATR FTIR concentration measurements performed during cooling crystallization operations in impure solvent. In order to better investigate the kinetics, the batch operations were also carried out with varying cooling rates [START_REF] Gherras | On the use of process analytical technologies and population balance equations for the estimation of crystallization kinetics. A case study[END_REF] . Figure 8 shows two typical supersaturation vs. temperature profiles obtained with a cooling rate of -20°C/h, for two different concentrations of impurities. Supersaturation was computed with respect to the actual solubility (i.e., with dissolved impurity):

* * ) ( ) ( imp imp C C t C t - = σ (7)
The relative uncertainty of the "measured" supersaturation was estimated to be of the order of 3%. When compared to the crystallization advancement in pure solvent, the overall solute consumption is clearly slowed down so that the supersaturation remains high during the cooling process following the onset of primary nucleation. It should also be noted from Fig. 8 that the primary nucleation burst is delayed. The arrows indicate the maximal values of supersaturation corresponding to the so-called "limit of metastable zone". The supersaturation profiles exhibit the following maximal values measured during the primary nucleation peak: σ lim =0.4 for C imp =0.02 mol/L and σ lim =0.22 for C imp =0.006 mol/L. The first nucleation burst is therefore delayed for increasing concentrations of impurities, which confirms the overall inhibition of the crystallization process. The average value of σ lim in pure water was about 0.2 for batch experiments performed with the cooling rate dT/dt=-20°C/h 1 . The enlargement of the metastable zone obtained with C imp =0.006 mol/L is thus negligible, despite the importance of the effect of the dissolved impurities on the overall concentration trajectory.

At the end of the batch process, a significant level of supersaturation remains after a 3-hours isothermal maturation stage. In order to check the reality of the crystallization stopping, the suspension was kept at 20°C for 30 hours. Figure 9 shows that the concentration goes on decreasing slowly before reaching its equilibrium value. Moreover, the equilibrium finally corresponds to the solubility measured with dissolved NiSO 4 . Three important points can be outlined from these observations. At the end of the batch cooling process, the crystallization seems to be stopped (i.e., the solute concentration does not seem to evolve any more, at the time scale of the batch process.) Rather than plotting σ=f(T) in Fig. (8), plotting σ=f(t)

shows that the growth rate is not nil at the end of the batch process (i.e. during the final isothermal stage): the variations of the solute concentration suggest strong inhibition phenomena. The apparent steady-state supersaturation threshold characterizing the system (in the sense of inhibition models) requires a significant maturation time to be reached. It addition to these two observations, it should be outlined that the final supersaturation threshold coincides with the increase of solubility investigated above (i.e. the increase of solubility resulting from the presence of impurities in water/AO solution.) One could therefore conclude that the final supersaturation threshold is nil, provided we consider the supersaturation reference as defined with respect to the actual solubility (i.e., the solubility in impure solvent.) From these results, the undesirable phenomenon referred to as a supersaturation threshold should be called into question. This is why the notion of residual supersaturation and its theoretical basis are now examined. The physics of the adsorption of impurity species was studied using molecular simulation modeling which was expected to allow better understanding of the interactions between the AO crystal and adsorbates. In particular, as presented below, molecular simulation allowed estimating the adsorption energy of Ni and SO 4 ions. The results were then compared to the experimental data obtained above.

It should firstly be outlined that several papers dealt with the inhibiting effect of the adsorption of Ni 2+ during the growth of AO crystals while the possible adsorption of SO 4 2was not clearly evaluated in the past. Actually, dynamic simulation studies performed during this work suggested that sulfate ions, to a lesser extent than Ni 2+ , could also play a role in the inhibition of the crystallization of AO. To evaluate this assumption, cooling solution crystallization experiments were performed in the presence of 0.0019 mol.L -1 nickel sulfate and ammonium sulfate, with a cooling rate of -30°C/h. Some pictures taken during these two experiments are presented in comparison with AO suspensions obtained in pure water (Fig. 10a vs. 10b). The deterioration of the crystal surfaces and shapes, and the presence of crystal fragments in Fig. 11b suggest that the presence of sulfate ions also affects the overall crystal generation process. Growth inhibition caused by the presence of sulfate ions was also confirmed thanks to ATR FTIR AO concentration measurements during batch cooling experiments which are not presented here, and by further related experimental observations (see Fig. 9). The likely adsorption of both Ni 2+ and SO 4 2-ions on AO crystal faces was also confirmed by the molecular modeling simulations presented in the sequel.

The separate effect of nickel and sulfate ions was also evaluated through monitoring batch cooling crystallization experiments performed with 0.0019 mol/L dissolved NiCl 2 and NiSO 4 .

The corresponding solute concentration profiles are presented in Fig. 11

Molecular modeling simulation of the adsorption of ionic species on AO crystal faces.

To compare the interaction between AO crystal faces and both Ni 2+ and SO 4 2-ions, the competitive adsorption occurring on the different AO crystal faces through was evaluated thanks to molecular modeling simulations. First of all, the theory of attachment energies predicts an equilibrium AO crystal structure which is rather close to the structure described in the literature [START_REF] Küppers | Anisotropy of thermal expansion of ammonium and potassium Oxalates[END_REF][START_REF] Küppers | The crystal structure of ammonium hydrogen oxalate hemihydrate[END_REF] , but the energy values were "adapted" so as to improve the fit between the equilibrium and experimental crystal shape. The simulation was carried out with 4000 atoms, with a cut-off radius of 13 Å. With respect to their tendency to attract ions, the {010} and {110} faces present a particular interest. The {010} face is rather steep with sharp edges of exposed ammonium species. The face contains electropositive wells that can interact with SO 4 2-ions and electronegative edges that can interact with Ni 2+ . The {110} face is smoother than the {010}, with edges favorable to the adsorption of SO 4 2-. For the sake of conciseness, the representation of the crystal structure will not be given here.

Additional molecular simulations were again performed to bring information about the possible adsorption of ionic species. As far as the {101} face is concerned, the adsorbed Ni is linked through 4 dative covalent bond, with three neighbor oxalate ions, and one oxygen of the water molecule. However, in view of the deeper position of adsorbed Ni, the ammonium ions are more facing to the outside. Sulfate ions can also interact with two ammonium ions through electrostatic bonds, and through weak hydrogen bonds with two oxygen atoms of neighboring water molecules. Owing to its small size and to its interactions with the environment, Ni 2+ could incorporate the crystal lattice, at least at a shallow depth, while sulfate would rather adsorb at the surface. A representation of the structure of the {101} face and of its interactions with nickel and sulfate species is presented in Fig. 14. Due to the lack of place, the representation of the {110} face and of its interactions is not presented in this paper. Regarding the {110} face, Ni 2+ ions interact through dative covalent bonds with 3 oxygens of the oxalate groups and one oxygen of the water molecule located above. In a fashion similar to the {010} face, SO 4 2-can adsorb on the {110} face. Both nickel and sulfate ions are therefore likely to restrain the growth of AO crystal through adsorption, and some incorporation of Ni 2+ through the {010} face cannot be excluded.

Molecular simulation was also applied to the systems AO/NiSO 4 for estimating the adsorption energies of the two impurity ionic species on the different crystal faces. The main computation results are presented in Table 4. When the difference between the attachment energy and the adsorption energy of a given specie on a given surface is positive, the adsorbed specie has low impact of the crystal growth process. In contrast, when the difference is negative, the impurity slows down or even blocks the face growth, leading to a spreading of the latter face. Weaker differences of energies induce increased competition between the solute molecules and impurities, which can cause significant changes in the crystal shape. As explained above, the {010} face is mostly affected by the presence of impurities and, to a lesser extent, so is the {110} face. The Stranski-Krastanov model, which mainly applies to organic crystals where Van der Waals bonds predominate, was used for the computations. However, as electrostatic energies are mainly involved in the growth of AO, even though adsorption energies contribute significantly to the growth, it was important to evaluate the relevance of the computations. Such validation was carried out through the comparison between the calculated total energy of the crystal and its experimental estimate.

Using GenMol®, the calculated crystal energy is -443 kJ.mol -1 . The experimental evaluation of the crystal energy was also tried using DSC (Differential Scanning Calorimetry).

Unfortunately, in the case of AO, the technique is not fully reliable. Indeed, AO starts to decompose at 70°C, until 200°C, while water starts evaporating, with a maximal evaporation rate observed around 115°C. Between 200°C-235°C and 235°C-300°C, melting followed by vaporization and concomitant decomposition are observed. Due to the simultaneity of the involved thermal phenomena, estimating the energy of the AO crystals is therefore rather hazardous. The sublimation energy was still assessed assuming that the sublimation energy is given by the sum of the energy of vaporization of structural water and of the AO melting/vaporization energy. A "pseudo sublimation energy" was thus evaluated to be 229 kJ.mol -1 . Considering that the crystal energy is of the order of twice the sublimation energy, it can therefore reasonably be concluded that the adsorption energies evaluated using GenMol® are acceptable.

Growth inhibition models and supersaturation thresholds

Several models describing the inhibition of crystal growth by impurities were proposed in the literature [START_REF] Sizemore | A New Model for the Effect of Molecular Imposters on the Shape of Faceted Molecular Crystals[END_REF][START_REF] Martins | The Influence of Impurities on the Crystal Growth Kinetics According to a Competitive Adsorption Model[END_REF][START_REF] Martins | Unsteady-state inhibition of crystal growth caused by solution impurities[END_REF][START_REF] Kubota | Effect of Impurities on the Growth Kinetics of Crystals[END_REF] . According to Cabrera and Vermilyea 27 , the elementary growth steps are hindered by the adsorption of impurity species which results in the socalled "pinning effect". When the impurity spacing becomes narrower than the critical radius of step curvature, the model predicts the cancellation of the growth rate. As a consequence, when supersaturation σ defined by Eq.( 7) falls below a critical threshold σ c , the pinning effect leads to a "dead zone" where the growth steps are stopped. Cabrera and Vermylyea model was further reexamined by Kubota and Mullin who considered that the coverage of the crystal surface, expressed thanks to Langmuir's adsorption theory, should be corrected by the efficiency parameter α introduced in the following Eq.(8). To allow further crystal growth, the growth-step has to circumvent the adsorbed impurity, which obviously reduces the overall growth rate. Several models describing such "pinning mechanism" were early described and many convincing experimental observations of the pinning mechanism were reported in the literature, notably thanks to the use of advanced imaging techniques such as AFM (Atomic Force Microscopy, see e.g. [START_REF] Thomas | Molecular properties of adsorbates that affect the growth kinetics of archerite (KDP)[END_REF][START_REF] Ristic | Does Impurity-Induced Step-Bunching Invalidate Key Assumptions of the Cabrera-Vermilyea Model?[END_REF][START_REF] Weaver | Inhibition of calcium oxalate monohydrate growth by citrate and the effect of the background electrolyte[END_REF] ).

According to Cabrera-Vermylyea and KM theory, the ratio Γ between the step velocities in pure solvent (u 0 ) and in the presence of dissolved impurity (u) is expressed as a function of the coverage fraction of inhibited active sites θ and of an overall effectiveness factor α :

n T L kT a u u θ σ ω αθ θ σ γ         - = - =         - = = Γ 1 1 1 0 ( 8
)
where γ is the edge free energy, a is the size of the growth unit, T is the absolute temperature, k is the Boltzmann constant, θ is the coverage fraction of active growing crystal surface by adsorbed impurities and L is the average distance between the adsorbed impurity species. Exponent n depends on the hindered active site: n=1 when the adsorption of impurity occurs at kink sites and n=1/2 when random surface active sites on terrace are concerned. In the sequel, n was set to 1.

Assuming that the overall growth rate G of a crystal is proportional to the step velocity u, Eq.( 8) yields:

( ) αθ - = 1 0 ) t ( G ) t ( G (9)
According to Eq.( 8), α does not depend on properties of the crystal only, but also on temperature and supersaturation 2 . When the adsorption process reaches its steady-state, the equilibrium coverage is referred to as θ * . According to KM Model [START_REF] Kubota | A kinetic model for crystal growth from aqueous solution in the presence of impurity[END_REF][START_REF] Kubota | Supersaturation dependence of crystal growth in solutions in the presence of impurity[END_REF] the equilibrium coverage of the growing surface is simply estimated thanks to Langmuir adsorption theory where K is the Langmuir adsorption constant:

        + = imp mp i KC KC 1 * θ (10)
In the present case, both Ni and SO 4 species are likely to adsorb on active kink sites, and to hinder the growth of AO. The overall coverage of growing surface by the two ionic impurity species accounts for the two different adsorption constants and for the concentration of the two impurity species [START_REF] Sangwal | Additives and Crystallization Processes: From Fundamentals to Applications[END_REF] . Assuming negligible incorporation of nickel sulfate during the crystallization and equality between the two ionic concentrations we get the following equation expressing the competitive adsorption of the two ionic species:

imp imp SO SO SO imp Ni Ni Ni SO SO SO SO SO Ni Ni Ni Ni Ni C C K K C K K C K C K C K C K                 + - +         + - =         + - +         + - = 4 4 4 4 4 4 4 4 1 ) 1 ( 1 ) 1 ( 1 ) 1 ( 1 ) 1 ( * θ θ θ θ θ (11)
Average values of the adsorption energies of nickel and sulfate displayed given in Table 4, make it possible to calculate the adsorption constants:

      = kT exp SO4 or Ni , SO4 or Ni ads E K (12)
The coverage fractions of Ni and SO 4 impurities was computed using Eqs.( 11) and ( 12), for an impurity concentration of 0.02 mol.L -1 . The results are displayed in Figure (13). The adsorption of sulfate cannot be neglected, even though Ni 2+ is the main adsorbed impurity. 8) and ( 9), the growth rate is cancelled (i.e., G=0) when the supersaturation falls below σ c such that:

θ ω θ γ σ T kTL a c = = (13)
As already outlined, many studies dealing with the effect of impurities on crystallization kinetics consider the solubility of the pure system as the actual solubility, which is questionable. In particular, such an assumption implicitly means that the equilibrium concentration measured in the presence of impurities should actually be considered as a threshold of supersaturation. Assuming that the difference

* 0 * * C C C imp imp - = ∆
, computed from the solubility data displayed in Fig. 3, can be regarded as the critical supersaturation σ c , the critical supersaturation should amount to the following ratio:

* 0 * / C C imp c ∆ = σ (14)
Combining Eqs.( 11) to (13) and assuming that the supersaturation threshold is given by Eq.( 13) allows estimating ω as a function of temperature:

imp imp SO SO SO imp Ni Ni Ni imp c C C K K C K K C C T T kL a                 + - +         + - ∆ = = = 4 4 4 1 ) 1 ( 1 ) 1 ( θ θ θ σ γ ω (15)
A plot of ω=f(T) is displayed in Fig. 13 for C imp =0.02 mol.L -1 . The variations exhibit a very significant decrease for increasing temperature. Knowing that the liquid-solid interfacial tension is expected to increase for increasing temperature [START_REF] Mersmann | Crystallization Technology Handbook[END_REF][START_REF] Spaepen | The temperature dependence of the crystal-melt interfacial tension: a simple model[END_REF] , the results suggest that the average distance between the active sites L increases very strongly with temperature, which is quite unlikely. Furthermore, an estimate of ω was proposed in a recent paper [START_REF] Fevotte | On Multiple Nucleation Bursts During Solution Crystallization in Pure and Impure Solvent[END_REF] reporting partial results about the dynamic population equation balance modeling of the crystallization of AO in the presence of 0.006 mol/L dissolved NiSO 4 . In this paper, ω was assumed constant and, in the average, found to be equal to 1080 K. This latter value was obtained through parameter optimization based on data acquired during batch cooling experiments (i.e. using supersaturation profiles like the ones in Figs 8 and9). However, ω was obtained without accounting for the adsorption of sulfate and no estimate of the adsorption energy of Ni 2+ was used. Other unpublished estimates accounting for the adsorption of both impurity species and using the estimated values of the adsorption energy in Table 4 were further obtained. It was finally found that, in the average, ω ≅ 350 K. This estimate [START_REF] Gherras | Etude expérimentale et modélisation par bilans de populations des cinétiques de nucléation et de croissance d'opérations discontinues de cristallisation par refroidissement en absence et en présence d'impuretés[END_REF] is probably more reliable than the previous value ω=1080 K since it was computed from several experimental datasets and led to satisfactory simulation of the whole crystallization process..

Discussion.

Most studies dealing with the dynamic modeling of inhibiting effects of impurities on crystal growth remain rather vague about the appropriate definition of the solubility * imp C in the presence of impurities, and about its use for computing supersaturation in impure solvents. As already outlined, this vagueness is all the more surprising that possible drifts of the final solute concentration during crystallization processes have an obvious effect on the yield of industrial batch operations and, consequently, on their economical performances. Indeed, apart from the problem of possible changes of "crystal quality" through the incorporation of foreign species, the presence of impurities affects the solution speciation and is thus expected to affect the solubility * 0 C of the pure system. In this context two different stances can be adopted for analyzing the growth inhibition problem even though both viewpoints can hardly be considered as compatible. On the one hand, one can assume that apparent changes of solubility in impure solvent would result from supersaturation thresholds described by inhibition models (e.g. Kubota-Mullin's pining mechanism model,

Martins competitive adsorption model, Doherty model of molecular imposters …)

Consistently with this first point of view, the definition of solubility and supersaturation should be unchanged with respect to the case of the pure system. On the other hand, it seems logical to consider that the actual solubility in impure solvent should be the appropriate reference for computing supersaturation. For increasing concentration of impurities, supersaturation would therefore be lowered when the solubility increases in the presence of dissolved impurities.

It is not straightforward to be definitive about the correct viewpoint to adopt in the case of the AO/NISO 4 /water system investigated here. From a thermodynamic viewpoint, while no incorporation of impurities in the crystal lattice was observed, it is clear that the speciation of the liquid phase is changed by dissolved impurities, as well as the properties of the liquidsolid interface. It is therefore consistent to observe drifts of the solubility in the presence of dissolved impurities. In such conditions, there is no logical reason to go on considering the solubility in pure solvent as a correct reference. In other words, it is difficult to support the idea that the difference ∆C imp = * imp C -* 0 C is a "true" supersaturation threshold, as predicted by inhibition theories. Moreover, as outlined in the experimental part of this paper, different solubility curves of impure systems were determined through both cooling and heating procedures. The fact that the "heating" and "cooling" curves coincide suggests that the solute-solvent system returns to a "true" equilibrium in both cases. It should also be noted, even though it was not fully conclusive, that the speciation study presented above allowed explaining the main trends of the observed drifts of solubility, in particular as a function of temperature. This argues again in favor of considering that each solute-solvent-impurity system is characterized by its own solubility curve and, consequently, that the appropriate definition of supersaturation in impure systems should refer to the impurity-concentration dependent solubility (i.e., not to the unique curve obtained in pure solvent).

Moreover, the assumption that ∆C imp is a supersaturation barrier was analyzed in the light of Kubota-Mullin pinning model. This analysis was based on estimates of the adsorption energy of nickel and sulfate impurity species. Identifying ∆C imp to σ c (i.e., the critical supersaturation defined by the growth inhibition model) allowed estimating ω=γa/kL, a lumped parameter characterizing the crystal/solution system. Compared to estimates of ω obtained during a more thorough unpublished study, the estimates of ω seem to be significantly underestimated, which suggests that ∆C imp should not be considered as a supersaturation threshold.

Conclusions and perspectives.

The present study of the ternary AO/NiSO 4 /water solute/impurity/solvent system aims at investigating the inhibition effects of impurities on the development of batch cooling crystallization processes. New experimental and theoretical results demonstrate unambiguously the inhibition of the growth of crystals by dissolved impurities. In such conditions, the transient behavior of the crystallization process appears to be rather complex and, before analyzing and modeling the nucleation and growth kinetics involved, it was first necessary to characterize the equilibrium features of the system. As a wide as possible range of techniques was applied to analyze the significance of residual supersaturation thresholds observed at the end of batch cooling crystallization experiments performed in the presence of NiSO 4 dissolved impurities: dynamic continuous and steadystate multi-step equilibrium concentration measurements, calorimetric analysis, solid state analysis (Raman, XRS, etc.), and molecular simulation modeling.

The quantitative information provided thanks to the various experimental and theoretical approaches applied during this work was put in contrast with Kubota-Mullin impurity adsorption theory. The main conclusion drawn from the latter comparison is that increases of final equilibrium concentrations with increasing impurity concentrations should rather be attributed to solubility drifts than to supersaturation thresholds (in the sense of Kubota-

Mullin pinning model)

This latter conclusion is somewhat unexpected because, in the present case-study, it seems to call into question the notion of "supersaturation thresholds" developed in the framework of growth-inhibiting adsorption mechanisms. We thus suggest that most papers dealing with the effect of impurities during solution crystallization processes present blind spotssometimes erroneous interpretations-about the right way of defining supersaturation. In the presence of impurities, using the right reference solubility and defining correctly supersaturation requires accounting for solubility variations. This is why we finally propose that no reliable study dealing with the effect of impurities on the kinetics of crystallization should be carried out without a rigorous preliminary evaluation of the variations of solubility in impure solvents.

For table of content use only.

Title: Batch cooling solution crystallization of ammonium oxalate in the presence of impurities: study of solubility, supersaturation and steady-state inhibition.

By FEVOTTE Gilles, GHERRAS Nesrine & MOUTTE Jacques Synopsis.

Papers dealing with the inhibiting effects of impurities in solution crystallization are rather unclear about the appropriate way of defining supersaturation in impure solvents. Setting the correct reference for computing supersaturation in impure solutions is obviously a key-issue for understanding and modeling the dynamics of elementary crystallization phenomena (i.e., nucleation, crystal growth, etc.). This is the object of the present paper.

  refers to the solubility of AO in pure water, low concentrations of dissolved impurities usually increase the actual solubility such that >

2 .

 2 -Mullin pinning model.) Both viewpoints are examined in the following. The paper is organized as follows. After presenting the experimental means implemented during the work, Part 3 is focused on solubility measurements obtained in pure water (Part 3.1) and in the presence of varying concentrations of dissolved nickel sulfate (Part 3.2.) For interpreting these experimental results, a speciation study of the impure system is presented in Part 3.3. In order to get an additional point of view on the solubility study, Part 3 is completed by a brief analysis of the variations of the dissolution enthalpy with varying concentrations of impurities. Part 4 is then devoted to the evaluation of possible steady-state supersaturation thresholds in impure water. Experimental results are presented and analyzed in Part 4.1 and compared with molecular modeling estimates of the adsorption of the involved impurity species on the two main AO crystal faces in Part 4.2. Thanks to the adsorption energies computed using molecular simulation, Part 5 finally presents a basic quantitative approach of supersaturation thresholds. Experimental crystallization setup and continuous in situ measurements of supersaturation during the crystallization of ammonium oxalate monohydrate.

Figure 1 : 3 .

 13 Figure 1: Experimental set up: temperature controlled well-mixed batch crystallizer equipped with in situ ATR-FTIR diamond probe and in situ image acquisition device.

  27 10 -5 (± 6.7 10 -8 ) T 2 + 5.39 10 -4 (± 5.6 10 -6 ) T + 2.87 10 -2 (± 4.8 10 -3 ) concentration in pure water in kg AO/kg solution and T is the temperature in °C.

Fig 3 ,

 3 Fig 3, three different dissolved impurity concentrations were selected: C imp = [0.04; 0.02; 0.006] mol NiSO 4 .L -1 . Solubility data obtained at 30°C, published by Sangwal and Mielniczek-

Figure 3 .

 3 Figure 3. Measured solubility curves of AO (in kg AO/ kg solvent) in pure solvent and with varying concentrations of NiSO 4 . Insert: data from 11 . "smoothed" Cooling DHSCA data (Continuous line) and heating SHSCA data (full points)

Figure 4 .

 4 Figure 4. AO solubility variations ∆C* imp computed using the experimental results displayed in Fig.3 with temperature varying between 10 and 60°C.

Figure 5 .

 5 Figure 5. Variations with pH of the concentrations of the various solute species for a 1 molal AO solution at T=25°C.

  the concentration of species in solution were computed as a function of the concentration of NiSO 4 and of temperature. The results allowed evaluating the differences of solubility displayed in Fig.6. The amplitude of the experimental solubility variations is a little bit more than twice as high as the model predictions, but the main trends of the temperature dependence are consistently reproduced by the speciation model. Fig.6 alsoshows that according to the same model, the trend of the variations of solubility with the concentration of NiSO 4 reverses around 40°C. The same feature was also observed with the experimental results displayed in Fig.4 around 50°C.Actually, the program ARXIM, used for speciation calculations, implements the Helgeson-Kirkham-Flowers (HKF) equation of state for aqueous species. The computations thus rely on the availability of HKF parameters that allows the computation of a set of dissociation constants for all species of interest. The latest update of the HKF-oriented database (slop07.dat) includes all the aqueous species that should be taken in consideration for speciation calculations in the H 2 O -AO -NiSO 4 system, but it lacks data for the ammonium oxalate salt. We have thus derived the equilibrium constant of ammonium oxalate hydrate, at various temperatures, from experimental data on its solubility in the H 2 O -AO system, using the speciation program ARXIM with the Davies -Samson activity model, and the data on all aqueous species of interest (H + , OH -C 2 O 4 ) -)). It was then possible, using the same speciation program with the same activity model, and adding the Ni-related species (Ni +2 , NiOH + , NiO, HNiO -2 , NiO 2 -2 , Ni(C 2 O 4 )), to compute the solubility of the ammonium oxalate salt for systems containing various amounts of NiSO 4 and for different temperatures.

Figure 7 .

 7 Figure 7. Plot of the logarithm of the solubility of ammonium oxalate (molar fraction) vs. 1/T in pure water and in the presence of varying nickel sulfate concentrations.

Figure 8

 8 Figure 8 Typical supersaturation profiles measured during the cooling crystallization of AO performed with a cooling rate of -20°C/h, with different concentrations C imp of nickel sulfate.

Figure 9 .

 9 Figure 9. Solute concentration profiles measured during the cooling crystallization of AO performed with 2 cooling rates (dT/dt= -20 and -10 °C/h) and constant concentration of nickel sulfate C imp = 0.02 mol.L -1 .

  also recalling the time variations of concentration measured in pure water after applying identical operating conditions. The experiments demonstrate a stronger inhibiting effect of Cl -than SO 4 2-but do not allow comparing the inhibiting effect of the adsorption of Ni 2+ and SO 4 2-. It can therefore be assumed that the overall growth inhibition process results from a competitive inhibiting adsorption of the dissolved ionic species during the crystallization of AO.

Figure 10 .

 10 Figure 10. Pictures taken during the cooling crystallization of AO performed with a cooling rate of -30°C/h, in the presence of 0.0019 mol/L (a) Nickel sulfate; (b) Ammonium sulfate, and in (c) Pure water.

Figure 11 .

 11 Figure 11. Comparison between the concentration profiles measured during the AO cooling crystallization performed at a cooling rate of -30°C/h in the presence of 0.0019 mol.L -1 nickel chloride and nickel sulfate.

Figure 12 :

 12 Figure 12: (a) Representation of the {010} AO face in the presence of Ni 2+ . For the sake of visibility the atom of nickel was expanded 3 times. (b) Stick representation of the same face in the presence of SO 4 2-. (c) Representation of the {010} face with adsorbed sulfate anion.

Figure 13 .

 13 Figure 13. Estimated temperature variations of parameter ω and of the coverage fractions of Ni and SO 4 species for an impurity concentration C NiSO4 =0.02 mol.L -1 .

Table 1

 1

	C i (mol/L)	a	b	c	R 2
	0	2.27 10 -5	5.39 10 -4	2.87 10 -2	0.999
	0.0006	2.06 10 -5	6.06 10 -4	2.92 10 -2	0.999
	0.02	1.72 10 -5	7.68 10 -4	3.17 10 -2	0.999
	0.04	1.64 10 -5	6.73 10 -4	4.0 10 -2	0.999

. Parameter estimates for the polynomial approximation of the solubility curves of ammonium oxalate with different concentrations of dissolved nickel sulfate (see eq.2)

Table 2 .

 2 2, one can see that 3 species are predicted to coexist during the whole crystallization process at almost constant concentration: ammonium and oxalate ions and the complex anion [(NH 4 )C 2 O 4 ] -; the latter being involved in the crystallization process. Equilibrium constants of the aqueous species computed according to the HKF equation of state. The values are for the reactions of formation from the following set of primary species: H 2 O, H + , C 2 O 4

			-2 , NH 4	+ , Ni +2 , SO 4 -2 . The constants for ammonium oxalate salt
	were calculated from the solubility data, using the same activity model as used in
	subsequent equilibrium calculations
	T (°C)	10	20		25	30	40	50
	OH -	-14.52	-14.16	-13.99	-13.83	-13.53	-13.26
	HSO 4 -	1.81	1.92		1.98	2.04	2.17	2.3
	NH 3	-9.72	-9.4		-9.24	-9.09	-8.81	-8.54
	NiOH +	-11.31	-10.96	-10.8	-10.63	-10.31	-10.01
	NiO	-22.25	-21.2	-20.7	-20.22	-19.3	-18.43
	HNiO -2	-32.17	-31.37	-30.99	-30.62	-29.9	-29.23
	NiO 2 -2	-45.88	-44.62	-44.02	-43.45	-42.36	-41.35
	H 2 C 2 O 4	5.47	5.5		5.53	5.56	5.63	5.72
	HC 2 O 4 -	4.21	4.24		4.26	4.28	4.33	4.39
	Ni(C 2 O 4 )	5.16	5.16		5.16	5.17	5.21	5.25
	NH 4 (C 2 O 4 ) -	0.9	0.91		0.92	0.94	0.98	1.03
	(NH 4 ) 2 (C 2 O 4 ),s	2.57	2.24		2.09	1.95	1.71	1.47

As for the pure system, a study of speciation accounting for the effect of dissolved nickel sulfate was again performed using ARXIM®. The equilibrium constants which are not available in the literature were calculated as a function of temperature, using the experimental solubility data. The computation was performed with saturated oxalate solutions at temperatures varying between 25°C and 50°C and NiSO4 concentrations between 10-12 and 0.095 mol.L -1 . The speciation model shows that increasing the concentration of dissolved NiSO 4 decreases the concentrations of oxalate, H-oxalate and NH4-oxalate ions while the concentration of NH 4+ ions increases. For example, with a concentration of NiSO 4 of 0.095 mol.L -1 , pH=6.16 at 25°C and 5.84 at 50°C. Nickel in solution is present as aquo-complexes involved in acid-base reactions that release H 3 O + ions.

Consequently, the pH decreases and thus affects the speciation of oxalates and ammonium ions. Moreover, oxalate ions are trapped in soluble nickel complexes, their concentration is approximately equivalent to the concentration of dissolved nickel sulfate. The concentration of ammonium ions increases in the same proportions. 4 , H 2 O) with increasing nickel sulfate concentrations and with temperature varying between 10 and 50°C.

Table 3 :

 3 Estimated values of ΔH sol (KJ/mol), the enthalpy of dissolution of AO estimated from linear plots of lnC* against 1/T of experimental solubility data obtained with varying NiSO 4 concentrations.

	C imp (mol/m 3 )	ΔH sol (kJ/mol)	R 2
	0	32.52	0.99
	6 10 -6	31.31	1
	2 10 -5	30.06	0.99
	4 10 -5	26.43	0.99

Table 4 :

 4 Computed adsorption energies (in kJ/mol) of sulfate and nickel species on the crystal faces subject to interactions with both impurities..

	Crystal	Attachement	SO 4	-2 Adsorption	Ni 2+ Adsorption	(E att -E ads )	(E att -E ads )
	Face	Energy		Energy	Energy	SO4 -2	Ni 2+
	110	-20.5		-26.85	-23.2	-6.4	-2.7
	010	-0.5		-20.39	-23.7	-20.9	-24.2

Acknowledgements:

We greatly acknowledge the French Research Agency ANR (Agence Nationale de la Recherche) for the support granted to the "white project" IPAPI (Improving the Properties of Active Pharmaceutical Ingredients), ref.07-BLAN-0183.

We also thank the following researchers for their very kind and valuable assistance: Professor K. SANGWAL for his invaluable help and wise advices about the manuscript and Doctor Gerard PEPE (Research Director at UPR-CNRS 3118) for his expertise and his contribution to the molecular simulation study.


Gloria Heilbronn

Cécile Robin

Laurie Barrier

Sylvie Bourquin

Marc Jolivet

Yingying Jia

Laure Guerit

Bihong Fu

Evolution paléogéographique et paléotopographique du Tian Shan chimois au mesozoïque : une approche sédimentologique et thermochronologique

ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

ÉVOLUTION PALÉOGÉOGRAPHIQUE ET PALÉOTOPOGRAPHIQUE DU TIAN SHAN CHINOIS AU MÉSOZOIQUE : UNE APPROCHE SÉDIMENTOLOGIQUE ET THERMOCHRONOLOGIQUE

Gloria HEILBRONN (1,@) , Cécile ROBIN (1) , Laurie BARRIER (2) , Sylvie BOURQUIN (1) , Marc JOLIVET (1) , Yingying JIA (3) , Laure GUÉRIT (2) , Bihong FU (3) (1) Géosciences Rennes, UMR CNRS 6118, Université de Rennes 1, Rennes, France.

(2) Institut de Physique du Globe de Paris, Sorbonne Paris Cité, UMR CNRS 7154, Université Paris Diderot, Paris, France.

(3) Laboratory of Lithospheric Evolution, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing, China.

(@) gloria.heilbronn@gmail.com Le Tian Shan est une chaîne intracontinentale d'Asie Centrale dont la structure lithosphérique résulte de l'accrétion de divers blocs au cours du Paléozoïque. Son histoire tectonique récente est marquée par sa réactivation au Tertiaire, liée à la collision entre l'Inde et l'Asie. Plusieurs études ont montré l'influence des structures tectoniques héritées sur la localisation de la déformation et donc de la topographie cénozoïque. Toutefois, l'histoire topographique et tectonique de la chaine entre les deux orogenèses majeures, à savoir paléozoïque terminal et cénozoïque, reste peu contrainte. La synthèse des données bibliographiques montre des incohérences majeures entre l'histoire topographique déduite des données de thermochronologie basse température sur la chaine, et celle déduite des séries sédimentaires exposées dans les piémonts.

Afin de mieux contraindre cette évolution, nous avons étudié six coupes réparties entre les piémonts nord et sud, ainsi que dans les bassins internes à la chaine. L'objectif de cette étude est de déterminer les grandes variations spatiales et temporelles des profils de dépôt depuis le Trias supérieur jusqu'au Paléocène. Ces variations sont alors interprétées en termes de proximalitédistalité par rapport à la chaine et donc en termes de croissanceérosion des reliefs. Un modèle d'environnement de dépôts intégrant sédimentologie de faciès, analyse pétrographique et mesures de paléocourants a donc été établi. La corrélation des principales séquences de dépôts entre les différentes coupes a permis de déterminer le degré de connexion entre les bassins au cours du temps, et ainsi d'avoir une idée qualitative de la topographie mésozoïque. En combinant nos données avec celles de la bibliographie (thermochronologie et géochronologie détritiques et sur le socle, sédimentologie et géomorphologie), un modèle détaillé de l'évolution topographique mésozoïque a été construit réconciliant les interprétations sédimentologiques et géochronologiques.

Le démantèlement des reliefs associés à la chaine tardi-paléozoïque s'achève au Trias supérieur. Le Jurassique est caractérisé par une activité tectonique très faible, majoritairement transtensive et le développement d'une surface majeure de pénéplanation sous un climat relativement humide. Toutefois, la mise en place de cônes alluviaux éphémères associés à des dépôts éoliens au Jurassique supérieur -Crétacé inférieur indique une réactivation de la topographie et des conditions climatiques plus arides. Le Crétacé est ensuite marqué par une activité tectonique limitée, qui s'arrête totalement au Crétacé terminal permettant le développement au Paléocène d'un niveau majeur de calcrêtes.