TY - JOUR

T1 - (HCl)2 and (HF)2 in small helium clusters

T2 - Quantum solvation of hydrogen-bonded dimers

AU - Jiang, Hao

AU - Sarsa, Antonio

AU - Murdachaew, Garold

AU - Szalewicz, Krzysztof

AU - Bačić, Zlatko

N1 - Funding Information:
Two of the authors (Z.B. and H.J.) have been supported in part by the National Science Foundation Grant No. CHE-0315508. Another two of the authors (K.S. and G.M.) acknowledge support from the NSF through Grant No. CHE-0239611. One author (A.S.) acknowledges support from the Spanish Dirección de Investigación Cientifica y Technica, Grant No. ESPB2002–00200, and partial support from the Junta de Andalucia.

PY - 2005

Y1 - 2005

N2 - We present a rigorous theoretical study of the solvation of (HCl)2 and (HF)2 by small (He4) n clusters, with n=1-14 and 30. Pairwise-additive potential-energy surfaces of Hen (HX)2 (X=Cl and F) clusters are constructed from highly accurate four-dimensional (rigid monomer) HX-HX and two-dimensional (rigid monomer) He-HX potentials and a one-dimensional He-He potential. The minimum-energy geometries of these clusters, for n=1-6 in the case of (HCl)2 and n=1-5 for (HF)2, correspond to the He atoms in a ring perpendicular to and bisecting the HX-HX axis. The quantum-mechanical ground-state energies and vibrationally averaged structures of Hen (HCl)2 (n=1-14 and 30) and Hen (HF)2 (n=1-10) clusters are calculated exactly using the diffusion Monte Carlo (DMC) method. In addition, the interchange-tunneling splittings of Hen (HCl)2 clusters with n=1-14 are determined using the fixed-node DMC approach, which was employed by us previously to calculate the tunneling splittings for Hen (HF)2 clusters, n=1-10 [A. Sarsa, Phys. Rev. Lett. 88, 123401 (2002)]. The vibrationally averaged structures of Hen (HX)2 clusters with n=1-6 for (HCl)2 and n=1-5 for (HF)2 have the helium density localized in an effectively one-dimensional ring, or doughnut, perpendicular to and at the midpoint of the HX-HX axis. The rigidity of the solvent ring varies with n and reaches its maximum for the cluster size at which the ring is filled, n=6 and n=5 for (HCl)2 and (HF)2, respectively. Once the equatorial ring is full, the helium density spreads along the HX-HX axis, eventually solvating the entire HX dimer. The interchange-tunneling splitting of Hen (HCl)2 clusters hardly varies at all over the cluster size range considered, n=1-14, and is virtually identical to that of the free HCl dimer. This absence of the solvent effect is in sharp contrast with our earlier results for Hen (HF)2 clusters, which show a ∼30% reduction of the tunneling splitting for n=4. A tentative explanation for this difference is proposed. The implications of our results for the interchange-tunneling dynamics of (HCl)2 in helium nanodroplets are discussed.

AB - We present a rigorous theoretical study of the solvation of (HCl)2 and (HF)2 by small (He4) n clusters, with n=1-14 and 30. Pairwise-additive potential-energy surfaces of Hen (HX)2 (X=Cl and F) clusters are constructed from highly accurate four-dimensional (rigid monomer) HX-HX and two-dimensional (rigid monomer) He-HX potentials and a one-dimensional He-He potential. The minimum-energy geometries of these clusters, for n=1-6 in the case of (HCl)2 and n=1-5 for (HF)2, correspond to the He atoms in a ring perpendicular to and bisecting the HX-HX axis. The quantum-mechanical ground-state energies and vibrationally averaged structures of Hen (HCl)2 (n=1-14 and 30) and Hen (HF)2 (n=1-10) clusters are calculated exactly using the diffusion Monte Carlo (DMC) method. In addition, the interchange-tunneling splittings of Hen (HCl)2 clusters with n=1-14 are determined using the fixed-node DMC approach, which was employed by us previously to calculate the tunneling splittings for Hen (HF)2 clusters, n=1-10 [A. Sarsa, Phys. Rev. Lett. 88, 123401 (2002)]. The vibrationally averaged structures of Hen (HX)2 clusters with n=1-6 for (HCl)2 and n=1-5 for (HF)2 have the helium density localized in an effectively one-dimensional ring, or doughnut, perpendicular to and at the midpoint of the HX-HX axis. The rigidity of the solvent ring varies with n and reaches its maximum for the cluster size at which the ring is filled, n=6 and n=5 for (HCl)2 and (HF)2, respectively. Once the equatorial ring is full, the helium density spreads along the HX-HX axis, eventually solvating the entire HX dimer. The interchange-tunneling splitting of Hen (HCl)2 clusters hardly varies at all over the cluster size range considered, n=1-14, and is virtually identical to that of the free HCl dimer. This absence of the solvent effect is in sharp contrast with our earlier results for Hen (HF)2 clusters, which show a ∼30% reduction of the tunneling splitting for n=4. A tentative explanation for this difference is proposed. The implications of our results for the interchange-tunneling dynamics of (HCl)2 in helium nanodroplets are discussed.

UR - http://www.scopus.com/inward/record.url?scp=29144519314&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=29144519314&partnerID=8YFLogxK

U2 - 10.1063/1.2136358

DO - 10.1063/1.2136358

M3 - Article

C2 - 16375482

AN - SCOPUS:29144519314

VL - 123

JO - Journal of Chemical Physics

JF - Journal of Chemical Physics

SN - 0021-9606

IS - 22

M1 - 224313

ER -