Ground and Electronic Excited States from Pairing Matrix Fluctuation and Particle-Particle Random Phase Approximation

The accurate description of ground and electronic excited states is an important and challenging topic in quantum chemistry. The pairing matrix fluctuation, as a counterpart of the density fluctuation, is applied to this topic. From the pairing matrix fluctuation, the exact electron correlation energy as well as two electron addition/removal energies can be extracted. Therefore, both ground state and excited states energies can be obtained and they are in principle exact with a complete knowledge of the pairing matrix fluctuation. In practice, considering the exact pairing matrix fluctuation is unknown, we adopt its simple approximation --- the particle-particle random phase approximation (pp-RPA) --- for ground and excited states calculations. The algorithms for accelerating the pp-RPA calculation, including spin separation, spin adaptation, as well as an iterative Davidson method, are developed. For ground states correlation descriptions, the results obtained from pp-RPA are usually comparable to and can be more accurate than those from traditional particle-hole random phase approximation (ph-RPA). For excited states, the pp-RPA is able to describe double, Rydberg, and charge transfer excitations, which are challenging for conventional time-dependent density functional theory (TDDFT). Although the pp-RPA intrinsically cannot describe those excitations excited from the orbitals below the highest occupied molecular orbital (HOMO), its performances on those single excitations that can be captured are comparable to TDDFT. The pp-RPA for excitation calculation is further applied to challenging diradical problems and is used to unveil the nature of the ground and electronic excited states of higher acenes. The pp-RPA and the corresponding Tamm-Dancoff approximation (pp-TDA) are also applied to conical intersections, an important concept in nonadiabatic dynamics. Their good description of the double-cone feature of conical intersections is in sharp contrast to the failure of TDDFT. All in all, the pairing matrix fluctuation opens up new channel of thinking for quantum chemistry, and the pp-RPA is a promising method in describing ground and electronic excited states.

Magnetic properties of paramagnetic systems : density functional studies

Thèse numérisée par la Direction des bibliothèques de l'Université de Montréal.

Kinetics and coverage dependent reaction mechanisms of the copper atomic layer deposition from copper dimethylamino-2-propoxide and diethylzinc

Atomic layer deposition (ALD) has been recognized as a promising method to deposit conformal and uniform thin film of copper for future electronic devices. However, many aspects of the reaction mechanism and the surface chemistry of copper ALD remain unclear. In this paper, we employ plane wave density functional theory (DFT) to study the transmetalation ALD reaction of copper dimethylamino-2-propoxide [Cu(dmap)2] and diethylzinc [Et2Zn] that was realized experimentally by Lee et al. [ Angew. Chem., Int. Ed. 2009, 48, 4536−4539]. We find that the Cu(dmap)2 molecule adsorbs and dissociates through the scission of one or two Cu–O bonds into surface-bound dmap and Cu(dmap) fragments during the copper pulse. As Et2Zn adsorbs on the surface covered with Cu(dmap) and dmap fragments, butane formation and desorption was found to be facilitated by the surrounding ligands, which leads to one reaction mechanism, while the migration of ethyl groups to the surface leads to another reaction mechanism. During both reaction mechanisms, ligand diffusion and reordering are generally endothermic processes, which may result in residual ligands blocking the surface sites at the end of the Et2Zn pulse, and in residual Zn being reduced and incorporated as an impurity. We also find that the nearby ligands play a cooperative role in lowering the activation energy for formation and desorption of byproducts, which explains the advantage of using organometallic precursors and reducing agents in Cu ALD. The ALD growth rate estimated for the mechanism is consistent with the experimental value of 0.2 Å/cycle. The proposed reaction mechanisms provide insight into ALD processes for copper and other transition metals.

Using Surface Spectroscopy to Help Interpret STM Images of Chemisorbed Molecules: Application to the Study of the Hydrogenation of Ketones and a-Ketoesters on Pt

