Quantifying Uncertainties in Solvation Procedures for Modeling Aqueous Phase Reaction Mechanisms.

Computational quantum chemistry provides fundamental chemical and physical insights into solvated reaction mechanisms across many areas of chemistry, especially in homogeneous and heterogeneous renewable energy catalysis. Such reactions may depend on explicit interactions with ions and solvent molecules that are nontrivial to characterize. Rigorously modeling explicit solvent effects with molecular dynamics usually brings steep computational costs while the performance of continuum solvent models such as polarizable continuum model (PCM), charge-asymmetric nonlocally determined local-electric (CANDLE), conductor-like screening model for real solvents (COSMO-RS), and effective screening medium method with the reference interaction site model (ESM-RISM) are less well understood for reaction mechanisms. Here, we revisit a fundamental aqueous hydride transfer reaction-carbon dioxide (CO2) reduction by sodium borohydride (NaBH4)-as a test case to evaluate how different solvent models perform in aqueous phase charge migrations that would be relevant to renewable energy catalysis mechanisms. For this system, quantum mechanics/molecular mechanics (QM/MM) molecular dynamics simulations almost exactly reproduced energy profiles from QM simulations, and the Na+ counterion in the QM/MM simulations plays an insignificant role over ensemble averaged trajectories that describe the reaction pathway. However, solvent models used on static calculations gave much more variability in data depending on whether the system was modeled using explicit solvent shells and/or the counterion. We pinpoint this variability due to unphysical descriptions of charge-separated states in the gas phase (i.e., self-interaction errors), and we show that using more accurate hybrid functionals and/or explicit solvent shells lessens these errors. This work closes with recommended procedures for treating solvation in future computational efforts in studying renewable energy catalysis mechanisms.

Li deposition and desolvation with electron transfer at a silicon/propylene-carbonate interface: transition-state and free-energy profiles by large-scale first-principles molecular dynamics.

We report the result of a large-scale first-principles molecular dynamics simulation under different electric biases performed to understand the charge transfer process coupling with lithium deposition and desolvation processes. We applied the effective screening medium (ESM) method to control the bias across the electrode/solution interface, and simulated a series of Li de-solvation and Li-deposition reactions occurring under the bias. Solvated Li+ in the bulk propylene carbonate migrates to the Si electrode surface and gradually de-solvates through the transition state. Introducing the blue-moon ensemble method, we determined the possible structures and activation energies for the transition states.

Toward full simulation of the electrochemical oxygen reduction reaction on Pt using first-principles and kinetic calculations.

Present fuel cells must increase the activity of the oxygen reduction reaction (ORR) on platinum (or Pt alloy) electrodes. Detailed simulation analyses can direct future investigations by providing a better understanding of the ORR. We adopted a density functional theory (DFT)-based, first-principles molecular dynamics simulation for the elementary steps of the electrochemical ORR on Pt(111). The two-step process involves successive protonation of O and OH, which are adsorbed on Pt. The relevant redox potentials were estimated by changing the coverage of OH(ad) and O(ad). The reaction energy profiles were determined along the reaction coordinate using the Blue-Moon ensemble method and a constant-bias scheme in the DFT calculations. These profiles at different biases were then used to generate activation energies and symmetry factors. Cyclic voltammetry (CV) and linear sweep voltammetry profiles were then calculated from the Butler-Volmer rate, Nernst equilibrium, and mass diffusion equations using these obtained parameters, literature values and appropriate prefactors in the rate equations. The experimentally observed reversible and irreversible peaks in CV were obtained. The irreversibility of the protonation of O(ad), , attributed to its higher activation energy, affects the ORR potential and thus fuel cell performance. It is therefore necessary not only to tune the adsorption energy of the O(ad) and OH(ad) intermediates, which are the origin of the "volcano plot", but also to tune 's activation energy to elevate the performance above that of the volcano-top.

Polyanthraquinone as a Reliable Organic Electrode for Stable and Fast Lithium Storage.

