Ab initio molecular dynamics of atomic-scale surface reactions: insights into metal organic chemical vapor deposition of AlN on graphene

Study on the growth mechanism of monolayer and few-layer hexagonal boron nitride films on copper foil

During the process of synthesizing h-BN on Cu foil via chemical vapor deposition (CVD), low-pressure CVD (LPCVD) typically synthesizes monolayer h-BN films, whereas atmospheric pressure CVD (APCVD) yields few-layer h-BN films. Herein, a growth mechanism for monolayer and few-layer h-BN on Cu foil is proposed using first-principles calculations: Cu(111) passivated h-BN hinders the diffusion of B and N atoms at the subsurface of Cu(111), preventing sufficient transportation of B and N atoms to the existing h-BN layer, thereby leading to the formation of monolayer h-BN in LPCVD. For APCVD, the edges of h-BN are passivated by H, which decreases the barrier energy for the diffusion of B and N atoms on the Cu(111) subsurface, and B and N atoms can easily migrate from the subsurface of Cu(111) to its surface, resulting in the nucleation of h-BN between the existing h-BN and Cu(111), and leading to the formation of few-layer h-BN films. This work provides a theoretical basis at the atomic scale for further understanding the growth of monolayer and few-layer h-BN films on Cu foil.

PHOTH-graphene: a new 2D carbon allotrope with low barriers for Li-ion mobility

The continued interest in 2D carbon allotropes stems from their unique structural and electronic characteristics, which are crucial for diverse applications. This work theoretically introduces PHOTH-Graphene (PHOTH-G), a novel 2D planar carbon allotrope formed by 4-5-6-7-8 carbon rings. PHOTH-G emerges as a narrow band gap semiconducting material with low formation energy, demonstrating good stability under thermal and mechanical conditions. This material has slight mechanical anisotropy with Young modulus and Poisson ratios varying between 7.08-167.8 GPa and 0.21-0.96. PHOTH-G presents optical activity restricted to the visible range. Li atoms adsorbed on its surface have a migration barrier averaging 0.38 eV.

Proposing TODD-graphene as a novel porous 2D carbon allotrope designed for superior lithium-ion battery efficiency

The category of 2D carbon allotropes has gained considerable interest due to its outstanding optoelectronic and mechanical characteristics, which are crucial for various device applications, including energy storage. This study uses density functional theory calculations, ab initio molecular dynamics (AIMD), and classical reactive molecular dynamics (MD) simulations to introduce TODD-Graphene, an innovative 2D planar carbon allotrope with a distinctive porous arrangement comprising 3-8-10-12 carbon rings. TODD-G exhibits intrinsic metallic properties with a low formation energy and stability in thermal and mechanical behavior. Calculations indicate a substantial theoretical capacity for adsorbing Li atoms, revealing a low average diffusion barrier of 0.83 eV. The metallic framework boasts excellent conductivity and positioning TODD-G as an active layer for superior lithium-ion battery efficiency. Charge carrier mobility calculations for electrons and holes in TODD-G surpass those of graphene. Classical reactive MD simulation results affirm its structural integrity, maintaining stability without bond reconstructions at 2200 K.

Remote Epitaxy: Fundamentals, Challenges, and Opportunities

Advanced heterogeneous integration technologies are pivotal for next-generation electronics. Single-crystalline materials are one of the key building blocks for heterogeneous integration, although it is challenging to produce and integrate these materials. Remote epitaxy is recently introduced as a solution for growing single-crystalline thin films that can be exfoliated from host wafers and then transferred onto foreign platforms. This technology has quickly gained attention, as it can be applied to a wide variety of materials and can realize new functionalities and novel application platforms. Nevertheless, remote epitaxy is a delicate process, and thus, successful execution of remote epitaxy is often challenging. Here, we elucidate the mechanisms of remote epitaxy, summarize recent breakthroughs, and discuss the challenges and solutions in the remote epitaxy of various material systems. We also provide a vision for the future of remote epitaxy for studying fundamental materials science, as well as for functional applications.

Toward steering the motion of surface rolling molecular machines by straining graphene substrate

The surface rolling molecular machines are proposed to perform tasks and carrying molecular payloads on the substrates. As a result, controlling the surface motion of these molecular machines is of interest for the design of nano-transportation systems. In this study, we evaluate the motion of the nanocar on the graphene nanoribbons with strain gradient, through the molecular dynamics (MD) simulations, and theoretical relations. The nanocar indicates directed motion from the maximum strained part of the graphene to the unstrained end of the substrate. The strain gradient induced driving force and diffusion coefficients of nanocars are analyzed from the simulation and theoretical points of view. To obtain the optimum directed motion of nanocar, we consider the effects of temperature, strain average, and magnitude of strain gradient on the directionality of the motion. Moreover, the mechanism of the motion of nanocar is studied by evaluating the direction of the nanocar's chassis and the rotation of wheels around the axles. Ultimately, the programmable motion of nanocar is shown by adjusting the strain gradient of graphene substrate.

Covalent Triazine Framework C6N6 as an Electrochemical Sensor for Hydrogen-Containing Industrial Pollutants. A DFT Study

