Open-shell Density Functional Theory Research Articles

Observation of Single Molecule Plasmon-Driven Electron Transfer in Isotopically Edited 4,4'-Bipyridine Gold Nanosphere Oligomers.

We clarify mechanistic questions regarding plasmon-driven chemistry and nanoscale photocatalysis within optically confined near-field plasmonic systems. Using surface-enhanced Raman scattering (SERS), we directly monitor the photoinduced reaction dynamics of 4,4'-bipyridine molecules, localized in plasmonic hot spots within individual gold nanosphere oligomers. Our experiment generates surface electrons from the gold nanoparticle using an intense off-molecular resonance continuous wave pump field, and detects radical anion products via SERS. This is done by adopting a dual-wavelength spectroscopic approach. Empirical evidence of plasmon-driven electron transfer is provided for the first time by direct detection of the 4,4'-bipyridine radical anion species localized in the plasmonic hot spots of individual gold nanosphere oligomers, corroborated by open-shell density functional theory calculations. An isotopologue approach using both protonated and deuterated 4,4'-bipyridine molecules demonstrates the single molecule response of plasmon-driven electron transfer occurring in single nanosphere oligomer systems with a 3% yield, a phenomenon unobserved in ensemble measurements under analogous experimental conditions. This mechanism has broad applicability to using nanoscale chemical reactors for surface redox reactions on the subnanometer scale.

Real-time TD-DFT study on the dioxygen/superoxide radical charge transfer reaction

Superoxide (O2·−) is an oxygen radical anion that damages biochemical systems. Here I executed real-time computations on the dioxygen/superoxide charge transfer (CT) reaction that is a mode of superoxide radical formation and propagation: aO2+bO2·−→aO2·−+bO2. Dioxygen/superoxide CT involves a diatomic oxygen dimer encounter complex (EC). I used time-dependent open shell density functional theory (TD-ODFT)/CAM−B3LYP to show that the parallel quartet EC had an optimized contact distance of 2.62Å, and a fast oscillating free radical CT reaction time of 97.2 femtoseconds.

Elucidating Potential Energy Surfaces for Singlet O2 Reactions with Protonated, Deprotonated, and Di-Deprotonated Cystine Using a Combination of Approximately Spin-Projected Density Functional Theory and Guided-Ion-Beam Mass Spectrometry

The reactivity of cystine toward electronically excited singlet O2 (a1Δg) has been long debated, despite the fact that most organic disulfides are susceptible to oxidation by singlet O2. We report a combined experimental and computational study on reactions of singlet O2 with gas-phase cystine at different ionization and hydration states, aimed to determine reaction outcomes, mechanisms, and potential energy surfaces (PESs). Ion-molecule collisions of protonated and di-deprotonated cystine ions with singlet O2, in both the absence and the presence of a water ligand, were measured over a center-of-mass collision energy (Ecol) range from 0.1 to 1.0 eV, using a guided-ion-beam scattering tandem mass spectrometer. No oxidation was observed for these reactant ions except collision-induced dissociation at high energies. Guided by density functional theory (DFT)-calculated PESs, reaction coordinates were established to unravel the origin of the nonreactivity of cystine ions toward singlet O2. To account for mixed open- and closed-shell characters, singlet O2 and critical structures along reaction coordinates were evaluated using broken-symmetry, open-shell DFT with spin contamination errors removed by an approximate spin-projection method. It was found that collision of protonated cystine with singlet O2 follows a repulsive potential surface and possesses no chemically significant interaction and that collision-induced dissociation of protonated cystine is dominated by loss of water and CO. Collision of di-deprotonated cystine with singlet O2, on the other hand, forms a short-lived electrostatically bonded precursor complex at low Ecol. The latter may evolve to a covalently bonded persulfoxide, but the conversion is blocked by an activation barrier lying 0.39 eV above reactants. At high Ecol, C-S bond cleavage dominates the collision-induced dissociation of di-deprotonated cystine, leading to charge-separated fragmentation. Cross section for the ensuing fragment ion H2NCH(CO2-)CH2SS• was measured as a function of Ecol, and the mechanism of charge-separated fragmentation was discussed. It was also found that the reaction of deprotonated cystine with singlet O2 follows a similar mechanism as that of di-deprotonated cystine, but with an even higher activation barrier (0.72 eV).

Atomic order and cluster energetics of a 17 wt% Si-based glass versus the liquid phase

