Metastable de-excitation spectroscopy and density functional theory study of the selective oxidation of crotyl alcohol over Pd(111)

The extremely surface sensitive technique of metastable de-excitation spectroscopy (MDS) has been utilized to probe the bonding and reactivity of crotyl alcohol over Pd(111) and provide insight into the selective oxidation pathway to crotonaldehyde. Auger de-excitation (AD) of metastable He (23S) atoms reveals distinct features associated with the molecular orbitals of the adsorbed alcohol, corresponding to emission from the hydrocarbon skeleton, the O n nonbonding, and C═C π states. The O n and C═C π states of the alcohol are reversed when compared to those of the aldehyde. Density functional theory (DFT) calculations of the alcohol show that an adsorption mode with both C═C and O bonds aligned somewhat parallel to the surface is energetically favored at a substrate temperature below 200 K. Density of states calculations for such configurations are in excellent agreement with experimental MDS measurements. MDS revealed oxidative dehydrogenation of crotyl alcohol to crotonaldehyde between 200 and 250 K, resulting in small peak shifts to higher binding energy. Intramolecular changes lead to the opposite assignment of the first two MOs in the alcohol versus the aldehyde, in accordance with DFT and UPS studies of the free molecules. Subsequent crotonaldehyde decarbonylation and associated propylidyne formation above 260 K could also be identified by MDS and complementary theoretical calculations as the origin of deactivation and selectivity loss. Combining MDS and DFT in this way represents a novel approach to elucidating surface catalyzed reaction pathways associated with a “real-world” practical chemical transformation, namely the selective oxidation of alcohols to aldehydes.

Generalized atoms in molecules approach (GAIM approach)

This thesis involves two parts. The first is a new-proposed theoretical approach called generalized atoms in molecules (GAIM). The second is a computational study on the deamination reaction of adenine with OH⁻/nH₂O (n=0, 1, 2, 3) and 3H₂O. The GAIM approach aims to solve the energy of each atom variationally in the first step and then to build the energy of a molecule from each atom. Thus the energy of a diatomic molecule (A-B) is formulated as a sum of its atomic energies, EA and EB. Each of these atomic energies is expressed as, EA = Hᴬ + Vₑₑᴬᴬ + 1/2Vₑₑᴬ<>ᴮ EB = Hᴮ + Vₑₑᴮᴮ + 1/2Vₑₑᴬ<>ᴮ where; Hᴬ and Hᴮ are the kinetic and nuclear attraction energy of electrons of atoms A and B, respectively; Vₑₑᴬᴬ and Vₑₑᴮᴮ are the interaction energy between the electrons on atoms A and B, respectively; and Vₑₑᴬ<>ᴮ is the interaction energy between the electrons of atom A with the electrons of atom B. The energy of the molecule is then minimized subject to the following constraint, |ρA(r)dr + |ρB(r)dr = N where ρA(r) and ρB(r) are the electron densities of atoms A and B, respectively, and N is the number of electrons. The initial testing of the performance of GAIM was done through calculating dissociation curves for H₂, LiH, Li₂, BH, HF, HCl, N₂, F₂, and Cl₂. The numerical results show that GAIM performs very well with H₂, LiH, Li₂, BH, HF, and HCl. GAIM shows convergence problems with N₂, F₂, and Cl₂ due to difficulties in reordering the degenerate atomic orbitals Pₓ, Py, and Pz in N, F, and Cl atoms. Further work for the development of GAIM is required. Deamination of adenine results in one of several forms of premutagenic lesions occurring in DNA. In this thesis, mechanisms for the deamination reaction of adenine with OH⁻/nH₂O, (n = 0, 1, 2, 3) and 3H₂O were investigated. HF/6-31G(d), B3LYP/6-31G(d), MP2/6-31G(d), and B3LYP/6-31+G(d) levels of theory were employed to optimize all the geometries. Energies were calculated at the G3MP2B3 and CBS-QB3 levels of theory. The effect of solvent (water) was computed using the polarizable continuum model (PCM). Intrinsic reaction coordinate (IRC) calculations were performed for all transition states. Five pathways were investigated for the deamination reaction of adenine with OH⁻/nH₂O and 3H₂O. The first four pathways (A-D) begin with by deprotonation at the amino group of adenine by OH⁻, while pathway E is initiated by tautomerization of adenine. For all pathways, the next two steps involve the formation of a tetrahedral intermediate followed by dissociation to yield products via a 1,3-hydrogen shift. Deamination with a single OH⁻ has a high activation barrier (190 kJ mol⁻¹ using G3MP2B3 level) for the rate-determining step. Addition of one water molecule reduces this barrier by 68 kJ mol⁻¹ calculated at G3MP2B3 level. Adding more water molecules decreases the overall activation energy of the reaction, but the effect becomes smaller with each additional water molecule. The most plausible mechanism is pathway E, the deamination reaction of adenine with 3H₂O, which has an overall G3MP2B3 activation energy of 139 and 137 kJ mol⁻¹ in the gas phase and PCM, respectively. This barrier is lower than that for the deamination with OH⁻/3H₂O by 6 and 2 kJ mol⁻¹ in the gas phase and PCM, respectively.

