US20130270492A1 - Method for predicting the conductivity of a liquid mixture - Google Patents
Method for predicting the conductivity of a liquid mixture Download PDFInfo
- Publication number
- US20130270492A1 US20130270492A1 US13/994,916 US201113994916A US2013270492A1 US 20130270492 A1 US20130270492 A1 US 20130270492A1 US 201113994916 A US201113994916 A US 201113994916A US 2013270492 A1 US2013270492 A1 US 2013270492A1
- Authority
- US
- United States
- Prior art keywords
- conductivity
- species
- ionic species
- liquid solution
- values
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Abandoned
Links
Images
Classifications
-
- G06F19/70—
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/30—Prediction of properties of chemical compounds, compositions or mixtures
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
- B01D—SEPARATION
- B01D15/00—Separating processes involving the treatment of liquids with solid sorbents; Apparatus therefor
- B01D15/08—Selective adsorption, e.g. chromatography
- B01D15/10—Selective adsorption, e.g. chromatography characterised by constructional or operational features
- B01D15/12—Selective adsorption, e.g. chromatography characterised by constructional or operational features relating to the preparation of the feed
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
- B01F—MIXING, e.g. DISSOLVING, EMULSIFYING OR DISPERSING
- B01F35/00—Accessories for mixers; Auxiliary operations or auxiliary devices; Parts or details of general application
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
- G01N27/02—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance
- G01N27/04—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating resistance
- G01N27/06—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating resistance of a liquid
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
- G01N27/26—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating electrochemical variables; by using electrolysis or electrophoresis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
- G01N27/26—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating electrochemical variables; by using electrolysis or electrophoresis
- G01N27/416—Systems
- G01N27/4166—Systems measuring a particular property of an electrolyte
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C99/00—Subject matter not provided for in other groups of this subclass
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
- G01N30/26—Conditioning of the fluid carrier; Flow patterns
- G01N30/28—Control of physical parameters of the fluid carrier
- G01N30/34—Control of physical parameters of the fluid carrier of fluid composition, e.g. gradient
Definitions
- the present invention relates to the preparation of liquid mixtures, and more particularly to the preparation of a liquid mixture, such as a buffer, wherein the conductivity of the liquid mixture is predicted.
- a buffer solution is typically an aqueous solution consisting of a mixture of a weak acid and its conjugate base, or a weak base and its conjugate acid, and has the property that the pH of the solution changes very little when a small amount of strong acid or base is added to it. Buffer solutions are therefore used as a means of keeping pH at a nearly constant value in a wide variety of chemical applications, including e.g. chromatography, filtration, etc.
- a buffer solution may be made up of more than one weak acid and its conjugate base. For instance, a wider buffer region may be created by mixing two buffering agents with overlapping individual buffer regions.
- a buffer solution having a desired pH, and optionally also ionic strength may be prepared by calculating the necessary amounts of the ingredients of the buffer and mixing them. While it is often necessary to solve several different equations for calculating the pH of a relatively simple mixture of a weak acid (or base) and a strong base (or acid) depending on their relative concentrations, there is commercial software available for performing such calculations. For many applications, it would be desirable to also be able to predict the expected conductivity of the buffer to be prepared. However, due to the complex nature of buffers solutions, being combinations of strong electrolytes like salts and strong acids or bases that dissociate completely and weak electrolytes like buffer substances that ionize partially, this has so far generally not been possible.
- the prediction of the conductivity comprises solving the exact concentrations of the different ions present in solution including the equilibrium concentration of each of the charged species of a weak electrolyte, determining the molar conductivity of each of the charged species, calculating the corresponding conductivities, and summing up all the individual conductivities to obtain the total conductivity of the liquid mixture.
- calculating the molar conductivities of the subspecies of the weak electrolyte includes applying predetermined values for the Kohlrausch coefficient used when calculating the molar conductivities.
- the predetermined Kohlrausch coefficient values are obtained by fitting measured conductivities to the equations used for determining molar conductivities.
- the present invention therefore provides a method of preparing a liquid solution by mixing ingredients according to a predetermined recipe, wherein the liquid solution comprises a plurality of species, at least one pair of species being derived from a weak electrolyte and corresponding to an acid-base pair, and wherein the conductivity of the liquid solution is predicted by:
- ⁇ is the molar conductivity
- ⁇ 0 is the molar conductivity at infinite dilution
- c is the concentration of the ionic species
- K is the Kohlrausch coefficient
- Sqrt is square root
- K and ⁇ 0 are obtained from a data set comprising predetermined values for K and ⁇ 0 for each ionic species
- the data set comprises values for the Kohlrausch coefficient K obtained by measuring the conductivity of a solution containing the ionic species at a plurality of different concentrations of the ionic species and different pH values, and optionally different temperatures, and fitting the resulting data to the formula in step (ii) above to obtain K.
- Values for ⁇ 0 in the data set may include previously known values, such as literature values, and/or values be obtained by the fitting.
- the concentration of each ionic species is preferably calculated by an algorithm comprising the equation of Debye-Hückel, wherein the ionic strength of each species is used as a weighting parameter in the calculation of an average hydrated radius of the ionic species, as will be described in more detail below.
- the above-described method may advantageously be used for controlling a buffer formulation system or an in-line dilution system.
- the method may also be used in screening experiments wherein conductivity is used as a design of experiment (DoE) parameter.
- DoE design of experiment
- the present invention provides a computer program product comprising instructions for performing the steps of the method aspect above.
- the conductivity prediction steps in the method of the invention may further be used in a “backwards” fashion to indirectly determine pH of a solution by measuring its conductivity.
- the present invention therefore provides a device for measuring pH, comprising a conductivity sensor and means for calculating pH from measured conductivity using the conductivity prediction steps of the method aspect above in a backwards calculation mode.
- Yet another aspect of the invention provides a method for predicting conductivity of a liquid solution comprising steps (i) to (iv) of the method aspect above.
- FIG. 1 is a graph showing measured conductivity versus predicted conductivity for a buffer solution.
- a buffer solution is an aqueous solution typically consisting of a mixture of a weak acid and its conjugate base or a weak base and its conjugate acid. It has the property that the pH of the solution changes very little when a small amount of strong acid or base is added to it. Buffer solutions are used as a means of keeping pH at a nearly constant value in a wide variety of chemical applications
- a titrand is the solution (or other substance) into which a titrant is added.
- the titrand may, for example, be a weak acid or base, and the titrant a strong base or acid, respectively.
- a “strong” electrolyte is a solute hat completely, or almost completely, ionizes or dissociates in a solution. These ions are good conductors of electric current in the solution.
- strong electrolytes include, for instance, salts, such as sodium chloride; strong acids, such as hydrochloric acid, sulfuric acid, nitric acid; and strong bases, such as potassium hydroxide, sodium hydroxide, calcium hydroxide.
- a weak electrolyte dissociates into ions only to a small extent, typically much less than 10%, and the vast bulk of the weak electrolyte remains in its original unionized form in solution.
- most organic acids and their salts, and the salts of organic bases are weak electrolytes.
- weak electrolytes include, for example acetic acid and citric acid and salts thereof.
- the amount of weak electrolyte that dissociates depends on the number of counter-ions present in the solution.
- the conductivity (or specific conductance) of an electrolyte solution is a measure of its ability to conduct electricity.
- the SI unit of conductivity is “siemens” per meter (S/m).
- the ionic strength of a solution is a function of the concentration of all ions in the solution (half the sum of concentration multiplied by the square of ionic charge for all ions).
- the ionic strength is typically given in the unit mol/dm 3 .
- the present invention relates to the prediction of the conductivity of a buffer solution prepared according to a desired recipe or formula.
- a buffer having a desired pH To prepare a buffer having a desired pH, the added amounts of moles of titrand, titrant, and salt solution, when required, can be calculated. From these mole values, the volume needed from corresponding stock solutions can in turn be calculated. Typically, such calculations are performed by suitable computer software, which is commercially available. The use of such software permits the preparation of buffers up to concentrations of 1 M or higher, as will be described in more detail below.
- Atkins, P., and de Paula, J., 2002 describes calculation of molar conductivities for weak electrolytes for the ideal case where activity coefficients are ignored, which assumption leads in most cases to gross errors in the solution of the equilibrium equations.
- ⁇ is the measured conductivity
- c is the electrolyte concentration
- ⁇ 0 is the molar conductivity at infinite dilution, or “limiting molar conductivity”
- K is the “Kohlrausch coefficient”, which depends on the nature of the specific salt in solution.
- the constant K can be expressed as A+B ⁇ 0 (Debye, P., Hückel, H., Phys. Z. 1924, 25, 49) or A+B ⁇ w ⁇ 0 (Onsager, L., Phys. Z, 1926, 27, 388-392, and 1927, 28, 277-298), where A and B are constants depending on the temperature (T), dielectric constant ( ⁇ ) and charges of the ions (A also depends on the viscosity of the liquid), and w is a factor introduced by Onsager which contains the contribution from self-Brownian motion of the ion (Onsager, L., 1926, 1927 above).
- the molar conductivity A decreases as a function of the square root of the concentration in a linear fashion. This is an approximation that is valid in the low concentration regime. As the concentration increases, the rate of decrease of the molar conductivity A decreases so that it levels off, as can for instance be seen in FIG. 1 of the above cited reference Chandra, A., 1999. Besides, this reference also presents a microscopic theory that accounts for the constants A and B by treating the solvent as a dielectric continuum and calculating the concentration dependent ion contribution to microscopic friction.
- the limiting molar conductivity ⁇ 0 can be decomposed into contributions from the different ions (commonly known as law of independent migration of ions):
- ⁇ i is the molar ionic conductivities of ion i
- ⁇ i is the number of ions i in the formula unit of the electrolyte (e.g. 2 and 1 for Na + and SO 4 2 ⁇ , respectively, in Na 2 SO 4 ).
- the conductivity can then be calculated with Equations (2) and (3) above once the constants K and ⁇ 0 , or alternatively the constants A, B and ⁇ 0 (and optionally w), have been determined.
- This determination can be made experimentally by measuring the conductivities of a large number of samples of the buffer at different concentrations, temperatures, and pH values, and then fitting the obtained data to Equation (2).
- this includes solving the equilibrium equations for the buffering species.
- ⁇ becomes 1 and the activity of every ion will be equal to the corresponding concentration.
- the ionic strength is different from zero and the activity coefficients of the different species will be less than 1.
- ⁇ a and ⁇ b are the activity coefficients for the acidic and basic species, respectively, and pK a ′ is an apparent pK a value which allows the use of the measureable concentrations of the different buffer species.
- A is a constant, or rather a temperature dependent parameter ⁇ 0.51.
- pK a ′ pK a +( AZ a 2I 1 0.5 )/(1+0.33*10 8 a a I 0.5 ) ⁇ ( AZ b 2 I 0.5 )/(1+0.33*10 8 a b I 0.5 ) (Eq. 9)
- Z a charge of acidic species
- Z b charge of basic species
- a a ion size parameter of the acidic species
- a b ion size parameter of the basic species
- pK a ′ pK a +( AZ a 2 I 0.5 )/(1+ I 0.5 ) ⁇ ( AZ b 2 I 0.5 )/(1+ +I 0.5 ) (Eq. 10)
- buffers are polyprotic, i.e. their buffer molecules can accept and give away more than one proton corresponding to more than one pK a value.
- the number of species in such a buffer system is always one more than the number of pK a values.
- the calculation of the amount of moles in each of the protonation species is equivalent to solving the equilibrium equations of each of the species with the “neighbouring” species with one more and/or one less proton and with the concentration of hydrogen atoms (the pH).
- start_species is meant the species of the buffer substance before mixing, i.e. the protonation state of the buffer in the can or stock solution. This protonation state is determined by the amount of counter-ions per buffer molecule because the macroscopic object, can or stock solution, has to be electrically neutral.
- the minus sign in front of [OH ⁇ ] is due to the minus sign of the charge of the OH ions, whereas the minus sign in front of the two last terms is due to the charge of the counter-ions of the titrant and the start_species, respectively.
- Equations (11) to (14) above imply that there are six equations with six unknowns (the four [s i ], [OH ⁇ ] and [H + ]) for the case of three pK a values, and the equilibrium concentrations of the respective acidic and basic species may thus be calculated.
- s 1 corresponds to the acidic species
- s 2 to the basic species
- concentrations of s 3 and s 4 are set to zero.
- a is determined as the weighted mean ion size of all species contributing to the ionic strength of the liquid mixture, using the ionic strength as weighting parameter.
- the ion size parameter a may then be calculated as
- I i is the ionic strength
- a i is the ion size parameter of species i
- I is the total ionic strength defined by
- the ion size parameter a may also be approximated as
- shell is fixed at a specific value for a positively charged ionic species (typically in the range of 3.9-4.1), and at another fixed value for a negatively charged species (typically in the range of 0-0.1, or 0).
- WO 2009/131524 A1 there is also disclosed a computer program, or software, for calculating a recipe (i.e. amounts of titrand, titrant, water and salt) to achieve a solution of a given pH and ionic strength (or, alternatively, salt concentration instead of ionic strength).
- a computer program “Buffalo Plus” is disclosed which is written in Visual C++ and can be run in Windows for accurate calculation of buffer pH at ionic strengths as high as 1 M for buffer systems including phosphate, citrate, acetate and Tris.
- Such software may conveniently be used for calculating the equilibrium concentrations of acidic and basic species for the purposes of the present invention.
- the conductivity of a solution containing the weak electrolyte is measured at a plurality of different concentrations and different pH values, and optionally different temperatures.
- the resulting conductivities are then used together with equilibrium concentrations obtained by solving equilibrium equations as described above to find through a computational procedure of regression analysis type the values for K that give the best fit to a regression model provided by the above-mentioned equations (2) and (3), i.e.
- measurements are made for a large number a different weak electrolyte ionic species at different concentrations and pH to produce a catalogue or database of values for the coefficient K, optionally together with ⁇ 0 / ⁇ 0 values.
- the method of the invention may be implemented by software run on an electrical data processing device, such as a computer.
- Such software may be provided to the computer on any suitable computer-readable medium, including a record medium, a read-only memory, or an electrical or optical signal which may be conveyed via electrical or optical cable or by radio or other means.
- Predicting conductivities as described above may be used for several purposes. Exemplary applications include the use for controlling a buffer formulation system or an in-line dilution system, wherein a buffer solution of predetermined conductivity, pH (and potentially other parameters) is prepared by mixing relative component proportions of at least one each of:
- Such prediction of conductivities may also be used in screening experiments wherein conductivity is used as a design of experiment (DoE) parameter. Still another application is for determining the exact concentration of a stock solution of a weak electrolyte. Yet another application is for indirect measurement of pH by measuring conductivity and calculating pH therefrom.
- DoE design of experiment
- conductivities for a plurality of compositions of a formate buffer (formic acid—Na-formate) of different pH were predicted and compared with the measured conductivities.
- Buffer solutions of different pH and salt concentration according to recipes calculated by using “Buffalo Plus” software (GE Healthcare) were prepared by weighing the respective amounts of the components and mixing in a TECAN robot according to Table 1 below. pH and conductivities were measured. Conductivities were measured using a WTWLF 340 conductivity meter, and the results are presented in Table 1. The temperature was 22° C.
- the total conductivity C tot of the buffer solution is obtained by adding up the conductivities of all ionic species according to the formula:
- Equation (2) Using Equations (1) to (3) with the concentration values in Table 2 and the ⁇ 0 -values in Table 3, a “manual” fitting procedure was used to obtain the respective values for the Kohlrausch coefficient K in Equation (2) that gives the best fit with the measured conductivity values in Table 1. More specifically, the fit was done manually in MS Office Excel by changing the parameters one by one gradually until the best fit had been obtained for a straight line through the origin in a graph where measured conductivity is plotted against predicted conductivity.
Landscapes
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Immunology (AREA)
- Biochemistry (AREA)
- Pathology (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Electrochemistry (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Theoretical Computer Science (AREA)
- Bioinformatics & Computational Biology (AREA)
- Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Molecular Biology (AREA)
- Crystallography & Structural Chemistry (AREA)
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
Abstract
In a method of preparing a liquid solution by mixing ingredients according to a predetermined recipe, wherein at least one pair of species of the liquid solution is derived from a weak electrolyte and corresponds to an acid-base pair, the conductivity of the liquid solution is predicted by:
-
- (i) for each pair of species derived from a weak electrolyte, solving a respective equilibrium equation to calculate the actual molar concentration of each such species at equilibrium in the liquid solution,
- (ii) calculating for each ionic species of said plurality of species the molar conductivity by the formula:
Λ=Λ0 −K×Sqrt(c)
-
- wherein Λ is the molar conductivity, Λ0 is the molar conductivity at infinite dilution, c is the concentration, and K is the Kohlrausch coefficient, and wherein K and Λ0 are predetermined values for K and Λ0 for each ionic species,
- (iii) calculating the conductivity x for each ionic species by the formula:
κ=c×Λ
-
- and
- (iv) adding up the conductivities determined in step (iii) for the different ionic species to obtain a predicted conductivity of the liquid solution. A computer program product comprises instructions for causing a computer to perform the method steps.
Description
- The present invention relates to the preparation of liquid mixtures, and more particularly to the preparation of a liquid mixture, such as a buffer, wherein the conductivity of the liquid mixture is predicted.
- A buffer solution is typically an aqueous solution consisting of a mixture of a weak acid and its conjugate base, or a weak base and its conjugate acid, and has the property that the pH of the solution changes very little when a small amount of strong acid or base is added to it. Buffer solutions are therefore used as a means of keeping pH at a nearly constant value in a wide variety of chemical applications, including e.g. chromatography, filtration, etc. In general, a buffer solution may be made up of more than one weak acid and its conjugate base. For instance, a wider buffer region may be created by mixing two buffering agents with overlapping individual buffer regions.
- A buffer solution having a desired pH, and optionally also ionic strength, may be prepared by calculating the necessary amounts of the ingredients of the buffer and mixing them. While it is often necessary to solve several different equations for calculating the pH of a relatively simple mixture of a weak acid (or base) and a strong base (or acid) depending on their relative concentrations, there is commercial software available for performing such calculations. For many applications, it would be desirable to also be able to predict the expected conductivity of the buffer to be prepared. However, due to the complex nature of buffers solutions, being combinations of strong electrolytes like salts and strong acids or bases that dissociate completely and weak electrolytes like buffer substances that ionize partially, this has so far generally not been possible.
- Atkins, P., and de Paula, J., Atkins' Physical Chemistry, Oxford University Press, 2002, 7th Ed., p. 834-835 discloses calculation of the molar conductivity for a weak acid at equilibrium for the ideal case when activity coefficients are ignored.
- It is an object of the present invention to provide a method of preparing a liquid mixture which comprises predicting the conductivity of the prepared liquid mixture wherein the activities of the ions are taken into account.
- The above-mentioned object as well as other objects and advantages are achieved by the method according to the present invention for preparing a liquid mixture, especially a buffer solution, having a predicted conductivity.
- According to the present invention, the prediction of the conductivity comprises solving the exact concentrations of the different ions present in solution including the equilibrium concentration of each of the charged species of a weak electrolyte, determining the molar conductivity of each of the charged species, calculating the corresponding conductivities, and summing up all the individual conductivities to obtain the total conductivity of the liquid mixture.
- According to a basic feature of the invention, calculating the molar conductivities of the subspecies of the weak electrolyte includes applying predetermined values for the Kohlrausch coefficient used when calculating the molar conductivities. Typically, the predetermined Kohlrausch coefficient values are obtained by fitting measured conductivities to the equations used for determining molar conductivities.
- In one aspect, the present invention therefore provides a method of preparing a liquid solution by mixing ingredients according to a predetermined recipe, wherein the liquid solution comprises a plurality of species, at least one pair of species being derived from a weak electrolyte and corresponding to an acid-base pair, and wherein the conductivity of the liquid solution is predicted by:
-
- (i) for each pair of species derived from a weak electrolyte, solving a respective equilibrium equation to calculate from said predetermined recipe the actual molar concentration of each such species, including all ionic species derived from the weak electrolyte, at equilibrium in the liquid solution,
- (ii) calculating for each ionic species of said plurality of species the molar conductivity by the formula:
-
Λ=Λ0 −K×Sqrt(c) - wherein Λ is the molar conductivity, Λ0 is the molar conductivity at infinite dilution, c is the concentration of the ionic species, K is the Kohlrausch coefficient, and Sqrt is square root, and wherein K and Λ0 are obtained from a data set comprising predetermined values for K and Λ0 for each ionic species,
-
- (iii) calculating the conductivity κ for each ionic species by the formula:
-
κ=c×Λ -
- wherein c and Λ are as defined above, and
- (iv) adding up the conductivities determined in step (iii) for the different ionic species to obtain a predicted conductivity of the liquid solution.
- In a preferred embodiment, the data set comprises values for the Kohlrausch coefficient K obtained by measuring the conductivity of a solution containing the ionic species at a plurality of different concentrations of the ionic species and different pH values, and optionally different temperatures, and fitting the resulting data to the formula in step (ii) above to obtain K.
- Values for Λ0 in the data set may include previously known values, such as literature values, and/or values be obtained by the fitting.
- The Kohlrausch coefficient is typically expressed as K=A+B×Λ0, wherein A and B are temperature-dependent constants and Λ0 is as defined above. Rather than obtaining K directly in the above-mentioned fitting, values for A and B may be obtained in the fitting, and K is then calculated therefrom.
- Alternatively, the Kohlrausch coefficient may be expressed as K=A+B+w×Λ0, wherein A and B are temperature-dependent constants, w is the Onsager factor, and Λ0 is as defined above. Values for A, B and w may then be obtained in the fitting and K calculated therefrom.
- The concentration of each ionic species is preferably calculated by an algorithm comprising the equation of Debye-Hückel, wherein the ionic strength of each species is used as a weighting parameter in the calculation of an average hydrated radius of the ionic species, as will be described in more detail below.
- The above-described method may advantageously be used for controlling a buffer formulation system or an in-line dilution system. The method may also be used in screening experiments wherein conductivity is used as a design of experiment (DoE) parameter.
- Other preferred embodiments are set forth in the dependent claims.
- In another aspect, the present invention provides a computer program product comprising instructions for performing the steps of the method aspect above.
- The conductivity prediction steps in the method of the invention may further be used in a “backwards” fashion to indirectly determine pH of a solution by measuring its conductivity.
- In still another aspect, the present invention therefore provides a device for measuring pH, comprising a conductivity sensor and means for calculating pH from measured conductivity using the conductivity prediction steps of the method aspect above in a backwards calculation mode.
- Yet another aspect of the invention provides a method for predicting conductivity of a liquid solution comprising steps (i) to (iv) of the method aspect above.
- In the following, the invention will be described in more detail, by way of example only, reference being made to the accompanying drawing.
-
FIG. 1 is a graph showing measured conductivity versus predicted conductivity for a buffer solution. - Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by a person skilled in the art related to this invention. Also, the singular forms “a”, “an”, and “the” are meant to include plural reference unless it is stated otherwise.
- In order to facilitate an understanding of the disclosed invention, a number of terms will be defined below.
- Buffer
- As used herein, a buffer solution is an aqueous solution typically consisting of a mixture of a weak acid and its conjugate base or a weak base and its conjugate acid. It has the property that the pH of the solution changes very little when a small amount of strong acid or base is added to it. Buffer solutions are used as a means of keeping pH at a nearly constant value in a wide variety of chemical applications
- Titrand and Titrant
- As used herein in the context of buffer systems, a titrand is the solution (or other substance) into which a titrant is added. The titrand may, for example, be a weak acid or base, and the titrant a strong base or acid, respectively.
- Electrolyte
- A “strong” electrolyte is a solute hat completely, or almost completely, ionizes or dissociates in a solution. These ions are good conductors of electric current in the solution. Examples of strong electrolytes include, for instance, salts, such as sodium chloride; strong acids, such as hydrochloric acid, sulfuric acid, nitric acid; and strong bases, such as potassium hydroxide, sodium hydroxide, calcium hydroxide.
- A weak electrolyte, on the other hand, dissociates into ions only to a small extent, typically much less than 10%, and the vast bulk of the weak electrolyte remains in its original unionized form in solution. In general, most organic acids and their salts, and the salts of organic bases, are weak electrolytes. Examples of weak electrolytes include, for example acetic acid and citric acid and salts thereof. The amount of weak electrolyte that dissociates depends on the number of counter-ions present in the solution.
- Conductivity (Electrolytic)
- The conductivity (or specific conductance) of an electrolyte solution is a measure of its ability to conduct electricity. The SI unit of conductivity is “siemens” per meter (S/m).
- Ionic Strength
- The ionic strength of a solution is a function of the concentration of all ions in the solution (half the sum of concentration multiplied by the square of ionic charge for all ions). The ionic strength is typically given in the unit mol/dm3.
- Description of the Method of the Invention
- As mentioned above, the present invention relates to the prediction of the conductivity of a buffer solution prepared according to a desired recipe or formula.
- To prepare a buffer having a desired pH, the added amounts of moles of titrand, titrant, and salt solution, when required, can be calculated. From these mole values, the volume needed from corresponding stock solutions can in turn be calculated. Typically, such calculations are performed by suitable computer software, which is commercially available. The use of such software permits the preparation of buffers up to concentrations of 1 M or higher, as will be described in more detail below.
- However, to predict the expected conductivity of the buffer to be prepared is not trivial. While for simple electrolytes like NaCl theoretical accounts have been developed (see e.g. Chandra, A., Biswas, R., and Bagchi, B., J. Am. Chem. Soc. 1999, 121, 4082-4083), this is generally not applicable to buffers which typically are complex combinations of, on the one hand, strong electrolytes like strong acids or bases and salts which dissociate completely, and weak electrolytes like buffer substances (weak acids and bases) that ionize partially. As mentioned above in the background section, Atkins, P., and de Paula, J., 2002 describes calculation of molar conductivities for weak electrolytes for the ideal case where activity coefficients are ignored, which assumption leads in most cases to gross errors in the solution of the equilibrium equations. There are, however, currently no models addressing the general case of real solutions where the activity of the ions is taken into account and data to parametrise such models.
- Before describing the present invention in more detail, the theory behind electrolytic conductivity will be briefly described.
- From its definition, the molar conductivity is given by the equation
-
Λ=κ/c (Eq. 1) - where κ is the measured conductivity, and c is the electrolyte concentration.
- For strong electrolytes, such as salts, strong acids and strong bases, the molar conductivity is only weakly dependent on concentration and, to a good approximation, fits to the Kohlrausch equation
-
Λ=Λ0 −K×Sqrt(c) (Eq. 2) - where Λ0 is the molar conductivity at infinite dilution, or “limiting molar conductivity”, and K is the “Kohlrausch coefficient”, which depends on the nature of the specific salt in solution.
- The constant K can be expressed as A+B×Λ0 (Debye, P., Hückel, H., Phys. Z. 1924, 25, 49) or A+B×w×Λ0 (Onsager, L., Phys. Z, 1926, 27, 388-392, and 1927, 28, 277-298), where A and B are constants depending on the temperature (T), dielectric constant (ε) and charges of the ions (A also depends on the viscosity of the liquid), and w is a factor introduced by Onsager which contains the contribution from self-Brownian motion of the ion (Onsager, L., 1926, 1927 above).
- According to Equation 2 above, the molar conductivity A decreases as a function of the square root of the concentration in a linear fashion. This is an approximation that is valid in the low concentration regime. As the concentration increases, the rate of decrease of the molar conductivity A decreases so that it levels off, as can for instance be seen in FIG. 1 of the above cited reference Chandra, A., 1999. Besides, this reference also presents a microscopic theory that accounts for the constants A and B by treating the solvent as a dielectric continuum and calculating the concentration dependent ion contribution to microscopic friction.
- The limiting molar conductivity Λ0 can be decomposed into contributions from the different ions (commonly known as law of independent migration of ions):
-
Λ0=Σiνiλi (Eq. 3) - where λi is the molar ionic conductivities of ion i, and νi is the number of ions i in the formula unit of the electrolyte (e.g. 2 and 1 for Na+ and SO4 2−, respectively, in Na2SO4).
- A table of limiting molar conductivity of a limited number of different ions in solution can, for instance, be found in Harrison, R. D., Revised Nuffield Advanced Science Book of data. Longman, Harlow, 1984.
- According to the present invention there is now suggested a novel approach for predicting conductivity, which is based, on the one hand, on the assumption that even if the buffering substance itself is not a strong electrolyte, every and each one of its charged subspecies can be considered as one strong electrolyte at equilibrium, and, on the other hand, on using the equilibrium concentrations to calculate the molar conductivities of the weak electrolyte subspecies from the Kohlrausch equation above, wherein the Kohlrausch coefficient K has been determined beforehand, typically by fitting measured conductivity data at known pH values and concentrations thereto.
- Using this approach, in brief, by first determining the exact concentrations of the different ions present in solution including the concentration of each of the charged species of the buffer substance, the conductivity can then be calculated with Equations (2) and (3) above once the constants K and Λ0, or alternatively the constants A, B and Λ0 (and optionally w), have been determined. This determination can be made experimentally by measuring the conductivities of a large number of samples of the buffer at different concentrations, temperatures, and pH values, and then fitting the obtained data to Equation (2).
- With regard to the first step of determining the concentrations of the different ions present in solution, methods to do this are documented in the literature and well known to persons skilled in the art and will therefore be only briefly addressed below. Basically this includes solving the equilibrium equations for the buffering species.
- Solving Equilibrium Equations
- For a particular basic species (which may be a base B or a conjugate base A−) in equilibrium with a corresponding acidic species (which may be a conjugate acid BH+ or an acid HA, respectively) the following equation applies
-
pH=pKa+log {(basic species)/(acidic species)} (Eq. 4) - which is often referred to as the Henderson-Hasselbach equation. In this equation, the parentheses denote the activity of each species rather than the concentration, the reason being that the ions involved tend to become shielded from the environment. The activity of each ion is related to the corresponding concentration through the activity coefficient φ
-
(species)=φ[species] (Eq. 5) - In the ideal state of infinite dilution, φ becomes 1 and the activity of every ion will be equal to the corresponding concentration. However, in practice, the ionic strength is different from zero and the activity coefficients of the different species will be less than 1.
- Inserting Eq. 5 into Eq. 4 above gives the pH in terms of the concentrations instead of activities
-
pH=pKa′+log{(basic species)/(acidic species)} (Eq. 6) - where
-
pKa′=pKa+log φb−log φa (Eq. 7) - in which φa and φb are the activity coefficients for the acidic and basic species, respectively, and pKa′ is an apparent pKa value which allows the use of the measureable concentrations of the different buffer species.
- A model for these deviations is given in the so-called Debye-Hückel theory, known as
-
log φ=(AZ 2 I 0.5)/(1+0.33*108 aI 0.5) (Eq. 8) - wherein A is a constant, or rather a temperature dependent parameter ˜0.51. Using well known data, the value of A can accurately be calculated as A=0.4918+0.0007*T+0.000004*T2, where T is the temperature in degrees Celsius, Z is the charge of the ion and the quantity α, the radii of the hydrated ions (in Å), is described as the “mean distance of approach of the ions, positive or negative” in the original paper of Debye and Hückel.
- The value of pKa′ can be calculated by inserting Eq. 8 into Eq. 7 giving
-
pKa′=pKa+(AZ a 2I10.5)/(1+0.33*108 a a I 0.5)−(AZ b 2 I 0.5)/(1+0.33*108 a b I 0.5) (Eq. 9) - where the subscripts a and b specify the parameters corresponding to the acid and the base respectively, Za=charge of acidic species, Zb=charge of basic species, aa=ion size parameter of the acidic species, and ab=ion size parameter of the basic species.
- As to the parameter a, Guggenheim E. A. & Schindler, T. D., (1934) J. Phys. Chem. 33. 533 suggests an approximation of the parameter a set to 3 Å for all buffer molecules, leading to the simplified formula
-
pKa′=pKa+(AZ a 2 I 0.5)/(1+I 0.5)−(AZ b 2 I 0.5)/(1++I 0.5) (Eq. 10) - Eq. 10 above is the formula for ionic strength correction usually found in the literature.
- Using Eq. 7 (or 4) in combination with the equations for (i) the conservation of mass, (ii) the conservation of charge, and (iii) the water dissociation equilibrium, the concentrations of the acidic and basic species, respectively, of a monoprotic buffer at equilibrium may be calculated.
- However, many buffers are polyprotic, i.e. their buffer molecules can accept and give away more than one proton corresponding to more than one pKa value. The number of species in such a buffer system is always one more than the number of pKa values. The calculation of the amount of moles in each of the protonation species is equivalent to solving the equilibrium equations of each of the species with the “neighbouring” species with one more and/or one less proton and with the concentration of hydrogen atoms (the pH).
- Assume, for example, a tritropic buffer. Four protonation states or species will be defined (s1, s2, s3 and s4), and three pKa values. Three equations (corresponding to three pKa values) may then be derived directly from Eq. 6 above.
-
xx[i]=10(pH-pKa′[i]) (Eq. 11) - where each i corresponds to each pKa(i) value (i=1, 2, 3), and xx[i] are the ratios between the concentrations of the corresponding base and the corresponding acid i.e. xx[1]=[s2/s1], xx[2]=[s3/s2], xx[3]=[s4/s3].
- In addition to these three equations, an equation arises because of the conservation of mass
-
[s 1 ]+[s 2 ]+[s 3 ]+[s 4]=buffer concentration (Eq. 12) - and the conservation of charge
-
[H+]—[OH−]+Σ specific charge (s i)−titrant_charge*[titrant]−spec_charge(start_species)*[start_species]=0 (Eq. 13) - By “start_species” is meant the species of the buffer substance before mixing, i.e. the protonation state of the buffer in the can or stock solution. This protonation state is determined by the amount of counter-ions per buffer molecule because the macroscopic object, can or stock solution, has to be electrically neutral. The minus sign in front of [OH−] is due to the minus sign of the charge of the OH ions, whereas the minus sign in front of the two last terms is due to the charge of the counter-ions of the titrant and the start_species, respectively.
- Finally, there is the water dissociation equilibrium
-
[OH−][H+]=1014 (Eq. 14) - Equations (11) to (14) above imply that there are six equations with six unknowns (the four [si], [OH−] and [H+]) for the case of three pKa values, and the equilibrium concentrations of the respective acidic and basic species may thus be calculated.
- For a monoprotic buffer, for example, the calculations are simplified; s1 corresponds to the acidic species, s2 to the basic species and the concentrations of s3 and s4 are set to zero.
- A more accurate determination of the ion size parameter a in the Debye-Hückel equation 8 above is described in WO 2009/131524 A1 (the disclosure of which is incorporated by reference herein), where a is determined as the weighted mean ion size of all species contributing to the ionic strength of the liquid mixture, using the ionic strength as weighting parameter. The ion size parameter a may then be calculated as
-
- where Ii is the ionic strength, ai is the ion size parameter of species i, and I is the total ionic strength defined by
-
I=½Σ(C i Z i) (Eq. 16) - where Ci is the concentration and Zi is the charge of ion present in the solution (in units of electronic charge), which gives
-
- The ion size parameter a may also be approximated as
-
a=0.5*(mass)1/3+shell (Eq. 18). - where “shell” is fixed at a specific value for a positively charged ionic species (typically in the range of 3.9-4.1), and at another fixed value for a negatively charged species (typically in the range of 0-0.1, or 0).
- This permits the determination of the relative proportions of the constituents of buffers having salt concentrations up to or above 2 M or even 5 M.
- In the above-mentioned WO 2009/131524 A1 there is also disclosed a computer program, or software, for calculating a recipe (i.e. amounts of titrand, titrant, water and salt) to achieve a solution of a given pH and ionic strength (or, alternatively, salt concentration instead of ionic strength). Specifically, a computer program “Buffalo Plus” is disclosed which is written in Visual C++ and can be run in Windows for accurate calculation of buffer pH at ionic strengths as high as 1 M for buffer systems including phosphate, citrate, acetate and Tris.
- Such software may conveniently be used for calculating the equilibrium concentrations of acidic and basic species for the purposes of the present invention.
- Predicting Conductivity
- 1. Determination of the Kohlrausch Coefficient K
- For each pair of ionic species of a weak electrolyte to be present in a buffer of a desired composition, pH, and optionally ionic strength, the conductivity of a solution containing the weak electrolyte is measured at a plurality of different concentrations and different pH values, and optionally different temperatures. The resulting conductivities are then used together with equilibrium concentrations obtained by solving equilibrium equations as described above to find through a computational procedure of regression analysis type the values for K that give the best fit to a regression model provided by the above-mentioned equations (2) and (3), i.e.
-
Λ=Λ0 −K×Sqrt(c) (Eq. 2) -
Λ0=Σiνi (Eq. 3) - Values for limiting molar conductivities are available in the literature for a number of common ions, including those of weak electrolytes. In cases, however, where such λ0 values are not known, they may be determined in the fitting to equations (2) and (3).
- Preferably, measurements are made for a large number a different weak electrolyte ionic species at different concentrations and pH to produce a catalogue or database of values for the coefficient K, optionally together with Λ0/λ0 values.
- 2. Calculation of Conductivities
- Values for K and Λ0 obtained as described above are used to calculate the molar conductivity of each ion (strong electrolyte ions as well as weak electrolyte ions) present in the desired buffer solution using equation (2) above. The corresponding conductivity of each ion is then determined from the previously mentioned equation (1), i.e.
-
Λ=κ/c (Eq. 1), - Finally, the conductivities of all ions are added up to give the predicted conductivity of the desired buffer solution.
- The method of the invention may be implemented by software run on an electrical data processing device, such as a computer. Such software may be provided to the computer on any suitable computer-readable medium, including a record medium, a read-only memory, or an electrical or optical signal which may be conveyed via electrical or optical cable or by radio or other means.
- Predicting conductivities as described above may be used for several purposes. Exemplary applications include the use for controlling a buffer formulation system or an in-line dilution system, wherein a buffer solution of predetermined conductivity, pH (and potentially other parameters) is prepared by mixing relative component proportions of at least one each of:
-
- a buffer;
- an acid or a base;
- a solvent;
- and optionally
- a salt,
and wherein the relative component proportions are determined using predicted conductivity according to the present method and using predicted pH and Ionic strength in accordance with the method disclosed in WO 2009/131524. According to one embodiment the relative component proportions are determined using an iterative procedure wherein the conductivity and other control parameters are used to iteratively determine relative component proportions.
- Such prediction of conductivities may also be used in screening experiments wherein conductivity is used as a design of experiment (DoE) parameter. Still another application is for determining the exact concentration of a stock solution of a weak electrolyte. Yet another application is for indirect measurement of pH by measuring conductivity and calculating pH therefrom.
- The invention will now be described in more detail, by way of example only, by the following non-limiting Example.
- In this example, conductivities for a plurality of compositions of a formate buffer (formic acid—Na-formate) of different pH were predicted and compared with the measured conductivities.
- Buffer solutions of different pH and salt concentration according to recipes calculated by using “Buffalo Plus” software (GE Healthcare) were prepared by weighing the respective amounts of the components and mixing in a TECAN robot according to Table 1 below. pH and conductivities were measured. Conductivities were measured using a WTWLF 340 conductivity meter, and the results are presented in Table 1. The temperature was 22° C.
-
TABLE 1 NaCl— Buffer- Measured conc. conc. conduc- HCOOH Na+HCOO− (mole/ (mole/ tivity pH (mole) (mole) dm3) dm3) (mS/cm) 3.1 0.040791 0.009225 0 0.050016 1.32 3.97 0.0165223 0.0328 0 0.049322 3.17 2.995 0.20755 0.044107 0 0.251657 4.7 4.04 0.0651762 0.186358 0 0.251534 15.38 2.985 0.410651 0.093503 0 0.504154 9.1 4.06 0.115571 0.389142 0 0.504713 29.2 2.99 0.605233 0.148559 0 0.753792 13.62 4.095 0.15289 0.60327 0 0.75616 42.3 3.005 0.0410208 0.009768 0.251337 0.050789 27.3 3.9 0.0163126 0.034522 0.251337 0.050835 28.9 2.975 0.0406123 0.009823 0.502674 0.050435 51.1 3.885 0.015963 0.035611 0.503535 0.051574 52.5 3.035 0.0395468 0.012015 0.754871 0.051561 72.4 3.93 0.0140662 0.036802 0.754011 0.050868 73.8 3.005 0.200951 0.054151 0.250476 0.255102 30.4 4.04 0.0602869 0.192086 0.251337 0.252373 39.2 2.985 0.196017 0.054781 0.503535 0.250797 53.5 3.99 0.0628657 0.189574 0.502674 0.25244 61.3 2.995 0.193722 0.058587 0.754011 0.25231 74.3 4.035 0.0556598 0.194628 0.754011 0.250287 80.2 - Using the Buffalo software, pKa′ values and the equilibrium concentrations of HCOOH and Na+COO−, were then calculated, as well as the concentrations of Na+, H+, OH− and Cl−, and the results are given in Table 2 below.
-
TABLE 2 HCOOH COO− Na+ H+ OH− Cl− pKa′ (mole/dm3) (mole) (mole/dm3) (mole/dm3) (mole/dm3) (mole/dm3) 3.7012 0.0399966 0.010019 0.00739624 0.000794328 1.25893E−11 0 3.66796 0.0164151 0.0329068 0.03402269 0.000107152 0.000107152 0 3.65565 0.206539 0.0451185 0.04511705 0.001011579 9.88553E−12 0 3.58292 0.065085 0.186449 0.18268706 9.12011E−05 1.09648E−10 0 3.62177 0.409616 0.094538 0.09763033 0.001035142 9.66051E−12 0 3.53232 0.115484 0.389229 0.37794773 8.70964E−05 1.14815E−10 0 3.59631 0.60421 0.149582 0.15384174 0.001023293 9.77237E−12 0 3.49858 3.49858 0.60335 0.57912539 8.03526E−05 1.24451E−10 0 3.57575 0.0400323 0.0107562 0.26169173 0.000988553 1.01158E−11 0.251337 3.56948 0.0161867 0.0346478 0.28905781 0.000125893 7.94328E−11 0.251337 3.53546 0.039553 0.0108823 0.51376836 0.001059254 9.44061E−12 0.502674 3.53138 0.0158327 0.0357414 0.54199543 0.000130317 7.67361E−11 0.503535 3.51002 0.0386243 0.0129371 0.76818423 0.000922571 1.08393E−11 0.754871 3.50728 0.0139487 0.036919 0.79395068 0.00011749 8.51138E−11 0.754011 3.56449 0.199963 0.0551394 0.3066874 0.000988553 1.01158E−11 0.250476 3.53586 0.0601957 0.192177 0.4414203 9.12011E−05 1.09648E−10 0.251337 3.52824 0.194982 0.0558156 0.56418415 0.001035142 9.66051E−12 0.503535 3.50969 0.0627634 0.189676 0.69719504 0.000102329 9.77237E−11 0.502674 3.50467 0.192711 0.0595988 0.81835827 0.001011579 9.88553E−12 0.754011 3.49041 0.0555676 0.19472 0.95075092 9.22571E−05 1.08393E−10 0.754011 - The total conductivity Ctot of the buffer solution is obtained by adding up the conductivities of all ionic species according to the formula:
-
C tot =C s1 +C s2 +C Na+ +C H+ +C OH− +C Cl− - where the subscripts s1 and s2 denote HCOOH and Na+COO−, respectively.
- The conductivity values are in turn obtained through equations (1) to (3) above, i.e.
-
Λ=κ/c (Eq. 1) -
Λ=Λ0 −K×Sqrt(c) (Eq. 2) -
Λ0=Σiνiλi (Eq. 3) - Λ0-values (in mSm2mol−1) for the respective species are obtained from Atkins' Physical Chemistry, 7th Ed. (supra) and listed in Table 3 below.
-
TABLE 3 Species Λ0 (mSm2mol−1) HCOOH 0 HCOO−, 5.46 Na+ 5.01 H+ 34.96 OH− 19.91 Cl− 7.635 - Using Equations (1) to (3) with the concentration values in Table 2 and the Λ0-values in Table 3, a “manual” fitting procedure was used to obtain the respective values for the Kohlrausch coefficient K in Equation (2) that gives the best fit with the measured conductivity values in Table 1. More specifically, the fit was done manually in MS Office Excel by changing the parameters one by one gradually until the best fit had been obtained for a straight line through the origin in a graph where measured conductivity is plotted against predicted conductivity.
- The resulting K-values for HCOO−, Na+ and Cl− are indicated in Table 4 below together with the Λ0 values. The K-values for HCOOH, H+ and OH− were assumed to be zero due to low concentration.
-
TABLE 4 HCOOH HCOO− Na+ H+ OH− Cl− Λ o0 5.46 5.01 34.96 19.91 7.365 K 0 0.91 3.6 0 0 0.2 - Corresponding conductivity values for the different species concentrations of each species indicated in Table 2, and the predicted total conductivity of the buffer solution, together with the measured conductivity values from Table 1, are given in Table 5 below.
-
TABLE 5 Conductivity (mS/cm) CHCOOH CCOO− CNa+ CH+ CCl− Ctot predicted Ctot measured 0 0.537911 0.347652 0.277697 0 1.163261 1.32 0 1.74239 1.478616 0.03746 0 3.258467 3.17 0 2.376259 1.915369 0.353648 0 4.645275 4.7 0 9.447491 6.3416 0.031884 0 15.82097 15.38 0 4.897259 3.793084 0.361886 0 9.052229 9.1 0 19.04212 10.57049 0.030449 0 29.64306 29.2 0 7.640723 5.535201 0.357743 0 13.53367 13.62 0 28.67814 13.1484 0.028091 0 41.85464 42.3 0 0.577137 8.291416 0.345598 18.93757177 28.15172 27.3 0 1.833081 8.887061 0.044012 18.93757177 29.70173 28.9 0 0.583843 12.48254 0.370315 37.66637313 51.10307 51.1 0 1.889991 12.78929 0.045559 37.73027836 52.45512 52.5 0 0.692975 14.24782 0.322531 56.32268699 71.58601 72.4 0 1.951224 14.30905 0.041074 56.25926694 72.56062 73.8 0 2.892787 9.250744 0.345598 18.87312826 31.36226 30.4 0 9.726221 11.55717 0.031884 18.93757177 40.25285 39.2 0 2.927534 13.00987 0.361886 37.73027836 54.02956 53.5 0 9.604583 13.97224 0.035774 37.66637313 61.27897 61.3 0 3.121692 14.34849 0.353648 56.25926694 74.0831 74.3 0 9.849801 14.25905 0.032253 56.25926694 80.40037 80.2 - The corresponding graph mentioned above where measured conductivity is plotted against predicted conductivity is shown in
FIG. 1 . - As apparent from Table 5 and
FIG. 1 , a good fit was obtained, which demonstrates that the procedure outlined above may be used for predicting the conductivity that will be obtained following a predetermined recipe for a buffer composition. - While the above described fitting procedure was done manually, it may conveniently be performed by an algorithm designed therefor. Based on the teaching above, the design of such a suitable algorithm may readily be made by a person skilled in the art.
- The present invention is not limited to the above-described preferred embodiments. Various alternatives, modifications and equivalents may be used. Therefore, the above embodiments should not be taken as limiting the scope of the invention, which is defined by the appending claims.
Claims (12)
1. A method of preparing a liquid solution by mixing ingredients according to a predetermined recipe, wherein the liquid solution comprises a plurality of species, at least one pair of species being derived from a weak electrolyte and corresponding to an acid-base pair, and wherein the conductivity of the liquid solution is predicted by:
(i) for each pair of species derived from a weak electrolyte, solving a respective equilibrium equation to calculate from said predetermined recipe the actual molar concentration of each such species, including all ionic species derived from the weak electrolyte, at equilibrium in the liquid solution,
(ii) calculating for each ionic species of said plurality of species the molar conductivity by the formula:
Λ=Λ0 −K×Sqrt(c)
Λ=Λ0 −K×Sqrt(c)
wherein Λ is the molar conductivity, Λ0 is the molar conductivity at infinite dilution, c is the concentration of the ionic species, and K is the Kohlrausch coefficient, and wherein K and Λ0 are obtained from a data set comprising predetermined values for K and Λ0 for each ionic species,
(iii) calculating the conductivity K for each ionic species by the formula:
κ=×Λ
κ=×Λ
wherein c and Λ are as defined above, and
(iv) adding up the conductivities determined in step (iii) for the different ionic species to obtain a predicted conductivity of the liquid solution.
2. The method of claim 1 , wherein the data set comprises values for the Kohlrausch coefficient K obtained by measuring the conductivity of a solution containing the ionic species at a plurality of different concentrations of the ionic species and different pH values, and optionally different temperatures, and fitting the resulting data to the formula in step (ii) of claim 1 to obtain K.
3. The method of claim 2 , wherein said data set further comprises values for Λ0 obtained by said fitting.
4. The method of claim 2 , wherein the Kohlrausch coefficient is expressed as K=A+B×Λ0, wherein A and B are temperature-dependent constants and Λ0 is as previously defined, and wherein the predetermined values for K in the data set are calculated from values for A and B obtained by said fitting.
5. The method of claim 2 , wherein the Kohlrausch coefficient is expressed as K=A+B+w×Λ0, wherein A and B are temperature-dependent constants, w is the Onsager factor, and Λ0 is as previously defined, and wherein the predetermined values for K in the data set are calculated from values for A, B and w obtained by said fitting.
6. The method of claim 1 , wherein the concentration of each ionic species is calculated by an algorithm comprising the equation of Debye-Hückel, wherein the ionic strength of each species is used as a weighting parameter in the calculation of an average hydrated radius of the ionic species.
7. The method of claim 1 , wherein the liquid solution is a buffer.
8. The method of claim 1 , which is computer-implemented.
9-11. (canceled)
12. A computer program product comprising instructions for causing a computer to perform the method steps of claim 1 .
13. A device for measuring pH, comprising a conductivity sensor and means for calculating pH from measured conductivity using the conductivity prediction steps of claim 1 , in a backwards calculation mode.
14. A method of predicting the conductivity of a liquid solution containing ingredients according to a predetermined recipe, wherein the liquid solution comprises a plurality of species, at least one pair of species being derived from a weak electrolyte and corresponding to an acid-base pair, comprising the steps of:
(i) for each pair of species derived from a weak electrolyte, solving a respective equilibrium equation to calculate from said predetermined recipe the actual molar concentration of each such species, including all ionic species derived from the weak electrolyte, at equilibrium in the liquid solution,
(ii) calculating for each ionic species of said plurality of species the molar conductivity by the formula:
Λ=Λ0 −K×Sqrt(c)
Λ=Λ0 −K×Sqrt(c)
wherein Λ is the molar conductivity, Λ0 is the molar conductivity at infinite dilution, c is the concentration of the ionic species, and K is the Kohlrausch coefficient, and wherein K and Λ0 are obtained from a data set comprising predetermined values for K and Λ0 for each ionic species,
(iii) calculating the conductivity K for each ionic species by the formula:
κ=c×Λ
κ=c×Λ
wherein c and Λ are as defined above, and
(iv) adding up the conductivities determined in step (iii) for the different ionic species to obtain a predicted conductivity of the liquid solution.
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
SE1051344-8 | 2010-12-17 | ||
SE1051344 | 2010-12-17 | ||
PCT/SE2011/051513 WO2012082061A1 (en) | 2010-12-17 | 2011-12-15 | Method for predicting the conductivity of a liquid mixture |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/SE2011/051513 A-371-Of-International WO2012082061A1 (en) | 2010-12-17 | 2011-12-15 | Method for predicting the conductivity of a liquid mixture |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US16/130,092 Division US11670402B2 (en) | 2010-12-17 | 2018-09-13 | Method for predicting the conductivity of a liquid mixture |
Publications (1)
Publication Number | Publication Date |
---|---|
US20130270492A1 true US20130270492A1 (en) | 2013-10-17 |
Family
ID=46244974
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US13/994,916 Abandoned US20130270492A1 (en) | 2010-12-17 | 2011-12-15 | Method for predicting the conductivity of a liquid mixture |
US16/130,092 Active 2033-08-13 US11670402B2 (en) | 2010-12-17 | 2018-09-13 | Method for predicting the conductivity of a liquid mixture |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US16/130,092 Active 2033-08-13 US11670402B2 (en) | 2010-12-17 | 2018-09-13 | Method for predicting the conductivity of a liquid mixture |
Country Status (6)
Country | Link |
---|---|
US (2) | US20130270492A1 (en) |
EP (1) | EP2651547B1 (en) |
JP (1) | JP5923518B2 (en) |
KR (1) | KR101950524B1 (en) |
CN (2) | CN110057874A (en) |
WO (1) | WO2012082061A1 (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20160282318A1 (en) * | 2015-03-27 | 2016-09-29 | Ge Healthcare Bio-Sciences Ab | Method for baseline correction in a chromatogram |
US11670402B2 (en) | 2010-12-17 | 2023-06-06 | Cytiva Sweden Ab | Method for predicting the conductivity of a liquid mixture |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3052226B1 (en) * | 2013-09-30 | 2019-07-03 | GE Healthcare Bio-Sciences AB | Method for preparation of liquid mixtures |
CN106990141A (en) * | 2017-05-19 | 2017-07-28 | 内蒙古岱海发电有限责任公司 | A kind of stove water conductivity computational methods and system |
GB201709531D0 (en) * | 2017-06-15 | 2017-08-02 | Ge Healthcare Bio Sciences Ab | Method and apparatus for determining one or more buffer composition recipes |
CN107764870A (en) * | 2017-11-14 | 2018-03-06 | 东莞市循美环保科技有限公司 | Cutting fluid concentration method of testing |
CN108896498B (en) * | 2018-04-25 | 2020-07-21 | 中国农业大学 | Method for accounting ion EC contribution rate and conductivity of nutrient solution based on ion activity |
US20220128523A1 (en) * | 2019-02-28 | 2022-04-28 | Cytiva Sweden Ab | Control of a Buffer Preparation Process |
JP7489172B2 (en) | 2019-05-14 | 2024-05-23 | 三浦工業株式会社 | Burner |
CN113707225B (en) * | 2020-05-21 | 2024-03-15 | 中国科学院过程工程研究所 | Method for predicting metal separation capacity in solution based on metal ion morphology and application thereof |
CN111681713B (en) * | 2020-06-10 | 2023-03-31 | 重庆邮电大学 | Model construction method for polymer molecule electrical property containing decoherence and application thereof |
CN113465660B (en) * | 2021-05-25 | 2022-07-05 | 湖南大学 | Non-contact temperature measurement and material component detection device and method based on conductivity |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6221250B1 (en) * | 1996-02-29 | 2001-04-24 | Amersham Pharmacia Biotech Ab | Method of preparing a liquid mixture |
US20080279038A1 (en) * | 2003-10-17 | 2008-11-13 | Louis Bellafiore | Multi-stage accurate blending system and method |
WO2009131524A1 (en) * | 2008-04-21 | 2009-10-29 | Ge Healthcare Bio-Sciences Ab | Preparation of liquid mixtures |
WO2011037530A1 (en) * | 2009-09-25 | 2011-03-31 | Ge Healthcare Bio-Sciences Ab | Method and system for preparation of liquid mixtures |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0718831B2 (en) * | 1987-07-17 | 1995-03-06 | 財団法人電力中央研究所 | Method for controlling circulating liquid in liquid circulation plant, especially method and device for controlling primary cooling water in pressurized water reactor |
CN101587150B (en) * | 2006-08-30 | 2011-04-13 | 梅特勒-托利多仪器(上海)有限公司 | Solution conductivity measuring method |
SI22818A (en) * | 2008-05-13 | 2009-12-31 | KovaÄŤ Kralj Anita | Procedure of determination of weak electrolyte concentrations |
CN101726518B (en) * | 2009-12-24 | 2012-09-26 | 福建师范大学 | Method for measuring electrical conductivity of solution |
CN110057874A (en) | 2010-12-17 | 2019-07-26 | 通用电气健康护理生物科学股份公司 | The method for predicting the conductivity of liquid mixture |
-
2011
- 2011-12-15 CN CN201910026056.7A patent/CN110057874A/en active Pending
- 2011-12-15 EP EP11848855.0A patent/EP2651547B1/en active Active
- 2011-12-15 WO PCT/SE2011/051513 patent/WO2012082061A1/en active Application Filing
- 2011-12-15 JP JP2013544430A patent/JP5923518B2/en active Active
- 2011-12-15 KR KR1020137018564A patent/KR101950524B1/en active IP Right Grant
- 2011-12-15 CN CN2011800607061A patent/CN103249478A/en active Pending
- 2011-12-15 US US13/994,916 patent/US20130270492A1/en not_active Abandoned
-
2018
- 2018-09-13 US US16/130,092 patent/US11670402B2/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6221250B1 (en) * | 1996-02-29 | 2001-04-24 | Amersham Pharmacia Biotech Ab | Method of preparing a liquid mixture |
US20080279038A1 (en) * | 2003-10-17 | 2008-11-13 | Louis Bellafiore | Multi-stage accurate blending system and method |
WO2009131524A1 (en) * | 2008-04-21 | 2009-10-29 | Ge Healthcare Bio-Sciences Ab | Preparation of liquid mixtures |
WO2011037530A1 (en) * | 2009-09-25 | 2011-03-31 | Ge Healthcare Bio-Sciences Ab | Method and system for preparation of liquid mixtures |
Non-Patent Citations (2)
Title |
---|
Hirokazu Okamoto et al. "Theory and Computer Programs for Calculating Solution pH, Buffer, Formula, and Buffer Capacity for multiple Component System at a Given Ionic Strength and Temperature" Pharmaceutical Research Volume 14 No. 3 Pages 299-302 1997 * |
Peter Atkins and Julio de Paula "Physical Chemistry" 2002 Freeman 7th Ed. * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11670402B2 (en) | 2010-12-17 | 2023-06-06 | Cytiva Sweden Ab | Method for predicting the conductivity of a liquid mixture |
US20160282318A1 (en) * | 2015-03-27 | 2016-09-29 | Ge Healthcare Bio-Sciences Ab | Method for baseline correction in a chromatogram |
US10598640B2 (en) * | 2015-03-27 | 2020-03-24 | Ge Healthcare Bio-Sciences Ab | Method for baseline correction in a chromatogram |
Also Published As
Publication number | Publication date |
---|---|
WO2012082061A1 (en) | 2012-06-21 |
KR101950524B1 (en) | 2019-02-20 |
EP2651547A4 (en) | 2016-05-11 |
EP2651547B1 (en) | 2021-02-17 |
US11670402B2 (en) | 2023-06-06 |
CN110057874A (en) | 2019-07-26 |
KR20140003466A (en) | 2014-01-09 |
US20190042704A1 (en) | 2019-02-07 |
JP2014501385A (en) | 2014-01-20 |
EP2651547A1 (en) | 2013-10-23 |
JP5923518B2 (en) | 2016-05-24 |
CN103249478A (en) | 2013-08-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11670402B2 (en) | Method for predicting the conductivity of a liquid mixture | |
Snyder et al. | Estimating the electrical conductivity of cement paste pore solutions from OH−, K+ and Na+ concentrations | |
Fuoss et al. | Properties of electrolytic solutions. II. The evaluations of Λ0 and of K for incompletely dissociated electrolytes | |
Tuttle et al. | A new cluster pair method of determining absolute single ion solvation energies demonstrated in water and applied to ammonia | |
EP3201618B1 (en) | Method for predicting the solubility of a buffer | |
Gomaa et al. | Thermodynamics of the solvation of CaSO4 in mixed DMF-H2O at 301.15 K | |
EP3052226B1 (en) | Method for preparation of liquid mixtures | |
Ribeiro et al. | Diffusion coefficients and electrical conductivities for calcium chloride aqueous solutions at 298.15 K and 310.15 K | |
Hernández-Luis et al. | Activity coefficients of KCl in PEG 4000+ water mixtures at 288.15, 298.15 and 308.15 K | |
EP3212300B1 (en) | Method for predicting the dynamic ph range of a buffer | |
Partanen et al. | Comparison of different methods for calculation of the stoichiometric dissociation constant of acetic acid from results of potentiometric titrations at 298.15 K in aqueous sodium or potassium chloride solutions | |
Partanen et al. | Determination of stoichiometric dissociation constant of n-butyric acid in aqueous sodium or potassium chloride solutions at 298.15 K | |
Michałowska-Kaczmarczyk et al. | Dynamic buffer capacity versus alkalinity. Formulation in terms of Simms constants idea | |
Rafols et al. | Ionic equilibria in aqueous organic solvent mixtures. Speciation of hydrofluoric acid in several ethanol/water solutions | |
Bijster et al. | Direct potentiometric determination of sodium ion in blood. I. Potentiometric response in simple sodium chloride solutions. | |
Gurzyński et al. | Potentiometric and ab initio studies of acid–base interactions of substituted 4-halo (Cl, Br) pyridine N-oxide systems | |
Wu et al. | Binding of Halide Ions to Bovine Serum Albumin and Hemoglobin: Studied with Ion-Selective Electrodes | |
Partanen et al. | Re-evaluation of the first and second stoichiometric dissociation constants of oxalic acid at temperatures from 0 to 60 C in aqueous oxalate buffer solutions with or without sodium or potassium chloride | |
Néher-Neumann | The Liquid Junction Potential in Potentiometric Titrations. IX. The Calculation of Potentials Across Liquid Junctions of the Type, AY| AY+ BY z (B)+ HY+ A y L, for Emf Cells where Strong or Weak Complexes are Formed, Respectively, at [Y−]= C mol⋅ L− 1, Constant and− log 10 [H+]≤ 7 | |
Prokunin et al. | Metrological support of ion-selective electrodes | |
Zevatskii et al. | Protolytic equilibria and density of solutions of weak electrolytes | |
Nechepurenko et al. | Standardization of acidity scale in a toluene–ethanol medium | |
Néher-Neumann | The Liquid Junction Potential in Potentiometric Titrations. X. The Calculation of Potentials Across Liquid Junctions of the Type AY| AY+ HY+ BY z (B)+ A y L for Emf Cells Where Strong or Weak Complexes are Formed at Constant Ionic Strength and− log 10 [H+]≤ 7 | |
Sergievskii et al. | On the adequacy of modeling the concentration dependences of the activity coefficients for the components of solutions | |
Domonkos et al. | Thermal Diffusion of Weak Electrolytes. I. A Model for the Determination of the Soret Coefficient of Weak Electrolytes |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: GE HEALTHCARE BIO-SCIENCES AB, SWEDEN Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:CARREDANO, ENRIQUE;REEL/FRAME:030624/0412 Effective date: 20120208 |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: FINAL REJECTION MAILED |
|
STCB | Information on status: application discontinuation |
Free format text: ABANDONED -- FAILURE TO RESPOND TO AN OFFICE ACTION |