research papers\(\def\hfill{\hskip 5em}\def\hfil{\hskip 3em}\def\eqno#1{\hfil {#1}}\)

Journal logoJOURNAL OF
APPLIED
CRYSTALLOGRAPHY
ISSN: 1600-5767

Raoult's law revisited: accurately predicting equilibrium relative humidity points for humidity control experiments

CROSSMARK_Color_square_no_text.svg

aDepartment of Physics, University of Oxford, Keble Road, Oxford OX1 3RH, UK, bDepartment of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, UK, cEuropean Molecular Biology Laboratory, Grenoble Outstation, 71 avenue des Martyrs, CS 90181, Grenoble F-38042, France, and dUnit for Virus Host Cell Interactions, Université Grenoble Alpes–EMBL–CNRS, 71 avenue des Martyrs, CS 90181, Grenoble F-38042, France
*Correspondence e-mail: michael.bowler@physics.ox.ac.uk, mbowler@embl.fr

Edited by A. R. Pearson, Universität Hamburg, Germany (Received 22 November 2016; accepted 7 March 2017; online 29 March 2017)

The humidity surrounding a sample is an important variable in scientific experiments. Biological samples in particular require not just a humid atmosphere but often a relative humidity (RH) that is in equilibrium with a stabilizing solution required to maintain the sample in the same state during measurements. The controlled dehydration of macromolecular crystals can lead to significant increases in crystal order, leading to higher diffraction quality. Devices that can accurately control the humidity surrounding crystals while monitoring diffraction have led to this technique being increasingly adopted, as the experiments become easier and more reproducible. Matching the RH to the mother liquor is the first step in allowing the stable mounting of a crystal. In previous work [Wheeler, Russi, Bowler & Bowler (2012[Wheeler, M. J., Russi, S., Bowler, M. G. & Bowler, M. W. (2012). Acta Cryst. F68, 111-114.]). Acta Cryst. F68, 111–114], the equilibrium RHs were measured for a range of concentrations of the most commonly used precipitants in macromolecular crystallography and it was shown how these related to Raoult's law for the equilibrium vapour pressure of water above a solution. However, a discrepancy between the measured values and those predicted by theory could not be explained. Here, a more precise humidity control device has been used to determine equilibrium RH points. The new results are in agreement with Raoult's law. A simple argument in statistical mechanics is also presented, demonstrating that the equilibrium vapour pressure of a solvent is proportional to its mole fraction in an ideal solution: Raoult's law. The same argument can be extended to the case where the solvent and solute molecules are of different sizes, as is the case with polymers. The results provide a framework for the correct maintenance of the RH surrounding a sample.

1. Introduction

Sample environments that control relative humidity (RH) are important in many experiments where a wide variety of samples require specific RH values to maintain sample integrity or RH is a parameter to be varied. Humidity control has been an important parameter in the study of lipid bilayers (Lin et al., 2007[Lin, J.-M., Lin, T.-L. & Jeng, U.-S. (2007). J. Appl. Cryst. 40, s367-s372.]) and amyloid fibres (McDonald et al., 2008[McDonald, M., Kendall, A., Tanaka, M., Weissman, J. S. & Stubbs, G. (2008). J. Appl. Cryst. 41, 206-209.]), and in small-molecule crystallography (Mo & Ramsøskar, 2009[Mo, F. & Ramsøskar, K. (2009). J. Appl. Cryst. 42, 531-534.]), coherent X-ray diffraction microscopy of cells (Takayama & Nakasako, 2012[Takayama, Y. & Nakasako, M. (2012). Rev. Sci. Instrum. 83, 054301.]) and serial crystallography (Roedig et al., 2016[Roedig, P., Duman, R., Sanchez-Weatherby, J., Vartiainen, I., Burkhardt, A., Warmer, M., David, C., Wagner, A. & Meents, A. (2016). J. Appl. Cryst. 49, 968-975.]). In biological crystallography, changing the RH can sometimes induce phase changes in crystals of macromolecules with a concomitant improvement in the quality of observed diffraction. This has been observed since the earliest days of macromolecular crystallography (Berthou et al., 1972[Berthou, J., Cesbron, F. & Laurent, A. (1972). J. Mol. Biol. 71, 809-813.]; Einstein & Low, 1962[Einstein, J. R. & Low, B. W. (1962). Acta Cryst. 15, 32-34.]; Huxley & Kendrew, 1953[Huxley, H. E. & Kendrew, J. C. (1953). Acta Cryst. 6, 76-80.]; Perutz, 1946[Perutz, M. F. (1946). Trans. Faraday Soc. 42, B187.]) and is most easily effected by altering the molar fraction of water in the crystal solution or by changing the RH of the air surrounding a crystal. Many successful examples are given in the literature (Adachi et al., 2009[Adachi, H., Umena, Y., Enami, I., Henmi, T., Kamiya, N. & Shen, J.-R. (2009). BBA Bioenergetics, 1787, 121-128.]; Bowler et al., 2006[Bowler, M. W., Montgomery, M. G., Leslie, A. G. W. & Walker, J. E. (2006). Acta Cryst. D62, 991-995.]; Cramer et al., 2000[Cramer, P., Bushnell, D. A., Fu, J. H., Gnatt, A. L., Maier-Davis, B., Thompson, N. E., Burgess, R. R., Edwards, A. M., David, P. R. & Kornberg, R. D. (2000). Science, 288, 640-649.]; Fratini et al., 1982[Fratini, A. V., Kopka, M. L., Drew, H. R. & Dickerson, R. E. (1982). J. Biol. Chem. 257, 4686-4707.]; Gupta et al., 2010[Gupta, V., Gupta, R. K., Khare, G., Salunke, D. M., Surolia, A. & Tyagi, A. K. (2010). PLoS One, 5, e9222.]; Heras et al., 2003[Heras, B., Edeling, M. A., Byriel, K. A., Jones, A., Raina, S. & Martin, J. L. (2003). Structure, 11, 139-145.]; Hu et al., 2011[Hu, N. J., Iwata, S., Cameron, A. D. & Drew, D. (2011). Nature, 478, 408-411.]; Kadlec et al., 2011[Kadlec, J., Hallacli, E., Lipp, M., Holz, H., Sanchez-Weatherby, J., Cusack, S. & Akhtar, A. (2011). Nat. Struct. Mol. Biol. 18, 142-149.]; Kuo et al., 2003[Kuo, A., Bowler, M. W., Zimmer, J., Antcliff, J. F. & Doyle, D. A. (2003). J. Struct. Biol. 141, 97-102.]; Nakamura et al., 2007[Nakamura, A., Wada, C. & Miki, K. (2007). Acta Cryst. F63, 346-349.]; Sam et al., 2006[Sam, M. D., Abbani, M. A., Cascio, D., Johnson, R. C. & Clubb, R. T. (2006). Acta Cryst. F62, 825-828.]; Vijayalakshmi et al., 2008[Vijayalakshmi, L., Krishna, R., Sankaranarayanan, R. & Vijayan, M. (2008). Proteins Struct. Function Bioinform. 71, 241-249.]; Yap et al., 2007[Yap, T. L., Chen, Y. L., Xu, T., Wen, D., Vasudevan, S. G. & Lescar, J. (2007). Acta Cryst. F63, 78-83.]; Zerrad et al., 2011[Zerrad, L., Merli, A., Schröder, G. F., Varga, A., Gráczer, E., Pernot, P., Round, A., Vas, M. & Bowler, M. W. (2011). J. Biol. Chem. 286, 14040-14048.]). Several specific devices have been developed to control the humidity surrounding a crystal (Einstein, 1961[Einstein, J. R. (1961). J. Sci. Instrum. 38, 449-451.]; Sjögren et al., 2002[Sjögren, T., Carlsson, G., Larsson, G., Hajdu, A., Andersson, C., Pettersson, H. & Hajdu, J. (2002). J. Appl. Cryst. 35, 113-116.]; Pickford et al., 1993[Pickford, M. G., Garman, E. F., Jones, E. Y. & Stuart, D. I. (1993). J. Appl. Cryst. 26, 465-466.]) with modern devices mounted at X-ray sources or synchrotron beamlines (Kiefersauer et al., 2000[Kiefersauer, R., Than, M. E., Dobbek, H., Gremer, L., Melero, M., Strobl, S., Dias, J. M., Soulimane, T. & Huber, R. (2000). J. Appl. Cryst. 33, 1223-1230.]; Russi et al., 2011[Russi, S., Juers, D. H., Sanchez-Weatherby, J., Pellegrini, E., Mossou, E., Forsyth, V. T., Huet, J., Gobbo, A., Felisaz, F., Moya, R., McSweeney, S. M., Cusack, S., Cipriani, F. & Bowler, M. W. (2011). J. Struct. Biol. 175, 236-243.]; Sanchez-Weatherby et al., 2009[Sanchez-Weatherby, J., Bowler, M. W., Huet, J., Gobbo, A., Felisaz, F., Lavault, B., Moya, R., Kadlec, J., Ravelli, R. B. G. & Cipriani, F. (2009). Acta Cryst. D65, 1237-1246.]). The ability to change the RH while characterizing changes via diffraction allows any changes undergone by the crystal to be observed in real time and increases the chances of characterizing a beneficial phase change.

The HC1 humidity control device was developed at the EMBL Grenoble to be a user-friendly device compatible with a complex beamline environment (Sanchez-Weatherby et al., 2009[Sanchez-Weatherby, J., Bowler, M. W., Huet, J., Gobbo, A., Felisaz, F., Lavault, B., Moya, R., Kadlec, J., Ravelli, R. B. G. & Cipriani, F. (2009). Acta Cryst. D65, 1237-1246.]). It produces an air stream with a controlled RH using a dispensing nozzle, in the same manner as cryostream devices produce a nitro­gen flow at 100 K, and is therefore easy to integrate with most diffractometers. It supplies a stream of humid air at an RH determined by a dew point controller acting on a water-saturated air supply. The device is now installed at laboratories and synchrotrons across the world (Bowler, Mueller et al., 2015[Bowler, M. W., Mueller, U., Weiss, M. S., Sanchez-Weatherby, J., Sorensen, T. L. M., Thunnissen, M. M. G. M., Ursby, T., Gobbo, A., Russi, S., Bowler, M. G., Brockhauser, S., Svensson, O. & Cipriani, F. (2015). Cryst. Growth Des. 15, 1043-1054.]), resulting in many successful experiments (Hu et al., 2011[Hu, N. J., Iwata, S., Cameron, A. D. & Drew, D. (2011). Nature, 478, 408-411.]; Kadlec et al., 2011[Kadlec, J., Hallacli, E., Lipp, M., Holz, H., Sanchez-Weatherby, J., Cusack, S. & Akhtar, A. (2011). Nat. Struct. Mol. Biol. 18, 142-149.]; Malinauskaite et al., 2014[Malinauskaite, L., Quick, M., Reinhard, L., Lyons, J. A., Yano, H., Javitch, J. A. & Nissen, P. (2014). Nat. Struct. Mol. Biol. 21, 1006-1012.]; Oliete et al., 2013[Oliete, R., Pous, J., Rodríguez-Puente, S., Abad-Zapatero, C. & Guasch, A. (2013). Acta Cryst. D69, 194-212.]). The device can also be used for ambient-temperature data collection (Bowler, Mueller et al., 2015[Bowler, M. W., Mueller, U., Weiss, M. S., Sanchez-Weatherby, J., Sorensen, T. L. M., Thunnissen, M. M. G. M., Ursby, T., Gobbo, A., Russi, S., Bowler, M. G., Brockhauser, S., Svensson, O. & Cipriani, F. (2015). Cryst. Growth Des. 15, 1043-1054.]; Russi et al., 2011[Russi, S., Juers, D. H., Sanchez-Weatherby, J., Pellegrini, E., Mossou, E., Forsyth, V. T., Huet, J., Gobbo, A., Felisaz, F., Moya, R., McSweeney, S. M., Cusack, S., Cipriani, F. & Bowler, M. W. (2011). J. Struct. Biol. 175, 236-243.]) where the RH must be matched to the mother liquor to prevent dehydration of the crystal. The first step in these experiments is to define the equilibrium point between the RH and the mother liquor of the sample. This is an essential step as it defines the starting point for the experiments and maintains the crystal in a stable environment when the mother liquor is removed. In order to facilitate this stage we measured the equilibrium RH points for a variety of solutions commonly used for the crystallization of proteins and nucleic acids (Wheeler et al., 2012[Wheeler, M. J., Russi, S., Bowler, M. G. & Bowler, M. W. (2012). Acta Cryst. F68, 111-114.]). This provided a starting point for most experiments and the results obtained were compared with Raoult's law (Raoult, 1887[Raoult, F.-M. (1887). C. R. Acad. Sci. Paris, 104, 1430-1433.]) for the equilibrium vapour pressure of water above a solution [and for solutions of polymers, with a generalization (Bowler, Mueller et al., 2015[Bowler, M. W., Mueller, U., Weiss, M. S., Sanchez-Weatherby, J., Sorensen, T. L. M., Thunnissen, M. M. G. M., Ursby, T., Gobbo, A., Russi, S., Bowler, M. G., Brockhauser, S., Svensson, O. & Cipriani, F. (2015). Cryst. Growth Des. 15, 1043-1054.])]. The measurements made were consistently higher than those predicted by Raoult's law and a satisfactory explanation for the discrepancy could not be found. Here, we have repeated the measurements using a device based on the HC1 but with higher precision in the control of RH. The new measurements are in very good agreement with Raoult's law. Because of its importance, we present a simple explanation for Raoult's law using statistical mechanics and also show how this treatment can be extended to polymer solutions, where Raoult's law breaks down. These results illuminate the machinery underlying a long-observed phenomenon and allow the accurate prediction of humid atmospheres for specific sample requirements, applicable to a wide variety of fields.

2. Experimental procedures

2.1. RH measurements

Solutions of polyethylene glycol (PEG) were made gravimetrically at concentrations between 50 and 10%(w/w). Stock solutions of salts at 3 M were made and then diluted to reach the desired concentration. A round 600 µm Micromount (MiTeGen, Ithica, New York, USA) was mounted on either the BM14 or MASSIF-1 (Bowler, Nurizzo et al., 2015[Bowler, M. W., Nurizzo, D. et al. (2015). J. Synchrotron Rad. 22, 1540-1547.]; Nurizzo et al., 2016[Nurizzo, D., Bowler, M. W., Caserotto, H., Dobias, F., Giraud, T., Surr, J., Guichard, N., Papp, G., Guijarro, M., Mueller-Dieckmann, C., Flot, D., McSweeney, S., Cipriani, F., Theveneau, P. & Leonard, G. A. (2016). Acta Cryst. D72, 966-975.]) diffractometers with an HC-Lab device (Arinax, Moirons, France) mounted at a distance of 5 mm from the loop. The HC-Lab is based on the original HC1 developed at the EMBL, Grenoble (Sanchez-Weatherby et al., 2009[Sanchez-Weatherby, J., Bowler, M. W., Huet, J., Gobbo, A., Felisaz, F., Lavault, B., Moya, R., Kadlec, J., Ravelli, R. B. G. & Cipriani, F. (2009). Acta Cryst. D65, 1237-1246.]), but with improvements in the dew point controller, temperature measurement and calculation of RH. These developments have led to a device with superior control and stability of RH levels. In order to determine the equilibrium RH, 2 µl of solution were taken and a small drop placed on the loop with a pipette. The diameter of the drop was measured using specific image analysis software. The humidity was adjusted until the drop diameter was stable. This was repeated a few times until the drop diameter was stable upon initial placement on the loop. Each measurement was then repeated three times at ambient temperature.

3. Results

3.1. Comparison of measured equilibria and predicted values

In previous work we measured the RH equilibrium points for a range of solutions commonly used in protein crystallization and examined the results in terms of Raoult's law and the Flory–Huggins model for the entropy of mixing of polymers (Bowler, Mueller et al., 2015[Bowler, M. W., Mueller, U., Weiss, M. S., Sanchez-Weatherby, J., Sorensen, T. L. M., Thunnissen, M. M. G. M., Ursby, T., Gobbo, A., Russi, S., Bowler, M. G., Brockhauser, S., Svensson, O. & Cipriani, F. (2015). Cryst. Growth Des. 15, 1043-1054.]; Wheeler et al., 2012[Wheeler, M. J., Russi, S., Bowler, M. G. & Bowler, M. W. (2012). Acta Cryst. F68, 111-114.]). While the measured values provided a starting point for humidity control experiments and Raoult's law should be a good explanation for the observed results, there was a considerable discrepancy between the two (Wheeler et al., 2012[Wheeler, M. J., Russi, S., Bowler, M. G. & Bowler, M. W. (2012). Acta Cryst. F68, 111-114.]). Measured values were consistently 1–3% higher than those predicted, which was attributed to the condenser used in the device being rather inaccurate at humidity values above 96%. Repeating these measurements using the new humidity control device, the HC-Lab, the discrepancy is no longer significant (Figs. 1[link]a and 2[link]a). The results obtained from the HC-Lab are also in agreement with detailed studies of the activity of water above salt (Robinson, 1945[Robinson, R. A. (1945). Trans. Proc. R. Soc. NZ, 75, 203-217.]; Wishaw & Stokes, 1954[Wishaw, B. F. & Stokes, R. H. (1954). Trans. Faraday Soc. 50, 952-954.]) and polymer solutions (Sadeghi & Shahebrahimi, 2011[Sadeghi, R. & Shahebrahimi, Y. (2011). J. Chem. Eng. Data, 56, 789-799.]; Sadeghi & Ziamajidi, 2006[Sadeghi, R. & Ziamajidi, F. (2006). J. Chem. Thermodyn. 38, 1335-1343.]) (Figs. 1[link]b and 2[link]b), with the salt solution measurements made in this study appearing to be more accurate. This now brings the control of RH surrounding crystals into line with measurements made using dedicated and accurate devices, as well as with theoretical calculations.

[Figure 1]
Figure 1
(a) Plots showing the equilibrium RH for salt solutions commonly used as precipitants or additives in macromolecular crystallogenesis measured using the HC-Lab. (b) The measured vapour pressures above solutions of ammonium sulfate (Wishaw & Stokes, 1954[Wishaw, B. F. & Stokes, R. H. (1954). Trans. Faraday Soc. 50, 952-954.]) and sodium chloride (Robinson, 1945[Robinson, R. A. (1945). Trans. Proc. R. Soc. NZ, 75, 203-217.]). The lines show the calculated RH from Raoult's law (Wheeler et al., 2012[Wheeler, M. J., Russi, S., Bowler, M. G. & Bowler, M. W. (2012). Acta Cryst. F68, 111-114.]). The measurements made using the HC-Lab [panel (a)] more accurately reflect the predicted values from Raoult's law.
[Figure 2]
Figure 2
(a) Plots showing the equilibrium RH for PEG solutions commonly used as precipitants or additives in macromolecular crystallogenesis measured using the HC-Lab. (b) The measured vapour pressures above PEG solutions from Sadeghi & Shahebrahimi (2011[Sadeghi, R. & Shahebrahimi, Y. (2011). J. Chem. Eng. Data, 56, 789-799.]) and Sadeghi & Ziamajidi (2006[Sadeghi, R. & Ziamajidi, F. (2006). J. Chem. Thermodyn. 38, 1335-1343.]). The lines show the calculated RH from Raoult's law modified for polymer solutions (Bowler, Mueller et al., 2015[Bowler, M. W., Mueller, U., Weiss, M. S., Sanchez-Weatherby, J., Sorensen, T. L. M., Thunnissen, M. M. G. M., Ursby, T., Gobbo, A., Russi, S., Bowler, M. G., Brockhauser, S., Svensson, O. & Cipriani, F. (2015). Cryst. Growth Des. 15, 1043-1054.]).

3.2. Derivation of the origin of Raoult's law

Raoult's law (Raoult, 1887[Raoult, F.-M. (1887). C. R. Acad. Sci. Paris, 104, 1430-1433.]) describes the reduction in the saturated vapour pressure above a solvent when a mole fraction x of some solute is dissolved within it. If the vapour pressure above the pure solvent is p0 then the vapour pressure of the solvent above the solution is given by

[p = p_0 (1 - x) . \eqno(1)]

This is of course an idealization, but it is remarkably good, particularly at low mole fractions of the solute. Originally empirical, from what principles can it be derived? Any such derivations depend on the assumption of an ideal solution, meaning that within the body of the solution the elements of the solute are nearly identical to the elements of the solvent (and yet for a non-volatile solute the solute cannot enter the vapour phase). In thermodynamics, equilibrium at constant temperature and pressure corresponds to a minimum of the Gibbs' function G and hence liquid–vapour equilibrium requires equal chemical potentials. The chemical potential of the solvent vapour phase is the same as that of the solvent, both above the pure liquid solvent and above a solution. The chemical potential in the solution is reduced by mixing; thermo­dynamic arguments are used to turn an entropy of mixing into a change in chemical potential. Thermodynamics does not deal with the mechanisms underlying these steps and it seems reasonable to ask, first, how the vapour pressure can be affected by the number of ways of arranging fixed numbers of two kinds of molecule and, secondly, why is there no apparent role for a work function related to the latent heat of vaporization?

Raoult's law is the direct result of the dilution of the solvent by the solute and can be extracted by applying elementary statistical mechanics. The machinery involves the energy levels the confined components can occupy and, in the simplest case of non-ideal solutions, differences in work functions are both important and easily calculated.

3.2.1. Statistical mechanics

It is a truth universally acknowledged that any system (such as an atom in a box) that has energy levels i and is in thermal equilibrium at temperature T has a probability of occupying a given level proportional to [\exp(-\varepsilon_i / k_{\rm B} T)], where kB is Boltzmann's constant, for in an ensemble the vast majority of possible configurations have this distribution and for macroscopic phenomena we are concerned with sums or averages over very many individual microscopic systems (here atoms, ions or molecules). For pure solvent we divide the energy levels into two classes, those in the liquid and those in the vapour phases. They are separated by a step in energy, a work function W, and so the number niv, from a total of N atoms, found in the ith vapour state of energy [\varepsilon_i^{\rm v}] + W is given by

[\eqalignno{n_i^{\rm v} = & \, N \exp\left (-{{\varepsilon_i^{\rm v} + W} \over {k_{\rm B} T}} \right) \bigg / \Biggl [\sum \limits_{j} \exp \left (-{{\varepsilon_j^{\rm v} + W} \over {k_{\rm B} T}} \right) \cr & \, + \sum \limits_{k} \exp \left (-{{\varepsilon_k^{\rm l}} \over {k_{\rm B} T}} \right) \Biggr] . &(2)}]

Here, the factor following the total number N is the probability of finding a solvent molecule in a vapour state of energy [\varepsilon_i^{\rm v}] above energy W, and [\varepsilon_k^{\rm l}] is the energy of the liquid state k. The sum over the index j in equation (2)[link] is over the vapour states and the index k over the liquid states.

For a given temperature, the total number of atoms in the vapour is found by summing the numerator of equation (2)[link] over the index i, yielding a fraction y of the total number N. The vapour energy levels start raised above the energy levels in the liquid by the work function W (closely related to the latent heat) and so the fraction of atoms in the vapour contains a suppression factor of exp(- W / kBT). We are not yet concerned with this factor, nor with the details of the structure of the energy levels. It suffices that, for a given temperature and container, the number of atoms in the vapour phase is the fraction y of the total number of solvent atoms N. This fraction is determined by the work function, the temperature and the detailed structure of the energy levels, in turn determined by the volumes available. If a fraction x of the solvent atoms are removed and replaced by Nx units of solute, changing nothing else, the volume of the container does not change and neither the detailed structure of the energy levels nor the work function for solvent atoms changes because of the (close) identity of the solvent and solute units in an ideal solution. The fraction of solvent atoms in the vapour phase does not change and, because there are now only (1 − x)N atoms of solvent, the number of atoms of solvent in the vapour phase is reduced by a factor (1 − x). Hence the reduced vapour pressure and Raoult's law.

This simple argument is indubitably correct, given the assumptions of an ideal solution. The flux of solvent molecules leaving the surface is reduced by a factor (1 − x), and for equilibrium both the returning flux and the number density of solvent molecules in the vapour phase are also reduced by a factor (1 − x), as a direct result of the lower concentration of solvent molecules. This approach can be extended to non-ideal solutions (such as solutions of polymers), but this is more complicated because of the need to calculate differences in work functions.

3.2.2. Some technical details concerning volume

A second result from elementary statistical mechanics removes a potential objection to the above argument. What if the volume of pure solvent is reduced? If the volumes of liquid and vapour are held constant, the number of vapour atoms is (for a fixed temperature) a definite fraction of the number of atoms in the liquid phase. The more general result is that the concentration of atoms in the vapour phase is a definite fraction of the concentration of atoms in the liquid phase. The vapour pressure above a liquid in a sealed container does not, in equilibrium, depend on the volume of liquid in the container. Thus (1 − x)N atoms of solvent in the container without xN atoms of dissolved solute would not (and does not) result in a pressure reduced by (1 − x). The reason is as follows. The energy levels for atoms in the vapour are those of particle waves confined within the volume between the liquid surface and the walls of the container. For an ideal gas, the number of energy levels in a given interval of energy is proportional to the volume – the spacing goes down as the volume goes up. If the volume available to vapour doubles, the number of levels in some interval Δ at also doubles and hence so does the number of molecules in the vapour. Thus the concentration of atoms in the vapour phase is constant as the volume increases – the pressure remains the same. Similarly, the molecules in the liquid roam throughout the liquid volume and their wavefunctions are constrained by the walls and the liquid surface. If the volume of liquid is reduced, the sum over the populations of liquid energy levels is reduced because there are fewer of them. The spacing between energy levels in the liquid goes up with the reduction in volume and the concentration in the liquid remains the same. Thus the saturated vapour pressure above the liquid remains constant as the ratio of vapour volume to liquid volume is increased, until of course all the atoms originally in the liquid are in the vapour phase. Thereafter, as the volume is increased (by pulling back on a piston perhaps) the vapour density, and so the pressure along the isotherm, falls.

When extracted solvent molecules are replaced by solute, the solute molecules make up the missing liquid volume. This makes available to the reduced number of solvent molecules the same energy level structures in both the liquid and vapour phases. This dependence of the energy-level density on the free-range volume results in the concentration of atoms in the vapour phase being a definite fraction of the concentration in the solution. This is important for considering the vapour pressure above solutions that are not ideal, for example polymers. Finally, it is essential for understanding the thermo­dynamic treatment and entropy of mixing.

3.2.3. Solutions of molecules of different sizes

Suppose now that, instead of replacing a fraction of molecules of solvent with molecules of solute pre-empting the same volume, the solute molecules require a different volume. For the case of polymers, such as polyethyl­ene glycol (PEG), the specific volume will be larger, very substantially larger for the heavier long-chain polymers. Let there be N1 molecules of solvent of specific volume v1; similarly for the solute N2, v2. The volume occupied by the liquid solution is N1v1 + N2v2 and the concentration of solvent molecules is less than for pure solvent occupying the same volume. The ratio of concentrations of the solvent molecules in the solution to pure solvent gives a factor in the vapour pressure ratio of

[{{N_1 v_1} \over {N_1 v_1 + N_2 v_2}}. \eqno(3)]

This factor reduces to Raoult's law as the specific volumes of solvent v1 and solute v2 approach equality. This is not the whole story because the work needed to remove a solvent molecule from solution is not equal to that required to remove a solvent molecule from pure solvent, except in this limit. The following simple calculation yields the requisite difference in work functions. The work function is the work that has to be done when removing a molecule against the cohesive forces in the liquid and any contribution from ambient pressure. Because the forces are cohesive, the removal of a volume ΔV of liquid to the vapour state requires energy −PcΔV, where the quantity Pc is the potential energy density associated with the cohesive forces. It is a contribution to the pressure in the liquid and is negative. If a volume ΔV is instead added, it acquires negative potential energy and the work done is PcΔV, where Pc is again negative. Consider the operation of replacing a molecule of solvent by one of solute. The liquid volume increases by ΔV = (v2v1) and this volume contains a negative potential energy density. The cohesive pressure term −Pc must balance that from the thermal energy density (both are of the order of 1000 atmospheres and ambient pressures permitting the liquid state are perhaps 1 atmosphere) and so is given by

[-P_{\rm c} \left (N_1 v_1 + N_2 v_2 \right) = \left (N_1 + N_2 \right) \, k_{\rm B} T . \eqno(4)]

Thus the work that has to be done to make the replacement is given by

[- {{k_{\rm B} T \,\left (N_1 + N_2 \right) \, \left (v_2 - v_1 \right)} \over {N_1 v_1 + N_2 v_2}} . \eqno(5)]

This is made up of two parts, the work necessary to insert a molecule of solute (a contribution to the chemical potential μ2) and the work necessary to extract a molecule of solvent (−μ1). The question is, what part of equation (5)[link] is to be identified with the component of μ1 and, because equation (5)[link] is a difference, what is the origin? As N2 [\to] 0, the solution approaches pure solvent and the term in the numerator involving N2 goes to zero, thus suggesting that the difference in the work that has to be done to deliver one molecule of solvent to the solution as opposed to pure solvent is

[\Delta W^{\uparrow} = {{k_{\rm B} T N_2 \left (v_2 - v_1 \right)} \over {N_1 v_1 + N_2 v_2}} . \eqno(6)]

This can be verified by calculating the work done against cohesive pressure to insert a solvent molecule into the solution as opposed to the same volume of pure solvent: calculate the (pressure-related) work done inserting an atom of solvent 1 into a solution and also calculate the work done inserting an additional atom into a volume of pure species 1. In both cases Δv = v1.

The magnitude of the cohesive pressure in a solution is given by

[{{\left (N_1 + N_2 \right) \, k_{\rm B} T} \over {N_1 v_1 + N_2 v_2}} \eqno(7)]

[from equation (4)[link]] and the pressure in a pure solvent is kBT/v1. Then

[- \Delta (p \Delta v) = - \left \{ {{\left (N_1 + N_2 \right) v_1} \over {N_1 v_1 + N_2 v_2}} - 1 \right \} k_{\rm B} T . \eqno(8)]

This also yields equation (6)[link].

The difference in work functions for removing atoms to the vapour phase, [\Delta {W}^{\downarrow}], is the negative of equation (6)[link]. The effect on RH is an exponential factor

[\exp \left ( -{{\Delta W^{\downarrow}} \over {k_{\rm B} T}} \right ) = \exp \left ( {{\Delta W^{\uparrow}} \over {k_{\rm B} T}} \right ) . \eqno(9)]

The concentration ratio of equation (3)[link] multiplied by this factor yields the RH of the solvent:

[{p \over {p_0}} = {{N_1 v_1} \over {N_1 v_1 + N_2 v_2}} \exp \left [ {{N_2 \left ( v_2 - v_1 \right )} \over {N_1 v_1 + N_2 v_2}} \right ] . \eqno(10)]

The first factor on the right-hand side is the volume fraction of solvent in the solution and reduces to Raoult's law as the specific volumes become equal. The second factor goes to unity in this same limit. It is less obvious that equation (10)[link] also reduces to Raoult's law in the limit of extreme dilution, regardless of the ratio of specific volumes, but it is so.

This expression [equation (10)[link]], derived using elementary notions from statistical mechanics, is the same as that derived using thermodynamics and the Flory–Huggins entropy of mixing devised for polymer solutions (Flory, 1942[Flory, P. J. (1942). J. Chem. Phys. 10, 51-61.], 1970[Flory, P. J. (1970). Discuss. Faraday Soc. 49, 7-29.]) or, equivalently, Hildebrand's entropy of solution of molecules of different sizes (Hildebrand, 1947[Hildebrand, J. H. (1947). J. Chem. Phys. 15, 225-228.]). In such treatments both factors in equation (10)[link] emerge from matching chemical potentials. Our treatment clarifies the physical meaning of the factors – the first factor is the concentration ratio, while the second (exponential) factor embodies the difference in work functions arising from different specific volumes. In Appendix A[link] we discuss the relationship between simple statistical mechanics and thermodynamic arguments, addressing in particular the significance of the entropy of mixing.

3.2.4. Relationship between observations and theory

We have shown that there is good agreement between measured values of the RH above a solution and the theoretical basis for vapour pressure above a solution. How do the curves shown in Figs. 1[link] and 2[link] relate to equations (1)[link] and (10)[link]? The RH given by Raoult's law is a linear function of the fraction of molecules that are solvent, the mole fraction, but the curves shown in the figures are not linear for two reasons: (i) the concentrations are shown in units that are commonly used in biochemistry, not the mole fraction, and (ii) the specific volume of the solute changes the relationship. This section explains the observed curves, first for salts and then for polymers.

Writing Raoult's law as

[{p \over {p_0}} = 1 - {{N_2} \over {N_1 + N_2}}, \eqno(11)]

it is obvious that the RH is a linear function of the fraction of more or less freely propagating components of species 2. However, for practical reasons solutions are not usually prepared as a mole fraction, but rather of a specified molarity, the number of moles of solute in a litre of solution (the graphs in Fig. 1[link] are plotted as a function of molarity). If the solute molarity is to be specified, there are two complications in the translation of Raoult's law. The first is that (ionic) salts when dissolved dissociate into freely drifting ions (such as Na+ and Cl). The second is that the volume of water is reduced below 1 l by the volume of the salt. Thus, if M is the solute molarity, the RH is given by

[{p \over {p_0}} = 1 - {{xM} \over {xM + {{1000} \over {18}} (1 - yM)}}, \eqno(12)]

where x is the number of independent ions into which the salt dissociates and y accounts for the specific volume of the salt. For sodium chloride the quantities x (y) are 2 (0.027), for ammonium sulfate 2 (0.074) and for sodium malonate 3 (0.095). The above equation is in fact equation (3) of Wheeler et al. (2012[Wheeler, M. J., Russi, S., Bowler, M. G. & Bowler, M. W. (2012). Acta Cryst. F68, 111-114.]). It is clear that the Raoult's law RH is not a linear function of molarity and also that the slope at low molarity depends directly on the degree of dissociation of the salt, clearly seen in Fig. 1[link](a).

Comparison of equations (1)[link] and (10)[link] makes it clear that if the specific volumes are not the same, and the solution is not very dilute, the RH is not a linear function of mole fraction. There is a further complication: equation (10)[link] is primarily used for solutions of polymers where mass fraction (w/w), rather than mole fraction, is the most commonly used expression of concentration. Thus, to obtain the RH in the form quoted as equation (1) of Bowler, Mueller et al. (2015[Bowler, M. W., Mueller, U., Weiss, M. S., Sanchez-Weatherby, J., Sorensen, T. L. M., Thunnissen, M. M. G. M., Ursby, T., Gobbo, A., Russi, S., Bowler, M. G., Brockhauser, S., Svensson, O. & Cipriani, F. (2015). Cryst. Growth Des. 15, 1043-1054.]) the following steps are taken. Equation (10)[link] is written as

[{p \over {p_0}} = {1 \over {1 + f}} \exp \left [\left ({1 - {1 \over r}} \right) {f \over {1 + f}} \right] , \eqno(13)]

where f = rN2/N1 and r = v2/v1. The mole ratio is written in terms of the ratio of polymer mass to solution mass, here called x but in Fig. 2[link], and in other usage, weight per weight (w/w):

[{{N_2} \over {N_1}} = {x \over {1 - x}} \, {{18} \over n} , \eqno(14)]

where n is the molecular weight of the polymer (as in PEG n) and 18 that of water, assumed to be the solvent.

The remaining problem is the value to be adopted for the volume ratio r. The expression for the RH above a polymer solution was originally worked out using the Flory–Huggins entropy of mixing. These early calculations supposed a lattice, with water molecules each occupying one site and each monomer unit of a polymer likewise occupying one site. Then, if the molecular weight of the monomer is m the quantity r is n/m and

[f = {x \over {1 - x}} \, {{18} \over m}. \eqno(15)]

Substitution yields equation (1) of Bowler, Mueller et al. (2015[Bowler, M. W., Mueller, U., Weiss, M. S., Sanchez-Weatherby, J., Sorensen, T. L. M., Thunnissen, M. M. G. M., Ursby, T., Gobbo, A., Russi, S., Bowler, M. G., Brockhauser, S., Svensson, O. & Cipriani, F. (2015). Cryst. Growth Des. 15, 1043-1054.]), used in the construction of the curves in Fig. 2[link]. (We have found that the best value of m for PEGs is 38.) Expressed as a function of mass fraction, the RH becomes independent of the polymer molecular mass n as n becomes very large, i.e. for very long chain polymers.

4. Discussion

The control of the RH surrounding samples is important to maintain their integrity and study the effects of increased or decreased humidity. Here we have established that the theoretical RH values we previously calculated (Bowler, Mueller et al., 2015[Bowler, M. W., Mueller, U., Weiss, M. S., Sanchez-Weatherby, J., Sorensen, T. L. M., Thunnissen, M. M. G. M., Ursby, T., Gobbo, A., Russi, S., Bowler, M. G., Brockhauser, S., Svensson, O. & Cipriani, F. (2015). Cryst. Growth Des. 15, 1043-1054.]; Wheeler et al., 2012[Wheeler, M. J., Russi, S., Bowler, M. G. & Bowler, M. W. (2012). Acta Cryst. F68, 111-114.]) are in satisfactory agreement with a humidity control device used on protein crystallography beamlines. As the predicted values are also in complete agreement with measurements made using specific devices, the previous discrepancies can be ascribed to shortcomings in the control of RH in the HC1c device used. We have also determined the origin of the observed vapour pressure changes above solutions of solutes. If N units of a liquid solvent are in an equilibrium where liquid and vapour phases coexist, a fixed fraction are (for a given temperature) in the vapour phase. If the number of units is reduced to N(1 − x) and if all else remains unchanged, because of the presence of Nx units of the solute in an ideal solution, then the number of units in the vapour phase (and hence the pressure) is reduced by the same factor (1 − x), Raoult's law. For unequal sizes of solvent and solute components, the dilution factor has to be multiplied by an exponentiated work function. These results provide a solid basis on which to predict the RHs required to maintain a wide variety of samples and solutions in homeostasis.

APPENDIX A

Statistical mechanics, entropy and chemical potentials

The origin of Raoult's law lies in the freedom of the units of solvent and solute to roam throughout the volume of liquid. For the assumptions of an ideal solution, access of both solvent and solute to the whole volume results in energy levels available (to the solvent) for a given volume, unchanged from those in the pure solvent, and the density of (energy) states is proportional to the volume. Entropy of mixing expresses these same ideas in the language of thermodynamics.

Suppose that we can decompose a system into many identical parts having energy levels i. This complex system is in thermal equilibrium at some temperature T, with an exponential distribution in the energies of the components. Let the total energy of our complex system be U. The following relations apply:

[n_i = {{N \exp \left (-{{\varepsilon_i} / {k_{\rm B} T}} \right)} \over {z}} \eqno(16)]

and

[U = \sum \limits_{i} n_i \varepsilon_i = N {{\sum _i \varepsilon_i \exp \left (- {{\varepsilon_i} / {k_{\rm B} T}} \right)} \over {z}}, \eqno(17)]

where

[z = \textstyle\sum \limits_i \exp \left (- \varepsilon_i/k_{\rm B} T \right). \eqno(18)]

A small change in the internal energy U can be written as

[\Delta U = \textstyle\sum \limits_i \varepsilon_i \Delta n_i + \sum \limits_i n_i \Delta \varepsilon_i + \mu \sum \limits_i \Delta n_i. \eqno(19)]

The last term vanishes if there is no change in the number of components in the system. The first term is the result of slightly redistributing the population over the energy levels i. It represents the addition of heat. The second term corresponds to the energy levels changing with no change in population – doing it very slowly. If the volume slowly increases the energy levels slowly sag as the wavelengths of standing waves increase and the system does work. Thus the equivalent expression in thermodynamics is the first law in the form

[\Delta U = \Delta Q + \Delta W + \mu \Delta N, \eqno(20)]

where the last term is called chemical work and μ is the chemical potential. (Generally, each species of atom has its own chemical potential.)

If everything is done very slowly and reversibly,

[\Delta Q = T \Delta S, \eqno(21)]

where S is the thermodynamic entropy going back to Carnot. We are identifying the heat term in the first law with

[\Delta Q = \textstyle\sum \limits_i \varepsilon_i \Delta n_i . \eqno(22)]

Express the energy of the ith level in terms of its population

[\varepsilon_i = -k_{\rm B} T \ln \left (P_i z \right), \quad P_i = n_i/N . \eqno(23)]

For fixed N the sum of Δni is zero and so the term in ln(z) drops out and

[\Delta Q = - k_{\rm B} TN \Delta \left ( \textstyle\sum \limits_i P_i \ln P_i \right ) = T \Delta S, \eqno(24)]

demonstrating the equivalence of the Carnot and Boltzmann entropies.

Thus the entropy associated with N units (atoms, molecules, ions…) distributed over these energy levels, in equilibrium at temperature T, is

[S = - N \, k_{\rm B} \textstyle\sum \limits_i P_i \ln P_i . \eqno(25)]

The probabilities Pi involve the normalizing factor z, a sum over all energy levels. The value depends on the level density. The more levels the components are spread over, the smaller the individual Pi and the larger the entropy.

Substitution of the expressions for the exponential probabilities yields for the entropy

[S = {U \over T} + N \, k_{\rm B} \ln z . \eqno(26)]

Rewriting the sum for z as an integral

[z = \int \exp \left ({{- \varepsilon} \over {k_{\rm B} T}} \right) {{{\rm d}n} \over {{\rm d}\varepsilon}} \, {\rm d} \varepsilon. \eqno(27)]

The density of states factor dn/d depends on (which integrates out) and is linearly proportional to the volume available to the wandering molecules. Consider taking a volume V1 of solvent and a volume V2 of solute. (This is most easily envisaged if the solute is also a liquid; otherwise, pretend that there is a solute liquid with the properties the solute will display in the solution.) Before mixing the two together each has its entropy, appropriate to volumes V1 and V2, respectively. After mixing, both solvent and solute have access to a total volume V1 + V2. For an ideal solution nothing else has changed and so, taking the difference in entropy after and before the mixing,

[\eqalignno{S_{12} - S_1 - S_2 = & \, k_{\rm B} [ \left ( N_1 + N_2 \right ) \ln \left ( V_1 + V_2 \right ) \cr & \, - N_1 \ln V_1 - N_2 \ln V_2 ] . &(28)}]

This is the entropy of mixing and arises entirely from the increased density of energy levels as more volume is made available for units of both solvent and solute to roam at random. (These volumes are defined by the boundaries confining the liquids, setting boundary conditions and hence determining the quantized energy levels.)

Since we are looking at mixing of two forms of condensed matter, each with the same specific volume (an ideal solution again), the above expression for the entropy of mixing can also be written with V1 and V2 replaced by N1 and N2, respectively. The result is essentially identical to the product of kB (Boltzmann's constant) and the logarithm of the number of different ways of arranging N1 and N2 units (for large N1 and N2, using Stirling's theorem). This is a purely combinatorial problem and the number of perceptibly different ways is given by

[{{\left ( N_1 + N_2 \right )!} \over {N_1 ! \, N_2!}} . \eqno(29)]

What can this have to do with the vapour pressure above a solution? We now see that N1 and N2 are (for an ideal solution) proxies for V1 and V2 and these volumes control the energy levels available to the components of the solution before and after mixing.

More generally, suppose that the solvent molecules are each associated with a free volume v1 and the solute molecules with v2. Then the entropy of mixing [equation (28)[link] above] is

[k_{\rm B} \left [ \left ( N_1 + N_2 \right ) \ln \left ( N_1 v_1 + N_2 v_2 \right ) - N_1 \ln \left ( N_1 v_1 \right ) - N_2 \ln \left ( N_2 v_2 \right ) \right ] . \eqno(30)]

This is essentially the expression for the entropy of mixing for solvent and solute molecules of different free volumes to be found in equation (3) of Hildebrand (1947[Hildebrand, J. H. (1947). J. Chem. Phys. 15, 225-228.]), where the volumes are introduced through a classical argument concerning uncertainty of location. It is also equivalent to the Flory–Huggins entropy for polymer solutions, most clearly discussed by Flory (1970[Flory, P. J. (1970). Discuss. Faraday Soc. 49, 7-29.]).

The derivative of the entropy of mixing with respect to the number of solvent molecules within the solution (N1) yields the difference in chemical potentials that must match the difference in chemical potential of the vapours above the solution and the pure solvent:

[\Delta \mu = - T \, {{\partial \Delta S} \over {\partial N_1}} . \eqno(31)]

For the solution, the derivative of S12 in equation (30)[link] with respect to N1 is

[\mu _{12}^1 = - k_{\rm B} T \, {\partial \over {\partial N_1}} \left ( N_1 + N_2 \right ) \ln \left ( N_1 v_1 + N_2 v_2 \right ) , \eqno(32)]

and for the pure solvent before mixing the derivative of S1 is

[\mu _1^1 = - k_{\rm B} T \, {\partial \over {\partial N_1}} N_1 \ln \left ( N_1 v_1 \right ) . \eqno(33)]

In the standard thermodynamic argument (e.g. Hildebrand, 1947[Hildebrand, J. H. (1947). J. Chem. Phys. 15, 225-228.]), taking the difference in chemical potentials and matching to the vapour phase eventually yields

[\ln \left ( {p \over {p_0}} \right ) = \ln \left [ {{N_1 v_1} \over {N_1 v_1 + N_2 v_2}} \right]+ {{N_2 \left ( v_2 - v_1 \right )} \over {N_1 v_1 + N_2 v_2}} , \eqno(34)]

which is equation (6) given by Hildebrand (1947[Hildebrand, J. H. (1947). J. Chem. Phys. 15, 225-228.]). Then the RH of the solvent above such a solution is

[{p \over {p_0}} = {{N_1 v_1} \over {N_1 v_1 + N_2 v_2}} \exp \left [ {{N_2 \left ( v_2 - v_1 \right )} \over {N_1 v_1 + N_2 v_2}} \right ] . \eqno(35)]

Equation (35)[link] is identical to equation (10)[link].

In §3.2.3[link] we calculated the difference in work functions for the solvent in a solution of volume V and for the pure solvent in the same volume. This result can also be obtained from the differential of the difference in entropies of the solution and the pure solvent in equal volumes. The only terms that survive in the difference are (N1 + N2)kBlnV for the solution and N10 kB lnV for the pure solvent. The volume V is given by

[V = N_1 v_1 + N_2 v_2 = N_1^0 v_1 . \eqno(36)]

The relevant difference in chemical potentials is then

[k_{\rm B} T \, {{N_2 \left ( v_2 - v_1 \right )} \over {N_1 v_1 + N_2 v_2}} . \eqno(37)]

The negative of this is the difference in work functions, needed to complete the ratio of vapour pressures at the end of §3.2.3[link]. The result given in equation (37)[link] above of course agrees with equation (6)[link].

Acknowledgements

MWB thanks Babu Manjasetty and Hassan Belrhali (EMBL, Grenoble) for access to the HC-Lab on beamline BM14, and Ralf Siebrecht and Jaouhar Nasri (Arinax) for stimulating discussions.

References

First citationAdachi, H., Umena, Y., Enami, I., Henmi, T., Kamiya, N. & Shen, J.-R. (2009). BBA Bioenergetics, 1787, 121–128.  CrossRef PubMed CAS Google Scholar
First citationBerthou, J., Cesbron, F. & Laurent, A. (1972). J. Mol. Biol. 71, 809–813.  CrossRef CAS PubMed Web of Science Google Scholar
First citationBowler, M. W., Montgomery, M. G., Leslie, A. G. W. & Walker, J. E. (2006). Acta Cryst. D62, 991–995.  Web of Science CrossRef CAS IUCr Journals Google Scholar
First citationBowler, M. W., Mueller, U., Weiss, M. S., Sanchez-Weatherby, J., Sorensen, T. L. M., Thunnissen, M. M. G. M., Ursby, T., Gobbo, A., Russi, S., Bowler, M. G., Brockhauser, S., Svensson, O. & Cipriani, F. (2015). Cryst. Growth Des. 15, 1043–1054.  Web of Science CrossRef CAS Google Scholar
First citationBowler, M. W., Nurizzo, D. et al. (2015). J. Synchrotron Rad. 22, 1540–1547.  CrossRef IUCr Journals Google Scholar
First citationCramer, P., Bushnell, D. A., Fu, J. H., Gnatt, A. L., Maier-Davis, B., Thompson, N. E., Burgess, R. R., Edwards, A. M., David, P. R. & Kornberg, R. D. (2000). Science, 288, 640–649.  CrossRef PubMed CAS Google Scholar
First citationEinstein, J. R. (1961). J. Sci. Instrum. 38, 449–451.  CrossRef CAS Web of Science Google Scholar
First citationEinstein, J. R. & Low, B. W. (1962). Acta Cryst. 15, 32–34.  CrossRef IUCr Journals Web of Science Google Scholar
First citationFlory, P. J. (1942). J. Chem. Phys. 10, 51–61.  CrossRef CAS Google Scholar
First citationFlory, P. J. (1970). Discuss. Faraday Soc. 49, 7–29.  CrossRef Google Scholar
First citationFratini, A. V., Kopka, M. L., Drew, H. R. & Dickerson, R. E. (1982). J. Biol. Chem. 257, 4686–4707.  Google Scholar
First citationGupta, V., Gupta, R. K., Khare, G., Salunke, D. M., Surolia, A. & Tyagi, A. K. (2010). PLoS One, 5, e9222.  Web of Science CrossRef PubMed Google Scholar
First citationHeras, B., Edeling, M. A., Byriel, K. A., Jones, A., Raina, S. & Martin, J. L. (2003). Structure, 11, 139–145.  Web of Science CrossRef PubMed CAS Google Scholar
First citationHildebrand, J. H. (1947). J. Chem. Phys. 15, 225–228.  CrossRef CAS Google Scholar
First citationHu, N. J., Iwata, S., Cameron, A. D. & Drew, D. (2011). Nature, 478, 408–411.  CrossRef CAS PubMed Google Scholar
First citationHuxley, H. E. & Kendrew, J. C. (1953). Acta Cryst. 6, 76–80.  CrossRef IUCr Journals Web of Science Google Scholar
First citationKadlec, J., Hallacli, E., Lipp, M., Holz, H., Sanchez-Weatherby, J., Cusack, S. & Akhtar, A. (2011). Nat. Struct. Mol. Biol. 18, 142–149.  CrossRef CAS PubMed Google Scholar
First citationKiefersauer, R., Than, M. E., Dobbek, H., Gremer, L., Melero, M., Strobl, S., Dias, J. M., Soulimane, T. & Huber, R. (2000). J. Appl. Cryst. 33, 1223–1230.  Web of Science CrossRef CAS IUCr Journals Google Scholar
First citationKuo, A., Bowler, M. W., Zimmer, J., Antcliff, J. F. & Doyle, D. A. (2003). J. Struct. Biol. 141, 97–102.  Web of Science CrossRef PubMed CAS Google Scholar
First citationLin, J.-M., Lin, T.-L. & Jeng, U.-S. (2007). J. Appl. Cryst. 40, s367–s372.  CrossRef CAS IUCr Journals Google Scholar
First citationMalinauskaite, L., Quick, M., Reinhard, L., Lyons, J. A., Yano, H., Javitch, J. A. & Nissen, P. (2014). Nat. Struct. Mol. Biol. 21, 1006–1012.  CrossRef CAS PubMed Google Scholar
First citationMcDonald, M., Kendall, A., Tanaka, M., Weissman, J. S. & Stubbs, G. (2008). J. Appl. Cryst. 41, 206–209.  Web of Science CrossRef CAS IUCr Journals Google Scholar
First citationMo, F. & Ramsøskar, K. (2009). J. Appl. Cryst. 42, 531–534.  Web of Science CrossRef CAS IUCr Journals Google Scholar
First citationNakamura, A., Wada, C. & Miki, K. (2007). Acta Cryst. F63, 346–349.  Web of Science CrossRef IUCr Journals Google Scholar
First citationNurizzo, D., Bowler, M. W., Caserotto, H., Dobias, F., Giraud, T., Surr, J., Guichard, N., Papp, G., Guijarro, M., Mueller-Dieckmann, C., Flot, D., McSweeney, S., Cipriani, F., Theveneau, P. & Leonard, G. A. (2016). Acta Cryst. D72, 966–975.  Web of Science CrossRef IUCr Journals Google Scholar
First citationOliete, R., Pous, J., Rodríguez-Puente, S., Abad-Zapatero, C. & Guasch, A. (2013). Acta Cryst. D69, 194–212.  Web of Science CrossRef IUCr Journals Google Scholar
First citationPerutz, M. F. (1946). Trans. Faraday Soc. 42, B187.  CrossRef Google Scholar
First citationPickford, M. G., Garman, E. F., Jones, E. Y. & Stuart, D. I. (1993). J. Appl. Cryst. 26, 465–466.  CrossRef Web of Science IUCr Journals Google Scholar
First citationRaoult, F.-M. (1887). C. R. Acad. Sci. Paris, 104, 1430–1433.  Google Scholar
First citationRobinson, R. A. (1945). Trans. Proc. R. Soc. NZ, 75, 203–217.  CAS Google Scholar
First citationRoedig, P., Duman, R., Sanchez-Weatherby, J., Vartiainen, I., Burkhardt, A., Warmer, M., David, C., Wagner, A. & Meents, A. (2016). J. Appl. Cryst. 49, 968–975.  Web of Science CrossRef CAS IUCr Journals Google Scholar
First citationRussi, S., Juers, D. H., Sanchez-Weatherby, J., Pellegrini, E., Mossou, E., Forsyth, V. T., Huet, J., Gobbo, A., Felisaz, F., Moya, R., McSweeney, S. M., Cusack, S., Cipriani, F. & Bowler, M. W. (2011). J. Struct. Biol. 175, 236–243.  Web of Science CrossRef CAS PubMed Google Scholar
First citationSadeghi, R. & Shahebrahimi, Y. (2011). J. Chem. Eng. Data, 56, 789–799.  Web of Science CrossRef CAS Google Scholar
First citationSadeghi, R. & Ziamajidi, F. (2006). J. Chem. Thermodyn. 38, 1335–1343.  Web of Science CrossRef CAS Google Scholar
First citationSam, M. D., Abbani, M. A., Cascio, D., Johnson, R. C. & Clubb, R. T. (2006). Acta Cryst. F62, 825–828.  Web of Science CrossRef CAS IUCr Journals Google Scholar
First citationSanchez-Weatherby, J., Bowler, M. W., Huet, J., Gobbo, A., Felisaz, F., Lavault, B., Moya, R., Kadlec, J., Ravelli, R. B. G. & Cipriani, F. (2009). Acta Cryst. D65, 1237–1246.  Web of Science CrossRef CAS IUCr Journals Google Scholar
First citationSjögren, T., Carlsson, G., Larsson, G., Hajdu, A., Andersson, C., Pettersson, H. & Hajdu, J. (2002). J. Appl. Cryst. 35, 113–116.  Web of Science CrossRef IUCr Journals Google Scholar
First citationTakayama, Y. & Nakasako, M. (2012). Rev. Sci. Instrum. 83, 054301.  Web of Science CrossRef PubMed Google Scholar
First citationVijayalakshmi, L., Krishna, R., Sankaranarayanan, R. & Vijayan, M. (2008). Proteins Struct. Function Bioinform. 71, 241–249.  CrossRef CAS Google Scholar
First citationWheeler, M. J., Russi, S., Bowler, M. G. & Bowler, M. W. (2012). Acta Cryst. F68, 111–114.  Web of Science CrossRef CAS IUCr Journals Google Scholar
First citationWishaw, B. F. & Stokes, R. H. (1954). Trans. Faraday Soc. 50, 952–954.  CrossRef CAS Web of Science Google Scholar
First citationYap, T. L., Chen, Y. L., Xu, T., Wen, D., Vasudevan, S. G. & Lescar, J. (2007). Acta Cryst. F63, 78–83.  Web of Science CrossRef CAS IUCr Journals Google Scholar
First citationZerrad, L., Merli, A., Schröder, G. F., Varga, A., Gráczer, E., Pernot, P., Round, A., Vas, M. & Bowler, M. W. (2011). J. Biol. Chem. 286, 14040–14048.  CrossRef CAS PubMed Google Scholar

This is an open-access article distributed under the terms of the Creative Commons Attribution (CC-BY) Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original authors and source are cited.

Journal logoJOURNAL OF
APPLIED
CRYSTALLOGRAPHY
ISSN: 1600-5767
Follow J. Appl. Cryst.
Sign up for e-alerts
Follow J. Appl. Cryst. on Twitter
Follow us on facebook
Sign up for RSS feeds