Abstract: Bifunctionality is a key feature of many industrial catalysts, supported metal clusters and particles in particular, and the development of such catalysts for the CO2 reduction reaction (CO2RR) to hydrocarbons and alcohols is gaining traction in light of recent advancements in the field. Carbon-supported Cu nanoparticles are suitable candidates for integration in the state-of-the-art reaction interfaces, and here, we propose, synthesize, and evaluate a bifunctional Ni–N-doped-C-supported Cu electrocatalyst, in which the support possesses active sites for selective CO2 conversion to CO and Cu nanoparticles catalyze either the direct CO2 or CO reduction to hydrocarbons. In this work, we introduce the scientific rationale behind the concept, its applicability, and the challenges with regard to the catalyst. From the practical aspect, the deposition of Cu nanoparticles onto carbon black and Ni–N–C supports via an ammonia-driven deposition precipitation method is reported and explored in more detail using X-ray diffraction, thermogravimetric analysis, and hydrogen temperature-programmed reduction. High-angle annular dark-field scanning transmission electron microscopy (HAADF-STEM) and energy-dispersive X-ray spectroscopy (EDXS) give further evidence of the presence of Cu-containing nanoparticles on the Ni–N–C supports while revealing an additional relationship between the nanoparticle’s composition and the electrode’s electrocatalytic performance. Compared to the benchmark carbon black-supported Cu catalysts, Ni–N–C-supported Cu delivers up to a 2-fold increase in the partial C2H4 current density at −1.05 VRHE (C1/C2 = 0.67) and a concomitant 10-fold increase of the CO partial current density. The enhanced ethylene production metrics, obtained by virtue of the higher intrinsic activity of the Ni–N–C support, point out toward a synergistic action between the two catalytic functionalities.

Abstract: Rapid development of experimental techniques has enabled real time studies of solid gas reactions at the level reaching the atomic scale. In the present paper, we focus on a combination of atomic resolution STEM/EDX, which visualizes the reaction zone, and near ambient pressure (NAP) XPS, which collects information for a surface layer of variable thickness under reaction conditions. We compare the behavior of two affined topological insulators, Bi2Te3 and Sb2Te3. We used a simple reaction with molecular oxygen occurring at 298 K, which is of practical importance to avoid material degradation. Despite certain limitations, a combination of in situ XPS and ex situ cross-sectional STEM/EDX allowed us to obtain a self-consistent picture of the solid gas reaction mechanism for oxidation of Sb2Te3 and Bi2Te3 crystals, which includes component redistribution between the oxide and the subsurface layer and Te segregation with formation of a thin ordered layer at the interface. The process is multistep in case of both compounds. At the very beginning of the oxidation process the reactivity is determined by the energy benefit of the corresponding element oxygen bond formation. Further in the oxidation process, the behavior of these two compounds becomes similar and features component redistribution between the oxide and the subsurface layer.

Abstract: The surface termination of MXenes greatly determines the electrochemical properties and ion kinetics on their surfaces. So far, hydroxyl-, oxygen-, and fluorine-terminated MXenes have been widely studied for energy storage applications. Recently, sulfur-functionalized MXene structures, which possess low diffusion barriers, have been proposed as candidate materials to enhance battery performance. We performed first-principles calculations on the structural, stability, electrochemical, and ion dynamic properties of Li-adsorbed sulfur-functionalized groups 3B, 4B, 5B, and 6B transition-metal (M)-based MXenes (i.e., M2CS2 with M = Sc, Ti, Zr, Hf, V, Nb, Ta, Cr, Mo, and W). We performed phonon calculations, which indicated that all of the above M2CS2 MXenes, except for Sc, are dynamically stable at T = 0 K. The ground-state structure of each M2CS2 monolayer depends on the type of M atom. For instance, while sulfur prefers to sit at the FCC site on Ti2CS2, it occupies the HCP site of Cr-based MXene. We determined the Li adsorption configurations at different concentrations using the cluster expansion method. The highest maximum open-circuit voltages were computed for the group 4B element (i.e., Ti, Zr, and Hf)-based M2CS2, which are larger than 2.1 V, while their average voltages are approximately 1 V. The maximum voltage for the group 6B element (i.e., Cr, Mo, W)-based M2CS2 is less than 1 V, and the average voltage is less than 0.71 V. We found that S functionalization is helpful for capacity improvements over the O-terminated MXenes. In this respect, the computed storage gravimetric capacity may reach up to 417.4 mAh/g for Ti2CS2 and 404.5 mAh/g for V2CS2. Ta-, Cr-, Mo-, and W-based M2CS2 MXenes show very low capacities, which are less than 100 mAh/g. The Li surface diffusion energy barriers for all of the considered MXenes are less than 0.22 eV, which is favorable for high charging and discharging rates. Finally, ab initio molecular dynamic simulations performed at 400 K and bond-length analysis with respect to Li concentration verify that selected promising systems are robust against thermally induced perturbations that may induce structural transformations or distortions and undesirable Li release.

