Tip Growth of Quasi-Metallic Bilayer Graphene Nanoribbons with Armchair Chirality

Graphene nanoribbons (GNRs), quasi one-dimensional (1D) narrow strips of graphene, have shown promise for high-performance nanoelectronics due to their exceptionally high carrier mobility and structurally tunable bandgaps. However, producing chirality-uniform GNRs on insulating substrates remains a big challenge. Here, we report the successful growth of bilayer GNRs with predominantly armchair chirality and ultranarrow widths (<5 nm) on insulating hexagonal boron nitride (h-BN) substrates using chemical vapor deposition (CVD). The growth of GNRs is catalyzed by transition metal nanoparticles, including Fe, Co, and Ni, through a unique tip-growth mechanism. Notably, GNRs catalyzed by Ni exhibit a high purity (97.3%) of armchair chirality. Electron transport measurements indicate that the ultrathin bilayer armchair GNRs exhibit quasi-metallic behavior. This quasi-metallicity is further supported by density functional theory (DFT) calculations, which reveal a significantly reduced bandgap in bilayer armchair GNRs. The chirality-specific GNRs reported here offer promising advancements for the application of graphene in nanoelectronics.

A mechanism for thickness-controllable single crystalline 2D materials growth

Recent efforts in growing two-dimensional (2D) multilayers have enabled the synthesis of single crystalline 2D multilayers in a wafer scale through the seamless stitching of multiple epitaxial 2D islands. Unlike previously observed wedding-cake or inverted-wedding-cake structures, these multilayer islands have the same size and shape in each layer with aligned edges. In this study, we investigated the underlying growth mechanisms of synchronic 2D multilayers growth and have showed that a heterogenous layer on a crystalline substrate is critical for maintaining the synchronic growth of 2D multilayers. During growth, the heterogenous layer passivates the edges of multilayer 2D island and thus prevents the coalescence of these active edges, while the high interfacial energy between the heterogenous surface layer and the substrate stabilizes the synchronic structure. Based on this model, we have successfully explained the previously observed synchronic growth of graphene and hexagonal boron nitride multilayers (Nat Nanotech 2020, 15: 861; Nature 2022, 606: 88). The deep understanding on the mechanism paves a way towards the synthesis of wafer-scale single-crystal 2D multilayers with a uniform thickness.

Substrate Engineering for Chemical Vapor Deposition Growth of Large-Scale 2D Transition Metal Dichalcogenides

The large-scale production of 2D transition metal dichalcogenides (TMDs) is essential to realize their industrial applications. Chemical vapor deposition (CVD) has been considered as a promising method for the controlled growth of high-quality and large-scale 2D TMDs. During a CVD process, the substrate plays a crucial role in anchoring the source materials, promoting the nucleation and stimulating the epitaxial growth. It thus significantly affects the thickness, microstructure, and crystal quality of the products, which are particularly important for obtaining 2D TMDs with expected morphology and size. Here, an insightful review is provided by focusing on the recent development associated with the substrate engineering strategies for CVD preparation of large-scale 2D TMDs. First, the interaction between 2D TMDs and substrates, a key factor for the growth of high-quality materials, is systematically discussed by combining the latest theoretical calculations. Based on this, the effect of various substrate engineering approaches on the growth of large-area 2D TMDs is summarized in detail. Finally, the opportunities and challenges of substrate engineering for the future development of 2D TMDs are discussed. This review might provide deep insight into the controllable growth of high-quality 2D TMDs toward their industrial-scale practical applications.

Growth mechanisms of monolayer hexagonal boron nitride (h-BN) on metal surfaces: theoretical perspectives