In spite of recent progress, there is still a lack of reliable organic electrodes for Li storage with high comprehensive performance, especially in terms of long-term cycling stability. Herein, we report an ideal polymer electrode based on anthraquinone, namely, polyanthraquinone (PAQ), or specifically, poly(1,4-anthraquinone) (P14AQ) and poly(1,5-anthraquinone) (P15AQ). As a lithium-storage cathode, P14AQ showed exceptional performance, including reversible capacity almost equal to the theoretical value (260 mA h g(-1); >257 mA h g(-1) for AQ), a very small voltage gap between the charge and discharge curves (2.18-2.14=0.04 V), stable cycling performance (99.4% capacity retention after 1000 cycles), and fast-discharge/charge ability (release of 69% of the low-rate capacity or 64% of the energy in just 2 min). Exploration of the structure-performance relationship between P14AQ and related materials also provided us with deeper understanding for the design of organic electrodes.

Effect of thermal motion on catalytic activity of nanoparticles in polar solvent.

In this study, we propose that electrode potential fluctuations due to the thermal motion of the solvent may serve to enhance the catalytic activity of nanostructures. The proposed model uses a simple, Marcus-type treatment of the statistical behavior of the solvent and the Butler-Volmer law for the instantaneous catalytic rate as a function of the electrode potential. The rapid development of probing techniques with high spatial and temporal resolution will help to further confirm and characterize the dynamical properties of nanostructures.

A method of orbital analysis for large-scale first-principles simulations.

An efficient method of calculating the natural bond orbitals (NBOs) based on a truncation of the entire density matrix of a whole system is presented for large-scale density functional theory calculations. The method recovers an orbital picture for O(N) electronic structure methods which directly evaluate the density matrix without using Kohn-Sham orbitals, thus enabling quantitative analysis of chemical reactions in large-scale systems in the language of localized Lewis-type chemical bonds. With the density matrix calculated by either an exact diagonalization or O(N) method, the computational cost is O(1) for the calculation of NBOs associated with a local region where a chemical reaction takes place. As an illustration of the method, we demonstrate how an electronic structure in a local region of interest can be analyzed by NBOs in a large-scale first-principles molecular dynamics simulation for a liquid electrolyte bulk model (propylene carbonate + LiBF4).

Charge manipulation in molecules encapsulated inside single-wall carbon nanotubes.

We report clear experimental evidence for the charge manipulation of molecules encapsulated inside single-wall carbon nanotubes (SWCNTs) using electrochemical doping techniques. We encapsulated ß-carotene (Car) inside SWCNTs and clarified electrochemical doping characteristics of their Raman spectra. C=C streching modes of encapsulated Car and a G band of SWCNTs showed clearly different doping behaviors as the electrochemical potentials were shifted. Electron extraction from encapsulated Car was clearly achieved. However, electrochemical characteristics of Car inside SWCNTs and doping mechanisms elucidated by calculations based on density-functional theory indicate the difficulty of charge manipulation of molecules inside SWCNTs due to the presence of strong on-site Coulomb repulsion energy at the molecules.

First-Principles Investigation of Charged Germagraphene as a Cathode Material for Dual-Carbon Batteries.

As part of the concerted effort for development of energy storage technologies, dual-ion batteries (DIBs) or dual-carbon batteries (DCBs) are attracting interest, owing primarily to their eco-friendly active materials. The use of carbon as the active materials of DCBs brings about several challenges involving capacity and stability. This contribution aims to provide an in-depth understanding of the structural and electronic properties of Ge-doped graphene (Germagraphene) as a novel cathode material for DCBs. Density functional theory (DFT) calculations are combined with the effective screening medium (ESM) method for analyzing the electronic and band structure of PF6 - anion-adsorbed Germagraphene under a potential bias. These theoretical investigations indicate that the use of Ge as a dopant for graphene has a positive impact on the adsorption of the anion on the cathode under both neutral and electrically biased conditions.

Electrocatalytic Mechanisms for an Oxygen Evolution Reaction at a Rhombohedral Boron Monosulfide Electrode/Alkaline Medium Interface.

