Structure and spectral features of H+(H2O)7: Eigen versus Zundel forms

The two dimensional (2D) to three dimensional (3D) transition for the protonated water cluster has been controversial, in particular, for H(+)(H(2)O)(7). For H(+)(H(2)O)(7) the 3D structure is predicted to be lower in energy than the 2D structure at most levels of theory without zero-point energy (ZPE) correction. On the other hand, with ZPE correction it is predicted to be either 2D or 3D depending on the calculational levels. Although the ZPE correction favors the 3D structure at the level of coupled cluster theory with singles, doubles, and perturbative triples excitations [CCSD(T)] using the aug-cc-pVDZ basis set, the result based on the anharmonic zero-point vibrational energy correction favors the 2D structure. Therefore, the authors investigated the energies based on the complete basis set limit scheme (which we devised in an unbiased way) at the resolution of the identity approximation Moller-Plesset second order perturbation theory and CCSD(T) levels, and found that the 2D structure has the lowest energy for H(+)(H(2)O)(7) [though nearly isoenergetic to the 3D structure for D(+)(D(2)O)(7)]. This structure has the Zundel-type configuration, but it shows the quantum probabilistic distribution including some of the Eigen-type configuration. The vibrational spectra of MP2/aug-cc-pVDZ calculations and Car-Parrinello molecular dynamics simulations, taking into account the thermal and dynamic effects, show that the 2D Zundel-type form is in good agreement with experiments.     

Structure, stability, thermodynamic properties, and infrared spectra of the protonated water octamer H(+)(H2O)8

We investigated various two-dimensional (2D) and three-dimensional (3D) structures of H (+)(H 2O) 8, using density functional theory (DFT), Moller-Plesset second-order perturbation theory (MP2), and coupled cluster theory with single, double, and perturbative triple excitations (CCSD(T)). The 3D structure is more stable than the 2D structure at all levels of theory on the Born-Oppenheimer surface. With the zero-point energy (ZPE) correction, the predicted structure varies depending on the level of theory. The DFT employing Becke's three parameters with Lee-Yang-Parr functionals (B3LYP) favors the 2D structure. At the complete basis set (CBS) limit, the MP2 calculation favors the 3D structure by 0.29 kcal/mol, and the CCSD(T) calculation favors the 3D structure by 0.27 kcal/mol. It is thus expected that both 2D and 3D structures are nearly isoenergetic near 0 K. At 100 K, all the calculations show that the 2D structure is much more stable in free binding energy than the 3D structure. The DFT and MP2 vibrational spectra of the 2D structure are consistent with the experimental spectra. First-principles Car-Parrinello molecular dynamics (CPMD) simulations show that the 2D Zundel-type vibrational spectra are in good agreement with the experiment.     

Accuracy and limitations of second-order many-body perturbation theory for predicting vertical detachment energies of solvated-electron clusters

Vertical electron detachment energies (VDEs) are calculated for a variety of (H(2)O)(n)(-) and (HF)(n)(-) isomers, using different electronic structure methodologies but focusing in particular on a comparison between second-order M?ller-Plesset perturbation theory (MP2) and coupled-cluster theory with noniterative triples, CCSD(T). For the surface-bound electrons that characterize small (H(2)O)(n)(-) clusters (n< or = 7), the correlation energy associated with the unpaired electron grows linearly as a function of the VDE but is unrelated to the number of monomers, n. In every example considered here, including strongly-bound "cavity" isomers of (H(2)O)(24)(-), the correlation energy associated with the unpaired electron is significantly smaller than that associated with typical valence electrons. As a result, the error in the MP2 detachment energy, as a fraction of the CCSD(T) value, approaches a limit of about -7% for (H(2)O)(n)(-) clusters with VDEs larger than about 0.4 eV. CCSD(T) detachment energies are bounded from below by MP2 values and from above by VDEs calculated using second-order many-body perturbation theory with molecular orbitals obtained from density functional theory. For a variety of both strongly- and weakly-bound isomers of (H(2)O)(20)(-) and (H(2)O)(24)(-), including both surface states and cavity states, these bounds afford typical error bars of +/-0.1 eV. We have found only one case where the Hartree-Fock and density functional orbitals differ qualitatively; in this case the aforementioned bounds lie 0.4 eV apart, and second-order perturbation theory may not be reliable.     

Ab initio electron correlated studies on the intracluster reaction of NO+ (H2O)(n) → H3O+ (H2O)(n-2) (HONO) (n = 4 and 5)

