Evaluation of ePC-SAFT for pH Calculation in Aqueous Itaconic Acid Solutions at High Ionic Strengths

In the context of the raw material change for sustainable production of chemicals, the selected bio-based amino acids and carboxylic acids are considered as promising platform chemicals. After fermentation, the acids are present in aqueous solutions with many side components and elevated ionic strength. The ionic strength is even further increased when pH-shift operations are applied for the separation of the target compounds. Since high ionic strengths strongly affect the solution properties, particularly the solid–liquid-equilibrium and the dissociation equilibrium in the solution, the high ionic strengths and the resulting effects on the solutions must also be taken into account in process modeling and the design of downstream processes. Various models have been reported in the literature but the majority cannot be applied for predicting the solution composition and pH at high ionic strengths. In this work, a procedure for the calculation of the composition, i.e. the distribution of the present species and pH, of aqueous itaconic acid solutions based on ePC-SAFT is developed and evaluated at different levels of ionic strengths. The ePC-SAFT parameters of itaconic acid are determined based on experimental solubility data from literature. The resulting model is validated with experimentally measured titration curves and compared with the ideal model and the Davies approximation to the Debye–Hückel model. It is demonstrated that the ePC-SAFT approach provides reliable results at high ionic strengths with lower deviations from experimental data than current literature models.


Introduction
Biotechnological processes represent an important basis for the sustainable production of chemicals from renewable raw materials. In these processes, mostly aqueous media with complex composition are produced, which have to be handled in downstream process steps. pH values in electrolyte solutions containing weak diprotic acids at high ionic strengths above 1.0 mol ⋅ L −1 .
In this work, the dynamical approach by Glaser et al. [1] is combined with the ePC-SAFT equation of state [14] to obtain a model for the calculation of species concentrations at equilibrium, including pH, in aqueous itaconic acid solutions, with validity at higher ionic strengths. The diprotic carboxylic acid is chosen as a model substance, since it represents a potential bio-based platform chemical. Thus, the model developed in this work can be used in future studies on the development of suitable downstream processes for this target molecule. The relevant ePC-SAFT parameters for itaconic acid are determined applying a procedure based on the approaches by Ruether and Sadowski [15] and Lange et al. [13]. The model is compared with the ideal model and with the Davies approximation to the Debye-Hückel theory at different ionic strengths. Additionally, the shifts in the dissociation equilibria of the itaconic acid species at increasing ionic strength are investigated.

Chemicals
The chemicals used in the experiments are specified in Table 1.

Titrations
Titrations of aqueous 2-methylidenebutanedioic acid (trivial name: itaconic acid (IA)) solutions with different sodium chloride (NaCl) concentrations were conducted with 20 wt% sodium hydroxide (NaOH) solution at 298.15 K. For two separate titrations, the following initial compositions were selected: 70 g ⋅ L −1 IA and 0 g ⋅ L −1 NaCl, and 55 g ⋅ L −1 IA and 55 g ⋅ L −1 NaCl. The concentrations are given in grams substance per liters solvent (water). The chosen IA concentrations are in the range of IA solubility in water at low pH and temperatures [16] and are representative for crystallization processes. The difference in the IA concentration in the initial solutions was chosen due to the reduced solubility of IA in water at increased sodium chloride concentration [16]. The solutions were stirred over night on a heated magnetic stirrer plate IKA C-MAG HS10 digital to ensure complete dissolution of solids. The titrations were carried out by adding the sodium hydroxide solution with a 50 mL burette to the continuously stirred aqueous IA solution. The temperature control was realized using a Pt100 sensor and a water bath heated by the heated magnetic