Industrial pollutants pose a serious threat to ecosystems. Hence, there is a need to search for new efficient sensor materials for the detection of pollutants. In the current study, we explored the electrochemical sensing potential of a C₆N₆ sheet for H-containing industrial pollutants (HCN, H₂S, NH₃ and PH₃) through DFT simulations. The adsorption of industrial pollutants over C₆N₆ occurs through physisorption, with adsorption energies ranging from -9.36 kcal/mol to -16.46 kcal/mol. The non-covalent interactions of analyte@C₆N₆ complexes are quantified by symmetry adapted perturbation theory (SAPT0), quantum theory of atoms in molecules (QTAIM) and non-covalent interaction (NCI) analyses. SAPT0 analyses show that electrostatic and dispersion forces play a dominant role in the stabilization of analytes over C₆N₆ sheets. Similarly, NCI and QTAIM analyses also verified the results of SAPT0 and interaction energy analyses. The electronic properties of analyte@C₆N₆ complexes are investigated by electron density difference (EDD), natural bond orbital analyses (NBO) and frontier molecular orbital analyses (FMO). Charge is transferred from the C₆N₆ sheet to HCN, H₂S, NH₃ and PH₃. The highest exchange of charge is noted for H₂S (-0.026 e^-). The results of FMO analyses show that the interaction of all analytes results in changes in the E_H-L gap of the C₆N₆ sheet. However, the highest decrease in the E_H-L gap (2.58 eV) is observed for the NH₃@C₆N₆ complex among all studied analyte@C₆N₆ complexes. The orbital density pattern shows that the HOMO density is completely concentrated on NH₃, while the LUMO density is centred on the C₆N₆ surface. Such a type of electronic transition results in a significant change in the E_H-L gap. Thus, it is concluded that C₆N₆ is highly selective towards NH₃ compared to the other studied analytes.

Discovering atomistic pathways for supply of metal atoms from methyl-based precursors to graphene surface

Conceptual 2D group III nitrides and oxides (e.g., 2D InN and 2D InO) in heterostructures with graphene have been realized by metal-organic chemical vapor deposition (MOCVD). MOCVD is expected to bring forth the same impact in the advancement of 2D semiconductor materials as in the fabrication of established semiconductor materials and device heterostructures. MOCVD employs metal-organic precursors such as trimethyl-indium, -gallium, and -aluminum, with (strong) metal-carbon bonds. Mechanisms that regulate MOCVD processes at the atomic scale are largely unknown. Here, we employ density-functional molecular dynamics - accounting for van der Waals interactions - to identify the reaction pathways responsible for dissociation of the trimethylindium (TMIn) precursor in the gas phase as well as on top-layer and zero-layer graphene. The simulations reveal how collisions with hydrogen molecules, intramolecular or surface-mediated proton transfer, and direct TMIn/graphene reactions assist TMIn transformations, which ultimately enables delivery of In monomers or InH and CH₃In admolecules, on graphene. This work provides knowledge for understanding the nucleation and intercalation mechanisms at the atomic scale and for carrying out epitaxial growth of 2D materials and graphene heterostructures.

Machine Learning Applications for Chemical Reactions

Machine learning (ML) approaches have enabled rapid and efficient molecular property predictions as well as the design of new novel materials. In addition to great success for molecular problems, ML techniques are applied to various chemical reaction problems that require huge costs to solve with the existing experimental and simulation methods. In this review, starting with basic representations of chemical reactions, we summarized recent achievements of ML studies on two different problems; predicting reaction properties and synthetic routes. The various ML models are used to predict physical properties related to chemical reaction properties (e. g. thermodynamic changes, activation barriers, and reaction rates). Furthermore, the predictions of reactivity, self-optimization of reaction, and designing retrosynthetic reaction paths are also tackled by ML approaches. Herein we illustrate various ML strategies utilized in the various context of chemical reaction studies.

Collective movement and thermal stability of fullerene clusters on the graphene layer

Understanding the motion characteristics of fullerene clusters on the graphene surface is critical for designing surface manipulation systems. Toward this purpose, using the molecular dynamics method, we evaluated six clusters of fullerenes including 1, 2, 3, 5, 10, and 25 molecules on the graphene surface, in the temperature range of 25 to 500 K. First, the surface motion of clusters is studied at 200 K and lower temperatures, in which fullerenes remain as a single group. The trajectories of the motion as well as the diffusion coefficients indicate the reduction of surface mobility as a response to the increase of the fullerene number. The clusters show normal diffusion at the temperature of 25 K, while they follow the super-diffusion regime at higher temperatures. The separation of fullerenes occurs at 300 K and higher temperatures. Due to the increase of vdW attraction with the increase of the fullerene number, the separation of fullerenes in larger clusters occurs at higher temperatures. The thermal energy at 500 K is sufficient to divide the large C60 clusters into smaller clusters. This energy level is related to the saturation of the interaction energy experienced by individual fullerenes, which can be estimated from the potential energy analysis. The results of simulations reveal that the separation occurs at the edge of clusters. Moreover, we studied the thermal stability of multilayer fullerene clusters on graphene. The simulation results indicate the tendency of multilayer clusters to locate on the surface, which implies the wetting property of C60s on the graphene layer.