Hydrated nitrosonium ion clusters NO(+)(H(2)O)(n) (n = 4 and 5) were investigated by using MP2/aug-cc-pVTZ level of theory to clarify isomeric reaction pathways for formation of HONO and fully hydrated hydride ions. We found some new isomers and transition state structures in each hydration number, whose lowest activation energies of the intracluster reactions were found to be 4.1 and 3.4 kcal mol(-1) for n = 4 and n = 5, respectively. These thermodynamic properties and full quantum mechanical molecular dynamics simulation suggest that product isomers with HONO and fully hydrated hydride ions can be obtained at n = 4 and n = 5 in terms of excess hydration binding energies which can overcome these activation barriers.     

Structural rearrangements and magic numbers in reactions between pyridine-containing water clusters and ammonia

Molecular cluster ions H(+)(H(2)O)(n), H(+)(pyridine)(H(2)O)(n), H(+)(pyridine)(2)(H(2)O)(n), and H(+)(NH(3))(pyridine)(H(2)O)(n) (n = 16-27) and their reactions with ammonia have been studied experimentally using a quadrupole-time-of-flight mass spectrometer. Abundance spectra, evaporation spectra, and reaction branching ratios display magic numbers for H(+)(NH(3))(pyridine)(H(2)O)(n) and H(+)(NH(3))(pyridine)(2)(H(2)O)(n) at n = 18, 20, and 27. The reactions between H(+)(pyridine)(m)(H(2)O)(n) and ammonia all seem to involve intracluster proton transfer to ammonia, thus giving clusters of high stability as evident from the loss of several water molecules from the reacting cluster. The pattern of the observed magic numbers suggest that H(+)(NH(3))(pyridine)(H(2)O)(n) have structures consisting of a NH(4)(+)(H(2)O)(n) core with the pyridine molecule hydrogen-bonded to the surface of the core. This is consistent with the results of high-level ab initio calculations of small protonated pyridine/ammonia/water clusters.     

Exploring the rich energy landscape of sulfate-water clusters SO4(2-) (H2O)(n=3-7): an electronic structure approach

We present a reinvestigation of sulfate-water clusters SO4(2-) (H2O)(n=3-7), which involves several new aspects. Using a joint molecular mechanics/first principles approach, we perform exhaustive searches for stable cluster geometries, showing that the sulfate-water landscape is much richer than anticipated previously. We check the compatibility of the new structures with experiment by comparing vertical detachment energies (VDEs) calculated at the B3LYP/6-311++G** level of theory and determine the energetic ordering of the isomers at the RI-MP2/aug-cc-pVTZ level. Our results are bench-marked carefully against reference energies of estimated CCSD(T)/aug-cc-VTZ quality and VDEs of CCSD(T)/aug-cc-pVDZ quality. Furthermore, we calculate anharmonic vibrational corrections for up to the n = 6 clusters, which are shown to be significant for isomer energy ordering. We use energy decomposition analysis (EDA) based on the absolutely localized fragment (ALMO) expansion to gain chemical insight into the binding motifs.     

Theoretical Study of Nascent Hydration in the Fe+(H2O)n system

The interactions of the iron monocation with water molecules and argon atoms in the gas phase were studied computationally to elucidate recent infrared vibrational spectroscopy on this system. These calculations employ first-principles all-electron methods performed with B3LYP/DZVP density functional theory. The ground state of Fe(+)(H(2)O) is found to be a quartet (M = 2S + 1 = 4, S is the total spin). Different binding sites for the addition of one or two argon atoms produce several low-lying states of different geometry and multiplicity in a relatively small energy range for Fe(+)(H(2)O)-Ar(2) and Fe(+)(H(2)O)(2)-Ar. In both species, quartet states are lowest in energy, and sextets and doublets lie at higher energies from the respective ground states. These results are consistent with the conclusion that the experimentally determined infrared photodissociation spectra (IRPD) of Fe(+)(H(2)O)-Ar(2) and Fe(+)(H(2)O)(2)-Ar are complicated because of the presence of multiple isomeric structures. The estimated IR bands for the symmetric and asymmetric O-H stretches from different isomers provide new insight into the observed IRPD spectra.     

Development of a 3-body:many-body integrated fragmentation method for weakly bound clusters and application to water clusters (H2O)(n = 3-10, 16, 17)