Modeling of Aqueous Electrolyte Solutions and pH Calculation
In this work, a model for aqueous electrolyte solutions containing IA, the pH modifying agents hydrochloric acid (HCl) and sodium hydroxide (NaOH) and the corresponding ions is developed. For the mathematical description of this system, a dynamical approach, describing all dissociation reactions, is chosen. First, the equilibrium state is defined for the dissociation of a monoprotic acid HA, shown in Eq. 1.
The thermodynamic equilibrium constant K a describes the dissociation equilibrium. In the definition in Eq. 2, the nonidealities in the solution are taken into account by considering activities instead of molar concentrations. The nonidealities cover all interactions between molecules present in the solution.
In this work, K a is divided into two terms, shown in Eq. 2. Q represents the quotient of the concentrations of the products and reactants, K represents the corresponding quotient of the activity coefficients. The conversion of activities and activity coefficients of different units is described in the appendix (Eqs. 46-48).
In dilute solutions, activity coefficients converge to unity and the equilibrium state can be described by the acid constant K acid [3]. This component specific property is measured by titration methods and extrapolated to infinite dilution [17,18]. Thus, it is not valid for high concentrations.
In Eqs. 4 to 9, the reaction system of interest for aqueous electrolyte solutions containing IA is defined along with the corresponding equilibrium constants.
The composition of the solution can be calculated on the basis of the rates of change of each component, applying the general mass action kinetics [19][20][21]. The rates of change of each component of the investigated reaction system are listed in Eqs. 10 to 20. The nonidealities are considered by replacing the molar concentrations with activities in the reaction kinetic [22]. Velocity constants are commonly determined in experimental investigations, but they are not documented for a wide range of substances. As suggested by Schell [3], k f is arbitrarily set to 100 and k b is calculated by the equilibrium constant.
In this work, ePC-SAFT is used to calculate the activity coefficients. This approach is compared to the ideal calculation ( i = 1 ) and the Davies approximation [6] to the Debye-Hückel theory [23]. In the Davies Eq. (eq 22), the activity coefficients z of ions mainly depend on the ionic strength I.
The coefficient A = e 2 B∕(2.3038π k B T) is calculated with the electron charge e , the static dielectric constant of water ϵ , the Boltzmann constant k B , the temperature T, B = (2e 2 N A ∕ k B T) 0.5 and the Avogadro constant N A [24,25]. At 298.15 K , A has an approximate value of 0.5108 kg 0. 5 ⋅ mol −0.5 and B is approximately 0.3287 ⋅ 10 8 kg 0. 5 ⋅ cm −1 ⋅ mol −0.5 [25][26][27]. The Davies approximation includes the empirical parameter b with a constant value for all ions. Davies original work assigns b = 0.2 , and this value was shown to give improved activity coefficients for large anions at low ionic strength based on conductivity measurements [6]. Schell et al. [3] showed, that b = 0.1 leads to more accurate results in the dihydrogen phosphate buffer system, so this value for b was chosen in this work. The ionic strength is calculated with the ion charge z i and the molar concentration c i , as shown in Eq. 23.
The derivation of the Debye-Hückel Theory and thus also the Davies approximation are based on mole fractions as concentration unit [6,23]. Hence, the ionic strength is calculated with the mole fraction of the ions. In diluted solutions, the mole fractions can be replaced by molar or molal concentration as a simplification [1,3,5]. In this work, we use molar concentrations to calculate the ionic strength (Eq. 23) and the activity coefficients (Eq. 22), even if we investigate concentrated solutions. The calculations with the Davies approximation are used as references to the ideal calculations and the ePC-SAFT model, to demonstrate the limits of the Davies approximation and the potential of ePC-SAFT at high concentrations.
In the cases of ideal calculation and the Davies approximation, the presented system of ordinary differential equations is solved. The equilibrium constants for the dissociation of IA are obtained from the acid constants K acid at 298.15 K from literature [28]. Hydrochloric acid, sodium hydroxide and sodium chloride fully dissociate in water, but are considered in the defined reaction system for the purpose of a general model. The values of K NaOH and K NaCl are set to 10 7 , which is the reported value of the acid constant of hydrochloric acid [28], to fulfill the condition of full dissociation. K water is calculated with the ion product of water K w with the reported value of 10 14 mol 2 ⋅ L −2 [28]. In the approach with ePC-SAFT, the ordinary differential equation system is simplified. The dissociation equilibria of hydrochloric acid, sodium hydroxide and sodium chloride are neglected due to missing ePC-SAFT parameters. The simplified ordinary differential equation system is given by the Eqs. 24 to 29.
The equilibrium constants of the dissociation of IA are fitted simultaneously with the ePC-SAFT parameters, as suggested by Lange [13]. The fitting procedure is presented in detail in Sect. 3.3.
The presented systems of ordinary differential equations (ODE) are solved with the MAT-LAB-solver ode23s. The results are the molar concentrations of each component. Thus the pH value in the solution can be determined as defined in Eq. 30 [29].
In this work, the dynamical implementation is only used for the determination of the species concentrations at steady state as proposed by Glaser et al. [1]. For this purpose, the evaluation time for the ODE system was set to 0.1 s , which ensured that the steady state concentration of each component and thus of the complete specified system was reached.