Structural dynamics of Ru clusters during nitrogen dissociation in ammonia synthesis

The dynamic evolution of catalyst structures greatly influences the reactivity, especially sub-nanometer clusters, exhibiting complex configurational fluctuation. In the present work, we study the structural dynamics of a Ru₁₉ cluster during the dissociation of N₂ and calculate the reaction free energies using ab initio molecular dynamics (AIMD). Our AIMD calculation predicts a peak-shaped reaction entropy curve due to the adsorption-induced phase transition of the Ru₁₉ cluster. The low melting points of sub-nanometer clusters make it possible to activate N₂ at low temperatures. This work demonstrates that the dynamic changes of cluster structures have a non-negligible effect on reaction free energy and offer an opportunity for achieving ammonia synthesis under mild conditions.

Transition metal decorated phthalocyanine as a potential host material for lithium polysulfides: a first-principles study

The shuttle effect caused by the soluble long-chain lithium polysulfides greatly hinders the practical application of lithium–sulfur (Li–S) batteries. Therefore, the introduction of suitable anchoring materials is more effective to mitigate this problem. Transition metal phthalocyanines (TMPc) are regarded as a new class of sulfur host materials. Here, 4d transition metal (Y, Zr, Nb, Mo, Tc, Ru, Rh, Pd, Ag, Cd) decorated phthalocyanines are designed and systematically researched for the performance analysis of anchoring S₈/LiPSs by first-principles calculations. The results reveal that the bonding strength of LiPSs can be well adjusted by introducing suitable 4d transition metals into the phthalocyanine structure. The electronic structure analysis indicates the formation of TM–S bonds between the TMPc substrate materials and the LiPSs, which is essential to weaken the Li–S bonds and hence slow down the shuttle effect of LiPSs. ZrPc and NbPc both exhibit excellent potential and thermal stability for facilitating the conversion of LiPSs, as well as a better promoting effect for the sulfur reduction reactions (SRR) with a reduced Gibbs free energy in the rate-determining step (*Li₂S₂ → *Li₂S) during the discharge reaction process. These findings in our work may encourage further experimental and theoretical research for anchoring LiPSs with TMPc as a host material.

DFT calculations reveal that TMPc, especially ZrPc and NbPc exhibit the best anchoring and catalytic effects for lithium polysulfides.

Understanding the CH4 Conversion over Metal Dimers from First Principles

Inspired by the advantages of bi-atom catalysts and recent exciting progresses of nanozymes, by means of density functional theory (DFT) computations, we explored the potential of metal dimers embedded in phthalocyanine monolayers (M₂-Pc), which mimics the binuclear centers of methane monooxygenase, as catalysts for methane conversion using H₂O₂ as an oxidant. In total, 26 transition metal (from group IB to VIIIB) and four main group metal (M = Al, Ga, Sn and Bi) dimers were considered, and two methane conversion routes, namely *O-assisted and *OH-assisted mechanisms were systematically studied. The results show that methane conversion proceeds via an *OH-assisted mechanism on the Ti₂-Pc, Zr₂-Pc and Ta₂-Pc, a combination of *O- and *OH-assisted mechanism on the surface of Sc₂-Pc, respectively. Our theoretical work may provide impetus to developing new catalysts for methane conversion and help stimulate further studies on metal dimer catalysts for other catalytic reactions.

Adsorption Free Energy of Cellulose Nanocrystal on Water–Oil Interface

To investigate the amphiphilicity of cellulose, a series of molecular dynamics simulations were performed with a cellulose nanocrystal and a water–octane interfacial system. Assuming that the axis of cellulose is parallel to the water–octane interface, the freedoms of motion of the nanocrystal were restricted to two, the distance from the interface and the orientation around the axis. The mean force and the mean torque on the nanocrystal were evaluated with sufficiently long simulation at each crystal configuration, and their numerical integration gave a smooth free energy surface as the potential of mean force. The cellulose sample used here was found to be much more hydrophilic than oleophilic with the free energy difference ΔFw→o=318 kcal/mol. Three adsorption states with local minimum of adsorption free energy are distinguished in the free energy surface—the direct contact type which is similar to previously reported one, the hydrophilic-surface/water/octane type where a thin water layer is sandwiched between the surface and the octane phase, and the oleophilic/water/octane type where a thin water layer also exists. Water molecules in these water layers contribute to stabilize the adsorption states by taking a special orientational order and slow self-diffusion.

Discovering surface reaction pathways using accelerated molecular dynamics and network analysis tools