Abstract: Water permeation between stacked layers of hBN sheets forming 2D nanochannels is investigated using large-scale ab initio-quality molecular dynamics simulations. A high-dimensional neural network potential trained on density-functional theory calculations is employed. We simulate water in van der Waals nanocapillaries and study the impact of nanometric confinement on the structure and dynamics of water using both equilibrium and nonequilibrium methods. At an interlayer distance of 10.2 A confinement induces a first-order phase transition resulting in a well-defined AA-stacked bilayer of hexagonal ice. In contrast, for h < 9 A, the 2D water monolayer consists of a mixture of different locally ordered patterns of squares, pentagons, and hexagons. We found a significant change in the transport properties of confined water, particularly for monolayer water where the water-solid friction coefficient decreases to half and the diffusion coefficient increases by a factor of 4 as compared to bulk water. Accordingly, the slip-velocity is found to increase under confinement and we found that the overall permeation is dominated by monolayer water adjacent to the hBN membranes at extreme confinements. We conclude that monolayer water in addition to bilayer ice has a major contribution to water transport through 2D nanochannels.

Abstract: The recent advantages of the fabrication of artificial nanochannels enabled new research on the molecular transport, permeance, and selectivity of various gases and molecules. However, the physisorption/chemisorption of the unwanted molecules (usually hydrocarbons) inside nanochannels results in the alteration of the functionality of the nanochannels. We investigated contamination due to hydrocarbon molecules, nanochannels made of graphene, hexagonal boron nitride, BC2N, and molybdenum disulfide using molecular dynamics simulations. We found that for a certain size of nanochannel (i.e., h = 0.7 nm), as a result of the anomalous hydrophilic nature of nanochannels made of graphene, the hydrocarbons are fully adsorbed in the nanochannel, giving rise to full uptake. An increasing temperature plays an important role in unclogging, while pressure does not have a significant role. The results of our pioneering work contribute to a better understanding and highlight the important factors in alleviating the contamination and unclogging of nanochannels, which are in good agreement with the results of recent experiments.

Abstract: A multiscale modeling and simulation approach, including first-principles calculations, ab initio molecular dynamics simulations, and a tight binding approach, is employed to study band flattening of the electronic band structure of oxidized monolayer graphene. The width offlat bands can be tuned by strain, the external electric field, and the density of functional groups and their distribution. A transition to a conducting state is found for monolayer graphene with impurities when it is subjected to an electric field of similar to 1.0 V/angstrom. Several parallel impurity-induced flat bands appear in the low-energy spectrum of monolayer graphene when the number of epoxy groups is changed. The width of the flat band decreases with an increase in tensile strain but is independent of the electric field strength. Here an alternative and easy route for obtaining band flattening in thermodynamically stable functionalized monolayer graphene is introduced. Our work discloses a new avenue for research on band flattening in monolayer graphene.

Abstract: Doping two single transition-metal (TM) atoms on a substrate host opens numerous possibilities for catalyst design. However, what if the substrate contains more than one vacancy site? Then, the combination of two TMs along with their distribution patterns becomes a design parameter potentially complementary to the substrate itself and the bimetal composition. In this study, we investigate ammonia synthesis under mild electrocatalytic conditions on a transition-metal-doped porous C24N24 catalyst using density functional theory (DFT). The TMs studied include Ti, Mn, and Cu in a 2:4 dopant ratio (Ti2Mn4@C24N24 and Ti2Cu4@N-24(24)). Our computations show that a single Ti atom in both catalysts exhibits the highest selectivity for N-2 fixation at ambient conditions. This work is a good theoretical model to establish the structure-activity relationship, and the knowledge earned from the metal-N-4 moieties may help studies of related nanomaterials, especially those with curved structures.

