Revised self-consistent continuum solvation in electronic-structure calculations

Toward Realistic Models of the Electrocatalytic Oxygen Evolution Reaction

The electrocatalytic oxygen evolution reaction (OER) supplies the protons and electrons needed to transform renewable electricity into chemicals and fuels. However, the OER is kinetically sluggish; it operates at significant rates only when the applied potential far exceeds the reversible voltage. The origin of this overpotential is hidden in a complex mechanism involving multiple electron transfers and chemical bond making/breaking steps. Our desire to improve catalytic performance has then made mechanistic studies of the OER an area of major scientific inquiry, though the complexity of the reaction has made understanding difficult. While historically, mechanistic studies have relied solely on experiment and phenomenological models, over the past twenty years ab initio simulation has been playing an increasingly important role in developing our understanding of the electrocatalytic OER and its reaction mechanisms. In this Review we cover advances in our mechanistic understanding of the OER, organized by increasing complexity in the way through which the OER is modeled. We begin with phenomenological models built using experimental data before reviewing early efforts to incorporate ab initio methods into mechanistic studies. We go on to cover how the assumptions in these early ab initio simulations─no electric field, electrolyte, or explicit kinetics─have been relaxed. Through comparison with experimental literature, we explore the veracity of these different assumptions. We summarize by discussing the most critical open challenges in developing models to understand the mechanisms of the OER.

Hybrid Cluster-Continuum Method for Single-Ion Solvation Free Energy in Acetonitrile Solvent

A new hybrid discrete-continuum approach named the cluster-continuum static approximation (CCSA) has been proposed for acetonitrile solvent. The continuum part uses the conductor-like polarizable continuum model for electrostatic and a surface area-dependent term for nonelectrostatic solvation. The CCSA includes only one explicit acetonitrile solvent molecule and a damping function, which makes the CCSA method reduce to pure continuum solvation in the case of weaker potential of mean force for solute-solvent interaction. The performance of the model was tested for 22 anions and 22 cations, including challenge species that cannot be adequately described by pure continuum solvation. A comparison was done with the widely used solvent model density (SMD) model. For anions, the CCSA reduces to pure continuum solvation and the method has the same performance as the SMD model, with a standard deviation of the mean signed error (SD-MSE) of 2.7 kcal mol-1 for both models. However, the CCSA method for cations considerably outperforms the SMD model, with an SD-MSE of 3.3 kcal mol-1 for the former and 8.4 kcal mol-1 for the latter. The method can be automated, and the present study suggests that continuum solvation models could be parameterized taking into account the explicit solvation as proposed in this work.

Correcting implicit solvation at metal/water interfaces through the incorporation of competitive water adsorption

Conventional continuum solvation models are ubiquitous in computational catalysis, including for describing metal/water interfaces, which are relevant to both solution-phase heterogeneous catalysis and electrocatalysis. Nonetheless, we find that such continuum models qualitatively fail to describe both the adsorption free energy and conformational preference for many organic molecules at such interfaces, largely due to the failure of continuum models to incorporate the role of competitive water adsorption. We develop a simple phenomenological model that accounts for competitive water adsorption and show that the model, when used in conjunction with continuum solvation, provides a dramatic improvement in the description of both adsorption and conformational preference. The model is also extended to additionally incorporate the influence of applied potential at the electrode surface, thus facilitating computationally efficient applications to scenarios including electrocatalysis.

Emerging Atomistic Modeling Methods for Heterogeneous Electrocatalysis

Heterogeneous electrocatalysis lies at the center of various technologies that could help enable a sustainable future. However, its complexity makes it challenging to accurately and efficiently model at an atomic level. Here, we review emerging atomistic methods to simulate the electrocatalytic interface with special attention devoted to the components/effects that have been challenging to model, such as solvation, electrolyte ions, electrode potential, reaction kinetics, and pH. Additionally, we review relevant computational spectroscopy methods. Then, we showcase several examples of applying these methods to understand and design catalysts relevant to green hydrogen. We also offer experimental views on how to bridge the gap between theory and experiments. Finally, we provide some perspectives on opportunities to advance the field.

PW-SMD: A Plane-Wave Implicit Solvation Model Based on Electron Density for Surface Chemistry and Crystalline Systems in Aqueous Solution

Electron density-based implicit solvation models are a class of techniques for quantifying solvation effects and calculating free energies of solvation without an explicit representation of solvent molecules. Integral to the accuracy of solvation modeling is the proper definition of the solvation shell separating the solute molecule from the solvent environment, allowing for a physical partitioning of the free energies of solvation. Unlike state-of-the-art implicit solvation models for molecular quantum chemistry calculations, e.g., the solvation model based on solute electron density (SMD), solvation models for systems under periodic boundary conditions with plane-wave (PW) basis sets have been limited in their accuracy. Furthermore, a unified implicit solvation model with both homogeneous solution-phase and heterogeneous interfacial structures treated on equal footing is needed. In order to address this challenge, we developed a high-accuracy solvation model for periodic PW calculations that is applicable to molecular, ionic, interfacial, and bulk-phase chemistry. Our model, PW-SMD, is an extension of the SMD molecular solvation model to periodic systems in water. The free energy of solvation is partitioned into the electrostatic and cavity-dispersion-solvent structure (CDS) contributions. The electrostatic contributions of the solvation shell surrounding solute structures are parametrized based on their geometric and physical properties. In addition, the nonelectrostatic contribution to the solvation energy is accounted for by extending the CDS formalism of SMD to incorporate periodic boundary conditions. We validate the accuracy and robustness of our solvation model by comparing predicted solvation free energies against experimental data for molecular and ionic systems, carved-cluster composite energetic models of solvated reaction energies and barriers on surface systems, and deep-learning-accelerated ab initio molecular dynamics (AIMD). Our developed periodic implicit solvation model shows significantly improved accuracy compared to previous work (namely, solvation models in aqueous solution) and can be applied to simulate solvent effects in a wide range of surface and crystalline materials.

In-silico mechanistic analysis of adsorption of Iodinated Contrast Media agents on graphene surface

The study aims at assessing the potential of graphene-based adsorbents to reduce environmental impacts of Iodinated Contrast Media Agents (ICMs). We analyze an extensive collection of ICMs. A modeling approach resting on molecular docking and Density Functional Theory simulations is employed to examine the adsorption process at the molecular level. The study also relies on a Quantitative Structure-Activity Relationship (QSAR) modeling framework to correlate molecular properties with the adsorption energy (Ead) of ICMs, thus enabling identification of the key mechanisms underpinning adsorption and of the key factors contributing to it. A collection of distinct QSAR-based models is developed upon relying on Multiple Linear Regression and a standard genetic algorithm method. Having at our disposal multiple models enables us to take into account the uncertainty associated with model formulation. Maximum Likelihood and formal model identification/discrimination criteria (such as Bayesian and/or information theoretic criteria) are then employed to complement the traditional QSAR modeling phase. This has the advantage of (a) providing a rigorous ranking of the alternative models included in the selected set and (b) quantifying the relative degree of likelihood of each of these models through a weight or posterior probability. The resulting workflow of analysis enables one to seamlessly embed DFT and QSAR studies within a theoretical framework of analysis that explicitly takes into account model and parameter uncertainty. Our results suggest that graphene-based surfaces constitute a promising adsorbent for ICMs removal, π-π stacking being the primary mechanism behind ICM adsorption. Furthermore, our findings offer valuable insights into the potential of graphene-based adsorbent materials for effectively removing ICMs from water systems. They contribute to ascertain the significance of various factors (such as, e.g., the distribution of atomic van der Waals volumes, overall molecular complexity, the presence and arrangement of Iodine atoms, and the presence of polar functional groups) on the adsorption process.

