A critical assessment of theoretical methods for finding reaction pathways and transition states of surface processes

Plasma Chemical Looping: Unlocking High-Efficiency CO2 Conversion to Clean CO at Mild Temperatures

We propose a plasma chemical looping CO2 splitting (PCLCS) approach that enables highly efficient CO2 conversion into O2-free CO at mild temperatures. PCLCS achieves an impressive 84% CO2 conversion and a 1.3 mmol g-1 CO yield, with no O2 detected. Crucially, this strategy significantly lowers the temperature required for conventional chemical looping processes from 650 to 1000 °C to only 320 °C, demonstrating a robust synergy between plasma and the Ce0.7Zr0.3O2 oxygen carrier (OC). Systematic experiments and density functional theory (DFT) calculations unveil the pivotal role of plasma in activating and partially decomposing CO2, yielding a mixture of CO, O2/O, and electronically/vibrationally excited CO2*. Notably, these excited CO2* species then efficiently decompose over the oxygen vacancies of the OCs, with a substantially reduced activation barrier (0.86 eV) compared to ground-state CO2 (1.63 eV), contributing to the synergy. This work offers a promising and energy-efficient pathway for producing O2-free CO from inert CO2 through the tailored interplay of plasma and OCs.

Kondo Effect of Co-Porphyrin: Remarkable Sensitivity to Adsorption Sites and Orientations

We investigated the Kondo effect of cobalt(II)-5-15-bis(4'-bromophenyl)-10,20-bis(4'-iodophenyl)porphyrin (CoTPPBr2I2) molecules on Au(111) with low-temperature scanning tunneling microscopy under ultrahigh vacuum conditions. The molecules exhibit four adsorption configurations at the top and bridge sites of the surface with different molecular orientations. The Kondo resonance shows extraordinary sensitivity to the adsorption configuration. By switching the molecule between different configurations, the Kondo temperature is varied over a wide range from ≈8 up to ≈250 K. Density functional theory calculations reveal that changes of the adsorption configuration lead to distinct variations of the hybridization between the molecule and the surface. Furthermore, we show that surface reconstruction plays a significant role for the molecular Kondo effect.

Surface Coverage as an Important Parameter for Predicting Selectivity Trends in Electrochemical CO2 Reduction

The electrochemical CO₂ reduction reaction (CO₂RR) is important for a sustainable future. Key insights into the reaction pathways have been obtained by density functional theory (DFT) analysis, but so far, DFT has been unable to give an overall understanding of selectivity trends without important caveats. We show that an unconsidered parameter in DFT models of electrocatalysts-the surface coverage of reacting species-is crucial for understanding the CO₂RR selectivities for different surfaces. Surface coverage is a parameter that must be assumed in most DFT studies of CO₂RR electrocatalysts, but so far, only the coverage of nonreacting adsorbates has been treated. Explicitly treating the surface coverage of reacting adsorbates allows for an investigation that can more closely mimic operating conditions. Furthermore, and of more immediate importance, the use of surface coverage-dependent adsorption energies allows for the extraction of ratios of adsorption energies of CO₂RR intermediates (COOH_ads and HCOO_ads) that are shown to be predictive of selectivity and are not susceptible to systematic errors. This approach allows for categorization of the selectivity of several monometallic catalysts (Pt, Pd, Au, Ag, Zn, Cu, Rh, W, Pb, Sn, In, Cd, and Tl), even problematic ones such as Ag or Zn, and does so by only considering the adsorption energies of known intermediates. The selectivity of the further reduction of COOH_ads can now be explained by a preference for Tafel or Heyrovsky reactions, recontextualizing the nature of selectivity of some catalysts. In summary, this work resolves differences between DFT and experimental studies of the CO₂RR and underlines the importance of surface coverage.

Kinetic Monte Carlo simulations for heterogeneous catalysis: Fundamentals, current status, and challenges

Kinetic Monte Carlo (KMC) simulations in combination with first-principles (1p)-based calculations are rapidly becoming the gold-standard computational framework for bridging the gap between the wide range of length scales and time scales over which heterogeneous catalysis unfolds. 1p-KMC simulations provide accurate insights into reactions over surfaces, a vital step toward the rational design of novel catalysts. In this Perspective, we briefly outline basic principles, computational challenges, successful applications, as well as future directions and opportunities of this promising and ever more popular kinetic modeling approach.

Formation of Graphene Nanoscrolls and Their Electronic Structures Based on Ab Initio Calculations

