COMPUTER METHOD AND SYSTEM FOR PREDICTING PHYSICAL PROPERTIES USING A CONCEPTUAL SEGMENT-BASED IONIC ACTIVITY COEFFICIENT MODEL

Abstract

In the present invention the NonRandom Two-Liquid segment activity coefficient model system of the parent application is extended for computation of ionic activity coefficients and solubilities of electrolytes, organic and inorganic, in common solvents and solvent mixtures. The invention method and system may be applied to the chemical and/or pharmaceutical design process. In addition to the three types of molecular parameters defined for organic nonelectrolytes, i.e., hydrophobicity X, polarity Y, and hydrophilicity Z, an electrolyte parameter, E, is introduced to characterize both local and long-range ion-ion and ion-molecule interactions attributed to ionized segments of electrolytes. Successful representations of mean ionic activity coefficients and solubilities of electrolytes, inorganic and organic, in aqueous and nonaqueous solvents are presented.

Description

BACKGROUND OF THE INVENTION

The chemical and pharmaceutical industries screen and develop hundreds of new chemicals and drug candidates each year. Chemists and engineers are tasked to develop process recipes for these new molecules and the recipes often involve multiple reaction steps coupled with separation steps such as crystallization or extraction. A critical consideration in the chemical process design is the choice of solvents and solvent mixtures, from among hundreds of solvent candidates, for reaction, separation, and purification (Frank, T. C. et al., “Quickly Screen Solvents for Organic Solids,” Chemical Engineering Progress, December 1999, 41). Phase behavior, especially solubility, of the new molecules in solvents or solvent mixtures weighs heavily in the choice of solvents for recipe development while little if any such experimental data is available for these new molecules. Although limited solubility experiments may be taken as part of the process research and development, during the early stages of development it is typical for limited experimental resource and drug substance availability, to restrict experimental program of solvent selection for process development. Where the process requires a mixed solvent system it is practically impossible to cover the full range of potential solvent combinations with sufficient detail to find the optimal solution, even with modern high throughput techniques. Consequently, solvent selection today is largely dictated by researchers' preferences or prior experiences. This often leads to a sub optimal manufacturing process reaching the pilot plant and potential manufacturability issues at various steps of the process. To overcome these obstacles, it is highly desirable to have a predictive model of chemical solubility in single and mixed solvent systems, based on a small initial set of measured solubilities.

