Heterostructures coupling ultrathin metal carbides and chalcogenides.

Non-layered transition metal carbides (TMCs) and layered transition metal dichalcogenides (TMDs) are two well-studied material families that have individually received considerable attention over the past century. In recent years, with the shift towards two-dimensional materials and heterostructures, a field has emerged that is focused on the structure and properties of TMC/TMD heterostructures, which through chemical conversion exhibit diverse types of heterostructure configuration that host coupled 2D-3D interfaces, giving rise to exotic properties. In this Review, we highlight experimental and computational efforts to understand the routes to fabricate TMC/TMD heterostructures. Furthermore, we showcase how controlling these heterostructures can lead to emergent electronic transport, optical properties and improved catalytic properties.

Predicted Separation of Acid Gases from Gas Mixtures by Functionalized Porous Aromatic Frameworks.

The selective adsorption of target acid gas molecules from binary gas mixtures by porous aromatic frameworks (PAFs) with two identical functional groups per aromatic ring (PAF-R2) was computationally investigated using grand canonical Monte Carlo simulations. PAF-R2 adsorption was considered for three binary mixtures of small molecular concentrations of acid gas and abundant nitrogen gas (CO2/N2, SO2/N2, and H2S/N2). The results indicate that additional functional groups enhance acid gas loadings and selectivity, compared with pristine PAF and single-functionalized PAFs. Low pressures yield linearly increasing gas loadings and constant selectivity, while high pressures yield much higher adsorption and selectivity. In particular, SO2 loading and selectivity under high pressures are heavily influenced by the PAF's maximum adsorption limit, which can be linked back to the functional groups and their configuration. In summary, PAF-(3,5)-(COOH)2 (nomenclature of PAFs is provided in the Appendix in the Supporting Information) and many other PAF with the same two electron-withdrawing groups are predicted to have great acid gas adsorption and selectivity from gas mixtures, while PAF-(3,5)-(OH)2 (one of PAFs with two identical electron-donating groups) is predicted to have good adsorption and selectivity, especially under elevated pressures. The results of this work can provide insights into various types of PAFs with great selective adsorption ability and their corresponding conditions. The simulation procedures and results may inspire the exploration and screening of other types of PAFs or porous materials, for acid gas absorption.

Superconductivity enhancement in phase-engineered molybdenum carbide/disulfide vertical heterostructures.

Stacking layers of atomically thin transition-metal carbides and two-dimensional (2D) semiconducting transition-metal dichalcogenides, could lead to nontrivial superconductivity and other unprecedented phenomena yet to be studied. In this work, superconducting α-phase thin molybdenum carbide flakes were first synthesized, and a subsequent sulfurization treatment induced the formation of vertical heterolayer systems consisting of different phases of molybdenum carbide-ranging from α to γ' and γ phases-in conjunction with molybdenum sulfide layers. These transition-metal carbide/disulfide heterostructures exhibited critical superconducting temperatures as high as 6 K, higher than that of the starting single-phased α-Mo2C (4 K). We analyzed possible interface configurations to explain the observed moiré patterns resulting from the vertical heterostacks. Our density-functional theory (DFT) calculations indicate that epitaxial strain and moiré patterns lead to a higher interfacial density of states, which favors superconductivity. Such engineered heterostructures might allow the coupling of superconductivity to the topologically nontrivial surface states featured by transition-metal carbide phases composing these heterostructures potentially leading to unconventional superconductivity. Moreover, we envisage that our approach could also be generalized to other metal carbide and nitride systems that could exhibit high-temperature superconductivity.

Optimized utilization of COMB3 reactive potentials in LAMMPS.

An investigation to optimize the application of the third-generation charge optimized many-body (COMB3) interatomic potential and associated input parameters was carried out through the study of solid-liquid interactions in classical molecular dynamics simulations. The rates of these molecular interactions are understood through the wetting rates of water nano-droplets on a bare copper (111) surface. Implementing the Langevin thermostat, the influence of simulation time step, the number of atoms in the system, the frequency at which charge equilibration is performed, and the temperature relaxation rate are all examined. The results indicate that time steps of 0.4 fs are possible when using longer relaxation times for the system temperature, which is almost double the typical time step used for reactive potentials. The use of the charge equilibration allows for a fewer atomic layers to be used in the Cu slab. In addition, charge equilibrium schemes do not need to be performed every time step to ensure accurate charge transfer. Interestingly, the rate of wetting for the nanodroplets is dominantly dependent on the temperature relaxation time, which is predicted to significantly change the viscosity of the water droplets. This work provides a pathway for optimizing simulations using the COMB3 reactive interatomic potential.