Coupling molecular density functional theory with converged selected configuration interaction methods to study excited states in aqueous solution

This paper presents the first implementation of a coupling between advanced wavefunction theories and molecular density functional theory (MDFT). This method enables the modeling of solvent effect into quantum mechanical (QM) calculations by incorporating an electrostatic potential generated by solvent charges into the electronic Hamiltonian. Solvent charges are deduced from the spatially and angularly dependent solvent particle density. Such a density is obtained through the minimization of the functional associated with the molecular mechanics (MM) Hamiltonian describing the interaction between the fluid particles. The introduced QM/MDFT framework belongs to QM/MM family of methods, but its originality lies in the use of MDFT as the MM solver, offering two main advantages. First, its functional formulation makes it competitive with respect to sampling-based molecular mechanics. Second, it preserves a molecular-level description lost in macroscopic continuum approaches. The excited state properties of water and formaldehyde molecules solvated into water have been computed at the selected configuration interaction (SCI) level. The excitation energies and dipole moments have been compared with experimental data and previous theoretical work. A key finding is that using the Hartree-Fock method to describe the solute allows for predicting the solvent charge around the ground state with sufficient precision for the subsequent SCI calculations of excited states. This significantly reduces the computational cost of the described procedure, paving the way for the study of more complex molecules.

Searching for the Thinnest Metallic Wire

One-dimensional materials have gained much attention in the last decades: from carbon nanotubes to ultrathin nanowires to few-atom atomic chains, these can all display unique electronic properties and great potential for next-generation applications. Exfoliable bulk materials could naturally provide a source for one-dimensional wires with a well-defined structure and electronics. Here, we explore a database of one-dimensional materials that could be exfoliated from experimentally known three-dimensional van der Waals compounds, searching for metallic wires that are resilient to Peierls distortions and could act as vias or interconnects for future downscaled electronic devices. As the one-dimensional nature makes these wires particularly susceptible to dynamical instabilities, we carefully characterize vibrational properties to identify stable phases and characterize electronic and dynamical properties. Our search discovers several stable wires; notably, we identify what could be the thinnest possible exfoliable metallic wire, CuC2, coming a step closer to the ultimate limit in material downscaling.

Solvation Free Energies from Machine Learning Molecular Dynamics

The present work proposes an extension to the approach of [Xi, C; et al. J. Chem. Theory Comput. 2022, 18, 6878] to calculate ion solvation free energies from first-principles (FP) molecular dynamics (MD) simulations of a hybrid solvation model. The approach is first re-expressed within the quasi-chemical theory of solvation. Then, to allow for longer simulation times than the original first-principles molecular dynamics approach and thus improve the convergence of statistical averages at a fraction of the original computational cost, a machine-learned (ML) energy function is trained on FP energies and forces and used in the MD simulations. The ML workflow and MD simulation times (≈200 ps) are adjusted to converge the predicted solvation energies within a chemical accuracy of 0.04 eV. The extension is successfully benchmarked on the same set of alkaline and alkaline-earth ions.

Constant potential energetics of metallic and semiconducting electrodes: A benchmark study on 2D materials

Grand-canonical (GC) constant-potential methods within an implicit solvent environment provide a general approach to compute the potential-dependent energetics at electrified solid-liquid interfaces with first-principles density-functional theory. Here, we use a mindfully chosen set of 27 isostructural 2D metal halides MX2 to analyze the variation of this energetics when the electronic structure changes from metallic to semiconducting and insulating state. Apart from expectable changes due to the opening up of the electronic bandgap, the calculations also show an increasing sensitivity to the numerical Brillouin zone integration and electronic smearing, which imposes computational burdens in practice. We rationalize these findings within the picture of the total interfacial capacitance arising from a series connection of the electrochemical double-layer capacitance and the so-called quantum capacitance resulting from the filling of electronic states inside the electrode. For metals, the electrochemical double-layer capacitance dominates at all potentials, and the entire potential drop takes place in the electrolyte. For semiconductors, the potential drop occurs instead fully or partially inside the electrode at potentials within or just outside the bandgap. For 2D semiconductors, the increased sensitivity to numerical parameters then results from the concomitantly increased contribution of the quantum capacitance that is harder to converge. Fortunately, this understanding motivates a simple extension of the CHE + DL approximation for metals, which provides the approximate GC energetics of 2D semiconductors using only quantities that can be obtained from computationally undemanding calculations at the point of zero charge and a generic double-layer capacitance.

Elucidating Solvatochromic Shifts in Two-Dimensional Photocatalysts by Solving the Bethe-Salpeter Equation Coupled with Implicit Solvation Method

Many studies have focused on tailoring the photophysical properties of two-dimensional (2D) materials for photocatalytic (PC) or photoelectrochemical (PEC) applications. To understand the optical properties of 2D materials in solution, we established a computational method that combined the Bethe-Salpeter equation (BSE) calculations with our GW-GPE method, allowing for GW/BSE-level calculations with implicit solvation described using the generalized Poisson equation (GPE). We applied this method to MoS2, phosphorene (PP), and g-C3N4 and found that when the solvent dielectric increased, it reduced the exciton binding energy and quasiparticle bandgap, resulting in almost no solvatochromic shift in the excitonic peaks of MoS2 and PP, which is consistent with previous experiments. However, our calculations predicted that the solvent dielectric had a significant impact on the excitonic properties of g-C3N4, exhibiting a large solvatochromic shift. We expect that our GW/BSE-GPE method will offer insights into the design of 2D materials for PC and PEC applications.

Soft-sphere continuum solvation models for nonaqueous solvents

Solvation effects profoundly influence the characteristics and behavior of chemical systems in liquid solutions. The interaction between solute and solvent molecules intricately impacts solubility, reactivity, stability, and various chemical processes. Continuum solvation models gained prominence in quantum chemistry by implicitly capturing these interactions and enabling efficient investigations of diverse chemical systems in solution. In comparison, continuum solvation models in condensed matter simulation are very recent. Among these, the self-consistent continuum solvation (SCCS) and the soft-sphere continuum solvation models (SSCS) have been among the first to be successfully parameterized and extended to model periodic systems in aqueous solutions and electrolytes. As most continuum approaches, these models depend on a number of parameters that are linked to experimental or theoretical properties of the solvent, or that can be tuned based on reference data. Here, we present a systematic parameterization of the SSCS model for over 100 nonaqueous solvents. We validate the model's efficacy across diverse solvent environments by leveraging experimental solvation-free energies and partition coefficients from comprehensive databases. The average root means square error over all the solvents was calculated as 0.85 kcal/mol which is below the chemical accuracy (1 kcal/mol). Similarly to what has been reported by Hille et al. (J. Chem. Phys. 2019, 150, 041710.) for the SCCS model, a single-parameter model accurately reproduces experimental solvation energies, showcasing the transferability and predictive power of these continuum approaches. Our findings underscore the potential for a unified approach to predict solvation properties, paving the way for enhanced computational studies across various chemical environments.