A 3-body:many-body integrated quantum mechanical (QM) fragmentation method for non-covalent clusters is introduced within the ONIOM formalism. The technique captures all 1-, 2-, and 3-body interactions with a high-level electronic structure method, while a less demanding low-level method is employed to recover 4-body and higher-order interactions. When systematically applied to 40 low-lying (H(2)O)(n) isomers ranging in size from n = 3 to 10, the CCSD(T):MP2 3-body:many-body fragmentation scheme deviates from the full CCSD(T) interaction energy by no more than 0.07 kcal mol(-1) (or <0.01 kcal mol(-1) per water). The errors for this QM:QM method increase only slightly for various low-lying isomers of (H(2)O)(16) and (H(2)O)(17) (always within 0.13 kcal mol(-1) of the recently reported canonical CCSD(T)/aug-cc-pVTZ energies). The 3-body:many-body CCSD(T):MP2 procedure is also very efficient because the CCSD(T) computations only need to be performed on subsets of the cluster containing 1, 2, or 3 monomers, which in the current context means the largest CCSD(T) calculations are for 3 water molecules, regardless of the cluster size.     

On the contribution of vibrational anharmonicity to the binding energies of water clusters

The second-order vibrational perturbation theory method has been used together with the B3LYP and MP2 electronic structure methods to investigate the effects of anharmonicity on the vibrational zero-point energy (ZPE) contributions to the binding energies of (H2O)n, n = 2-6, clusters. For the low-lying isomers of (H2O)6, the anharmonicity correction to the binding energy is calculated to range from -248 to -355 cm(-1). It is also demonstrated that although high-order electron correlation effects are important for the individual vibrational frequencies, they are relatively unimportant for the net ZPE contributions to the binding energies of water clusters.     

Tuning of the internal energy and isomer distribution in small protonated water clusters H(+)(H2O)(4-8): an application of the inert gas messenger technique

Infrared spectroscopy of gas-phase hydrated clusters provides us much information on structures and dynamics of water networks. However, interpretation of spectra is often difficult because of high internal energy (vibrational temperature) of clusters and coexistence of many isomers. Here we report an approach to vary these factors by using the inert gas (so-called "messenger")-mediated cooling technique. Protonated water clusters with a messenger (M), H(+)(H(2)O)(4-8)·M (M = Ne, Ar, (H(2))(2)), are formed in a molecular beam and probed with infrared photodissociation spectroscopy in the OH stretch region. Observed spectra are compared with each other and with bare H(+)(H(2)O)(n). They show clear messenger dependence in their bandwidths and relative band intensities, reflecting different internal energy and isomer distribution, respectively. It is shown that the internal energy follows the order H(+)(H(2)O)(n) > H(+)(H(2)O)(n)·(H(2))(2) > H(+)(H(2)O)(n)·Ar > H(+)(H(2)O)(n)·Ne, while the isomer-selectivity, which changes the isomer distribution in the bare system, follows the order H(+)(H(2)O)(n)·Ar > H(+)(H(2)O)(n)·(H(2))(2) > H(+)(H(2)O)(n)·Ne ~ (H(+)(H(2)O)(n)). Although the origin of the isomer-selectivity is unclear, comparison among spectra measured with different messengers is very powerful in spectral analyses and makes it possible to easily assign spectral features of each isomer.     

Benchmark structures and binding energies of small water clusters with anharmonicity corrections

For (H(2)O)(n) where n = 1-10, we used a scheme combining molecular dynamics sampling with high level ab initio calculations to locate the global and many low lying local minima for each cluster. For each isomer, we extrapolated the RI-MP2 energies to their complete basis set limit, included a CCSD(T) correction using a smaller basis set and added finite temperature corrections within the rigid-rotor-harmonic-oscillator (RRHO) model using scaled and unscaled harmonic vibrational frequencies. The vibrational scaling factors were determined specifically for water clusters by comparing harmonic frequencies with VPT2 fundamental frequencies. We find the CCSD(T) correction to the RI-MP2 binding energy to be small (<1%) but still important in determining accurate conformational energies. Anharmonic corrections are found to be non-negligble; they do not alter the energetic ordering of isomers, but they do lower the free energies of formation of the water clusters by as much as 4 kcal/mol at 298.15 K.     

Nonstandard cages in the formation process of methane clathrate: stability, structure, and spectroscopic implications from first-principles

