Voronoi deformation density (VDD) charges: Assessment of the Mulliken, Bader, Hirshfeld, Weinhold, and VDD methods for charge analysis

We present the Voronoi Deformation Density (VDD) method for computing atomic charges. The VDD method does not explicitly use the basis functions but calculates the amount of electronic density that flows to or from a certain atom due to bond formation by spatial integration of the deformation density over the atomic Voronoi cell. We compare our method to the well-known Mulliken, Hirshfeld, Bader, and Weinhold [Natural Population Analysis (NPA)] charges for a variety of biological, organic, and inorganic molecules. The Mulliken charges are (again) shown to be useless due to heavy basis set dependency, and the Bader charges (and often also the NPA charges) are not realistic, yielding too extreme values that suggest much ionic character even in the case of covalent bonds. The Hirshfeld and VDD charges, which prove to be numerically very similar, are to be recommended because they yield chemically meaningful charges. We stress the need to use spatial integration over an atomic domain to get rid of basis set dependency, and the need to integrate the deformation density in order to obtain a realistic picture of the charge rearrangement upon bonding. An asset of the VDD charges is the transparency of the approach owing to the simple geometric partitioning of space. The deformation density based charges prove to conform to chemical experience.     

Complex Reaction Environments and Competing Reaction Mechanisms in Zeolite Catalysis: Insights from Advanced Molecular Dynamics

The methanol‐to‐olefin process is a showcase example of complex zeolite‐catalyzed chemistry. At real operating conditions, many factors affect the reactivity, such as framework flexibility, adsorption of various guest molecules, and competitive reaction pathways. In this study, the strength of first principle molecular dynamics techniques to capture this complexity is shown by means of two case studies. Firstly, the adsorption behavior of methanol and water in H‐SAPO‐34 at 350 °C is investigated. Hereby an important degree of framework flexibility and proton mobility was observed. Secondly, the methylation of benzene by methanol through a competitive direct and stepwise pathway in the AFI topology was studied. Both case studies clearly show that a first‐principle molecular dynamics approach enables unprecedented insights into zeolite‐catalyzed reactions at the nanometer scale to be obtained.     

Electronic spectrum of UO2(2+) and [UO2Cl4]2- calculated with time-dependent density functional theory

The electronic spectra of UO(2) (2+) and [UO(2)Cl(4)](2-) are calculated with a recently proposed relativistic time-dependent density functional theory method based on the two-component zeroth-order regular approximation for the inclusion of spin-orbit coupling and a noncollinear exchange-correlation functional. All excitations out of the bonding sigma(u) (+) orbital into the nonbonding delta(u) or phi(u) orbitals for UO(2) (2+) and the corresponding excitations for [UO(2)Cl(4)](2-) are considered. Scalar relativistic vertical excitation energies are compared to values from previous calculations with the CASPT2 method. Two-component adiabatic excitation energies, U-O equilibrium distances, and symmetric stretching frequencies are compared to CASPT2 and combined configuration-interaction and spin-orbit coupling results, as well as to experimental data. The composition of the excited states in terms of the spin-orbit free states is analyzed. The results point to a significant effect of the chlorine ligands on the electronic spectrum, thereby confirming the CASPT2 results: The excitation energies are shifted and a different luminescent state is found.     

Realistic modeling of ruthenium-catalyzed transfer hydrogenation

We report the first computational study of a fully atomistic model of the ruthenium-catalyzed transfer hydrogenation of formaldehyde and the reverse reaction in an explicit methanol solution. Using ab initio molecular dynamics techniques, we determined the thermodynamics, mechanism, and electronic structure along the reaction path. To assess the effect of the solvent quantitatively, we make a direct comparison with the gas-phase reaction. We find that the energy profile in solution bears little resemblance to the profile in the gas phase and a distinct solvation barrier is found: the activation barriers in both directions are lowered and the concerted hydride and proton transfer in the gas phase are converted into a sequential mechanism in solution with the substrate appearing as methoxide-like intermediate. Our results indicate that besides the metal-ligand bifunctional mechanism, as proposed by Noyori, also a concerted solvent-mediated mechanism is feasible. Our study gives a new perspective of the active role a solvent can have in transition-metal-catalyzed reactions.     

Polyoxometalates with internal cavities: redox activity, basicity, and cation encapsulation in [Xn+P5W30O110](15-n)- Preyssler complexes, with X = Na+, Ca2+, Y3+, La3+, Ce3+, and Th4+

