MolSym: A Python package for handling symmetry in molecular quantum chemistry.

A consideration of the point group symmetry of molecules is often advantageous from a computational efficiency standpoint and sometimes necessary for the correct treatment of chemical physics problems. Many modern electronic structure software packages include a treatment of symmetry, but these are sometimes incomplete or unusable outside of that program's environment. Therefore, we have developed the MolSym package for handling molecular symmetry and its associated functionalities to provide a platform for including symmetry in the implementation and development of other methods. Features include point group detection, molecule symmetrization, arbitrary generation of symmetry element sets and character tables, and symmetry adapted linear combinations of real spherical harmonic basis functions, Cartesian displacement coordinates, and internal coordinates. We present some of the advantages of using molecular symmetry as achieved by MolSym, particularly with respect to Hartree-Fock theory, and the reduction of finite difference displacements in gradient/Hessian computations. This package is designed to be easily integrated into other software development efforts and may be extended to further symmetry applications.

Concordant Mode Approach for Molecular Vibrations.

The Concordant Mode Approach (CMA) is advanced as a novel hierarchy for increasing the system size and level of theory feasible for quantum chemical computations of harmonic vibrational frequencies. The key concept behind CMA is that transferrable, internal-coordinate normal modes computed at an appropriate lower level of theory (B) comprise a superb basis for converging to vibrational frequencies given by a higher level of theory (A). Accordingly, high-level harmonic frequencies can be evaluated via CMA from a collection of single-point energies that essentially scales linearly in the number of atoms, providing nearly order-of-magnitude CPU time speedups. The accuracy of CMA methods was established by comprehensive tests on over 120 molecules for target Level A = CCSD(T)/cc-pVTZ with auxiliary Level B choices of both CCSD(T)/cc-pVDZ and B3LYP/6-31G(2df,p). Remarkably, the frequency residuals given by the diagonal CMA-0A(nc) scheme exhibit mean absolute deviations (MADs) of only 0.2 cm-1 and standard deviations less than 0.5 cm-1; the corresponding zero-point vibrational energies (ZPVEs) have negligible errors in the vicinity of 0.3 cm-1.

Toward the Observation of the Tin and Lead Analogs of Formaldehyde.

Heavy aldehyde and ketone analogues, R2X═O (X = Si, Ge, Sn, or Pb), differ from their R2C═O counterparts due to their greater tendency to oligeramize as the X═O bond polarity increases as one goes down the periodic table. To date, H2Sn═O and H2Pb═O have eluded experimental detection. Herein we present the most rigorous theoretical study to date on these structures, providing CCSD(T)/pwCVTZ fundamental frequencies computed on CCSD(T)/CBS optimized structures for the H2X═O (X = Sn, Pb) potential energy surface. The focal point approach is employed to produce the CCSDTQ/CBS relative energies. For the Sn and Pb structures, the carbene-like cis-HXOH was the global minima, with the trans species being less than 0.6 and 1.1 kcal mol-1 above the cis structures, respectively. The formaldehyde-like H2X═O structure is in an energy well of at least 34.8 and 25.4 kcal mol-1 for Sn and Pb, respectively. Our results provide guidance for future work that may detect H2Sn═O or H2Pb═O for the first time.

The HOX⋯SO2 (X=F, Cl, Br, I) Binary Complexes: Implications for Atmospheric Chemistry.

Sulfur dioxide and hypohalous acids (HOX, X=F, Cl, Br, I) are ubiquitous molecules in the atmosphere that are central to important processes like seasonal ozone depletion, acid rain, and cloud nucleation. We present the first theoretical examination of the HOX⋯SO2 binary complexes and the associated trends due to halogen substitution. Reliable geometries were optimized at the CCSD(T)/aug-cc-pV(T+d)Z level of theory for HOF and HOCl complexes. The HOBr and HOI complexes were optimized at the CCSD(T)/aug-cc-pV(D+d)Z level of theory with the exception of the Br and I atoms which were modeled with an aug-cc-pwCVDZ-PP pseudopotential. 27 HOX⋯SO2 complexes were characterized and the focal point method was employed to produce CCSDT(Q)/CBS interaction energies. Natural Bond Orbital analysis and Symmetry Adapted Perturbation Theory were used to classify the nature of each principle interaction. The interaction energies of all HOX⋯SO2 complexes in this study ranged from 1.35 to 3.81 kcal mol-1 . The single-interaction hydrogen bonded complexes spanned a range of 2.62 to 3.07 kcal mol-1 , while the single-interaction halogen bonded complexes were far more sensitive to halogen substitution ranging from 1.35 to 3.06 kcal mol-1 , indicating that the two types of interactions are extremely competitive for heavier halogens. Our results provide insight into the interactions between HOX and SO2 which may guide further research of related systems.