PATH-INTEGRAL SIMULATIONS OF MIXED PARA-D-2 AND ORTHO-D-2 CLUSTERS - THE ORIENTATIONAL EFFECTS

A recently suggested computational scheme [Buch and Devlin, J. Chem. P hys. 98, 4195 (1993)] is extended to describe orientational phenomena in systems containing assemblies of odd-J and even-J hydrogen and deut erium molecules. Nuclear symmetry and rotational effects are incorpora ted in the path integral Monte Carlo algorithm. The scheme is employed in a study of the (para-D-2)(3)(ortho-D-2)(10) and (para-D-2)(13) clu sters in the 1-3 K temperature range. In the (para-D-2)(3)(ortho-D-2)( 10) cluster at 1 K, the three J=1 molecules form a ''cluster within cl uster,'' held. together by the anisotropic quadrupole-quadrupole inter actions. At 3 K the three molecule J=1 cluster is dissolved within the (D-2)(13) cluster. The orientational structure of the (para-D-2)(13) cluster at 1 K is composed of two distinct distorted-T pair configurat ions of neighboring J=1 molecules. Upon heating to 3 K the orientation al structure changes, and a clear separation between the two pair geom etries is no longer obtained. Moreover at 3 K thermal distortion of th e cluster structure from icosahedral relaxes the extent of orientation al frustration of the system, and thus contributes to the lowering of the rotational energy.