Rolling up two-dimensional (2D) materials can form quasi-one-dimensional nanoscrolls, which are expected to have novel properties due to their larger space of structural parameters. In this Letter, the structural dependence of formation energy was investigated based on more than 90 different graphene nanoscrolls (GNSs) through ab initio calculations. A quantified relationship between formation energy and structural parameters is discovered, which could provide universal description of rolling up 2D materials beyond graphene. Further calculations on electronic structures show the opening of bandgap in GNSs with ultrahigh carrier mobilities up to 10⁷ cm² V^-1 s^-1. The structural stability under room temperature was also testified by using molecular dynamic simulations. This work provides general insights into the rolling-up strategy and demonstrates the tunable properties of GNSs, thus extending the scope of the research field for 2D materials.

SnS2 Nanosheets Anchored on Nitrogen and Sulfur Co-Doped MXene Sheets for High-Performance Potassium-Ion Batteries

Potassium-ion batteries (KIBs) are emerging as the prospective alternatives to lithium-ion batteries in energy storage systems owing to the sufficient resources and relatively low cost of K-related materials. However, serious volume expansion and low specific capacity are found in most materials systems resulting from the large intrinsic radius of K⁺. Herein, SnS₂ nanosheets anchored on nitrogen and sulfur co-doped MXene (SnS₂ NSs/MXene) are creatively designed as advanced anode materials for KIBs. SnS₂ NSs/MXene with a unique hierarchical structure can not only provide fast transmission channels for K⁺ but also avoid the accumulation of K⁺ and volume expansion. These novel features make SnS₂ NSs/MXene electrodes exhibit a superior reversible specific capacity of 342.4 mA h g^-1 under 50 mA g^-1. Also, they maintain 206.1 mA h g^-1 at an even higher current density of 0.5 A g^-1 over 800 cycles almost without capacity decay. Moreover, the multistep alloying reaction mechanism of SnS₂ NSs/MXene composites and K⁺ is revealed by the ex situ X-ray diffraction measurement. In addition, the density functional theory calculations confirm the existence of Ti-S bonds between SnS₂ nanosheets and MXene, which significantly enhance the structural stability and cycling electrochemical performance of SnS₂ NSs/MXene composites.

Unravelling phosphate adsorption on hydrous ferric oxide surfaces at the molecular level

The thorough understanding of the adsorption mechanism of phosphate on hydrous ferric oxides is necessary to deal with the environmental issues related to high phosphate concentrations in soils and open water. In this work, we consider three different adsorption geometries (monodentate and bidentate chemisorption and physisorption) and calculate the adsorption geometries and related adsorption energies at optPBE-vdW level. Using the Maxwell-Boltzmann distribution, it is estimated that about 83% of the phosphate molecules is in a monodentate chemisorption configuration, while 17% is physisorbed. Furthermore, theoretical infra-red spectra are obtained and compared to equivalent experimental spectra, supporting the conclusion that mainly monodentate chemisorption and physisorption occur. Most interestingly, a weighed infra-red spectrum is then calculated, using the weights from the Maxwell-Boltzmann distribution, showing a very good comparison with the experimental spectra.

Exploring dynamic interactions of single nanoparticles at interfaces for surface-confined electrochemical behavior and size measurement

With the development of new instruments and methodologies, the highly dynamic behaviors of nanoparticle at the liquid-solid interface have been studied. However, the dynamic nature of the electrochemical behavior of individual nanoparticles on the electrode interface is still poorly understood. Here, we generalize scaling relations to predict nanoparticle-electrode interactions by examining the adsorption energy of nanoparticles at an ultramicroelectrode interface. Based on the theoretical predictions, we investigate the interaction-modulated dynamic electrochemical behaviors for the oxidation of individual Ag nanoparticles. Typically, significantly distinct current traces are observed owing to the adsorption-mediated motion of Ag nanoparticles. Inspired by restraining the stochastic paths of particles in the vicinity of the electrode interface to produce surface-confined current traces, we successfully realize high-resolution size measurements of Ag nanoparticles in mixed-sample systems. This work offers a better understanding of dynamic interactions of nanoparticles at the electrochemical interface and displays highly valuable applications of single-entity electrochemistry.

Single-entity electrochemistry has been proposed for studying properties of single nanoparticles (NPs). Here, the authors make use of adsorption-mediated motion of Ag NPs to take individual NP size measurements using electrochemical impacts with excellent agreement to standard imaging techniques.

Highly Promoted Carrier Mobility and Intrinsic Stability by Rolling Up Monolayer Black Phosphorus into Nanoscrolls

