Monte Carlo computations have been used to study the behavior of the TIPS2 model of water next to various surfaces. This paper reports on the model next to inert plane surfaces at 399 K. This study serves as a preliminary to studies of the model next to other surfaces especially interesting for electrochemistry; these are reported in an accompanying paper. Meanwhile the surface structure and polarization are examined and compared with those previously reported for other water models. There is a mild tendency to form an oriented ice structure near the surface, as had been found [C. Y. Lee, J. A. McCammon, and P. J. Rossky, J. Chem. Phys. 80, 4448 (1984)] for the ST2 model. Extremely vexing convergence problems were encountered, and as well an unexpected anisotropy of the angular distributions. The latter problem turns out to involve questions of the sample shape and boundary conditions and seems to have worrying implications even for simulations of homogeneous fluids of nonspherical particles.

1.
For an excellent review, see
S. L.
Carnie
and
G. M.
Torrie
,
Adv. Chem. Phys.
56
,
141
(
1984
).
2.
See, for example, work cited in P. Delahay, Double Layer and Electrode Kinetics (Wiley, New York, 1965), and in M. J. Sparnaay, The Electrical Double Layer (Pergamon, Oxford, 1972);
also for example
R.
Parsons
,
Electroanal. Chem. Interfac. Electrochem.
59
,
229
(
1975
), and
W. R.
Fawcett
,
J. Chem. Phys.
82
,
1385
(
1978
).
3.
D. Y. C.
Chan
,
D. J.
Mitchell
, and
B. W.
Ninham
,
J. Chem. Phys.
72
,
5159
(
1980
);
S. L.
Carnie
and
D. Y. C.
Chan
,
J. Chem. Phys.
75
,
3485
(
1981
), and other references to the work of Chan, Carnie et al. cited there; ,
J. Chem. Phys.
L.
Blum
and
D.
Henderson
,
J. Chem. Phys.
74
,
1902
(
1981
); ,
J. Chem. Phys.
D.
Henderson
,
L.
Blum
, and
M.
Lozada‐Cassou
,
J. Electroanal. Chem.
150
,
291
(
1983
);
F.
Vericat
,
L.
Blum
, and
D.
Henderson
,
J. Chem. Phys.
77
,
5808
(
1982
);
G. N. Patey and G. M. Torrie (unpublished).
4.
W. L.
Jorgensen
,
J. Chem. Phys.
77
,
4156
(
1982
).
5.
W. L.
Jorgensen
,
J.
Chandrasekhar
,
J. D.
Madura
,
R.
Impey
, and
M. L.
Klein
,
J. Chem. Phys.
79
,
926
(
1983
).
6.
B.
Jo/nsson
,
Chem. Phys. Lett.
82
,
520
(
1981
).
7.
N. I.
Christou
,
J. S.
Whitehouse
,
D.
Nicholson
, and
N. G.
Parsonage
,
Chem. Soc. Faraday Symp.
16
,
139
(
1981
).
8.
M.
Marchesi
,
Chem. Phys. Lett.
97
,
224
(
1983
).
9.
N.
Anastasiou
,
D.
Fincham
, and
K.
Singer
,
Chem. Soc. Faraday Trans. 2
79
,
1639
(
1983
).
10.
R.
Sonnenschein
and
K.
Heinziger
,
Chem. Phys. Lett.
102
,
550
(
1983
).
11.
G.
Barabino
,
C.
Gavotti
, and
M.
Marchesi
,
Chem. Phys. Lett.
104
,
478
(
1984
).
12.
C. Y.
Lee
,
J. A.
McCammon
, and
P. J.
Rossky
,
J. Chem. Phys.
80
,
4448
(
1984
).
13.
N.
Metropolis
,
A. W.
Metropolis
,
M. N.
Rosenbluth
,
A. H.
Teller
, and
E.
Teller
,
J. Chem. Phys.
21
,
1087
(
1953
).
14.
J. P. Valleau and S. G. Whittington, Statistical Mechanics, Part A: Equilibrium Techniques, edited by B. J. Berne, Vol. 5 of Modern Theoretical Chemistry (Plenum, New York, 1977), p. 137.
15.
G.
Marsaglia
,
Ann. Math. Stat.
43
,
645
(
1972
).
16.
G. M.
Torrie
and
J. P.
Valleau
,
J. Chem. Phys.
73
,
5807
(
1980
);
J. P. Valleau and G. M. Torrie, in The Problems of Long‐Range Forces in Computer Simulation of Condensed Media, NRCC Proceedings No. 9 (Lawrence Berkeley Laboratory, California, 1980), p. 65.
17.
The energy of interaction Uds of the charges qi on a TIPS2 particle, with a set of n such sheets carrying dipole densities ρ1(l = 1,2,…n), is given explicitly by
with
where (xi, yi, zil) are the coordinates of the charge ρi with respect to the center of the square hole in the l th sheet. The sum over i is over the three charges on the particle: M and the two H’s.
18.
Reference 14 and
D. N.
Card
and
J. P.
Valleau
,
J. Chem. Phys.
52
,
6232
(
1970
).
This content is only available via PDF.
You do not currently have access to this content.