Linear scaling computation of the Fock matrix. VII. Periodic density functional theory at the Gamma point

Linear scaling quantum chemical methods for density functional theory are extended to the condensed phase at the Gamma point. For the two-electron Coulomb matrix, this is achieved with a tree-code algorithm for fast Coulomb summation [M. Challacombe and E. Schwegler, J. Chem. Phys. 106, 5526 (1997)], together with multipole representation of the crystal field [M. Challacombe, C. White, and M. Head-Gordon, J. Chem. Phys. 107, 10131 (1997)]. A periodic version of the hierarchical cubature algorithm [M. Challacombe, J. Chem. Phys. 113, 10037 (2000)], which builds a telescoping adaptive grid for numerical integration of the exchange-correlation matrix, is shown to be efficient when the problem is posed as integration over the unit cell. Commonalities between the Coulomb and exchange-correlation algorithms are discussed, with an emphasis on achieving linear scaling through the use of modern data structures. With these developments, convergence of the Gamma-point supercell approximation to the k-space integration limit is demonstrated for MgO and NaCl. Linear scaling construction of the Fockian and control of error is demonstrated for RBLYP6-21G* diamond up to 512 atoms.     

Linear scaling computation of the Fock matrix. VIII. Periodic boundaries for exact exchange at the Gamma point

A translationally invariant formulation of the Hartree-Fock (HF) Gamma-point approximation is presented. This formulation is achieved through introduction of the minimum image convention (MIC) at the level of primitive two-electron integrals, and implemented in a periodic version of the ONX algorithm [E. Schwegler, M. Challacombe, and M. Head-Gordon, J. Chem. Phys. 106, 9708 (1997)] for linear scaling computation of the exchange matrix. Convergence of the HF-MIC Gamma-point model to the HF k-space limit is demonstrated for fully periodic magnesium oxide, ice, and diamond. Computation of the diamond lattice constant using the HF-MIC model together with the hybrid PBE0 density functional [C. Adamo, M. Cossi, and V. Barone, THEOCHEM 493, 145 (1999)] yields a0=3.569 A with the 6-21G* basis set and a 3x3x3 supercell. Linear scaling computation of the HF-MIC exchange matrix is demonstrated for diamond and ice in the condensed phase.     

Monte Carlo simulations of critical cluster sizes and nucleation rates of water

We have calculated the critical cluster sizes and homogeneous nucleation rates of water at temperatures and vapor densities corresponding to experiments by Wolk and Strey [J. Phys. Chem B 105, 11683 (2001)]. The calculations have been done with an expanded version of a Monte Carlo method originally developed by Vehkamaki and Ford [J. Chem. Phys. 112, 4193 (2000)]. Their method calculates the statistical growth and decay probabilities of molecular clusters. We have derived a connection between these probabilities and kinetic condensation and evaporation rates, and introduce a new way for the calculation of the work of formation of clusters. Three different interaction potential models of water have been used in the simulations. These include the unpolarizable SPC/E [J. Phys. Chem. 91, 6269 (1987)] and TIP4P [J. Chem. Phys. 79, 926 (1983)] models and a polarizable model by Guillot and Guissani [J. Chem. Phys. 114, 6720 (2001)]. We show that TIP4P produces critical cluster sizes and a temperature and vapor density dependence for the nucleation rate that agree well with the experimental data, although the magnitude of nucleation rate is constantly overestimated by a factor of 2 x 10(4). Guissani and Guillot's model is somewhat less successful, but both the TIP4P and Guillot and Guissani models are able to reproduce a much better experimental temperature dependency of the nucleation rate than the classical nucleation theory. Using SPC/E results in dramatically too small critical clusters and high nucleation rates. The water models give different average binding energies for clusters. We show that stronger binding between cluster molecules suppresses the decay probability of a cluster, while the growth probability is not affected. This explains the differences in results from different water models.     

Molecular nitrogen-N2 properties: the intermolecular potential and the equation of state