We present an automated method that maps surface reaction pathways with no experimental data and with minimal human interventions. In this method, bias potentials promoting surface reactions are applied to enable statistical samplings of the surface reaction within the timescale of ab initio molecular dynamics (MD) simulations, and elementary reactions are extracted automatically using an extension of the method constructed for gas- or liquid-phase reactions. By converting the extracted elementary reaction data to directed graph data, MD trajectories can be efficiently mapped onto reaction pathways using a network analysis tool. To demonstrate the power of the method, it was applied to the steam reforming of methane on the Rh(111) surface and to propane reforming on the Pt(111) and Pt₃Sn(111) surfaces. We discover new energetically favorable pathways for both reactions and reproduce the experimentally-observed materials-dependence of the surface reaction activity and the selectivity for the propane reforming reactions.

We present an automated method that maps surface reaction pathways with no experimental data and with minimal human interventions.

Interfacial acidity on the strontium titanate surface: a scaling paradigm and the role of the hydrogen bond

A fundamental understanding of acidity at an interface, as mediated by structure and molecule-surface interactions, is essential to elucidate the mechanisms of a range of chemical transformations. While the strength of an acid in homogeneous gas and solution phases is conceptually well understood, acid-base chemistry at heterogeneous interfaces is notoriously more complicated. Using density functional theory and nonlinear vibrational spectroscopy, we present a method to determine the interfacial Brønsted-Lowry acidity of aliphatic alcohols adsorbed on the (100) surface of the model perovskite, strontium titanate. While shorter and less branched alkanols are known to be less acidic in the gas phase and more acidic in solution, here we show that shorter alcohols are less acidic whereas less substituted alkanols are more acidic at the gas-oxide interface. Hydrogen bonding plays a critical role in defining acidity, whereas structure-acidity relationships are dominated by van der Waals interactions between the alcohol and the surface.

Adsorption of water in Na-LTA zeolites: an ab initio molecular dynamics investigation

The very wide range of applications of LTA zeolites, including the storage of tritiated water, implies that a detailed and accurate atomic-scale description of the adsorption processes taking place in their structure is crucial. To unravel with an unprecedented accuracy the mechanisms behind the water filling in NaA, we have conducted a systematic ab initio molecular dynamics investigation. Two LTA structural models, the conventional Z4A and the reduced one ZK4, have been used for static and dynamic ab initio calculations, respectively. After assessing this reduced model with comparative static DFT calculations, we start the filling of the α and β cages by water, molecule by molecule. This allowed us to thoroughly study the interaction of water molecules with the zeolite structure and between water molecules, progressively forming H-bond chains and ring patterns as the cage is being filled. The adsorption energies could then be calculated with an unprecedented accuracy, which showed that the interaction of the molecules with the zeolite weakens as their number increases. By these methods, we have been able to highlight the primary role of Na⁺ cations in the interaction of water with zeolite, and inversely, the role of water in the displacement of cations when it is sufficiently solvated, allowing the passage between the α and β cages. This phenomenon is possible thanks to the inhomogeneous distribution of water molecules on the cationic sites, as shown by our AIMD simulations, which allows the formation of water clusters. These results are important because they help in understanding how the coverage of cationic sites by water will affect the adsorption of other molecules inside the Na-LTA zeolite.

Oxidative etching of S-vacancy defective MoS2 monolayer upon reaction with O2

The reactions of O2 with S vacancy sites within a MoS2 monolayer were investigated using density functional theory calculations. We considered the following defects: single S vacancy, double S vacancy, two adjacent S vacancies and two S vacancies separated by a sulphur atom. We found that the surface distribution of S vacancy sites plays a key role in determining the surface reactivity towards O2. We observed the desorption of SO2 only for the last vacancy distribution. For the other cases, the surface becomes passivated with very stable structures having O atoms on the original vacancy sites and in some cases an SO group in an adjacent position. The ab initio molecular dynamics simulations showed that the impingement of the O2 molecule on an S vacancy site produces a stable chemisorbed O2 molecule with an upright configuration. The surface reactions initiate after the O2 molecule switches to the lying-down configuration which favours the breakage of the O-O bond and the concurrent formation of S-O bonds. In the most reactive vacancy site configuration, the dissociation of the first O2 molecule produces an SO intermediate which finally leads to desorption of SO2 after oxygen abstraction from the other adjacent O2 molecule. The formation of a MoO3 moiety within the monolayer was also observed in the molecular dynamic simulations at higher oxidation levels.

Reaction mechanism of atomic layer deposition of aluminum sulfide using trimethylaluminum and hydrogen sulfide

Atomic layer deposition (ALD) is a nanopreparation technique for materials and is widely used in the fields of microelectronics, energy and catalysis. ALD methods for metal sulfides, such as Al2S3 and Li2S, have been developed for lithium-ion batteries and solid-state electrolytes. In this work, using density functional theory calculations, the possible reaction pathways of the ALD of Al2S3 using trimethylaluminum (TMA) and H2S were investigated at the M06-2X/6-311G(d, p) level. Al2S3 ALD can be divided into two consecutive and complementary half-reactions involving TMA and H2S, respectively. In the TMA half-reaction, the methyl group can be eliminated through the reaction with the sulfhydryl group on the surface. This process is a ligand exchange reaction between the methyl and sulfhydryl groups via a four-membered ring transition state. TMA half-reaction with the sulfhydrylated surface is more difficult than that with the hydroxylated surface. When the temperature increases, the reaction requires more energy, owing to the contribution of the entropy. In the H2S half-reaction, the methyl group on the surface can further react with the H2S precursor via a four-membered ring transition state. The orientation of H2S and more molecules have minimal effect on the H2S half-reaction. The reaction involving H2S through a six-membered ring transition state is unfavorable. In addition, the methyl and sulfhydryl groups on the surface can both react with the adjacent sulfhydryl group on the subsurface to form and release CH4 or H2S in the two half-reactions. Furthermore, sulfhydryl elimination occurs more easily than methyl elimination on the surface. These findings for the TMA and H2S half-reactions of Al2S3 ALD may be used for studying precursor chemistry and improvements in the preparation of other metal sulfides for emerging applications.