Rhombohedral boron monosulfide (r-BS) with a layer stacking structure is a promising electrocatalyst for an oxygen evolution reaction (OER) within an alkaline solution. We investigated the catalytic mechanisms at the r-BS electrode/alkaline medium interface for an OER using hybrid solvation theory based on the first-principles method combined with classical solution theory. In this study, we elucidate the activities of the OER at the outermost r-BS sheet with and without various surface defects. The Gibbs free energies along the OER path indicate that the boron vacancies at the first and second layers of the r-BS surface (VB1 and VB2) can promote the OER. However, we found that the VB1 is easily occupied by the oxygen atom during the OER, degrading its electrocatalytic performance. In contrast, VB2 is suitable for the active site of the OER due to its structure stability. Next, we applied a bias voltage with the OER potential to the r-BS electrode. The bias voltage incorporates the positive excess surface charge into pristine r-BS and VB2, which can be understood by the relationship between the OER potential and potentials of zero charge at the r-BS electrode. Because the OH- ions are the starting point of the OER, the positively charged surface is kinetically favorable for the electrocatalyst owing to the attractive interaction with the OH- ions. Finally, we qualitatively discuss the flat-band potential at a semiconductor/alkaline solution interface. It suggests that p-type carrier doping could promote the catalytic performance of r-BS. These results explain the previous measurement of the OER performance with the r-BS-based electrode and provide valuable insights into developing a semiconductor electrode/water interface.

A Combined Reaction Path Search and Hybrid Solvation Method for the Systematic Exploration of Elementary Reactions at the Solid-Liquid Interface.

We present a combined simulation method of single-component artificial force induced reaction (SC-AFIR) and effective screening medium combined with the reference interaction site model (ESM-RISM), termed SC-AFIR+ESM-RISM. SC-AFIR automatically and systematically explores the chemical reaction pathway, and ESM-RISM directly simulates the precise electronic structure at the solid-liquid interface. Hence, SC-AFIR+ESM-RISM enables us to explore reliable reaction pathways at the solid-liquid interface. We applied it to explore the dissociation pathway of an H2O molecule at the Cu(111)/water interface. The reaction path networks of the whole reaction and the minimum energy paths from H2O to H2 + O depend on the interfacial environment. The qualitative difference in the energy diagrams and the resulting change in the kinematically favored dissociation pathway upon changing the solvation environments are discussed. We believe that SC-AFIR+ESM-RISM will be a powerful tool to reveal the details of chemical reactions in surface catalysis and electrochemistry.

First-principles molecular dynamics at a constant electrode potential.

A simulation scheme for performing first-principles molecular dynamics at a constant electrode potential is presented, opening the way for a more realistic modeling of voltage-driven devices. The system is allowed to exchange electrons with a reservoir at fixed potential, and dynamical equations for the total electronic charge are derived by using the potential energy of the extended system. In combination with a thermostat, this potentiostat scheme reproduces thermal fluctuations of the charge with the correct statistics, implying a realistic treatment of the potential as a control variable. Practically, the dynamics of the charge are decoupled from the electronic structure calculations, making the scheme easily implementable in existing first-principles molecular dynamics codes. Our approach is demonstrated on a test system by considering various test cases.

Large-scale first-principles molecular dynamics for electrochemical systems with O(N) methods.

A method for large-scale first-principles molecular dynamics (MD) simulations on electrochemical systems has been developed by combining the effective screening medium (ESM) method with O(N) density functional theory (DFT). This implementation has been significantly simplified by the introduction of neutral atom potentials, which minimizes the modifications to existing DFT code. In order to demonstrate ability of this implementation, it has been applied to an electrochemical system consisting of a H-Si(111) electrode, which is a candidate anode for high-capacity Li-ion secondary batteries, and a propylene carbonate (PC) solvent to simulate how PC molecules in the vicinity of the electrode surface respond to an imposed electric field. The large-scale MD simulation clearly demonstrates that the combination of the ESM and O(N) DFT methods provides a useful tool for first-principles investigation of complicated electrochemical systems such as high-capacity batteries.

