Functionals for TDDFT description of electron transport in nanostructures

La Teoria di Densità Funzionale (DFT) e la sua versione dipendente dal tempo (TDDFT) sono strumenti largamente usati per simulare e calcolare le proprietà statiche e dinamiche di sistemi con elettroni interagenti. La precisione del metodo si basa su una serie di approssimazioni degli effetti di exchange correlation fra gli elettroni, descritti da un funzionale della sola densità di carica. Nella presente tesi viene testata l'affidabilità del funzionale Mixed Localization Potential (MLP), una media pesata fra Single Orbital Approximation (SOA) e un potenziale di riferimento, ad esempio Local Density Approximation (LDA). I risultati mostrano capacità simulative superiori a LDA per i sistemi statici (curando anche un limite di LDA noto in letteratura come fractional dissociation) e dei progressi per sistemi dinamici quando si sviluppano correnti di carica. Il livello di localizzazione del sistema, inteso come la capacità di un elettrone di tenere lontani da sé altri elettroni, è descritto dalla funzione Electron Localization Function (ELF). Viene studiato il suo ruolo come guida nella costruzione e ottimizzazione del funzionale MLP.

Characterization of High-Performance Organic Dyes for Dye-Sensitized Solar Cell: A DFT/TDDFT Study

Dye-sensitized solar cell (DSSC) is currently a promising technology that makes solar energy efficient and cost-effective to harness. In DSSC, metal free dyes, such indoline-containing D149 and D205, are proved to be potential alternatives for traditional metal organic dyes. In this work, a DFT/TDDFT characterization for D149 and D205 were carried out using different functionals, including B3LYP, MPW1K, CAM-B3LYP and PBE0. Three different conformers for D149 and four different conformers for D205 were identified and calculated in vacuum. The performance of different functionals on calculating the maximum absorbance of the dyes in vacuum and five common solvents (acetonitrile, chloroform, ethanol, methanol, and THF) were examined and compared to determine the suitable computational setting for predicting properties of these two dyes. Furthermore, deprotonated D149 and D205 in solvents were also considered, and the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) were calculated, which elucidates the substitution effect on the rhodanine ring of D149 and D205 dyes on their efficiency. Finally, D149 and D205 molecules were confirmed to be firmly anchored on ZnO surface by periodic DFT calculations. These results would shed light on the design of new highly efficiency metal-free dyes.

Spin gaps and spin-flip energies in density-functional theory

Energy gaps are crucial aspects of the electronic structure of finite and extended systems. Whereas much is known about how to define and calculate charge gaps in density-functional theory (DFT), and about the relation between these gaps and derivative discontinuities of the exchange-correlation functional, much less is known about spin gaps. In this paper we give density-functional definitions of spin-conserving gaps, spin-flip gaps and the spin stiffness in terms of many-body energies and in terms of single-particle (Kohn-Sham) energies. Our definitions are as analogous as possible to those commonly made in the charge case, but important differences between spin and charge gaps emerge already on the single-particle level because unlike the fundamental charge gap spin gaps involve excited-state energies. Kohn-Sham and many-body spin gaps are predicted to differ, and the difference is related to derivative discontinuities that are similar to, but distinct from, those usually considered in the case of charge gaps. Both ensemble DFT and time-dependent DFT (TDDFT) can be used to calculate these spin discontinuities from a suitable functional. We illustrate our findings by evaluating our definitions for the Lithium atom, for which we calculate spin gaps and spin discontinuities by making use of near-exact Kohn-Sham eigenvalues and, independently, from the single-pole approximation to TDDFT. The many-body corrections to the Kohn-Sham spin gaps are found to be negative, i.e., single-particle calculations tend to overestimate spin gaps while they underestimate charge gaps.

Physical signatures of discontinuities of the time-dependent exchange-correlation potential

The exact exchange-correlation (XC) potential in time-dependent density-functional theory (TDDFT) is known to develop steps and discontinuities upon change of the particle number in spatially confined regions or isolated subsystems. We demonstrate that the self-interaction corrected adiabatic local-density approximation for the XC potential has this property, using the example of electron loss of a model quantum well system. We then study the influence of the XC potential discontinuity in a real-time simulation of a dissociation process of an asymmetric double quantum well system, and show that it dramatically affects the population of the resulting isolated single quantum wells. This indicates the importance of a proper account of the discontinuities in TDDFT descriptions of ionization, dissociation or charge transfer processes.

Assignment of the absolute configuration at the sulfur atom of thioridazine metabolites by the analysis of their chiroptical properties: The case of thioridazine 2-sulfoxide