Depuis que la haute énantiopureté est nécessaire dans l’industrie pharmaceutique, les études visant à découvrir les mécanismes pour l’hydrogénation énantiosélective de cétones ou céto-esters sur les surfaces, et à rechercher de nouveaux et plus performants catalyseurs asymétriques, sont d’une grande importance. La microscopie à effet tunnel (STM), la spectroscopie infrarouge de réflexion-absorption, la spectroscopie de désorption à température programmée et la spectrométrie de photoélectrons induits par rayons X sont des méthodes performantes facilitant la compréhension des mécanismes de réaction. En plus de nous permettre de comprendre les mécanismes réactionnels, les études peuvent fournir des informations sur la dynamique des réactions en catalyse hétérogène ainsi que sur le développement de la théorie de la fonctionnelle de la densité (DFT) afin de calculer des interactions faibles dans les processus de surface. D’autres parts, les calculs DFT fournissent une aide essentielle à l’interprétation des données de STM et spectroscopie de surface. Dans cette thèse, certains cétones et céto-esters sur la surface de platine sont étudiées par les techniques sophistiquées mentionnées ci-dessus. Mes études démontrent que la combinaison de l’utilisation de la spectroscopie de routine, des nanotechnologies et de nombreux calculs élaborés, est une méthode efficace pour étudier les réactions à la surface car ces techniques explorent les différents aspects de la surface ainsi que s’entraident mutuellement lors de certaines interprétations.

Random path sampling applied to vector models of magnetism

The recently reported Monte Carlo Random Path Sampling method (RPS) is here improved and its application is expanded to the study of the 2D and 3D Ising and discrete Heisenberg models. The methodology was implemented to allow use in both CPU-based high-performance computing infrastructures (C/MPI) and GPU-based (CUDA) parallel computation, with significant computational performance gains. Convergence is discussed, both in terms of free energy and magnetization dependence on field/temperature. From the calculated magnetization-energy joint density of states, fast calculations of field and temperature dependent thermodynamic properties are performed, including the effects of anisotropy on coercivity, and the magnetocaloric effect. The emergence of first-order magneto-volume transitions in the compressible Ising model is interpreted using the Landau theory of phase transitions. Using metallic Gadolinium as a real-world example, the possibility of using RPS as a tool for computational magnetic materials design is discussed. Experimental magnetic and structural properties of a Gadolinium single crystal are compared to RPS-based calculations using microscopic parameters obtained from Density Functional Theory.

Benchmarking Density Functional Based Tight-Binding for Silver and Gold Materials: From Small Clusters to Bulk

International audience

Theoretical Rationalization of the Singlet–Triplet Gap in OLEDs Materials: Impact of Charge-Transfer Character

New materials for OLED applications with low singlet–triplet energy splitting have been recently synthesized in order to allow for the conversion of triplet into singlet excitons (emitting light) via a Thermally Activated Delayed Fluorescence (TADF) process, which involves excited-states with a non-negligible amount of Charge-Transfer (CT). The accurate modeling of these states with Time-Dependent Density Functional Theory (TD-DFT), the most used method so far because of the favorable trade-off between accuracy and computational cost, is however particularly challenging. We carefully address this issue here by considering materials with small (high) singlet–triplet gap acting as emitter (host) in OLEDs and by comparing the accuracy of TD-DFT and the corresponding Tamm-Dancoff Approximation (TDA), which is found to greatly reduce error bars with respect to experiments thanks to better estimates for the lowest singlet–triplet transition. Finally, we quantitatively correlate the singlet–triplet splitting values with the extent of CT, using for it a simple metric extracted from calculations with double-hybrid functionals, that might be applied in further molecular engineering studies.

Designing intrinsically photostable low band gap polymers:a smart tool combining EPR spectroscopy and DFT calculations

A rapid and efficient method to identify the weak points of the complex chemical structure of low band gap (LBG) polymers, designed for efficient solar cells, when submitted to light exposure is reported. This tool combines Electron Paramagnetic Resonance (EPR) using the 'spin trapping method' coupled with density functional theory modelling (DFT). First, the nature of the short life-time radicals formed during the early-stages of photo-degradation processes are determined by a spin-trapping technique. Two kinds of short life-time radical (R and R′O) are formed after 'short-duration' illumination in an inert atmosphere and in ambient air, respectively. Second, simulation allows the identification of the chemical structures of these radicals revealing the most probable photochemical process, namely homolytical scission between the Si atom of the conjugated skeleton and its pendent side-chains. Finally, DFT calculations confirm the homolytical cleavage observed by EPR, as well as the presence of a group that is highly susceptible to photooxidative attack. Therefore, the synergetic coupling of a spin trapping method with DFT calculations is shown to be a rapid and efficient method for providing unprecedented information on photochemical mechanisms. This approach will allow the design of LBG polymers without the need to trial the material within actual solar cell devices, an often long and costly screening procedure.

