Ring currents in azulene

We propose a self consistent polarisable ion tight binding theory for the study of push-pull processes in aromatic molecules. We find that the method quantitatively reproduces ab initio calculations of dipole moments and polarisability. We apply the scheme in a simulation which solves the time dependent Schroedinger equation to follow the relaxation of azulene from the second excited to the ground states. We observe rather spectacular oscillating ring currents which we explain in terms of interference between the HOMO and LUMO states.

Lowest energy potentials of Li2 + F accessible at ultracold temperatures

Ab initio calculations for the strongly exoergic Li2+F harpoon reaction are presented using density-functional theory, complete active space self-consistent field, and multireference configuration interaction methods to argue that this reaction would be an ideal candidate for investigation with ultracold molecules. The lowest six states are calculated with the aug-correlation-consistent polarized valence triple-zeta basis set and at least two can be accessed by a ground rovibronic Li2 molecule with zero collision energy at all reaction geometries. The large reactive cross section (characteristic of harpoon reactions) and chemiluminescent products are additional attractive features of these reactions.

Prospects for ultracold carbon via charge exchange reactions and laser cooled carbides

Strategies to produce an ultracold sample of carbon atoms are explored and assessed with the help of quantum chemistry. After a brief discussion of the experimental difficulties using conventional methods, two strategies are investigated. The first attempts to exploit charge exchange reactions between ultracold metal atoms and sympathetically cooled C+ ions. Ab initio calculations including electron correlation have been conducted on the molecular ions [LiC]+ and [BeC]+ to determine whether alkali or alkaline earth metals are a suitable buffer gas for the formation of C atoms but strong spontaneous radiative charge exchange ensure they are not ideal. The second technique involves the stimulated production of ultracold C atoms from a gas of laser cooled carbides. Calculations on LiC suggest that the alkali carbides are not suitable but the CH radical is a possible laser cooling candidate thanks to very favourable Frank-Condon factors. A scheme based on a four pulse STIRAP excitation pathway to a Feshbach resonance is outlined for the production of atomic fragments with near zero centre of mass velocity.

Singly, doubly and triply resonant multiphoton processes involving autoionizing states in magnesium

Using the R-matrix Floquet theory we have carried out non-perturbative, ab initio one- and two-colour calculations of the multiphoton ionization of magnesium with the laser frequencies chosen such that the initial state of the atom is resonantly coupled with autoionizing resonances of the atom. Good agreement is obtained with previous calculations in the low-intensity regimes. The single-photon ionization from the 3s3p P excited state of magnesium has been studied in the vicinity of the 3p S autoionizing resonance at non-perturbative laser intensities. Laser-induced degenerate states (LIDS) are observed for modest laser intensities. By adding a second laser which resonantly couples the 3p S = and 3p3d P autoionizing levels, we show that, due to the small width of the 3p3d P state, LIDS occur between this state and the 3s3p P state at intensities of the first laser below 10 W cm . We next investigate the case in which the first laser induces a resonant two-photon coupling between the ground state and the 3p S autoionizing state, while the second laser again resonantly couples the respective 3p S and 3p3d P autoionizing states. At weak intensities, our calculations compare favourably with recent experimental data and calculations. We show that when the intensity of the first laser is increased, the effect of an additional autoionizing state, the 4s5s S state, becomes significant. This state is coupled to the 3p3d P autoionizing level by one photon, inducing a triply resonant processes. We show that LIDS occur among the three autoionizing levels and we discuss their effect on the decay rate of the ground state. We consider dressed two- and three-level atoms which can be used to model the results of our calculations.

Exciton-Plasmon States in Nanoscale Materials: Breakdown of the Tamm-Dancoff Approximation

