Electron-Nuclear Correlation in Laser-Driven Diatomic Molecules

The interaction of short intense laser pulses with atoms/molecules produces a multitude of highly nonlinear processes requiring a non-perturbative treatment. Detailed study of these highly nonlinear processes by numerically solving the time-dependent Schrodinger equation becomes a daunting task when the number of degrees of freedom is large. Also the coupling between the electronic and nuclear degrees of freedom further aggravates the computational problems. In the present work we show that the time-dependent Hartree (TDH) approximation, which neglects the correlation effects, gives unreliable description of the system dynamics both in the absence and presence of an external field. A theoretical framework is required that treats the electrons and nuclei on equal footing and fully quantum mechanically. To address this issue we discuss two approaches, namely the multicomponent density functional theory (MCDFT) and the multiconfiguration time-dependent Hartree (MCTDH) method, that go beyond the TDH approximation and describe the correlated electron-nuclear dynamics accurately. In the MCDFT framework, where the time-dependent electronic and nuclear densities are the basic variables, we discuss an algorithm to calculate the exact Kohn-Sham (KS) potentials for small model systems. By simulating the photodissociation process in a model hydrogen molecular ion, we show that the exact KS potentials contain all the many-body effects and give an insight into the system dynamics. In the MCTDH approach, the wave function is expanded as a sum of products of single-particle functions (SPFs). The MCTDH method is able to describe the electron-nuclear correlation effects as the SPFs and the expansion coefficients evolve in time and give an accurate description of the system dynamics. We show that the MCTDH method is suitable to study a variety of processes such as the fragmentation of molecules, high-order harmonic generation, the two-center interference effect, and the lochfrass effect. We discuss these phenomena in a model hydrogen molecular ion and a model hydrogen molecule. Inclusion of absorbing boundaries in the mean-field approximation and its consequences are discussed using the model hydrogen molecular ion. To this end, two types of calculations are considered: (i) a variational approach with a complex absorbing potential included in the full many-particle Hamiltonian and (ii) an approach in the spirit of time-dependent density functional theory (TDDFT), including complex absorbing potentials in the single-particle equations. It is elucidated that for small grids the TDDFT approach is superior to the variational approach.

Many-body theory of laser-induced ultrafast demagnetization and angular momentum transfer in ferromagnetic transition metals

An electronic theory is developed, which describes the ultrafast demagnetization in itinerant ferromagnets following the absorption of a femtosecond laser pulse. The present work intends to elucidate the microscopic physics of this ultrafast phenomenon by identifying its fundamental mechanisms. In particular, it aims to reveal the nature of the involved spin excitations and angular-momentum transfer between spin and lattice, which are still subjects of intensive debate. In the first preliminary part of the thesis the initial stage of the laser-induced demagnetization process is considered. In this stage the electronic system is highly excited by spin-conserving elementary excitations involved in the laser-pulse absorption, while the spin or magnon degrees of freedom remain very weakly excited. The role of electron-hole excitations on the stability of the magnetic order of one- and two-dimensional 3d transition metals (TMs) is investigated by using ab initio density-functional theory. The results show that the local magnetic moments are remarkably stable even at very high levels of local energy density and, therefore, indicate that these moments preserve their identity throughout the entire demagnetization process. In the second main part of the thesis a many-body theory is proposed, which takes into account these local magnetic moments and the local character of the involved spin excitations such as spin fluctuations from the very beginning. In this approach the relevant valence 3d and 4p electrons are described in terms of a multiband model Hamiltonian which includes Coulomb interactions, interatomic hybridizations, spin-orbit interactions, as well as the coupling to the time-dependent laser field on the same footing. An exact numerical time evolution is performed for small ferromagnetic TM clusters. The dynamical simulations show that after ultra-short laser pulse absorption the magnetization of these clusters decreases on a time scale of hundred femtoseconds. In particular, the results reproduce the experimentally observed laser-induced demagnetization in ferromagnets and demonstrate that this effect can be explained in terms of the following purely electronic non-adiabatic mechanism: First, on a time scale of 10–100 fs after laser excitation the spin-orbit coupling yields local angular-momentum transfer between the spins and the electron orbits, while subsequently the orbital angular momentum is very rapidly quenched in the lattice on the time scale of one femtosecond due to interatomic electron hoppings. In combination, these two processes result in a demagnetization within hundred or a few hundred femtoseconds after laser-pulse absorption.

