On the Unusual Synclinal Conformations of Hexafluorobutadiene and Structurally Similar Molecules.

An explanation is presented for the unusual conformations of some molecules that contain the C═C-C═C core, namely, butadienes, biphenyls, and styrenes. Small substituents often induce a synclinal conformation, which brings the substituents into close proximity, and sometimes, there is no anticlinal minimum at all. This would not be predicted from steric repulsion arguments nor would it be expected that atoms that are nonbonded in a Lewis structure would approach closer than the sum of their van der Waals radii. Atomic energies calculated according to the quantum theory of atoms in molecules (QTAIM) do not show a consistent pattern for these structurally similar molecules, nor are intersubstituent bond paths consistently found, nor favorable diatomic interaction energies calculated using the interacting quantum atoms (IQA) scheme. Instead, the synclinal conformations are found to be driven by the attraction energy of the electron distribution of the carbon atoms and the nuclei of the molecule.

The barrier to the methyl rotation in Cis-2-butene and its isomerization energy to Trans-2-butene, revisited.

We respond to the two questions posed by Weinhold, Schleyer, and McKee (WSM) in their study of cis-2-butene (Weinhold et al., J Comput Chem 2014, 35, 1499), in which they solicit explanations for the relative conformational energies of this molecule in terms of the Quantum Theory of Atoms in Molecules (QTAIM). WSM requested answers to the questions: (1) why is cis-2-butene less stable than trans-2-butene despite the presence of a hydrogen-hydrogen (H⋯H) bond path in the former but not in the latter if the H⋯H bond path is stabilizing? (2) Why is the potential well of the conformational global minimum of cis-2-butene only 0.8 kcal/mol deep when the H⋯H bonding is stabilizing by 5 kcal/mol? Both questions raised by WSM are answered by considering the changes in the energies of all atoms as a function of the rotation of one of the two methyl groups from the minimum-energy structure, which exhibits the H⋯H bond path, to the transition state, which is devoid of this bond path. It is found that the stability gained by the H⋯H bonding interaction is cancelled by the destabilization of one of the ethylenic carbon atoms which, alone, destabilizes the system by as much as 5 kcal/mol in the global minimum conformation. Further, it is found that the 1.1 kcal/mol stability of trans-2-butene with respect to the cis-isomer is driven by the considerable destabilization of the ethylenic carbons by 11 kcal/mol, while the changes in the atomic energies of the other corresponding atoms in the two isomers account for the observed different stabilities. The error introduced into QTAIM atomic energies by neglecting the virials of the forces on the nuclei for partially optimized structures is discussed.

Aqueous coordination chemistry of H2: why is coordinated H2 inert to substitution by water in trans-Ru(P2)2(H2)H+-type complexes (P2 = a chelating phosphine)?

The reactivity of a series of trans-Ru(P(2))(2)Cl(2) complexes with H(2) was explored. The complexes reacted with H(2) via a stepwise H(2) addition/heterolysis pathway to form the trans-[Ru(P(2))(2)(H(2))H](+) dihydrogen complexes. Some of the resulting eta(2)-H(2) complexes were surprisingly inert to substitution by water, even at concentrations as high as 55 M; however, the identity of the bidentate phosphine ligand greatly influenced the lability of the coordinated eta(2)-H(2) ligand. With less electron-donating phosphine ligands, the H(2) ligand was susceptible to substitution by H(2)O, whereas with more electron-rich phosphine ligands, the H(2) ligand was inert to substitution by water. Density functional theory (DFT) calculations of the ligand substitution reactions showed that the Ru-H(2) and Ru-H(2)O complexes are very close in energy, and therefore slight changes in the donor properties of the bidentate phosphine ligand can inhibit or promote the substitution of H(2)O for H(2).

Pseudospectral time-dependent density functional theory.

Time-dependent density functional theory (TDDFT) is implemented within the Tamm-Dancoff approximation (TDA) using a pseudospectral approach to evaluate two-electron repulsion integrals. The pseudospectral approximation uses a split representation with both spectral basis functions and a physical space grid to achieve a reduction in the scaling behavior of electronic structure methods. We demonstrate here that exceptionally sparse grids may be used in the excitation energy calculation, following earlier work employing the pseudospectral approximation for determining correlation energies in wavefunction-based methods with similar conclusions. The pseudospectral TDA-TDDFT method is shown to be up to ten times faster than a conventional algorithm for hybrid functionals without sacrificing chemical accuracy.

