Abinitiocalculations and high resolution infrared investigation on XeF4

Accelerating Anharmonic Spectroscopy Simulations via Local-Mode, Multilevel Methods

Ab initio computer simulations of anharmonic vibrational spectra provide nuanced insight into the vibrational behavior of molecules and complexes. The computational bottleneck in such simulations, particularly for ab initio potentials, is often the generation of mode-coupling potentials. Focusing specifically on two-mode couplings in this analysis, the combination of a local-mode representation and multilevel methods is demonstrated to be particularly symbiotic. In this approach, a low-level quantum chemistry method is employed to predict the pairwise couplings that should be included at the target level of theory in vibrational self-consistent field (and similar) calculations. Pairs that are excluded by this approach are "recycled" at the low level of theory. Furthermore, because this low-level pre-screening will eventually become the computational bottleneck for sufficiently large chemical systems, distance-based truncation is applied to these low-level predictions without substantive loss of accuracy. This combination is demonstrated to yield sub-wavenumber fidelity with reference vibrational transitions when including only a small fraction of target-level couplings; the overhead of predicting these couplings, particularly when employing distance-based, local-mode cutoffs, is a trivial added cost. This combined approach is assessed on a series of test cases, including ethylene, hexatriene, and the alanine dipeptide. Vibrational self-consistent field (VSCF) spectra were obtained with an RI-MP2/cc-pVTZ potential for the dipeptide, at approximately a 5-fold reduction in computational cost. Considerable optimism for increased accelerations for larger systems and higher-order couplings is also justified, based on this investigation.

Heats of formation of XeF(3)(+), XeF(3)(-), XeF(5)(+), XeF(7)(+), XeF(7)(-), and XeF(8) from high level electronic structure calculations

Atomization energies at 0 K and heats of formation at 0 and 298 K are predicted for XeF(3)(+), XeF(3)(-), XeF(5)(+), XeF(7)(+), XeF(7)(-), and XeF(8) from coupled cluster theory (CCSD(T)) calculations with effective core potential correlation-consistent basis sets for Xe and including correlation of the nearest core electrons. Additional corrections are included to achieve near chemical accuracy of +/-1 kcal/mol. Vibrational zero point energies were computed at the MP2 level of theory. Unlike the other neutral xenon fluorides, XeF(8) is predicted to be thermodynamically unstable with respect to loss of F(2) with the reaction calculated to be exothermic by 22.3 kcal/mol at 0 K. XeF(7)(+) is also predicted to be thermodynamically unstable with respect to the loss of F(2) by 24.1 kcal/mol at 0 K. For XeF(3)(+), XeF(5)(+), XeF(3)(-), XeF(5)(-), and XeF(7)(-), the reactions for loss of F(2) are endothermic by 14.8, 37.8, 38.2, 59.6, and 31.9 kcal/mol at 0 K, respectively. The F(+) affinities of Xe, XeF(2), XeF(4), and XeF(6) are predicted to be 165.1, 155.3, 172.7, and 132.5 kcal/mol, and the corresponding F(-) affinities are 6.3, 19.9, 59.1, and 75.0 kcal/mol at 0 K, respectively.

Molecular single-bond covalent radii for elements 1-118

A self-consistent system of additive covalent radii, R(AB)=r(A) + r(B), is set up for the entire periodic table, Groups 1-18, Z=1-118. The primary bond lengths, R, are taken from experimental or theoretical data corresponding to chosen group valencies. All r(E) values are obtained from the same fit. Both E-E, E-H, and E-CH(3) data are incorporated for most elements, E. Many E-E' data inside the same group are included. For the late main groups, the system is close to that of Pauling. For other elements it is close to the methyl-based one of Suresh and Koga [J. Phys. Chem. A 2001, 105, 5940] and its predecessors. For the diatomic alkalis MM' and halides XX', separate fits give a very high accuracy. These primary data are then absorbed with the rest. The most notable exclusion are the transition-metal halides and chalcogenides which are regarded as partial multiple bonds. Other anomalies include H(2) and F(2). The standard deviation for the 410 included data points is 2.8 pm.

Xe129 chemical shift by the perturbational relativistic method: Xenon fluorides