A chemical Hamiltonian approach study of the basis set superposition error changes on electron densities and one- and two-center energy components

A comparision of the local effects of the basis set superposition error (BSSE) on the electron densities and energy components of three representative H-bonded complexes was carried out. The electron densities were obtained with Hartee-Fock and density functional theory versions of the chemical Hamiltonian approach (CHA) methodology. It was shown that the effects of the BSSE were common for all complexes studied. The electron density difference maps and the chemical energy component analysis (CECA) analysis confirmed that the local effects of the BSSE were different when diffuse functions were present in the calculations

A new all-round density functional based on spin states and SN2 barriers

We report here a new empirical density functional that is constructed based on the performance of OPBE and PBE for spin states and SN 2 reaction barriers and how these are affected by different regions of the reduced gradient expansion. In a previous study [Swart, Sol̀, and Bickelhaupt, J. Comput. Methods Sci. Eng. 9, 69 (2009)] we already reported how, by switching between OPBE and PBE, one could obtain both the good performance of OPBE for spin states and reaction barriers and that of PBE for weak interactions within one and the same (SSB-sw) functional. Here we fine tuned this functional and include a portion of the KT functional and Grimme's dispersion correction to account for π- π stacking. Our new SSB-D functional is found to be a clear improvement and functions very well for biological applications (hydrogen bonding, π -π stacking, spin-state splittings, accuracy of geometries, reaction barriers)

Infrared intensities of furan, pyrrole and thiophene: beyond the double harmonic approximation

Infrared intensities of the fundamental, overtone and combination transitions in furan, pyrrole and thiophene have been calculated using the variational normal coordinate code MULTIMODE. We use pure vibrational wavefunctions, and quartic force fields and cubic dipole moment vector surfaces, generated by density functional theory. The results are compared graphically with second-order perturbation calculations and with relative intensities from experiment for furan and pyrrole.

A DFT study of free radicals formed from artemisinin and related compounds

Ab initio calculations using density functional theory have shown that the reactions that occur between artemisinin, 1, a cyclic trioxane active against malaria, and some metal ions and complexes lead to a series of radicals which are probably responsible for its therapeutic activity. In particular it has been shown that the interaction of Fe(H) with artemisinin causes the O-O bond to be broken as indeed does Fe(III) and Cu(I), while Zn(II) does not. Calculations were carried out with Fe(II) in several different forms including the bare ion, [Fe(H2O)(5)](2+) and [FeP(Im)] (P, porphyrin; Im, imadazole) and similar results were obtained. The resulting oxygen-based radicals are readily converted to more stable carbon-based radicals and/or. stable products. Similar radicals and products are also formed from two simple model trioxanes 2 and 3 that show little or no therapeutic action against malaria although some subtle differences were obtained. This suggests that the scaffold surrounding the pharmacophore may be involved in molecular recognition events allowing efficient uptake of this trioxane warhead into the parasite. (C) 2004 Elsevier B.V. All rights reserved.

A pulse radiolysis study of free radicals formed by one-electron oxidation of the antimalarial drug pyronaridine

Free radicals from one-electron oxidation of the antimalarial drug pyronaridine have been studied by pulse radiolysis. The results show that pyronaridine is readily oxidised to an intermediate semi-iminoquine radical by inorganic and organic free radicals, including those derived from tryptophan and acetaminophen. The pyronaridine radical is rapidly reduced by both ascorbate and caffeic acid. The results indicate that the one-electron reduction potential of the pyronaridine radical at neutral pH lies between those of acetaminophen (707 mV) and caffeic acid (534 mV). The pyronaridine radical decays to produce the iminoquinone, detected by electrospray mass spectrometry, in a second-order process that density functional theory (DFT) calculations (UB3LYP/6-31+G*) suggest is a disproportionation reaction. Important calculated dimensions of pyronaridine, its phenoxyl and aminyl radical, as well as the iminoquinone, are presented.

UV absorption spectra of methyl-substituted hydroxy-cyclohexadienyl radicals in the gas phase