Existing solubility parameter models such as that of Hansen (Hansen, C. M., Hansen Solubility Parameters: A User's Handbook; CRC Press, 2000) offer very limited predictive power while group contribution models such as UNIFAC (Fredenslund, A. et al., “Group-Contribution Estimation of Activity Coefficients in Nonideal Liquid Mixtures,” AIChE J. 21:1086, 1975) are rather inadequate due to missing functional groups or collapse of functional group additivity rule with large, complex molecules.

Recently Chen and Song (Chen, C.-C. and Y. Song, “Solubility Modeling with NonRandom Two-Liquid Segment Activity Coefficient Model,” Industrial and Engineering Chemistry Research, 43:8354, 2004a) and the related U.S. patent application Ser. No. 10/785,925, proposed a NonRandom Two-Liquid segment activity coefficient (NRTL-SAC) model for fast, qualitative correlation and estimation of solubility of organic nonelectrolytes in common solvents and solvent mixtures. Conceptually, the approach suggests that one could account for the liquid phase nonideality of mixtures of small solvent molecules and complex chemical molecules in terms of pre-defined conceptual segments with pre-determined binary interaction characteristics. Examples of the conceptual segments are hydrophobic segment, polar segment, and hydrophilic segment. The numbers of conceptual segments for each molecule, solvent or solute, reflect the characteristic surface interaction area and nature of the surface interactions. While loosely correlated with molecular structure, they are identified from true behavior of the molecules in solution, i.e., available experimental phase equilibrium data. The molecular make-up in terms of numbers of conceptual segments, i.e. hydrophobicity X, polarity types Y− and Y+, and hydrophilicity Z, constitutes the molecular parameters for the solvent and solute molecules. Given the molecular parameters for solvent and solute molecules, the model offers a thermodynamically consistent expression for estimation of phase behavior, including solubilities, for organic nonelectrolytes in chemical process design.

SUMMARY OF THE INVENTION

The recently developed NRTL-SAC model (parent patent application) proves to be an effective tool for the correlation and prediction of nonelectrolyte molecule solubility in single and mixed solvent systems. The present invention extends the NRTL-SAC model to salts as much of chemical synthesis and medicinal therapy are derived from or administered as salts. The present invention provides a simple and practical qualitative ionic activity coefficient model to support the early stage solvent selection of manufacturing process development for salts and help define the impact of polymorphism from a theoretical perspective.

In some embodiments, the present invention features a method of modeling at least one physical property of a mixture of at least one electrolyte dissolved in one or more solvents. In one embodiment, the method comprises computer implemented steps of determining a conceptual electrolyte segment for the electrolyte, using the determined conceptual electrolyte segment, computing at least one physical property of the mixture, and providing an analysis of the computed physical property. The step of determining a conceptual electrolyte segment includes defining a segment number. The analysis forms a model of at least one physical property of the mixture.

The electrolyte can be organic or inorganic, symmetrical or unsymmetrical, or univalent or multivalent. The electrolyte can be a pharmaceutical compound, a nonpolymeric compound, a polymer or an oligomer. The electrolyte may include two or more ionic species.

In one embodiment, the conceptual electrolyte segment includes a cationic segment and an anionic segment. Both segments are of unity charge. The step of computing at least one physical property includes calculating an ionic activity coefficient of the electrolyte.

In further embodiment, the present invention includes a method of modeling at least one physical property of a mixture of at least one electrolyte dissolved in one or more solvents. The method comprises the computer implemented steps of determining a combination of conceptual segments including at least one conceptual electrolyte segment for the electrolyte, using one or more of the conceptual electrolyte segments, computing at least one physical property of the mixture, and providing an analysis of the computed physical property. For each conceptual segment the step of determining a combination of conceptual segments includes defining a segment number. The analysis forms a model of at least one physical property of the mixture. The conceptual segments can include a polar segment, a hydrophilic segment, a hydrophobic segment and an electrolyte segment.

In another embodiment, this invention features a method of modeling at least one physical property of a mixture of at least one chemical compound dissolved in one or more solvents. The method comprises the computer implemented steps of: determining a conceptual electrolyte segment for the chemical component, using the determined conceptual electrolyte segment, computing at least one physical property of the mixture, and providing an analysis of the computed physical property. The step of determining a conceptual electrolyte segment includes defining a segment number. The analysis forms a model of at least one physical property of the mixture.

In further embodiments, the present invention provides a computer program product or computer system carrying out the above invention methods.

In one embodiment, the computer system comprises a user input means for obtaining or defining chemical data from a user, a digital processor coupled to receive the chemical data input from the input means, and an output means coupled to the digital processor. The output means provides to the user the formed model of the physical property of a mixture. The digital processor executes a modeling system in working memory. The modeling system: (i) uses the chemical data to determine a conceptual electrolyte segment for a subject electrolyte, including defining a segment number; (ii) uses the determined conceptual electrolyte segment to compute at least one physical property of a mixture formed of the subject electrolyte; and (iii) provides an analysis of the computed physical property. The analysis forms a model of at least one physical property of the mixture.

In some embodiments, the present invention features a chemical and/or pharmaceutical compound manufactured by a process that employs the above modeling methods.

In some embodiments, the present invention features a chemical or pharmaceutical compound comprising at least one electrolyte. The electrolyte is selected by a model of the modeling methods and systems above. The model determines a conceptual electrolyte segment of the electrolyte and computes physical properties of the electrolyte.

Accordingly, the present invention provides for the fast, practical modeling of physical properties or behaviors of mixtures of at least one electrolyte dissolved in one or more solvents, even when there is little or no experimental data to which the behavior of the mixture can be correlated. The formed models offer improved accuracy over most or all prior modeling methods. For example, this invention offers a simple and practical tool for practitioners to estimate solubility of various components of a mixture (e.g., a mixture including a chemical compound), even when there is little or no phase equilibrium data available for the mixture.

This invention provides for modeling of mixtures having significant electrolyte interactions, hydrophobic interactions, polar interactions, and/or hydrophilic interactions. This invention eliminates the need to characterize mixture constituents with sets of pre-defined functional groups and provides for the modeling of mixtures comprising large, complex electrolytes for which a functional group additivity rule becomes invalid and/or for which there are a number of un-defined functional groups.

BRIEF DESCRIPTION OF THE DRAWINGS

The foregoing and other objects, features and advantages of the invention will be apparent from the following more particular description of preferred embodiments of the invention, as illustrated in the accompanying drawings in which like reference characters refer to the same parts throughout the different views. The drawings are not necessarily to scale, emphasis instead being placed upon illustrating the principles of the invention.

FIG. 1 is a graph illustrating the effect of hydrophobicity parameter X on natural logarithm of mean ionic activity coefficient of aqueous electrolytes with E=1.

FIG. 2 is a graph illustrating the effect of polarity parameter Y− on natural logarithm of mean ionic activity coefficient of aqueous electrolytes with E=1.

FIG. 3 is a graph illustrating the effect of polarity parameter Y+on natural logarithm of mean ionic activity coefficient of aqueous electrolytes with E=1.

FIG. 4 is a graph illustrating the effect of hydrophilicity parameter Z on natural logarithm of mean ionic activity coefficient of aqueous electrolytes with E=1.

FIG. 5 is a graph illustrating the effect of electrolyte parameter E on natural logarithm of mean ionic activity coefficient of aqueous electrolytes.

FIG. 6 is a graph illustrating comparison of experimental and calculated molality scale mean ionic activity coefficients of representative aqueous electrolytes at 298.15 K.

FIG. 7 is a graph illustrating the present invention model results for sodium chloride solubility at 298.15 K.

FIG. 8 is a graph illustrating the present invention model results for sodium acetate solubility at 298.15 K.

FIG. 9 a is a graph illustrating the present invention model results for benzoic acid solubility at 298.15 K.

FIG. 9 b is a graph illustrating the present invention model results for sodium benzoate solubility at 298.15 K.

FIG. 10 a is a graph illustrating the present invention model results for salicylic acid solubility at 298.15 K.

FIG. 10 b is a graph illustrating the present invention model results for sodium salicylate solubility at 298.15 K.

FIG. 11 a is a graph illustrating the present invention model results for p-aminobenzoic acid solubility at 298.15 K.

FIG. 11 b is a graph illustrating the present invention model results for sodium p-aminobenzoate solubility at 298.15 K.

FIG. 12 a is a graph illustrating the present invention model results for ibuprofen solubility at 298.15 K FIG. 12b is a graph illustrating the present invention model results for sodium ibuprofen solubility at 298.15 K.

FIG. 13 a is a graph illustrating the present invention model results for diclofenac solubility at 298.15 K.

FIG. 13 b is a graph illustrating the present invention model results for sodium diclofenac solubility at 298.15 K.

FIG. 14 is a schematic view of a computer network in which the present invention may be implemented.

FIG. 15 is a block diagram of a computer of the network of FIG. 14.

FIGS. 16-17 are flow diagrams of one embodiment of the present invention employed in the computer network environment of FIGS. 14 and 15.

DETAILED DESCRIPTION OF THE INVENTION

Electrolytes are substances which exist as free ions or dissociate into free ions when dissolved to produce an electrically conductive medium. Electrolytes generally exist as salts. Salts have a vast array of chemical and medicinal applications. Predicting the phase behaviors of the electrolytes/salts can assist chemists in the design and in the synthesis of a molecule. For example, it is estimated that half of all drug molecules used in medicinal therapy are administered as salts (Stahl, P. H. and C. F. Wermuth, (Eds.), Handbook of Pharmaceutical Salts: Properties, Selection and Use, Verlag Helvetica Chimica Acta, Zürich, Switzerland (2002)). The NRTL-SAC model (parent application/invention) for organic nonelectrolytes is extended to account for the liquid phase nonideality due to presence of ionic charges in organic electrolytes. Applicants have developed a simple and practical qualitative ionic activity coefficient model to aid in designing the chemical process for salts. The extension represents a major advance to the scope of the NRTL-SAC thermodynamic framework. In the absence of organic moiety, organic electrolytes become inorganic electrolytes. In the absence of ionized moiety, organic electrolytes become organic nonelectrolytes. The electrolyte extension of NRTL-SAC model (present invention eNRTL-SAC) provides a consistent and comprehensive thermodynamic framework for systems involving all types of electrolytes and nonelectrolytes. In other words, the NRTL-SAC model of the parent invention/application becomes a limiting case for the present invention eNRTL-SAC model.

NRTL Segment Activity Coefficient Model

The NRTL-SAC activity coefficient model for component I is composed of the combinatorial term γ_I^C and the residual term γ_I^R:

lnγ_I=lnγ_I^C+lnγ_I^R  (1)

Here the combinatorial term γ_I^C is calculated from the Flory-Huggins equation for the combinatorial entropy of mixing. The residual term γ_I^R is calculated from the local composition (lc) interaction contribution γ_I^lc of Polymer NRTL (Chen, C.-C., “A Segment-Based Local Composition Model for the Gibbs Energy of Polymer Solutions,” Fluid Phase Equilibria, 83:301, 1993) (herein “Chen 1993”). The Polymer NRTL equation incorporates the segment interaction concept and computes activity coefficient for component I in a solution by summing up contributions to activity coefficient from all segments that make up component I. The equation is given as follows:

$\begin{array}{cc}\ln {\gamma }_I^R=\ln {\gamma }_I^{\mathrm{lc}}=\sum _mr_{m,I}\left[\ln {\Gamma }_m^{\mathrm{lc}}-\ln {\Gamma }_m^{\mathrm{lc},I}\right]\mathrm{with}& \left(2\right)\\ \ln {\Gamma }_m^{\mathrm{lc}}=\frac{\sum _jx_jG_{\mathrm{jm}}{\tau }_{\mathrm{jm}}}{\sum _kx_kG_{km}}+\sum _{m^{\prime }}\frac{x_{m^{\prime }}G_{mm^{\prime }}}{\sum _kx_kG_{km^{\prime }}}\left(\begin{array}{c}{\tau }_{{\mathrm{mm}}^{\prime }}-\\ \frac{\sum _jx_jG_{{\mathrm{jm}}^{\prime }}{\tau }_{{\mathrm{jm}}^{\prime }}}{\sum _kx_kG_{{\mathrm{km}}^{\prime }}}\end{array}\right)& \left(3\right)\\ \ln {\Gamma }_m^{\mathrm{lc},I}=\frac{\sum _jx_{j,I}G_{\mathrm{jm}}{\tau }_{\mathrm{jm}}}{\sum _kx_{k,I}G_{km}}+\sum _{m^{\prime }}\frac{x_{m^{\prime },I}G_{mm^{\prime }}}{\sum _kx_{k,I}G_{km^{\prime }}}\left(\begin{array}{c}{\tau }_{{\mathrm{mm}}^{\prime }}-\\ \frac{\sum _jx_{j,I}G_{{\mathrm{jm}}^{\prime }}{\tau }_{{\mathrm{jm}}^{\prime }}}{\sum _kx_{k,I}G_{{\mathrm{km}}^{\prime }}}\end{array}\right)& \left(4\right)\\ x_j=\frac{\sum _Ix_Ir_{j,I}}{\sum _I\sum _ix_Ir_{i,I}}x_{j,I}=\frac{r_{j,I}}{\sum _ir_{i,I}}& \left(5\right)\end{array}$

where I is the component index, i, j, k, m, m′ are the segment species index, x_l is the mole fraction of component I, x_j is the segment-based mole fraction of segment species j, r_m,l is the number of segment species m contained only in component I, Γ_m^lc is the activity coefficient of segment species m, and Γ_m^lc,I is the activity coefficient of segment species m contained only in component I. G and τ in Eqs. 3 and 4 are local binary quantities related to each other by the NRTL non-random factor parameter α:

G=exp(−ατ).  (6)

Four pre-defined conceptual segments were suggested by Chen and Song (2004 above and in parent patent application): one hydrophobic (x), two polar (y− and y+), and one hydrophilic (z). The model molecular parameters, i.e., hydrophobicity X, polarity types Y− and Y+, and hydrophilicity Z, correspond to r_m,l (m=x, y−, y+, z), numbers of various conceptual segments in component I.

In the notation used throughout this disclosure, subscript I (upper case) refers to components while subscript i (lower case) refers to segments.

eNRTL Segment Activity Coefficient Model

The extension of NRTL-SAC model for electrolytes is based on the generalized eNRTL model as summarized by Chen and Song (Chen, C.-C. and Y. Song, “Generalized Electrolyte NRTL Model for Mixed-Solvent Electrolyte Systems,” AIChE J., 50:1928, 2004b herein incorporated by reference) (herein “Chen, 2004b). Here Applicants briefly present the generalized eNRTL model followed by details of the extended NRTL-SAC model of the present invention.

The generalized eNRTL model is applied to correlate mean ionic activity coefficient of mixed-solvent electrolyte systems. The segment interaction concept provides the framework to explicitly account for the attractive interaction of ions with the hydrophilic segments of organic solvents and the repulsive interaction of ions with the hydrophobic segments of organic solvents. In the generalized eNRTL model, any component, electrolyte or solvent, can be defined as an oligomer consisting of various segment species. For instance, an organic electrolyte species can be defined as an oligomer consisting of cationic segment, anionic segment and molecular segment. An organic solvent can be also defined as an oligomer consisting of multiple molecular segments of different nature. Accordingly, with the conventional activity coefficient accounting for the local interaction (Chen, 1993) and the long-range interaction, the model that uses the unsymmetric Pitzer-Debye-HUckel (PDH) formula (Pitzer, K. S., “Electrolytes: From Dilute Solutions to Fused Salts,” J. Am. Chem. Soc., 102, 2902 (1980)) (herein “Pitzer, 1980”) is calculated as follows:

$\begin{array}{cc}\begin{array}{c}\ln {\gamma }_I^{*}=\frac{1}{\mathrm{RT}}{\left(\frac{\partial G_m^{*\mathrm{ex}}}{\partial n_I}\right)}_{T,P,n_{J\ne I}}\\ =\frac{1}{\mathrm{RT}}{\left(\frac{\partial G_m^{*\mathrm{ex},\mathrm{lc}}}{\partial n_I}\right)}_{T,P,n_{J\ne I}}+\frac{1}{\mathrm{RT}}{\left(\frac{\partial G_m^{*\mathrm{ex},\mathrm{PDH}}}{\partial n_I}\right)}_{T,P,n_{J\ne I}}\end{array}\mathrm{Or}& \left(7\right)\\ \ln {\gamma }_I^{*}=\ln {\gamma }_I^{*\mathrm{lc}}+\ln {\gamma }_I^{*\mathrm{PDH}}& \left(8\right)\end{array}$

where I is the component index, “*” denotes the unsymmetric convention, γ_I is the activity coefficient of the component I in the mixture; R is the gas constant; T is the temperature; P is the pressure; and n_I is the mole number of the component I in the mixture. The unsymmetric PDH formula, G*_m^ex,PDH, is obtained by normalization to mole fractions of unity for solvents and zero for electrolytes (Pitzer, K. S., “Thermodynamics of Electrolytes. I: Theoretical and General Equations,” J. Phys. Chem., 77, 268 (1973)). The local interaction NRTL model, G*_m^ex,lc, is developed as a symmetric model (Chen, C.-C., “A Segment-Based Local Composition Model for the Gibbs Energy of Polymer Solutions,” Fluid Phase Equilib., 83, 301 (1993); and Chen, C.-C., C. P. Bokis, and P. M. Mathias, “A Segment-Based Excess Gibbs Energy Model for Aqueous Organic Electrolyte Systems,” AIChE J., 47, 2593 (2001)), based on the symmetrical reference state so that the derived activity coefficient, is γ_I^lc=1 as x→1 for any component (species). The model is then normalized by the unsymmetric reference state (that is, the infinite-dilution activity coefficient in an aqueous or mixed-solvent solution) to obtain the unsymmetric model, G*_m^ex,lc. Accordingly, the unsymmetric convention activity coefficient is calculated as follows:

$\begin{array}{cc}\ln {\gamma }_I^{*\mathrm{lc}}=\ln {\gamma }_I^{\mathrm{lc}}-\ln {\gamma }_I^{\infty \mathrm{lc}}& \left(9\right)\\ \ln {\gamma }_I^c=\frac{1}{\mathrm{RT}}{\left(\frac{\partial G_m^{\mathrm{ex},\mathrm{lc}}}{\partial n_I}\right)}_{T,P,n_{J\ne I}}& \left(10\right)\end{array}$

where γ_I^∞ is the infinite-dilution activity coefficient of the ionic component I in an aqueous or mixed-solvent solution as calculated by Equation 10. A more detailed description on the generalized electrolyte-NRTL model is depicted in Chen, 2004b.

This generalized segment interaction concept is advantageous when one must exactly account for the different interaction characteristics that may be attributed to different molecules, solvents or solutes. The ability to exactly account for such different segment-segment interactions between different species in a system is shown to be key for quantitative correlation of mean ionic activity coefficients in mixed-solvent electrolyte systems. In the generalized eNRTL model, however, it is necessary to account for an electrolyte segment for each and every species separately. Therefore, in a system that involves multiple components, there could be tens of different segments to consider and hundreds of segment-segment interactions to account for, and the computation for activity coefficients becomes much more complicated.

Derived from and improved upon the generalized eNRTL model, the electrolyte extension of NRTL model of the present invention provides one conceptual electrolyte segment. A “conceptual electrolyte segment” herein is one predefined electrolyte segment that characterizes the prominent interaction mechanisms between molecules in the liquid phase, that account for the liquid phase nonideality. This pre-defined electrolyte segment is used as a reference against which all electrolyte segments are measured in terms of their liquid phase interaction characteristics. Unlike the generalized eNRTL model, which has no such “conceptual electrolyte segment” as a reference point, surface interaction characteristics of electrolyte segments of the present invention are normalized against the “conceptual electrolyte segment” (in a preferred embodiment, one with interaction characteristics of NaCl) and mathematically expressed as an equivalent number of the reference one. Having a point of reference for the calculation of the electrolyte segment provides a unified and consistent description of liquid phase nonideality of all electrolyte segments and a more intuitive and powerful predictive tool in modeling physical properties including solubility. Together with the numbers of “conceptual” hydrophobic segment, hydrophilic segment and polar segment, the number of “conceptual electrolyte segment” reflects the nature of the surface interactions and their characteristic surface interaction areas that determine their phase behavior.

Electrolytes dissociate to ionic species in solutions. For “strong” electrolytes, the dissociation is “completely” to ionic species. For “weak” electrolytes, the dissociation is partially to ionic species while undissociated electrolytes, similar to nonelectrolytes, remain as neutral molecular species. Complexation of ionic species with solvent molecules or other ionic species may also occur. An implication of the electrolyte solution chemistry is that the extended model should provide a thermodynamically consistent framework to compute activity coefficients for both molecular species and ionic species.

In the simplest case of a strong electrolyte CA, one may use the following chemical reaction to describe the complete dissociation of the electrolyte:

CA→ν _C C ^{z C} +ν_AA^zA  (11)

with

ν_CZ_V=ν_AZ_A  (12)

where ν_C is the cationic stoichiometric coefficient, ν_A is the anionic stoichiometric coefficient, Z_C is the absolute charge number for cation C, and Z_A is the absolute charge number for anion A.

In applying the segment contribution concept to electrolytes, Applicants introduce a new conceptual electrolyte segment e. This conceptual segment e would completely dissociate to a cationic segment (c) and an anionic segment (a), both of unity charge. Applicants then follow the like-ion repulsion and the electroneutrality constraints imposed by the generalized eNRTL model to derive the activity coefficient equations for ionic segments c and a. All electrolytes, organic or inorganic, symmetric or unsymmetric, univalent or multivalent, are to be represented with this conceptual uni-univalent electrolyte segment e together with previously defined hydrophobic segment, x, polar segments, y− and y+, and hydrophilic segment, z. Due to the fact that Applicants introduce only one (a universally useable one) conceptual electrolyte segment e, the resulting eNRTL-SAC model of the present invention is much simpler than the generalized eNRTL model proposed earlier.

In some embodiments, this invention includes methods and apparatus/systems of modeling at least one physical property of a mixture of at least one electrolyte dissolved in one or more solvents. Preferably, the method, apparatus and system comprise the computer implemented steps of: (a) determining a conceptual electrolyte segment for the electrolyte, including defining a segment number; (b) using the determined conceptual electrolyte segment, computing at least one physical property of the mixture; and (c) providing an analysis of the computed physical property. The analysis forms a model of the at least one physical property of the mixture.

The invention method, apparatus/system models mixtures of one or more electrolytes that are organic or inorganic, symmetrical or unsymmetrical, or univalent or multivalent. The electrolyte may be a chemical compound, a pharmaceutical compound, a nonpolymeric compound, a polymer, an oligomer, an inorganic compound and an organic compound. The electrolyte may include two or more ionic species.

As used herein, a “pharmaceutical compound” can include drugs, therapeutic agents, salts or a precursor thereof (i.e., a compound used as an ingredient in a pharmaceutical compound production process) and all manner of compounds residing within an organism or system, such as nutrients, metabolites, endogenous hormones, toxins, and the likes. In some embodiments, the chemical compound has a weight greater than about 900 daltons, at least one molecule having a molecular weight in the range of between about 100 daltons and about 900 daltons, and/or at least one molecule having a molecular weight in the range of between about 200 daltons and about 600 daltons.

In accordance with one aspect of the present invention, at least one conceptual segment (e.g., at least 1, 2, 3, 4, or more than 4 conceptual segments) is determined or defined for each of the chemical species of a subject mixture. The conceptual segments are molecular descriptors of the various molecular species in the mixture. An identity and an equivalent number are determined for each of the conceptual segments. Examples of identities for conceptual segments include a hydrophobic segment, a polar segment, a hydrophilic segment, an electrolyte segment, and the like. Experimental phase equilibrium data can be used to determine the equivalent number of the conceptual segment(s).

The determined conceptual segments are used to compute at least one physical property of the mixture, and an analysis of the computed physical property is provided to form a model of at least one physical property of the mixture. The present invention enables a wide variety of physical properties to be modeled. Examples of physical properties include an activity coefficient, vapor pressure, solubility (e.g., the equilibrium concentration of one or more chemical species in one or more phases of the mixture), boiling point, freezing point, octanol/water partition coefficient, lipophilicity of the electrolyte and other physical properties that are measured or determined for use in the chemical processes.

The models of the physical property or properties of the mixture are produced by determining the interaction characteristics of the conceptual segments. In some embodiments, the segment-segment interaction characteristics of the conceptual segments are represented by their corresponding binary eNRTL parameters. Given the eNRTL parameters for the conceptual segments and the molecular descriptors for the molecules, the eNRTL-SAC model computes activity coefficients for the segments and then for the various electrolytes in the mixture. In other words, the physical properties or behavior of the mixture will be accounted for based on the segment compositions of the electrolytes and their mutual interactions. In further embodiments, the step of computing at least one physical property includes calculating an activity coefficient of the electrolyte.

For example, the solubility of an electrolyte is described well by the expression:

$\begin{array}{cc}{K}_{\mathrm{sp}}\left(T\right)=\prod _Cx_C^{v_c,\mathrm{SAT}}{\gamma }_C^{*v_c,\mathrm{SAT}}\prod _Ax_A^{v_A,\mathrm{SAT}}{\gamma }_A^{*v_A,\mathrm{SAT}}\prod _Mx_M^{\mathrm{SAT}}{\gamma }_M^{\mathrm{SAT}},& \left(13\right)\end{array}$

where K_sp is the solubility product constant for the electrolyte, T is the temperature of the mixture, x_C^vCSAT is the mole fraction of a cation derived from the electrolyte at saturation point of the electrolyte, x^vASAT is the mole fraction of an anion derived from the electrolyte at saturation point of the electrolyte, x_M^vMSAT is the mole fraction of a neutral molecule derived from the electrolyte at saturation point of the electrolyte, γ*_C^vC,SAT is the activity coefficient of a cation derived from the electrolyte at the saturation concentration, γ*_A^vA,SAT is the activity coefficient of an anion derived from the electrolyte at the saturation concentration, γ*_M^vM,SAT is the activity coefficient of a neutral molecule derived from the electrolyte at the saturation concentration, C is the cation, γ* is the unsymmetric activity coefficient of a species in solution, SAT is saturation concentration, v_C is the cationic stoichiometric coefficient, D_A is the anionic stoichiometric coefficient, and v_m is the neutral stoichiometric coefficient.

A major consideration in the extension of NRTL-SAC for electrolytes is the treatment of reference state for activity coefficient calculations. While the conventional reference state for nonelectrolyte systems is the pure liquid component, the conventional reference state for electrolytes in solution is the infinite-dilution aqueous solution and the corresponding activity coefficient is “unsymmetric.”

Following the generalized eNRTL model, the logarithm of unsymmetric activity coefficient of an ionic species, lnγ*_I, is the sum of three terms: the local composition term, lnγ*_I^lc, the Pitzer-Debye-Hückel term, lnγ*_I^PDH, and the Flory-Huggins term, lnγ*_I^FH.

lnγ*_I=lnγ*_I^lc+lnγ*_I^PDH+lnγ*_I^FH  (14)

Eq. 14 applies to aqueous electrolyte systems where water is a sole solvent within the solution. For mixed-solvent solutions, the Born term, Δlnγ_I^Born, is used to correct the change of the infinite dilution reference state from the mixed-solvent composition to the aqueous solution for the Pitzer-Debye-Hückel term:

lnγ*_I=lnγ*_I^lc+lnγ*_I^PDH+lnγ*_I^FH+Δlnγ_I^Born  (15)

Since Applicants adopt the aqueous phase infinite dilution reference state for γ*_I, the Born term correction is required for non-aqueous systems.

With the introduction of the conceptual electrolyte segment e and the corresponding conceptual ionic segments c and a, one can rewrite Eq. 15 in terms of contributions from all conceptual segments:

$\begin{array}{cc}\begin{array}{c}\ln {\gamma }_I^{*}=\ln {\gamma }_I^{*\mathrm{lc}}+\ln {\gamma }_I^{*\mathrm{PDH}}+\ln {\gamma }_I^{*\mathrm{FH}}+\Delta \ln {\gamma }_I^{\mathrm{Born}}\\ =\sum _mr_{m,I}\left(\ln {\Gamma }_m^{*\mathrm{lc}}+\ln {\Gamma }_m^{*\mathrm{PDH}}\right)+\\ r_{c,I}\left(\ln {\Gamma }_c^{*\mathrm{lc}}+\ln {\Gamma }_c^{*\mathrm{PDH}}+{\mathrm{\Delta ln\Gamma }}_c^{\mathrm{Born}}\right)+\\ r_{a,I}\left(\ln {\Gamma }_a^{*\mathrm{lc}}+\ln {\Gamma }_a^{*\mathrm{PDH}}+{\mathrm{\Delta ln\Gamma }}_a^{\mathrm{Born}}\right)+\ln {\gamma }_I^{*\mathrm{FH}}\end{array}& \left(16\right)\end{array}$

where r is the segment number, m is the conceptual molecular segment index y−, y+, z), c and a are cationic and anionic segments, respectively, resulting from the dissociation of the conceptual electrolyte segment e. Also notice that in Eq. 16, unlike the local composition term and the long range ion-ion interaction terms, the Flory-Huggins term remains as the component-based contribution.