Quantum mechanical (QM) high precision calculations were used to determine N(2)-N(2) intermolecular interaction potential. Using QM numerical data the anisotropic potential energy surface was obtained for all orientations of the pair of the nitrogen molecules in the rotation invariant form. The new N(2)-N(2) potential is in reasonably good agreement with the scaled potential obtained by van der Avoird et al. using the results of Hartree-Fock calculations [J. Chem. Phys. 84, 1629 (1986)]. The molecular dynamics (MD) of the N(2) molecules has been used to determine nitrogen equation of state. The classical motion of N(2) molecules was integrated in rigid rotor approximation, i.e., it accounted only translational and rotational degrees of freedom. Fincham [Mol. Simul. 11, 79 (1993)] algorithm was shown to be superior in terms of precision and energy stability to other algorithms, including Singer [Mol. Phys. 33, 1757 (1977)], fifth order predictor-corrector, or Runge-Kutta, and was therefore used in the MD modeling of the nitrogen pressure [S. Krukowski and P. Strak, J. Chem. Phys. 124, 134501 (2006)]. Nitrogen equation of state at pressures up to 30 GPa (300 kbars) and temperatures from the room temperature to 2000 K was obtained using MD simulation results. Results of MD simulations are in very good agreement (the error below 1%) with the experimental data on nitrogen equation of state at pressures below 1 GPa (10 kbars) for temperatures below 1800 K [R. T. Jacobsen et al., J. Phys. Chem. Ref. Data 15, 735 (1986)]. For higher temperatures, the deviation is slightly larger, about 2.5% which still is a very good agreement. The slightly larger difference may be attributed to the vibrational motion not accounted explicitly by rigid rotor approximation, which may be especially important at high temperatures. These results allow to obtain reliable equation of state of nitrogen for pressures up to 30 GPa (300 kbars), i.e., close to molecular nitrogen stability limit, determined by Nellis et al. [Phys. Rev. Lett. 53, 1661 (1984)].     

Polarizable and flexible model for ethanol

A polarizable, flexible model for ethanol is obtained based on an extensive series of B3LYP/6-311++G(d,p) calculations and molecular dynamics simulations. The ethanol model includes electric-field dependence in both the atomic charges and the intramolecular degrees of freedom. Field-dependent intramolecular potentials have been attempted only once previously, for OH and HH stretches in water [P. Cicu et al., J. Chem. Phys. 112, 8267 (2000)]. The torsional potential involving the hydrogen-bonding hydrogen in ethanol is found to be particularly field sensitive. The methodology for developing field-dependent potentials can be readily generalized to other molecules and is discussed in detail. Molecular dynamics simulations of bulk ethanol are performed and the results are assessed based on comparisons with the self-diffusion coefficient [N. Karger et al., J. Chem. Phys. 93, 3437 (1990)], dielectric constant [J. T. Kindt and C. A. Schmuttenmaer, J. Phys. Chem. 100, 10373 (1996)], enthalpy of vaporization [R. C. Wilhoit and B. J. Zwolinski, J. Phys. Chem. Ref. Data, Suppl. 2, 2 (1973)], and experimental interatomic distributions [C. J. Benmore and Y. L. Loh, J. Chem. Phys. 112, 5877 (2000)]. The simultaneous variation of the atomic charges and the intramolecular potentials requires modified equations of motion and a multiple time step algorithm has been implemented to solve these equations. The article concludes with a discussion of the bulk structure and properties with an emphasis on the hydrogen bonding network.     

The flexible, polarizable, thole-type interaction potential for water (TTM2-F) revisited

We present a revision of the flexible, polarizable, Thole-type interaction potential for water [J. Chem. Phys.2002, 116, 5115], which allows for condensed-phase simulations. The revised version (TTM2.1-F) of the potential correctly describes the individual water molecular dipole moment and alleviates problems arising at short intermolecular separations that can be sampled in the course of molecular dynamics and Monte Carlo simulations of condensed environments. Furthermore, its parallel implementation under periodic boundary conditions enables the efficient calculation of the macroscopic structural and thermodynamic properties of liquid water, as its performance scales superlinearly with up to a number of 64 processors for a simulation box of 512 molecules. We report the radial distribution functions, average energy, internal geometry, and dipole moment in the liquid as well as the density, dielectric constant, and self-diffusion coefficient at T = 300 K from (NVT) and (NPT) classical molecular dynamics simulations by using the revised version of the potential.     

Analytical angular solutions for the atom–diatom interaction potential in a basis set of products of two spherical harmonics: two approaches