Two-dimensional hexagonal boron nitride (h-BN) has appeared as a promising material in diverse areas of applications, including as an excellent substrate for graphene devices, deep-ultraviolet emitters, and tunneling barriers, thanks to its outstanding stability, flat surface, and wide-bandgap. However, for achieving such exciting applications, controllable mass synthesis of high-quality and large-scale h-BN is a precondition. The synthesis of h-BN on metal surfaces using chemical vapor deposition (CVD) has been extensively studied, aiming to obtain large-scale and high-quality materials. The atomic-scale growth process, which is a prerequisite for rationally optimizing growth circumstances, is a key topic in these investigations. Although theoretical investigations on h-BN growth mechanisms are expected to reveal numerous new insights and understandings, different growth methods have completely dissimilar mechanisms, making theoretical research extremely challenging. In this article, we have summarized the recent cutting-edge theoretical research on the growth mechanisms of h-BN on different metal substrates. On the frequently utilized Cu substrate, h-BN development was shown to be more challenging than a simple adsorption-dehydrogenation-growth scenario. Controlling the number of surface layers is also an important challenge. Growth on the Ni surface is controlled by precipitation. An unusual reaction-limited aggregation growth behavior has been seen on interfaces having a significant lattice mismatch to h-BN. With intensive theoretical investigations employing advanced simulation approaches, further progress in understanding h-BN growth processes is predicted, paving the way for guided growth protocol design.

Tunable growth of one-dimensional graphitic materials: graphene nanoribbons, carbon nanotubes, and nanoribbon/nanotube junctions

Graphene nanoribbons (GNRs) and carbon nanotubes (CNTs), two representative one-dimensional (1D) graphitic materials, have attracted tremendous research interests due to their promising applications for future high-performance nanoelectronics. Although various methods have been developed for fabrication of GNRs or CNTs, a unified method allowing controllable synthesis of both of them, as well as their heterojunctions, which could largely benefit their nano-electronic applications, is still lacking. Here, we report on a generic growth of 1D carbon using nanoparticles catalyzed chemical vapor deposition (CVD) on atomically flat hexagonal boron nitride (h-BN) substrates. Relative ratio of the yielded GNRs and CNTs is able to be arbitrarily tuned by varying the growth temperature or feeding gas pressures. The tunability of the generic growth is quantitatively explained by a competing nucleation theory: nucleation into either GNRs or CNTs by the catalysts is determined by the free energy of their formation, which is controlled by the growth conditions. Under the guidance of the theory, we further realized growth of GNR/CNT intramolecular junctions through changing H₂ partial pressure during a single growth process. Our study provides not only a universal and controllable method for growing 1D carbon nanostructures, but also a deep understanding of their growth mechanism, which would largely benefit future carbon-based electronics and optoelectronics.

Edge Reconstruction-Dependent Growth Kinetics of MoS2

Understanding the growth mechanisms of multielement two-dimensional (2D) crystals is challenging because of the unbalanced stoichiometry and possible reconstruction of their edges. Here, we present a systematic theoretical study on the chemical vapor deposition (CVD) growth mechanism of MoS₂. We found that the growth kinetics of MoS₂ highly depends on its edge reconstruction determined by concentrations of Mo and S in the growth environment. Based on the calculated energies of nucleation and propagation of various MoS₂ edges, we predicted the transition of a MoS₂ island growth from a regime of a triangle enclosed by Mo-terminated zigzag edges that are passivated by 50% S (Mo-II edges), to a regime of continuous evolution within a triangle, hexagon, and inverted triangle with 75%-S-terminated edges (S-III edges) and Mo-II edges, and finally to a regime of triangles with Mo-terminated zigzag edges that are passivated by 100% S (Mo-III edges) by tuning the growth condition from Mo-rich to S-rich, which provides a reasonable explanation to many experimental observations. This study provides a general guideline on theoretical studies of 2D crystals' growth mechanisms, deepens our understanding on the growth mechanism of multielement 2D crystals, and is beneficial for the controllable synthesis of various 2D crystals.

Catalytic Growth of Ultralong Graphene Nanoribbons on Insulating Substrates

