Abstract: Ripples in pristine freestanding graphene naturally orient themselves in an array that is alternately curved-up and curved-down; maintaining an average height of zero. Using scanning tunneling microscopy (STM) to apply a local force, the graphene sheet will reversibly rise and fall in height until the height reaches 60%-70% of its maximum at which point a sudden, permanent jump occurs. We successfully model the ripples as a spin-half Ising magnetic system, where the height of the graphene plays the role of the spin. The permanent jump in height, controlled by the tunneling current, is found to be equivalent to an antiferromagnetic-to-ferromagnetic phase transition. The thermal load underneath the STM tip alters the local tension and is identified as the responsible mechanism for the phase transition. Four universal critical exponents are measured from our STM data, and the model provides insight into the statistical role of graphene's unusual negative thermal expansion coefficient.

Abstract: Using first-principle calculations, we systematically study strain tuning of the electronic properties of two- dimensional graphitic carbon nitride nanosheets with empirical formula CnNm. We found the following: (i) the ferromagnetic ordered state in the metal-free systems (n, m) = (4,3), (10,9), and (14,12) remains stable in the presence of strain of about 6%. However, the system (9,7) loses its ferromagnetic ordering when increasing strain. This is due to the presence of topological defects in the (9,7) system, which eliminates the asymmetry between spin up and spin down of the p(z) orbitals when strain is applied. (ii) By applying uniaxial strain, a band gap opens in systems which are initially gapless. (iii) In semiconducting systems which have an initial gap of about 1 eV, the band gap is closed with applying uniaxial strain.

Abstract: Integration of graphene-based two-dimensional materials is essential for nanoelectronics applications. Using density-functional theory, we systematically investigate the electronic properties of vertically stacked graphene-graphitic carbon nitrides (GE/GCN). We also studied the covalently lateral stitched graphene-graphitic carbon nitrides (GE-GCN heterojunctions). The effects of perpendicular electric field on the electronic properties of six different heterostructures, i.e., (i) one layer of GE on top of a layer of CnNm with (n, m) = (3,1), (3,4), and (4,3) and (ii) three heterostructures CnNm/Cn'Nm', where (n, m) not equal (n', m') are elucidated. The most important calculated features are (i) the systems GE/C3N4, C3N/C3N4, GE-C3N, GE-C4N3, and C3N-C3N4 exhibit semiconducting characteristics having small band gaps of Delta(0)=20, 250, 100, 100, 80 meV, respectively while (ii) the systems GE/C4N3, C3N/C4N3, and C3N-C4N3 show ferromagnetic-metallic properties. In particular, we found that, in semiconducting heterostructures, the band gap increases nontrivially with increasing the absolute value of the applied perpendicular electric field. This work is useful for designing heterojunctions and heterostructures made of graphene and other two-dimensional materials such as those proposed in recent experiments [X. Liu and M. C. Hersam Sci. Adv. 5, 6444 (2019)].

Abstract: We investigate the stability of circular monolayer graphene subjected to a radial load using non-equilibrium molecular dynamics simulations. When monolayer graphene is radially stressed, after some small circular strain (~0.4%) it buckles and bends into a new bowl-like shape. Young's modulus is calculated from the linear relation between stress and strain before the buckling threshold, which is in agreement with experimental results. The prediction of elasticity theory for the buckling threshold of a radially stressed plate is presented and its results are compared to the one of our atomistic simulation. The Jarzynski equality is used to estimate the difference between the free energy of the non-compressed states and the buckled states. From a calculation of the free energy we obtain the optimum radius for which the system feels the minimum boundary stress.

Abstract: Using the valence force field model of Perebeinos and Tersoff (2009 Phys. Rev. B 79 241409(R)), different energy modes of suspended graphene subjected to tensile or compressive strain are studied. By carrying out Monte Carlo simulations it is found that: (i) only for small strains (vertical bar epsilon vertical bar (sic) 0.02) is the total energy symmetrical in the strain, while it behaves completely differently beyond this threshold; (ii) the important energy contributions in stretching experiments are stretching, angle bending, an out-of-plane term, and a term that provides repulsion against pi-pi misalignment; (iii) in compressing experiments the two latter terms increase rapidly, and beyond the buckling transition stretching and bending energies are found to be constant; (iv) from stretching-compressing simulations we calculated the Young's modulus at room temperature 350 +/- 3.15 N m(-1), which is in good agreement with experimental results (340 +/- 50 N m(-1)) and with ab initio results (322-353) N m(-1); (v) molar heat capacity is estimated to be 24.64 J mol(-1) K-1 which is comparable with the Dulong-Petit value, i. e. 24.94 J mol(-1) K-1, and is almost independent of the strain; (vi) nonlinear scaling properties are obtained from height-height correlations at finite temperature; (vii) the used valence force field model results in a temperature independent bending modulus for graphene, and (viii) the Gruneisen parameter is estimated to be 0.64.

Abstract: The strain-induced pseudomagnetic field in supported graphene deposited on top of a nanostructured substrate is investigated by using atomistic simulations. A step, an elongated trench, a one-dimensional barrier, a spherical bubble, a Gaussian bump, and a Gaussian depression are considered as support structures for graphene. From the obtained optimum configurations we found very strong induced pseudomagnetic fields which can reach up to similar to 1000 T due to the strain-induced deformations in the supported graphene. Different magnetic confinements with controllable geometries are found by tuning the pattern of the substrate. The resulting induced magnetic fields for graphene on top of a step, barrier, and trench are calculated. In contrast to the step and trench the middle part of graphene on top of a barrier has zero pseudomagnetic field. This study provides a theoretical background for designing magnetic structures in graphene by nanostructuring substrates. We found that altering the radial symmetry of the deformation changes the sixfold symmetry of the induced pseudomagnetic field.