Electronic Instability and Anharmonicity in SnSe

The binary compound SnSe exhibits record high thermoelectric performance, largely because of its very low thermal conductivity. The origin of the strong phonon anharmonicity leading to the low thermal conductivity of SnSe is investigated through first-principles calculations of the electronic structure and phonons. It is shown that a Jahn-Teller instability of the electronic structure is responsible for the high-temperature lattice distortion between the Cmcm and Pnma phases. The coupling of phonon modes and the phase transition mechanism are elucidated, emphasizing the connection with hybrid improper ferroelectrics. This coupled instability of electronic orbitals and lattice dynamics is the origin of the strong anharmonicity causing the ultralow thermal conductivity in SnSe. Exploiting such bonding instabilities to generate strong anharmonicity may provide a new rational to design efficient thermoelectric materials.

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.

Level-resolved quantum statistical theory of electron capture into many-electron compound resonances in highly charged ions

The strong mixing of many-electron basis states in excited atoms and ions with open f shells results in very large numbers of complex, chaotic eigenstates that cannot be computed to any degree of accuracy. Describing the processes which involve such states requires the use of a statistical theory. Electron capture into these “compound resonances” leads to electron-ion recombination rates that are orders of magnitude greater than those of direct, radiative recombination and cannot be described by standard theories of dielectronic recombination. Previous statistical theories considered this as a two-electron capture process which populates a pair of single-particle orbitals, followed by “spreading” of the two-electron states into chaotically mixed eigenstates. This method is similar to a configuration-average approach because it neglects potentially important effects of spectator electrons and conservation of total angular momentum. In this work we develop a statistical theory which considers electron capture into “doorway” states with definite angular momentum obtained by the configuration interaction method. We apply this approach to electron recombination with W²⁰⁺, considering 2×10⁶ doorway states. Despite strong effects from the spectator electrons, we find that the results of the earlier theories largely hold. Finally, we extract the fluorescence yield (the probability of photoemission and hence recombination) by comparison with experiment.

Morphology and material stability in polymer solar cells

Polymer solar cells are promising in that they are inexpensive to produce, and due to their mechanical flexibility have the potential for use in applications not possible for more traditional types of solar cells. The performance of polymer solar cells depends strongly on the distribution of electron donor and acceptor material in the active layer. Understanding the connection between morphology and performance as well as how to control the morphology, is therefore of great importance. Furthermore, improving the lifetime of polymer solar cells has become at least as important as improving the efficiency.   In this thesis, the relation between morphology and solar cell performance is studied, and the material stability for blend films of the thiophene-quinoxaline copolymer TQ1 and the fullerene derivatives PCBM and PC70BM. Atomic force microscopy (AFM) and scanning transmission X-ray microscopy (STXM) are used to investigate the lateral morphology, secondary ion mass spectrometry (SIMS) to measure the vertical morphology and near-edge X-ray absorption fine structure (NEXAFS) spectroscopy to determine the surface composition. Lateral phase-separated domains are observed whose size is correlated to the solar cell performance, while the observed TQ1 surface enrichment does not affect the performance. Changes to the unoccupied molecular orbitals as a result of illumination in ambient air are observed by NEXAFS spectroscopy for PCBM, but not for TQ1. The NEXAFS spectrum of PCBM in a blend with TQ1 changes more than that of pristine PCBM. Solar cells in which the active layer has been illuminated in air prior to the deposition of the top electrode exhibit greatly reduced electrical performance. The valence band and absorption spectrum of TQ1 is affected by illumination in air, but the effects are not large enough to account for losses in solar cell performance, which are mainly attributed to PCBM degradation at the active layer surface.