Thioridazine (THD) is a commonly prescribed phenotiazine neuroleptic drug, which is extensively biotransformed in the organism producing as main metabolites sulfoxides and a sulfone by sulfur oxidation Significant differences have been observed in the activity of the THD enantiomers as well as for its main metabolites, and enantioselectivity phenomena have been proved in the metabolic pathway. Here the assignment of the absolute configuration at the sulfur atom of enantiomeric THD-2-sulfoxide (THD-2-SO) has been carried out by circular dichroism (CD) spectroscopy The stereoisomers were separated by HPLC on Chiralpak AS column, recording the CD spectra for the two collected enantiomeric fractions The theoretical electronic CD spectrum has been obtained by the TDDFT/B3LYP/6-31G*. as Boltzmann averaging of the contributions calculated for the most stable conformations of the drug The comparison of the simulated and experimental spectra allowed the absolute configuration at the sulfur atom of the four THD-2-SO stereoisomers to be assigned The developed method should be useful for a reliable correlation between stereochemistry and activity and/or toxicity

A fluorescente bicyclic Calix[4] Arene-Oxacyclophane with planar chirality: Resolution, chiroptical, properties, and absolute configuration.

A new inherently chiral calix[4]arene ICC 1 has been disclosed. The dissymmetry of 1 is generated from a chirality plane in the quinol moiety of a 1,3-bridged bicyclic calix[4]arene. ICC 1 has been resolved by enantioselective HPLC, and the chiroptical properties of both isolated antipodes (pS)-1 and (pR)-1 confirm their enantiomeric nature. The absolute configuration of the (pS)-1/(pR)-1 enantiomeric pair was established through time-dependent density functional theory (TDDFT) calculations of electronic circular dichroism (CD) spectra. (C) 2014 Elsevier Ltd. All rights reserved.

Inherently chiral calix[4]arenes with planar chirality: two new entries to the family

The synthesis of two new inherently chiral calix[4]arenes (ICCs, 1 and 2), endowed with electron-rich concave surfaces, has been achieved through the desymmetrization of a lower rim distal-bridged oxacyclophane (OCP) macrocycle. The new highly emissive ICCs were resolved by chiral HPLC, and the enantiomeric nature of the isolated antipodes proved by electronic circular dichroism (CD). Using time-dependent density functional calculations of CD spectra, their absolute configurations were established. NMR studies with (S)-Pirkle's alcohol unequivocally showed that the host-guest interactions occur in the chiral pocket comprehending the calix-OCP exo cavities and the carbazole moieties.

Mono(Eta(5)-cyclopentadienyl)metal(II) complexes with thienyl acetylide chromophores: synthesis, electrochemical studies, and first hyperpolarizabilities

A series of mono(eta(5)-cyclopentadienyl)metal-(II) complexes with nitro-substituted thienyl acetylide ligands of general formula [M(eta(5)-C5H5)(L)(C C{C4H2S}(n)NO2)] (M = Fe, L = kappa(2)-DPPE, n = 1,2; M = Ru, L = kappa(2)-DPPE, 2 PPh3, n = 1, 2; M = Ni, L = PPh3, n = 1, 2) has been synthesized and fully characterized by NMR, FT-IR, and UV-Vis spectroscopy. The electrochemical behavior of the complexes was explored by cyclic voltammetry. Quadratic hyperpolarizabilities (beta) of the complexes have been determined by hyper-Rayleigh scattering (HRS) measurements at 1500 nm. The effect of donor abilities of different organometallic fragments on the quadratic hyperpolarizabilities was studied and correlated with spectroscopic and electrochemical data. Density functional theory (DFT) and time-dependent DFT (TDDFT) calculations were employed to get a better understanding of the second-order nonlinear optical properties in these complexes. In this series, the complexity of the push pull systems is revealed; even so, several trends in the second-order hyperpolarizability can still be recognized. In particular, the overall data seem to indicate that the existence of other electronic transitions in addition to the main MLCT clearly controls the effectiveness of the organometallic donor ability on the second-order NLO properties of these push pull systems.

Termocromismo em soluções de alcóxidos de Vanádio(IV): uma abordagem pela modelagem molecular

