Paul R. Horn

Unrestricted absolutely localized molecular orbitals for energy decomposition analysis: Theory and applications to intermolecular interactions involving radicals

Radical-closed shell and radical-radical intermolecular interactions are less well-understood than those between closed shell species. With the objective of gaining additional insight, this work reports a generalization of the absolutely localized molecular orbital (ALMO) energy decomposition analysis (EDA) to open shell fragments, described by self-consistent field methods, such as standard density functional theory. The ALMO-EDA variationally partitions an intermolecular interaction energy into three separate contributions; frozen orbital interactions, polarization, and charge transfer. The first examples involve comparison of the interactions of alkanes and alkyl radicals (methyl radical, methane, tertiary butyl radical, and isobutane) with sodium, potassium, hydronium, and ammonium cations. A second series of examples involve benzene cation interacting with a series of nucleophiles in both on-top and side-on geometries. The ALMO-EDA yields a variety of interesting insights into the relative roles of its component contributions as the interacting partners and their geometries are changed.

Useful lower limits to polarization contributions to intermolecular interactions using a minimal basis of localized orthogonal orbitals: Theory and analysis of the water dimer

The problem of describing the energy-lowering associated with polarization of interacting molecules is considered in the overlapping regime for self-consistent field wavefunctions. The existing approach of solving for absolutely localized molecular orbital (ALMO) coefficients that are block-diagonal in the fragments is shown based on formal grounds and practical calculations to often overestimate the strength of polarization effects. A new approach using a minimal basis of polarized orthogonal local MOs (polMOs) is developed as an alternative. The polMO basis is minimal in the sense that one polarization function is provided for each unpolarized orbital that is occupied; such an approach is exact in second-order perturbation theory. Based on formal grounds and practical calculations, the polMO approach is shown to underestimate the strength of polarization effects. In contrast to the ALMO method, however, the polMO approach yields results that are very stable to improvements in the underlying AO basis expansion. Combining the ALMO and polMO approaches allows an estimate of the range of energy-lowering due to polarization. Extensive numerical calculations on the water dimer using a large range of basis sets with Hartree-Fock theory and a variety of different density functionals illustrate the key considerations. Results are also presented for the polarization-dominated Na(+)CH4 complex. Implications for energy decomposition analysis of intermolecular interactions are discussed.

Dissociative Photoionization of Glycerol and its Dimer Occurs Predominantly via a Ternary Hydrogen-Bridged Ion–Molecule Complex

The photoionization and dissociative photoionization of glycerol are studied experimentally and theoretically. Time-of-flight mass spectrometry combined with vacuum ultraviolet synchrotron radiation ranging from 8 to 15 eV is used to investigate the nature of the major fragments and their corresponding appearance energies. Deuterium (1,1,2,3,3-D5) and (13)C (2-(13)C) labeling is employed to narrow down the possible dissociation mechanisms leading to the major fragment ions (C3H(x)O2(+), C2H(x)O2(+), C2H(x)O(+), CH(x)O(+)). We find that the primary fragmentation of the glycerol radical cation (m/z 92) occurs only via two routes. The first channel proceeds via a six-membered hydrogen-transfer transition state, leading to a common stable ternary intermediate, comprised of neutral water, neutral formaldehyde, and a vinyl alcohol radical cation, which exhibits a binding energy of ≈42 kcal/mol and a very short (1.4 Å) hydrogen bond. Fragmentation of this intermediate gives rise to experimentally observed m/z 74, 62, 44, and 45. Fragments m/z 74 and 62 both consist of hydrogen-bridged ion-molecule complexes with binding energy >25 kcal/mol, whereas the m/z 44 species lacks such stabilization. This explains why water- or formaldehyde-loss products are observed first. The second primary fragmentation route arises from cleaving the elongated C-C bond. Also for this channel, intermediates comprised of hydrogen-bridged ion-molecule complexes exhibiting binding energies >24 kcal/mol are observed. Energy decomposition analysis reveals that electrostatic and charge-transfer interactions are equally important in hydrogen-bridged ion-molecule complexes. Furthermore, the dissociative photoionization of the glycerol dimer is investigated and compared to the main pathways for the monomeric species. To a first approximation, the glycerol dimer radical cation can be described as a monomeric glycerol radical cation in the presence of a spectator glycerol, thus giving rise to a dissociation pattern similar to that of the monomer.

A systematic study on Pt based, subnanometer-sized alloy cluster catalysts for alkane dehydrogenation: effects of intermetallic interaction

Platinum-based bimetallic nanoparticles are analyzed by the application of density functional theory to a series of tetrahedral Pt3X cluster models, with element X taken from the P-block, preferably group 14, or from the D-block around group 10. Almost identical cluster geometries allow a systematic investigation of electronic effects induced by different elements X. Choosing the propane-to-propene conversion as the desired dehydrogenation reaction, we provide estimates for the activity and selectivity of the various catalysts based on transition state theory. No significant Brønsted-Evans-Polanyi-relation could be found for the given reaction. A new descriptor, derived from an energy decomposition analysis, captures the effect of element X on the rate-determining step of the first hydrogen abstraction. Higher activities than obtained for pure Pt4 clusters are predicted for Pt alloys containing Ir, Sn, Ge and Si, with Pt3Ir showing particularly high selectivity.

Approaching the basis set limit for DFT calculations using an environment-adapted minimal basis with perturbation theory: Formulation, proof of concept, and a pilot implementation

Recently developed density functionals have good accuracy for both thermochemistry (TC) and non-covalent interactions (NC) if very large atomic orbital basis sets are used. To approach the basis set limit with potentially lower computational cost, a new self-consistent field (SCF) scheme is presented that employs minimal adaptive basis (MAB) functions. The MAB functions are optimized on each atomic site by minimizing a surrogate function. High accuracy is obtained by applying a perturbative correction (PC) to the MAB calculation, similar to dual basis approaches. Compared to exact SCF results, using this MAB-SCF (PC) approach with the same large target basis set produces <0.15 kcal/mol root-mean-square deviations for most of the tested TC datasets, and <0.1 kcal/mol for most of the NC datasets. The performance of density functionals near the basis set limit can be even better reproduced. With further improvement to its implementation, MAB-SCF (PC) is a promising lower-cost substitute for conventional large-basis calculations as a method to approach the basis set limit of modern density functionals.