ePC-SAFT Equation of State
In this work, the dynamical approach for pH calculation by Glaser et al. [1] is combined with the ePC-SAFT equation of state for the calculation of activity coefficients. ePC-SAFT is extensively described in literature [11,12]. Based on the residual Helmholtz energy a res , different thermodynamic quantities are calculated with ePC-SAFT. a res takes into account the repulsion of the reference system (hc), attractive forces (disp), hydrogen-bonding interactions (assoc), as well as Coulomb interactions (ion) and is calculated with Eq. 31.
To describe non-associating molecules, three pure component parameters, the segment diameter , the segment number m and the dispersion-energy parameter ∕k B are necessary [11]. Associating molecules require two additional pure component parameters, the association-energy parameter A i B j ∕k B and the association-volume parameter A i B j [12]. The dispersion and association energy parameters are commonly divided by the Boltzmann constant k B . The pure component parameters for non-associating components in mixtures are combined with the mixing rule of Berthelot-Lorentz [30], shown in Eqs. 32 and 33.
k ij is the binary interaction parameter, which can be considered to adjust the dispersionenergy in mixtures. In this work, a linear dependency of k ij on the temperature is assumed Eq. 34).
The mixing rules from Wolbach and Sandler [31] allow the description of associating components in mixtures, shown in the Eqs. 35 and 36.
The ion-ion interaction is a central aspect for the description of electrolyte solutions, which is taken into account by the Helmholtz-energy contribution a ion . For this, Cameretti [14] applied a Debye-Hückel [23] term (Eq. 37) that requires no additional parameter, only the charge of the ions z i .
i is defined by Eq. 38 and is the Debye screening length.
Besides the contribution of ion-ion interactions to the Helmholtz-energy a ion , Held et al. [10] introduce an effect on the dispersion energy ij between two ions. The dispersion interaction of ions with the same charge is neglected ( ij = 0). The activity coefficients are calculated with the fugacity coefficients obtained from ePC-SAFT. Two definitions are relevant for the description of electrolyte solutions. They differ in their reference state, the pure substance as shown in Eq. 39 or the infinite dilution, shown in Eq. 40 [32].
The pure component reference state is used to calculate the activity coefficients of the solvent and the target component. The activity coefficients of the solutes that only exist in solution, in this study HIA − , IA 2− , OH − , H 3 O + , require the infinite dilution reference state.