Calculated ionization potentials of the linear alkanes

The equivalent orbital (EO) method is investigated and used for predicting outer and inner ionization potentials of the linear alkanes. The calculated ionization potentials are in good agreement with those observed in photoelectron spectra provided that a set of 12 parameters is used in the theory. An optimization technique is used to find the best values for thle parameters and a single transferable parameter set can be found which is applicable to all the n-alkanes. A good fit to the experimental results can only be obtained if the uppermost molecular orbital of the n-alkanes is an antisymmetrical orbital built up from CH equivalent orbitals.

Coupling of magnetic, strain and electric polarization fields in the structure of multiferroic material

The comprehensive study on the coupling of magnetism, electrical polarization and the crystalline lattice with the off-stoichiometric effects in self-doped multiferroic hexagonal h-LuMnxO3±δ (0.92≤x≤1.12) ceramic oxides was carried out for the PhD work. There is a complex coupling of the three ferroic degrees. The cancelation of the magnetic moments of ions in the antiferromagnetic order, electric polarization with specific vortex/antivortex topology and lattice properties have pushed researchers to find out ways to disclose the underlying physics and chemistry of magneto-electric and magneto-elastic couplings of h-RMnO3 multiferroic materials. In this research work, self-doping of Lu-sites or Mn-sites of h-LuMnxO3±δ ceramics prepared via solid state route was done to pave a way for deeper understanding of the antiferromagnetic transition, the weak ferromagnetism often reported in the same crystalline lattices and the ferroelectric properties coupled to the imposed lattice changes. Accordingly to the aim of the PhD thesis, the objectives set for the sintering study in the first chapter on experimental results were two. First, study of sintering off-stoichiometric samples within conditions reported in the bibliography and also extracted from the phase diagrams of the LuMnxO3±δ, with a multiple firings ending with a last high temperature step at 1300ºC for 24 hours. Second, explore longer annealing times of up to 240 hours at the fixed temperature of 1300 ºC in a search for improving the properties of the solid solution under study. All series of LuMnxO3±δ ceramics for each annealing time were characterized to tentatively build a framework enabling comparison of measured properties with results of others available in literature. XRD and Rietveld refinement of data give the evolution the lattice parameters as a function to x. Shrinkage of the lattice parameters with increasing x values was observed, the stability limit of the solid solution being determined by analysis of lattice parameters. The evolution of grain size and presence of secondary phases have been investigated by means of TEM, SEM, EDS and EBSD techniques. The dependencies of grain growth and regression of secondary phases on composition x and time were further characterized. Magnetic susceptibility of samples and magnetic irreversibility were extensively examined in the present work. The dependency of magnetic susceptibility, Neel ordering transition and important magnetic parameters are determined and compared to observation in other multiferroics in the following chapter of the thesis. As a tool of high sensitivity to detect minor traces of the secondary phase hausmannite, magnetic measurements are suggested for cross-checking of phase diagrams. Difficulty of previous studies on interpreting the magnetic anomaly below 43 K in h-RMnO3 oxides was discussed and assigned to the Mn3O4 phase, with supported of the electron microscopy. Magneto-electric coupling where AFM ordering is coupled to dielectric polarization is investigated as a function of x and of sintering condition via frequency and temperature dependent complex dielectric constant measurements in the final chapter of the thesis. Within the limits of solid solubility, the crystalline lattice of off-stoichiometric ceramics was shown to preserve the magneto-electric coupling at TN. It represents the first research work on magneto-electric coupling modified by vacancy doping to author’s knowledge. Studied lattices would reveal distortions at the atomic scale imposed by local changes of x dependent on sintering conditions which were widely inspected by using TEM/STEM methods, complemented with EDS and EELS spectroscopy all together to provide comprehensive information on cross coupling of distortions, inhomogeneity and electronic structure assembled and discussed in a specific chapter. Internal interfaces inside crystalline grains were examined. Qualitative explanations of the measured magnetic and ferroelectric properties were established in relation to observed nanoscale features of h-LuMnxO3±δ ceramics. Ferroelectric domains and topological defects are displayed both in TEM and AFM/PFM images, the later technique being used to look at size, distribution and switching of ferroelectric domains influenced by vacancy doping at the micron scale bridging to complementary TEM studies on the atomic structure of ferroelectric domains. In support to experimental study, DFT simulations using Wien2K code have been carried out in order to interpret the results of EELS spectra of O K-edge and to obtain information on the cation hybridization to oxygen ions. The L3,2 edges of Mn is used to access the oxidation state of the Mn ions inside crystalline grains. In addition, rehybridization driven ferroelectricity is also evaluated by comparing the partial density of states of the orbitals of all ions of the samples, also the polarization was calculated and correlated to the off-stoichiometric effect.