Conformational preference of dipeptide zwitterions in aqueous solvents

Proper description of solvent effects is challenging for theoretical methods, particularly if the solute is a zwitterion. Here, a series of theoretical procedures are used to determine the preferred solvated conformations of twelve hydrophobic dipeptides (Leu-Leu, Leu-Phe, Phe-Leu, Ile-Leu, Phe-Phe, Ala-Val, Val-Ala, Ala-Ile, Ile-Ala, Ile-Val, Val-Ile and Val-Val) in the zwitterionic state. First, the accuracy of density functional theory (DFT), combined with different implicit solvent models, for describing zwitterions in aqueous solvent is assessed by comparing the predicted against the experimental glycine tautomerization energy, i.e., the energetic difference between canonical and zwitterionic glycine in aqueous solvents. It is found that among the tested solvation schemes, the charge-asymmetric nonlocally determined local-electric solvation model (CANDLE) predicts an energetic difference in excellent agreement with the experimental value. Next, DFT-CANDLE is used to determine the most favorable solvated conformation for each of the investigated dipeptide zwitterions. The CANDLE-solvated structures are obtained by exploring the conformational space of each dipeptide zwitterion concatenating DFT calculations, in vacuum, with classical molecular dynamics simulations, in explicit solvents, and DFT calculations including explicit water molecules. It is found that the energetically most favorable conformations are similar to those of the dipeptide zwitterions in their respective crystal structures. Such structural agreement is indicative of the DFT-CANDLE accomplishment of the description of solvated zwitterions, and suggests that these biomolecules self-assemble as quasi-rigid objects.

Incorporating Electrolyte Correlation Effects into Variational Models of Electrochemical Interfaces

We propose a way for obtaining a classical free energy density functional for electrolytes based on a first-principle many-body partition function. Via a one-loop expansion, we include coulombic correlations beyond the conventional mean-field approximation. To examine electrochemical interfaces, we integrate the electrolyte free energy functional into a hybrid quantum-classical model. This scheme self-consistently couples electronic, ionic, and solvent degrees of freedom and incorporates electrolyte correlation effects. The derived free energy functional causes a correlation-induced enhancement in interfacial counterion density and leads to an overall increase in capacitance. This effect is partially compensated by a reduction of the dielectric permittivity of interfacial water. At larger surface charge densities, ion crowding at the interface stifles these correlation effects. While scientifically intriguing already at planar interfaces, we anticipate these correlation effects to play an essential role for electrolytes in nanoconfinement.

Formic acid stability in different solvents by DFT calculations

CONTEXT

New technologies have been developed toward the use of green energies. The production of formic acid (FA) from carbon dioxide (CO[Formula: see text]) hydrogenation with H[Formula: see text] is a sustainable process for H[Formula: see text] storage. However, the FA adduct stabilization is thermodynamically dependent on the type of solvent and thermodynamic conditions. The results suggest a wide range of dielectric permittivity values between the dimethyl sulfoxide (DMSO) and water solvents to stabilize the FA in the absence of base. The thermodynamics analysis and the infrared and charge density difference results show that the formation of the FA complex with H[Formula: see text]O is temperature dependent and has a major influence on aqueous solvents compared to the FA adduct with amine, in good agreement with the experiment. In these conditions, the stability thermodynamic of the FA molecule may be favorable at non-organic solvents and dielectric permittivity values closer to water. Therefore, a mixture of aqueous solvents with possible ionic composition could be used to increase the thermodynamic stability of H[Formula: see text] storage in CO[Formula: see text] conversion processes.

METHODS

Using the Quantum ESPRESSO package, density functional theory (DFT) calculations were performed with periodic boundary conditions, and the electronic wave functions were expanded in plane waves. For the exchange-correlation functional, we use the vdW-DF functional with the inclusion of van der Waals (vdW) forces. Electron-ion interactions are treated by the projector augmented wave (PAW) method with pseudopotentials available in the PSlibrary repository. The wave functions and the electronic densities were expanded employing accurate cut-off energies of 6.80[Formula: see text]10[Formula: see text] and 5.44[Formula: see text]10[Formula: see text] eV, respectively. The electronic density was computed from the wave functions calculated at the [Formula: see text]-point in the first Brillouin-zone. Each structural optimization was minimized according to the Broyden-Fletcher-Goldfarb-Shanno (BFGS) algorithm, with force and energy convergence criteria of 25 meV[Formula: see text]Å[Formula: see text] and 1.36 meV, respectively. The electrostatic solvation effects were performed by the [Formula: see text] package with the Self-Consistent Continuum Solvation (SCCS) approach.

Microfluidic Printing of Vertically-Oriented Nanosheets/MOFs Hetero-Interface for Intensive Pseudocapacitive Storage

Efficient manufacture of electroactive vertically-oriented nanosheets with enhanced electrolyte mass diffusion and strong interfacial redox dynamics is critical for realizing high energy density of miniature supercapacitor (SC), but still challenging. Herein, microfluidic droplet printing is developed to controllably construct vertically-oriented graphene/ZIF-67 hetero-microsphere (VAGS/ZIF-67), where the ZIF-67 is coordinately grown on vertically-oriented graphene framework via Co─O─C bonds. The VAGS/ZIF-67 shows ordered porous channel, high electroactivity and strong interfacial interaction, providing rapid electrolyte diffusion dynamics and high faradaic capacitance in KOH solution (1674 F g-1 , 1004 C g-1 ), which are verified by computational fluid dynamics (CFD) and density functional theory (DFT). Moreover, the VAGS/ZIF-67 based SC exhibits large energy density (100 Wh kg-1 ), excellent durability (10 000 cycles and high/low temperature), and robust power-supply applications in portable electronics.

Differences in interaction of graphene/graphene oxide with bacterial and mammalian cell membranes

