RESUMO
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.
RESUMO
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.
RESUMO
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.
RESUMO
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.
RESUMO
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.
RESUMO
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.
RESUMO
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.
RESUMO
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.
RESUMO
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.
RESUMO
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).
RESUMO
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.
RESUMO
Octahedral tilts are the most ubiquitous distortions in perovskite-related structures that can dramatically influence ferroelectric, magnetic, and electronic properties; yet the paradigm of tilt epitaxy in thin films is barely explored. Non-destructively characterizing such epitaxy in three-dimensions for low symmetry complex tilt systems composed of light anions is a formidable challenge. Here we demonstrate that the interfacial tilt epitaxy can transform ultrathin calcium titanate, a non-polar earth-abundant mineral, into high-temperature polar oxides that last above 900 K. The comprehensive picture of octahedral tilts and polar distortions is revealed by reconstructing the three-dimensional electron density maps across film-substrate interfaces with atomic resolution using coherent Bragg rod analysis. The results are complemented with aberration-corrected transmission electron microscopy, film superstructure reflections, and are in excellent agreement with density functional theory. The study could serve as a broader template for non-destructive, three-dimensional atomic resolution probing of complex low symmetry functional interfaces.
RESUMO
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.
RESUMO
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.
RESUMO
The friction of adjacent Tin+1Cn (n = 1, 2, and 3) MXene layers is investigated using density functional theory (DFT) calculations and classical molecular dynamics simulations with ReaxFF potentials. The calculations reveal the sliding pathways in all three MXene systems with low energy barriers. The friction coefficients for interlayer sliding are evaluated using static calculations. Both DFT and ReaxFF methods predict friction coefficients between 0.24 and 0.27 for normal loads less than 1.2 GPa. The effect of titanium (Ti) vacancies in sublayers and terminal oxygen (O) vacancies at surfaces on the interlayer friction is further investigated using the ReaxFF potential. These defects are found to increase the friction coefficients by increasing surface roughness and creating additional attractive forces between adjacent layers. However, these defective MXenes still maintain friction coefficients below 0.31. We also consider functionalized Ti3C2 MXene terminated with -OH and -OCH3 and find that compared to the -O-terminated surface both groups further reduce the interlayer friction coefficient to 0.10-0.14.
RESUMO
Unraveling the physical and chemical properties of graphene-metal contacts is a key step toward the development of graphitic electronic nanodevices. Although many studies have revealed the way that various metals interact with graphene, few have described the structure and behavior of large pieces of graphene-metal nanostructures under different conditions. Here, we present the first classical molecular dynamics study of graphene-titanium (G-Ti) structures, with and without substrates. Physical and chemical properties of equilibrium structures of G-Ti interfaces with different amounts of titanium coverage are investigated. Adhesion of Ti films on graphene is shown to be enhanced by the vacancies in graphene or the electrostatic influence of substrates. The dynamics of pristine G-Ti structures at different temperatures on planar and nonplanar substrates are investigated, and the results show that G-Ti interfaces are thermally stable, that is, not prone to any reaction toward the formation of titanium carbide.
RESUMO
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.
RESUMO
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).
RESUMO
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.