Aerodynamic levitation of a multicomponent 17 wt% Si glass formed by rapid quenching of the melt phase was studied by high resolution x-ray diffraction (XRD) and reverse Monte Carlo (RMC) modelling. The main local atomic order features comprised interactions between Si, Fe and Mg polyhedra, the stereochemistry of which was on a par with the literature. Both the glass and the liquid state appeared to consist of the same fundamental Si–O, Fe–O and Mg–O clusters, with only the relative number of each varying between the two. Transition from liquid to glass involved a three-fold decrease in uncoordinated O (to within the first minimum of the total g(r)) and a marked increase of Fe–Si–Mg polyhedra bridging O. Octahedral Fe coordination was not suggested by the RMC data. All-electron open-shell density functional theory (DFT) calculations of the most prominent clusters suggested independence between the Fe oxidation state and its polyhedra O-coordination. Of secondary thermodynamic importance were indications of network-forming Fe2+ and Fe3+ distorted trigonal and tetrahedral polyhedra. In all occasions, the Fe ferrous and ferric states involved comparable binding energies within similar clusters which indicate a dynamic equilibrium between the two.

Theoretical study on vibrational circular dichroism spectra of tris(acetylacetonato)metal(III) complexes: Anharmonic effects and low-lying excited states

The open-shell density functional theory calculations with M06 exchange-correlation functional and all-electron Douglas-Kroll second order scalar relativistic correction were performed to interpret the vibrational circular dichroism (VCD) spectra of four kinds of tris(acetylacetonato)metal(III), [M(III)(acac)(3)] (acac = acetylacetonato, M = Ru, Cr, Co, and Rh). It was deduced that the experimental spectra were well reproduced by the calculation with harmonic approximation in case of [Co(III)(acac)(3)] (d(6); S = 0), [Rh(III)(acac)(3)] (d(6); S = 0), and [Ru(III)(acac)(3)] (d(5); S = 1/2). In case of [Cr(III)(acac)(3)] (d(3); S = 3/2), anharmonic effects should be taken into account to predict the accurate vibrational frequencies of closely located modes. Time-dependent density functional theory calculations were performed to estimate the contribution of excited states in the VCD spectra. As a consequence, the presence of the low-lying excited states was predicted for [Ru(III)(acac)(3)] alone, which agreed with the experimental observation.

Synthesis and Electronic Structure of Cationic, Neutral, and Anionic Bis(imino)pyridine Iron Alkyl Complexes: Evaluation of Redox Activity in Single-Component Ethylene Polymerization Catalysts

A family of cationic, neutral, and anionic bis(imino)pyridine iron alkyl complexes has been prepared, and their electronic and molecular structures have been established by a combination of X-ray diffraction, Mössbauer spectroscopy, magnetochemistry, and open-shell density functional theory. For the cationic complexes, [((iPr)PDI)Fe-R][BPh(4)] ((iPr)PDI = 2,6-(2,6-(i)Pr(2)-C(6)H(3)N═CMe)(2)C(5)H(3)N; R = CH(2)SiMe(3), CH(2)CMe(3), or CH(3)), which are known single-component ethylene polymerization catalysts, the data establish high spin ferrous compounds (S(Fe) = 2) with neutral, redox-innocent bis(imino)pyridine chelates. One-electron reduction to the corresponding neutral alkyls, ((iPr)PDI)Fe(CH(2)SiMe(3)) or ((iPr)PDI)Fe(CH(2)CMe(3)), is chelate-based, resulting in a bis(imino)pyridine radical anion (S(PDI) = 1/2) antiferromagnetically coupled to a high spin ferrous ion (S(Fe) = 2). The neutral neopentyl derivative was reduced by an additional electron and furnished the corresponding anion, [Li(Et(2)O)(3)][((iPr)PDI)Fe(CH(2)CMe(3))N(2)], with concomitant coordination of dinitrogen. The experimental and computational data establish that this S = 0 compound is best described as a low spin ferrous compound (S(Fe) = 0) with a closed-shell singlet bis(imino)pyridine dianion (S(PDI) = 0), demonstrating that the reduction is ligand-based. The change in field strength of the bis(imino)pyridine coupled with the placement of the alkyl ligand into the apical position of the molecule induced a spin state change at the iron center from high to low spin. The relevance of the compounds and their electronic structures to olefin polymerization catalysis is also presented.

Modelling the visible absorption spectra of copper(II) acetylacetonate by Density Functional Theory