UV absorption spectra of five methyl-substituted hydroxy-cyclohexadienyl radicals, formed by the addition of the hydroxyl radical (OH) to toluene (methyl benzene), o-, m- and p-xylene (1,2-, 1,3- and 1,4-dimethyl benzene, respectively) and mesitylene (1,3,5-trimethylbenzene), have been determined at 298 K, 1 atm pressure (N-2 + O-2), and the corresponding absolute absorption cross-sections measured, using laser flash photolysis and time-resolved UV absorption detection. As observed for other cyclohexadienyl-type radicals, a strong absorption band is present in the 260-340 nm spectral region, with maximum cross-sections in the range (0.9-2.2) x 10(-17) cm(2) molecule(-1). The shape of the band varies significantly from one radical to the next for the series of aromatic precursors investigated. The nature and yields of hydroxylated ring-retaining oxidation products, identified in previous studies of the OH-initiated oxidation of aromatic hydrocarbons, and the results of theoretical density functional theory (DFT) calculations indicate that one or more possible isomers of the various OH-adducts may contribute to the observed spectra. Isomers where the OH-group is ortho- (or both ortho- and ipso-) to a substituent methyl-group are likely to be the most abundant but other isomers may also be formed to a significant extent. Nonetheless, the present study provides absorption spectra of the adduct radicals formed from the gas phase addition of OH to the aromatic hydrocarbons considered, near room temperature and I atm pressure. (c) 2005 Elsevier B.V. All rights reserved.

Surface geometry of Cu{531}

We present a combined quantitative low-energy electron diffraction (LEED) and density-functional theory (DFT) study of the chiral Cu{531} surface. The surface shows large inward relaxations with respect to the bulk interlayer distance of the first two layers and a large expansion of the distance between the fourth and fifth layers. (The latter is the first layer having the same coordination as the Cu atoms in the bulk.) Additional calculations have been performed to study the likelihood of faceting by comparing surface energies of possible facet terminations. No overall significant reduction in energy with respect to planar {531} could be found for any of the tested combinations of facets, which is in agreement with the experimental findings.

The local adsorption geometry and electronic structure of alanine on Cu{110}

The adsorption of alanine on Cu {110} was studied by a combination of near edge X-ray absorption fine structure (NEXAFS) spectroscopy, X-ray photoelectron spectroscopy (XPS) and density functional theory (DFT). Large chemical shifts in the C 1s, N 1s, and O 1s XP spectra were found between the alanine multilayer and the chemisorbed and pseudo-(3 x 2) alaninate layers. From C, N, and O K-shell NEXAFS spectra the tilt angles of the carboxylate group (approximate to 26 degrees in plane with respect to [1 (1) over bar0] and approximate to 45 degrees out of plane) and the C-N bond angle with respect to [1 (1) over bar0] could be determined for the pseudo-(3 x 2) overlayer. Using this information three adsorption geometries could be eliminated from five p(3 x 2) structures which lead to almost identical heats of adsorption in the DFT calculations between 1.40 and 1.47 eV/molecule. Due to the small energy difference between the remaining two structures it is not unlikely that these coexist on the surface at room temperature. (c) 2006 Elsevier B.V. All rights reserved.

Ag(I) and Cu(I) complexes of tetramethyldiphosphinedisulfide: synthesis and structure

Reaction of [M(NCCH3)(4)][PF6] (M = Ag, Cu) with the S2P2Me4 ligand in dichloromethane solution led to substitution of all the nitrile ligands by two molecules of the sulfur ligand, affording the new species [Ag(S2P2Me4)(2)][PF6] (1) and [Cu(S2P2Me4)(2)][PF6] (2). The structures of these complexes were determined by single crystal X-ray diffraction. showing the expected tetrahedral coordination around each metal. Density functional theory (DFT) calculations confirmed the different geometries and energies of the free and coordinated ligand, and provided a very good reproduction of the experimental structures, both for Ag and Cu. The lengths of the S=P bonds are barely affected by coordination, indicating that the pi bond is not important in binding to the metal. (C) 2002 Elsevier Science B.V. All rights reserved.

The Structure of the chiral Pt{531} surface: a combined LEED and DFT study