The thermochromic behavior exhibited by vanadium(IV) alkoxides, [V2(μ-OPr i)2(OPr i) 6] and [V2(μ-ONep)2(ONep)6 ], OPr i = isopropoxide and ONep = neopentoxide, was studied by molecular modeling using DFT, TDDFT and INDO/S methods. The vibrational and electronic spectra calculated for [V2(μ-OPr i)2(OPr i) 6] were very similar to the experimental data registered for crystalline samples of the complex and for its solutions at low temperature (< 210 K), while spectra recorded at high temperature (> 315 K) were compatible with those calculated for the monomeric form, [V(OPr i)4]. These results consistently point to a monomer/dimer equilibrium as an explanation for the solution thermochromism of {V(OPr i)4}n. In spite of the structural similarity between [V2(μ-ONep)2(ONep)6 ] and [V2(μ-OPr i)2(OPr i) 6] in the solid state, the thermochromic behavior of the former could not be explained by the same model, and the possibility of tetranuclear aggregation at low temperatures was also investigated.

Design de nouvelles fonctionnelles en théorie de la fonctionnelle de la densité et conception de polymères pour application à la photovoltaïque organique

La présente thèse porte sur les calculs utilisant la théorie de la fonctionnelle de la densité (DFT) pour simuler des systèmes dans lesquels les effets à longue portée sont importants. Une emphase particulière est mise sur les calculs des énergies d’excitations, tout particulièrement dans le cadre des applications photovoltaïques. Cette thèse aborde ces calculs sous deux angles. Tout d’abord, des outils DFT déjà bien établis seront utilisés pour simuler des systèmes d’intérêt expérimental. Par la suite, la théorie sous-jacente à la DFT sera explorée, ses limites seront identifiées et de nouveaux développements théoriques remédiant à ceux-ci seront proposés. Ainsi, dans la première partie de cette thèse, des calculs numériques utilisant la DFT et la théorie de la fonctionnelle de la densité dépendante du temps (TDDFT) telles qu’implémentées dans le logiciel Gaussian [1] sont faits avec des fonctionnelles courantes sur des molécules et des polymères d’intérêt expérimental. En particulier, le projet présenté dans le chapitre 2 explore l’utilisation de chaînes latérales pour optimiser les propriétés électroniques de polymères déjà couramment utilisés en photovoltaïque organique. Les résultats obtenus montrent qu’un choix judicieux de chaînes latérales permet de contrôler les propriétés électroniques de ces polymères et d’augmenter l’efficacité des cellules photovoltaïques les utilisant. Par la suite, le projet présenté dans le chapitre 3 utilise la TDDFT pour explorer les propriétés optiques de deux polymères, le poly-3-hexyl-thiophène (P3HT) et le poly-3-hexyl- sélénophène (P3HS), ainsi que leur mélange, dans le but d’appuyer les observations expérimentales indiquant la formation d’exciplexe dans ces derniers. Les calculs numériques effectués dans la première partie de cette thèse permettent de tirer plusieurs conclusions intéressantes, mais mettent également en évidence certaines limites de la DFT et de la TDDFT pour le traitement des états excités, dues au traitement approximatif de l’interaction coulombienne à longue portée. Ainsi, la deuxième partie de cette thèse revient aux fondements théoriques de la DFT. Plus précisément, dans le chapitre 4, une série de fonctionnelles modélisant plus précisément l’interaction coulombienne à longue portée grâce à une approche non-locale est élaborée. Ces fonctionnelles sont basées sur la WDA (weighted density approximation), qui est modifiée afin d’imposer plusieurs conditions exactes qui devraient être satisfaites par le trou d’échange. Ces fonctionnelles sont ensuite implémentées dans le logiciel Gaussian [1] et leurs performances sont évaluées grâce à des tests effectués sur une série de molécules et d’atomes. Les résultats obtenus indiquent que plusieurs de ces fonctionnelles donnent de meilleurs résultats que la WDA. De plus, ils permettrent de discuter de l’importance relative de satisfaire chacune des conditions exactes.

Structure Absorption Spectra Correlation in a Series of 2,6-Dimethyl-4-arylpyrylium Salts

A combined experimental and theoretical study of the absorption spectra of a group of closely related pyrylium perchlorates 1-11 are presented. Minor changes in the position of the substituents lead to drastic changes in the absorption spectra in this series of compounds. We have attempted to explain the observed changes using the x,y-band notation developed by Balaban and co-workers. Absorption spectra of all compounds are compared with results from time-dependent density functional theory (TDDFT) and Zerner’s intermediate neglect of differential overlap (ZINDO/S) level calculations. Results of the calculations are in good agreement with experimental observations and an interesting correlation between Balaban’s notations and the MO transitions are obtained for simple derivatives. It is suggested that for more complex systems such as R- and â-naphthyl substituted systems, the empirical method is not appropriate.

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.

Nickel(II) and copper(II) complexes of Schiff base ligands containing N4O2 and N4S2 donors with pyrrole terminal binding groups: synthesis, characterization, X-ray structures, DFT and electrochemical studies