Theoretical Analysis on the Stability of 1-Pyrenebutanoic Acid Succinimidyl Ester Adsorbed on Graphene.

The adsorbed structure of 1-pyrenebutanoic acid succinimidyl ester (PASE) on graphene was investigated based on density functional theory. We found two locally stable structures: a straight structure with the chainlike part of butanoic acid succinimidyl ester (BSE) lying down and a bent structure with the BSE part directed away from graphene, keeping the pyrene (Py) part adsorbed on graphene. Then, to elucidate the adsorption mechanism, we separately estimated the contributions of the Py and BSE parts to the entire PASE adsorption, and the adsorption effect of the BSE part was found to be secondary in comparison to the contribution of the Py. Next, the mobility of the BSE part at room temperature was confirmed by the activation energy barrier between straight and bent structures. To take account of the external environment, we considered the presence of amino acids and the hydration effect by a three-dimensional reference interaction site model. The contributions of glycine molecules and the solvent environment to stabilizing the bent PASE structure relative to the straight PASE structure were found. Therefore, the effect of the external environment around PASE is of importance when the standing-up process of the BSE part from graphene is considered.

Semiconducting electronic property of graphene adsorbed on (0001) surfaces of SiO2.

First-principles total energy calculations are performed to investigate the energetics and electronic structures of graphene adsorbed on both an oxygen-terminated SiO2 (0001) surface and a fully hydroxylated SiO2 (0001) surface. We find that there are several stable adsorption sites for graphene on both O-terminated and hydroxylated SiO2 surfaces. The binding energy in the most stable geometry is found to be 15 meV per C atom, indicating a weak interaction between graphene and SiO2 (0001) surfaces. We also find that the graphene adsorbed on SiO2 is a semiconductor irrespective of the adsorption arrangement due to the variation of on-site energy induced by the SiO2 substrate.

Reversible redox reaction and water configuration on a positively charged platinum surface: first principles molecular dynamics simulation.

The water dissociation reaction and water molecule configuration on a positively charged platinum (111) surface were investigated by means of first principles molecular dynamics under periodic boundary conditions. Water molecules on the Pt surface were mostly in the O-down orientation but some H-down structures were also found. OH(-) ion, generated by removing H from H(2)O in the bulk region, moved to the Pt surface, on which a positive charge is induced, by a Grotthuss-like proton-relay mechanism and adsorbed on it as OH(Pt). Hydrogen atom exchange between OH(Pt) and a near-by water molecule frequently occurred on the Pt surface and had a low activation energy of the same order as room temperature energy. When a positive charge (7 µC cm(-2)) was added to the Pt surface, H(3)O(+) and OH(Pt) were generated from 2H(2)O on the Pt. This may be coupled with an electron transfer to the Pt electrode [2H(2)O → H(3)O(+) + OH(Pt) + e(-)]. The opposite reaction was also observed on the same charged surface during a simulation of duration about 10 ps; it is a reversible redox reaction. When further positive charge (14 µC cm(-2)) was added, the reaction shifted to the right hand side completely. Thus, this one-electron transfer reaction, which is a part of the oxygen electrode reaction in fuel cells and water electrolysis, was confirmed to be a low activation energy process.

Influence of disorder on conductance in bilayer graphene under perpendicular electric field.

Electron transport in bilayer graphene placed under a perpendicular electric field is revealed experimentally. Steep increase of the resistance is observed under high electric field; however, the resistance does not diverge even at low temperatures. The observed temperature dependence of the conductance consists of two contributions: the thermally activated (TA) conduction and the variable range hopping (VRH) conduction. We find that for the measured electric field range (0-1.3 V/nm) the mobility gap extracted from the TA behavior agrees well with the theoretical prediction for the band gap opening in bilayer graphene, although the VRH conduction deteriorates the insulating state more seriously in bilayer graphene with smaller mobility. These results show that the improvement of the mobility is crucial for the successful operation of the bilayer graphene field effect transistor.

Toward Engineering of Solution Microenvironments for the CO2 Reduction Reaction: Unraveling pH and Voltage Effects from a Combined Density-Functional-Continuum Theory.