For systems of single electrolyte CA with a segment number r_e, r_c and r_a must satisfy electroneutrality and they can be computed from r_e, Z_C, and Z_A.

r_c,C=r_e,CAZ_C  (17)

r_a,A=r_e,CAZ_A  (18)

For systems of multiple electrolytes, the mixing rule is needed to compute segment number r_c and r_a for each cation C and anion A.

$\begin{array}{cc}{r}_{c,C}=\sum _Ar_{e,\mathrm{CA}}Z_C\left(x_AZ_A/\sum _{A^{\prime }}x_{A^{\prime }}Z_{A^{\prime }}\right)& \left(19\right)\\ r_{a,C}=\sum _Cr_{e,\mathrm{CA}}Z_A\left(x_CZ_C/\sum _{C^{\prime }}x_{C^{\prime }}Z_{C^{\prime }}\right)& \left(20\right)\end{array}$

r_e,CA, the number of conceptual electrolyte segment e in electrolyte CA, becomes the new model parameter for electrolytes. For the sake of brevity, Applicants call r_e,CA parameter E, the electrolyte segment number.

Local Composition Interaction Contribution

To derive the expression for the local composition interaction contribution, Applicants simplify the generalized excess Gibbs energy expression of the prior Chen and Song model (Chen, 2004b) for systems with multiple molecular segments m and single electrolyte segment e. The single electrolyte segment e is then decomposed into a cationic segment c and an anionic segment a:

$\begin{array}{cc}\frac{G^{\mathrm{ex},\mathrm{lc}}}{\mathrm{nRT}}=\sum _I\left[\begin{array}{c}\sum _mr_{m,I}x_I\left(\frac{\sum _jx_jG_{\mathrm{jm}}{\tau }_{\mathrm{jm}}}{\sum _kx_kG_{km}}\right)+\\ r_{c,I}x_I\left(\frac{\sum _jx_jG_{\mathrm{jc},\mathrm{ac}}{\tau }_{\mathrm{jc},\mathrm{ac}}}{\sum _kx_kG_{\mathrm{kc},\mathrm{ac}}}\right)+r_{a,I}x_I\left(\frac{\sum _jx_jG_{\mathrm{ja},\mathrm{ca}}{\tau }_{\mathrm{ja},\mathrm{ca}}}{\sum _kx_kG_{\mathrm{ka},\mathrm{ca}}}\right)\end{array}\right]\mathrm{with}& \left(21\right)\\ x_j=\frac{\sum _Ix_Ir_{j,I}}{\sum _I\sum _ix_Ir_{i,I}}i,j=m,c,a& \left(22\right)\end{array}$

where G^ex,lc is the excess Gibbs energy from local composition interactions, n is the total mole number, R is the gas constant and T is the temperature.

To derive the segment activity coefficient, one can rewrite Eq. 21 as follows:

$\begin{array}{cc}\frac{G^{\mathrm{ex},\mathrm{lc}}}{n_S\mathrm{RT}}=\sum _mx_m\left(\frac{\sum _jx_jG_{\mathrm{jm}}{\tau }_{\mathrm{jm}}}{\sum _kx_kG_{\mathrm{km}}}\right)+x_c\left(\frac{\sum _jx_jG_{\mathrm{jc},\mathrm{ac}}{\tau }_{\mathrm{jc},\mathrm{ac}}}{\sum _kx_kG_{\mathrm{kc},\mathrm{ac}}}\right)+x_a\left(\frac{\sum _jx_jG_{\mathrm{ja},\mathrm{ca}}{\tau }_{\mathrm{ja},\mathrm{ca}}}{\sum _kx_kG_{\mathrm{ka},\mathrm{ca}}}\right)& \left(23\right)\end{array}$

where n_S is the total number of all segments. Accordingly, the segment activity coefficient can be calculated as follows:

$\begin{array}{cc}\begin{array}{cc}\ln {\Gamma }_j^{\mathrm{lc}}=\frac{1}{\mathrm{RT}}{\left(\frac{\partial G^{\mathrm{ex},\mathrm{lc}}}{\partial n_j}\right)}_{T,P,n_{i\ne j}}& i,j=m,c,a\end{array}& \left(24\right)\end{array}$

Specifically, the activity coefficients from Eq. 24 for molecular segments, cationic segment, and anionic segment can be carried out as follows:

$\begin{array}{cc}\ln {\Gamma }_m^{\mathrm{lc}}=\frac{\sum _jx_jG_{\mathrm{jm}}{\tau }_{\mathrm{jm}}}{\sum _kx_kG_{\mathrm{km}}}+\sum _{m^{\prime }}\frac{x_{m^{\prime }}G_{{\mathrm{mm}}^{\prime }}}{\sum _kx_kG_{{\mathrm{km}}^{\prime }}}\left({\tau }_{{\mathrm{mm}}^{\prime }}-\frac{\sum _jx_jG_{{\mathrm{jm}}^{\prime }}{\tau }_{{\mathrm{jm}}^{\prime }}}{\sum _kx_kG_{{\mathrm{km}}^{\prime }}}\right)+\frac{x_cG_{\mathrm{mc},\mathrm{ac}}}{\sum _kx_kG_{\mathrm{kc},\mathrm{ac}}}\left({\tau }_{\mathrm{mc},\mathrm{ac}}-\frac{\sum _jx_jG_{\mathrm{jc},\mathrm{ac}}{\tau }_{\mathrm{jc},\mathrm{ac}}}{\sum _kx_kG_{\mathrm{kc},\mathrm{ac}}}\right)+\frac{x_aG_{\mathrm{ma},\mathrm{ca}}}{\sum _kx_kG_{\mathrm{ka},\mathrm{ca}}}\left({\tau }_{\mathrm{ma},\mathrm{ca}}-\frac{\sum _jx_jG_{\mathrm{ja},\mathrm{ca}}{\tau }_{\mathrm{ja},\mathrm{ca}}}{\sum _kx_kG_{\mathrm{ka},\mathrm{ca}}}\right)& \left(25\right)\\ \ln {\Gamma }_c^{\mathrm{lc}}=\sum _m\frac{x_mG_{\mathrm{cm}}}{\sum _kx_kG_{\mathrm{km}}}\left({\tau }_{\mathrm{cm}}-\frac{\sum _jx_jG_{\mathrm{jm}}{\tau }_{\mathrm{jm}}}{\sum _kx_kG_{\mathrm{km}}}\right)+\frac{\sum _jx_jG_{\mathrm{jc},\mathrm{ac}}{\tau }_{\mathrm{jc},\mathrm{ac}}}{\sum _kx_kG_{\mathrm{kc},\mathrm{ac}}}-\frac{x_a}{\sum _kx_kG_{\mathrm{ka},\mathrm{ca}}}\left(\frac{\sum _jx_jG_{\mathrm{ja},\mathrm{ca}}{\tau }_{\mathrm{ja},\mathrm{ca}}}{\sum _kx_kG_{\mathrm{ka},\mathrm{ca}}}\right)& \left(26\right)\\ \ln {\Gamma }_a^{\mathrm{lc}}=\sum _m\frac{x_mG_{\mathrm{am}}}{\sum _kx_kG_{\mathrm{km}}}\left({\tau }_{\mathrm{am}}-\frac{\sum _jx_jG_{\mathrm{jm}}{\tau }_{\mathrm{jm}}}{\sum _kx_kG_{\mathrm{km}}}\right)+\frac{\sum _jx_jG_{\mathrm{ja},\mathrm{ca}}{\tau }_{\mathrm{ja},\mathrm{ca}}}{\sum _kx_kG_{\mathrm{ka},\mathrm{ca}}}-\frac{x_c}{\sum _kx_kG_{\mathrm{kc},\mathrm{ac}}}\left(\frac{\sum _mx_mG_{\mathrm{mc},\mathrm{ac}}{\tau }_{\mathrm{mc},\mathrm{ac}}}{\sum _kx_kG_{\mathrm{kc},\mathrm{ac}}}\right)& \left(27\right)\end{array}$

The local composition term for the logarithm of activity coefficient of component I is computed as the sum of the individual segment contributions.

$\begin{array}{cc}\begin{array}{c}\ln {\gamma }_I^{\mathrm{lc}}=\sum _ir_{i,I}\ln {\Gamma }_i^{\mathrm{lc}}\mathrm{ }i=m,c,a\\ \sum _mr_{m,I}\ln {\Gamma }_m^{\mathrm{lc}}+r_{c,I}\ln {\Gamma }_c^{\mathrm{lc}}+r_{a,I}\ln {\Gamma }_a^{\mathrm{lc}}\end{array}& \left(28\right)\end{array}$

However, the activity coefficient by Eq. 28 needs to be further normalized so that γ_I^lc=1 as x_I→1 for any component; this is the so-called symmetric reference state. The normalization can be done as follows:

$\begin{array}{cc}\begin{array}{c}\ln {\gamma }_I^{\mathrm{lc}}=\sum _ir_{i,I}\lfloor \ln {\Gamma }_i^{\mathrm{lc}}-\ln {\Gamma }_i^{\mathrm{lc},I}\rfloor i=m,c,a\\ =\sum _mr_{m,I}\left[\ln {\Gamma }_m^{\mathrm{lc}}-\ln {\Gamma }_m^{\mathrm{lc},I}\right]+r_{c,I}\left[\ln {\Gamma }_c^{\mathrm{lc}}-\ln {\Gamma }_{c}^{\mathrm{lc},I}\right]+\\ r_{a,I}\left[\ln {\Gamma }_a^{\mathrm{lc}}-\ln {\Gamma }_a^{\mathrm{lc},I}\right]\end{array}& \left(29\right)\end{array}$

Here Γ_i^lc,I is the activity coefficient of the segment i contained in the symmetric reference state of component I; it can be calculated from Eqs. 25-27 by setting x_I=1:

lnΓ_i^lc,I=lnΓ_i^lc(x_I=1) i=m,c,a  (30)

Finally, the unsymmetric convention in Eq. 15 requires us to compute the infinite-dilution activity coefficient, γ_I^∞lc, for a component:

lnγ*_I^lc=lnγ_I^lc−lnγ_I^∞lc  (31)

with

$\begin{array}{cc}\begin{array}{c}\ln {\gamma }_I^{\infty \mathrm{lc}}=\sum _ir_{i,I}\lfloor \ln {\Gamma }_i^{\infty \mathrm{lc}}-\ln {\Gamma }_i^{\mathrm{lc},I}\rfloor i=m,c,a\\ =\sum _mr_{m,I}\left[\ln {\Gamma }_m^{\infty \mathrm{lc}}-\ln {\Gamma }_m^{\mathrm{lc},I}\right]+r_{c,I}\left[\ln {\Gamma }_c^{\infty \mathrm{lc}}-\ln {\Gamma }_c^{\mathrm{lc},I}\right]+\\ r_{a,I}\left[\ln {\Gamma }_a^{\infty \mathrm{lc}}-\ln {\Gamma }_a^{\mathrm{lc},I}\right]\end{array}& \left(32\right)\end{array}$

Combining Eqs. 29 and 32, one can obtain:

$\begin{array}{cc}\begin{array}{c}\ln {\gamma }_I^{*\mathrm{lc}}=\ln {\gamma }_I^{\mathrm{lc}}-\ln {\gamma }_I^{\infty \mathrm{lc}}\\ =\sum _ir_{i,I}\left[\ln {\Gamma }_i^{\mathrm{lc}}-\ln {\Gamma }_i^{\infty \mathrm{lc}}\right]i=m,c,a\\ =\sum _mr_{m,I}\left[\ln {\Gamma }_m^{\mathrm{lc}}-\ln {\Gamma }_m^{\infty \mathrm{lc}}\right]+r_{c,I}\left[\ln {\Gamma }_c^{\mathrm{lc}}-\ln {\Gamma }_c^{\infty \mathrm{lc}}\right]+\\ r_{a,I}\left[\ln {\Gamma }_a^{\mathrm{lc}}-\ln {\Gamma }_a^{\infty \mathrm{lc}}\right]\\ =\sum _mr_{m,I}\ln {\Gamma }_m^{*\mathrm{lc}}+r_{c,I}\ln {\Gamma }_c^{*\mathrm{lc}}+r_{a,I}\ln {\Gamma }_a^{*\mathrm{lc}}\end{array}\mathrm{with}& \left(33\right)\\ \ln {\Gamma }_m^{*\mathrm{lc}}=\ln {\Gamma }_m^{\mathrm{lc}}-\ln {\Gamma }_m^{\infty \mathrm{lc}}& \left(34\right)\\ \ln {\Gamma }_c^{*\mathrm{lc}}=\ln {\Gamma }_c^{\mathrm{lc}}-\ln {\Gamma }_c^{\infty \mathrm{lc}}& \left(35\right)\\ \ln {\Gamma }_a^{*\mathrm{lc}}=\ln {\Gamma }_a^{\mathrm{lc}}-\ln {\Gamma }_a^{\infty \mathrm{lc}}& \left(36\right)\end{array}$

Because Applicants adopt the aqueous phase infinite dilution reference state, the infinite-dilution activity coefficients of conceptual segments can be calculated from Eqs. 25-28 by setting x_W=1:

lnΓ_i^∞lc=lcΓ_i^lc(x_W=1) i=m, c, a  (37)

where x_W is the mole fraction of water in the solution.Long-Range Interaction Contribution from Pitzer-Debye-Hückel (PDH) Model

To account for the long-range ion-ion interactions, the present invention eNRTL-SAC model uses the unsymmetric Pitzer-Debye-Hückel (PDH) formula (Pitzer, 1980) on the segment basis:

$\begin{array}{cc}\frac{G^{*\mathrm{ex},\mathrm{PDH}}}{n_S\mathrm{RT}}=-{\left(\frac{1000}{{\stackrel{\_}{M}}_S}\right)}^{\frac{1}{2}}\left(\frac{4A_{\varphi }I_x}{\rho }\right)\ln \left(1+\rho I_x^{\frac{1}{2}}\right)\mathrm{with}& \left(38\right)\\ A_{\varphi }=\frac{1}{3}{\left(\frac{2\pi N_A{\stackrel{\_}{d}}_S}{1000}\right)}^{\frac{1}{2}}{\left(\frac{Q_e^2}{{\stackrel{\_}{\varepsilon }}_Sk_BT}\right)}^{\frac{3}{2}}& \left(39\right)\\ I_x=\frac{1}{2}\sum _ix_iz_i^2& \left(40\right)\end{array}$

where A_φ is the Debye-Hückel parameter, I_x is the ionic strength (segment mole fraction scale), M_S is the average molecular weight of the mixed-solvents, ρ is the closest approach parameter, N_A is the Avogadro's number, d_S is the average density of the mixed-solvents, Q_e is the electron charge, ∈_S is the average dielectric constant of the mixed-solvents, k_B is the Boltzmann constant, and z_i(z_m=0;z_c=z_a=1) is the charge number of segment-based species i.

Applying the PDH model to the conceptual segments, the activity coefficient of segment species i can be derived as follows:

$\begin{array}{cc}\begin{array}{c}\ln {\Gamma }_i^{*\mathrm{PDH}}=\frac{1}{\mathrm{RT}}{\left(\frac{\partial G^{*\mathrm{ex},\mathrm{PDH}}}{\partial n_i}\right)}_{T,P,n_{j\ne i}}\\ =-{\left(\frac{1000}{{\stackrel{\_}{M}}_S}\right)}^{\frac{1}{2}}A_{\varphi }\left[\left(\frac{2z_i^2}{\rho }\right)\ln \left(1+\rho I_x^{\frac{1}{2}}\right)+\frac{z_i^2I_x^{\frac{1}{2}}-2I_x^{\frac{3}{2}}}{1+\rho I_x^{\frac{1}{2}}}\right]\end{array}i,j=m,c,a& \left(41\right)\end{array}$

The unsymmetric long range term for the logarithm of activity coefficient of component I is the sum of contributions from its various segments:

$\begin{array}{cc}\ln {\gamma }_I^{*\mathrm{PDH}}=\sum _mr_{m,I}\ln {\Gamma }_m^{*\mathrm{PDH}}+r_{c,I}\ln {\Gamma }_c^{*\mathrm{PDH}}+r_{a,I}\ln {\Gamma }_a^{*\mathrm{PDH}}\mathrm{where}& \left(42\right)\\ \ln {\Gamma }_m^{*\mathrm{PDH}}=2{\left(\frac{1000}{{\stackrel{\_}{M}}_S}\right)}^{\frac{1}{2}}\frac{A_{\varphi }I_x^{\frac{3}{2}}}{1+\rho I_x^{\frac{1}{2}}}& \left(43\right)\\ \ln {\Gamma }_c^{*\mathrm{PDH}}=\ln {\Gamma }_a^{*\mathrm{PDH}}=-{\left(\frac{1000}{{\stackrel{\_}{M}}_S}\right)}^{\frac{1}{2}}A_{\varphi }\left[\left(\frac{2}{\rho }\right)\ln \left(1+\rho I_x^{\frac{1}{2}}\right)+\frac{I_x^{\frac{1}{2}}-2I_x^{\frac{3}{2}}}{1+\rho I_x^{\frac{1}{2}}}\right]\mathrm{With}& \left(44\right)\\ A_{\varphi }=\frac{1}{3}{\left(\frac{2\pi N_A{\stackrel{\_}{d}}_S}{1000}\right)}^{\frac{1}{2}}{\left(\frac{Q_e^2}{{\stackrel{\_}{\varepsilon }}_Sk_BT}\right)}^{\frac{3}{2}}& \left(45\right)\\ I_x=\frac{1}{2}\left(x_c+x_a\right)& \left(46\right)\end{array}$