Graphene, a single layer, hexagonally packed two-dimensional carbon sheet is an attractive candidate for diverse applications including antibacterial potential and drug delivery. One of the knowledge gaps in biomedical application of graphene is the interaction of these materials with the cells. To address this, we investigated the interaction between graphene materials (graphene and graphene oxide) and plasma membranes of cells (bacterial and mammalian cells). The interactions of four of the most abundant phospholipids in bacteria and mammalian plasma membranes with graphene materials were studied using density functional theory (DFT) at the atomic level. The calculations showed that the mammalian phospholipids have stronger bonding to each other compared to bacterial phospholipids. When the graphene/graphene oxide sheet is approaching the phospholipid pairs, the bacterial pairs exhibit less repulsive interactions, thereby a more stable system with the sheets was found. We also assembled bacterial and mammalian phospholipids into liposomes. We further observed that the bacterial liposomes and cells let the graphene flakes penetrate the membrane. The differential scanning calorimetry measurements of liposomes revealed that the bacterial liposomes have the lowest heat capacity; this strengthens the theoretical predictions of weaker interaction between the bacterial phospholipids compared to the mammalian phospholipids. We further demonstrated that graphene oxide could be internalized into the mammalian liposomes without disrupting the membrane integrity. The results suggest that the weak bonding among bacteria phospholipids and less repulsive force when graphene materials approach, result in graphene materials interacting differently with the bacteria compared to mammalian cells.

Prediction of Strong Solvatochromism in a Molecular Photocatalyst

Based on quantum chemical calculations, we predict strong solvatochromism in a light-driven molecular photocatalyst for hydrogen generation, that is we show that the electronic and optical properties of the photocatalyst strongly depend on the solvent it is dissolved in. Our calculations in particular indicate a solvent-dependent relocation of the highest occupied molecular orbital (HOMO). Ground-state density functional theory and linear response time-dependent density functional theory calculations were applied in order to investigate the influence of implicit solvents on the structural, electronic and optical properties of a molecular photocatalyst. Only at high dielectric constants of the solvent, is the HOMO located at the metal center of the photosensitizer, whereas at low dielectric constants the HOMO is centered at the metal atom of the catalytically active complex. We elucidate the electronic origins of this strong solvatochromic effect and sketch the consequences of these insights for the use of photocatalysts in different environments.

An implicit electrolyte model for plane wave density functional theory exhibiting nonlinear response and a nonlocal cavity definition

We have developed and implemented an implicit electrolyte model in the Vienna Ab initio Simulation Package (VASP) that includes nonlinear dielectric and ionic responses as well as a nonlocal definition of the cavities defining the spatial regions where these responses can occur. The implementation into the existing VASPsol code is numerically efficient and exhibits robust convergence, requiring computational effort only slightly higher than the original linear polarizable continuum model. The nonlinear + nonlocal model is able to reproduce the characteristic "double hump" shape observed experimentally for the differential capacitance of an electrified metal interface while preventing "leakage" of the electrolyte into regions of space too small to contain a single water molecule or solvated ion. The model also gives a reasonable prediction of molecular solvation free energies as well as the self-ionization free energy of water and the absolute electron chemical potential of the standard hydrogen electrode. All of this, combined with the additional ability to run constant potential density functional theory calculations, should enable the routine computation of activation barriers for electrocatalytic processes.

Ab Initio-Based Modeling of Thermodynamic Cyclic Voltammograms: A Benchmark Study on Ag(100) in Bromide Solutions

Experimental cyclic voltammograms (CVs) measured in the slow scan rate limit can be entirely described in terms of the thermodynamic equilibrium quantities of the electrified solid-liquid interface. They correspondingly serve as an important benchmark for the quality of first-principles calculations of interfacial thermodynamics. Here, we investigate the partially drastic approximations made presently in computationally efficient calculations for the well-defined showcase of an Ag(100) model electrode in Br-containing electrolytes, where the nontrivial part of the CV stems from the electrosorption of Br ions. We specifically study the entanglement of common approximations in the treatment of solvation and field effects, as well as in the way macroscopic averages of the two key quantities, namely, the potential-dependent adsorbate coverage and electrosorption valency, are derived from the first-principles energetics. We demonstrate that the combination of energetics obtained within an implicit solvation model and a perturbative second order account of capacitive double layer effects with a constant-potential grand-canonical Monte Carlo sampling of the adsorbate layer provides an accurate description of the experimental CV. However, our analysis also shows that error cancellation at lower levels of theory may equally lead to good descriptions even though key underlying physics such as the disorder-order transition of the Br adlayer at increasing coverages is inadequately treated.

Cavity formation at metal-water interfaces

The free energy cost of forming a cavity in a solvent is a fundamental concept in rationalizing the solvation of molecules and ions. A detailed understanding of the factors governing cavity formation in bulk solutions has inter alia enabled the formulation of models that account for this contribution in coarse-grained implicit solvation methods. Here, we employ classical molecular dynamics simulations and multistate Bennett acceptance ratio free energy sampling to systematically study cavity formation at a wide range of metal-water interfaces. We demonstrate that the obtained size- and position-dependence of cavitation energies can be fully rationalized by a geometric Gibbs model, which considers that the creation of the metal-cavity interface necessarily involves the removal of interfacial solvent. This so-called competitive adsorption effect introduces a substrate dependence to the interfacial cavity formation energy that is missed in existing bulk cavitation models. Using expressions from scaled particle theory, this substrate dependence is quantitatively reproduced by the Gibbs model through simple linear relations with the adsorption energy of a single water molecule. Besides providing a better general understanding of interfacial solvation, this paves the way for the derivation and efficient parametrization of more accurate interface-aware implicit solvation models needed for reliable high-throughput calculations toward improved electrocatalysts.

Stay Hydrated! Impact of Solvation Phenomena on the CO2 Reduction Reaction at Pb(100) and Ag(100) surfaces

Herein, a comprehensive computational study of the impact of solvation on the reduction reaction of CO2 to formic acid (HCOOH) and carbon monoxide on Pb(100) and Ag(100) surfaces is presented. Results further the understanding of how solvation phenomena influence the adsorption energies of reaction intermediates. We applied an explicit solvation scheme harnessing a combined density functional theory (DFT)/microkinetic modeling approach for the CO2 reduction reaction. This approach reveals high selectivities for CO formation at Ag and HCOOH formation on Pb, resolving the prior disparity between ab initio calculations and experimental observations. Furthermore, the detailed analysis of adsorption energies of relevant reaction intermediates shows that the total number of hydrogen bonds formed by HCOO plays a primary role for the adsorption strength of intermediates and the electrocatalytic activity. Results emphasize the importance of explicit solvation for adsorption and electrochemical reaction phenomena on metal surfaces.

Coordination of Mercury(II) in Water Promoted over Hydrolysis in Solvated Clusters [Hg(H2O)1-6](aq)2+: Insights from Relativistic Effects and Free Energy Analysis