Graphene nanoribbons (GNRs) with widths of a few nanometers are promising candidates for future nanoelectronic applications due to their structurally tunable bandgaps, ultrahigh carrier mobilities, and exceptional stability. However, the direct growth of micrometer-long GNRs on insulating substrates, which is essential for the fabrication of nanoelectronic devices, remains an immense challenge. Here, the epitaxial growth of GNRs on an insulating hexagonal boron nitride (h-BN) substrate through nanoparticle-catalyzed chemical vapor deposition is reported. Ultranarrow GNRs with lengths of up to 10 µm are synthesized. Remarkably, the as-grown GNRs are crystallographically aligned with the h-BN substrate, forming 1D moiré superlattices. Scanning tunneling microscopy reveals an average width of 2 nm and a typical bandgap of ≈1 eV for similar GNRs grown on conducting graphite substrates. Fully atomistic computational simulations support the experimental results and reveal a competition between the formation of GNRs and carbon nanotubes during the nucleation stage, and van der Waals sliding of the GNRs on the h-BN substrate throughout the growth stage. This study provides a scalable, single-step method for growing micrometer-long narrow GNRs on insulating substrates, thus opening a route to explore the performance of high-quality GNR devices and the fundamental physics of 1D moiré superlattices.

Epitaxy of 2D Materials toward Single Crystals

Two-dimensional (2D) materials exhibit unique electronic, optical, magnetic, mechanical, and thermal properties due to their special crystal structure and thus have promising potential in many fields, such as in electronics and optoelectronics. To realize their real applications, especially in integrated devices, the growth of large-size single crystal is a prerequisite. Up to now, the most feasible way to achieve 2D single crystal growth is the epitaxy: growth of 2D materials of one or more specific orientations with single-crystal substrate. Only when the 2D domains have the same orientation, they can stitch together seamlessly and single-crystal 2D films can be obtained. In this view, four different epitaxy modes of 2D materials on various substrates are presented, including van der Waals epitaxy, edge epitaxy, step-guided epitaxy, and in-plane epitaxy focusing on the growth of graphene, hexagonal boron nitride (h-BN), and transition metal dichalcogenide (TMDC). The lattice symmetry relation and the interaction between 2D materials and the substrate are the key factors determining the epitaxy behaviors and thus are systematically discussed. Finally, the opportunities and challenges about the epitaxy of 2D single crystals in the future are summarized.

Sub-4 nm Nanodiamonds from Graphene-Oxide and Nitrated Polycyclic Aromatic Hydrocarbons at 423 K

Nanodiamonds are interesting materials from the point of view of their biocompatibility and their chemical, spectroscopic, and mechanical properties. Current synthetic methods for nanodiamonds involve harsh environments, which are potentially hazardous in addition to being expensive. We report a low-temperature (423 K) hydrothermal approach to form nanodiamonds by using graphene-oxide or nitrated polycyclic aromatic hydrocarbons (naphthalene, anthracene, phenanthrene, or pyrene) as a starting material. The reaction products contain single-crystalline or twinned nanodiamonds with average diameters in the 2-3 nm range. Theoretical calculations prove that, at the nanoscale, sub-4 nm nanodiamonds may adopt a structure that is more stable than graphene-oxide and nitrated polycyclic aromatic hydrocarbons. Our findings show that sp² carbon in the polycyclic aromatic precursor can be converted to sp³ carbon under unexpectedly moderate temperature conditions by using nanoscale precursors and thus offer a low-temperature approach for the synthesis of sub-4 nm nanodiamonds.

Universal Zigzag Edge Reconstruction of an α-Phase Puckered Monolayer and Its Resulting Robust Spatial Charge Separation

Edges are important, because they dictate the stability and properties of nanoribbons. Here, we reveal a universal reconstruction of the ZZ edge into a (2 × 1) tubed [ZZ(Tube)] edge, enabling an ultimate narrow nanotube to terminate nanoribbons for α-puckered group-V elemental and compound monolayers (GeS/Se and SnS/Se). The reconstructed edge formations are confirmed by CALYPSO. The ZZ(Tube) edge forms easily, is highly stable, and is semiconducting. Remarkably, the ZZ(Tube) edge always exhibits a type-II band structure and robust spatial charge separation. For a compound monolayer monochalcogenide, mild (2 × 1) ZZ(S-R) occurs at the chalcogenide-terminated edge. TDDFT simulations indicate that charge separation occurs only at 672 fs, while the lifetime is over 5 ns, thus facilitating robust spatial charge accumulation. These remarkable features of ZZ(Tube) edge-terminated α-puckered nanoribbons are ideal for optoelectronic and photocatalytic applications.

