The efficient computation of dynamic (time-dependent) molecular properties is a broad field with numerous applications in aiding molecular synthesis and design, with a particular preva- lence in spectroscopic predictions. Typical methods for computing the response of a molecu- lar system to an electromagnetic field (EMF) considers a quantum mechanical description of the molecule and a classical approximation for the EMF. Methods for describing light-matter interactions with high-accuracy electronic structure methods, such as coupled cluster (CC), are discussed, with a focus on improving the efficiency of such methods. The CC method suffers from high-degree polynomial scaling. In addition to the ground-state calculation, computing dynamic properties requires the description of sensitive excited-state effects. The cost of such methods often prohibits the accurate calculation of response prop- erties for systems of significant importance, such as large-molecule drug candidates or chiral species present in biological systems. While the literature is ripe with reduced-scaling meth- ods for CC ground-state calculations, considerably fewer approaches have been applied to excited-state properties, with even fewer still providing adequate results for realistic systems. This work presents three studies on the reduction of the cost of molecular property evalu- ations, in the hopes of closing this gap in the literature and widening the scope of current theoretical methods. There are two main ways of simulating time-dependent light-matter interactions: one may consider these effects in the frequency domain, where the response of the system to an EMF is computed directly; or, the response may be considered explicitly in the time domain, where wave function (or density) parameters can be propagated in time and examined in detail. Each methodology has unique advantages and computational bottlenecks. The first two studies focus on frequency-domain calculations, and employ fragmentation and machine- learning techniques to reduce the cost of single-molecule calculations or sets of calculations across a series of geometric conformations. The third study presents a novel application of the local correlation technique to real-time CC calculations, and highlights deficiencies and possible solutions to the approach. / Doctor of Philosophy / Theoretical chemistry plays a key role in connecting experimental results with physical inter- pretation. Paramount to the success of theoretical methods is the ability to predict molecular properties without the need for costly high-throughput synthesis, aiding in the determina- tion of molecular structure and the design of new materials. Light-matter interactions, which govern spectroscopic techniques, are particularly complicated, and sensitive to the theoreti- cal tools employed in their prediction. Compounding the issue of accuracy is one of efficiency — accurate theoretical methods typically incur steep scaling of computational cost (memory and processor time) with respect to the size of the system. An important aspect in improving the efficiency of these methods is understanding the nature of light-matter interactions at a quantum level. Many unanswered questions still remain, such as, "Can light-matter interactions be thought of as a sum of interactions be- tween smaller fragments of the system?" and "Can conventional methods of accelerating ground-state calculations be expected to perform well for spectroscopic properties?" The present work seeks to answer these questions through three studies, focusing on improving the efficiency of these techniques, while simultaneously addressing their fundamental flaws and providing reasonable alternatives.