Understanding the nature of the interaction between mercury(II) ions, Hg2+, and water molecules is crucial to describe the stability and chemical behavior of structures formed during solvation, as well as the conditions that favor the Hg2+ coordination or inducing water hydrolysis. In our study, we explored exhaustively the potential energy surface of Hg2+ with up to six water molecules. We analyzed electronic and Gibbs free energies, binding, and nuclear magnetic resonance parameters. We used the zeroth-order regular approximation Hamiltonian, including scalar and spin-orbit relativistic corrections for free energy calculations and geometry optimizations to explore the interplay between electron correlation and relativistic effects. We analyzed intermolecular interactions with energy decomposition analysis, quantum theory of atoms in molecules, and natural bond orbital. Additionally, we used the four-component Dirac Hamiltonian to compute solvent effect on the magnetic shielding and J-coupling constants. Our results revealed that the water hydrolysis by Hg2+ requires a minimum of three water molecules. We found that the interaction between Hg2+ and water molecules is an orbital interaction due to relativistic effects and the most stable structures are opened-shape clusters, reducing the number of oxygen-mercury contacts and maximizing the formation of hydrogen bonds among water molecules. In these types of clusters, Hg2+ promotes the water hydrolysis over coordination with oxygen atoms. However, when we considered the change associated with the transfer of a cluster from the ideal gas to a solvated system, our solvation free energy analysis revealed that closed-shape clusters are more favorable, maximizing the number of oxygen-mercury contacts and reducing the formation of hydrogen bonds among water molecules. This finding suggests that, under room conditions, the coordination of Hg2+ is more favorable than hydrolysis. Our results have significant implications for understanding Hg2+ behavior in water, helping to develop targeted strategies for mercury remediation and management, and contributing to advancements in the broader field of environmental chemistry.

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.

Density functional theory study of the corrosion inhibition performance of 6-mercaptopurine and 6-thioguanine expired drugs toward the aluminium (111) surface

The potentiality of the 6-mercaptopurine (MP) and 6-thioguanine (TG) expired drugs toward the corrosion inhibition of the aluminium (Al) (111) surface was widely investigated using a series of density functional theory (DFT) calculations. A competition between the anti-corrosive features of the studied drugs in the gas and aqueous phases was conducted on both neutral and protonated forms by means of quantum mechanical descriptors. The results of the electrostatic potential analysis demonstrated the prominent nucleophilic nature of the sulfur and nitrogen atoms over the structures of the examined drugs. The frontier molecular orbital theory findings outlined the higher preferability of TG over MP as a corrosion inhibitor. Upon determining the most beneficial configurations of the MP/TG⋯Al (111) complexes, first-principles molecular dynamics simulations were executed. Interestingly, the competence of the TG drug in the corrosion inhibition process of Al (111) was more extensive than that of the MP one, which was confirmed by the interaction energy values of -1.79 and -1.64 eV, respectively. Upon obtaining the relaxed complexes, the effect of the presence of water solvent on the adsorption process was studied. These findings provide a foundation for developing green anti-corrosive inhibitors for the aluminium surface.

Thermodynamic and kinetic modeling of electrocatalytic reactions using a first-principles approach

The computational modeling of electrochemical interfaces and their applications in electrocatalysis has attracted great attention in recent years. While tremendous progress has been made in this area, however, the accurate atomistic descriptions at the electrode/electrolyte interfaces remain a great challenge. The Computational Hydrogen Electrode (CHE) method and continuum modeling of the solvent and electrolyte interactions form the basis for most of these methodological developments. Several posterior corrections have been added to the CHE method to improve its accuracy and widen its applications. The most recently developed grand canonical potential approaches with the embedded diffuse layer models have shown considerable improvement in defining interfacial interactions at electrode/electrolyte interfaces over the state-of-the-art computational models for electrocatalysis. In this Review, we present an overview of these different computational models developed over the years to quantitatively probe the thermodynamics and kinetics of electrochemical reactions in the presence of an electrified catalyst surface under various electrochemical environments. We begin our discussion by giving a brief picture of the different continuum solvation approaches, implemented within the ab initio method to effectively model the solvent and electrolyte interactions. Next, we present the thermodynamic and kinetic modeling approaches to determine the activity and stability of the electrocatalysts. A few applications to these approaches are also discussed. We conclude by giving an outlook on the different machine learning models that have been integrated with the thermodynamic approaches to improve their efficiency and widen their applicability.

Role of graphene in scavenging methyl cations: a DFT study

CONTEXT

It is known that methylating agents methylate DNA by transferring a methyl cation (CH3+) to the nucleophilic sites in DNA bases and DNA methylation is implicated in cancer and other pathological conditions. Therefore, it is important to scavenge CH3+ ion in order to protect DNA from methylation. Graphene is considered to be a versatile material for use in a wide variety of fields including sensors, antioxidants, drug delivery and DNA sequencing. In this work, we have theoretically investigated the interaction of CH3+ ions with graphene surface with an aim to understand if pristine graphene can be used as a substrate to adsorb CH3+ cations generated from harmful methylating agents. The computed adsorption energies show that adsorption of one, two and three CH3+ ions on graphene is favourable as the adducts thus formed are found to be substantially stable in both gas phase and aqueous media. The Bader charge transfer analysis and density of states (DOS) calculation also indicate a strong interaction between graphene and CH3+ ions. Thus, our results show that pristine graphene can be used as a substrate to scavenge CH3+ ions.

METHODS

The spin polarised density functional theory (DFT) calculations employing PBE functional, ultrasoft pseudopotentials and plane wave basis set having kinetic energy cut-offs of 40 Ry and 400 Ry, respectively, for wave functions and charge densities were carried out to study the adsorption of CH3+ ion(s) on the pristine graphene surface. The Grimme's DFT-D2 method was used for the estimation of van der Waals interactions. The 'dipole correction' along z-direction was also applied for adsorption study. The Marzari-Vanderbilt smearing and Monkhorst-Pack k-point grid were employed for the Brillouin zone sampling. A 6 × 6 graphene supercell with a vertical cell dimension of 18 Å was considered for the adsorption study. The charge transfer between the CH3+ ion(s) and graphene was estimated using Bader charge analysis. The implicit solvation model (SCCS) was used to estimate the solvent effect of aqueous media. All the calculations were performed using QUANTUM ESPRESSO package.

Amyloid-β Tetramers and Divalent Cations at the Membrane/Water Interface: Simple Models Support a Functional Role

Charge polarization at the membrane interface is a fundamental process in biology. Despite the lower concentration compared to the abundant monovalent ions, the relative abundance of divalent cations (Ca2+, Mg2+, Zn2+, Fe2+, Cu2+) in particular spaces, such as the neuron synapse, raised many questions on the possible effects of free multivalent ions and of the required protection of membranes by the eventual defects caused by the free forms of the cations. In this work, we first applied a recent realistic model of divalent cations to a well-investigated model of a polar lipid bilayer, di-myristoyl phosphatidyl choline (DMPC). The full atomistic model allows a fairly good description of changes in the hydration of charged and polar groups upon the association of cations to lipid atoms. The lipid-bound configurations were analyzed in detail. In parallel, amyloid-β 1-42 (Aβ42) peptides assembled into tetramers were modeled at the surface of the same bilayer. Two of the protein tetramers' models were loaded with four Cu2+ ions, the latter bound as in DMPC-free Aβ42 oligomers. The two Cu-bound models differ in the binding topology: one with each Cu ion binding each of the monomers in the tetramer; one with pairs of Cu ions linking two monomers into dimers, forming tetramers as dimers of dimers. The models here described provide hints on the possible role of Cu ions in synaptic plasticity and of Aβ42 oligomers in storing the same ions away from lipids. The release of structurally disordered peptides in the synapse can be a mechanism to recover ion homeostasis and lipid membranes from changes in the divalent cation concentration.