On the rotation of co-orbital bodies in eccentric orbits

We investigate the resonant rotation of co-orbital bodies in eccentric and planar orbits. We develop a simple analytical model to study the impact of the eccentricity and orbital perturbations on the spin dynamics. This model is relevant in the entire domain of horseshoe and tadpole orbit, for moderate eccentricities. We show that there are three different families of spin-orbit resonances, one depending on the eccentricity, one depending on the orbital libration frequency, and another depending on the pericenter's dynamics. We can estimate the width and the location of the different resonant islands in the phase space, predicting which are the more likely to capture the spin of the rotating body. In some regions of the phase space the resonant islands may overlap, giving rise to chaotic rotation.

Interactions of pyridinium, pyrrolidinium or piperidinium based ionic liquids with water: measurements and COSMO-RS modelling

Looking for a better knowledge concerning water and ionic liquids (ILs) interactions, a systematic study of the activity coefficients of water in pyridinium, pyrrolidinium and piperidinium-based ILs at 298.2 K is here presented based on water activity measurements. Additionally, the study of the structural effects of the pyridinium-based cation is also pursued. The results show that non-aromatic ILs are interacting more with water than aromatic ones, and among the ortho, meta and para isomers of 1-butyl-methylpyridinium chloride, the ortho position confers a more hydrophilic character to that specific IL. The physicalchemistry of the solutions was interpreted based on dissociation constants, natural bond orbitals and excess enthalpies providing a sound basis for the interpretation of the experimental observations. These results show that hydrogen bonding controls the behavior of these systems, being the anion-water one of the most relevant interactions, but modulated by the anionecation interactions.

Direct experimental evidence for a multiparticle-hole ground state configuration of deformed Mg-33

The first direct experimental evidence of a multiparticle-hole ground state configuration of the neutron-rich Mg-33 isotope has been obtained via intermediate energy (400 A MeV) Coulomb dissociation measurement. The major part similar to(70 +/- 13)% of the cross section is observed to populate the excited states of Mg-32 after the Coulomb breakup of Mg-33. The shapes of the differential Coulomb dissociation cross sections in coincidence with different core excited states favor that the valence neutron occupies both the s(1/2) and p(3/2) orbitals. These experimental findings suggest a significant reduction and merging of sd-pf shell gaps at N similar to 20 and 28. The ground state configuration of Mg-33 is predominantly a combination of Mg-32(3.0,3.5MeV; 2(-), 1(-)) circle times nu(s1/2), Mg-32(2.5MeV; 2(+)) circle times nu(p3/2), and Mg-32(0; 0(+)) circle times nu(p3/2). The experimentally obtained quantitative spectroscopic information for the valence neutron occupation of the s and p orbitals, coupled with different core states, is in agreement with Monte Carlo shell model (MCSM) calculation using 3 MeV as the shell gap at N = 20.

Double Hybrid Functionals and the Π-System Bond Length Alternation Challenge: Rivaling Accuracy of Post-HF Methods