Effects of surface charge and cluster size on the electrochemical dissolution of platinum nanoparticles using COMB3 and continuum electrolyte models.

We study the site-dependent dissolution of platinum nanoparticles under electrochemical conditions to assess their thermodynamic stability as a function of shape and size using empirical molecular dynamics and electronic-structure models. The third-generation charge optimized many-body potential is employed to determine the validity of uniform spherical representations of the nanoparticles in predicting dissolution potentials (the Kelvin model). To understand the early stages of catalyst dissolution, implicit solvation techniques based on the self-consistent continuum solvation method are applied. It is demonstrated that interfacial charge and polarization can shift the dissolution energies by amounts on the order of 0.74 eV depending on the surface site and nanoparticle shape, leading to the unexpected preferential removal of platinum cations from highly coordinated sites in some cases.

Multi-Step Topochemical Pathway to Metastable Mo2AlB2 and Related Two-Dimensional Nanosheet Heterostructures.

The rational synthesis of metastable inorganic solids, which is a grand challenge in solid-state chemistry, requires the development of kinetically controlled reaction pathways. Topotactic strategies can achieve this goal by chemically modifying reactive components of a parent structure under mild conditions to produce a closely related analogue that has otherwise inaccessible structures and/or compositions. Refractory materials, such as transition metal borides, are difficult to structurally manipulate at low temperatures because they generally are chemically inert and held together by strong covalent bonds. Here, we report a multistep low-temperature topotactic pathway to bulk-scale Mo2AlB2, which is a metastable phase that has been predicted to be the precursor needed to access a synthetically elusive family of 2-D metal boride (MBene) nanosheets. Room-temperature chemical deintercalation of Al from the stable compound MoAlB (synthesized as a bulk powder at 1400 °C) formed highly strained and destabilized MoAl1-xB, which was size-selectively precipitated to isolate the most reactive submicron grains and then annealed at 600 °C to deintercalate additional Al and crystallize Mo2AlB2. Further heating resulted in topotactic decomposition into bulk-scale Mo2AlB2-AlOx nanolaminates that contain Mo2AlB2 nanosheets with thickness of 1-3 nm interleaved by 1-3 nm of amorphous aluminum oxide. The combination of chemical destabilization, size-selective precipitation, and low-temperature annealing provides a potentially generalizable kinetic pathway to metastable variants of refractory compounds, including bulk Mo2AlB2 and Mo2AlB2-AlOx nanosheet heterostructures, and opens the door to other previously elusive 2-D materials such as 2-D MoB (MBene).

The structure of graphene on graphene/C60/Cu interfaces: a molecular dynamics study.

Two experimental studies reported the spontaneous formation of amorphous and crystalline structures of C60 molecules intercalated between graphene and a surface. The findings observed included interesting phenomena ranging from reaction between fullerene C60s ('C60s' stands for plural of C60) under graphene to graphene sheets sagging between C60s and control of strain in these sheets. Motivated by this work, we performed fully atomistic reactive molecular dynamics simulations to investigate the formation and thermal stability of graphene sheet wrinkles as well as graphene attachment to and detachment from a surface when the sheet is laid over a previously distributed array of C60 molecules on a copper surface at different temperatures. As graphene compresses the C60s against the surface, and graphene attachment to the surface in between C60s depends on the height of the wrinkles in the graphene sheet, configurations with both frozen and non-frozen fullerenes were investigated in the simulations in order to examine the experimental result of stable, sagged graphene sheets when the distance between C60s is about 4 nm and the height of the wrinkles in the sheet is about 0.8 nm. Below a distance of 4 nm between fullerenes, the graphene is predicted to become locally suspended and less strained. The simulations predict that this happens when the fullerenes can deform under the compressive action of the graphene sheet. If the fullerenes are kept frozen, spontaneous 'blanketing' of graphene is predicted only when the distance between neighbouring C60s is equal to or great than about 7 nm. These predictions agree with a mechanical model relating the rigidity of a graphene sheet to the energy of graphene-surface adhesion. This work further reveals the structure of intercalated molecules and the role of stability and sheet wrinkling on the preferred configuration of graphene. This study thus might assist in the development of two-dimensional confined nanoreactors for chemical reactions.