Adsorption Features of Tetrahalomethanes (CX4; X = F, Cl, and Br) on β12 Borophene and Pristine Graphene Nanosheets: A Comparative DFT Study

The potentiality of the β₁₂ borophene (β₁₂) and pristine graphene (GN) nanosheets to adsorb tetrahalomethanes (CX₄; X = F, Cl, and Br) were investigated using density functional theory (DFT) methods. To provide a thorough understanding of the adsorption process, tetrel (XC-X₃∙∙∙β₁₂/GN)- and halogen (X₃C-X∙∙∙β₁₂/GN)-oriented configurations were characterized at various adsorption sites. According to the energetic manifestations, the adsorption process of the CX₄∙∙∙β₁₂/GN complexes within the tetrel-oriented configuration led to more desirable negative adsorption energy (E_ads) values than that within the halogen-oriented analogs. Numerically, E_ads values of the CBr₄∙∙∙Br1@β₁₂ and T@GN complexes within tetrel-/halogen-oriented configurations were -12.33/-8.91 and -10.03/-6.00 kcal/mol, respectively. Frontier molecular orbital (FMO) results exhibited changes in the E_HOMO, E_LUMO, and E_gap values of the pure β₁₂ and GN nanosheets following the adsorption of CX₄ molecules. Bader charge transfer findings outlined the electron-donating property for the CX₄ molecules after adsorbing on the β₁₂ and GN nanosheets within the two modeled configurations, except the adsorbed CBr₄ molecule on the GN sheet within the tetrel-oriented configuration. Following the adsorption process, new bands and peaks were observed in the band structure and density of state (DOS) plots, respectively, with a larger number in the case of the tetrel-oriented configuration than in the halogen-oriented one. According to the solvent effect affirmations, adsorption energies of the CX₄∙∙∙β₁₂/GN complexes increased in the presence of a water medium. The results of this study will serve as a focal point for experimentalists to better comprehend the adsorption behavior of β₁₂ and GN nanosheets toward small toxic molecules.

The Oxygen Evolution Reaction at MoS2 Edge Sites: The Role of a Solvent Environment in DFT-Based Molecular Simulations

Density functional theory (DFT) calculations are employed to study the oxygen evolution reaction (OER) on the edges of stripes of monolayer molybdenum disulfide. Experimentally, this material has been shown to evolve oxygen, albeit with low efficiency. Previous DFT studies have traced this low catalytic performance to the unfavourable adsorption energies of some reaction intermediates on the MoS₂ edge sites. In this work, we study the effects of the aqueous liquid surrounding the active sites. A computational approach is used, where the solvent is modeled as a continuous medium providing a dielectric embedding of the catalyst and the reaction intermediates. A description at this level of theory can have a profound impact on the studied reactions: the calculated overpotential for the OER is lowered from 1.15 eV to 0.77 eV. It is shown that such variations in the reaction energetics are linked to the polar nature of the adsorbed intermediates, which leads to changes in the calculated electronic charge density when surrounded by water. These results underline the necessity to computationally account for solvation effects, especially in aqueous environments and when highly polar intermediates are present.

Solv: An Alternative Continuum Model Implementation Based on Fixed Atomic Charges, Scaled Particle Theory, and the Atom-Atom Potential Method

An alternative continuum model implementation is reported. The electrostatic contribution to the solvation Gibbs free energy utilizes the noniterative conductor-like screening model of Vyboishchikov and Voityuk (DOI: 10.1002/jcc.26531) based on the fixed partial atomic charges. The nonelectrostatic solute-solvent dispersion-repulsion energy is calculated through the Caillet-Claverie atom-atom potential method employing the grid-based approach. The nonelectrostatic cavitation energy is computed within the scaled particle theory (SPT) formalism with the solute hard-sphere radius obtained via the Pierotti-Claverie (PC) scheme, from the solute molecular surface (SPT-S) or volume (SPT-V). The solvent hard-sphere radius is derived through the fitting to the experimental total solvation free energies of 2530 neutral species in 92 solvents. Application of the model to reproduce both the absolute and relative (reaction net) solvation free energies indicates that the SPT-V approach based on the CM5 charges is the best performer. The method is suggested for the solvation free energy calculation in the nonaqueous solvents.

On the Use of Graphene Nanosheets for Drug Delivery: A Case Study of Cisplatin and Some of Its Analogs

Graphene (GN) nanosheets have been widely exploited in biomedical applications as potential nanocarriers for various drugs due to their distinct physical and chemical properties. In this regard, the adsorption behavior of cisplatin (cisPtCl₂) and some of its analogs on a GN nanosheet was investigated in perpendicular and parallel configurations by using density functional theory (DFT). According to the findings, the most significant negative adsorption energies (E_ads) within the cisPtX₂⋯GN complexes (where X = Cl, Br, and I) were observed for the parallel configuration, with values up to -25.67 kcal/mol at the H@GN site. Within the perpendicular configuration of the cisPtX₂⋯GN complexes, three orientations were investigated for the adsorption process, namely, X/X, X/NH₃, and NH₃/NH₃. The negative E_ads values of the cisPtX₂⋯GN complexes increased with the increasing atomic weight of the halogen atom. The Br@GN site showed the largest negative E_ads values for the cisPtX₂⋯GN complexes in the perpendicular configuration. The Bader charge transfer outcomes highlighted the electron-accepting properties of cisPtI₂ within the cisPtI₂⋯GN complexes in both configurations. The electron-donating character of the GN nanosheet increased as the electronegativity of the halogen atom increased. The band structure and density of state plots revealed the occurrence of the physical adsorption of the cisPtX₂ on the GN nanosheet, which was indicated by the appearance of new bands and peaks. Based on the solvent effect outlines, the negative E_ads values generally decreased after the adsorption process in a water medium. The recovery time results were in line with the E_ads findings, where the cisPtI₂ in the parallel configuration took the longest time to be desorbed from the GN nanosheet with values of 61.6 × 10⁸ ms at 298.15 K. The findings of this study provide better insights into the utilization of GN nanosheets in drug delivery applications.

Concerted Proton-Coupled Electron Transfer to a Graphite Adsorbed Metalloporphyrin Occurs by Band to Bond Electron Redistribution

Surface immobilized catalysts are highly promising candidates for a range of energy conversion reactions, and atomistic mechanistic understanding is essential for their rational design. Cobalt tetraphenylporphyrin (CoTPP) nonspecifically adsorbed on a graphitic surface has been shown to undergo concerted proton-coupled electron transfer (PCET) in aqueous solution. Herein, density functional theory calculations on both cluster and periodic models representing π-stacked interactions or axial ligation to a surface oxygenate are performed. As the electrode surface is charged due to applied potential, the adsorbed molecule experiences the electrical polarization of the interface and nearly the same electrostatic potential as the electrode, regardless of the adsorption mode. PCET occurs by electron abstraction from the surface to the CoTPP concerted with protonation to form a cobalt hydride, thereby circumventing Co(II/I) redox. Specifically, the Co(II) d-state localized orbital interacts with a proton from solution and an electron from the delocalized graphitic band states to produce a Co(III)-H bonding orbital below the Fermi level, corresponding to redistribution of electrons from the band states to the bonding states. These insights have broad implications for electrocatalysis by chemically modified electrodes and surface immobilized catalysts.