The Debye-Hückel theory is based on the infinite dilution reference state for ionic species in the actual solvent media. For systems with water as the only solvent, the reference state is the infinite dilution aqueous solution. For mixed-solvent systems, the reference state for which the Pitzer-Debye-Hückel formula remains valid is the infinite dilution solution with the corresponding mixed-solvent composition. Consequently, the molecular quantities for the single solvent need to be extended for mixed-solvents; simple composition average mixing rules are adequate to calculate them as follows:

$\begin{array}{cc}{\stackrel{\_}{M}}_S=\sum _Sx_S^{\prime }M_S& \left(47\right)\\ \frac{1}{{\stackrel{\_}{d}}_S}=\sum _S\frac{x_S^{\prime }}{d_S}& \left(48\right)\\ {\stackrel{\_}{\varepsilon }}_S=\sum _Sw_S^{\prime }{\varepsilon }_S\mathrm{with}& \left(49\right)\\ x_S^{\prime }=\frac{x_S}{\sum _Sx_S}& \left(50\right)\\ w_S^{\prime }=\frac{M_Sx_S}{\sum _SM_Sx_S}& \left(51\right)\end{array}$

where S is a solvent component in the mixture, and M_S is the molecular weight of the solvent S. It should be pointed out that Eqs. 47-51 should be used only in Eq. 41 and M_S, d_S, and ∈_S were already assumed as constants in Eqs. 38 and 39 when deriving Eq. 41 for mixed-solvent systems. Table 1 shows the values of dielectric constant at 298.15 K used in this study for the same sixty-two solvents investigated by Chen and Song (Chen, 2004a and U.S. application Ser. No. 10/785,925) above. These values were compiled from various sources including internet websites and commercial software Aspen Properties v2004.1 (by Aspen Technology, Inc. of Cambridge, Mass., assignee of the present invention).

Born Term Correction to Activity Coefficient

Given that the infinite dilution aqueous solution is chosen as the reference state, one needs to correct the change of the reference state from the mixed-solvent composition to aqueous solution for the Pitzer-Debye-Hüekel term. The Born term (Robinson, R. A. and R. H. Stokes, Electrolyte Solutions, 2^nd ed., Butterworths (1970), Rashin, A. A. and B. Honig, “Reevaluation of the Born Model of Ion Hydration, J. Phys. Chem., 89: 5588 (1985)) on the segment basis is used for this purpose:

$\begin{array}{cc}\frac{\Delta G^{\mathrm{Born}}}{n_S\mathrm{RT}}=\frac{Q_e^2}{2k_BT}\left(\frac{1}{{\stackrel{\_}{\varepsilon }}_S}-\frac{1}{{\varepsilon }_W}\right)\sum _i\frac{x_iz_i^2}{r_i}{10}^{-2}& \left(52\right)\end{array}$

ΔG^Born is the Born term correction to the unsymmetric Pitzer-Debye-Hückel formula, G*^ex,PDH, ∈_W is the dielectric constant of water, and r_i is the Born radius of segment specie i.

Applying Eq. 52 to all conceptual segments, the corresponding expression for the activity coefficient of segment species i can be derived as follows:

$\begin{array}{cc}\Delta \ln {\Gamma }_m^{\mathrm{Born}}=\frac{1}{\mathrm{RT}}{\left(\frac{\partial \Delta G^{\mathrm{Born}}}{\partial n_m}\right)}_{T,P,n_{j\ne m}}=0m=x,y-,y+,z& \left(53\right)\\ \begin{array}{c}\mathrm{\Delta ln}{\Gamma }_i^{\mathrm{Born}}=\frac{1}{\mathrm{RT}}{\left(\frac{\partial \Delta G^{\mathrm{Born}}}{\partial n_i}\right)}_{T,P,n_{j\ne i}}\\ =\frac{Q_e^2}{2k_BT}\left(\frac{1}{{\stackrel{\_}{\varepsilon }}_S}-\frac{1}{{\varepsilon }_W}\right)\frac{z_i^2}{r_i}{10}^{-2}\end{array}i=c,a& \left(54\right)\end{array}$

The Born correction term on the logarithm of activity coefficient of component I is the sum of contributions from its various segments:

$\begin{array}{cc}\Delta \ln {\gamma }_I^{\mathrm{Born}}=r_{c,I}\mathrm{\Delta ln}{\Gamma }_c^{\mathrm{Born}}+r_{a,I}\mathrm{\Delta ln}{\Gamma }_a^{\mathrm{Born}}& \left(55\right)\\ \mathrm{\Delta ln}{\Gamma }_c^{\mathrm{Born}}=\frac{Q_e^2}{2k_BT}\left(\frac{1}{{\stackrel{\_}{\varepsilon }}_S}-\frac{1}{{\varepsilon }_W}\right)\frac{1}{r_c}{10}^{-2}& \left(56\right)\\ \mathrm{\Delta ln}{\Gamma }_a^{\mathrm{Born}}=\frac{Q_e^2}{2k_BT}\left(\frac{1}{{\stackrel{\_}{\varepsilon }}_S}-\frac{1}{{\varepsilon }_W}\right)\frac{1}{r_a}{10}^{-2}& \left(57\right)\end{array}$

Flory-Huggins Term Correction to Activity Coefficient

Although in most common electrolyte systems, the combinatorial entropy of mixing term is much smaller than the residual term, one may still want to include it in a general model. Applicants follow the Polymer NRTL model (Chen 1993 above) and use the Flory-Huggins term to describe the combinatorial term:

$\begin{array}{cc}\frac{G^{\mathrm{ex},\mathrm{FH}}}{\mathrm{nRT}}=\sum _Ix_I\ln \left(\frac{{\phi }_I}{x_I}\right)\mathrm{with}& \left(58\right)\\ {\phi }_I=\frac{x_Ir_I}{\sum _Jx_Jr_J}& \left(59\right)\end{array}$

where G^ex,FH is the Flory-Huggins term for the excess Gibbs energy, φ_I is the segment fraction of component I, and r_I is the number of all conceptual segments in component I:

$\begin{array}{cc}{r}_I=\sum _mr_{m,I}+r_{c,I}+r_{a,I}& \left(60\right)\end{array}$

The activity coefficient of component I from the combinatorial term can be derived from Eq. 60:

$\begin{array}{cc}\ln {\gamma }_I^{\mathrm{FH}}=\ln \left(\frac{{\phi }_I}{x_I}\right)+1-r_I\sum _J\frac{{\phi }_J}{r_J}=\ln \left(\frac{r_I}{\sum _Jx_Jr_J}\right)+1-\frac{r_I}{\sum _Jx_Jr_J}& \left(61\right)\end{array}$

The infinite-dilution activity coefficient of a component in water is:

$\begin{array}{cc}\ln {\gamma }_I^{\infty \mathrm{FH}}=\ln \left(\frac{r_I}{r_W}\right)+1-\frac{r_I}{r_W}& \left(62\right)\end{array}$

In both NRTL-SAC (parent patent application) and present invention eNRTL-SAC, water is selected as the reference for the hydrophilic segment z. Therefore, one can set r_W=1. Thus, one has:

lnγ_I^∞FH=lnr_I+1−r_I  (63)

One can then compute the unsymmetric activity coefficient from the Flory-Huggins term as follows:

$\begin{array}{cc}\ln {\gamma }_I^{*,\mathrm{FH}}=\ln {\gamma }_I^{\mathrm{FH}}-\ln {\gamma }_I^{\infty \mathrm{FH}}=r_I-\ln \left(\sum _Jx_Jr_J\right)-\frac{r_I}{\sum _Jx_Jr_J}& \left(64\right)\end{array}$

NRTL Binary Parameters

In Eqs. 3 and 4 for NRTL-SAC, the model formulation requires the asymmetric interaction energy parameters, τ, and the symmetric nonrandom factor parameters, α, for each binary pair of the conceptual segments. In Eqs. 25-27 for eNRTL-SAC of the present invention, one needs additional binary parameters of τ and α between conceptual molecular segments, m and ionic segments, c or a. In practice, Applicants fix the values of α's for the binary pairs of molecular segment and ionic segment to the single value of 0.2 while the values of τ for the binary pairs of molecular segment and ionic segment are calculated from the τ′s for the binary pairs of molecular segment and electrolyte segment. Following the same scheme in generalized eNRTL (Chen and Song, 2004b above), one can calculate these binary interaction energy parameters as follows:

τ_cm=τ_am=τ_em  (65)

τ_mc,ac=τ_ma,ca=τ_me  (66)

Following the treatment of NRTL-SAC (disclosed in U.S. application Ser. No. 10/785,925), Applicants identify a reference electrolyte for the conceptual electrolyte segment e. In searching for the reference electrolyte, Applicants choose one elemental electrolyte that has abundant literature data. In one example study, NaCl is used as the reference electrolyte for e. The ionic radii for sodium ion and chloride ion are 1.680×10⁻¹⁰ m and 1.937×10⁻¹⁰ m, respectively. With NaCl as the reference electrolyte, the energy parameters for the z-e pair are set to (8.885, −4.549) for the water-NaCl pair. The energy parameters for the x-e pair are set to (15, 5), in line with the parameters identified for C₂H₄-NaCl pair earlier by Chen and Song (Chen, 2004b). The energy parameters for the y-e pairs are set to (12, −3) after limited trials to optimize the performance of the model in this study. The complete set of NRTL binary interaction energy parameters are given in Table 2. Other choices of the reference electrolyte and parameter values may be suitable. The below reports the general behavior of the present invention eNRTL-SAC model based on the parameters reported in Table 2.

The electrolyte segment e is the only extra molecular descriptor and the electrolyte parameter E is the only extra molecular parameter for all electrolytes, inorganic or organic. All local and long range interactions derived from the existence of cationic and anionic species of various ionic charge valence, radius, chemical make-up, etc., are to be accounted for with this extra molecular descriptor for electrolytes together with combinations of conceptual molecular segments, i.e., hydrophobicity, polarity and hydrophilicity. In other words, every electrolyte, organic or inorganic, are modeled as combinations of E, X, Y, and Z. As such, electrolytes are recognized as “hydrophobic” electrolytes, “polar” electrolytes, “hydrophilic” electrolytes, and their various combinations. Likewise, ionic activity coefficient of each ionic species will be computed from its share of E, X, Y, and Z. The ions are to be considered as “hydrophobic” ions, “polar” ions, or “hydrophilic” ions.

FIGS. 1 to 5 show effects of the molecular parameters on mean ionic activity coefficients (mole fraction scale) of the reference electrolyte, i.e., electrolyte with E=1. As shown in FIGS. 1 to 5, hydrophobicity parameter X brings down the mean ionic activity coefficient at low electrolyte concentration but in a rather nonlinear way. Polarity parameter Y− raises the mean ionic activity coefficient while polarity parameter Y+ lowers the mean ionic activity coefficient. Hydrophilicity parameter Z has a relatively slight downshift effect on the mean ionic activity coefficient. Electrolyte parameter E brings down the mean ionic activity coefficient at low electrolyte concentration and pushes up the mean ionic activity coefficient at high electrolyte concentration.

Experimental data for ionic activity coefficients are not readily available though emerging (Wilczek-Vera, G. et al, “On the Activity of Ions and the Junction Potential: Revised Values for All Data,” AIChE J., 50:445, 2004). Given the fact that existing experimental data are limited to mean ionic activity coefficient for neutral electrolytes, Applicants are not able to directly identify the molecular parameters for ionic species. In preparing FIGS. 1 to 5 discussed above and the subsequent studies reported in the Model Applications section below, Applicants use Eqs. 17-18 to determine from electrolyte parameter E the ionic segment numbers for the ions and Applicants arbitrarily assign molecular segment parameters (X, Y−, Y+, and Z) only to the anion. This practice is acceptable since virtually all electrolytes investigated in this study are electrolytes with elemental cations.

