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

Metal organic chemical vapor deposition (MOCVD) of group III nitrides on graphene heterostructures offers new opportunities for the development of flexible optoelectronic devices and for the stabilization of conceptually-new two-dimensional materials. However, the MOCVD of group III nitrides is regulated by an intricate interplay of gas-phase and surface reactions that are beyond the resolution of experimental techniques. We use density-functional ab initio molecular dynamics (AIMD) with van der Waals corrections to identify atomistic pathways and associated electronic mechanisms driving precursor/surface reactions during metal organic vapor phase epitaxy at elevated temperatures of aluminum nitride on graphene, considered here as model case study. The results presented provide plausible interpretations of atomistic and electronic processes responsible for delivery of Al, C adatoms, and C-Al, CHx, AlNH2 admolecules on pristine graphene via precursor/surface reactions. In addition, the simulations reveal C adatom permeation across defect-free graphene, as well as exchange of C monomers with graphene carbon atoms, for which we obtain rates of ∼0.3 THz at typical experimental temperatures (1500 K), and extract activation energies E = 0.28 ± 0.13 eV and attempt frequencies Aexc = 2.1 (×1.7±1) THz via Arrhenius linear regression. The results demonstrate that AIMD simulations enable understanding complex precursor/surface reaction mechanisms, and thus propose AIMD to become an indispensable routine prediction-tool toward more effective exploitation of chemical precursors and better control of MOCVD processes during synthesis of functional materials.

Van der Waals stacks of few-layer h-AlN with graphene: an ab initio study of structural, interaction and electronic properties.

Graphite-like hexagonal AlN (h-AlN) multilayers have been experimentally manifested and theoretically modeled. The development of any functional electronics applications of h-AlN would most certainly require its integration with other layered materials, particularly graphene. Here, by employing vdW-corrected density functional theory calculations, we investigate structure, interaction energy, and electronic properties of van der Waals stacking sequences of few-layer h-AlN with graphene. We find that the presence of a template such as graphene induces enough interlayer charge separation in h-AlN, favoring a graphite-like stacking formation. We also find that the interface dipole, calculated per unit cell of the stacks, tends to increase with the number of stacked layers of h-AlN and graphene.

Tuning band inversion symmetry of buckled III-Bi sheets by halogenation.

First-principles calculations are employed to investigate structural, electronic and topological insulating properties of XBi (X = B, Al, Ga, and In) monolayers upon halogenation. It is known that Y-XBi (X = Ga, In, Tl; Y = F, Cl, Br, I) can originate inversion-asymmetric topological insulators with large bulk band gaps. Our results suggest that Y-XBi (X = B, Al; Y = F, Cl, Br, I) may also result in nontrivial topological insulating phases. Despite the lower atomic number of B and Al, the spin-orbit coupling opens a band gap of about 400 meV in Y-XBi (X = B, Al), exhibiting an unusual electronic behavior for practical applications in spintronics. The nature of the bulk band gap and Dirac-cone edge states in their nanoribbons depends on the group-III elements and Y chemical species. They lead to a chemical tunability, giving rise to distinct band inversion symmetries and exhibiting Rashba-type spin splitting in the valence band of these systems. These findings indicate that a large family of Y-XBi sheets can exhibit nontrivial topological characteristics, by a proper tuning, and open a new possibility for viable applications at room temperature.

Spin-orbit-induced gap modification in buckled honeycomb XBi and XBi3 (X = B, Al, Ga, and In) sheets.

The band structure and stability of XBi and XBi3 (X = B, Al, Ga, and In) single sheets are predicted using first-principles calculations. It is demonstrated that the band gap values of these new classes of two-dimensional (2D) materials depend on both the spin-orbit coupling (SOC) and type of group-III elements in these hetero-sheets. Thus, topological properties can be achieved, allowing for viable applications based on coherent spin transport at room temperature. The spin-orbit effects are proved to be essential to explain the tunability by group-III atoms. A clear effect of including SOC in the calculations is lifting the spin degeneracy of the bands at the Γ point of the Brillouin zone. The nature of the band gaps, direct or indirect, is also tuned by SOC, and by the appropriate X element involved. It is observed that, in the case of XBi single sheets, band inversions naturally occur for GaBi and InBi, which exhibit band gap values around 172 meV. This indicates that these 2D materials are potential candidates for topological insulators. On the contrary, a similar type of band inversion, as obtained for the XBi, was not observed in the XBi3 band structure. In general, the calculations, taking into account SOC, reveal that some of these buckled sheets exhibit sizable gaps, making them suitable for applications in room-temperature spintronic devices.

Feasibility of novel (H3C)nX(SiH3)3-n compounds (X = B, Al, Ga, In): structure, stability, reactivity, and Raman characterization from ab initio calculations.

We employ ab initio calculations to predict the equilibrium structure, stability, reactivity, and Raman scattering properties of sixteen different (H3C)nX(SiH3)3-n compounds (X = B, Al, Ga, In) with n = 0-3. Among this methylsilylmetal family, only the (H3C)3X members, i.e., trimethylboron (TMB), trimethylaluminum (TMA), trimethylgallium (TMG), and trimethylindium (TMI), are currently well-studied. The remaining twelve compounds proposed here open up a two-dimensional array of new possibilities for precursors in various deposition processes, and evoke potential applications in the chemical synthesis of other compounds. We infer that within the (H3C)nX(SiH3)3-n family, the compounds with fewer silyl groups (and consequently with more methyl groups) are less reactive and more stable. This trend is verified from the calculated cohesive energy, Gibbs free energy of formation, bond strength, and global chemical indices. Furthermore, we propose sequential reaction routes for the synthesis of (H3C)nX(SiH3)3-n by substitution of methyl by silyl groups, where the silicon source is the silane gas. The corresponding reaction barriers for these chemical transformations lie in the usual energy range typical for MOCVD processes. We also report the Raman spectra and light scattering properties of the newly proposed (H3C)nX(SiH3)3-n compounds, in comparison with available data of known members of this family. Thus, our computational experiment provides useful information for a systematic understanding of the stability/reactivity and for the identification of these compounds.

Quantum size effects in the polarizability of carbon fullerenes.

We investigate the size-dependent dielectric response of carbon fullerenes with up to 3840 atoms in the framework of the linear response theory. Our results suggest a significant polarizability enhancement due to quantum size effects with respect to classical or semiclassical calculations. The accuracy of our results, based on a parametrized Hamiltonian, is verified by ab initio time dependent density functional calculations for smaller fullerenes. Our findings underline the importance of quantum effects in the electronic response of nano- and mesoscopic systems.

Structural identification of metcars.

The ground-state structure of the metcar Ti8C12 is investigated using first-principles computer simulations. Comparison with recent experimental data on the vibrational spectrum of gas phase Ti8C12 allows one to identify the geometrical structure of the clusters studied in the experiment. The present combination of predictive first-principles computer simulations and detailed experimental measurements of the vibrational spectra of clusters offers the first viable tool for structural identification of cluster shapes.