Bromination of 2D materials.

The adsorption, reaction and thermal stability of bromine on Rh(111)-supported hexagonal boron nitride (h-BN) and graphene were investigated. Synchrotron radiation-based high-resolution x-ray photoelectron spectroscopy (XPS) and temperature-programmed XPS allowed us to follow the adsorption process and the thermal evolutionin situon the molecular scale. Onh-BN/Rh(111), bromine adsorbs exclusively in the pores of the nanomesh while we observe no such selectivity for graphene/Rh(111). Upon heating, bromine undergoes an on-surface reaction onh-BN to form polybromides (170-240 K), which subsequently decompose to bromide (240-640 K). The high thermal stability of Br/h-BN/Rh(111) suggests strong/covalent bonding. Bromine on graphene/Rh(111), on the other hand, reveals no distinct reactivity except for intercalation of small amounts of bromine underneath the 2D layer at high temperatures. In both cases, adsorption is reversible upon heating. Our experiments are supported by a comprehensive theoretical study. DFT calculations were used to describe the nature of theh-BN nanomesh and the graphene moiré in detail and to study the adsorption energetics and substrate interaction of bromine. In addition, the adsorption of bromine onh-BN/Rh(111) was simulated by molecular dynamics using a machine-learning force field.

Caracal: A Versatile Ring Polymer Molecular Dynamics Simulation Package.

A new open-source program package named Caracal covering simulations of molecular systems with ring polymer molecular dynamics (RPMD) is presented. It combines a powerful RPMD implementation including chemical reaction rate calculations and biased periodic and nonperiodic samplings with a collection of easy to set up potential energy surface (PES) methodologies, thus delivering an all-inclusive approach. Most implemented PESs are based on the QMDFF and EVB-QMDFF methods. Where the quantum mechanically derived force field (QMDFF) can be set up for an arbitrary molecular system in a black-box fashion, the empirical valence bond (EVB)-QMDFF connects two QMDFFs and is able to represent the PES of a chemical reaction. With our previously published flavors of this composite method, PESs for almost arbitrary gas phase thermal ground state reactions can be set up. Given an optimized reaction path, the mechanism of the reaction can be classified and RPMD rate constants can be obtained via umbrella sampling and recrossing calculations on an EVB-QMDFF PES. Further, QMDFFs can be polymerized for the description of liquid systems. In this paper, the internal structure as well as the handling philosophy of Caracal are outlined. Further, examples of the different possible kinds of calculations are given.

Bromine Adsorption and Thermal Stability on Rh(111): A Combined XPS, LEED and DFT Study.

This study addresses a fundamental question in surface science: the adsorption of halogens on metal surfaces. Using synchrotron radiation-based high-resolution X-ray photoelectron spectroscopy (XPS), temperature-programmed XPS, low-energy electron diffraction (LEED) and density functional theory (DFT) calculations, we investigated the adsorption and thermal stability of bromine on Rh(111) in detail. The adsorption of elemental bromine on Rh(111) at 170 K was followed in situ by XPS in the Br 3d region, revealing two individual, coverage-dependent species, which we assign to fcc hollow- and bridge-bound atomic bromine. In addition, we find a significant shift in binding energy upon increasing coverage due to adsorbate-adsorbate interactions. Subsequent heating shows a high thermal stability of bromine on Rh(111) up to above 1000 K, indicating strong covalent bonding. To complement the XPS data, LEED was used to study the long-range order of bromine on Rh(111): we observe a (√3×√3)R30° structure for low coverages (≤0.33 ML) and a star-shaped compression structure for higher coverages (0.33-0.43 ML). Combining LEED and DFT calculations, we were able to visualize bromine adsorption on Rh(111) in real space for varying coverages.

Structure Formation in an Ionic Liquid Wetting Layer: A Combined STM, IRAS, DFT and MD Study of [C2 C1 Im][OTf] on Au(111).