Hybrid Nanostructures Obtained by Transport and Condensation of Tungsten Oxide Vapours onto CNW Templates

We present hybrid nanomaterial architectures, consisting of carbon nanowalls (CNW) templates decorated with tungsten oxide nanoparticles, synthesized using a mechanism based on tungsten oxide sublimation, vapor transport, followed by vapor condensation, in the absence or presence of plasma. The key steps in the decoration mechanism are the sublimation of tungsten oxides, when are exposed in vacuum at high temperature (800 °C), and their redeposition on colder surfaces (400-600 °C). The morphology and chemical composition of the hybrid architectures, as obtained from Scanning Electron Microscopy and X-ray Photoelectron Spectroscopy, are discussed with respect to substrate nature and the physical conditions of synthesis. We pointed out that the decoration process is strongly dependent on the temperature of the CNW templates and plasma presence. Thus, the decoration process performed with plasma was effective for a wider range of template temperatures, in contrast with the decoration process performed without plasma. The results are useful for applications using the sensing and photochemical properties of tungsten oxides, and have also relevance for fusion technology, tungsten walls erosion and material redeposition being widely observed in fusion machines.

Characterization of the mechanical properties of van der Waals heterostructures of stanene adsorbed on graphene, hexagonal boron-nitride and silicon carbide

Stanene has revealed a new horizon in the field of quantum condensed matter and energy conversion devices but its significantly lower tensile strength limits its further applications and effective operation in these nanodevices. Van der Waals heterostructures have given substantial flexibility to integrate different two-dimensional (2D) layered materials over the past few years and have proven highly functional with exceptional features, appealing applications, and innovative physics. Considerable efforts have been made for the preparation, thorough understanding, and applications of van der Waals heterostructures in the fields of electronics and optoelectronics. In this paper, we have executed Molecular Dynamics (MD) simulations to predict the tensile strength of van der Waals heterostructures of stanene (Sn) adsorbed on graphene (Gr), hexagonal boron nitride (hBN), and silicon carbide (SiC) (Sn/Gr, Sn/hBN, and Sn/SiC, respectively) subjected to both armchair and zigzag directional loading at different strain rates for the first time, which has enticing applications in electronic, optoelectronic, energy storage and bio-engineered devices. Among all the van der Waals heterostructures, the Sn/SiC heterostructure exhibits the lowest tensile strength and tensile strain. Furthermore, it has been found that zigzag directional loading could endure more tensile strain before fracture. Besides, it has been disclosed that though the rule of mixtures may accurately reproduce the Young's modulus of these heterostructures, it has limitations to predict the tensile strength. Fracture analysis suggests that for the Sn/hBN heterostructure the fracture initiates from the stanene layer while for the Sn/Gr and Sn/SiC heterostructures the fracture initiates from the Gr and SiC layer, respectively, for both armchair and zigzag directional loading. Overall, this study would aid in the design and efficient operation of Sn/Gr, Sn/hBN, and Sn/SiC heterostructures when subjected to mechanical force.

Mass production of metal-doped graphene from the agriculture waste of Quercus ilex leaves for supercapacitors: inclusive DFT study

This work reports a facile, eco-friendly, and cost-effective mass-scale synthesis of metal-doped graphene sheets (MDGs) using agriculture waste of Quercus ilex leaves for supercapacitor applications. A single step-degradation catalyst-based pyrolysis route was used for the manufacture of MDGs. Obtained MDGs were further evaluated via advanced spectroscopy and microscopic techniques including Raman spectroscopy, FT-IR, XRD, SEM/EDX, and TEM imaging. The Raman spectrum showed D and G bands at 1300 cm^-1 and 1590 cm^-1, respectively, followed by a 2D band at 2770 cm^-1, which confirmed the synthesis of few-layered MDGs. The SEM/EDX data confirmed the presence of 6.15%, 3.17%, and 2.36% of potassium, calcium and magnesium in the obtained MDGs, respectively. Additionally, the FT-IR, XRD, TEM, and SEM data including the plot profile diagrams confirmed the synthesis of MDGs. Further, a computational study was performed for the structural validation of MDGs using Gaussian 09. The density functional theory (DFT) results showed a chemisorption/decoration pattern of doping for metal ions on the few-layered graphene nanosheets, rather than a substitutional pattern. Further, resulting MDGs were used as an active material for the fabrication of a supercapacitor electrode using the polymer gel of PVA-H₃PO₄ as the electrolyte. The fabricated device showed a decent specific capacitance of 18.2 F g^-1 at a scan rate of 5 mV s^-1 with a power density of 1000 W kg^-1 at 5 A g^-1.