Endohedral CH(4)@(H(2)O)(n) (n = 16, 18, 20, 22, 24) clusters with standard and nonstandard cage configurations containing four-, five-, six-, seven-membered rings were generated by spiral algorithm and were systematically explored using DFT-D methods. The geometries of all isomers were optimized in vacuum and aqueous solution. In vacuum, encapsulation of methane molecules can stabilize the hollow (H(2)O)(n) cage by 2.31~5.44 kcal/mol; but the endohedral CH(4)@(H(2)O)(n) cages are still less stable than the pure (H(2)O)(n) clusters. Aqueous environment could promote the stabilities of the hollow (H(2)O)(n) cages as well as the CH(4)@(H(2)O)(n) clusters, and the CH(4)@(H(2)O)(n) clusters possess larger stabilization energies with regard to the pure (H(2)O)(n) clusters except for n = 24. The lowest energy structures of the CH(4)@(H(2)O)(20) and CH(4)@(H(2)O)(24) cages are identical to the building units in the crystalline sI clathrate hydrate. All of the low-energy cages (including both regular and irregular ones) have large structural similarity and can be connected by "dimer-insertion" operation and Stone-Wales transformation. Our calculation also showed that in the range of cluster size n = 16-24, the relative energies of cage isomers tend to decrease with increasing number of the adjacent pentagons in the oxygen skeleton structures. In addition to the regular endohedral CH(4)@(H(2)O)(20) and CH(4)@(H(2)O)(24) cage structures, some nonstandard CH(4)@(H(2)O)(n) (n = 18, 20, 22, 24) cages have lower energies and might appear during nucleation process of methane hydrate. For the methane molecules in these low-energy cage isomers, we found that the C-H symmetric stretching frequencies show a red-shift trend and the (13)C NMR chemical shifts generally move toward negative values as the cavity size increases. These theoretical results are comparable to the available experimental data and might help experimental identification of the endohedral water cages during nucleation.     

Sequential bond energies and barrier heights for the water loss and charge separation dissociation pathways of Cd(2+)(H2O)n, n = 3-11

The bond dissociation energies for losing one water from Cd(2+)(H(2)O)(n) complexes, n = 3-11, are measured using threshold collision-induced dissociation in a guided ion beam tandem mass spectrometer coupled with a thermal electrospray ionization source. Kinetic energy dependent cross sections are obtained for n = 4-11 complexes and analyzed to yield 0 K threshold measurements for loss of one, two, and three water ligands after accounting for multiple collisions, kinetic shifts, and energy distributions. The threshold measurements are converted from 0 to 298 K values to give the hydration enthalpies and free energies for sequentially losing one water from each complex. Theoretical geometry optimizations and single point energy calculations are performed on reactant and product complexes using several levels of theory and basis sets to obtain thermochemistry for comparison to experiment. The charge separation process, Cd(2+)(H(2)O)(n) → CdOH(+)(H(2)O)(m) + H(+)(H(2)O)(n-m-1), is also observed for n = 4 and 5 and the competition between this process and water loss is analyzed. Rate-limiting transition states for the charge separation process at n = 3-6 are calculated and compared to experimental threshold measurements resulting in the conclusion that the critical size for this dissociation pathway of hydrated cadmium is n(crit) = 4.     

Isotope exchange and structural rearrangements in reactions between size-selected ionic water clusters, H3O(+)(H2O)(n) and NH4(+)(H2O)(n), and D2O

Hydrogen/deuterium exchange in reactions of H3O(+)(H2O)n and NH4(+)(H2O)n (1 < or = n < or = 30) with D2O has been studied experimentally at center-of-mass collisions energies of < or = 0.2 eV. For a given cluster size, the cross-sections for H3O(+)(H2O)n and NH4(+)(H2O)n are similar, indicating a structural resemblance and energetics of binding. For protonated pure water clusters, H3O(+)(H2O)n, reacting with D2O the main H/D exchange mechanism is found to be proton catalyzed. In addition the H/D scrambling becomes close to statistically randomized for the larger clusters. For NH4(+)(H2O)n clusters reacting with D2O, the main mechanism is a D2O/H2O swap reaction. The lifetimes of H3O(+)(H2O)n clusters have been estimated using RRKM theory and a plateau in lifetime vs. cluster size is found already at n = 10.     

Hydration of the bisulfate ion: atmospheric implications