Nanoscale Structure and Dynamics of Water on Pt and Cu Surfaces from MD Simulations.

The interaction of liquid water with Pt(111) is investigated with classical molecular dynamics (MD) simulations, where the forces are determined using the third-generation charge optimized many-body (COMB3) interatomic potential. In cases of sub-monolayer water coverage, the parameterized empirical potential predicts experimentally observed and energetically favorable √37 and √39 reconstructed water structures with "575757" di-interstitial defects. At both sub-monolayer and multilayer water coverages, the structure of the first wetting layer of liquid water on Pt(111) exhibits a characteristic distribution where the molecules form two distinct buckled layers as a result of the interplay between water-metal adsorption and water-water hydrogen bonds. The dynamic spreading rate of water nanodroplets on large Pt surfaces (>200 nm2) characterized by molecular kinetic spreading theory is an order of magnitude slower than the molecular kinetic rate of the same droplet on close-packed Cu surfaces due to variation in molecular distributions at the water-metal interface. These nanoscale MD simulation predictions using the COMB3 interatomic potential demonstrate the capability of capturing both many-body interactions between H2O and Pt or Cu and hydrogen bonding in liquid water.

Applied Potentials in Variable-Charge Reactive Force Fields for Electrochemical Systems.

An atomic description of water dynamics and electrochemical properties at electrode-electrolyte interfaces is presented using molecular dynamics with the third generation of the charge-optimized many-body (COMB3) potential framework. Externally applied potentials in electrochemical applications were simulated by offsetting electronegativity on electrode atoms. This approach is incorporated into the variable charge scheme within COMB3 and is used to investigate electrochemical systems consisting of two Cu electrodes and a water electrolyte with varying concentrations of hydroxyls (OH-) and protons (H+). The interactions between the electronegativity offset method and the charge equilibration method in a variable charge scheme are analyzed. In addition, a charge equilibration method for electrochemical applications is proposed, where the externally applied potentials are treated by the electronegativity offset on the electrodes thus enforcing charge neutrality on the electrolyte. This method is able to qualitatively capture the relevant electrochemistry and predict consistently correct voltages with precalibration.

Two-Dimensional Intrinsic Half-Metals With Large Spin Gaps.

Through a systematic search of all layered bulk compounds combined with density functional calculations employing hybrid exchange-correlation functionals, we predict a family of three magnetic two-dimensional (2D) materials with half-metallic band structures. The 2D materials, FeCl2, FeBr2, and FeI2, are all sufficiently stable to be exfoliated from bulk layered compounds. The Fe2+ ions in these materials are in a high-spin octahedral d6 configuration leading to a large magnetic moment of 4 µB. Calculations of the magnetic anisotropy show an easy-plane for the magnetic moment. A classical XY model with nearest neighbor coupling estimates critical temperatures, Tc, for the Berezinskii-Kosterlitz-Thouless transition ranging from 122 K for FeI2 to 210 K for FeBr2. The quantum confinement of these 2D materials results in unusually large spin gaps, ranging from 4.0 eV for FeI2 to 6.4 eV for FeCl2, which should defend against spin current leakage even at small device length scales. Their purely spin-polarized currents and dispersive interlayer interactions should make these materials useful for 2D spin valves and other spintronic applications.

Topology-Scaling Identification of Layered Solids and Stable Exfoliated 2D Materials.

The Materials Project crystal structure database has been searched for materials possessing layered motifs in their crystal structures using a topology-scaling algorithm. The algorithm identifies and measures the sizes of bonded atomic clusters in a structure's unit cell, and determines their scaling with cell size. The search yielded 826 stable layered materials that are considered as candidates for the formation of two-dimensional monolayers via exfoliation. Density-functional theory was used to calculate the exfoliation energy of each material and 680 monolayers emerge with exfoliation energies below those of already-existent two-dimensional materials. The crystal structures of these two-dimensional materials provide templates for future theoretical searches of stable two-dimensional materials. The optimized structures and other calculated data for all 826 monolayers are provided at our database (https://materialsweb.org).

Effect of Surface Chemistry on Water Interaction with Cu(111).