Estimation of ePC-SAFT Parameters
For the application of the ePC-SAFT equation of state for the calculation of activity coefficients in aqueous solutions containing IA, the ePC-SAFT parameters for IA need to be determined. In this section, the determination procedure is described.
The pure component parameters are usually fitted to the liquid density and the vapor pressure of the pure component. However, pure IA exists only in the solid state under normal conditions. For this reason, the parameters are fitted to solubility data in different solvents, as suggested by Ruether and Sadowski [15]. The solid-liquid equilibrium is described with the equation according to Prausnitz et al. (Eq. 41) [33].
x L i is the intrinsic solubility, the mole fraction of the fully protonated IA species in solution at the limit of solubility. x,L i is the corresponding activity coefficient in the liquid phase, Δh SL 0i the enthalpy of fusion at the melting temperature T SL 0i and R the universal gas constant. In the first step of the determination procedure, the pure component parameters of IA are determined by fitting to the solubility of IA in seven organic solvents: methanol, ethanol, ethylacetate, 1-propanol, 2-propanol, acetone and acetonitrile in the temperature range from 283.15 K to 328.15 K. The experimental data are taken from Yang et al. [34]. Due to the use of organic solvents, dissociation of IA is avoided and therefore not considered. Besides the pure component parameters, the temperature-dependent binary interaction parameters of IA and the solvents are fitted. The binary interaction parameters are implemented with a temperature dependency, as shown in Eq. 34. The fitting procedure of this step is shown in Fig. 1.
First, the input data for the parameter determination procedure is listed. Temperature, pressure and the temperature-dependent solubility are given by the experimental data from Yang et al. [34]. The ePC-SAFT pure component parameters of the solvents were taken from literature and are listed in the appendix (Table 6).
At the beginning of each iteration step the fitting variables are set by the MATLAB solver. These are the pure component ePC-SAFT parameters of IA and the binary interaction parameters of IA for each solvent. With the solubility x IA = x H 2 IA [34], the activity coefficients ePC−SAFT (ePC-SAFT) and exp (Eq. 41) are calculated. The objective function, which represents the difference between the calculated activity coefficients, is formulated as a least-squares problem and solved with the MATLAB solver lsqnonlin [35]. The variable v in Fig. 2 represents the vector that contains the fitted ePC-SAFT parameters.
In the second step, the binary interaction parameters of IA and water are determined by fitting the pH-dependent solubility of IA in water at 298.15 K using the experimental data from Holtz and Goertz [16]. In this step, dissociation of IA is considered. The fitting procedure of this step is shown in Fig. 2.
Again, the input data for the parameter determination procedure is listed first. Temperature, pressure and the pH-dependent solubility of IA in water are given by the experimental data from Holtz and Goertz [16]. The ePC-SAFT pure component parameters of undissociated IA are already known. The same parameters are used for the dissociated IA species, only the charge and association sites are adjusted, as suggested by Lange et al. [13]. At the beginning of each iteration step, the fitting variables are set. These are the binary interaction parameters k ij of the IA species and water, as well as the thermodynamic equilibrium constants of both dissociation steps K a,1 and K a,2 . The thermodynamic equilibrium constants are integrated into the fitting procedure, as the literature values [28] were determined assuming infinite dilution. Due to the high ionic strength at the solubility limit, the literature values do not describe the dissociation equilibrium correctly. Binary interaction parameters k ij of zero and values from literature for the acid constants are chosen as reasonable starting points. The inner iteration loop is used to determine the solution composition, including the species distribution, and activity coefficients at each experimental data point of the pH-dependent solubility [16]. The loop starts with an estimated species distribution so that the experimental pH value is reached by the dissociation of the fully protonated IA species and initial values of i = 1 (ideal case). The solution composition is calculated with the ODE reaction system defined in Eq. 24-29. The resulting composition is then used to calculate the activity coefficients with ePC-SAFT using the current parameter values for the next iteration step. If the change in composition falls below the termination criterion , the calculation of the composition and the corresponding activity coefficients is finished. For the termination criterion, a value = 10 −6 was chosen. The fitting is performed with the functions f 1 and f 2 in Eqs. 42 and 43, which represent the difference between the calculated pH value pH calc and the experimentally measured pH value pH exp and the difference between the calculated activity coefficient ePC−SAFT for the fully protonated species and the activity coefficient obtained for the fully protonated species at solid-liquid equilibrium exp = L i (Eq. 41). The functions f 1 and f 2 are evaluated for each data point. The objective function is formulated as a least-squares problem, as shown in Eq. 44 and solved with the MATLAB solver lsqnonlin. v represents again the vector that contains the fitted ePC-SAFT parameters. (42)