Within the Tamm-Dancoff approximation, ab initio approaches describe excitons as packets of electron-hole pairs propagating only forward in time. However, we show that in nanoscale materials excitons and plasmons hybridize, creating exciton-plasmon states where the electron-hole pairs oscillate back and forth in time. Then, as exemplified by the trans-azobenzene molecule and the carbon nanotubes, the Tamm-Dancoff approximation yields errors larger than the accuracy claimed in ab initio calculations. Instead, we propose a general and efficient approach that avoids the Tamm-Dancoff approximation, correctly describes excitons, plasmons, and exciton-plasmon states, and provides a good agreement with experimental results.

Strong electronic correlation in the hydrogen chain: A variational Monte Carlo study

In this paper, we report a fully ab initio variational Monte Carlo study of the linear and periodic chain of hydrogen atoms, a prototype system providing the simplest example of strong electronic correlation in low dimensions. In particular, we prove that numerical accuracy comparable to that of benchmark density-matrix renormalization-group calculations can be achieved by using a highly correlated Jastrow-antisymmetrized geminal power variational wave function. Furthermore, by using the so-called "modern theory of polarization" and by studying the spin-spin and dimer-dimer correlations functions, we have characterized in detail the crossover between the weakly and strongly correlated regimes of this atomic chain. Our results show that variational Monte Carlo provides an accurate and flexible alternative to highly correlated methods of quantum chemistry which, at variance with these methods, can be also applied to a strongly correlated solid in low dimensions close to a crossover or a phase transition.

Universal tight binding model for chemical reactions in solution and at surfaces. II. Water

A revised water model intended for use in condensed phase simulations in the framework of the self consistent polarizable ion tight binding theory is constructed. The model is applied to water monomer, dimer, hexamers, ice, and liquid, where it demonstrates good agreement with theoretical results obtained by more accurate methods, such as DFT and CCSD(T), and with experiment. In particular, the temperature dependence of the self diffusion coefficient in liquid water predicted by the model, closely reproduces experimental curves in the temperature interval between 230 K and 350 K. In addition, and in contrast to standard DFT, the model properly orders the relative densities of liquid water and ice. A notable, but inevitable, shortcoming of the model is underestimation of the static dielectric constant by a factor of two. We demonstrate that the description of inter and intramolecular forces embodied in the tight binding approximation in quantum mechanics leads to a number of valuable insights which can be missing from ab initio quantum chemistry and classical force fields. These include a discussion of the origin of the enhanced molecular electric dipole moment in the condensed phases, and a detailed explanation for the increase of coordination number in liquid water as a function of temperature and compared with ice-leading to insights into the anomalous expansion on freezing. The theory holds out the prospect of an understanding of the currently unexplained density maximum of water near the freezing point. 

Site symmetry dependence of repulsive interactions between chemisorbed oxygen atoms on Pt{100}-(1x1)

Ab initio total energy calculations using density functional theory with the generalized gradient approximation have been performed for the chemisorption of oxygen atoms on a Pt{100}-(1 x 1) slab. Binding energies for the adsorption of oxygen on different high-symmetry sites are presented. The bridge site is the most stable at a coverage of 0.5 ML, followed by the fourfold hollow site. The atop site is the least stable. This finding is rationalized by analyzing the ''local structures'' formed upon oxygen chemisorption. The binding energies and heats of adsorption at different oxygen coverages show that pairwise repulsive interactions are considerably stronger between oxygen atoms occupying fourfold sites than those occupying bridge sites. Analysis of the partial charge densities associated with Bloch states demonstrates that the O-Pt bond is considerably more localized at the bridge site. These effects cause a sharp drop in the heats of adsorption for oxygen on hollow sites when the coverage is increased from 0.25 to 0.5 ML. Mixing between oxygen p orbitals and Pt d orbitals can be observed over the whole metal d-band energy range.

ORBITAL MIXING IN CO CHEMISORPTION ON TRANSITION-METAL SURFACES