The structure of the chiral kinked Pt{531} surface has been determined by low-energy electron diffraction intensity-versus-energy (LEED-IV) analysis and density functional theory (DFT). Large contractions and expansions of the vertical interlayer distances with respect to the bulk-terminated surface geometry were found for the first six layers (LEED: d(12) = 0.44 angstrom, d(23) = 0.69 angstrom, d(34) = 0.49 angstrom, d(45) = 0.95 angstrom, d(56) = 0.56 angstrom; DFT: d(12) = 0.51 angstrom, d(23) = 0.55 angstrom, d(34) = 0.74 angstrom, d(45) = 0.78 angstrom, d(56) = 0.63 angstrom; d(bulk) = 0.66 angstrom). Energy-dependent cancellations of LEED spots over unusually large energy ranges, up to 100 eV, can be explained by surface roughness and reproduced by applying a model involving 0.25 ML of vacancies and adatoms in the scattering calculations. The agreement between the results from LEED and DFT is not as good as in other cases, which could be due to this roughness of the real surface.

Electronic transitions and bonding properties in a series of five-coordinate “16-electron” complexes [Mn(CO)3(L2)]− (L2=chelating redox-active π-donor ligand)

In this article we present for the first time accurate density functional theory (DFT) and time-dependent (TD) DFT data for a series of electronically unsaturated five-coordinate complexes [Mn(CO)(3)(L-2)](-), where L-2 stands for a chelating strong pi-donor ligand represented by catecholate, dithiolate, amidothiolate, reduced alpha-diimine (1,4-dialkyl-1,4-diazabutadiene (R-DAB), 2,2'-bipyridine) and reduced 2,2'-biphosphinine types. The single-crystal X-ray structure of the unusual compound [Na(BPY)][Mn(CO)(3)(BPY)]center dot Et2O and the electronic absorption spectrum of the anion [Mn(CO)(3)(BPY)](-) are new in the literature. The nature of the bidentate ligand determines the bonding in the complexes, which varies between two limiting forms: from completely pi-delocalized diamagnetic {(CO)(3)Mn-L-2}(-) for L-2 = alpha-diimine or biphosphinine, to largely valence-trapped {(CO)(3)Mn-1-L-2(2-)}(-) for L-2(2-) = catecholate, where the formal oxidation states of Mn and L-2 can be assigned. The variable degree of the pi-delocalization in the Mn(L-2) chelate ring is indicated by experimental resonance Raman spectra of [Mn(CO)(3)(L-2)](-) (L-2=3,5-di-tBu-catecholate and iPr-DAB), where accurate assignments of the diagnostically important Raman bands have been aided by vibrational analysis. The L-2 = catecholate type of complexes is known to react with Lewis bases (CO substitution, formation of six-coordinate adducts) while the strongly pi-delocalized complexes are inert. The five-coordinate complexes adopt usually a distorted square pyramidal geometry in the solid state, even though transitions to a trigonal bipyramid are also not rare. The experimental structural data and the corresponding DFT-computed values of bond lengths and angles are in a very good agreement. TD-DFT calculations of electronic absorption spectra of the studied Mn complexes and the strongly pi-delocalized reference compound [Fe(CO)(3)(Me-DAB)] have reproduced qualitatively well the experimental spectra. Analyses of the computed electronic transitions in the visible spectroscopic region show that the lowest-energy absorption band always contains a dominant (in some cases almost exclusive) contribution from a pi(HOMO) -> pi*(LUMO) transition within the MnL2 metallacycle. The character of this optical excitation depends strongly on the composition of the frontier orbitals, varying from a partial L-2 -> Mn charge transfer (LMCT) through a fully delocalized pi(MnL2) -> pi*(MnL2) situation to a mixed (CO)Mn -> L-2 charge transfer (LLCT/MLCT). The latter character is most apparent in the case of the reference complex [Fe(CO)(3)(Me-DAB)]. The higher-lying, usually strongly mixed electronic transitions in the visible absorption region originate in the three lower-lying occupied orbitals, HOMO - 1 to HOMO - 3, with significant metal-d contributions. Assignment of these optical excitations to electronic transitions of a specific type is difficult. A partial LLCT/MLCT character is encountered most frequently. The electronic absorption spectra become more complex when the chelating ligand L-2, such as 2,2'-bipyridine, features two or more closely spaced low-lying empty pi* orbitals.

Heterosite effects in novel heteronuclear clusters [OS2Ru(CO)(11)(PPh3)] and Os2Ru(CO)(10)(2-acetylpyridine-N-isopropylimine)]