Reference is now made to a preferred embodiment of the present invention as illustrated in FIGS. 14-17. FIG. 14 illustrates a computer network or similar digital processing environment in which the present invention may be implemented.

Referring to FIG. 14, client computer(s)/devices 50 and server computer(s) 60 provide processing, storage, and input/output devices executing application programs and the like. Client computer(s)/devices 50 can also be linked through communications network 70 to other computing devices, including other client devices/processes 50 and server computer(s) 60. Communications network 70 can be part of a remote access network, a global network (e.g., the Internet), a worldwide collection of computers, Local area or Wide area networks, and gateways that currently use respective protocols (TCP/IP, Bluetooth, etc.) to communicate with one another. Other electronic device/computer network architectures are suitable.

FIG. 15 is a diagram of the internal structure of a computer (e.g., client processor/device 50 or server computers 60) in the computer system of FIG. 14. Each computer 50, 60 contains system bus 79, where a bus is a set of hardware lines used for data transfer among the components of a computer or processing system. Bus 79 is essentially a shared conduit that connects different elements of a computer system (e.g., processor, disk storage, memory, input/output ports, network ports, etc.) that enables the transfer of information between the elements. Attached to system bus 79 is I/O device interface 82 for connecting various input and output devices (e.g., keyboard, mouse, displays, printers, speakers, etc.) to the computer 50, 60. Network interface 86 allows the computer to connect to various other devices attached to a network (e.g., network 70 of FIG. 14). Memory 90 provides volatile storage for computer software instructions 92 and data 94 used to implement an embodiment of the present invention (e.g., eNRTL-SAC modeler 20 in FIGS. 16-17). Disk storage 95 provides non-volatile storage for computer software instructions 92 and data 94 used to implement an embodiment of the present invention. Central processor unit 84 is also attached to system bus 79 and provides for the execution of computer instructions.

In one embodiment, the processor routines 92 and data 94 are a computer program product (generally referenced 92 or 20), including a computer readable medium (e.g., a removable storage medium such as one or more DVD-ROM's, CD-ROM's, diskettes, tapes, etc.) that provides at least a portion of the software instructions for the invention system 20. Computer program product 92 can be installed by any suitable software installation procedure, as is well known in the art. In another embodiment, at least a portion of the software instructions may also be downloaded over a cable, communication and/or wireless connection. In other embodiments, the invention programs are a computer program propagated signal product 107 embodied on a propagated signal on a propagation medium (e.g., a radio wave, an infrared wave, a laser wave, a sound wave, or an electrical wave propagated over a global network such as the Internet, or other network(s)). Such carrier medium or signals provide at least a portion of the software instructions for the present invention routines/program 92.

In alternate embodiments, the propagated signal is an analog carrier wave or digital signal carried on the propagated medium. For example, the propagated signal may be a digitized signal propagated over a global network (e.g., the Internet), a telecommunications network, or other network. In one embodiment, the propagated signal is a signal that is transmitted over the propagation medium over a period of time, such as the instructions for a software application sent in packets over a network over a period of milliseconds, seconds, minutes, or longer. In another embodiment, the computer readable medium of computer program product 92 (e.g., eNRTL-SAC modeler 20) is a propagation medium that the computer system 50 may receive and read, such as by receiving the propagation medium and identifying a propagated signal embodied in the propagation medium, as described above for computer program propagated signal product.

Generally speaking, the term “carrier medium” or transient carrier encompasses the foregoing transient signals, propagated signals, propagated medium, storage medium and the like.

FIGS. 16 and 17 illustrate data flow and process steps for a modeler 20 performing the methods of the present invention. With reference to FIG. 16, chemical data describing one or more chemical species (e.g., an electrolyte and solvent) of the mixture and/or environmental conditions (e.g., pressure and/or temperature) is entered at step 105 of the modeler process. Step 110 uses that data and Equation 16 to determine at least one conceptual segment including a conceptual electrolyte segment for each of the chemical species of the mixture. The determined electrolyte conceptual segment and other determined conceptual segments are used to compute at least one physical property of the mixture during step 115. The computed physical properties are analyzed to form a model of at least one physical property of the mixture (e.g., solubility of one or more chemical species in one or more phases of the mixture) in step 120. The model information is then given as output at step 125. The output can take the form of data or an analysis appearing on a computer monitor, data or instructions sent to a process control system or device, data entered into a data storage device, and/or data or instructions relayed to additional computer systems or programs.

FIG. 17 illustrates in more detail the computation at step 115 in FIG. 16. Step 115 begins with the receipt of determined conceptual electrolyte and other segments for each of the chemical species (e.g., electrolyte, solvent, etc.) of the mixture. The determined conceptual electrolyte and other segments and Equation 16 are used to compute at least one physical property of the mixture during step 215. The computed physical properties are provided as output 220 from computation step 215. In step 220, the computed physical properties are passed to step 120 of FIG. 16 for forming a model of the physical property of the mixture as described above.