Predicting accurate bond length alternations (BLAs) in long conjugated oligomers has been a significant challenge for electronic-structure methods for many decades, made particularly important by the close relationships between BLA and the rich optoelectronic properties of π-delocalized systems. Here, we test the accuracy of recently developed, and increasingly popular, double hybrid (DH) functionals, positioned at the top of Jacobs Ladder of DFT methods of increasing sophistication, computational cost, and accuracy, due to incorporation of MP2 correlation energy. Our test systems comprise oligomeric series of polyacetylene, polymethineimine, and polysilaacetylene up to six units long. MP2 calculations reveal a pronounced shift in BLAs between the 6-31G(d) basis set used in many studies of BLA to date and the larger cc-pVTZ basis set, but only modest shifts between cc-pVTZ and aug-cc-pVQZ results. We hence perform new reference CCSD(T)/cc-pVTZ calculations for all three series of oligomers against which we assess the performance of several families of DH functionals based on BLYP, PBE, and TPSS, along with lower-rung relatives including global- and range-separated hybrids. Our results show that DH functionals systematically improve the accuracy of BLAs relative to single hybrid functionals. xDH-PBE0 (N4 scaling using SOS-MP2) emerges as a DH functional rivaling the BLA accuracy of SCS-MP2 (N5 scaling), which was found to offer the best compromise between computational cost and accuracy the last time the BLA accuracy of DFT- and wave function-based methods was systematically investigated. Interestingly, xDH-PBE0 (XYG3), which differs to other DHs in that its MP2 term uses PBE0 (B3LYP) orbitals that are not self-consistent with the DH functional, is an outlier of trends of decreasing average BLA errors with increasing fractions of MP2 correlation and HF exchange.

Materials design for ambipolar devices: tuning orbital energetics in oligothiophene-naphthalimides semiconductors

Ambipolar organic field-effect transistors (OFETs), which can efficiently transport both holes and electrons, using a single type of electrode, are currently of great interest due to their possible applications in complementary metal oxide semiconductor (CMOS)-like circuits, sensors, and in light-emitting transistors. Several theoretical and experimental studies have argued that most organic semiconductors should be able to transport both types of carrier, although typically unipolar behavior is observed. One factor that can compromise ambipolar transport in organic semiconductors is poor solid state overlap between the HOMO (p-type) or LUMO (n-type) orbitals of neighboring molecules in the semiconductor thin film. In the search of low-bandgap ambipolar materials, where the absence of skeletal distortions allows closer intermolecular π-π stacking and enhanced intramolecular π-conjugation, a new family of oligothiophene-naphthalimide assemblies have been synthesized and characterized, in which both donor and acceptor moieties are directly conjugated through rigid linkers. In previous works we found that oligothiophene-napthalimide assemblies connected through amidine linkers (NDI derivates) exhibit skeletal distortions (50-60º) arising from steric hindrance between the carbonyl group of the arylene core and the sulphur atom of the neighbored thiophene ring (see Figure 1). In the present work we report novel oligo- and polythiophene–naphthalimide analogues NAI-3T, NAI-5T and poly-NAI-8C-3T, in which the connections of the amidine linkage have been inverted in order to prevent steric interactions. Thus, the nitrogen atoms are directly connected to the naphthalene moiety in NAI derivatives while they were attached directly to the thiophene moiety in the previously investigated NDI-3T and NDI-5T. In Figure 1 is depicted the calculated molecular structure of NAI-3T together with that of NDI-3T showing how the steric interactions are not present in the novel NAI derivative. The planar skeletons in these new family induce higher degree of crystallinity and the carrier charge transport can be switched from n-type to ambipolar behaviour. The highest FET performance is achieved for vapor-deposited films of NAI-3T with mobilities of 1.95x10-4cm2V-1s-1 and 2.00x10-4cm2V-1s-1 for electrons and holes, respectively. Finally, these planar semiconductors are compared with their NDI derivates analogues, which exhibit only n-type mobility, in order to understand the origin of the ambipolarity in this new series of molecular semiconductors.

Modern Computational Chemistry Methods for Prediction of Ground- and Excited-State Properties in Open-Shell Systems

The main goal of the research presented in this work is to provide some important insights about computational modeling of open-shell species. Such projects are: the investigation of the size-extensivity error in Equation-of-Motion Coupled Cluster methods, the analysis of the Long-Range corrected scheme in predicting UV-Vis spectra of Cu(II) complexes with the 4-imidazole acetate and its ethylated derivative, and the exploration of the importance of choosing a proper basis set for the description of systems such as the lithium monoxide anion. The most significant findings of this research are: (i) The contribution of the left operator to the size-extensivity error of the CR-EOMCC(2,3) approach, (ii) The cause of d-d shifts when varying the range-separation parameter and the amount of the exact exchange arising from the imbalanced treatment of localized vs. delocalized orbitals via the "tuned" CAM-B3LYP* functional, (iii) The proper acidity trend of the first-row hydrides and their lithiated analogs that may be reversed if the basis sets are not correctly selected.