We present general analytical expressions for the matrix elements of the atom–diatom interaction potential, expanded in terms of Legendre polynomials, in a basis set of products of two spherical harmonics, especially significant to the recently developed adiabatic variational theory for cold molecular collision experiments [J. Chem. Phys. 143, 074114 (2015); J. Phys. Chem. A 121, 2194 (2017)]. We used two approaches in our studies. The first involves the evaluation of the integral containing trigonometric functions with arbitrary powers. The second approach is based on the theorem of addition of spherical harmonics.

     

Recovering the Crooks equation for dynamical systems in the isothermal-isobaric ensemble: a strategy based on the equations of motion

The Crooks equation [Eq. (10) in J. Stat. Phys. 90, 1481 (1998)] relates the work done on a system during a nonequilibrium transformation to the free energy difference between the final and the initial state of the transformation. Recently, the authors have derived the Crooks equation for systems in the canonical ensemble thermostatted by the Nose-Hoover or Nose-Hoover chain method [P. Procacci et al., J. Chem. Phys. 125, 164101 (2006)]. That proof is essentially based on the fluctuation theorem by Evans and Searles [Adv. Phys. 51, 1529 (2002)] and on the equations of motion. Following an analogous approach, the authors derive here the Crooks equation in the context of molecular dynamics simulations of systems in the isothermal-isobaric (NPT) ensemble, whose dynamics is regulated by the Martyna-Tobias-Klein algorithm [J. Chem. Phys. 101, 4177 (1994)]. Their present derivation of the Crooks equation correlates to the demonstration of the Jarzynski identity for NPT systems recently proposed by Cuendet [J. Chem. Phys. 125, 144109 (2006)].     

Vibronic coupling in asymmetric bichromophores: Theory and application to diphenylmethane

The theory for modeling vibronic interactions in bichromophores was introduced in sixties by Witkowski and Moffitt [J. Chem. Phys. 33, 872 (1960)] and extended by Fulton and Gouterman [J. Chem. Phys. 35, 1059 (1961)]. The present work describes extension of this vibronic model to describe bichromophores with broken vibrational symmetry such as partly deuterated molecules. Additionally, the model is extended to include inter-chromophore vibrational modes. The model can treat multiple vibrational modes by employing Lanczos diagonalization procedure of sparse matrices. The developed vibronic model is applied to simulation of vibronic spectra of flexible bichromophore diphenylmethane and compared to high-resolution experimental spectra [J. A. Stearns, N. R. Pillsbury, K. O. Douglass, C. W. Mu?ller, T. S. Zwier, and D. F. Plusquellic, J. Chem. Phys. 129, 224305 (2008)].     

Diffusion in linear porous media with periodic entropy barriers: A tube formed by contacting spheres

The problem of transport in quasi-one-dimensional periodic structures has been studied recently by several groups [D. Reguera et al., Phys. Rev. Lett.96, 130603 (2006); P. S. Burada et al., Phys. Rev. E75, 051111 (2007); B. Q. Ai and L. G. Liu, ibid.74, 051114 (2006); B. Q. Ai et al., ibid.75, 061126 (2007); B. Q. Ai and L. G. Liu, J. Chem. Phys.126, 204706 (2007); 128, 024706 (2008); E. Yariv and K. D. Dorfman, Phys. Fluids19, 037101 (2007); N. Laachi et al., Europhys. Lett.80, 50009 (2007); A. M. Berezhkovskii et al., J. Chem. Phys.118, 7146 (2003); 119, 6991 (2003)]. Using the concept of "entropy barrier" [R. Zwanzig, J. Phys. Chem.96, 3926 (1992)] one can classify such structures based on the height of the entropy barrier. Structures with high barriers are formed by chambers, which are weakly connected with each other because they are connected by small apertures. To escape from such a chamber a diffusing particle has to climb a high entropy barrier to find an exit that takes a lot of time [I. V. Grigoriev et al., J. Chem. Phys.116, 9574 (2002)]. As a consequence, the particle intrachamber lifetime tau(esc) is much larger than its intrachamber equilibration time, tau(rel), tau(esc)>tau(rel). When the aperture is not small enough, the intrachamber escape and relaxation times are of the same order and the hierarchy fails. This is the case of low entropy barriers. Transport in this case is analyzed in the works of Schmid and co-workers, Liu and co-workers, and Dorfman and co-workers, while the work of Berezhkovskii et al. is devoted to diffusion in the case of high entropy barriers.     

Insights in quantum dynamical effects in the infrared spectroscopy of liquid water from a semiclassical study with an ab initio-based flexible and polarizable force field