According to the foregoing, in some embodiments, the present invention 20 features a computer system for modeling at least one physical property of a mixture of at least two chemical species (e.g., an electrolyte dissolved in a liquid component/solvent. The computer system is formed of a user input means for determining chemical data from a user, a digital processor coupled to receive input from the input means, and an output means coupled to the digital processor. The digit processor hosts and executes a modeling system 20 in working memory. The modeling system 20 (i) uses the chemical data to determine one or more conceptual segments including a conceptual electrolyte segment for each of the chemical species; (ii) uses the determined conceptual segments to compute at least one physical property of the chemical mixture; and (iii) provides an analysis of the computed physical property. The analysis forms a model of the at least one physical property of the mixture. The output means provides to the user the formed model of the physical property of the chemical mixture. In some embodiments, the present invention features a chemical compound manufactured by a process that includes a modeling method 20. The modeling method 20 models at least one physical property of a mixture of at least one electrolyte and one or more solvents. The invention method 20 comprises the computer implemented steps of (i) determining at least a conceptual electrolyte segment for the electrolyte solute, (ii) using the determined conceptual electrolyte segment computing at least one physical property of the mixture; and (iii) providing an analysis of the computed physical property. The step of determining at least a conceptual electrolyte segment may further include determining other conceptual segments and defining an identity and an equivalent number of each conceptual segment. The provided analysis forms a model of at least one physical property of the mixture.

Model 20 Applications

The following Examples are illustrative of the invention, and are not meant to be limiting in any way.

Limited amount of mean ionic activity coefficient data are available in the public literature for aqueous electrolytes. Applicants test the present eNRTL-SAC model 20 against mean ionic activity coefficient data of aqueous electrolyte systems. In addition, Applicants test the model 20 against salt solubility data in multiple solvents for a number of inorganic electrolytes and organic electrolytes. To the best of Applicants' knowledge, public literature data is very scarce for such salt solubility data. Proprietary solubility data from industrial collaborators was also used to test the applicability of the present invention model 20. However, results with such proprietary solubility data are not included in this discussion.

Mean Ionic Activity Coefficients in Aqueous Systems

For an electrolyte CA that dissociates to cation C and anion A, the mean ionic activity coefficients γ*_± is related to individual ionic activity coefficients as follows:

$\begin{array}{cc}\ln {\gamma }_{\pm }^{*}=\frac{1}{v}\left(v_C\ln {\gamma }_C^{*}+v_A\ln {\gamma }_A^{*}\right)& \left(67\right)\end{array}$

where v=v_C+v_A.

Equation 67 gives the mean ionic activity coefficient on the mole fraction scale and it can be converted to the molality scale:

lnγ*_±m=lnγ*_±+ln(1+vmM_S/1000)  (68)

where γ*_±m is the mean ionic activity coefficient on the molality scale, m is the molality of the salt (mol/kg-solvent), and M_S is the molecular weight of the solvent (g/mol).

Table 3 shows the fit to molality scale mean ionic activity coefficient data and the identified electrolyte and molecular parameters for the aqueous inorganic and organic electrolytes at 298.15 K as compiled by of Robinson and Stokes (1970) cited above. All mean ionic activity coefficient data are assumed to have standard deviation of 5%. The data for C5 and higher sodium carboxylates were excluded from the fit because these organic electrolytes were known to form micelles at high electrolyte concentrations (Chen, C.-C. et al., “Segment-Based Excess Gibbs Energy Model for Aqueous Organic Electrolytes, AIChE J., 47:2593, 2001). With a few exceptions such as LiBr, most uni-univalent and uni-bivalent electrolytes are well represented as combinations of E and Y− or Y+ parameters. Most uni-univalent electrolytes have E parameter around unity while higher E values are found for higher valent electrolytes. Applicants also found that the fit seems to deteriorate for electrolytes with higher E values. This observation is consistent with the understanding that higher valent electrolytes are known to prone to the formation of hydrated species or other complexation species. The relatively poor representation of these electrolytes with the model reflects the inadequate assumption of complete dissociation for such electrolytes (Chen, C.-C.; et al., “Unification of Hydration and Dissociation Chemistries with the Electrolyte NRTL Model,” AIChE Journal, 45:1576, 1999). As a derived property, mean ionic activity coefficient becomes meaningless if the complete dissociation assumption of electrolytes does not hold true.

To illustrate the quality of the fit, FIG. 6 shows the comparison of experimental and calculated molality scale mean ionic activity coefficients for five aqueous electrolytes at 298.15 K. The solid lines are the calculated values from the model. It shows that the present invention eNRTL-SAC model 20 provides reasonable qualitative representation of the data while the original eNRTL model (Chen, C.-C. et al., “Local Composition Model for Excess Gibbs Energy of Electrolyte Systems,” AIChE J., 28:588, 1982) achieves excellent quantitative representation of the data.

Salt Solubility in Mixed Solvent Systems

At the solubility limit of a nonelectrolytes, the solubility product constant, K_sp, can be written in terms of the product of the solute concentration and the solute activity coefficient at the saturation concentration:

K_sp=x_Iγ_I  (69)

At the solubility limit of an electrolyte, ionic species precipitate to form salt.

v_CC^ZC+v_AA^ZA−>C_vCA_vA(s)  (70)

The corresponding solubility product constant can be defined as follows.

K_sp=x_C^v_Cγ*_C^vCx_A^vAγ*_A^vA  (71)

Eqs. 70 and 71 can be expanded to include solvent molecules and other species if the solid polymorph involves hydrates, other solvent-containing salts, double salts, triple salts, and others.

Applicants tested the applicability of present invention eNRTL-SAC model 20 with the very limited public literature data and some proprietary data on solubilities of a number of inorganic and organic electrolytes in various solvents. This description presents the results with solubility data from public literature. To bring certain consistency to the data treatment, Applicants convert all solute solubility data to mole fraction (except for sodium chloride and sodium acetate). Applicants also assign standard deviation of 10% to all solute solubility data within range of 1 to 0.1, standard deviation of 20% to all solute solubility data with range of 0.1 to 0.01, standard deviation of 30% to data with range of 0.01 to 0.001, and so on.

Solubility data of sodium chloride in twelve different solvents at 298.15 K were successfully fitted with the present invention model 20. (Note that the temperature for the acetone data is 291.15 K and the temperature for the ethyl acetate data is 292.15 K. However, they are included as if they were data at 298.15 K.) The sodium chloride solubilities in the twelve solvents vary by six orders of magnitude. The satisfactory fit of the data for ten solvents (formic acid and ethyl acetate excluded) is shown in FIG. 7. The present invention model 20 predicts one order-of-magnitude higher solubility for sodium chloride in formic acid and virtually no solubility for sodium chloride in ethyl acetate while the data suggests very low but measurable solubility. The molecular parameters and the solubility product constant were adjusted simultaneously to provide the best fit to the data and the identified values are given in Table 4. In Table 4, the last column to the right quantifies goodness of fit to data. It is worth noting that the electrolyte parameter E for sodium chloride is near unity, similar to the parameters reported in Table 3 for sodium chloride.

Solubility data of sodium acetate in five different solvents was also fitted successfully with the present invention model 20. The solubilities in the five solvents vary by four orders of magnitude. The fit of the data is shown in FIG. 8. The solid phase for the solubility measurements is anhydrous sodium acetate. Note that the data for methanol and acetone was taken at 291.15 K while the data for water and ethylene glycol was taken at 298.15 K. The temperature for the 1-propanol data is not known. In fitting the data, Applicants treated all data as if it was 298.15 K data. The identified molecular parameters and the solubility product constant are given in Table 4. As an organic electrolyte, the electrolyte parameter E for sodium acetate is found to be significantly less than unity.

FIGS. 9 a and 9b show satisfactory representations of the solubility data of benzoic acid in twenty-six solvents (Beerbower, A. et al., “Expanded Solubility Parameter Approach. I. Naphthalene and Benzoic Acid in Individual Solvents,” J. Pharm. Sci., 73:179, 1984) and the solubility data of sodium benzoate in ten solvents (Bustamante, P. et al., “The Modified Extended Hansen Method to Determine Partial Solubility Parameters of Drugs Containing a Single Hydrogen Bonding Group and Their Sodium Derivatives: Benzoic Acid/Na and Ibuprofen/Na,” Int. J. of Pharmaceutics, 194:117, 2000). These solvents are chosen in this study because of the availability of the NRTL-SAC parameters for the solvents from Applicants' prior work. The identified molecular parameters for the two solutes were given in Table 4. It is interesting that the molecular parameters identified for benzoic acid with twenty-six solvents in this study are quite similar to the molecular parameters identified for benzoic acid with seven solvents in Applicants' earlier study. Applicants also noted that the solubility range expands as benzoic acid is converted to sodium benzoate. Furthermore, the molecular parameters have changed from a hydrophobic/polar/hydrophilic combination (benzoic acid) to a polar/hydrophilic/electrolytic combination (sodium benzoate). Solubility data of sodium benzoate in seven other solvents (chloroform, benzene, dioxane, cyclohexane, ethyl acetate, heptane and chlorobenzene) is excluded from FIG. 9b because the present invention model 20 predicts virtually no solubility for sodium benzoate in these solvents while the data suggests very low but measurable solubility. It is probable that the molecular form of sodium benzoate may be present in such highly hydrophobic solvents. However, due to their low concentrations, Applicants chose to ignore these low solubility solvents in this study although the current thermodynamic framework can be used to account for the two solubility routes, i.e., Eqs. 69 and 71, individually or simultaneously.

FIGS. 10 a and 10b show successful representations of the solubility data of salicylic acid in eighteen solvents and the solubility data of sodium salicylate in thirteen solvents (Barra, J. et al., “Proposition of Group Molar Constants for Sodium to Calculate the Partial Solubility Parameters of Sodium Salts Using the van Krevelen Group Contribution Method,” Eur. J. of Pharm. Sci., 10:153, 2000). Their molecular parameters were given in Table 4. Like the molecular parameters for benzoic acid and the sodium salt, the molecular parameters have changed from a hydrophobic/polar/hydrophilic combination (salicylic acid) to a polar/hydrophilic/electrolytic combination (sodium salicylate). Solubility data of sodium salicylate in benzene, cyclohexane, and heptane is excluded from FIG. 10b, again because the present invention model 20 predicts virtually no solubility of sodium salicylate in these three solvents although the data suggests very low but measurable solubility. Acetic acid is the only outlier among solvents with significant solubility for sodium salicylate. The present invention model 20 prediction for the solubility of sodium salicylate in acetic acid is about one order of magnitude too high. Acetic acid is not included in the thirteen solvents shown in FIG. 10b.

The present invention model 20 results for the solubility data of p-aminobenzoic acid in nineteen solvents and sodium p-aminobenzoate in twelve solvents (Barra et al., 2000, above) are given in FIGS. 11a and 11b. Again, low solubility solvents (benzene, cyclohexane and heptane) are excluded from FIG. 11b for sodium aminobenzoate. Acetone and DMF are two outliers for sodium aminobenzoate and they are also excluded from FIG. 11b. The present invention model 20 predicts two orders of magnitude higher solubilities in these two solvents.

The solubility data and model calculations for ibuprofen in nineteen solvents and sodium ibuprofen in eleven solvents (Bustamante et al., 2000 above) are given in FIGS. 12a and 12b. In comparison to other organic solutes, one embodiment of the model 20 provides a rather poor fit to the ibuprofen data albeit a better fit than in prior nonelectrolyte models. Applicants did notice that the ibuprofen solubility data from Bustamante et al. are significantly different from those reported by Gracin and Rasmuson (Gracin, S. and A. C. Rasmuson, “Solubility of Phenylacetic Acid, p-Hydroxyphenylacetic Acid, p-Aminophenylacetic acid, p-Hydroxybenzoic acid, and Ibuprofen in Pure Solvents,” J. Chem. Eng. Data, 47:1379, 2002) for certain common solvents including methanol, ethanol, acetone and ethyl acetate. No attempt was made to reconcile the differences between the Bustamante data and the Gracin and Rasmuson data. The present invention model 20 fit to the sodium ibuprofen solubility data appears to be more satisfactory. Again, the eleven solvents reported in FIG. 12b do not include low solubility solvents (benzene, cyclohexane, heptane, and chlorobenzene). Similarly, acetone and DMF are two outliers for sodium ibuprofen and they are also excluded from FIG. 12b. The present invention model 20 predicts two orders of magnitude higher solubilities in these two solvents than the available data. Bustamante et al. (2000, above) reported high water content of the ibuprofen sample (3.3 wt % water) and the sodium ibuprofen sample (13 wt % water). It is not clear how such high water contents in the samples could impact on the solubility measurements.

The solubility data for diclofenac in sixteen solvents and sodium diclofenac in ten solvents (Barra et al., 2000 above) are fitted and reported in FIGS. 13a and 13b. The present invention model 20 significantly overestimates the solubilities of diclofenac in acetic acid, formamide and ethylene glycol. These three solvents are excluded from the sixteen solvents shown in FIG. 13a. Data for low solubility solvents (benzene, cyclohexane, ethyl acetate, heptane and chlorobenzene) for sodium diclofenac are excluded from FIG. 13b. Acetic acid and acetone are two outliers with the model estimations one to three orders of magnitude higher solubilities for sodium diclofenac. The two solvents are not included in FIG. 13b.

The solubility data treatment above assumes complete dissociation of electrolytes and considers the solubility problem as formation of salts from ionized species of electrolytes, i.e., Eq. 71. One may argue that electrolytes do not dissociate completely into ionic species especially in organic solvents of low dielectric constant. In the absence of dissociation to ionic species, the solubility relationship can be described by Eq. 69 and the eNRTL-SAC model of the present invention reduces to the NRTL-SAC model of the parent patent application. Applicants have treated the electrolyte systems above as nonelectrolytes (i.e., no dissociation to ionic species) with NRTL-SAC and the model results are also included in Table 4. With the absence of electrolyte parameter, the representation of the solubility data deteriorates substantially. Applicants also noted that the identified molecular parameters (X, Y−, Y+, and Z) with the complete dissociation treatment are roughly twice as large as those reported with the non-dissociation treatment. This finding is consistent with the fact that Applicants only assign the molecular parameters (X, Y−, Y+, and Z) to the anion.

TABLE 1
Dielectric Constant of Solvents at 298.15 K
solvent name                  dielectric constant at 298.15 K
Acetic acid                   6.13
Acetone                       20.83
Acetonitrile                  36.97
Anisole                       4.3
Benzene                       2.27
1-Butanol                     17.7
2-Butanol                     15.8
n-Butyl-acetate               5.1
Methyl-tert-butyl-ether       2.6
Carbon-tetrachloride          2.23
Chlorobenzene                 5.56
Chloroform                    4.7
Cumene                        2.22
Cyclohexane                   2.02
1,2-Dichloroethane            10.19
1,1-Dichloroethylene          4.6
1,2-Dichloroethylene          4.6
Dichloromethane               8.9
1,2-Dimethoxyethane           not available
N,N-Dimethylacetamide         not available
N,N-Dimethylformamide         38.3
Dimethyl-sulfoxide            47.2
1,4-Dioxane                   2.21
Ethanol                       24.11
2-Ethoxyethanol               not available
Ethyl-acetate                 6.02
Ethylene-glycol               41.2
Diethyl-ether                 4.26
Ethyl-formate                 7.16
Formamide                     109.5
Formic-acid                   58.5
n-Heptane                     1.92
n-Hexane                      1.89
Isobutyl-acetate              5.6
Isopropyl-acetate             not available
Methanol                      32.62
2-Methoxyethanol              not available
Methyl-acetate                6.68
3-Methyl-1-butanol            14.7
2-Hexanone                    14.6
Methylcyclohexane             2.02
Methyl-ethyl-ketone           18.5
Methyl-isobutyl-ketone        13.1
Isobutanol                    17.9
N-Methyl-2-pyrrolidone        33
Nitromethane                  6.26
n-Pentane                     1.84
1-Pentanol                    13.9
1-Propanol                    20.1
Isopropyl-alcohol             19.9
n-propyl-acetate              6
Pyridine                      2.3
Sulfolane                     43.3
Tetrahydrofuran               7.52
1,2,3,4-Tetrahydronaphthalene not available
Toluene                       2.36
1,1,1-Trichloroethane         7.5
Trichloroethylene             3.42
m-Xylene                      2.24
Water                         78.54
Triethylamine                 2.44
1-Octanol                     10.3

TABLE 2
NRTL Binary Interaction Parameters
	Segment (1)
	x	x	y−	y+	x
	Segment (2)
		y−	z	z	z	y+
	τ12	1.643	6.547	−2.000	2.000	1.643
	τ21	1.834	10.949	1.787	1.787	1.834
	α12 = α21	0.2	0.2	0.3	0.3	0.2
	Segment (1)
	x	y−	y+	z
	Segment (2)
		e	e	e	e
	τ12	15	12	12	8.885
	τ21	5	−3	−3	−4.549
	α12 = α21	0.2	0.2	0.2	0.2

TABLE 3
Results of Fit for Molality Scale Mean Ionic Activity
Coefficient Data of Aqueous Electrolytes
at 298.15K (Data from Robinson and Stokes, 1970)
					max.
	E	Y−	Y+	σ□	molality
1-1 Electrolytes
AgNO3	0.738		1.758	0.050	6.0
CsAc	1.002	0.438		0.011	3.5
CsBr	0.950		0.678	0.013	5.0
CsCl	0.948		0.643	0.014	6.0
CsI	0.956		0.719	0.012	3.0
CsNO3	0.981		1.328	0.005	1.4
CsOH	0.942	0.354		0.002	1.0
HBr	1.135	0.654		0.034	3.0
HCl	1.324	0.524		0.087	6.0
HClO4	1.476	0.569		0.136	6.0
HI	1.117	0.824		0.035	3.0
HNO3	0.971	0.211		0.005	3.0
KAc	0.998	0.386		0.009	3.5
KBr	0.910		0.311	0.011	5.5
KBrO3	0.968		1.141	0.002	0.5
KCl	0.920		0.370	0.010	4.5
KClO3	0.958		1.053	0.003	0.7
KCNS	0.876		0.477	0.019	5.0
KF	0.987		0.042	0.004	4.0
KH Malonate	0.846		0.920	0.022	5.0
KH Succinate	0.912		0.665	0.011	4.5
KH2PO4	0.970		1.362	0.006	1.8
KI	0.903		0.168	0.011	4.5
KNO3	0.856		1.461	0.027	3.5
KOH	1.236	0.344		0.058	6.0
K Tol	0.750		1.296	0.026	3.5
LiAc	0.962	0.097		0.002	4.0
LiBr	1.422	0.526		0.116	6.0
LiCl	1.282	0.436		0.084	6.0
LiClO4	1.145	0.681		0.047	4.0
LiI	1.058	0.712		0.033	3.0
LiNO3	1.050	0.294		0.022	6.0
LiOH	1.028		0.652	0.022	4.0
LiTol	0.881		0.392	0.014	4.5
NaAc	0.978	0.301		0.005	3.5
NaBr	0.992	0.115		0.008	4.0
NaBrO3	0.923		0.802	0.010	2.5
Na Butyrate	0.989	0.566		0.009	3.5
NaCl	1.000			0.017	6.0
NaClO3	0.891		0.507	0.011	3.5
NaClO4	0.894		0.267	0.010	6.0
NaCNS	0.925	0.128		0.006	4.0
NaF	0.976		0.425	0.002	1.0
Na Formate	0.905		0.094	0.013	3.5
NaH Malonate	0.878		0.664	0.019	5.0
NaH Succinate	0.924		0.495	0.010	5.0
NaH2PO4	0.864		1.256	0.020	6.0
NaI	1.009	0.266		0.012	3.5
NaNO3	0.825		0.842	0.029	6.0
NaOH	1.080	0.109		0.039	6.0
Na Propionate	0.992	0.448		0.006	3.0
Na Tol	0.793		0.920	0.026	4.0
NH4Cl	0.884		0.424	0.019	6.0
NH4NO3	0.813		1.128	0.043	6.0
RbAc	1.012	0.416		0.011	3.5
RbBr	0.914		0.519	0.016	5.0
RbCl	0.929		0.466	0.012	5.0
RbI	0.925		0.520	0.014	5.0
RbNO3	0.815		1.611	0.038	4.5
TlAc	0.864		0.952	0.033	6.0
TlClO4	1.020		1.231	0.000	0.5
TlNO3	1.069		1.692	0.003	0.4
1-2 Electrolytes
Cs2SO4	1.161		2.568	0.050	1.8
K2CrO4	1.048		2.738	0.075	3.5
K2SO4	1.386		2.475	0.021	0.7
Li2SO4	1.138		2.177	0.051	3.0
Na2CrO4	1.091		2.443	0.051	4.0
Na2 Fumarate	1.259		1.770	0.041	2.0
Na2 Maleate	1.202		2.699	0.075	3.0
Na2SO4	0.988		3.273	0.090	4.0
Na2S2O3	1.071		2.709	0.064	3.5
(NH4)2SO4	1.006		3.477	0.118	4.0
Rb2SO4	1.150		2.743	0.052	1.8
1-3 Electrolytes
K3Fe(CN)6	1.328		4.996	0.101	1.4
1-4 Electrolytes
K4Fe(CN)6	1.449		9.448	0.146	0.9
2-1 Electrolytes
BaAc2	1.016		0.997	0.128	3.5
BaBr2	1.267		0.358	0.018	2.0
BaCl2	1.227		0.585	0.029	1.8
Ba(ClO4)2	1.305		0.261	0.049	5.0
BaI2	1.354	0.028		0.017	2.0
Ba(NO3)2	1.435		1.268	0.008	0.4
CaBr2	1.969		0.171	0.495	6.0
CaCl2	1.701		0.309	0.283	6.0
Ca(ClO4)2	2.021			0.431	6.0
CaI2	1.419	0.131		0.036	2.0
Ca(NO3)2	1.108		0.875	0.053	6.0
CdBr2	1.324		3.164	0.294	4.0
CdCl2	1.052		3.047	0.315	6.0
CdI2	1.780		3.820	0.337	2.5
Cd(NO3)2	1.176		0.500	0.037	2.5
CoBr2	1.779			0.218	5.0
CoCl2	1.397		0.194	0.046	4.0
CoI2	2.260			0.488	6.0
Co(NO3)2	1.444		0.296	0.113	5.0
CuCl2	1.033	0.425	1.217	0.069	6.0
Cu(NO3)2	1.409		0.416	0.117	6.0
FeCl2	1.319		0.255	0.011	2.0
MgAc2	1.192		0.946	0.059	4.0
MgBr2	1.941			0.347	5.0
MgCl2	1.745		0.144	0.275	5.0
Mg(ClO4)2	1.988	0.162		0.303	4.0
MgI2	2.237			0.470	5.0
Mg(NO3)2	1.493		0.198	0.140	5.0
MnCl2	1.273		0.343	0.020	6.0
NiCl2	1.533		0.189	0.123	5.0
Pb(ClO4)2	1.549		0.236	0.184	6.0
Pb(NO3)2	1.129		1.964	0.083	2.0
SrBr2	1.330		0.183	0.023	2.0
SrCl2	1.401		0.357	0.082	4.0
Sr(ClO4)2	1.742		0.034	0.261	6.0
SrI2	1.384	0.076		0.030	2.0
Sr(NO3)2	0.978		1.250	0.091	4.0
UO2Cl2	1.277		0.024	0.017	3.0
UO2(ClO4)2	2.854			0.883	5.5
UO2(NO3)2	1.392	0.372	0.490	0.036	5.5
ZnBr2	0.906		0.337	0.088	6.0
ZnCl2	0.953		0.971	0.065	6.0
Zn(ClO4)2	2.045	0.130		0.318	4.0
ZnI2	0.868	0.132		0.116	6.0
Zn(NO3)2	1.518		0.214	0.176	6.0
2-2 Electrolytes
BeSO4	1.376		4.077	0.233	4.0
MgSO4	1.380		4.206	0.238	3.0
MnSO4	1.287		4.460	0.271	4.0
NiSO4	1.398		4.381	0.220	2.5
CuSO4	1.587		4.114	0.154	1.4
ZnSO4	1.339		4.417	0.242	3.5
CdSO4	1.295		4.547	0.271	3.5
UO2SO4	1.215		4.528	0.309	6.0
3-1 Electrolytes
AlCl3	1.730		0.579	0.087	1.8
CeCl3	1.562		0.883	0.047	1.8
CrCl3	1.589		0.641	0.022	1.2
Cr(NO3)3	1.551		0.761	0.036	1.4
EuCl3	1.586		0.820	0.049	2.0
LaCl3	1.553		0.877	0.042	2.0
NdCl3	1.575		0.882	0.045	2.0
PrCl3	1.562		0.892	0.042	2.0
ScCl3	1.636		0.709	0.041	1.8
SmCl3	1.581		0.843	0.046	2.0
YCl3	1.629		0.807	0.057	2.0
3-2 Electrolytes
Al2(SO4)3	1.354		4.886	0.222	1.0
Cr2(SO4)3	1.257		4.549	0.218	1.2
4-1 Electrolytes
Th(NO3)4	1.273		1.251	0.056	5.0
1σ is defined to be
${\left[\sum _i^N{\left(\frac{{\gamma }_{\pm i}^{*\exp }-{\gamma }_{\pm i}^{*\mathrm{cal}}}{{\gamma }_{\pm i}^{*\exp }}\right)}^2/N\right]}^{1/2}$where γ±* is the mean ionic acitivity acoefficient of electrolyte and N is the number of data used in correlations

TABLE 4
eNRTL-SAC Model Parameters for Solutes
solute	no. of solvents	X	Y−	Y+	Z	E	ln Ksp	σ4
benzoic acid1	26	0.494		0.336	0.468		−1.714	0.292
salicylic acid1	18	0.726	0.176		0.749		−1.624	0.774
p-aminobenzoic acid1	19	0.552	0.423	0.594	0.881		−3.348	1.206
Ibuprofen1	19	1.038	0.051	0.028	0.318		−1.423	1.055
Diclofenac1	16	0.158		1.678	0.451		−3.560	0.991
sodium chloride2	10			1.444		0.994	−6.252	0.783
sodium acetate2	5			1.417		0.521	−6.355	0.241
sodium benzoate2	10		0.750	1.685	2.201	0.539	−7.312	0.493
sodium salicylate2	13			0.845	2.417	0.090	−4.889	0.771
sodium p-	12			2.299	2.387	0.192	−8.293	1.258
aminobenzoate2
sodium ibuprofen2	11	1.819	1.743		2.362	0.150	−17.844	0.886
sodium diclofenac2	10	0.409		3.558	3.486	0.161	−14.202	0.858
sodium chloride3	10			1.060	2.200		−3.540	0.923
sodium acetate3	5			0.249	0.679		−2.277	0.281
sodium benzoate3	10			0.179	1.825		−2.978	0.699
sodium salicylate3	13			0.373	1.572		−2.153	1.058
sodium p-	12		0.125	0.649	1.895		−3.247	1.904
aminobenzoate3
sodium ibuprofen3	11	0.270	0.394		0.823		−2.364	1.685
sodium diclofenac3	10	0.454		0.124	2.493		−4.405	1.473
1nonelectrolytes
2electrolytes
3treated as nonelectrolytes
4σ is defined to be
${\left(\sum _i^N{\left(\ln x_i^{\exp }-\ln x_i^{\mathrm{cal}}\right)}^2/N\right)}^{1/2}$where x is the solubility of solute, i.e., mole fraction (note that mass fraction for sodium chloride and sodium acetate), and N is the number of data used in correlations.

CONCLUSIONS

The NRTL-SAC model of related U.S. patent application Ser. No. 10/785,925 (by assignee), a practical thermodynamic framework for solubility modeling of organic nonelectrolytes, has been extended for electrolytes. The electrolyte NRTL-SAC model 20 of the present invention requires only one additional component-specific electrolyte parameter over the three types of molecular parameters associated with the NRTL-SAC model. For solute molecules, these parameters are identified from solubility measurements of the solute in a few representative solvents, i.e., hydrophobic, hydrophilic and polar solvents. While scarcity of public literature solubility data on organic electrolytes has hampered extensive testing, Applicants have shown the extended model to be a promising tool for qualitative correlation and estimation of solubility of electrolyte systems including systems with large, complex organic electrolytes.

While this invention has been particularly shown and described with references to preferred embodiments thereof, it will be understood by those skilled in the art that various changes in form and details may be made therein without departing from the scope of the invention encompassed by the appended claims.

Claims

1. A chemical compound manufactured by a process that includes a modeling method, wherein the modeling method models at least one physical property of a mixture of at least one electrolyte dissolved in at least one solvent and comprises the computer implemented steps of: a) determining a conceptual electrolyte segment for the electrolyte, including defining a segment number;b) using the determined conceptual electrolyte segment, computing at least one physical property of the mixture; andc) providing an analysis of the computed physical property, wherein the analysis forms a model of the at least one physical property of the mixture.