Theoretical Study of Chemical Vapor Deposition Synthesis of Graphene and Beyond: Challenges and Perspectives

Two-dimensional (2D) materials have attracted great attention in recent years because of their unique dimensionality and related properties. Chemical vapor deposition (CVD), a crucial technique for thin-film epitaxial growth, has become the most promising method of synthesizing 2D materials. Different from traditional thin-film growth, where strong chemical bonds are involved in both thin films and substrates, the interaction in 2D materials and substrates involves the van der Waals force and is highly anisotropic, and therefore, traditional thin-film growth theories cannot be applied to 2D material CVD synthesis. During the last 15 years, extensive theoretical studies were devoted to the CVD synthesis of 2D materials. This Perspective attempts to present a theoretical framework for 2D material CVD synthesis as well as the challenges and opportunities in exploring CVD mechanisms. We hope that this Perspective can provide an in-depth understanding of 2D material CVD synthesis and can further stimulate 2D material synthesis.

Strategies, Status, and Challenges in Wafer Scale Single Crystalline Two-Dimensional Materials Synthesis

The successful exfoliation of graphene has given a tremendous boost to research on various two-dimensional (2D) materials in the last 15 years. Different from traditional thin films, a 2D material is composed of one to a few atomic layers. While atoms within a layer are chemically bonded, interactions between layers are generally weak van der Waals (vdW) interactions. Due to their particular dimensionality, 2D materials exhibit special electronic, magnetic, mechanical, and thermal properties, not found in their 3D counterparts, and therefore they have great potential in various applications, such as 2D materials-based devices. To fully realize their large-scale practical applications, especially in devices, wafer scale single crystalline (WSSC) 2D materials are indispensable. In this review, we present a detailed overview on strategies toward the synthesis of WSSC 2D materials while highlighting the recent progress on WSSC graphene, hexagonal boron nitride (hBN), and transition metal dichalcogenide (TMDC) synthesis. The challenges that need to be addressed in future studies have also been described. In general, there have been two distinct routes to synthesize WSSC 2D materials: (i) allowing only one nucleus on a wafer scale substrate to be formed and developed into a large single crystal and (ii) seamlessly stitching a large number of unidirectionally aligned 2D islands on a wafer scale substrate, which is generally single crystalline. Currently, the synthesis of WSSC graphene has been realized by both routes, and WSSC hBN and MoS₂ have been synthesized by route (ii). On the other hand, the growth of other WSSC 2D materials and WSSC multilayer 2D materials still remains a big challenge. In the last section, we wrap up this review by summarizing the future challenges and opportunities in the synthesis of various WSSC 2D materials.

The epitaxy of 2D materials growth

Two dimensional (2D) materials consist of one to a few atomic layers, where the intra-layer atoms are chemically bonded and the atomic layers are weakly bonded. The high bonding anisotropicity in 2D materials make their growth on a substrate substantially different from the conventional thin film growth. Here, we proposed a general theoretical framework for the epitaxial growth of a 2D material on an arbitrary substrate. Our extensive density functional theory (DFT) calculations show that the propagating edge of a 2D material tends to align along a high symmetry direction of the substrate and, as a conclusion, the interplay between the symmetries of the 2D material and the substrate plays a critical role in the epitaxial growth of the 2D material. Based on our results, we have outlined that orientational uniformity of 2D material islands on a substrate can be realized only if the symmetry group of the substrate is a subgroup of that of the 2D material. Our predictions are in perfect agreement with most experimental observations on 2D materials' growth on various substrates known up to now. We believe that this general guideline will lead to the large-scale synthesis of wafer-scale single crystals of various 2D materials in the near future.

Kinetics of Graphene and 2D Materials Growth