The dynamical properties of liquid water play an important role in many processes in nature. In this paper, we focus on the infrared (IR) absorption spectrum of liquid water based on the linearized semiclassical initial value representation (LSC-IVR) with the local Gaussian approximation (LGA) [J. Liu and W. H. Miller, J. Chem. Phys. 131, 074113 (2009)] and an ab initio based, flexible, polarizable Thole-type model (TTM3-F) [G. S. Fanourgakis and S. S. Xantheas, J. Chem. Phys. 128, 074506 (2008)]. Although the LSC-IVR (LGA) gives the exact result for the isolated three-dimensional shifted harmonic stretching model, it yields a blueshifted peak position for the more realistic anharmonic stretching potential. By using the short-time information of the LSC-IVR correlation function; however, it is shown how one can obtain more accurate results for the position of the stretching peak. Due to the physical decay in the condensed phase system, the LSC-IVR (LGA) is a good and practical approximate quantum approach for the IR spectrum of liquid water. The present results offer valuable insight into future attempts to improve the accuracy of the TTM3-F potential or other ab initio-based models in reproducing the IR spectrum of liquid water.     

Exchange energy gradients with respect to atomic positions and cell parameters within the Hartree-Fock Gamma-point approximation

Recently, linear scaling construction of the periodic exact Hartree-Fock exchange matrix within the Gamma-point approximation has been introduced [J. Chem. Phys. 122, 124105 (2005)]. In this article, a formalism for evaluation of analytical Hartree-Fock exchange energy gradients with respect to atomic positions and cell parameters at the Gamma-point approximation is presented. While the evaluation of exchange gradients with respect to atomic positions is similar to those in the gas phase limit, the gradients with respect to cell parameters involve the accumulation of atomic gradients multiplied by appropriate factors and a modified electron repulsion integral (ERI). This latter integral arises from use of the minimum image convention in the definition of the Gamma-point Hartree-Fock approximation. We demonstrate how this new ERI can be computed with the help of a modified vertical recurrence relation in the frame of the Obara-Saika and Head-Gordon-Pople algorithm. As an illustration, the analytical gradients have been used in conjunction with the QUICCA algorithm [K. Nemeth and M. Challacombe, J. Chem. Phys. 121, 2877 (2004)] to optimize periodic systems at the Hartree-Fock level of theory.     

Properties of the H-bond network for two-dimensional lattice water model

A microscopic Hamiltonian of the hydrogen-bond network in two-dimensional lattice water is proposed, which describes the formation and disruption of the H bonds, their bending, and which satisfies the Bernal-Fowler rules [J. D. Bernal and R. H. Fowler, J. Chem. Phys. 1, 515 (1933)]. The thermodynamic properties of the H-bond network are studied using the method of many-particle irreducible distribution functions, which is a generalization of the Kikuchi cluster approach [R. Kikuchi, Phys. Rev. 81, 988 (1951)] and the Bethe-Peierls quasiactivities method [H. A. Bethe, Prog. R. Soc. A 150, 552 (1935)]. The temperature dependencies of the average number of H bonds per molecules, the contribution of the H bonds into the heat capacity of the system, and the parameters describing the correlations between the states of molecules on the neighboring sites are investigated. It is shown that depending on the magnitude of the interaction between the H bonds in the H-bond subsystem either smooth or sharp first-order phase transition can occur. The role of different factors in the formation of the properties of the H-bond network is discussed.     

State-to-state reactive differential cross sections for the H+H2-->H2+H reaction on five different potential energy surfaces employing a new quantum wavepacket computer code: DIFFREALWAVE

State-to-state differential cross sections have been calculated for the hydrogen exchange reaction, H+H2-->H2+H, using five different high quality potential energy surfaces with the objective of examining the sensitivity of these detailed cross sections to the underlying potential energy surfaces. The calculations were performed using a new parallel computer code, DIFFREALWAVE. The code is based on the real wavepacket approach of Gray and Balint-Kurti [J. Chem. Phys. 108, 950 (1998)]. The calculations are parallelized over the helicity quantum number Omega' (i.e., the quantum number for the body-fixed z component of the total angular momentum) and wavepackets for each J,Omega' set are assigned to different processors, similar in spirit to the Coriolis-coupled processors approach of Goldfield and Gray [Comput. Phys. Commun. 84, 1 (1996)]. Calculations for J=0-24 have been performed to obtain converged state-to-state differential cross sections in the energy range from 0.4 to 1.2 eV. The calculations employ five different potential energy surfaces, the BKMP2 surface and a hierarchical family of four new ab initio surfaces [S. L. Mielke, et al., J. Chem. Phys. 116, 4142 (2002)]. This family of four surfaces has been calculated using three different hierarchical sets of basis functions and also an extrapolation to the complete basis set limit, the so called CCI surface. The CCI surface is the most accurate surface for the H3 system reported to date. Our calculations of differential cross sections are the first to be reported for the A2, A3, A4, and CCI surfaces. They show that there are some small differences in the cross sections obtained from the five different surfaces, particularly at higher energies. The calculations also show that the BKMP2 performs well and gives cross sections in very good agreement with the results from the CCI surface, displaying only small divergences at higher energies.     