The interfacial dynamics of water in contact with bare, oxidized, and hydroxylated copper surfaces are examined using classical molecular dynamics (MD) simulations. A third-generation charge-optimized many-body (COMB3) potential is used in the MD simulations to investigate the adsorption of water molecules on Cu(111), and the results are compared to the findings of density functional theory (DFT) calculations. The adsorption energies and structures predicted by COMB3 are generally consistent with those determined with DFT. The COMB3 potential is then used to investigate the wetting behavior of water nanodroplets on Cu(111) at 20, 130, and 300 K. At room temperature, the simulations predict that the spreading rate of the base radius, R0, of a water droplet with a diameter of about 1.5 nm exhibits a spreading rate of R0 ≈ t(0.16) and a final base radius of 3.5 nm. At 20 and 130 K, water droplets are predicted to retain their structure after adsorption on Cu(111) and to undergo minimal spreading in agreement with scanning tunneling microscopy data. When the same water droplet encounters a reconstructed, oxidized Cu(111) surface, the classical MD simulations predict wetting with a spreading rate of R ≈ t(0.14) and a final base radius of 3.0 nm. Similarly, our MD simulations predict a spreading rate of R ≈ t(0.14) and a final base radius of 2.5 nm when water encounters OH-covered Cu(111). These results indicate that oxidation and hydroxylation cause a reduction in the degree of spreading and final base radius that is directly associated with a decreased spreading rate for water nanodroplets on copper.

Chemical Environment and Structural Variations in High Entropy Oxide Thin Film Probed with Electron Microscopy.

We employ analytical transmission electron microscopy (TEM) to correlate the structural and chemical environment variations within a stacked epitaxial thin film of the high entropy oxide (HEO) Mg0.2Co0.2Ni0.2Cu0.2Zn0.2O (J14), with two layers grown at different substrate temperatures (500 and 200 °C) using pulsed laser deposition (PLD). Electron diffraction and atomically resolved STEM imaging reveal the difference in out-of-plane lattice parameters in the stacked thin film, which is further quantified on a larger scale using four-dimensional STEM (4D-STEM). In the layer deposited at a lower temperature, electron energy loss spectroscopy (EELS) mapping indicates drastic changes in the oxidation states and bonding environment for Co ions, and energy-dispersive X-ray spectroscopy (EDX) mapping detects more significant cation deficiency. Ab initio density functional theory (DFT) calculations validate that vacancies on the cation sublattice of J14 result in significant electronic and structural changes. The experimental and computational analyses indicate that low temperatures during film growth result in cation deficiency, an altered chemical environment, and reduced lattice parameters while maintaining a single phase. Our results demonstrate that the complex correlation of configurational entropy, kinetics, and thermodynamics can be utilized for accessing a range of metastable configurations in HEO materials without altering cation proportions, enabling further engineering of functional properties of HEO materials.

Quasi-Van der Waals Epitaxial Growth of γ'-GaSe Nanometer-Thick Films on GaAs(111)B Substrates.

GaSe is an important member of the post-transition-metal chalcogenide family and is an emerging two-dimensional (2D) semiconductor material. Because it is a van der Waals material, it can be fabricated into atomic-scale ultrathin films, making it suitable for the preparation of compact, heterostructure devices. In addition, GaSe possesses unusual optical and electronic properties, such as a shift from an indirect-bandgap single-layer film to a direct-bandgap bulk material, rare intrinsic p-type conduction, and nonlinear optical behaviors. These properties make GaSe an appealing candidate for the fabrication of field-effect transistors, photodetectors, and photovoltaics. However, the wafer-scale production of pure GaSe single-crystal thin films remains challenging. This study develops an approach for the direct growth of nanometer-thick GaSe films on GaAs substrates by using molecular beam epitaxy. It yields smooth thin GaSe films with a rare γ'-polymorph. We analyze the formation mechanism of γ'-GaSe using density-functional theory and speculate that it is stabilized by Ga vacancies since the formation enthalpy of γ'-GaSe tends to become lower than that of other polymorphs when the Ga vacancy concentration increases. Finally, we investigate the growth conditions of GaSe, providing valuable insights for exploring 2D/three-dimensional (3D) quasi-van der Waals epitaxial growth.

Variable charge reactive potential for hydrocarbons to simulate organic-copper interactions.