During the last 10 years, remarkable achievements on the chemical vapor deposition (CVD) growth of 2D materials have been made, but the understanding of the underlying mechanisms is still relatively limited. Here, the current progress on the understanding of the growth kinetics of 2D materials, especially for their CVD synthesis, is reviewed. In order to present a complete picture of 2D materials' growth kinetics, the following factors are discussed: i) two types of growth modes, namely attachment-limited growth and diffusion-limited growth; ii) the etching of 2D materials, which offers an additional degree of freedom for growth control; iii) a number of experimental factors in graphene CVD synthesis, such as structure of the substrate, pressure of hydrogen or oxygen, temperature, etc., which are found to have profound effects on the growth kinetics; iv) double-layer and few-layer 2D materials' growth, which has distinct features different from the growth of single-layer 2D materials; and v) the growth of polycrystalline 2D materials by the coalescence of a few single crystalline domains. Finally, the current challenges and opportunities in future 2D materials' synthesis are summarized.

Catalytic CVD synthesis of boron nitride and carbon nanomaterials - synergies between experiment and theory

Low-dimensional carbon and boron nitride nanomaterials - hexagonal boron nitride, graphene, boron nitride nanotubes and carbon nanotubes - remain at the forefront of advanced materials research. Catalytic chemical vapour deposition has become an invaluable technique for reliably and cost-effectively synthesising these materials. In this review, we will emphasise how a synergy between experimental and theoretical methods has enhanced the understanding and optimisation of this synthetic technique. This review examines recent advances in the application of CVD to synthesising boron nitride and carbon nanomaterials and highlights where, in many cases, molecular simulations and quantum chemistry have provided key insights complementary to experimental investigation. This synergy is particularly prominent in the field of carbon nanotube and graphene CVD synthesis, and we propose here it will be the key to future advances in optimisation of CVD synthesis of boron nitride nanomaterials, boron nitride - carbon composite materials, and other nanomaterials generally.

Formation mechanism of overlapping grain boundaries in graphene chemical vapor deposition growth

The formation mechanisms of two different types of grain boundaries (GBs), the weakly bound overlapping GB and the covalent bound GB, during graphene domain coalescence are revealed by both theoretical modeling and experimental observations.

The formation of grain boundaries (GBs) in graphene films is both fundamentally interesting and practically important for many applications. A GB in graphene is known as a linear defect and is formed during the coalescence of two single crystalline graphene domains. The covalent binding between domains is broadly known as the mechanism of GB formation during graphene chemical vapor deposition (CVD) growth. Here, we demonstrate another GB formation mechanism, where two graphene domains are connected by weak van der Waals interactions between overlapping graphene layers. The formation mechanism of the overlapping GBs (OLGBs) is systematically explored theoretically and the proposed conditions for forming OLGBs are validated by experimental observations. This discovery leads to a deep understanding of the mechanism of graphene CVD growth and reveals potential means for graphene quality control in CVD synthesis.

Stacking sequence and interlayer coupling in few-layer graphene revealed by in situ imaging

In the transition from graphene to graphite, the addition of each individual graphene layer modifies the electronic structure and produces a different material with unique properties. Controlled growth of few-layer graphene is therefore of fundamental interest and will provide access to materials with engineered electronic structure. Here we combine isothermal growth and etching experiments with in situ scanning electron microscopy to reveal the stacking sequence and interlayer coupling strength in few-layer graphene. The observed layer-dependent etching rates reveal the relative strength of the graphene-graphene and graphene-substrate interaction and the resulting mode of adlayer growth. Scanning tunnelling microscopy and density functional theory calculations confirm a strong coupling between graphene edge atoms and platinum. Simulated etching confirms that etching can be viewed as reversed growth. This work demonstrates that real-time imaging under controlled atmosphere is a powerful method for designing synthesis protocols for sp2 carbon nanostructures in between graphene and graphite.

How Graphene Islands Are Unidirectionally Aligned on the Ge(110) Surface