Wave function based electronic structure methods have became a robust and reliable tool for the prediction and interpretation of the results of chemical experiments. However, they suffer from very steep scaling behavior with respect to an increase in the size of the system as well as very slow convergence of the correlation energy with respect to the basis set size. Thus these methods are limited to small systems of up to a dozen atoms. The first of these issues can be efficiently resolved by exploiting the local nature of electron correlation effects while the second problem is alleviated by the use of explicitly correlated R12/F12 methods. Since R12/F12 methods are central to this work, we start by reviewing their modern formulation. Next, we present the explicitly correlated second-order Mo ller-Plesset (MP2-F12) method in which all nontrivial post-mean-field steps are formulated with linear computational complexity in system size [Pavov{s}evi'c et al., {em J. Chem. Phys.} {bf 144}, 144109 (2016)]. The two key ideas are the use of pair-natural orbitals for compact representation of wave function amplitudes and the use of domain approximation to impose the block sparsity. This development utilizes the concepts for sparse representation of tensors described in the context of the DLPNO-MP2 method by Neese, Valeev and co-workers [Pinski et al., {em J. Chem. Phys.} {bf 143}, 034108 (2015)]. Novel developments reported here include the use of domains not only for the projected atomic orbitals, but also for the complementary auxiliary basis set (CABS) used to approximate the three- and four-electron integrals of the F12 theory, and a simplification of the standard B intermediate of the F12 theory that avoids computation of four-index two-electron integrals that involve two CABS indices. For quasi-1-dimensional systems (n-alkanes) the bigO{N} DLPNO-MP2-F12 method becomes less expensive than the conventional bigO{N^{5}} MP2-F12 for $n$ between 10 and 15, for double- and triple-zeta basis sets; for the largest alkane, C$_{200}$H$_{402}$, in def2-TZVP basis the observed computational complexity is $N^{sim1.6}$, largely due to the cubic cost of computing the mean-field operators. The method reproduces the canonical MP2-F12 energy with high precision: 99.9% of the canonical correlation energy is recovered with the default truncation parameters. Although its cost is significantly higher than that of DLPNO-MP2 method, the cost increase is compensated by the great reduction of the basis set error due to explicit correlation. We extend this formalism to develop a linear-scaling coupled-cluster singles and doubles with perturbative inclusion of triples and explicitly correlated geminals [Pavov{s}evi'c et al., {em J. Chem. Phys.} {bf 146}, 174108 (2017)]. Even for conservative truncation levels, the method rapidly reaches near-linear complexity in realistic basis sets; e.g., an effective scaling exponent of 1.49 was obtained for n-alkanes with up to 200 carbon atoms in a def2-TZVP basis set. The robustness of the method is benchmarked against the massively parallel implementation of the conventional explicitly correlated coupled-cluster for a 20-water cluster; the total dissociation energy of the cluster ($sim$186 kcal/mol) is affected by the reduced-scaling approximations by only $sim$0.4 kcal/mol. The reduced-scaling explicitly correlated CCSD(T) method is used to examine the binding energies of several systems in the L7 benchmark data set of noncovalent interactions. Additionally, we discuss a massively parallel implementation of the Laplace transform perturbative triple correction (T) to the DF-CCSD energy within density fitting framework. This work is closely related to the work by Scuseria and co-workers [Constans et al., {em J. Chem. Phys.} {bf 113}, 10451 (2000)]. The accuracy of quadrature with respect to the number of quadrature points has been investigated on systems of the 18-water cluster, uracil dimer and pentacene dimer. In the case of the 18-water cluster, the $mu text{E}_{text{h}}$ accuracy is achieved with only 3 quadrature points. For the uracil dimer and pentacene dimer, 6 or more quadrature points are required to achieve $mu text{E}_{text{h}}$ accuracy; however, binding energy of $<$1 kcal/mol is obtained with 4 quadrature points. We observe an excellent strong scaling behavior on distributed-memory commodity cluster for the 18-water cluster. Furthermore, the Laplace transform formulation of (T) performs faster than the canonical (T) in the case of studied systems. The efficiency of the method has been furthermore tested on a DNA base-pair, a system with more than one thousand basis functions. Lastly, we discuss an explicitly correlated formalism for the second-order single-particle Green's function method (GF2-F12) that does not assume the popular diagonal approximation, and describes the energy dependence of the explicitly correlated terms [Pavov{s}evi'c et al., {em J. Chem. Phys.} {bf 147}, 121101 (2017)]. For small and medium organic molecules the basis set errors of ionization potentials of GF2-F12 are radically improved relative to GF2: the performance of GF2-F12/aug-cc-pVDZ is better than that of GF2/aug-cc-pVQZ, at a significantly lower cost. / Ph. D.