Rolling up two-dimensional (2D) materials into nanoscrolls could not only retain the excellent properties of their 2D hosts but also display intriguing physical and chemical properties that arise from their 1D tubular structures. Here, we report a new class of black phosphorus nanoscrolls (bPNSs), which are stable at room-temperature and energetically more favorable than 2D bP. Most strikingly, these bPNSs hold tunable direct band gaps and extremely high mobilities (e.g., the mobility of the double-layer bPNS is about 20-fold higher than that of 2D bP monolayer). Their unique self-encapsulation structure and layer-dependent conduction band minimum can largely prevent the entrance of O₂ and the production of O₂^- and thereby suppress the possible environmental degradation as well. The enhanced intrinsic stability and promoted electronic properties render bPNSs great promise in many advanced electronics or optoelectronics applications.

Exciting H2 Molecules for Graphene Functionalization

Hydrogen functionalization of graphene by exposure to vibrationally excited H₂ molecules is investigated by combined scanning tunneling microscopy, high-resolution electron energy loss spectroscopy, X-ray photoelectron spectroscopy measurements, and density functional theory calculations. The measurements reveal that vibrationally excited H₂ molecules dissociatively adsorb on graphene on Ir(111) resulting in nanopatterned hydrogen functionalization structures. Calculations demonstrate that the presence of the Ir surface below the graphene lowers the H₂ dissociative adsorption barrier and allows for the adsorption reaction at energies well below the dissociation threshold of the H-H bond. The first reacting H₂ molecule must contain considerable vibrational energy to overcome the dissociative adsorption barrier. However, this initial adsorption further activates the surface resulting in reduced barriers for dissociative adsorption of subsequent H₂ molecules. This enables functionalization by H₂ molecules with lower vibrational energy, yielding an avalanche effect for the hydrogenation reaction. These results provide an example of a catalytically active graphene-coated surface and additionally set the stage for a re-interpretation of previous experimental work involving elevated H₂ background gas pressures in the presence of hot filaments.

Spectroscopic Identification of the Au-C Bond Formation upon Electroreduction of an Aryl Diazonium Salt on Gold

Electroreduction of aryl diazonium salts on gold can produce organic films that are more robust than their analogous self-assembled monolayers formed from chemical adsorption of organic thiols on gold. However, whether the enhanced stability is due to the Au-C bond formation remains debated. In this work, we report the electroreduction of an aryl diazonium salt of 4,4'-disulfanediyldibenzenediazonium on gold forming a multilayer of Au-(Ar-S-S-Ar)_n, which can be further degraded to a monolayer of Au-Ar-S^- by electrochemical cleavage of the S-S moieties within the multilayer. By conducting an in situ surface-enhanced Raman spectroscopic study of both the multilayer formation/degradation and the monolayer reduction/oxidation processes, coupled to density functional theory calculations, we provide compelling evidence that an Au-C bond does form upon electroreduction of aryl diazonium salts on gold and that the enhanced stability of the electrografted organic films is due to the Au-C bond being intrinsically stronger than the Au-S bond for a given phenylthiolate compound by ca. 0.4 eV.

Probing Carrier Transport and Structure-Property Relationship of Highly Ordered Organic Semiconductors at the Two-Dimensional Limit

One of the basic assumptions in organic field-effect transistors, the most fundamental device unit in organic electronics, is that charge transport occurs two dimensionally in the first few molecular layers near the dielectric interface. Although the mobility of bulk organic semiconductors has increased dramatically, direct probing of intrinsic charge transport in the two-dimensional limit has not been possible due to excessive disorders and traps in ultrathin organic thin films. Here, highly ordered single-crystalline mono- to tetralayer pentacene crystals are realized by van der Waals (vdW) epitaxy on hexagonal BN. We find that the charge transport is dominated by hopping in the first conductive layer, but transforms to bandlike in subsequent layers. Such an abrupt phase transition is attributed to strong modulation of the molecular packing by interfacial vdW interactions, as corroborated by quantitative structural characterization and density functional theory calculations. The structural modulation becomes negligible beyond the second conductive layer, leading to a mobility saturation thickness of only ∼3  nm. Highly ordered organic ultrathin films provide a platform for new physics and device structures (such as heterostructures and quantum wells) that are not possible in conventional bulk crystals.

Efficiently computing pathway free energies: New approaches based on chain-of-replica and Non-Boltzmann Bennett reweighting schemes

BACKGROUND

Accurately modeling condensed phase processes is one of computation's most difficult challenges. Include the possibility that conformational dynamics may be coupled to chemical reactions, where multiscale (i.e., QM/MM) methods are needed, and this task becomes even more daunting.

METHODS

Free energy simulations (i.e., molecular dynamics), multiscale modeling, and reweighting schemes.

RESULTS