A new post-metallocene-ti catalyst with maltolate bidentade ligand: an investigation in heterogeneous polymerization reactions in different mesoporous supports

A new titanium catalyst easily synthesized from ethylmaltol bidentate chelator ligand was studied in homogeneous and heterogeneous ethylene polymerization. The dichlorobis(3-hydroxy-2-ethyl-4-pyrone)titanium(IV) complex was characterized by 1H and 13C NMR (nuclear magnetic resonance), UV-Vis and elemental analysis. Theoretical study by density functional theory (DFT) showed that the complex chlorines exhibit cis configuration, which is important for the activity in olefin polymerization. The complex was supported by two methods, direct impregnation or methylaluminoxane (MAO) pre-treatment, in five mesoporous supports: MCM-41 (micro and nano), SBA-15 and also the corresponding modified Al species. All the catalytic systems were active in ethylene polymerization and the catalytic activity was strongly influenced by the method of immobilization of the catalyst and the type of support.

First-Principles Calculations of the Electronic and Structural Properties of GaSb

In this paper, we carried out first-principles calculations in order to investigate the structural and electronic properties of the binary compound gallium antimonide (GaSb). This theoretical study was carried out using the Density Functional Theory within the plane-wave pseudopotential method. The effects ofexchange and correlation (XC) were treated using the functional Local Density Approximation (LDA), generalized gradient approximation (GGA): Perdew–Burke–Ernzerhof (PBE), Perdew-Burke-Ernzerhof revised for solids (PBEsol), Perdew-Wang91 (PW91), revised Perdew–Burke–Ernzerhof (rPBE), Armiento–Mattson 2005 (AM05) and meta-generalized gradient approximation (meta-GGA): Tao–Perdew– Staroverov–Scuseria (TPSS) and revised Tao–Perdew–Staroverov–Scuseria (RTPSS) and modified Becke-Johnson (MBJ). We calculated the densities of state (DOS) and band structure with different XC potentials identified and compared them with the theoretical and experimental results reported in the literature. It was discovered that functional: LDA, PBEsol, AM05 and RTPSS provide the best results to calculate the lattice parameters (a) and bulk modulus (B0); while for the cohesive energy (Ecoh), functional: AM05, RTPSS and PW91 are closer to the values obtained experimentally. The MBJ, Rtpss and AM05 values found for the band gap energy is slightly underestimated with those values reported experimentally.

Selective Sonogashira Couplings in 1,2,4,5-tetrazine derivatives

1,2,4,5-Tetrazines are six-membered heterocyclic compounds in which the four nitrogen atoms are displayed in a symmetric fashion. Their reactivity is quite different from other heterocyclic aromatic systems due to its unique electron-withdrawing character, comparable to tetra-nitrobenzene. 1 In particular, 1,2,4,5- tetrazines are known to take part in [4+2] inverse-Diels–Alder cycloaddition processes which efficiently lead to the construction of substituted pyridazine systems that are important in drug development and biomarker applications. 2 However, the electronic character of 1,2,4,5-tetrazines hampered the development of 3- ethynyl- and 3,6-diethynyl-1,2,4,5-tetrazine derivatives for molecular electronic applications, proved by the scarcity of examples found in the literature. 3 Herein, we describe the synthesis and characterization of two novel ethynyl-based 1,2,4,5-tetrazine derivatives. Synthesis of 3,6-(4-bromophenyl)-1,2,4,5-tetrazine precursor (1) was achieved in good yield by Pinner’s method, starting from 4-bromobenzonitrile. Despite its low solubility in common organic solvents, this precursor was found to react smoothly under typical Sonogashira coupling conditions to selectively afford the 3-ethynyl (2) and 3,6-diethynyl (3) protected derivatives (Figure 1). Reaction conditions were evaluated in order to provide the best yields and to promote selectivity of the mono- or disubstituted ethynyl derivatives. Finally, deprotection was achieved affording, in the case of compound 3, an unprecedented 3,6- diethynyl-1,2,4,5-tetrazine compound. Time-Dependent Density Functional Theory (TDDFT) calculations for both deprotected ethynyl derivatives were used to simulate electronic spectra. A deep knowledge of the relevant electronic transitions involved and quantitatively satisfactory results of the calculated electronic excitations in comparison with experimental data were obtained.

Theoretical and experimental studies of gold catalysts for glucose oxidation