Spectroscopic and ab initio characterization of the [ReH9]2- ion.

The dynamics and bonding of the hydrido complex Ba[ReH9], containing the D3h face-capped trigonal prismatic [ReH9]2- ion, have been investigated by vibrational spectroscopy and density functional theory (DFT). The combination of infrared, Raman, and inelastic neutron-scattering (INS) spectroscopies has enabled observation of all the modes of the [ReH9]2- ion for the first time. We demonstrate that calculations of the isolated [ReH9]2- ion are unable to reproduce the INS spectrum and that the complete unit cell must be considered with periodic DFT to have reliable results. This is shown to be a consequence of the long-range Coulomb potential present. Analysis of the electronic structure shows that the bonding between the rhenium and the hydrogen is largely covalent. There is a small degree of covalency between the prism hydrides and the barium. The counterion is crucial to the stability of the materials; hence, variation of it potentially offers a method to fine-tune the properties of the material.

Nuclear-magnetic-resonance shielding constants calculated by pseudospectral methods.

We have developed an algorithm based upon pseudospectral (PS) ab initio electronic structure methods for evaluating nuclear magnetic shielding constants using gauge-including atomic orbitals (GIAOs) in the spin-restricted and spin-unrestricted formalisms of Hartree-Fock (HF) theory and density-functional theory (DFT). The nuclear magnetic shielding constants for both 1H and 13C calculated using PS methodology for 21 small molecules have absolute mean errors of less than 0.3 ppm in comparison with analytic integral results. CPU timing comparisons between PS methods and conventional methods carried out for seven large molecules ranging from 510 to 1285 basis functions demonstrate that the PS methods are an order of magnitude more efficient than the conventional methods. PS-HF was between 9 and 26 times faster than conventional integral technology, and PS-DFT (Becke three-parameter Lee-Yang-Parr) was between 6 and 21 times faster.

Pseudospectral Local Second-Order Møller-Plesset Methods for Computation of Hydrogen Bonding Energies of Molecular Pairs.

We present a methodology for computing the binding energy of molecular dimers based on extrapolation of pseudospectral local second-order Moller-Plesset (MP2), or PS-LMP2, energies to the basis set limit. The extrapolation protocol is based on carrying out PS-LMP2 calculations with the Dunning cc-pVTZ (-f) and cc-pVQZ (-g) basis sets and then using a simple two-parameter function to compute the final basis set limit results. The function is parametrized to ultralarge basis set MP2 calculations for 5 molecular pairs taken from the literature and then tested by calculating results for a set of formamide dimers for which such calculations have also been carried out. The results agree to within ca. 0.2 kcal/mol with the conventional MP2 large basis set calculations. A specialized, but relatively simple, protocol is described for eliminating noise due to overcompleteness of the basis set. Timing results are presented for the LMP2 calculations, and comparisons are made with the LMP2 methodology of the QChem program. CPU time required by each of the methods scales as N(3), where N is the number of the basis functions, with the PS-LMP2 approach displaying a 2- to 3-fold advantage in the prefactor. We also discuss one set of test cases for which the PS-LMP2 results disagree with those obtained from an alternative type of MP2 calculation, N-methyl acetamide (NMA) dimers, and show that the results for liquid-state simulations using polarizable parameters derived by fitting to the PS-LMP2 binding energies appear to produce better results when compared with experimental data. The convergence issues associated with the alternative MP2 formulation remain to be investigated.

Photochemical studies as a function of solvent viscosity. A new photochemical pathway in the reaction of (eta5-C5H4Me)2Mo2(CO)6 with CCl4.

This paper reports the results of a study that used systematic changes in the solvent viscosity to probe the photochemical reactivity of the Cp'2Mo2(CO)6 (Cp' = eta5-C5H4CH3) molecule. The quantum yields for photolysis of Cp'2Mo2(CO)6 in the presence of CCl4 were studied as a function of solvent viscosity. The quantum yields did not decrease to zero with increasing solvent viscosity but rather leveled off at a constant, non-zero value. This result cannot be explained by any of the previously reported radical or Mo-CO dissociation photochemical pathways for this molecule, and therefore an additional photochemical pathway is suggested to be operating in the reaction. The new pathway may involve a reactive isomer of Cp'2Mo2(CO)6 or possibly electron transfer between the excited state of Cp'2Mo2(CO)6 and CCl4.