Using molecular dynamics configurational sampling combined with ab initio energy calculations, we determined the low energy isomers of the bisulfate hydrates. We calculated the CCSD(T) complete basis set (CBS) binding electronic and Gibbs free energies for 53 low energy isomers of HSO(4)(-)(H(2)O)(n=1-6) and derived the thermodynamics of adding waters sequentially to the bisulfate ion and its hydrates. Comparing the HSO(4)(-)/H(2)O system to the neutral H(2)SO(4)/H(2)O cluster, water binds more strongly to the anion than it does to the neutral molecules. The difference in the binding thermodynamics of HSO(4)(-)/H(2)O and H(2)SO(4)/H(2)O systems decreases with increasing number of waters. The thermodynamics for the formation of HSO(4)(-)(H(2)O)(n=1-5) is favorable at 298.15 K, and that of HSO(4)(-)(H(2)O)(n=1-6) is favorable for T < 273.15 K. The HSO(4)(-) ion is almost always hydrated at temperatures and relative humidity values encountered in the troposphere. Because the bisulfate ion binds more strongly to sulfuric acid than it does to water, it is expected to play a role in ion-induced nucleation by forming a strong complex with sulfuric acid and water, thus facilitating the formation of a critical nucleus.     

The effect of basis set superposition error on the convergence of intermolecular interaction energies for deprotonated complexes

The intermolecular interaction energies of the deprotonated hydrogen-bonded complexes F(-)(HF), F(-)(H(2)O), F(-)(NH(3)), Cl(-)(HF), SH(-)(HF), H(2)P(-)(HF), OH(-)(H(2)O), OH(-)(H(2)O)(2), OH(-)(NH(3)), Cl(-)(H(2)O), SH(-)(H(2)O), H(2)P(-)(H(2)O), Cl(-)(NH(3)), SH(-)(NH(3)), H(2)P(-)(NH(3)), Cl(-)(HCl), Cl(-)(H(2)S), Cl(-)(PH(3)), SH(-)(H(2)S), SH(-)(PH(3)), and H(2)P(-)(PH(3)) were calculated with correlation consistent basis sets at the MP2, MP4, QCISD(T), and CCSD(T) levels. When the basis set is smaller, the counterpoise-uncorrected intermolecular interaction energies are closer to the complete basis set limit than the counterpoise-corrected intermolecular interaction energies. The counterpoise-uncorrected intermolecular interaction energies obtained at the MP2/aug-cc-pVDZ level of theory are close to the interaction energies obtained at the extrapolated complete basis set limit in most of the complexes. Also, we investigate the accuracy of the other levels.     

Infrared spectroscopy of Cu+(H2O)(n) and Ag+(H2O)(n): coordination and solvation of noble-metal ions

M(+)(H(2)O)(n) and M(+)(H(2)O)(n)Ar ions (M=Cu and Ag) are studied for exploring coordination and solvation structures of noble-metal ions. These species are produced in a laser-vaporization cluster source and probed with infrared (IR) photodissociation spectroscopy in the OH-stretch region using a triple quadrupole mass spectrometer. Density functional theory calculations are also carried out for analyzing the experimental IR spectra. Partially resolved rotational structure observed in the spectrum of Ag(+)(H(2)O)(1) x Ar indicates that the complex is quasilinear in an Ar-Ag(+)-O configuration with the H atoms symmetrically displaced off axis. The spectra of the Ar-tagged M(+)(H(2)O)(2) are consistent with twofold coordination with a linear O-M(+)-O arrangement for these ions, which is stabilized by the s-d hybridization in M(+). Hydrogen bonding between H(2)O molecules is absent in Ag(+)(H(2)O)(3) x Ar but detected in Cu(+)(H(2)O)(3) x Ar through characteristic changes in the position and intensity of the OH-stretch transitions. The third H(2)O attaches directly to Ag(+) in a tricoordinated form, while it occupies a hydrogen-bonding site in the second shell of the dicoordinated Cu(+). The preference of the tricoordination is attributable to the inefficient 5s-4d hybridization in Ag(+), in contrast to the extensive 4s-3d hybridization in Cu(+) which retains the dicoordination. This is most likely because the s-d energy gap of Ag(+) is much larger than that of Cu(+). The fourth H(2)O occupies the second shells of the tricoordinated Ag(+) and the dicoordinated Cu(+), as extensive hydrogen bonding is observed in M(+)(H(2)O)(4) x Ar. Interestingly, the Ag(+)(H(2)O)(4) x Ar ions adopt not only the tricoordinated form but also the dicoordinated forms, which are absent in Ag(+)(H(2)O)(3) x Ar but revived at n=4. Size dependent variations in the spectra of Cu(+)(H(2)O)(n) for n=5-7 provide evidence for the completion of the second shell at n=6, where the dicoordinated Cu(+)(H(2)O)(2) subunit is surrounded by four H(2)O molecules. The gas-phase coordination number of Cu(+) is 2 and the resulting linearly coordinated structure acts as the core of further solvation processes.     