Oxide dielectric breakdown is an old problem that has been studied over decades. It causes power dissipations and irreversible damage to the electronic devices. The aggressive downscaling of the device size exponentially increases the leakage current density, which also raises the risk of dielectric breakdown. It has been proposed that point defects, current leakages, impurity diffusions, etc. all contribute to the change of oxide chemical composition and ultimately lead to the dielectric breakdown. However, the conclusive cause and a clear understanding of the entire process of dielectric breakdown are still under debate. In this research, the electronic structure at metal/oxide interfaces is studied using first-principle calculations within the framework of Density Functional Theory (DFT) to investigate any possible key signature that would trigger the dielectric breakdown. A classical band alignment method, the Van de Walle method, is applied to the case study of the Al/crystal-SiO2 (Al/c-SiO2) interface. Point defects, such as oxygen vacancy (VO) and hydrogen impurity (IH), are introduced into the Al/c-SiO2 interface to study the effects on band offset and electronic structure caused by point defects at metal/oxide interfaces. It is shown that the bonding chemistry at metal/oxide interfaces, which is mainly ionic bond, polarizes the interface. It results in many interface effects such as the interface dipole, built-in voltage, band bending, etc. Charge density analysis also indicates that the interface can localize charge due to such ionic bonding. It is also found that VO at the interface traps metal electrons which closes the open -sp3 orbital. The analysis on local potential shows that the metal potential penetrates through a few layers of oxide starting from the interface, which metalizes the interfacial region and induces unoccupied states in the oxide band gap. In addition, it is shown that higher oxygen content at metal/oxide interfaces minimizes such metal potential invasion. In addition, an oxygen vacancy is created at multiple sites through the Al/c-SiO2 and Al/a-SiO2 interface systems, separately. The oxygen local pressure is also calculated before its removal using Quantum Stress Density theory. Correlations among electronic structure, stress density, and vacancy formation energy are found, which provide informative insights into the defect generation controlling and dielectric breakdown analysis. A new band alignment approach based on the projection of plane-waves (PWs) into the space-dependent atomic orbital (LCAO) basis is presented and tested against classical band offset methods -- the Van de Walle method. It is found that the new band alignment approach can provide a quantitative and reliable band alignment and can be applied to the heterojunctions consisting of amorphous materials. The new band alignment approach reveals the real-space dependency of the electronic structure at interfaces. In addition, it includes all interface effects, such as the interface dipole, built-in voltage, virtual oxide thinning, and band deformation, which cannot be derived using classical band offset methods. This new band alignment approach is applied to the case study of both the Al/amorphous-SiO2 (Al/a-SiO2) interface and the Al/c-SiO2. We have found that at extremely low dimensions, the reduction of the insulator character due to the virtual oxide thinning is a pure quantum effect. I highlight that the quantum tunneling current leakage is more critical than the decrease of the potential barrier height on the failure of the devices. / PHD

京都大学 / 新制・課程博士 / 博士(人間・環境学) / 甲第23975号 / 人博第1027号 / 新制||人||242(附属図書館) / 2022||人博||1027(吉田南総合図書館) / 京都大学大学院人間・環境学研究科相関環境学専攻 / (主査)教授 田部 勢津久, 教授 吉田 寿雄, 教授 中村 敏浩, 教授 田中 勝久 / 学位規則第4条第1項該当 / Doctor of Human and Environmental Studies / Kyoto University / DFAM