Parameter Estimation and Fitting Results
The ePC-SAFT equation of state enables the calculation of activity coefficients of the components in electrolyte solutions. The ePC-SAFT parameters for the IA species obtained from the parameter estimation procedure (Sect. 3.3) are shown in Table 2. In addition, the parameters of other components considered in this work are listed in Table 2, including their literature reference.
As mentioned in Sect. 3.3, the parameters of the different IA species only differ in the number of association sites and charge. The starting values for IA used in the determination procedure were set to the parameter values determined for succinic acid in the work by Lange at al. [13] and are listed in Appendix ( Table 5). The value of the association volume is (A i B i ) ∕k B = 0.020069 and therefore lies between 0.01 and 0.03, which is considered typical for organic substances [15]. The binary interaction parameters k ij of IA and the organic solvents are temperature-dependent, as presented in Eq. 34, and are listed in the appendix (Table 7).
In Fig. 3, the modeling results for the solubility of IA in organic solvents based on the determined ePC-SAFT parameter set are compared to the experimental data. The predicted solubility is in good agreement with the experimental data, which demonstrates the validity of the parameter set for the protonated acid in the temperature range from 283.15 K to 328.15 K.
The binary interaction parameters of IA and water, as well as the equilibrium constants of both dissociation steps were determined based on the pH-dependent solubility data of IA [16]. The resulting parameters and dissociation constants (Eqs. 4 and 5) are listed in Tables 3 and 4, respectively.
The modeling results for the solubility of IA in water at different pH and the experimental data are compared in Fig. 4. The measured solubility in the pH range from 2.2 to 2.8 has a constant value of x IA = 0.0128 . Above a pH of 2.8, the solubility increases nonlinearly with increasing pH up to x IA = 0.0301 at a pH of 3.77. The calculated solubility of IA is in good   [10] agreement with the experimental data, in particular in the pH range from 2.2 to 3.5. The slope of the calculated curve increases continuously with increasing pH, in accordance with the dissociation of IA. At pH values above 3.5, the deviation between calculated solubility and measured solubility increases. In the same ranges, the errors of the experimental data increase as well. The averaged relative deviation between the calculated and experimental solubility data is 4.65 % and therefore in the range of the investigation carried out by Lange et al. [13] for the pH-dependent solubility of succinic acid.  x IA x IA,ePC-SAFT  The presence of NaCl leads to a decrease of the measured pH value, in agreement with Hartono et al. [37]. Comparing the data of the titration curves, the measured pH is decreased by approximately 0.2 pH units at 55 g ⋅ L −1 NaCl. Since in ideal solutions the additional NaCl should not affect the measured pH, the decrease of the pH can be attributed to the interactions between the solutes and demonstrates the nonideal behavior of the system. Thus, the nonidealities influence the dissociation equilibria of IA, which was also reported for other components, like methyldiethanolamin (MDEA) [37] and 2,2,3,4,5,5-hexamethylimidazolidin-1-oxyl (HMI) [38].
The experimental data is used for the comparison and validation of different models for the calculation of activity coefficients at different levels of the ionic strength. The ideal calculation with activity coefficients equal to 1 is compared with the Davies approximation (Eq. 22) and the ePC-SAFT model (Eq. 39 and 40). A comparison and an evaluation of these models on the basis of resulting pH was chosen, since it enables the comparison with the experimental data. In addition, the obtained activity coefficients by itself are regarded not informative. In Fig. 6, the calculation results for the different models are shown for the two solutions with different NaCl concentrations. Figure 6a shows the titration in the absence of NaCl. The ideal calculation ( = 1 ) leads to good results in the pH range from 2.5 to 4.0. In this range, the deviations are smaller than 0.16 pH units. At higher pH values, the ionic strength increases due to dissociation and the difference between calculated and measured pH value increases in the same way. The pH value is overestimated with the ideal calculation. The difference can be explained with the intensified interactions between the increased number of ions, which is neglected in the ideal calculation. The calculation with ePC-SAFT leads to qualitative similar results. In contrast to the ideal calculation, the pH value is underestimated over the entire pH range. Compared to the ideal calculation and the calculation with ePC-SAFT, the Davies approximation for the activity coefficients leads to better results in the pH range from 2.5 to 3.5. The ionic strength exceeds 0.1 mol ⋅ L −1 at the third titration step, representing the reported limit for the validity of the model [9], and reaches a value of 0.25 mol ⋅ L −1 at pH 3.5. Above pH 3.5, the difference between the measured and calculated pH values determined with the Davies approximation exceeds the difference determined with the ePC-SAFT approach. At a salt concentration of 55 g ⋅ L −1 NaCl, shown in Fig. 6b, the deviation between the ideal calculation and the experimental data increases. The deviation is at least   Fig. 6 Comparison of the pH determined using different models for the calculation of activity coefficients at different NaCl concentrations: (a) 70 g ⋅ L −1 IA, 0 g ⋅ L −1 NaCl; (b) 55 g ⋅ L −1 IA, 55 g ⋅ L −1 NaCl. Experimental data is reproduced from Fig. 5. The ionic strength is calculated neglecting the activity coefficients 0.25 pH units. In this case, the consideration of the activity coefficients with ePC-SAFT leads to excellent results in the pH range between 2.5 and 4.5. Deviations in this pH range are lower than 0.1 pH units. At pH above 4.5, the pH value is again underestimated, but the deviation to the experimental data is in the same range, as for the ideal calculations. The Davies approximation and ePC-SAFT lead to similar results in the pH range from 2.5 to 3.7. Between pH 3.7 and 5.0 the deviations of the measured and calculated pH with the Davies approximation are between 0.1 and 0.2 pH units and are therefore greater than with the ePC-SAFT calculation. Figure 6a and b demonstrate that the activity model ePC-SAFT is able to calculate the pH value of aqueous solutions with IA up to ionic strengths of 1.75 mol ⋅ L −1 accurately. Besides the more accurate pH calculation at high ionic strength, the advantage of the ePC-SAFT calculation over the ideal calculation and the Davies approximation is the knowledge of the activity coefficients of non charged components that can be used in further models, for example in crystallization or adsorption models.