The chemisorption of CO on metal surfaces is widely accepted as a model for understanding chemical bonding between molecules and solid surfaces, but is nevertheless still a controversial subject. Ab initio total energy calculations using density functional theory with gradient corrections for CO chemisorption on an extended Pd{110} slab yield good agreement with experimental adsorption energies. Examination of the spatial distribution of individual Bloch states demonstrates that the conventional model for CO chemisorption involving charge donation from CO 5 sigma states to metal states and back-donation from metal states into CO 2 pi states is too simplistic, but the computational results provide direct insight into the chemical bonding within the framework of orbital mixing (or hybridisation). The results provide a sound basis for understanding the bonding between molecules and metal surfaces.

Cement As a Waste Form for Nuclear Fission Products: The Case of 90Sr and Its Daughters

One of the main challenges faced by the nuclear industry is the long-term confinement of nuclear waste. Because it is inexpensive and easy to manufacture, cement is the material of choice to store large volumes of radioactive materials, in particular the low-level medium-lived fission products. It is therefore of utmost importance to assess the chemical and structural stability of cement containing radioactive species. Here, we use ab initio calculations based on density functional theory (DFT) to study the effects of ⁹⁰Sr insertion and decay in C-S-H (calcium-silicate-hydrate) in order to test the ability of cement to trap and hold this radioactive fission product and to investigate the consequences of its β-decay on the cement paste structure. We show that ⁹⁰Sr is stable when it substitutes the Ca²⁺ ions in C-S-H, and so is its daughter nucleus ⁹⁰Y after β-decay. Interestingly, ⁹⁰Zr, daughter of ⁹⁰Y and final product in the decay sequence, is found to be unstable compared to the bulk phase of the element at zero K but stable when compared to the solvated ion in water. Therefore, cement appears as a suitable waste form for ⁹⁰Sr storage.

Towards a spectroscopically accurate set of potentials for heavy hydride laser cooling candidates: effective core potential calculations of BaH

BaH (and its isotopomers) is an attractive molecular candidate for laser cooling to ultracold temperatures and a potential precursor for the production of ultracold gases of hydrogen and deuterium. The theoretical challenge is to simulate the laser cooling cycle as reliably as possible and this paper addresses the generation of a highly accurate ab initio $^{2}\Sigma^+$ potential for such studies. The performance of various basis sets within the multi-reference configuration-interaction (MRCI) approximation with the Davidson correction (MRCI+Q)is tested and taken to the Complete Basis Set (CBS) limit. It is shown that the calculated molecular constants using a 46 electron Effective Core-Potential (ECP) and even-tempered augmented polarized core-valence basis sets (aug-pCV$n$Z-PP, n= 4 and 5) but only including three active electrons in the MRCI calculation are in excellent agreement with the available experimental values. The predicted dissociation energy De for the X$^2\Sigma^+$ state (extrapolated to the CBS limit) is 16895.12 cm$^{-1}$ (2.094 eV), which agrees within 0.1$\%$ of a revised experimental value of <16910.6 cm$^{-1}$, while the calculated re is within 0.03 pm of the experimental result.

Ultrafast charge-transfer in organic photovoltaic interfaces: geometrical and functionalization effects

Understanding the microscopic mechanisms of electronic excitation in organic photovoltaic cells is a challenging problem in the design of efficient devices capable of performing sunlight harvesting. Here we develop and apply an ab initio approach based on time-dependent density functional theory and Ehrenfest dynamics to investigate photoinduced charge transfer in small organic molecules. Our calculations include mixed quantum–classical dynamics with ions moving classically and electrons quantum mechanically, where no experimental external parameter other than the material geometry is required. We show that the behavior of photocarriers in zinc phthalocyanine (ZnPc) and C60 systems, an effective prototype system for organic solar cells, is sensitive to the atomic orientation of the donor and the acceptor units as well as the functionalization of covalent molecules at the interface. In particular, configurations with the ZnPc molecules facing on C60 facilitate charge transfer between substrate and molecules that occurs within 200 fs. In contrast, configurations where ZnPc is tilted above C60 present extremely low carrier injection efficiency even at longer times as an effect of the larger interfacial potential level offset and higher energetic barrier between the donor and acceptor molecules. An enhancement of charge injection into C60 at shorter times is observed as binding groups connect ZnPc and C60 in a dyad system. Our results demonstrate a promising way of designing and controlling photoinduced charge transfer on the atomic level in organic devices that would lead to efficient carrier separation and maximize device performance.