A variable charge reactive empirical potential for carbon-based materials, hydrocarbons, organometallics, and their interfaces is developed within the framework of charge optimized many-body (COMB) potentials. The resulting potential contains improved expressions for the bond order and self-energy, which gives a flexible, robust, and integrated treatment of different bond types in multicomponent and multifunctional systems. It furthermore captures the dissociation and formation of the chemical bonds and appropriately and dynamically determines the associated charge transfer, thus providing a powerful method to simulate the complex chemistry of many-atom systems in changing environments. The resulting COMB potential is used in a classical molecular dynamics simulation of the room temperature, low energy deposition of ethyl radicals on the Cu (111) surface (a system with ∼5000 atoms) to demonstrate its capabilities at describing organic-metal interactions in a dynamically changing environment.

Effect of fluorocarbon molecules confined between sliding self-mated PTFE surfaces.

We report on the frictional response and atomic process that occur when molecular fluorocarbon molecules of varying lengths are sheared between two polytetrafluoroethylene (PTFE) surfaces. The thicknesses of the molecular layers are also varied. The approach is classical molecular dynamics simulations using a reactive bond-order potential parametrized for fluorocarbons. The results indicate that the presence of the molecules has a significant impact on the measured friction and wear of the surfaces, and that this impact depends on the nature of the fluorocarbon molecules and the thickness of the molecular film. The molecular mechanisms responsible for these differences are presented.

Single-Step Direct Laser Writing of Multimetal Oxygen Evolution Catalysts from Liquid Precursors.

We investigate a laser direct-write method to synthesize and deposit metastable, mixed transition metal oxides and evaluate their performance as oxygen evolution reaction catalysts. This laser processing method enabled the rapid synthesis of diverse heterogeneous alloy and oxide catalysts directly from cost-effective solution precursors, including catalysts with a high density of nanocrystalline metal alloy inclusions within an amorphous oxide matrix. The nanoscale heterogeneous structures of the synthesized catalysts were consistent with reactive force-field Monte Carlo calculations. By evaluating the impact of varying transition metal oxide composition ratios, we created a stable Fe0.63Co0.19Ni0.18Ox/C catalyst with a Tafel slope of 38.23 mV dec-1 and overpotential of 247 mV, a performance similar to that of IrO2. Synthesized Fe0.63Co0.19Ni0.18Ox/C and Fe0.14Co0.46Ni0.40Ox/C catalysts were experimentally compared in terms of catalytic performance and structural characteristics to determine that higher iron content and a less crystalline structure in the secondary matrix decrease the charge transfer resistance and thus is beneficial for electrocatalytic activity. This conclusion is supported by density-functional theory calculations showing distorted active sites in ternary metal catalysts are key for lowering overpotentials for the oxygen evolution reaction.

The effect of normal load on polytetrafluoroethylene tribology.

The tribological behavior of oriented poly(tetrafluoroethylene) (PTFE) sliding surfaces is examined as a function of sliding direction and applied normal load in classical molecular dynamics (MD) simulations. The forces are calculated with the second-generation reactive empirical bond-order potential for short-range interactions, and with a Lennard-Jones potential for long-range interactions. The range of applied normal loads considered is 5-30 nN. The displacement of interfacial atoms from their initial positions during sliding is found to vary by a factor of seven, depending on the relative orientation of the sliding chains. However, within each sliding configuration the magnitude of the interfacial atomic displacements exhibits little dependence on load over the range considered. The predicted friction coefficients are also found to vary with chain orientation and are in excellent quantitative agreement with experimental measurements.

Carbon doping of WS2 monolayers: Bandgap reduction and p-type doping transport.

Chemical doping constitutes an effective route to alter the electronic, chemical, and optical properties of two-dimensional transition metal dichalcogenides (2D-TMDs). We used a plasma-assisted method to introduce carbon-hydrogen (CH) units into WS2 monolayers. We found CH-groups to be the most stable dopant to introduce carbon into WS2, which led to a reduction of the optical bandgap from 1.98 to 1.83 eV, as revealed by photoluminescence spectroscopy. Aberration corrected high-resolution scanning transmission electron microscopy (AC-HRSTEM) observations in conjunction with first-principle calculations confirm that CH-groups incorporate into S vacancies within WS2. According to our electronic transport measurements, undoped WS2 exhibits a unipolar n-type conduction. Nevertheless, the CH-WS2 monolayers show the emergence of a p-branch and gradually become entirely p-type, as the carbon doping level increases. Therefore, CH-groups embedded into the WS2 lattice tailor its electronic and optical characteristics. This route could be used to dope other 2D-TMDs for more efficient electronic devices.