Spin restricted open-shell Density Functional Theory calculations have been carried out by means of linear response theory to investigate the visible absorption spectrum of copper(II) acetylacetonate complex, Cu(acac) 2. The 3 d → 3 d transition energies and the influence of molecular structure and non-coordinating solvent on the spectra have been investigated. The obtained four 3 d → 3 d transition energies accord well with the experimental data in the crystal phase. The presented results indicate that the experimentally observed four band structure of Cu(acac) 2 is of molecular nature, and not caused by factor-group splitting in the crystal environment as previously suggested.

Time-dependent closed and open-shell density functional theory from the perspective of partitioning techniques and projections

The formalism of time-dependent density functional theory and density functional response functions is reviewed from the perspective of partitioning techniques and projections, the hallmark of the school of quantum chemistry represented by P.O. Löwdin. Extension to open-shell density functional theory is described, reviewing some of its common problems. Sample calculations on hyper-Raleigh scattering in the diamagnetic and paramagnetic complexes of copper bound to (glycyl)glycine dipeptide are presented.

Structural Origin of Copper Ion Containing Artificial DNA: A Density Functional Study

In order to investigate an origin of structural stability of a copper ion containing artificial DNA, we evaluated the stacking energy between [H-Cu(2+)-H] (H: hydroxypyridone) dimer by means of density functional theory (DFT) with an Anderson-Langreth-Lindqvist van der Waals (vdW) functional. The calculated distance between the copper ions is about 3.6 angstroms, which agrees well with the experimental data. Evaluated stacking energy is about 8-10 kcal/mol, which is slightly smaller than that of two base pairs in a natural B-DNA. This tendency does not change in [H-2H(+)-H], which does not contain copper ions. These results indicate that the vdW interaction dominates the inter-base-pair interaction over spin-spin interaction, in contrast to a conjecture by an experimental group. According to the results by the open-shell DFT, antiferromagnetic (singlet) and ferromagnetic (triplet) states are almost degenerated when the two bases are vertically located and both bases have a planar structure as found in the B-DNA.

Spin−Spin Contributions to the Zero-Field Splitting Tensor in Organic Triplets, Carbenes and BiradicalsA Density Functional and Ab Initio Study

An evaluation study for the direct dipolar electron spin-spin (SS) contribution to the zero-field splitting (ZFS) tensor in electron paramagnetic resonance (EPR) spectroscopy is presented. Calculations were performed on a wide variety of organic systems where the SS contribution to the ZFS dominates over the second-order spin-orbit coupling (SOC) contribution. Calculations were performed using (hybrid) density functional theory (DFT), as well as complete active space self-consistent field (CASSCF) wave functions. In the former case, our implementation is an approximation, because we use the two-particle reduced spin-density matrix of the noninteracting reference system. In the latter case, the SS contribution is approximated by a mean-field method which, nevertheless, gives accurate results, compared to the approximation free computation of the SS part in a CASSCF framework. For the case of the triplet dioxygen molecule, it was shown that restricted open-shell density functional theory (RODFT), as well as CASSCF, can provide accurate spin-spin couplings while spin-unrestricted DFT leads to much larger errors. Furthermore, 15 organic radicals, including several 1,3 and 1,5 diradicals, dinitroxide biradicals, and even a chlorophyll a model system, were examined as test cases to demonstrate the accuracy and efficiency of our approach within a DFT framework. Accurate D values with root-mean-square deviations of 0.0035 cm(-1) were obtained. Furthermore, all trends, including those due to substituent effects, were correctly reproduced. In a different set of calculations, the polyacenes benzene, naphthalene, anthracene, and tetracene were studied. Applying DFT, the absolute D values were noticeably underestimated, but it was possible to correctly reproduce the trend to smaller D values with larger size of the systems. Finally, it was demonstrated that our approach is also well-suited for the study of carbenes. The smaller organic radicals of this work were also studied, through the use of CASSCF wave functions. This was a special advantage in the case of the triplet polyacenes, where the CASSCF approach gave better results than the DFT method. In comparing spin-restricted and spin-unrestricted results, it was shown through a natural orbital analysis and comparison to high-level ab initio calculations that even small amounts of spin polarization introduced by the unrestricted calculations lead to large deviations between the unrestricted Kohn-Sham (UKS) and restricted open-shell Kohn-Sham (ROKS) approaches. It is challenging to understand why the ROKS results show much better correlation with the experimental data.