In a solid catalyst with ionic liquid layer (SCILL), ionic liquid (IL) coatings are used to improve the selectivity of noble metal catalysts. To understand the origins of this selectivity control, we performed model studies by surface science methods in ultrahigh vacuum (UHV). We investigated the growth and thermal stability of ultrathin IL films by infrared reflection absorption spectroscopy (IRAS). We combined these experiments with scanning tunneling microscopy (STM) to obtain information on the orientation of the ions, the interactions with the surface, the intermolecular interactions, and the structure formation. Additionally, we performed DFT calculations and molecular dynamics (MD) simulations to interpret the experimental data. We studied the IL 1-ethyl-3-methylimidazolium trifluoromethanesulfonate [C2 C1 Im][OTf] on Au(111) surfaces. We observe a weakly bound multilayer of [C2 C1 Im][OTf], which is stable up to 390 K, while the monolayer desorbs at ∼450 K. [C2 C1 Im][OTf] preferentially adsorbs at the step edges and elbows of the herringbone reconstruction of Au(111). The anion adsorbs via the SO3 group with the molecular axis perpendicular to the surface. At low coverage, the [C2 C1 Im][OTf] crystallizes in a glass-like 2D phase with short-range order. At higher coverage, we observe a phase transition to a 6-membered ring structure with long-range order.

Comparison of Implicit and Explicit Solvent Approaches in Ab Initio Evaluation of Thermochemistry in Solution: Application in Studying Boron Isotope Fractionation in Water.

Evaluation of thermochemistry in solution plays a key role in numerous fields. For this task, the solvent effects are commonly included in theoretical computations based on either implicit or explicit solvent approaches. In the present study, we evaluate and compare the performance of some of the most widely applied methods based on these two approaches. For studying the solvent effect on thermochemistry with an explicit solvent, we demonstrate that partial normal mode analysis with frozen geometry of solvent molecules for multiple solute-solvent configurations can yield quite accurate and reliable results for a drastically reduced computational cost. As a case study, we consider the evaluation of the equilibrium constant for the boron isotope exchange between boric acid and borate (k3-4) in pure and saline water which is of high geochemical importance. Employing three different rigorous and high-precision theoretical approaches, we provide a reliable estimation of k3-4 which is a value within 1.028 to 1.030 for both pure and saline water which is in excellent agreement with experimental data.

A new class of reaction path based potential energy surfaces enabling accurate black box chemical rate constant calculations.

A new method for constructing a full-dimensional potential energy surface representation in black-box fashion for an arbitrary reaction is presented. With limited knowledge of the system and with a limited number of reference-level data points, it is possible to calculate reaction rate constants with high quality. Building on our recently published application of Grimme's quantum-mechanically derived force field (QMDFF) and its empirical valence bond extension EVB-QMDFF to rate constant calculations, an improved EVB coupling method with local corrections was developed in order to avoid spurious problems for certain systems and hence to achieve an even wider range of applicability. A given reaction path (RP) is modeled as a parametric curve via cubic spline interpolation; regions offside this path are then extrapolated with quadratic Taylor series, and regions around the transition state are corrected by introduction of direct reference interpolation; the method is named transition region corrected RP-EVB-QMDFF (TREQ). To verify the quality of TREQ, six reactions were chosen for which full-dimensional analytical potential surfaces are available in the literature. Chemical reaction rates were calculated with ring polymer molecular dynamics on the reference surfaces as well as on the TREQ surfaces resulting in excellent agreement.

Cheap but accurate calculation of chemical reaction rate constants from ab initio data, via system-specific, black-box force fields.

Building on the recently published quantum-mechanically derived force field (QMDFF) and its empirical valence bond extension, EVB-QMDFF, it is now possible to generate a reliable potential energy surface for any given elementary reaction step in an essentially black box manner. This requires a limited and pre-defined set of reference data near the reaction path and generates an accurate approximation of the reference potential energy surface, on and off the reaction path. This intermediate representation can be used to generate reaction rate data, with far better accuracy and reliability than with traditional approaches based on transition state theory (TST) or variational extensions thereof (VTST), even if those include sophisticated tunneling corrections. However, the additional expense at the reference level remains very modest. We demonstrate all this for three arbitrarily chosen example reactions.