The Preyssler anion, of general formula [Xn+P5W30O110](15-n)-, is the smallest polyoxometalate (POM) with an internal cavity allowing cation exchange with the solution. The Preyssler anion features a rich chemistry evidenced by its ability to accept electrons at low potentials, to selectively capture various metal cations, and to undergo acid-base reactions. A deep understanding of these topics is herein provided by means of DFT calculations on the title series of compounds. We evaluate the energetics of the release/encapsulation process for several Xn+ cations and identify the effect of the encapsulated ion on the properties of the Preyssler anion. We revisit the relationship between the internal cation charge and the electrochemical behavior of the POM. A linear dependence between the first one-electron reduction energies and the encapsulated Xn+ charge is found, with a slope of 48 mV per unit charge. The protonation also shifts the reduction potential to more positive values, but the effect is much larger. In connection to this, the last proton's pKa = 2 for the Na+ derivative was estimated to be in reasonable agreement with experiment. The electronic structure of lanthanide derivatives is more complex than conventional POM structures. The reduction energy for the CeIV-Preyssler + 1e- --> CeIII-Preyssler process was computed to be more exothermic than that of very oxidant species such as S2Mo18O624-.     

Oxidative properties of FeO2+: electronic structure and solvation effects

An electronic structure analysis is provided of the action of solvated FeO(2+), [FeO(H(2)O)(5)](2+), as a hydroxylation catalyst. It is emphasized that the oxo end of FeO(2+) does not form hydrogen bonds (as electron donor and H-bond acceptor) with H-bond donors nor with aliphatic C-H bonds, but it activates C-H bonds as an electron acceptor. It is extremely electrophilic, to the extent that it can activate even such poor electron donors as aliphatic C-H bonds, the C-H bond orbital acting as electron donor in a charge transfer type of interaction. Lower lying O-H bonding orbitals are less easily activated. The primary electron accepting orbital in a water environment is the 3sigma*alpha orbital, an antibonding combination of Fe-3d(z(2)) and O-2p(z), which is very low-lying relative to the pi*alpha compared with, for example, the sigma* orbital in O(2) relative to its pi*. This is ascribed to relatively small Fe-3d(z(2)) with O-2p(z) overlap, due to the nodal structure of the 3d(z(2)).The H-abstraction barrier is very low in the gas phase, but it is considerably enhanced in water solvent. This is shown to be due to strong screening effects of the dielectric medium, leading to relative destabilization of the levels of the charged [FeO(H(2)O)(5)](2+) species compared to those of the neutral substrate molecules, making it a less effective electron acceptor. The solvent directly affects the orbital interactions responsible for the catalytic reaction.     

Gas phase rate coefficients and activation energies for the reaction of butanal and 2‐methyl‐propanal with nitrate radicals

Rate coefficients have been determined for the reaction of butanal and 2‐methyl‐propanal with NO₃ using relative and absolute methods. The relative measurements were accomplished by using a static reactor with long‐path FTIR spectroscopy as the analytical tool. The absolute measurements were made using fast‐flow–discharge technique with detection of NO₃ by optical absorption. The resulting average coefficients from the relative rate experiments were k = (1.0 ± 0.1) × 10⁻¹⁴ and k = (1.2 ± 0.2) × 10⁻¹⁴ (cm³ molecule⁻¹ s⁻¹) for butanal and 2‐methyl‐propanal, respectively. The results from the absolute measurements indicated secondary reactions involving NO₃ radicals and the primary formed acyl radicals. The prospect of secondary reactions was investigated by means of mathematical modeling. Calculations indicated that the unwanted NO₃ radical reactions could be suppressed by introducing molecular oxygen into the flow tube. The rate coefficients from the absolute rate experiments with oxygen added were and k = (1.2 ± 0.1) × 10⁻¹⁴ and = (0.9 ± 0.1) × 10⁻¹⁴ (cm³ molecule⁻¹ s⁻¹) for butanal and 2‐methyl‐propanal. The temperature dependence of the reactions was studied in the range between 263 and 364 K. Activation energies for the reactions were determined to 12 ± 2 kJ mole⁻¹ and 14 ± 1 kJ mole⁻¹ for butanal and 2‐methyl‐propanal, respectively. © 2000 John Wiley & Sons, Inc. Int J Chem Kinet 32: 294–303, 2000     