Influence of the Nonidealities on the Dissociation Equilibrium
To investigate the influence of the nonidealities on the dissociation equilibrium at the different NaCl concentrations, K a,i , Q i and K γ,i , as defined in Eq. 2, are plotted in Fig. 7 as function of the ionic strength. For each titration step of each titration, the dissociation constants K a,i (Eqs. 4 and 5) are calculated as products of Q i and K γ,i , which are obtained from the with ePC-SAFT calculated species concentrations and activity coefficients. The concentration quotients Q 1,Davies and Q 2,Davies are calculated with the Eqs. 2 and 22 as comparison to the ePC-SAFT approach. Using the Davies equation, the activity coefficients and therefore K γ,i only depend on the ionic strength and Q i,Davies can be calculated with the knowledge of K a,i and the ionic strength as independent variable. The equilibrium constants K a,i only depend on the temperature and remain constant for each dissociation step and NaCl concentration. The negative logarithm of the concentration quotient pQ i for each dissociation step of IA is below the corresponding pK a,i value. Schell et al. [3] and Zars et al. [5] showed for different buffer systems, using the Davies approximation, a nonlinear decrease of pQ i with increasing ionic strength up to I = 0.5 mol ⋅ L −1 . This behavior can be observed for both dissociation steps of IA ( pQ 1 and pQ 2 ) calculated with ePC-SAFT up to an ionic strength of 0.9 mol ⋅ L −1 as well. Above this ionic strength, the concentration quotient of the first dissociation ( pQ 1 ) converges against the equilibrium constant pK a,1 in contrast to the concentrations quotient calculated with the Davies approximation pQ 1,Davies . It can be calculated that pQ 1,Davies reaches pK a,1 not until an ionic strength of 7.3 mol ⋅ L −1 . Concurrently, the difference between pQ 2 and pQ 2,Davies increases above I = 0.9 mol ⋅ L −1 . Therefore, the calculated impact on the dissociation equilibrium differs strongly at high ionic strength ( I > 0.9 mol ⋅ L −1 ) between the Davies and the ePC-SAFT approach.
The impact of the increased ionic strength on the dissociation equilibrium is associated with a change of the distribution of the IA species. Figure 8 shows the relative share of each acid species as a function of pH for each of the two titrations. The black lines indicate the ideal calculation neglecting the activity coefficients, the grey lines represent the calculation applying ePC-SAFT.
The acid constant (Eq. 3) describes the dissociation equilibrium in diluted solutions, in which nonidealities can be neglected ( i = 1 ). Therefore, the intersections of the -curves in the Fig. 8a and 8b mark the acid constants at pH = pK acid,1 = 3.84 and pH = pK acid,2 = 5.45 in case of the ideal calculation method. The acid constants are not valid in concentrated solutions, because the increased ionic strength leads to a shift of these constants [37][38][39], which can be calculated with activity coefficients. The acid constants calculated with ePC-SAFT (grey lines) are shifted to lower pH values. In Fig. 8a, pK acid,1 is marked at pH 3.58 and pK acid,2 at pH 5.40. Therefore, the first dissociation equilibrium is more shifted than the second one. With increasing NaCl concentration (Fig. 8b), the influence of the ionic strength on pK acid,1 decreases and on pK acid,2 increases. The -curves intersect at pH = 3.67 and pH = 5.26. The investigations demonstrate that ePC-SAFT enables the consideration of the ionic strength and the influence on the dissociation equilibrium. As the dissociation equilibrium determines the actual species distribution in solution, the correct description of the dissociation equilibrium is mandatory in the design of downstream processes, e.g. in crystallization or adsorption processes. This can be demonstrated by the pK acid,1 -shift in Fig. 8a: Pure IA can be separated from fermentation broths by crystallization of the protonated form. For the correct calculation of the supersaturation and expected yield in the crystallization process, the correct concentration of the protonated form must be used. Comparing the species distributions of the ideal model and the ePC-SAFT model, it becomes apparent that lower supersaturation and yield is obtained with the ePC-SAFT calculation.