The importance of a charge transfer descriptor for screening potential CO2 reduction electrocatalysts

It has been over twenty years since the linear scaling of reaction intermediate adsorption energies started to coin the fields of heterogeneous and electrocatalysis as a blessing and a curse at the same time. It has established the possibility to construct activity volcano plots as a function of a single or two readily accessible adsorption energies as descriptors, but also limited the maximal catalytic conversion rate. In this work, it is found that these established adsorption energy-based descriptor spaces are not applicable to electrochemistry, because they are lacking an important additional dimension, the potential of zero charge. This extra dimension arises from the interaction of the electric double layer with reaction intermediates which does not scale with adsorption energies. At the example of the electrochemical reduction of CO₂ it is shown that the addition of this descriptor breaks the scaling relations, opening up a huge chemical space that is readily accessible via potential of zero charge-based material design. The potential of zero charge also explains product selectivity trends of electrochemical CO₂ reduction in close agreement with reported experimental data highlighting its importance for electrocatalyst design.

Characterizing Grain Boundary Effects on Mg2+ Conduction in Metal-Organic Frameworks

Next-generation materials for fast ion conduction have the potential to revolutionize battery technology. Metal-organic frameworks (MOFs) are promising candidates for achieving this goal. Given their structural diversity, the design of efficient MOF-based conductors can be accelerated by a detailed understanding and accurate prediction of ion conductivity. However, the polycrystalline nature of solid-state materials requires consideration of grain boundary effects, which is complicated by challenges in characterizing grain boundary structures and simulating ensemble transport processes. To address this, we have developed an approach for modeling ion transport at grain boundaries and predicting their contribution to conductivity. Mg²⁺ conduction in the Mg-MOF-74 thin film was studied as a representative system. Using computational techniques and guided by experiments, we investigated the structural details of MOF grain boundary interfaces to determine accessible Mg²⁺ transport pathways. Computed transport kinetics were input into a simplified MOF nanocrystal model, which combines ion transport in the bulk structure and at grain boundaries. The model predicts Mg²⁺ conductivity in the MOF-74 film within chemical accuracy (<1 kcal/mol activation energy difference), validating our approach. Physically, Mg²⁺ conduction in MOF-74 is inhibited by strong Mg²⁺ binding at grain boundaries, such that only a small fraction of grain boundary alignments allow for fast Mg²⁺ transport. This results in a 2-3 order-of-magnitude reduction in conductivity, illustrating the critical impact of the grain boundary contribution. Overall, our work provides a computation-aided platform for molecular-level understanding of grain boundary effects and quantitative prediction of ion conductivity. Combined with experimental measurements, it can serve as a synergistic tool for characterizing the grain boundary composition of MOF-based conductors.

Inductive Effect Alone Cannot Explain Lewis Adduct Formation and Dissociation at Electrode Interfaces

Understanding breaking and formation of Lewis bonds at an electrified interface is relevant to a large range of phenomena, including electrocatalysis and electroadsorption. The complexities of interfacial environments and associated reactions often impede a systematic understanding of this type of bond at interfaces. To address this challenge, we report the creation of a main group classic Lewis acid-base adduct on an electrode surface and its behavior under varying electrode potentials. The Lewis base is a self-assembled monolayer of mercaptopyridine and the Lewis acid is BF₃, forming a Lewis bond between nitrogen and boron. The bond is stable at positive potentials but cleaves at potentials more negative of approximately -0.3 V vs Ag/AgCl without an associated current. We also show that if the Lewis acid BF₃ is supplied from a reservoir of Li⁺BF₄^- electrolyte, the cleavage is completely reversible. We propose that the N-B Lewis bond is affected both by the field-induced intramolecular polarization (electroinduction) and by the ionic structures and ionic equilibria near the electrode. Our results indicate that the second effect is responsible for the Lewis bond cleavage at negative potentials. This work is relevant to understanding the fundamentals of electrocatalytic and electroadsorption processes.

Adaptive Electric Fields Embedded Electrochemical Barrier Calculations

Electrochemical interfaces are grand canonical ensembles of varying electrons. Simulating them by standard first-principles methods is a challenging task, since the number of electrons (or charge) is fixed in the calculation. Under the constant charge framework, we developed a constant potential simulation method realized by adding an adaptive electric field to a charge neutral system. Electric field is the controlling variable. In addition, we defined an internal reversible hydrogen electrode potential (ϕ_IRHE), which can ensure the model independence of our method. To validate our method, the reaction energies of some electrochemical reactions are calculated, the results are comparable with the computational hydrogen electrode model and experiments. At last, the evolution of transition state structures and charge transfer coefficients of some electrochemical reactions on Ag(111) surface were discussed by our method.

Theoretical Screening of CO2 Electroreduction over MOF-808-Supported Self-Adaptive Dual-Metal-Site Pairs

Electrochemical CO₂ reduction to transportation fuels and valuable platform chemicals provides a sustainable avenue for renewable energy storage and realizes an artificially closed carbon loop. However, the rational design of highly active and selective CO₂ reduction electrocatalysts remains a challenging task. Herein, a series of metal-organic framework (MOF)-supported flexible, self-adaptive dual-metal-site pairs (DMSPs) including 21 pairwise combinations of six transition metal single sites (MOF-808-EDTA-M₁M₂, M₁/M₂ = Fe, Cu, Ni, Pd, Pt, Au) for the CO₂ reduction reaction (CO₂RR) were theoretically screened using density functional theory calculations. Against the competitive hydrogen evolution reaction, MOF-808-EDTA-FeFe and MOF-808-EDTA-FePt were identified as the promising CO₂RR electrocatalysts toward C₁ and C₂ products. The calculated limiting potential for CO₂ electroreduction to C₂H₆ and C₂H₅OH over MOF-808-EDTA-FeFe is -0.87 V. Compared with an applied potential of -0.56 eV toward CH₄ production over MOF-808-EDTA-FeFe, MOF-808-EDTA-FePt exhibits an even better activity for CO₂ reduction to C₁ products at a limiting potential of -0.35 V. The present work not only identifies promising candidates for highly selective CO₂RR electrocatalysts leading to C₁ and C₂ products but also provides mechanistic insights into the dynamic nature of DMSPs for stabilizing various reaction intermediates in the CO₂RR process.

Physical and numerical aspects of sodium ion solvation free energies via the cluster-continuum model