A photoelectron spectroscopy and ab initio study of B21-: negatively charged boron clusters continue to be planar at 21

The structures and chemical bonding of the B(21)(-) cluster have been investigated by a combined photoelectron spectroscopy and ab initio study. The photoelectron spectrum at 193 nm revealed a very high adiabatic electron binding energy of 4.38 eV for B(21)(-) and a congested spectral pattern. Extensive global minimum searches were conducted using two different methods, followed by high-level calculations of the low-lying isomers. The global minimum of B(21)(-) was found to be a quasiplanar structure with the next low-lying planar isomer only 1.9 kcal/mol higher in energy at the CCSD(T)/6-311-G* level of theory. The calculated vertical detachment energies for the two isomers were found to be in good agreement with the experimental spectrum, suggesting that they were both present experimentally and contributed to the observed spectrum. Chemical bonding analyses showed that both isomers consist of a 14-atom periphery, which is bonded by classical two-center two-electron bonds, and seven interior atoms in the planar structures. A localized two-center two-electron bond is found in the interior of the two planar isomers, in addition to delocalized multi-center σ and π bonds. The structures and the delocalized bonding of the two lowest lying isomers of B(21)(-) were found to be similar to those in the two lowest energy isomers in B(19)(-).     

Dissolution nature of cesium fluoride by water molecules

The structures, stabilities, thermodynamic quantities, dissociation energies, infrared spectra, and electronic properties of CsF hydrated by water molecules are investigated by using density functional theory, M?ller-Plesset second-order perturbation theory (MP2), coupled cluster theory with singles, doubles, and perturbative triples excitations (CCSD(T)), and ab initio molecular dynamic (AIMD) simulations. It is revealed that at 0 K three water molecules (as a global minimum structure) begin to half-dissociate the Cs-F, and six water molecules (though not a global minimum energy structure) can dissociate it. By the combination of the accurate CCSD(T) conformational energies for Cs(H2O)6 at 0 K with the AIMD thermal energy contribution, it reveals that the half-dissociated structure is the most stable at 0 K, but this structure (which is still the most stable) changes to the dissociated structure above 50 K. The spectra of CsF(H2O)(1-6) from MP2 calculations and the power spectra of CsF(H2O)6 from 50 and 100 K AIMD simulations are also reported.     

Heats of formation of xenon fluorides and the fluxionality of XeF(6) from high level electronic structure calculations

Atomization energies at 0 K and heats of formation at 0 and 298 K are predicted for XeF(+), XeF(-), XeF(2), XeF(4), XeF(5)(-), and XeF(6) from coupled cluster theory (CCSD(T)) calculations with new correlation-consistent basis sets for Xe. To achieve near chemical accuracy (+/-1 kcal/mol), up to four corrections were added to the complete basis set binding energies based on frozen core coupled cluster theory energies: a correction for core-valence effects, a correction for scalar relativistic effects, a correction for first-order atomic spin-orbit effects, and in some cases, a second-order spin-orbit correction. Vibrational zero-point energies were computed at the coupled cluster level of theory. The structure of XeF(6) is difficult to obtain with the C(3)(v)() and O(h)() structures having essentially the same energy. The O(h)() structure is only 0.19 kcal/mol below the C(3)(v)() one at the CCSD(T)/CBS level using an approximate geometry for the C(3)(v)() structure. With an optimized C(3)(v)() geometry, the C(3)(v)() structure would probably become slightly lower in energy than the O(h)() one. The calculated heats of formation for the neutral XeF(n)() fluorides are less negative than the experimental values from the equilibrium measurements by 2.0, 7.7, and 12.2 kcal/mol for n = 2, 4, and 6, respectively. For the experimental values, derived from the photoionization measurements, this discrepancy becomes even larger, suggesting a need for a redetermination of the experimental values. Evidence is presented for the fluxionality of XeF(6) caused by the presence of a sterically active, free valence electron pair on Xe.