2. The chemical compound in claim 1, wherein the chemical compound is any of a pharmaceutical compound, a nonpolymeric compound, a polymer, an oligomer, an inorganic compound and an organic compound.

3. A chemical compound comprising at least one electrolyte selected by a model, wherein the model determines a conceptual electrolyte segment of the electrolyte and computes at least one physical property of the electrolyte.

4. The chemical compound of claim 3, wherein the chemical compound is a pharmaceutical compound, a nonpolymeric compound, a polymer, an oligomer, an inorganic compound or an organic compound.

5. The chemical compound of claim 3, wherein the conceptual electrolyte segment includes a segment number, and the segment number includes a cationic segment number and an anionic segment number.

6. The chemical compound of claim 5, wherein the cationic segment number and the anionic segment number are represented by the following formulas: rc,C=re,CAZC and ra,A=re,CAZA,

7. The chemical compound of claim 5, wherein the cationic segment number and the anionic segment number are represented by the following formulas:

8. The chemical compound of claim 3, wherein the computed physical property includes at least one of activity coefficient, vapor pressure, solubility, boiling point, freezing point, octanol/water partition coefficient, and lipophilicity of the electrolyte.

9. The chemical compound of claim 8, wherein solubility is represented by the following formula:

10. The chemical compound of claim 3, wherein the electrolyte is symmetrical or unsymmetrical.

11. The chemical compound of claim 3, wherein the electrolyte is univalent or multivalent.

12. The chemical compound of claim 3, wherein the electrolyte includes two or more ionic species.

13. The chemical compound of claim 3, wherein the step of computing at least one physical property includes calculating the activity coefficient of the ionic species derived from the at least one electrolyte.

14. The chemical compound of claim 3, wherein the solvent is water, and the step of computing at least one physical property includes calculating: lnγ*I=lnγ*Ilc+lnγ*IPDH+lnγ*IFH,

15. The chemical compound of claim 3, wherein the one or more solvents include mixed-solvent solutions, and the step of computing at least one physical property including calculating: lnγ*I=lnγ*Ilc+lnγ*IPDH+lnγ*IFH+ΔlnγIBorn,

16. The chemical compound of claim 1, wherein the computed physical property includes at least one of activity coefficient, vapor pressure, solubility, boiling point, freezing point, octanol/water partition coefficient, and lipophilicity of the electrolyte.

17. The chemical compound of claim 1, wherein the electrolyte is symmetrical or unsymmetrical.

18. The chemical compound of claim 1, wherein the electrolyte is univalent or multivalent.

19. The chemical compound of claim 1, wherein the electrolyte includes two or more ionic species.

20. The chemical compound of claim 1, wherein the step of computing at least one physical property includes calculating the activity coefficient of the ionic species derived from the at least one electrolyte.