Epitaxially Grown Strained Pentacene Thin Film on Graphene Membrane

Organic-graphene system has emerged as a new platform for various applications such as flexible organic photovoltaics and organic light emitting diodes. Due to its important implication in charge transport, the study and reliable control of molecular packing structures at the graphene-molecule interface are of great importance for successful incorporation of graphene in related organic devices. Here, an ideal membrane of suspended graphene as a molecular assembly template is utilized to investigate thin-film epitaxial behaviors. Using transmission electron microscopy, two distinct molecular packing structures of pentacene on graphene are found. One observed packing structure is similar to the well-known bulk-phase, which adapts a face-on molecular orientation on graphene substrate. On the other hand, a rare polymorph of pentacene crystal, which shows significant strain along the c-axis, is identified. In particular, the strained film exhibits a specific molecular orientation and a strong azimuthal correlation with underlying graphene. Through ab initio electronic structure calculations, including van der Waals interactions, the unusual polymorph is attributed to the strong graphene-pentacene interaction. The observed strained organic film growth on graphene demonstrates the possibility to tune molecular packing via graphene-molecule interactions.

Multiscale Analysis for Field-Effect Penetration through Two-Dimensional Materials

Gate-tunable two-dimensional (2D) materials-based quantum capacitors (QCs) and van der Waals heterostructures involve tuning transport or optoelectronic characteristics by the field effect. Recent studies have attributed the observed gate-tunable characteristics to the change of the Fermi level in the first 2D layer adjacent to the dielectrics, whereas the penetration of the field effect through the one-molecule-thick material is often ignored or oversimplified. Here, we present a multiscale theoretical approach that combines first-principles electronic structure calculations and the Poisson–Boltzmann equation methods to model penetration of the field effect through graphene in a metal–oxide–graphene–semiconductor (MOGS) QC, including quantifying the degree of “transparency” for graphene two-dimensional electron gas (2DEG) to an electric displacement field. We find that the space charge density in the semiconductor layer can be modulated by gating in a nonlinear manner, forming an accumulation or inversion layer at the semiconductor/graphene interface. The degree of transparency is determined by the combined effect of graphene quantum capacitance and the semiconductor capacitance, which allows us to predict the ranking for a variety of monolayer 2D materials according to their transparency to an electric displacement field as follows: graphene > silicene > germanene > WS2 > WTe2 > WSe2 > MoS2 > phosphorene > MoSe2 > MoTe2, when the majority carrier is electron. Our findings reveal a general picture of operation modes and design rules for the 2D-materials-based QCs.

Structural and Electrical Investigation of C 60 –Graphene Vertical Heterostructures

Graphene, with its unique electronic and structural qualities, has become an important playground for studying adsorption and assembly of various materials including organic molecules. Moreover, organic/graphene vertical structures assembled by van der Waals interaction have potential for multifunctional device applications. Here, we investigate structural and electrical properties of vertical heterostructures composed of C60 thin film on graphene. The assembled film structure of C60 on graphene is investigated using transmission electron microscopy, which reveals a uniform morphology of C60 film on graphene with a grain size as large as 500 nm. The strong epitaxial relations between C60 crystal and graphene lattice directions are found, and van der Waals ab initio calculations support the observed phenomena. Moreover, using C60-graphene heterostructures, we fabricate vertical graphene transistors incorporating n-type organic semiconducting materials with an on/off ratio above 3 × 10(3). Our work demonstrates that graphene can serve as an excellent substrate for assembly of molecules, and attained organic/graphene heterostructures have great potential for electronics applications.