Conclusion
The dynamical approach by Glaser et al. [1] was combined with the ePC-SAFT equation of state to calculate the composition and pH in aqueous itaconic acid solutions at high ionic strengths. A two step determination procedure for the ePC-SAFT parameters of each IA species, based on the works from Ruether [15] et al. and Lange et al. [13], was applied. The calculated solubility of IA in different organic solvents are in very good agreement with the experimental data from Yang et al. [34]. The pH-dependent solubility of IA in water [16] was reproduced with an averaged relative deviation of 4.65 %. The resulting model was then compared with the ideal model and with the Davies approximation to the Debye-Hückel model and validated with titration experiments. A large influence of the ionic strength on measured pH was observed, indicating a shift in the dissociation equilibrium of itaconic acid. This impact on the dissociation equilibrium is similar to the calculations with the Davies approximation carried out by Schell et. al and Zars et al. [3,5] at low ionic strength ( I < 1 mol ⋅ L −1 ). The comparison of the calculation models demonstrated that results obtained with the ePC-SAFT model are in better agreement with experimental titration data at high ionic strengths than the other models. The resulting change of the species distribution of IA compared to the ideal calculation ( i = 1 ) can be relevant for the design of downstream processes such as pHshift-crystallization and adsorption.
Funding Open Access funding enabled and organized by Projekt DEAL. This work was funded by the Bundesministerium fuer Bildung und Forschung (BMBF, Federal Ministry of Education and Research) within the Project BioSorp (FKZ 031B0678A).
Data Availability All data generated or analysed during this study are included in this published article and its supplementary information files.