Sodium cation solvation Gibbs free energies (ΔG_solv(Na⁺)) have been obtained in water, dimethylformamide, dimethyl sulfoxide, ethanol, acetone, acetonitrile, and methanol through the "monomer cycle" cluster-continuum approach where a solvent reference state is described by infinitely separated molecules. The following steps are vital for obtaining reliable ΔG_solv(Na⁺) values: (a) a meticulous conformational search involving dispersion corrected density functional theory (DFT-D) and the continuum solvation model (CSM); (b) gas-phase DFT-D geometry optimization followed by single-point (SP) domain-based local pair natural orbital coupled clusters including single, double, and partly triple excitation (DLPNO-CCSD(T)) calculations in conjunction with the complete basis set extrapolation; (c) advanced statistical thermodynamic treatment of the low harmonic frequencies (<100 cm^-1) to obtain the robust gas-phase Gibbs free energy correction; (d) gas-phase and dielectric continuum SP with non-electrostatic contributions included in the CSM; (e) an evaluation of the relative thermodynamic stability of the Na⁺(S)_n clusters to identify the number of explicit solvent molecules n to be considered. Our refined computational protocol is promising with a Pearson correlation coefficient between the predicted and experimental data, ρ, of 0.82, and the mean signed and mean unsigned errors of 0.3 and 1.4 kcal mol^-1, respectively.

Unravelling the 2e- ORR Activity Induced by Distance Effect on Main-Group Metal InN4 Surface Based on First Principles

The p-electron-dominated main-group metals (Sb, Se, In, etc.) have recently been reported to possess excellent oxygen reduction reaction (ORR) activity by means of heteroatom doping into graphene. However, on these main group metal surfaces, other approaches especially the distance effect to modulate catalytic activity are rarely involved. In this work, the origin of excellent 2e^- ORR catalytic activity of graphene-supported InN₄ moiety by tuning the distance between metallic In atoms is thoroughly investigated by employing the first-principles calculations. Our DFT calculations show that the 2e^- ORR catalytic activity strongly depends on the crystal orbital Hamilton population (COHP) between In and O atoms. This work is useful for the rational design of main group metal single atom electrocatalysts.

Ion Solvation Free Energy Calculation Based on Ab Initio Molecular Dynamics Using a Hybrid Solvent Model

Free energy calculation of small molecules or ion species in aqueous solvent is one of the most important problems in electrochemistry study. Although there are many previous approaches to calculate such free energies, they are based on ab initio methods and suffer from various limitations and approximations. In the current work, we developed a hybrid approach based on ab initio molecular dynamics (AIMD) simulations to calculate the ion solvation energy. In this approach, a small water cluster surrounding the central ion is used, and implicit solvent model is applied outside the water cluster. A dynamic potential well is used during AIMD to keep the water cluster together. Quasi-harmonic approximation is used to calculate the entropy contribution, while the total energy average is used to calculate the enthalpy term. The obtained solvation voltages of the bulk metal agree with experiments within 0.3 eV, and the simulation results for the solvation energies of gaseous ions are close to the experimental observations. Besides the free energies, radial pair distribution functions and coordination numbers of hydrated cations are also obtained. The remaining challenges of this method are also discussed.

A multiscale approach for electronic transport simulation of carbon nanostructures in aqueous solvent

Theoretical works addressing electronic nano-devices operating in an aqueous environment often neglect solvent effects. In order to assess the role played by the polarization effects on the electronic transport properties of solvated graphene, for example in possible bio-sensing applications, we have used here a combination of polarizable force-field molecular dynamics, hybrid quantum mechanics/molecular mechanics (QM/MM) approach, density functional theory, and non-equilibrium Green's function method. We considered different solvation conditions, the presence of defects in graphene, as well as various choices for the partitions between the quantum and classical regions in QM/MM, in which we explicitly account for polarization effects. Our results show that the polarization effects on graphene lead to changes in the structure of interfacial water molecules which are more pronounced in the vicinity of defects. The presence of water leads to increased scattering due to the long-range charge interactions with graphene. At the same time, changes in the conductance due to polarization or salt concentration are found to be small, paving the way for robust electronic nano-devices operating in aqueous environments.

Multilevel Computational Studies Reveal the Importance of Axial Ligand for Oxygen Reduction Reaction on Fe-N-C Materials

The systematic improvement of Fe-N-C materials for fuel cell applications has proven challenging, due in part to an incomplete atomistic understanding of the oxygen reduction reaction (ORR) under electrochemical conditions. Herein, a multilevel computational approach, which combines ab initio molecular dynamics simulations and constant potential density functional theory calculations, is used to assess proton-coupled electron transfer (PCET) processes and adsorption thermodynamics of key ORR intermediates. These calculations indicate that the potential-limiting step for ORR on Fe-N-C materials is the formation of the Fe^III-OOH intermediate. They also show that an active site model with a water molecule axially ligated to the iron center throughout the catalytic cycle produces results that are consistent with the experimental measurements. In particular, reliable prediction of the ORR onset potential and the Fe(III/II) redox potential associated with the conversion of Fe^III-OH to Fe^II and desorbed H₂O requires an axial H₂O co-adsorbed to the iron center. The observation of a five-coordinate rather than four-coordinate active site has significant implications for the thermodynamics and mechanism of ORR. These findings highlight the importance of solvent-substrate interactions and surface charge effects for understanding the PCET reaction mechanisms and transition-metal redox couples under realistic electrochemical conditions.

GW Quasiparticle Energies and Bandgaps of Two-Dimensional Materials Immersed in Water

Computational simulations have become of major interest to screen potential photocatalysts for optimal band edge positions which straddle the redox potentials. Unfortunately, these methods suffer from a difficulty in resolving the dynamic solvent response on the band edge positions. We have developed a computational method based on the GW approximation coupled with an implicit solvation model that solves a generalized Poisson equation (GPE), that is, GW-GPE. Using GW-GPE, we have investigated the band edge locations of (quasi) 2D materials immersed in water and found a good agreement with experimental data. We identify two contributions of the solvent effect, termed a "polarization-field effect" and an "environmental screening effect", which are found to be highly sensitive to the atomic and charge distribution of the 2D materials. We believe that the GW-GPE scheme can pave the way to predict band edge positions in solvents, enabling design of 2D material-based photocatalysts and energy systems.

Predicting the Mechanisms for H2O2 Activation and Phenol Oxidation Catalyzed by Modified Graphene-Based Systems Using Density Functional Theory

The heterogeneous Fenton-like reaction on metal-free graphene-based catalysts attracts great attention. However, a systematic and comprehensive understanding of the mechanisms for H₂O₂ activation and pollutant oxidation is still lacking. In this study, the heterogeneous Fenton-like mechanisms on doped and oxygen-containing graphene are investigated using density functional theory. The H₂O₂ tends to form a surface oxygen and a water molecule on the doped graphene. For the oxygen-containing graphene-based systems, relative to the groups in the basal plane, the separated groups on the edge including hydroxyl, carbonyl, and carboxyl readily activate H₂O₂ to hydroxyls. However, when the groups are close to each other, more additional side reactions might occur upon H₂O₂ adsorption, which may inhibit catalyst retrieval. Phenol is selected as a model pollutant to study its oxidation reaction with the adsorbed oxygen formed from the dissociated H₂O₂. The thermodynamics of the reactions depends significantly on the co-adsorption strengths over different catalysts. Our work provides key fundamental insights into the catalytic performance of various modified graphene-based systems, which could guide the future design and applications of heterogeneous Fenton reactions.