Herein, we present two new approaches for mitigating the aforementioned challenges. The first is a new chain-of-replica method (off-path simulations, OPS) for computing potentials of mean force (PMFs) along an easily defined reaction coordinate. This development is coupled with a new distributed, highly-parallel replica framework (REPDstr) within the CHARMM package. Validation of these new schemes is carried out on two processes that undergo conformational changes. First is the simple torsional rotation of butane, while a much more challenging glycosidic rotation (in vacuo and solvated) is the second. Additionally, a new approach that greatly improves (i.e., possibly an order of magnitude) the efficiency of computing QM/MM PMFs is introduced and compared to standard schemes. Our efforts are grounded in the recently developed method for efficiently computing QM-based free energies (i.e., QM-Non-Boltzmann Bennett, QM-NBB). Again, we validate this new technique by computing the QM/MM PMF of butane's torsional rotation.

CONCLUSIONS

The OPS-REPDstr method is a promising new approach that overcomes many limitations of standard pathway simulations in CHARMM. The combination of QM-NBB with pathway techniques is very promising as it offers significant advantages over current procedures.

GENERAL SIGNIFICANCE

Efficiently computing potentials of mean force is a major, unresolved, area of interest. This article is part of a Special Issue entitled Recent developments of molecular dynamics.

Dispersion interactions in room-temperature ionic liquids: results from a non-empirical density functional

The role of dispersion or van de Waals (VDW) interactions in imidazolium-based room-temperature ionic liquids is studied within the framework of density functional theory, using a recently developed non-empirical functional [M. Dion, H. Rydberg, E. Schröder, D. C. Langreth, and B. I. Lundqvist, Phys. Rev. Lett. 92, 246401 (2004)], as efficiently implemented in the SIESTA code [G. Román-Pérez and J. M. Soler, Phys. Rev. Lett. 103, 096102 (2009)]. We present results for the equilibrium structure and lattice parameters of several crystalline phases, finding a general improvement with respect to both the local density (LDA) and the generalized gradient approximations (GGA). Similar to other systems characterized by VDW bonding, such as rare gas and benzene dimers as well as solid argon, equilibrium distances and volumes are consistently overestimated by ≈7%, compared to -11% within LDA and 11% within GGA. The intramolecular geometries are retained, while the intermolecular distances and orientations are significantly improved relative to LDA and GGA. The quality is superior to that achieved with tailor-made empirical VDW corrections ad hoc [M. G. Del Pópolo, C. Pinilla, and P. Ballone, J. Chem. Phys. 126, 144705 (2007)]. We also analyse the performance of an optimized version of this non-empirical functional, where the screening properties of the exchange have been tuned to reproduce high-level quantum chemical calculations [J. Klimes, D. Bowler, and A. Michaelides, J. Phys.: Condens. Matter 22, 074203 (2010)]. The results for solids are even better with volumes and geometries reproduced within 2% of experimental data. We provide some insight into the issue of polymorphism of [bmim][Cl] crystals, and we present results for the geometry and energetics of [bmim][Tf] and [mmim][Cl] neutral and charged clusters, which validate the use of empirical force fields.

An enhanced splined saddle method

We present modifications for the method recently developed by Granot and Baer [J. Chem. Phys. 128, 184111 (2008)]. These modifications significantly enhance the efficiency and reliability of the method. In addition, we discuss some specific features of this method. These features provide important flexibilities which are crucial for a double-ended saddle point search method in order to be applicable to complex reaction mechanisms. Furthermore, it is discussed under what circumstances this methods might fail to find the transition state and remedies to avoid such situations are provided. We demonstrate the performance of the enhanced splined saddle method on several examples with increasing complexity, isomerization of ammonia, ethane and cyclopropane molecules, tautomerization of cytosine, the ring opening of cyclobutene, the Stone-Wales transformation of the C(60) fullerene, and finally rolling a small NaCl cube on NaCl(001) surface. All of these calculations are based on density functional theory. The efficiency of the method is remarkable in regard to the reduction of the total computational time.

Real-world predictions from ab initio molecular dynamics simulations

In this review we present the techniques of ab initio molecular dynamics simulation improved to its current stage where the analysis of existing processes and the prediction of further chemical features and real-world processes are feasible. For this reason we describe the relevant developments in ab initio molecular dynamics leading to this stage. Among them, parallel implementations, different basis set functions, density functionals, and van der Waals corrections are reported. The chemical features accessible through AIMD are discussed. These are IR, NMR, as well as EXAFS spectra, sampling methods like metadynamics and others, Wannier functions, dipole moments of molecules in condensed phase, and many other properties. Electrochemical reactions investigated by ab initio molecular dynamics methods in solution, on surfaces as well as complex interfaces, are also presented.