The field of organic electronics, in which -conjugated, organic molecules and polymers are used as the active components (e.g., semiconductor, light emitter/harvester, etc.), has lead to a number a number of key technological developments that have been founded within fundamental research disciplines. In the Dissertation that follows, the research involves the use of quantum-chemical techniques to elucidate fundamental aspects of both intermolecular and intramolecular electron-transfer processes in organic, -conjugated molecules. The Dissertation begins with an introduction and brief review of organic molecular systems used as electron-transport semiconducting materials in device applications and/or in the fundamental studies of intramolecular mixed-valence processes. This introductory material is then followed by a brief review of the electronic-structure methods (e.g., Hartree-Fock theory and Density Functional Theory) and electron-transfer theory (i.e., semiclassical Marcus theory) employed throughout the investigations. The next three Chapters deal with investigations related to the characterization of non-rigid, -conjugated molecular systems that have amorphous solid-state properties used as the electron-transport layer in organic electronic and optoelectronic devices. Chapters 3 and 4 involve studies of silole- (silacyclopentadiene)-based materials that possess attractive electronic and optical properties in the solid state. Chapter 5 offers a preliminary study of dioxaborine-based molecular structures as electron-transport systems. In Chapters 6 8, the focus of the work shifts to investigations of organic mixed-valence systems. Chapter 6 centers on the examination of tetraanisylarylenediamine systems where the inter-redox site distances are approximately equal throughout the series. Chapter 7 examines the bridge-length dependence of the geometric structure, charge-(de)localization, and electronic coupling for a series of vinylene- and phenylene-vinylene-bridged bis-dianisylamines. In Chapter 8, the role of symmetric vibrations in the delocalization of the excess charge is studied in a dioxaborine radical-anion and a series of radical-cation bridged-bisdimethylamines. Finally, Chapter 9 provides a synopsis of the work and goals for future consideration.

Electronic structure calculations constitute a valuable tool to predict the properties of materials. In this study we propose an efficient scheme to study correlated electron systems with essentially only one free parameter, the screening length of the Coulomb potential. A general reformulation of the exchange energy of the correlated electron shell is combined with this method in order to analyze the calculations. The results are interpreted in terms of different polarization channels, due to different multipoles. The method is applied to various actinide compounds, in order to increase the understanding of the complicate behaviour of 5f electrons in these systems. We studied the non-magnetic phase of δ-Pu, where the spin polarization is taken over by a spin-orbit-like term that does not break the time reversal symmetry. We also find that a non-trivial high multipole of the magnetization density, the triakontadipole, constitutes the ordering parameter in the mysterious hidden order phase of the heavy-fermion superconductor URu2Si2. This type of multipolar ordering is also found to play an essential role in the hexagonal-based superconductors UPd2Al3,  UNi2Al3 and UPt3 and in the dioxide insulators UO2, NpO2 and PuO2. The triakontadipole moments are also present in all magnetic actinides we considered, except for Cm. These results led us to formulate a new set of rules for the ground state of a system, that are valid in presence of strong spin-orbit coupling interaction instead of those of Hund; the Katt's rules. Finally, we applied our method to a new class of high-Tc superconductors, the Fe-pnictides, where the Fe 3d electrons are moderately correlated. In these materials we obtain the stabilization of a low spin moment solution, in agreement with experiment, over a large moment solution, due to the gain in exchange energy in the formation of large multipoles of the spin magnetization density. / Felaktigt tryckt som Digital Comprehensive Summaries of Uppsala Dissertations from the Faculty of Science and Technology 705

