Electron transfer of functionalized quinones in acetonitrile.

Quinones are redox active organic molecules that have been proposed as an alternative choice to metal-based materials in electrochemical energy storage devices. Functionalization allows one to fine tune not only their chemical stability but also the redox potential and kinetics of the electron transfer reaction. However, the reaction rate constant is not only determined by the redox species but also impacted by solvent effects. In this work, we show how the functionalization of benzoquinone with different functional groups impacts the solvent reorganization free energies of electron transfer half-reactions in acetonitrile. The use of molecular density functional theory, whose computational cost for studying the electron transfer reaction is considerably reduced compared to the state-of-the-art molecular dynamics simulations, enables us to perform a systematic study. We validate the method by comparing the predictions of the solvation shell structure and the free energy profiles for electron transfer reaction to the reference classical molecular dynamics simulations in the case of anthraquinone solvated in acetonitrile. We show that all the studied electron transfer half-reactions follow the Marcus theory, regardless of functional groups. Consequently, the solvent reorganization free energy decreases as the molecular size increases.

Multi-scale simulation of the adsorption of lithium ion on graphite surface: From quantum Monte Carlo to molecular density functional theory.

The structure of the double-layer formed at the surface of carbon electrodes is governed by the interactions between the electrode and the electrolyte species. However, carbon is notoriously difficult to simulate accurately, even with well-established methods such as electronic density functional theory and molecular dynamics. Here, we focus on the important case of a lithium ion in contact with the surface of graphite, and we perform a series of reference quantum Monte Carlo calculations that allow us to benchmark various electronic density functional theory functionals. We then fit an accurate carbon-lithium pair potential, which is used in molecular density functional theory calculations to determine the free energy of the adsorption of the ion on the surface in the presence of water. The adsorption profile in aqueous solution differs markedly from the gas phase results, which emphasize the role of the solvent on the properties of the double-layer.

Solvation of anthraquinone and TEMPO redox-active species in acetonitrile using a polarizable force field.

Redox-active molecules are of interest in many fields, such as medicine, catalysis, or energy storage. In particular, in supercapacitor applications, they can be grafted to ionic liquids to form so-called biredox ionic liquids. To completely understand the structural and transport properties of such systems, an insight at the molecular scale is often required, but few force fields are developed ad hoc for these molecules. Moreover, they do not include polarization effects, which can lead to inaccurate solvation and dynamical properties. In this work, we developed polarizable force fields for redox-active species anthraquinone (AQ) and 2,2,6,6-tetra-methylpiperidinyl-1-oxyl (TEMPO) in their oxidized and reduced states as well as for acetonitrile. We validate the structural properties of AQ, AQ•-, AQ2-, TEMPO•, and TEMPO+ in acetonitrile against density functional theory-based molecular dynamics simulations and we study the solvation of these redox molecules in acetonitrile. This work is a first step toward the characterization of the role played by AQ and TEMPO in electrochemical and catalytic devices.

A semiclassical Thomas-Fermi model to tune the metallicity of electrodes in molecular simulations.

Spurred by the increasing needs in electrochemical energy storage devices, the electrode/electrolyte interface has received a lot of interest in recent years. Molecular dynamics simulations play a prominent role in this field since they provide a microscopic picture of the mechanisms involved. The current state-of-the-art consists of treating the electrode as a perfect conductor, precluding the possibility to analyze the effect of its metallicity on the interfacial properties. Here, we show that the Thomas-Fermi model provides a very convenient framework to account for the screening of the electric field at the interface and differentiating good metals such as gold from imperfect conductors such as graphite. All the interfacial properties are modified by screening within the metal: the capacitance decreases significantly and both the structure and dynamics of the adsorbed electrolyte are affected. The proposed model opens the door for quantitative predictions of the capacitive properties of materials for energy storage.

A first-principles investigation of the structural and electrochemical properties of biredox ionic species in acetonitrile.

Biredox ionic liquids are a new class of functionalized electrolytes that may play an important role in future capacitive energy storage devices. By allowing additional storage of electrons inside the liquids, they can improve device performance significantly. However current devices employ nanoporous carbons in which the diffusion of the liquid and the adsorption of the ions could be affected by the occurrence of electron-transfer reactions. It is therefore necessary to understand better the thermodynamics and the kinetics of such reactions in biredox ionic liquids. Here we perform ab initio molecular dynamics simulations of both the oxidized and reduced species of several redox-active ionic molecules (used in biredox ionic liquids) dissolved in acetonitrile solvent and compare them with the bare redox molecules. We show that in all the cases, it is necessary to introduce a two Gaussian state model to calculate the reaction free energies accurately. These reaction free energies are only slightly affected by the presence of the IL group on the molecule. We characterize the structure of the solvation shell around the redox active part of the molecules and show that in the case of TEMPO-based molecules strong reorientation effects occur during the oxidation reaction.

