Molecular Dynamics Study of Water Clathrates (P. L. M. Plummer

Whalley: Halos are usually caused by ice crystals in cirrus clouds, and presumably the octahedral crystals of cubic ice that appear to cause the 28' h...
12 downloads 0 Views 334KB Size
The Journal of Physical Chemlstry, Vol. 87, No. 27, 1983 4335

(“quasi-liquidlayer”). When the relative amount of the fluctuating mobile phase with respect to the bulk ice is higher, a cooperative fluctuation from the surface goes more easily into the interior of the bulk ice. The fluctuation on the surface must increase in this fashion. B. R. Stauffer: At what partial COz pressure are the monomolecular layer and the formation of clathrates observed? Ocampo: It depends on the temperature. In general, it appears for reduced pressures x = P/Psgreater than 0.2. In some cases a coverage of one monolayer is attained before the clathrate begins to nucleate. All the pressures in adsorption measurements were up to 1bar. We cannot use our results to compare with the COz in polar ice cores where the partial pressure is very low.

D. E. Barnaal: How did you establish that a C 0 2 clathrate appeared on the ice surface? Ocampo: In powder samples the pressure falls to the decomposition value.

Optimal Conditions of the Homogeneous Nucleation of Cubic Ice (IC)from Concentrated Solutions of LiCl 0 D20 (J. Dupuy) J. Perez: Do your results mean that the nucleation process implies the formation of only cubic ice? That seems evident with your results obtained at low temperature; for higher temperatures, growth might be too rapid, preventing the observation of cubic ice as, above a certain size, there is a transformation IC Ih. Dupuy: If the concentration “c” of LicP in D20 is 9.75 < c < 12 mol %, the nucleation process implies the formation of only cubic ice a t low temperatures near Tr If the size of cubic ice crystallites exceeds 180 A, there is a contamination of cubic ice by hexagonal ice, and the transformation IC Ih occurs.

-

-

Cubic Ice in Nature (E. Whalley) F. Stillinger: I suppose you have assumed that cubic ice in nature has the same refractive index as that of the hexagonal form. Since we cannot be certain to what extent that natural ice IC is proton disordered (or partially ordered), is it not possible that the refractive index is correspondingly uncertain? Whalley: Ice ICas made in the laboratory is certainly proton disordered, as is shown by its infrared’s2 and bans spectra and by its dielectric proper tie^.^ There is no reason to expect natural ice ICto differ from artificial ice in ita degree of order. The density of ice Ih and IC are identical within experimental uncertainty. The refractive index of ice Ih is almost isotropic, and so the refractive index of ice IC can be estimated to well within the accuracy needed for predicting halos. In any case, if ice ICwere partly ordered its symmetry would probably be degraded from face-centered cubic and it would not crystallize in octahedra. Bertie, J. E.; Whalley, E. J. Chem. Phys. 1964, 40, 1657. Bertie, J. E.; Whalley, E. Ibid. 1967, 46, 1271. Taylor, M. J.; Whalley, E. Ibid. 1964, 40, 1660. Gough, S. R.; Davidson, D. W. Ibid. 1970,52, 5442. J. W. Glen: Relative to Stillinger’s question, is ice IC were proton ordered its relative stability with ice Ih would vary strongly with temperature which does not seem to happen. As regards the tendency to observe it in Doriis, could this be due to the presence there of some mineral aerosol, possibly from industrial pollution, which is favorable for epitaxial nucleation of ice IC? Whalley: The industrial pollution would have to be in Paris from at least 1629 to 1920, with possibly a gap of 160 years. In any case, the cirrus clouds that usually form the halos may be too high to be much affected by emissions immediately below them. G . P. Johari: Do you think that the possibility of Scheiner’s halos as arising from octahedral clathrate hydrate crystals in the atmosphere could be completely eliminated? Whalley: The equilibrium vapor pressure of the guest molecule over any clathrate that might form from the atmosphere is much higher than its partial pressure in the atmosphere. A clathrate could, therefore, be formed only under highly metastable conditions. I do not know of any evidence that this might occur.

E. Hindman: What can you say about the atmosphere conditions which give rise to cubic ice, temperature, pressures, etc. Whalley: Halos are usually caused by ice crystals in cirrus clouds, and presumably the octahedral crystals of cubic ice that appear to cause the 28’ halo are formed there. The pressure and temperature appear to be roughly 0.2 bar and 230 K.

Molecular Dynamics Study of Water Clathrates (P. L. M. Plummer) R. Smoluchowski: (1) How do the hydrogens change their positions as function of time and temperature? (2) Do you think that a structure made of several cages would show different stability? Plummer: (1)At the lower temperatures (below 200 K), the hydrogens remain localized between the same two oxygens throughout the simulation although substantial deviation in 0-H- - -0angle occurs. As the temperahre increases, more proton “switching” reactions and double hydrogen bonds between a pair of oxygens occur. The lifetime of a specific hydrogen bond decreases. (2) I believe the basic trends would be the same. However, I would expect the temperature of the transition to be somewhat higher. With multiple cages, fewer of the molecules would be in the surface and the internal structure would persist longer.

Molecular Dynamics Studies of Ice ICand the Structure I Clathrate Hydrates of Methane (J. S. Tse) F. Stillinger: How did you decide which way to orient water molecules in the two networks, and was there any attempt to average over such canonical orientations? Tse: Since the water molecules in both the clathrate hydrate and ice ICare rotationally disordered, initial proton configurations were selected in a random fashion to conform to the Bernal-Fowler rules. Only one configuration was considered in each simulation and we make no attempt to average over different canonical orientations. P. Halleck: When you say that the empty hydrate is stable at high pressure, do you mean thermodynamically stable relative to high pressure ice (say ice VII)? Tse: No. I make no attempt to compare the stability of the hydrates with other forms of ice.

R. Smoluchowski: To my knowledge the empty clathrate has never been observed. How can you claim that it is stable below 230 K? With respect to what? Tse: Within the time scale of our molecular dynamic simulation, the fluctuations in the total energy and the pressure of the system from their mean values are very small. We observe no indication that the empty hydrate structure is transforming into other phases and it appears to be mechanically stable. G. P. Johari: I think you can determine the heat capacity of an empty clathrate lattice by extrapolating to zero occupancy the data obtained for partially occupied cages. Tse: Yes, in principle you can do so. Unfortunately, it is very difficult to prepare hydrate samples within a wide range of guest occupancies.

Ice-Rich Moons and the Physical Properties of Ice (C. J. Consolmagno) J. Klinger: Why did you exclude C02 when you enumerated primary molecules that formed in the primitive nebula? Consolmagno: In a gas of solar composition, rich in hydrogen and relatively poor in oxygen and carbon, it can be shown that CO is thermodynamically stable relative to COz. The relevant reaction is CO + H2 F? C02 + H2;so long as the [HI/ [H20]ratio is fixed, by cosmic abundances, to be much greater than unity, then [CO]K,![CO,] must also be much greater than unity. For the P-T conditions in the solar nebula, K , is always close enough to unity that [CO] will be greater than [CO,]. S. Warren: You showed pictures of bright white moons and suggested that the high reflectivity was an indication of ice. But ice is actually highly transparent at visible wavelengths and

4336

The Journal of Physical Chemistty, Vol. 87, No. 21, 1983

therefore very dark. It only becomes highly reflective when it is broken up into small separated particles. So the high reflectivity is actually evidence for snow. It may be possible to use the spectral reflectivity to estimate the sizes of these small particles, as well as their composition (Wiscombe; Warren J. Atmos. Sci. 1980,37, 2712). Consolmagno: Polarization studies can also reveal the small particle size. Since these are airless bodies, subjected to constant micrometeorite bombardment, we expect the surfaces to be broken up into small, snowlike fragments. These bodies also have a very strong radar reflectivity, which is still not well understood but which may be evidence for a complex surface structure full of voids and fractures.

T. Takahashi: May I ask why the ice coverage on planets is scattered in appearance? Consolmagno: We must remember that these are bodies made up primarily of ice; your question could be turned around and asked, why is the rocky coverage of the ice scattered? The dark material may be micrometeorites which have contaminated the surface, or may be a thin layer of rocky material from the primordial ice-rock mixture that made up these bodies which was left behind after some process (sublimation or ion sputtering) has stripped away the ice. We see craters with icy floors that appear to have “punctured” this layer of rocky material. Icy regions on Ganymede appear to be younger, resurfaced regions or associated with extensional “cracks”. Some transport of icy material to polar regions does occur, but since there is no atmosphere and the temperatures are quite low, it is not a very efficient process.

Extraterrestrial Ice. A Review (J. Klinger) R. Smoluchowski: There may be plenty of water on Mars. High temperatures occur on Mars every few 1OOOOO years. Enceladus surface temperature is approximately 75 K. Thus there is a need for considerable heat to cause its crystallization near 150 K. The heat evolved below 150 K is very small. Klinger: In fact we need plenty of heat to heat up Enceladus to the transition temperature, but I think that the model of tidal heating as it was proposed by Yoder (ref 31) will work much easier when phase transitions take place. Further it has been shown by Ghormley (ref 33) that amorphous ice will release some heat even before the phase transition. What may help too is resurfacing. J . C. Comiso: What additional ice observational data would you require to gain a better understanding of the origin of the solar system? Klinger: I think that the identification of different ice phases (for example by infrared reflection spectra) may give us some ideas the thermal history of the material. A . H. Delsemme: You mentioned that the OH band was observed in comets by radio astronomers. This is true. However, it was discovered in its 3050-A ultraviolet band in the spectrum of comet Cummingham (1941) by Swings and co-workers, and repeatedly observed in many comets in the ultraviolet. Klinger: OK!

Ice in Comets (A. H. Delsemme) J . Klinger: May the fact that Kohoutek was essentially outgassing COz be related to the finding that after perihelion the production of OH, in a wide range of heliocentric distances, showed no variation? Delsemme: Yes indeed; on its leg before perihelion, the surface temperature of the nucleus was controlled by C02 sublimation; apparently, COz became depleted near perihelion, but took a long time to restore a new steady state involving only water ice; this explains why the sublimation of water ice remained for a while a t a constant rate, the rise in temperature due to C02depleting being compensated by the larger and larger heliocentric distance of the comet.

Energetic Charged Particle Erosion of Ices in the Solar System (R. E. Johnson) J. M . Warman: You seemed to be attaching some significance to the fact that several molecules were being detached from the

surface per incident 1.5-MeV He+ ion. However, the total number of molecules undergoing excitation and dissociation will in fact be in the region of lo4 to lo6 per ion. Detachment therefore corresponds to a very small fraction of excitation events. A likely process of molecular detachment would seem to be molecular dissociation of surface molecules followed by escape of high velocity (“hot”)radical fragments. Would you be able tp detect such radical products and if so what is their ratio to detached water molecules? Johnson: It is true that an enormous number of bonds are broken as the electronic excited states relax. In fact decomposition apparently occurs fairly efficiently. The g value is of the order of 0.6 for a-particles. We see only a fraction of the O2 and Dz, those produced and ejected in the surface region. When I discussed the modification of the surface layer of a comet irradiated in the Oort cloud I was referring to the total effect you mentioned. I used the word large, I meant the number ejected from the surface was large compared to sputtering produced by direct collisions. Radical fragments are probably lost, but our estimates indicate they are not the dominant ejected species. We can detect these radicals by using water having isotopes of oxygen. Certainly if the events at the surface were produced by individual excitation events radicals may be dominant products of dissociative electronic relaxation as,for example, in desorption. However, the nonlinear dependences of yield on (dE/dn) suggests that the relaxation events might be producing “hot” overlapping regions which produces the observed molecular ejection.

R. Smoluchowski: You indicated that Saturnian rings are eroded photoelectrically rather than by particulate radiation. How did you estimate the electronic and proton flux in proximity of the rings? Is this based on Morfill, Grein et al. data? Johnson: The fluxes were obtained by modeling the diffusive loss in the vicinity of the main ring system. It is based on a recent paper by Cheng, Lanzerotti, and Pirronello. J. M. Greenberg: What is the depth of penetration of a 1.5-MeV He+ ion? If it is deep enough you are here dealing with bulk (well inside surface) phenomena differentiate. Johnson: A 1.5-MeV He+ ion will penetrate to micron depths. However, we believe the energy deposited only within a few 100 A at most would affect the single particle surface erosion. If the erosion continues or if the surface is subsequently heated the effects a t greater depth will eventually contribute.

R. Wolff: The energy of ejected H 2 0 from your experiment is generally