Hydroxylation catalysis by mononuclear and dinuclear iron oxo catalysts: a methane monooxygenase model system versus the Fenton reagent Fe(IV)O(H2O)5(2+)

Hydroxylation of aliphatic C-H bonds is a chemically and biologically important reaction, which is catalyzed by the oxidoiron group FeO(2+) in both mononuclear (heme and nonheme) and dinuclear complexes. We investigate the similarities and dissimilarities of the action of the FeO(2+) group in these two configurations, using the Fenton-type reagent [FeO(2+) in a water solution, FeO(H(2)O)(5)(2+)] and a model system for the methane monooxygenase (MMO) enzyme as representatives. The high-valent iron oxo intermediate MMOH(Q) (compound Q) is regarded as the active species in methane oxidation. We show that the electronic structure of compound Q can be understood as a dimer of two Fe(IV)O(2+) units. This implies that the insights from the past years in the oxidative action of this ubiquitous moiety in oxidation catalysis can be applied immediately to MMOH(Q). Electronically the dinuclear system is not fundamentally different from the mononuclear system. However, there is an important difference of MMOH(Q) from FeO(H(2)O)(5)(2+): the largest contribution to the transition state (TS) barrier in the case of MMOH(Q) is not the activation strain (which is in this case the energy for the C-H bond lengthening to the TS value), but it is the steric hindrance of the incoming CH(4) with the ligands representing glutamate residues. The importance of the steric factor in the dinuclear system suggests that it may be exploited, through variation in the ligand framework, to build a synthetic oxidation catalyst with the desired selectivity for the methane substrate.     

Assessment of density functional methods for reaction energetics: Iridium‐catalyzed water oxidation as case study

We investigate basis set convergence for a series of density functional theory (DFT) functionals (both hybrid and nonhybrid) and compare to coupled‐cluster with single and double excitations and perturbative triples [CCSD(T)] benchmark calculations. The case studied is the energetics of the water oxidation reaction by an iridium‐oxo complex. Complexation energies for the reactants and products complexes as well as the transition state (TS) energy are considered. Contrary to the expectation of relatively weak basis set dependence for DFT, the basis set effects are large, for example, more than 10 kcal mol^?1 difference from converged basis for the activation energy with “small” basis sets (DZ/6‐31G** for Ir/other atoms, or SVP) and still more than 6 kcal mol^?1 for def2‐TZVPP/6‐31G**. Inclusion of the dispersion correction in DFT‐D3 schemes affects the energies of reactant complex (RC), TS, and product complex (PC) by almost the same amount; it significantly improves the complexation energy (the formation of RC), but has little effect on the activation energy with respect to RC. With converged basis, some pure GGAs (PBE‐D3, BP86‐D3) as well as the hybrid functional B3LYP‐D3 are very accurate compared to benchmark CCSD(T) calculations. © 2012 Wiley Periodicals, Inc.     

Growth of sputter-deposited gold nanoparticles in ionic liquids

The growth of gold nanoparticles (NPs) synthesized by sputter deposition on an ionic liquid surface is studied in situ in the bulk phase of the ionic liquids (ILs) 1-butyl-3-methylimidazolium dicyanamide [C(1)C(4)Im][N(CN)(2)], 1-butyl-3-methylimidazolium bis(trifluoromethanesulfonyl)amide [C(1)C(4)Im][Tf(2)N], 1-butyl-3-methylimidazolium tetrafluoroborate [C(1)C(4)Im][BF(4)], 1-butyl-3-methylimidazolium hexafluorophosphate [C(1)C(4)Im][PF(6)] and 1-butyl-3-methylimidazolium triflate [C(1)C(4)Im][TfO]. It is found that primary nanoparticles with a diameter smaller than 2.5 nm are present in the sample immediately after sputtering. Growth of these primary particles proceeds after the end of the sputtering process and stops when the nanoparticles reach a certain size. Depending on the viscosity of the ionic liquid this growth process can proceed several hours to several days. The growth speed is fastest for the least viscous ionic liquid and follows the trend [C(1)C(4)Im][N(CN)(2)] > [C(1)C(4)Im][Tf(2)N] > [C(1)C(4)Im][TfO] > [C(1)C(4)Im][BF(4)] > [C(1)C(4)Im][PF(6)]. It is also found that a higher concentration of sputtered gold results in faster growth of the gold nanoparticles. A discussion on the growth mechanism of sputtered gold NPs is included.