The vibrational spectra of krypton and xenon difluoride: High‐resolution infrared studies and ab initio calculations

The Near-Sightedness of Many-Body Interactions in Anharmonic Vibrational Couplings

Couplings between vibrational motions are driven by electronic interactions, and these couplings carry special significance in vibrational energy transfer, multidimensional spectroscopy experiments, and simulations of vibrational spectra. In this investigation, the many-body contributions to these couplings are analyzed computationally in the context of clathrate-like alkali metal cation hydrates, including Cs+(H2O)20, Rb+(H2O)20, and K+(H2O)20, using both analytic and quantum-chemistry potential energy surfaces. Although the harmonic spectra and one-dimensional anharmonic spectra depend strongly on these many-body interactions, the mode-pair couplings were, perhaps surprisingly, found to be dominated by one-body effects, even in cases of couplings to low-frequency modes that involved the motion of multiple water molecules. The origin of this effect was traced mainly to geometric distortion within water monomers and cancellation of many-body effects in differential couplings, and the effect was also shown to be agnostic to the identity of the ion. These outcomes provide new understanding of vibrational couplings and suggest the possibility of improved computational methods for the simulation of infrared and Raman spectra.

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.

Noble-gas compounds: A general procedure of bonding analysis

This paper accounts for a general procedure of bonding analysis that is, expectedly, adequate to describe any type of interaction involving the noble-gas (Ng) atoms. Building on our recently proposed classification of the Ng-X bonds (X = binding partner) [New J. Chem. 44, 15536 (2020)], these contacts are first distinguished into three types, namely, A, B, or C, based on the topology of the electron energy density H(r) and on the shape of its plotted form. Bonds of type B or C are, then, further assigned as B-loose (B^l) or B-tight (B^t) and C-loose (C^l) or C-tight (C^t) depending on the sign that H(r) takes along the Ng-X bond path located from the topological analysis of ρ(r), particularly at around the bond critical point (BCP). Any bond of type A, B^l/B^t, or C^l/C^t is, finally, assayed in terms of contribution of covalency. This is accomplished by studying the maximum, minimum, and average value of H(r) over the volume enclosed by the low-density reduced density gradient (RDG) isosurface associated with the bond (typically, the RDG isosurface including the BCP) and the average ρ(r) over the same volume. The bond assignment is also corroborated by calculating the values of quantitative indices specifically defined for the various types of interactions (A, B, or C). The generality of our taken approach should encourage its wide application to the study of Ng compounds.

Freezing in resonance structures for better packing: XeF2 becomes (XeF+)(F-) at large compression

Recent high-pressure experiments conducted on xenon difluoride (XeF(2)) suggested that this compound undergoes several phase transitions up to 100 GPa, becoming metallic above 70 GPa. In this theoretical study, in contrast to experiment, we find that the ambient pressure molecular structure of xenon difluoride, of I4/mmm symmetry, remains the most stable one up to 105 GPa. In our computations, the structures suggested from experiment have either much higher enthalpies than the I4/mmm structure or converge to that structure upon geometry optimization. We discuss these discrepancies between experiment and calculation and point to an alternative interpretation of the measured cell vectors of XeF(2) at high pressure. At pressures exceeding those studied experimentally, above 105 GPa, the I4/mmm structure transforms to one of Pnma symmetry. The Pnma phase contains bent FXeF molecules, with unequal Xe-F distances, and begins to bring other fluorines into the coordination sphere of the Xe. Further compression of this structure up to 200 GPa essentially results in self-dissociation of XeF(2) into an ionic solid (i.e., [XeF](+)F(-)), similar to what is observed for nitrous oxide (N(2)O) at high pressure.

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.

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).

Microscopic structure of a vortex line in a dilute superfluid fermi gas

The microscopic properties of a single vortex in a dilute superfluid Fermi gas at zero temperature are examined within the framework of self-consistent Bogoliubov-de Gennes theory. Using only physical parameters as input, we study the pair potential, the density, the energy, and the current distribution. Comparison of the numerical results with analytical expressions clearly indicates that the energy of the vortex is governed by the zero-temperature BCS coherence length.