Computer simulation study of metastable ice VII and amorphous phases obtained by its melting

Molecular dynamics simulations of metastable ice VII and cubic ice Ic are carried out in order to examine (1) the ability of commonly used water interaction potentials to reproduce the properties of ices, and (2) the possibility of generating low-density amorphous (LDA) structures by heating ice VII, which is known to transform to LDA at approximately 135 K at normal pressure [S. Klotz, J. M. Besson, G. Hamel, R. J. Nelmes, J. S. Loveday, and W. G. Marshall, Nature (London) 398, 681 (1999)]. We test four simple empirical interaction potentials of water: TIP4P [W. L. Jorgensen, J. Chandrasekhar, J. D. Madura, R. W. Impey, and M. L. Klein, J. Chem. Phys. 79, 926 (1983)], SPC/E [H. J. C. Berendsen, J. R. Grigera, and T. P. Straatsma, J. Phys. Chem. B 91, 6269 (1987)], TIP5P [M. W. Mahoney and W. L. Jorgensen, J. Chem. Phys. 112, 8910 (2000)], and ST2 [F. H. Stillinger and A. Rahman, J. Chem. Phys. 60, 1545 (1974)]. We have found that TIP5P ice VII melts at 210 K, TIP4P at 90 K, and SPC/E at 70 K. Only TIP5P water after transition has a structure similar to that of LDA. TIP4P and SPC/E have almost identical structures, dissimilar to any known water or amorphous phases, but upon heating both slowly evolve towards LDA-like structure. ST2 ice VII is remarkably stable up to 430 K. TIP4P and SPC/E predict correctly the cubic ice collapse into a high-density amorphous ice (HDA) at approximately 1 GPa whereas TIP5P remains stable up to approximately 5 GPa. The densities of the simulated ice phases differ significantly, depending on the potential used, and are generally higher than experimental values. The importance of proper treatment of long-range electrostatic interactions is also discussed.     

Local response dispersion method. II. Generalized multicenter interactions

Recently introduced local response dispersion method [T. Sato and H. Nakai, J. Chem. Phys. 131, 224104 (2009)], which is a first-principles alternative to empirical dispersion corrections in density functional theory, is implemented with generalized multicenter interactions involving both atomic and atomic pair polarizabilities. The generalization improves the asymptote of intermolecular interactions, reducing the mean absolute percentage error from about 30% to 6% in the molecular C(6) coefficients of more than 1000 dimers, compared to experimental values. The method is also applied to calculations of potential energy curves of molecules in the S22 database [P. Jure?ka et al., Phys. Chem. Chem. Phys. 8, 1985 (2006)]. The calculated potential energy curves are in a good agreement with reliable benchmarks recently published by Molnar et al. [J. Chem. Phys. 131, 065102 (2009)]. These improvements are achieved at the price of increasing complexity in the implementation, but without losing the computational efficiency of the previous two-center (atom-atom) formulation. A set of different truncations of two-center and three- or four-center interactions is shown to be optimal in the cost-performance balance.     

Isotope effects in liquid water by infrared spectroscopy. II. Factor analysis of the temperature effect on H2O and D2O