An In-Depth Assessment of the Electronic and Magnetic Properties of a Highly Ordered Hybrid Interface: The Case of Nickel Tetra-Phenyl-Porphyrins on Fe(001)-p(1 × 1)O

In this paper we focus on the structural, electronic, and magnetic properties of Ni tetra-phenyl-porphyrins (NiTPP) grown on top of Fe(001)–p(1 × 1)O. Ordered thin films of metal TPP molecules are potentially interesting for organic electronic and spintronic applications, especially when they are coupled to a ferromagnetic substrate. Unfortunately, porphyrin layers deposited on top of ferromagnetic substrates do not generally show long-range order. In this work, we provide evidence of an ordered disposition of the organic film above the iron surface and we prove that the thin layer of iron oxide decouples the molecules from the substrate, thus preserving the molecular electronic features, especially the HOMO-LUMO gap, even when just a few organic layers are deposited. The effect of the exposure to molecular oxygen is also investigated and an increased robustness against oxidation with respect to the bare substrate is detected. Finally, we present our results for the magnetic analysis performed by spin resolved spectroscopy, finding a null magnetic coupling between the molecules and the substrate.

From Behavior of Water on Hydrophobic Graphene Surfaces to Ultra-Confinement of Water in Carbon Nanotubes

In recent years and with the achievement of nanotechnologies, the development of experiments based on carbon nanotubes has allowed to increase the ionic permeability and/or selectivity in nanodevices. However, this new technology opens the way to many questionable observations, to which theoretical work can answer using several approximations. One of them concerns the appearance of a negative charge on the carbon surface, when the latter is apparently neutral. Using first-principles density functional theory combined with molecular dynamics, we develop here several simulations on different systems in order to understand the reactivity of the carbon surface in low or ultra-high confinement. According to our calculations, there is high affinity of the carbon atom to the hydrogen ion in every situation, and to a lesser extent for the hydroxyl ion. The latter can only occur when the first hydrogen attack has been achieved. As a consequence, the functionalization of the carbon surface in the presence of an aqueous medium is activated by its protonation, then allowing the reactivity of the anion.

MOCVD of AlN on epitaxial graphene at extreme temperatures

Appearance of luminescent centers with narrow spectral emission at room temperature in nanometer thin AlN is reported.

Indium Nitride at the 2D Limit

The properties of 2D InN are predicted to substantially differ from the bulk crystal. The predicted appealing properties relate to strong in- and out-of-plane excitons, high electron mobility, efficient strain engineering of their electronic and optical properties, and strong application potential in gas sensing. Until now, the realization of 2D InN remained elusive. In this work, the formation of 2D InN and measurements of its bandgap are reported. Bilayer InN is formed between graphene and SiC by an intercalation process in metal-organic chemical vapor deposition (MOCVD). The thickness uniformity of the intercalated structure is investigated by conductive atomic force microscopy (C-AFM) and the structural properties by atomic resolution transmission electron microscopy (TEM). The coverage of the SiC surface is very high, above 90%, and a major part of the intercalated structure is represented by two sub-layers of indium (In) bonded to nitrogen (N). Scanning tunneling spectroscopy (STS) measurements give a bandgap value of 2 ± 0.1 eV for the 2D InN. The stabilization of 2D InN with a pragmatic wide bandgap and high lateral uniformity of intercalation is demonstrated.

Ab initio study for molecular-scale adsorption, decomposition and desorption on AlN surfaces during MOCVD growth

Since AlGaN offers new opportunities for the development of the solid state ultraviolet (UV) luminescence, detectors and high-power electronic devices, the growth of AlN buffer substrate is concerned. However, the growth of AlN buffer substrate during MOCVD is regulated by an intricate interplay of gas-phase and surface reactions that are beyond the resolution of experimental techniques, especially the surface growth process. We used density-functional ab initio calculations to analyze the adsorption, decomposition and desorption of group-III and group-V sources on AlN surfaces during MOCVD growth in molecular-scale. For AlCH₃ molecule the group-III source, the results indicate that AlCH₃ is more easily adsorbed on AlN (0001) than (000\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\overline{1}$$\end{document}1¯) surface on the top site. For the group-V source decomposition we found that NH₂ molecule is the most favorable adsorption source and adsorbed on the top site. We investigated the adsorption of group-III source on the reconstructed AlN (0001) surface which demonstrates that NH₂-rich condition has a repulsion effect to it. Furthermore, the desorption path of group-III and group-V radicals has been proposed. Our study explained the molecular-scale surface reaction mechanism of AlN during MOCVD and established the surface growth model on AlN (0001) surface.

A Cellular Automata Approach for the Modeling of a Polyamide and Carbon Aerogel Structure and Its Properties