(129)Xe nuclear shielding tensor is calculated at the leading-order, one-electron Breit-Pauli perturbation theory (BPPT) level for the xenon fluorides XeF(+), XeF(2), XeF(3) (+), and XeF(4) that cover the large nuclear magnetic resonance chemical shift range of this nucleus. BPPT is found to improve the shift range and relative shifts as compared to the nonrelativistic (NR) theory. While the full BPPT expansion consists of 16 relativistic terms, 5 of them are responsible for the entire chemical shift and shielding anisotropy. The remaining terms are practically isotropic, corelike contributions that are significant for the absolute shielding constant but cancel for the relative chemical shifts. The five principal terms are due to the spin-orbit-modified wave function allowing the Fermi contact and spin-dipole hyperfine interactions to be coupled to the orbital Zeeman interaction, as well as three distinct scalar relativistic modifications of the NR paramagnetic shielding: wave function change due to mass-velocity and Darwin interactions and the relativistic modification of the orbital hyperfine interaction. A very good agreement with the experimental shifts is obtained for XeF(2) and the particularly challenging XeF(+) species when both the NR and the five main relativistic terms are calculated at electron-correlated ab initio levels of theory. The performance of density-functional theory (DFT) with different pure and hybrid exchange-correlation functionals (with increasing exact exchange admixture) is tested against the ab initio data for each individual contribution. It is shown that DFT has difficulties in the description of paramagnetic shielding, already and especially in the NR contribution, which causes a large discrepancy of DFT results with experiment for xenon fluorides. In contrast, the DFT errors for the relativistic terms cancel out to the extent that a fairly good approximation of the total relativistic shift and anisotropy contributions may be obtained. A combination of high-level ab initio NR calculation with hybrid DFT estimates of the five main BPPT terms is proposed for reasonable estimates of xenon chemical shift in molecules. For the difficult cases such as the present XeF(+) and XeF(3) (+) cations, correlated ab initio calculations are unavoidable throughout. None of the other currently available relativistic methods, either at the fully relativistic or a variationally stable quasirelativistic levels of theory, surpasses the quality of the present approach for Xe shifts in these systems.

Accurate calculation of anharmonic vibrational frequencies of medium sized molecules using local coupled cluster methods

Local coupled cluster methods were applied for the automated generation of accurate multidimensional potential energy surfaces for a set of test molecules ranging from six to nine atoms. Based on these surfaces anharmonic fundamental frequencies were computed using vibrational self-consistent field and configuration interaction methods. The computed vibrational frequencies are compared to those obtained from similar calculations using conventional coupled cluster methods and to experimental values. The results from local and conventional methods are found to be of similar accuracy and in close agreement with experimental values. In addition, an efficient parallelization of the fully automated surface generation code is presented.

Vibrational analyses for CHFClBr and CDFClBr based on high levelab initiocalculations

Anharmonicity corrections to the harmonic vibrational spectra of CHFClBr and its deuterated isotopomer were computed by means of variational and perturbational approaches. A comparison of both methods is provided. Based on CCSD(T)/aug-cc-pVTZ electronic structure calculations excellent agreement with experimental data was obtained. Absolute mean deviations are in the range of about 4 cm(-1) for the fundamental modes, while slightly larger values of about 7 cm(-1) were found for the first vibrational overtones. In addition, vibrationally averaged structural parameters are provided for both molecules. The calculations will serve as a future starting point for parity-violation effects in vibrational transitions in these chiral molecules.

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.

Calculations of Vibrational Energy Levels by Using a Hybrid ab Initio and DFT Quartic Force Field: Application to Acetonitrile

A hybrid quartic force field with quadratic force constants calculated at the CCSD(T)/cc-pVTZ level and cubic and quartic terms determined by a B3LYP/cc-pVTZ treatment is proposed to compute the vibrational energy levels of acetonitrile from a variational method. Fundamentals and overtones calculated in the range of 300-3200 cm(-1) are in fair agreement with the 31 observed data, with an absolute mean deviation of less than 0.8%. These results allow us to explain several Fermi resonances.

Efficient calculation of potential energy surfaces for the generation of vibrational wave functions

An automatic procedure for the generation of potential energy surfaces based on high level ab initio calculations is described. It allows us to determine the vibrational wave functions for molecules of up to ten atoms. Speedups in computer time of about four orders of magnitude in comparison to standard implementations were achieved. Effects due to introduced approximations--within the computation of the potential--on fundamental modes obtained from vibrational self-consistent field and vibrational configuration interaction calculations are discussed. Benchmark calculations are provided for formaldehyde and 1,2,5-oxadiazole (furazan).