Abstract: Multimodal contrast agents, which take advantage of different imaging modalities, have emerged as an interesting approach to overcome the technical limitations of individual techniques. We developed hybrid nanoparticles comprising an iron oxide core and an outer gold spiky layer, stabilized by a biocompatible polymeric shell. The combined magnetic and optical properties of the different components provide the required functionalities for magnetic resonance imaging (MRI), surface-enhanced Raman scattering (SERS), and fluorescence imaging. The fabrication of such hybrid nanoprobes comprised the adsorption of small gold nanoparticles onto premade iron oxide cores, followed by controlled growth of spiky gold shells. The gold layer thickness and branching degree (tip sharpness) can be controlled by modifying both the density of Au nanoparticle seeds on the iron oxide cores and the subsequent nanostar growth conditions. We additionally demonstrated the performance of these hybrid multifunctional nanoparticles as multimodal contrast agents for correlative imaging of in vitro cell models and ex vivo tissues.

Abstract: Direct experimental measurement of thermal expansion coefficient without substrate effects is a challenging task for two-dimensional (2D) materials, and its accurate estimation with large-scale ab initio molecular dynamics is computationally very expensive. Machine learning-based interatomic potentials trained with ab initio data have been successfully used in molecular dynamics simulations to decrease the computational cost without compromising the accuracy. In this study, we investigated using Gaussian approximation potentials to reproduce the density functional theory-level accuracy for graphene within both lattice dynamical and molecular dynamical methods, and to extend their applicability to larger length and time scales. Two such potentials are considered, GAP17 and GAP20. GAP17, which was trained with pristine graphene structures, is found to give closer results to density functional theory calculations at different scales. Further vibrational and structural analyses verify that the same conclusions can be deduced with density functional theory level in terms of the reasoning of the thermal expansion behavior, and the negative thermal expansion behavior is associated with long-range out-of-plane phonon vibrations. Thus, it is argued that the enabled larger system sizes by machine learning potentials may even enhance the accuracy compared to small-size-limited ab initio molecular dynamics.

Abstract: Electrenes, an atomically thin form of layered electrides, are very recent members of the 2D materials family. In this work, we employed first principle calculations to determine stable, exfoliatable, and application-promising 2D electrene materials among possible M2X compounds, where M is a group II-A metal and X is a nonmetal element (C, N, P, As, and Sb). The promise of stable electrene compounds for battery applications is assessed via their exfoliation energy, adsorption properties, and migration energy barriers toward relevant Li, Na, K, and Ca atoms. Our calculations revealed five new stable electrene candidates in addition to previously known Ca2N and Sr2N. Among these seven dynamically stable electrenes, Ba2As, Ba2P, Ba2Sb, Ca2N, Sr2N, and Sr2P are found to be very promising for either K or Na ion batteries due to their extremely low migration energy barriers (5-16 meV), which roughly demonstrates 105 times higher mobility than graphene and two to four times higher mobility than other promising 2D materials such as MXene (Mo2C).

Abstract: The quantum electrostatic heterostructure (QEH) model allows for efficient computation of the dielectric screening properties of layered van der Waals (vdW)-bonded heterostructures in terms of the dielectric functions of the individual two-dimensional (2D) layers. Here, we extend the QEH model by including (1) contributions to the dielectric function from infrared active phonons in the 2D layers, (2) screening from homogeneous bulk substrates, and (3) intraband screening from free carriers in doped 2D semiconductor layers. We demonstrate the potential of the extended QEH model by calculating the dispersion of coupled phonons in multilayer stacks of hexagonal boron-nitride (hBN), the strong hybridization of plasmons and optical phonons in graphene/hBN heterostructures, the effect of substrate screening on the exciton series of monolayer MoS2, and the properties of hyperbolic plasmons in a doped phosphorene sheet. The new QEH code is distributed as a Python package with a simple command line interface and a comprehensive library of dielectric building blocks for the most common 2D materials, providing an efficient open platform for dielectric modeling of realistic vdW heterostructures.

Abstract: We evaluate the importance of quantum effects on the sticking of H2 on Al(110) for conditions that are close to those of molecular beam experiments that have been done on this system. Calculations with the quasi-classical trajectory (QCT) method and with quantum dynamics (QD) are performed using a model in which only motion in the six molecular degrees of freedom is allowed. The potential energy surface used has a minimum barrier height close to the value recently obtained with the quantum Monte Carlo method. Monte Carlo averaging over the initial rovibrational states allowed the QD calculations to be done with an order of magnitude smaller computational expense. The sticking probability curve computed with QD is shifted to lower energies relative to the QCT curve by 0.21 to 0.05 kcal/mol, with the highest shift obtained for the lowest incidence energy. Quantum effects are therefore expected to play a small role in calculations that would evaluate the accuracy of electronic structure methods for determining the minimum barrier height to dissociative chemisorption for H2 + Al(110) on the basis of the standard procedure for comparing results of theory with molecular beam experiments.