In this work, a cellular automata (CA) approach was used to generate 3D structures of polyamide and carbon aerogels. Experimental results are used as initial data for materials’ digital representations and to verify the developed CA models. Based on the generated digital structures, a computer study of aerogels’ mechanical properties was conducted. The offered CA models can be applied for the development of new nanoporous materials such as aerogels of different nature and allow for a reduction in the amount of required full-scale experiments, consequently decreasing development time and costs of new material formulations.

Nanoscale phenomena ruling deposition and intercalation of AlN at the graphene/SiC interface

The possibility for kinetic stabilization of prospective 2D AlN was explored by rationalizing metal organic chemical vapor deposition (MOCVD) processes of AlN on epitaxial graphene. From the wide range of temperatures which can be covered in the same MOCVD reactor, the deposition was performed at the selected temperatures of 700, 900, and 1240 °C. The characterization of the structures by atomic force microscopy, electron microscopy and Raman spectroscopy revealed a broad range of surface nucleation and intercalation phenomena. These phenomena included the abundant formation of nucleation sites on graphene, the fragmentation of the graphene layers which accelerated with the deposition temperature, the delivery of excess precursor-derived carbon adatoms to the surface, as well as intercalation of sub-layers of aluminum atoms at the graphene/SiC interface. The conceptual understanding of these nanoscale phenomena was supported by our previous comprehensive ab initio molecular dynamics (AIMD) simulations of the surface reaction of trimethylaluminum, (CH₃)₃Al, precursor with graphene. A case of applying trimethylindium, (CH₃)₃In, precursor to epitaxial graphene was considered in a comparative way.

Graphene Oxide Membrane Immobilized Aptamer as a Highly Selective Hormone Removal

Three-dimensional (3D) reduced graphene oxide (rGO) modified by polyethyleneimine (PEI) was prepared and functionalized by fluorophore-labeled dexamethasone-aptamer (Flu-DEX-apt) via π-π stacking interaction. The rGO/PEI/Flu-DEX-apt was used as a selective membrane for dexamethasone hormone removal from water. The prepared rGO/PEI/Flu-DEX-apt membranes were stable, insoluble, and easily removable from liquid media. The membrane was characterized by Raman spectroscopy, scanning electron spectroscopy, and FTIR spectroscopy. The rGO/PEI/Flu-DEX-apt membrane showed high sensitivity and specificity toward the dexamethasone hormone in the presence of other steroid hormone analogs, such as progesterone, estrone, estradiol, and 19-norethindrone. The fluorescence and UV-visible spectroscopy were used to confirm the membranes performance and the quantification of hormones removal. The resulting data clearly show that the graphene oxide concentration influence the aptamers and analytes interaction (π-π stacking interaction). It was found that by varying the graphene oxide concentration yields to different porosities of rGO/PEI/Flu-DEX-apt membranes affects the adsorption recovery rate, as well as the specificity and selectivity toward the dexamethasone hormone.

The C2N surface as a highly selective sensor for the detection of nitrogen iodide from a mixture of NX3 (X = Cl, Br, I) explosives

Explosives are quite toxic and destructive; therefore, it is necessary to not only detect them but also remove them. The adsorption behavior of NX₃ analytes (NCl₃, NBr₃ and NI₃) over the microporous C₂N surface was evaluated by DFT calculations. The nature of interactions between NX₃ and C₂N was characterized by adsorption energy, NCI, QTAIM, SAPT0, NBO, EDD and FMO analysis. The interaction energies of NX₃ with C₂N are in the range of −10.85 to −16.31 kcal mol⁻¹ and follow the order of NCl₃@C₂N > NBr₃@C₂N > NI₃@C₂N, respectively. The 3D isosurfaces and 2D-RGD graph of NCI analysis qualitatively confirmed the existence of halogen bonding interactions among the studied systems. Halogen bonding was quantified by SAPT0 component energy analysis. The SAPT0 results revealed that ΔE_disp (56.75%) is the dominant contributor towards interaction energy, whereas contributions from ΔE_elst and ΔE_ind are 29.41% and 14.34%, respectively. The QTAIM analysis also confirmed the presence of halogen bonding between atoms of NX₃ and C₂N surface. EDD analysis also validated NCI, QTAIM and NBO analysis. FMO analysis revealed that the adsorption of NI₃ on the C₂N surface caused the highest change in the E_HOMO–LUMO gap (from 5.71 to 4.15 eV), and resulted in high sensitivity and selectivity of the C₂N surface towards NI₃, as compared to other analytes. It is worth mentioning that in all complexes, a significant difference in the E_HOMO–LUMO gap was seen when electronic transitions occurred from the analyte to the C₂N surface.

Explosives are quite toxic and destructive; therefore, it is necessary to not only detect them but also remove them.

The Reaction Thermodynamics during Plating Al on Graphene Process