A first-principles computational comparison of defect-free and disordered, fluorinated anatase TiO2 (001) interfaces with water.

Chemical doping and other surface modifications have been used to engineer the bulk properties of materials, but their influence on the surface structure and consequently the surface chemistry are often unknown. Previous work has been successful in fluorinating anatase TiO2 with charge balance achieved via the introduction of Ti vacancies rather than the reduction of Ti. Our work here investigates the interface between this fluorinated titanate with cationic vacancies and a monolayer of water via density functional theory based molecular dynamics. We compute the projected density of states for only those atoms at the interface and for those states that fall within 1 eV of the Fermi level for various steps throughout the simulation, and we determine that the variation in this visualization of the density of states serves as a reasonable tool to anticipate where surfaces are most likely to be reactive. In particular, we conclude that water dissociation at the surface is the main mechanism that influences the anatase (001) surface whereas the change in the density of states at the surface of the fluorinated structure is influenced primarily through the adsorption of water molecules.

Electronic Excitation Dynamics in Liquid Water under Proton Irradiation.

Molecular behaviour of liquid water under proton irradiation is of great importance to a number of technological and medical applications. The highly energetic proton generates a time-varying field that is highly localized and heterogeneous at the molecular scale, and massive electronic excitations are produced as a result of the field-matter interaction. Using first-principles quantum dynamics simulations, we reveal details of how electrons are dynamically excited through non-equilibrium energy transfer from highly energetic protons in liquid water on the atto/femto-second time scale. Water molecules along the path of the energetic proton undergo ionization at individual molecular level, and the excitation primarily derives from lone pair electrons on the oxygen atom of water molecules. A reduced charge state on the energetic proton in the condensed phase of water results in the strongly suppressed electronic response when compared to water molecules in the gas phase. These molecular-level findings provide important insights into understanding the water radiolysis process under proton irradiation.

Diffusion quantum Monte Carlo study of martensitic phase transition energetics: The case of phosphorene.

Recent technical advances in dealing with finite-size errors make quantum Monte Carlo methods quite appealing for treating extended systems in electronic structure calculations, especially when commonly used density functional theory (DFT) methods might not be satisfactory. We present a theoretical study of martensitic phase transition energetics of a two-dimensional phosphorene by employing diffusion Monte Carlo (DMC) approach. The DMC calculation supports DFT prediction of having a rather diffusive barrier that is characterized by having two transition states, in addition to confirming that the so-called black and blue phases of phosphorene are essentially degenerate. At the same time, the DFT calculations do not provide the quantitative accuracy in describing the energy changes for the martensitic phase transition even when hybrid exchange-correlation functional is employed. We also discuss how mechanical strain influences the stabilities of the two phases of phosphorene.

Role of Surface Termination on Hot Electron Relaxation in Silicon Quantum Dots: A First-Principles Dynamics Simulation Study.

The role of surface termination on phonon-mediated relaxation of an excited electron in quantum dots was investigated using first-principles simulations. The surface terminations of a silicon quantum dot with hydrogen and fluorine atoms lead to distinctively different relaxation behaviors, and the fluorine termination shows a nontrivial relaxation process. The quantum confined electronic states are significantly affected by the surface of the quantum dot, and we find that a particular electronic state dictates the relaxation behavior through its infrequent coupling to neighboring electronic states. Dynamical fluctuation of this electronic state results in a slow shuttling behavior within the manifold of unoccupied electronic states, controlling the overall dynamics of the excited electron with its characteristic frequency of this shuttling behavior. The present work revealed a unique role of surface termination, dictating the hot electron relaxation process in quantum-confined systems in the way that has not been considered previously.

Theoretical oxidation state analysis of Ru-(bpy)3: influence of water solvation and Hubbard correction in first-principles calculations.

Oxidation state is a powerful concept that is widely used in chemistry and materials physics, although the concept itself is arguably ill-defined quantum mechanically. In this work, we present impartial comparison of four, well-recognized theoretical approaches based on Lowdin atomic orbital projection, Bader decomposition, maximally localized Wannier function, and occupation matrix diagonalization, for assessing how well transition metal oxidation states can be characterized. Here, we study a representative molecular complex, tris(bipyridine)ruthenium. We also consider the influence of water solvation through first-principles molecular dynamics as well as the improved electronic structure description for strongly correlated d-electrons by including Hubbard correction in density functional theory calculations.