In disordered alloys, atoms belonging to the same chemical element will exhibit different environments. This leads to variations in the atoms’ local electronic structures, which in turn leads to variations in the binding energies of their core levels. These binding energies can be measured experimentally using core level X-ray photoelectron spectroscopy (XPS). Therefore, in theory at least, core level XPS can be used to resolve different environments in alloys. However, to make this a reality one must understand how an atom’s local electronic structure, and hence the binding energies of its core levels, are affected by local environment. In this thesis, two simple phenomenological models are explored which purport to correctly describe the local electronic structure of disordered alloys. The first model which we consider has its roots in chemical intuition; specifically, the notion that pairs of unlike atoms, i.e. atoms belonging to different chemical elements, transfer a certain quantity of charge, while like atoms do not. Using this model - known as the optimised linear charge model (OLCM) - the relationship between an atom’s local electronic structure, core level binding energies, and its environment is explored in detail, both in the bulk of disordered alloys and near their surfaces. As well as ‘homogeneous’ disordered alloys, in which the concentrations of the alloy’s constituent elements are the same throughout the entire alloy, various ‘inhomogeneous’ disordered alloy systems are considered. These include alloys exhibiting surface segregation - in which the concentrations at the surface differ from those in the bulk - as well as interfaces between two metals with various levels of intermixing. The results of our investigation of bulk inhomogeneous alloys are compared to analogous ab initio results, which confirms the model’s viability as a tool for rationalising the relationship between local electronic structure, core level binding energies, and environment. More generally, our results also reveal a number of interesting new phenomena. Firstly, the widths of spectra in inhomogeneous disordered alloys are significantly larger in some cases than is possible in any analogous homogeneous disordered alloy. Secondly, differences between the concentrations of each element at the surface and deep within the bulk cause a shift in the work function of the alloy under consideration. The latter results in qualitatively different trends than one would expect if this phenomenon was ignored, and prompts an alternative interpretation of the results of a recent experimental study. The second model which we consider is a particular case of the charge-excess functional model, in which the realised charges on all atoms are those which minimise a particular expression for the total energy of the system, and whose accuracy has been well established. The underlying assumptions and properties of this model are explored in detail, adding insight into the nature of the screening and inter-atomic interactions in disordered alloys. The model is shown to be equivalent to the OLCM for the case of binary alloys, and can therefore be considered to be the generalisation of the OLCM for alloys containing more than two chemical elements. The model is also used to derive analytical expressions for various physical quantities for any alloy, including the width of core level XPS spectra and the Madelung energy. These expressions are then used to investigate how the physical quantities to which they pertain vary with the concentrations of each element in a homogeneous disordered alloy consisting of three elements. Among other things, it was observed that the width of the core level XPS spectra is maximised when the concentrations of the two elements in the alloy with the largest electronegativity difference have equal concentrations, while the remaining element has a vanishing concentration.

The localized model of a chemical bond has had a long and prominent role in chemistry, but situations of extended charge delocalization and dipole effects remain topics in need of greater understanding. Both orbital delocalization in isolated molecules and induced molecular dipoles in condensed phases serve to move electron density and influence the chemical and physical properties of a system. This dissertation studies these aspects of electronic structure for selected organic, inorganic, and organometallic systems by means of electronic structure calculations and photoelectron spectroscopy, which is well-suited for studying both intramolecular and intermolecular effects by providing a direct probe of orbital energies and characters. Photoelectron spectra of P₄ and AsP₃ reveal differences in the molecular symmetry and cationic state effects between the two molecules in Chapter 3. Despite these differences, AsP₃ is found to have electron delocalization and vibrational structures that are comparable to P₄. A similar study of the delocalized -system of 2H-1,2,3-triazole in Chapter 4 relates the vibrational structure in photoelectron spectroscopy data to a series of Rydberg excitations in the vacuum UV photoabsorption spectrum. Chapters 5 and 6 examine extended electronic structures in organometallic complexes. The electron delocalization and charge transfer between two Ru centers along a bridging ethynediyl ligand is studied in [CpRu(CO)₂]₂(μC≡C). Details of the Ru-alkynyl interaction were explored by comparing the spectra of CpRu(CO)₂C≡CMe with CpRu(CO)₂Cl, including the -backbonding ability of alkynyl ligands. Chapter 6 moves from the realm of intramolecular effects to intermolecular interactions to understand how surrounding media affect electronic properties of molecules. The reversal of ionization energies between the gas and solid phases of M(CO)₄dmpe and M(CO)₄dppe, where M = Mo, W, is explored with photoelectron spectroscopy. The surrounding molecular environment stabilizes the cation, resulting in this reversal that extends to core ionization energies. The variety of systems presented illustrates the wide applicability of photoelectron spectroscopy and computations to different electronic structure studies, including how gas phase results can be related to condensed phase studies. This work continues the progress of photoelectron spectroscopy from small molecules to larger molecular systems and even further to bulk systems.

No description available.

No description available.