This research explored a novel chemical reduction of organic aluminum for plating Al on a graphene surface. The thermodynamics of the Al plating reaction process were studied. The Al plating process consisted of two stages: the first was to prepare (C₂H₅)₃Al. In this reaction, the ΔH(enthalpy) was 10.64 kcal/mol, the ΔG(Gibbs free energy) was 19.87 kcal/mol and the ΔS(entropy) was 30.9 cal/(mol·K); this was an endothermic reaction. In the second stage, the (C₂H₅)₃Al decomposed into Al atoms, which were gradually deposited on the surface of the graphene and the Al plating formed. At 298.15 K, the ΔH was −20.21 kcal/mol, the ΔG was −54.822 kcal/mol, the ΔS was 116.08 cal/(mol·K) and the enthalpy change was negative, thus indicating an endothermic reaction.

Epitaxial Growth of Orthorhombic GaFeO₃ Thin Films on SrTiO₃ (111) Substrates by Simple Sol-Gel Method

A Sol-gel method assisted with spin-coating has been successfully used to grow orthorhombic GaFeO₃ epitaxial films on SrTiO₃ (111) substrates for the first time. The film with Pna2₁ crystal structure has been grown along the c-axis. The rocking curve of (004) reflection shows that the Full-Width at Half-Maximum (FWHM) value could be determined to be 0.230°, indicating good single crystallinity and high quality. X-ray Φ scan reveals a three-fold symmetry of the substrate and a six-fold symmetry of the film, respectively. The in-plane domains rotate 60° from each other in the film. Uniform film with dense structure, columnar grains with similar grain size was obtained. The thickness of the film was evaluated to be ~170 nm. The roughness value (RMS) measured by AFM was 4.5 nm, revealing a flat film. The in-plane Magnetization versus Magnetic field (M-H) curve at 5 K performs a typical ferri- or ferromagnetic hysteresis loop with a saturated magnetization (M_s) value of 136 emu/cm³. The Curie temperature could be determined to be 174 K. Compared to Pulsed Laser Deposition (PLD), the sol-gel method can prepare large area films with low cost. These new films show promising applications in multiferroic devices.

Molecular Dynamics Study on the Reverse Osmosis Using Multilayer Porous Graphene Membranes

In this study, the reverse osmosis (RO) of a salt solution was investigated using a molecular dynamics method to explore the performance of a multilayer porous graphene membrane. The effects of the salt solution concentration, pressure, layer separation and pore offset on the RO performance of the membrane were investigated and the influences of the number of layers and the gradient structure were determined. The results show that as the salt solution concentration increases, the energy barrier of the water molecules passing through the bilayer porous graphene membranes changes slightly, indicating that the effect of the water flux on the membrane can be ignored. The salt rejection performance of the membrane improves with an increase in the concentration of the salt solution. When the pressure is increased, the energy barrier decreases, the water flux increases and the salt rejection decreases. When the layer separation of the bilayer porous graphene membrane is the same as the equilibrium spacing of the graphene membrane, the energy barrier is the lowest and the membrane water flux is the largest. The energy barrier of the bilayer porous graphene membrane increases with increasing layer separation, resulting in a decrease in the water flux of the membrane. The salt rejection increases with increasing layer separation. The water flux of the membrane decreases as the energy barrier increases with increasing pore offset and the salt rejection increases. The energy barrier effect is more pronounced for a larger number of graphene layers and the water flux of the membrane decreases because it is more difficult for the water molecules to pass through the porous graphene membrane. However, the salt rejection performance improves with the increase in the number of layers. The gradient pore structure enhances the energy barrier effect of the water molecules permeating through the membrane and the water flux of the membrane decreases. The salt rejection performance is improved by the gradient pore structure. The research results provide theoretical guidance for research on the RO performance of porous graphene membranes and the design of porous graphene membranes.

Molecular Dynamics Investigation of Graphene Nanoplate Diffusion Behavior in Poly-α-Olefin Lubricating Oil

Graphene as a type of novel additive significantly enhanced the tribological performance of blended lubricating oil. However, the dispersibility of graphene with long-term stability in lubricating oil is still a challenge. Chemical modification for graphene, rather than using surfactants, provided a better method to improve the dispersibility of graphene in lubricants. In this study, the equilibrium molecular dynamics (EMD) simulations were carried out to investigate the diffusion behavior of graphene nanoplates in poly-α-olefin (PAO) lubricating oil. The effects of graphene-size, edge-functionalization, temperature, and pressure on the diffusion coefficient were studied. In order to understand the influence of edge-functionalization, three different functional groups were grafted to the edge of graphene nanoplates: COOH, COON(CH3)2, CONH(CH2)8CH3 (termed GO, MG, and AG, respectively). The EMD simulations results demonstrated that the relationships between diffusion coefficient and graphene-size and number of functional groups were linear while the temperature and pressure had a nonlinear influence on the diffusion coefficient. It was found that the larger dimension and more functional groups provided the lower diffusion coefficient. AG with eight CONH(CH2)8CH3 groups exhibited the lowest diffusion coefficient. Furthermore, the experimental results and radial distribution function for graphene-PAO illustrated that the diffusion coefficient reflected the dispersibility of nanoparticles in nanofluids to some degree. To our best knowledge, this study is the first time the diffusion behavior of graphene in PAO lubricating oil was investigated using EMD simulations.