Abstract: Atoms start behaving as waves rather than classical particles if confined in spaces commensurate with their de Broglie wavelength. At room temperature this length is only about one angstrom even for the lightest atom, hydrogen. This restricts quantum-confinement phenomena for atomic species to the realm of very low temperatures(1-5). Here, we show that van der Waals gaps between atomic planes of layered crystals provide angstrom-size channels that make quantum confinement of protons apparent even at room temperature. Our transport measurements show that thermal protons experience a notably higher barrier than deuterons when entering van der Waals gaps in hexagonal boron nitride and molybdenum disulfide. This is attributed to the difference in the de Broglie wavelengths of the isotopes. Once inside the crystals, transport of both isotopes can be described by classical diffusion, albeit with unexpectedly fast rates comparable to that of protons in water. The demonstrated angstrom-size channels can be exploited for further studies of atomistic quantum confinement and, if the technology can be scaled up, for sieving hydrogen isotopes.

Abstract: The influence of triaxial in-plane strain on the electronic properties of a hexagonal boron-nitride sheet is investigated using density functional theory. Different from graphene, the triaxial strain localizes the molecular orbitals of the boron-nitride flake in its center depending on the direction of the applied strain. The proposed technique for localizing the molecular orbitals that are close to the Fermi level in the center of boron nitride flakes can be used to actualize engineered nanosensors, for instance, to selectively detect gas molecules. We show that the central part of the strained flake adsorbs polar molecules more strongly as compared with an unstrained sheet.

Abstract: Large-scale atomistic simulations using the reactive empirical bond order force field approach is implemented to investigate thermal and mechanical properties of single-layer (SL) and multilayer (ML) molybdenum disulfide (MoS2). The amplitude of the intrinsic ripples of SL MoS2 are found to be smaller than those exhibited by graphene (GE). Furthermore, because of the van der Waals interaction between layers, the out-of-plane thermal fluctuations of ML MoS2 decreases rapidly with increasing number of layers. This trend is confirmed by the buckling transition due to uniaxial stress which occurs for a significantly larger applied tension as compared to graphene. For SL MoS2, the melting temperature is estimated to be 3700 K which occurs through dimerization followed by the formation of small molecules consisting of two to five atoms. When different types of vacancies are inserted in the SL MoS2 it results in a decrease of both the melting temperature as well as the stiffness.

Abstract: The electronic properties of boron nitride nanoribbons (BNNRs) doped with a line of carbon atoms are investigated using density functional calculations. By replacing a line of alternating B and N atoms with carbons, three different configurations are possible depending on the type of the atoms which bond to the carbons. We found very different electronic properties for these configurations: (i) the NCB arrangement is strongly polarized with a large dipole moment having an unexpected direction, (ii) the BCB and NCN arrangements are nonpolar with zero dipole moment, (iii) the doping by a carbon line reduces the band gap regardless of the local arrangement of the borons and the nitrogens around the carbon line, and (iv) the polarization and energy gap of the carbon-doped BNNRs can be tuned by an electric field applied parallel to the carbon line. Similar effects were found when either an armchair or zigzag line of carbon was introduced.

Abstract: The buckling of defected rectangular graphene nanoribbons when subjected to axial stress with supported boundary conditions is investigated using atomistic simulations. The buckling strain and mechanical stiffness of monolayer graphene decrease with the percentage of randomly distributed vacancies. The elasticity to plasticity transition in the stress-strain curve, at low percentage of vacancies, are found to be almost equal to the buckling strain thresholds and they decrease with increasing percentage of vacancies.

Abstract: The electronic properties of a triaxially strained hexagonal graphene flake with either armchair or zigzag edges are investigated using molecular dynamics simulations and tight-binding calculations. We found that (i) the pseudomagnetic field in strained graphene flakes is not uniform neither in the center nor at the edge of zigzag terminated flakes, (ii) the pseudomagnetic field is almost zero in the center of armchair terminated flakes but increases dramatically near the edges, (iii) the pseudomagnetic field increases linearly with strain, for strains lower than 15% but increases nonlinearly beyond it, (iv) the local density of states in the center of the zigzag hexagon exhibits pseudo-Landau levels with broken sublattice symmetry in the zeroth pseudo-Landau level, and in addition there is a shift in the Dirac cone due to strain induced scalar potentials, and (v) there is size effect in pseudomagnetic field. This study provides a realistic model of the electronic properties of inhomogeneously strained graphene where the relaxation of the atomic positions is correctly included together with strain induced modifications of the hopping terms up to next-nearest neighbors.

Abstract: We report subnanometer, high-bandwidth measurements of the out-of-plane (vertical) motion of atoms in freestanding graphene using scanning tunneling microscopy. By tracking the vertical position over a long time period, a 1000-fold increase in the ability to measure space-time dynamics of atomically thin membranes is achieved over the current state-of-the-art imaging technologies. We observe that the vertical motion of a graphene membrane exhibits rare long-scale excursions characterized by both anomalous mean-squared displacements and Cauchy-Lorentz power law jump distributions.

Abstract: The spatial dependence of the van der Waals (vdW) energy between graphene and hexagonal boron-nitride (h-BN) is investigated using atomistic simulations. The van der Waals energy between graphene and h-BN shows a hexagonal superlattice structure identical to the observed Moire pattern in the local density of states, which depends on the lattice mismatch and misorientation angle between graphene and h-BN. Our results provide atomistic features of the weak van der Waals interaction between graphene and BN which are in agreement with experiment and provide an analytical expression for the size of the spatial variation of the weak van der Waals interaction. We also found that the A-B-lattice symmetry of graphene is broken along the armchair direction. (C) 2014 AIP Publishing LLC.