The unidirectional alignment of graphene islands is essential to the synthesis of wafer-scale single-crystal graphene on Ge(110) surface, but the underlying mechanism is not well-understood. Here we report that the necessary coalignment of the nucleating graphene islands on Ge(110) surface is caused by the presence of step-pattern; we show that on the preannealed Ge(110) textureless surface the graphene islands appear nonpreferentially orientated, while on the Ge(110) surfaces with natural step pattern, all graphene islands emerge coaligned. First-principles calculations and theoretical analysis reveal this different alignment behaviors originate from the strong chemical binding formed between the graphene island edges and the atomic steps on the Ge(110) surface, and the lattice matching at edge-step interface dictates the alignment of graphene islands with the armchair direction of graphene along the [-110] direction of the Ge(110) substrate.

Large scale atomistic simulation of single-layer graphene growth on Ni(111) surface: molecular dynamics simulation based on a new generation of carbon-metal potential

To explore the mechanism of graphene chemical vapor deposition (CVD) growth on a catalyst surface, a molecular dynamics (MD) simulation of carbon atom self-assembly on a Ni(111) surface based on a well-designed empirical reactive bond order potential was performed. We simulated single layer graphene with recorded size (up to 300 atoms per super-cell) and reasonably good quality by MD trajectories up to 15 ns. Detailed processes of graphene CVD growth, such as carbon atom dissolution and precipitation, formation of carbon chains of various lengths, polygons and small graphene domains were observed during the initial process of the MD simulation. The atomistic processes of typical defect healing, such as the transformation from a pentagon into a hexagon and from a pentagon-heptagon pair (5|7) to two adjacent hexagons (6|6), were revealed as well. The study also showed that higher temperature and longer annealing time are essential to form high quality graphene layers, which is in agreement with experimental reports and previous theoretical results.

The reconstructed edges of the hexagonal BN

As an important two-dimensional material which shows exceptional mechanical and chemical stability, superior electronic properties, along with broad applications, the hexagonal-BN (h-BN) has drawn great attention recently. Here we report a systematic study on the structural stability, electronic and magnetic properties of various h-BN edges, including both bare and hydrogen-terminated ones. It is found that along the armchair (AC) direction, the pristine edge is the most stable one because of the formation of a triple B≡N bond, while, along the zigzag (ZZ) directions, the reconstructed ones, ZZB + N and ZZN57 are more stable. The pristine edges are more stable in bare BN in most cases if saturated with hydrogen. By applying the theory of Wulff construction, we predicted that an unpassivated BN domain prefers the hexagonal shape enclosed with bare AC edges i.e., AC-Ns, AC, AC-Bs if the feedstock varies from N-rich to B-rich. However, the evolution from ZZN edged triangular domain, to hexagonal domain enclosed with AC edges, and ZZB edged triangle may occur if the edges are terminated by hydrogen atoms. Further calculation shows that these edges present rich type-dependent properties and thus are important for various applications. This theoretical study showed that controlling the morphologies of BN domains and BN edges is crucial for various applications.

What are the active carbon species during graphene chemical vapor deposition growth?

The dissociation of carbon feedstock is a crucial step for understanding the mechanism of graphene chemical vapor deposition (CVD) growth. Using first-principles calculations, we performed a comprehensive theoretical study for the population of various active carbon species, including carbon monomers and various radicals, CHi (i = 1, 2, 3, 4), on four representative transition-metal surfaces, Cu(111), Ni(111), Ir(111) and Rh(111), under different experimental conditions. On the Cu surface, which is less active, the population of CH and C monomers at the subsurface is found to be very high and thus they are the most important precursors for graphene CVD growth. On the Ni surface, which is more active than Cu, C monomers at the subsurface dominate graphene CVD growth under most experimental conditions. In contrast, on the active Ir and Rh surfaces, C monomers on the surfaces are found to be very stable and thus are the main precursors for graphene growth. This study shows that the mechanism of graphene CVD growth depends on the activity of catalyst surfaces and the detailed graphene growth process at the atomic level can be controlled by varying the temperature or partial pressure of hydrogen.