A series of hexadentate ligands, H2Lm (m = 1−4), [1H-pyrrol-2-ylmethylene]{2-[2-(2-{[1H-pyrrol-2-ylmethylene]amino}phenoxy)ethoxy]phenyl}amine (H2L1), [1H-pyrrol-2-ylmethylene]{2-[4-(2-{[1H-pyrrol-2-ylmethylene]amino}phenoxy)butoxy]phenyl}amine (H2L2), [1H-pyrrol-2-ylmethylene][2-({2-[(2-{[1H-pyrrol-2-ylmethylene]amino}phenyl)thio]ethyl}thio)phenyl]amine (H2L3) and [1H-pyrrol-2-ylmethylene][2-({4-[(2-{[1H-pyrrol-2-lmethylene]amino}phenyl)thio]butyl}thio) phenyl]amine (H2L4) were prepared by condensation reaction of pyrrol-2-carboxaldehyde with {2-[2-(2-aminophenoxy)ethoxy]phenyl}amine, {2-[4-(2-aminophenoxy)butoxy]phenyl}amine, [2-({2-[(2-aminophenyl)thio]ethyl}thio)phenyl]amine and [2-({4-[(2-aminophenyl)thio]butyl}thio)phenyl]amine respectively. Reaction of these ligands with nickel(II) and copper(II) acetate gave complexes of the form MLm (m = 1−4), and the synthesized ligands and their complexes have been characterized by a variety of physico-chemical techniques. The solid and solution states investigations show that the complexes are neutral. The molecular structures of NiL3 and CuL2, which have been determined by single crystal X-ray diffraction, indicate that the NiL3 complex has a distorted octahedral coordination environment around the metal while the CuL2 complex has a seesaw coordination geometry. DFT calculations were used to analyse the electronic structure and simulation of the electronic absorption spectrum of the CuL2 complex using TDDFT gives results that are consistent with the measured spectroscopic behavior of the complex. Cyclic voltammetry indicates that all copper complexes are electrochemically inactive but the nickel complexes with softer thioethers are more easily oxidized than their oxygen analogs.

π-dimerization of pleiadiene radical cations at low temperatures revealed by UV–vis spectroelectrochemistry and quantum theory

One-electron oxidation of the non-alternant polycyclic aromatic hydrocarbon pleiadiene and related cyclohepta[ c,d]pyrene and cyclohepta[c,d]fluoranthene in THF produces corresponding radical cations detectable in the temperature range of 293–263 K only on the subsecond time scale of cyclic voltammetry. Although the EPR-active red-coloured pleiadiene radical cation is stable according to the literature in concentrated sulfuric acid, spectroelectrochemical measurements reported in this study provide convincing evidence for its facile conversion into the green-coloured, formally closed shell and, hence, EPRsilent π-bound dimer dication stable in THF at 253 K. The unexpected formation of the thermally unstable dimeric product featuring a characteristic intense low-energy absorption band at 673 nm (1.84 eV; logεmax=4.0) is substantiated by ab initio calculations on the parent pleiadiene molecule and the PF6 − salts of the corresponding radical cation and dimer dication. The latter is stabilized with respect to the radical cation by 14.40 kcal mol−1 (DFT B3LYP) [37.64 kcal mol−1 (CASPT2/DFT B3LYP)]. An excellent match has been obtained between the experimental and TDDFT- calculated UV–vis spectra of the PF6 − salt of the pleiadiene dimer dication, considering solvent (THF) effects.

Bis(azidophenyl)phosphole building block for extended π-conjugated systems

Novel bis(azidophenyl)phosphole sulfide building block 8 has been developed to give access to a plethora of phosphole-containing π-conjugated systems in a simple synthetic step. This was explored for the reaction of the two azido moieties with phenyl-, pyridyl- and thienylacetylenes, to give bis(aryltriazolyl)-extended π-systems, having either the phosphole sulfide (9) or the phosphole (10) group as central ring. These conjugated frameworks exhibit intriguing photophysical and electrochemical properties that vary with the nature of the aromatic end-group. The λ3-phospholes 10 display blue fluorescence (λem = 460–469 nm) with high quan-tum yield (ΦF = 0.134–0.309). The radical anion of pyridylsubstituted phosphole sulfide 9b was observed with UV/Vis spectroscopy. TDDFT calculations on the extended π-systems showed some variation in the shape of the HOMOs, which was found to have an effect on the extent of charge transfer, depending on the aromatic end-group. Some fine-tuning of the emission maxima was observed, albeit subtle, showing a decrease in conjugation in the order thienyl � phenyl � pyridyl. These results show that variations in the distal ends of such π-systems have a subtle but significant effect on photophysical properties.