Some 500 infrared (IR) spectra of light and heavy waters were obtained between 29 and 93 degrees C in order to identify the species present in liquid water. Factor analysis of these gives two species for each type of water with their IR spectra and abundance curves. Using an orthogonalization procedure, we obtained the temperature factor limits of -22 and +118 degrees C (+/-5 degrees C) that we coined cold and hot factors, respectively. Within experimental error, these limits are the same for light and heavy waters. The spectra of the orthogonalized factors presented show a decrease of the OH (OD) stretch band integrated intensities of almost 36% from the cold to the hot factors. No "free" OH (OD) group is present or formed in the temperature ramp. This indicates that all water molecules in the bulk are made of an oxygen atom surrounded with four hydrogen atoms, two covalently bonded, and two hydrogen bonded. This is consistent with the previous study of mixtures of H(2)O and D(2)O [part I: J.-J. Max and C. Chapados, J. Chem. Phys. 116, 4626 (2002)]. To maintain the ordinary liquid within the limits of 0 and 100 degrees C at atmospheric pressure, a fraction of the cold and hot factors are necessary. With the spectra of the cold and hot factors and the abundance curves, one can generate the spectrum at any temperature between -22 and +118 degrees C of light and heavy liquid waters.     

Rigid-body dynamics in the isothermal-isobaric ensemble: a test on the accuracy and computational efficiency

We have developed a time-reversible rigid-body (rRB) molecular dynamics algorithm in the isothermal-isobaric (NPT) ensemble. The algorithm is an extension of rigid-body dynamics [Matubayasi and Nakahara, J Chem Phys 1999, 110, 3291] to the NPT ensemble on the basis of non-Hamiltonian statistical mechanics [Martyna, G. J. et al., J Chem Phys 1994, 101, 4177]. A series of MD simulations of water as well as fully hydrated lipid bilayer systems have been undertaken to investigate the accuracy and efficiency of the algorithm. The rRB algorithm was shown to be superior to the state-of-the-art constraint-dynamics algorithm SHAKE/RATTLE/ROLL, with respect to computational efficiency. However, it was revealed that both algorithms produced accurate trajectories of molecules in the NPT as well as NVT ensembles, as long as a reasonably short time step was used. A couple of multiple time-step (MTS) integration schemes were also examined. The advantage of the rRB algorithm for computational efficiency increased when the MD simulation was carried out using MTS on parallel processing computer systems; total computer time for MTS-MD of a lipid bilayer using 64 processors was reduced by about 40% using rRB instead of SHAKE/RATTLE/ROLL.     

Density functional theory of adsorption of mixtures of charged chain particles and spherical counterions

We propose a microscopic density functional theory to describe nonuniform ionic fluids composed of chain molecules with charged "heads" and spherical counterions. The chain molecules are modeled as freely jointed chains of hard spheres, the counterions are oppositely charged spheres of the same diameter as all segments of chain molecules. The theory is based on the approach of Yu and Wu [J. Chem. Phys. 117, 2368 (2002)] of adsorption of chain molecules and on theory of adsorption of electrolytes [O. Pizio, A. Patrykiejew, and S. Sokolowski, J. Chem. Phys. 121, 11957 (2004)]. As an application of the proposed formalism we investigate the structure and adsorption of fluids containing segments of different length in a slitlike pore.     

Contribution of water dimer absorption to the millimeter and far infrared atmospheric water continuum

We present a rigorous calculation of the contribution of water dimers to the absorption coefficient alpha(nu,T) in the millimeter and far infrared domains, over a wide range (276-310 K) of temperatures. This calculation relies on the explicit consideration of all possible transitions within the entire rovibrational bound state manifold of the dimer. The water dimer is described by the flexible 12-dimensional potential energy surface previously fitted to far IR transitions [C. Leforestier et al., J. Chem. Phys. 117, 8710 (2002)], and which was recently further validated by the good agreement obtained for the calculated equilibrium constant Kp(T) with experimental data [Y. Scribano et al., J. Phys. Chem. A. 110, 5411 (2006)]. Transition dipole matrix elements were computed between all rovibrational states up to an excitation energy of 750 cm(-1), and J=K=5 rotational quantum numbers. It was shown by explicit calculations that these matrix elements could be extrapolated to much higher J values (J=30). Transitions to vibrational states located higher in energy were obtained from interpolation of computed matrix elements between a set of initial states spanning the 0-750 cm(-1) range and all vibrational states up to the dissociation limit (approximately 1200 cm(-1)). We compare our calculations with available experimental measurements of the water continuum absorption in the considered range. It appears that water dimers account for an important fraction of the observed continuum absorption in the millimeter region (0-10 cm(-1)). As frequency increases, their relative contribution decreases, becoming small (approximately 3%) at the highest frequency considered nu=944 cm(-1).