The glucaric acid (GLA) has been identified as a “top value-added chemical from biomass” that can be employed for many uses; for instance, it could be a precursor of adipic acid, a monomer of Nylon-6,6. GLA can be synthetized by the oxidation of glucose (GLU), passing through the intermediate gluconic acid (GLO). In recent years, a new process has been sought to obtain GLA in an economic and environmental sustainable way, in order to replace the current use of HNO3 as a stoichiometric oxidant, or electrocatalysis and biochemical synthesis, which show several disadvantages. Thereby, this work is focused on the study of catalysts based on gold nanoparticles supported on activated carbon for the oxidation reaction of GLU to GLA using O2 as an oxidant agent and NaOH as base. The sol-immobilization method leads us to obtain small and well dispersed nanoparticles, characterized by UV-Vis, XRD and TEM techniques. Repeating the reaction on different batches of catalyst, both the synthesis and the reaction were confirmed to be reproducible. The effect of the reaction time feeding GLO as reagent was studied: the results show that the conversion of GLO increases as the reaction time increases; however, the yields of GLA and others increase up to 1 hour, and then they remain constant. In order to obtain information on the catalytic mechanism at the atomistic level, a computational study based on density functional theory and atomistic modeling of the gold nano-catalyst were performed. Highly symmetric (icosahedral and cubo-octahedral) and distorted Au55 nanoparticles have been optimized along with Au(111) and Au(100) surfaces. Distorted structures were found to be more stable than symmetrical ones due to relativistic effects. On these various models the adsorptions of various species involved in the catalysis have been studied, including OH- species, GLU and GLO. The study carried out aims to provide a method for approaching to the study of nanoparticellary catalytic systems.

Species Fractionation In Atomic Chains From Mechanically Stretched Alloys.

Bettini et al (2006 Nat. Nanotechnol. 1 182-5) reported the first experimental realization of linear atomic chains (LACs) composed of different atoms (Au and Ag). The different contents of Au and Ag were observed in the chains from what was found in the bulk alloys, which raises the question of what the wire composition is, if it is in equilibrium with a bulk alloy. In this work we address the thermodynamic driving force for species fractionation in LACs under tension, and we present the density-functional theory results for Ag-Au chain alloys. A pronounced stabilization of the wires with an alternating Ag-Au sequence is observed, which could be behind the experimentally observed Au enrichment in LACs from alloys with high Ag content.

A New Platinum Complex With Tryptophan: Synthesis, Structural Characterization, Dft Studies And Biological Assays In Vitro Over Human Tumorigenic Cells.

A new platinum(II) complex with the amino acid L-tryptophan (trp), named Pt-trp, was synthesized and characterized. Elemental, thermogravimetric and ESI-QTOF mass spectrometric analyses led to the composition [Pt(C11H11N2O2)2]⋅6H2O. Infrared spectroscopic data indicate the coordination of trp to Pt(II) through the oxygen of the carboxylate group and also through the nitrogen atom of the amino group. The (13)C CP/MAS NMR spectroscopic data confirm coordination through the oxygen atom of the carboxylate group, while the (15)N CP/MAS NMR data confirm coordination of the nitrogen of the NH2 group to the metal. Density functional theory (DFT) studies were applied to evaluate the cis and trans coordination modes of trp to platinum(II). The trans isomer was shown to be energetically more stable than the cis one. The Pt-trp complex was evaluated as a cytotoxic agent against SK-Mel 103 (human melanoma) and Panc-1 (human pancreatic carcinoma) cell lines. The complex was shown to be cytotoxic over the considered cells.

Synthesis, characterization, DFT and Td-dfT study of the [Fe(mnt)(L)(t-BuNC) 2] octahedral complex (L = phen, bipy)

FeBr2 has reacted with an equivalent of mnt2- (mnt = cis-1,2-dicyanoethylene-1,2-dithiolate) and the α-diimine L (L = 1,10'-phenantroline, 2,2'-bipyridine) in THF solution, and followed by adding of t-butyl-isocyanide to give [Fe(mnt)(L)(t-BuNC)2] neutral compound. The products were characterized by infrared, UV-visible and Mössbauer spectroscopy, besides thermogravimetric and conductivity data. The geometry in the equilibrium was calculated by the density functional theory and the electronic spectrum by the time-dependent. The experimental and theoretical results in good agreement have defined an octahedral geometry with two isocyanide neighbours. The π→π* intraligand electronic transition was not observed for cis-isomers in the near-IR spectral region.