A new synthetic route towards the mixed-metal cluster [OS2Ru(CO)(12)] is described together with the syntheses of its PPh3 and iPr-AcPy (iPr-AcPy = 2-acetylpyridine-N-isopropylimine) derivatives. The molecular structures of the novel clusters [Os2Ru(CO)(11)(PPh3)] and [Os2Ru(CO)(10)(iPr-AcPy)] were determined on the basis of crystalline solid solutions of the Os2Ru and corresponding Os-3 species. The structures reveal that coordination of the Lewis bases occurs exclusively at the ruthenium site of [Os2Ru(CO)(12)], which is in agreement with density functional theory (DFT) calculations on several structural isomers of these compounds. According to the time-dependent DFT results, the lowest optically accessible excited state of [Os2Ru(CO)(10)(iPr-AcPy)] has a prevailing sigma(Ru-Os-2)pi*(iPr-AcPy) character, with a partial sigma sigma*(Ru-Os-2) contribution. In weakly coordinating 2-chlorobutane, the excited state has a lifetime tau = 10.4 +/- 1.2 ps and produces biradicals considerably faster than observed for [Os-3(CO)10(iPr-AcPy) (tau = 25.3 +/- 0.7ps)]. In coordinating acetonitrile, the excited state of [Os2Ru(CO)(10)(iPr-AcPy)] decays mono-exponentially with a lifetime tau = 2.1 +/- 0.2 ps. In contrast to [Os-3(CO)(10)(iPr-AcPy)] that forms biradicals as the main primary photoproduct even in strongly coordinating solvents, zwitterion formation from the solvated lowest excited state is observed for the heterometallic cluster. This is concluded from time-resolved absorption studies in the microsecond time domain. Due to the lower tendency of the coordinatively unsaturated Ru+(CO)(2)(iPr-AcPy-/0) moiety to bind a Lewis base, the heteronuclear biradical and zwitterionic photoproducts live significantly shorter than their triosmium counterparts. The influence of the weaker Os-2-Ru(iPr-AcPy) bond on the redox reactivity is clearly reflected in very reactive radical anions formed upon electrochemical reduction of [Os2Ru(CO)(10)(iPr-AcPy)]. The dimer [-OS(CO)(4)-Os(CO)(4)-Ru(CO)(2)(iPr-AcPy)](2)(2-) is the only IR-detectable intermediate reduction product. The dinuclear complex [Os-2(CO)(8)](2-) and insoluble [Ru(CO)(2)(iPr-AcPy)](n), are the ultimate reduction products, proving fragmentation of the OS2Ru core.

Magnetic coupling in trinuclear partial cubane copper(II) complexes with a hydroxo bridging core and peripheral phenoxo bridges from NNO donor Schiff base ligands

Three new trinuclear copper(II) complexes, [(CuL1)(3)(mu(3)-OH)](ClO4)(2)center dot 3.75H(2)O (1), [(CuL2)(3)(mu(3)-OH)](ClO4)(2) (2) and [(CuL3)(3)(mu(3)-OH)](BF4)(2)center dot 0.5CH(3)CN (3) have been synthesized from three tridentate Schiff bases HL1, HL2, and HL3 (HL1 = 2-[(2-amino-ethylimino)-methyl]-phenol, HL2 = 2-[(2-methylamino-ethylimino)-methyl]-phenol and HL3 = 2-[1-(2-dimethylamino-ethylimino)-ethyl]-phenol). The complexes are characterized by single-crystal X-ray diffraction analyses, IR, UV-vis and EPR spectroscopy, and variable-temperature magnetic measurements. All the compounds contain a partial cubane [Cu3O4] core consisting of the trinuclear unit [(CuL)(3)(mu(3)-OH)](2+) together with perchlorate or fluoroborate anions. In each of the complexes, the three copper atoms are five-coordinated with a distorted square-pyramidal geometry except in complex 1, in which one of the Cu-II ions of the trinuclear unit is six-coordinate being in addition weakly coordinated to one of the perchlorate anions. Variable-temperature magnetic measurements and EPR spectra indicate an antiferromagnetic exchange coupling between the CuII ions of complexes 1 and 2, while this turned out to be ferromagnetic for complex 3. Experimental values have been fitted according to an isotropic exchange Hamiltonian. Calculations based on Density Functional Theory have also been performed in order to estimate the exchange coupling constants in these three complexes. Both sets of values indicate similar trends and specially calculated J values establish a magneto-structural correlation between them and the Cu-O-Cu bond angle, in that the coupling is more ferromagnetic for smaller bond angle values.