Engineering the electrolyte microenvironment represents an attractive route to tuning the selectivity of electrocatalytic reactions beyond catalyst composition and morphology. However, harnessing the full potential of this approach requires understanding the interplay between voltage, electrolyte composition, and adsorbate binding within the electrical double layer, which is absent from the usual theoretical approaches. In this work, we apply a recently developed density functional theory (DFT)-continuum approach based on the effective screening medium method and reference interaction site model (ESM-RISM) to explore electrolyte effects with an enhanced description of the electrochemical interface. Applying this method to the binding of CO adsorbates in potassium-containing electrolytes on copper, a problem of direct relevance to CO2 electroreduction to value-added products, we show that the interdependence of voltage and pH leads to an unexpected change in adsorption site preference on Cu(001) terraces. Our findings highlight the often-overlooked importance of the electrical double-layer structure for predicting catalyst operation.

Two-Phase Reaction Mechanism for Fluorination and Defluorination in Fluoride-Shuttle Batteries: A First-Principles Study.

Fluoride-shuttle batteries (FSBs), which are based on fluoride-ion transfer, have attracted attention because of their high theoretical energy densities. The fluorination and defluorination reactions at the electrodes are the possible rate-determining steps in FSBs, and understanding the mechanism is important to achieve smooth charge/discharge. In this study, we discuss the thermodynamically favored pathways for the fluorination and defluorination reactions and compare the reactions through the solid-solution and two-phase-coexistent states by density functional theory (DFT) calculations. The free energies of the solid-solution and two-phase states approximate the energies calculated by DFT, and their accuracy was validated by comparison with experimental formation enthalpies and free energies. The relative formation enthalpies of typical, transition, and relativistic metal (Tl, Pb, and Bi) fluorides are well reproduced by DFT calculations within 0.1, 0.2, and 0.4 eV, respectively. We also show that the reaction pathway can be determined by comparing the formation enthalpies of the metal fluoride H, a fluorine vacancy HV, and an interstitial fluorine defect HI from the simple selection rule. The enthalpy relation of HI > H > -HV observed in all the calculations strongly suggests that fluorination and defluorination in FSB electrodes occur by a two-phase reaction. This fluorination and defluorination mechanism will be useful to clarify the rate-determining step in FSBs.

Negative dielectric constant of water confined in nanosheets.

Electric double-layer capacitors are efficient energy storage devices that have the potential to account for uneven power demand in sustainable energy systems. Earlier attempts to improve an unsatisfactory capacitance of electric double-layer capacitors have focused on meso- or nanostructuring to increase the accessible surface area and minimize the distance between the adsorbed ions and the electrode. However, the dielectric constant of the electrolyte solvent embedded between adsorbed ions and the electrode surface, which also governs the capacitance, has not been previously exploited to manipulate the capacitance. Here we show that the capacitance of electric double-layer capacitor electrodes can be enlarged when the water molecules are strongly confined into the two-dimensional slits of titanium carbide MXene nanosheets. Using electrochemical methods and theoretical modeling, we find that dipolar polarization of strongly confined water resonantly overscreens an external electric field and enhances capacitance with a characteristically negative dielectric constant of a water molecule.

Specific ion effects at graphitic interfaces.

Improved understanding of aqueous solutions at graphitic interfaces is critical for energy storage and water desalination. However, many mechanistic details remain unclear, including how interfacial structure and response are dictated by intrinsic properties of solvated ions under applied voltage. In this work, we combine hybrid first-principles/continuum simulations with electrochemical measurements to investigate adsorption of several alkali-metal cations at the interface with graphene and within graphene slit-pores. We confirm that adsorption energy increases with ionic radius, while being highly dependent on the pore size. In addition, in contrast with conventional electrochemical models, we find that interfacial charge transfer contributes non-negligibly to this interaction and can be further enhanced by confinement. We conclude that the measured interfacial capacitance trends result from a complex interplay between voltage, confinement, and specific ion effects-including ion hydration and charge transfer.