Abstract: Nitrogen-doped graphene was recently synthesized and was reported to be a catalyst for hydrogen dissociative adsorption under a perpendicular applied electric field (F). In this work, the diffusion of H atoms on N-doped graphene, in the presence and absence of an applied perpendicular electric field, is studied using density functional theory. We demonstrate that the applied field can significantly facilitate the binding of hydrogen molecules on N-doped graphene through dissociative adsorption and diffusion on the surface. By removing the applied field the absorbed H atoms can be released efficiently. Our theoretical calculation indicates that N-doped graphene is a promising hydrogen storage material with reversible hydrogen adsorption/desorption where the applied electric field can act as a switch for the uptake/release processes.

Abstract: Quantum thermal transport in armchair and zigzag graphene nanoribbons is investigated in the presence of single atomic vacancies and subject to different boundary conditions. We start with a full comparison of the phonon polarizations and energy dispersions as given by a fifth-nearest-neighbor force-constant model (5NNFCM) and by elasticity theory of continuum membranes (ETCM). For free-edge ribbons, we discuss the behavior of an additional acoustic edge-localized flexural mode, known as fourth acoustic branch (4ZA), which has a small gap when it is obtained by the 5NNFCM. Then, we show that ribbons with supported edges have a sample-size dependent energy gap in the phonon spectrum which is particularly large for in-plane modes. Irrespective to the calculation method and the boundary condition, the dependence of the energy gap for the low-energy optical phonon modes against the ribbon width W is found to be proportional to 1/W for in-plane, and 1/W-2 for out-of-plane phonon modes. Using the 5NNFCM, the ballistic thermal conductance and its contributions from every single phonon mode are then obtained by the nonequilibrium Green's function technique. We found that, while edge and central localized single atomic vacancies do not affect the low-energy transmission function of in-plane phonon modes, they reduce considerably the contributions of the flexural modes. On the other hand, in-plane modes contributions are strongly dependent on the boundary conditions and at low temperatures can be highly reduced in supported-edge samples. These findings could open a route to engineer graphene based devices where it is possible to discriminate the relative contribution of polarized phonons and to tune the thermal transport on the nanoscale.

Abstract: The vortex-vortex interaction potential in bulk superconductors is calculated within the Ginzburg-Landau (GL) theory and is obtained from a numerical solution of a set of two coupled nonlinear GL differential equations for the vector potential and the superconducting order parameter, where the merger of vortices into a giant vortex is allowed. Further, the interaction potentials between a vortex and a giant vortex and between a vortex and an antivortex are obtained for both type-I and type-II superconductors. Our numerical results agree asymptotically with the analytical expressions for large intervortex separations that are available in the literature. We propose empirical expressions valid over the full interaction range, which are fitted to our numerical data for different values of the GL parameter.

Abstract: The scattering of a Gaussian wave packet in armchair and zigzag graphene edges is theoretically investigated by numerically solving the time-dependent Schrodinger equation for the tight-binding model Hamiltonian. Our theory allows us to investigate scattering in reciprocal space, and depending on the type of graphene edge we observe scattering within the same valley, or between different valleys. In the presence of an external magnetic field, the well-known skipping orbits are observed. However, our results demonstrate that in the case of a pseudomagnetic field, induced by nonuniform strain, the scattering by an armchair edge results in a nonpropagating edge state.

Abstract: The time evolution of a low-energy two-dimensional Gaussian wave packet in ABC-stacked n-layer graphene (ABC-NLG) is investigated. Expectation values of the position (x, y) of center-of-mass and the total probability densities of the wave packet are calculated analytically using the Green's function method. These results are confirmed using an alternative numerical method based on the split-operator technique within the Dirac approach for ABC-NLG, which additionally allows to include external fields and potentials. The main features of the zitterbewegung (trembling motion) of wave packets in graphene are demonstrated and are found to depend not only on the wave packet width and initial pseudospin polarization, but also on the number of layers. Moreover, the analytical and numerical methods proposed here allow to investigate wave packet dynamics in graphene systems with an arbitrary number of layers and arbitrary potential landscapes.

Abstract: We describe a semianalytic approach to the two-band Ginzburg-Landau theory, which predicts the behavior of vortices in two-band superconductors. We show that the character of the short-range vortex-vortex interaction is determined by the sign of the normal domain-superconductor interface energy, in analogy with the conventional differentiation between type I and type II superconductors. However, we also show that the long-range interaction is determined by a modified Ginzburg-Landau parameter κ*, different from the standard κ of a bulk superconductor. This opens the possibility for nonmonotonic vortex-vortex interaction, which is temperature dependent, and can be further tuned by alterations of the material on the microscopic scale.

Abstract: Semiconductors are the basis of many vital technologies such as electronics, computing, communications, optoelectronics, and sensing. Modern semiconductor technology can trace its origins to the invention of the point contact transistor in 1947. This demonstration paved the way for the development of discrete and integrated semiconductor devices and circuits that has helped to build a modern society where semiconductors are ubiquitous components of everyday life. A key property that determines the semiconductor electrical and optical properties is the bandgap. Beyond graphene, recently discovered two-dimensional (2D) materials possess semiconducting bandgaps ranging from the terahertz and mid-infrared in bilayer graphene and black phosphorus, visible in transition metal dichalcogenides, to the ultraviolet in hexagonal boron nitride. In particular, these 2D materials were demonstrated to exhibit highly tunable bandgaps, achieved via the control of layers number, heterostructuring, strain engineering, chemical doping, alloying, intercalation, substrate engineering, as well as an external electric field. We provide a review of the basic physical principles of these various techniques on the engineering of quasi-particle and optical bandgaps, their bandgap tunability, potentials and limitations in practical realization in future 2D device technologies.

Abstract: Atomically thin nanoribbons (NRs) have been at the forefront of materials science and nanoelectronics in recent years. State-of-the-art research on nanoscale materials has revealed that electronic, magnetic, phononic, and optical properties may differ dramatically when their one-dimensional forms are synthesized. The present article aims to review the recent advances in synthesis techniques and theoretical studies on NRs. The structure of the review is organized as follows: After a brief introduction to low dimensional materials, we review different experimental techniques for the synthesis of graphene nanoribbons (GNRs) with their advantages and disadvantages. In addition, theoretical investigations on width and edge-shape-dependent electronic and magnetic properties, functionalization effects, and quantum transport properties of GNRs are reviewed. We then devote time to the NRs of the transition metal dichalcogenides (TMDs) family. First, various synthesis techniques, E-field-tunable electronic and magnetic properties, and edge-dependent thermoelectric performance of NRs of MoS2 and WS2 are discussed. Then, strongly anisotropic properties, growth-dependent morphology, and the weakly width-dependent bandgap of ReS2 NRs are summarized. Next we discuss TMDs having a T-phase morphology such as TiSe2 and stable single layer NRs of mono-chalcogenides. Strong edge-type dependence on characteristics of GaS NRs, width-dependent Seebeck coefficient of SnSe NRs, and experimental analysis on the stability of ZnSe NRs are reviewed. We then focus on the most recently emerging NRs belonging to the class of transition metal trichalcogenides which provide ultra-high electron mobility and highly anisotropic quasi-1D properties. In addition, width-, edge-shape-, and functionalization-dependent electronic and mechanical properties of blackphosphorus, a monoatomic anisotropic material, and studies on NRs of group IV elements (silicene, germanene, and stanene) are reviewed. Observation of substrate-independent quantum well states, edge and width dependent properties, the topological phase of silicene NRs are reviewed. In addition, H-2 concentration-dependent transport properties and anisotropic dielectric function of GeNRs and electric field and strain sensitive I-V characteristics of SnNRs are reviewed. We review both experimental and theoretical studies on the NRs of group III-V compounds. While defect and N-termination dependent conductance are highlighted for boron nitride NRs, aluminum nitride NRs are of importance due to their dangling bond, electric field, and strain dependent electronic and magnetic properties. Finally, superlattice structure of NRs of GaN/AlN, Si/Ge, G/BN, and MoS2/WS2 is reviewed. Published by AIP Publishing.

Abstract: The mechanical properties of monolayer GaS and GaSe crystals are investigated in terms of their elastic constants: in-plane stiffness (C), Poisson ratio (nu), and ultimate strength (sigma(U)) by means of first-principles calculations. The calculated elastic constants are compared with those of graphene and monolayer MoS2. Our results indicate that monolayer GaS is a stiffer material than monolayer GaSe crystals due to the more ionic character of the Ga-S bonds than the Ga-Se bonds. Although their Poisson ratio values are very close to each other, 0.26 and 0.25 for GaS and GaSe, respectively, monolayer GaS is a stronger material than monolayer GaSe due to its slightly higher sU value. However, GaS and GaSe crystals are found to be more ductile and flexible materials than graphene and MoS2. We have also analyzed the band-gap response of GaS and GaSe monolayers to biaxial tensile strain and predicted a semiconductor-metal crossover after 17% and 14% applied strain, respectively, for monolayer GaS and GaSe. In addition, we investigated how the mechanical properties are affected by charging. We found that the flexibility of single layer GaS and GaSe displays a sharp increase under 0.1e/cell charging due to the repulsive interactions between extra charges located on chalcogen atoms. These charging-controllable mechanical properties of single layers of GaS and GaSe can be of potential use for electromechanical applications.

Abstract: By means of density functional theory based first-principles calculations, the structural, vibrational, and electronic properties of 1H- and 1T-phases of single-layer CaX2 (X = F, Cl, Br, or I) structures are investigated. Our results reveal that both the 1H- and 1T-phases are dynamically stable in terms of their phonon band dispersions with the latter being the energetically favorable phase for all single-layers. In both phases of single-layer CaX2 structures, significant phonon softening occurs as the atomic radius increases. In addition, each structural phase exhibits distinctive Raman active modes that enable one to characterize either the phase or the structure via Raman spectroscopy. The electronic band dispersions of single-layer CaX2 structures reveal that all structures are indirect bandgap insulators with a decrease in bandgaps from fluorite to iodide crystals. Furthermore, the calculated linear elastic constants, in-plane stiffness, and Poisson ratio indicate the ultra-soft nature of CaX2 single-layers, which is quite important for their nanoelastic applications. Overall, our study reveals that with their dynamically stable 1T- and 1H-phases, single-layers of CaX2 crystals can be alternative ultra-thin insulators.

Abstract: Using first-principles calculations, we study the dependence of the electronic and vibrational properties of multilayered PbI2 crystals on the number of layers and focus on the electronic-band structure and the Raman spectrum. Electronic-band structure calculations reveal that the direct or indirect semiconducting behavior of PbI2 is strongly influenced by the number of layers. We find that at 3L thickness there is a direct-to-indirect band gap transition (from bulk-to-monolayer). It is shown that in the Raman spectrum two prominent peaks, A(1g) and E-g, exhibit phonon hardening with an increasing number of layers due to the interlayer van der Waals interaction. Moreover, the Raman activity of the A(1g) mode significantly increases with an increasing number of layers due to the enhanced out-of-plane dielectric constant in the few-layer case. We further characterize rigid-layer vibrations of low-frequency interlayer shear (C) and breathing (LB) modes in few-layer PbI2. A reduced monatomic (linear) chain model (LCM) provides a fairly accurate picture of the number of layers dependence of the low-frequency modes and it is shown also to be a powerful tool to study the interlayer coupling strength in layered PbI2.

Abstract: The surface termination of MXenes greatly determines the electrochemical properties and ion kinetics on their surfaces. So far, hydroxyl-, oxygen-, and fluorine-terminated MXenes have been widely studied for energy storage applications. Recently, sulfur-functionalized MXene structures, which possess low diffusion barriers, have been proposed as candidate materials to enhance battery performance. We performed first-principles calculations on the structural, stability, electrochemical, and ion dynamic properties of Li-adsorbed sulfur-functionalized groups 3B, 4B, 5B, and 6B transition-metal (M)-based MXenes (i.e., M2CS2 with M = Sc, Ti, Zr, Hf, V, Nb, Ta, Cr, Mo, and W). We performed phonon calculations, which indicated that all of the above M2CS2 MXenes, except for Sc, are dynamically stable at T = 0 K. The ground-state structure of each M2CS2 monolayer depends on the type of M atom. For instance, while sulfur prefers to sit at the FCC site on Ti2CS2, it occupies the HCP site of Cr-based MXene. We determined the Li adsorption configurations at different concentrations using the cluster expansion method. The highest maximum open-circuit voltages were computed for the group 4B element (i.e., Ti, Zr, and Hf)-based M2CS2, which are larger than 2.1 V, while their average voltages are approximately 1 V. The maximum voltage for the group 6B element (i.e., Cr, Mo, W)-based M2CS2 is less than 1 V, and the average voltage is less than 0.71 V. We found that S functionalization is helpful for capacity improvements over the O-terminated MXenes. In this respect, the computed storage gravimetric capacity may reach up to 417.4 mAh/g for Ti2CS2 and 404.5 mAh/g for V2CS2. Ta-, Cr-, Mo-, and W-based M2CS2 MXenes show very low capacities, which are less than 100 mAh/g. The Li surface diffusion energy barriers for all of the considered MXenes are less than 0.22 eV, which is favorable for high charging and discharging rates. Finally, ab initio molecular dynamic simulations performed at 400 K and bond-length analysis with respect to Li concentration verify that selected promising systems are robust against thermally induced perturbations that may induce structural transformations or distortions and undesirable Li release.

Abstract: Using density functional theory, we obtain the mechanical and thermal properties of MX2 monolayers (where M = Cr, Mo, W and X = O, S, Se, Te). The C-centered phonon frequencies (i.e., A(1), A(2)'', E ', and E ''), relative frequency values of A(1), and E ' modes, and mechanical properties (i.e., elastic constants, Young modulus, and Poisson's ratio) display a strong dependence on the type of metal and chalcogenide atoms. In each chalcogenide (metal) group, transition-metal dichalcogenides (TMDCs) with W (O) atom are found to be much stiffer. Consistent with their stability, the thermal expansion of lattice constants for TMDCs with O (Te) is much slower (faster). Furthermore, in a heterostructure of these materials, the difference of the thermal expansion of lattice constants between the individual components becomes quite tiny over the whole temperature range. The calculated mechanical and thermal properties show that TMDCs are promising materials for heterostructures. (C) 2014 AIP Publishing LLC.

Abstract: The phase transition processes induced by ultrashort, 100 fs pulsed laser irradiation of Au, Cu, and Ni are studied by means of a combined atomistic-continuum approach. A moderately low absorbed laser fluence range, from 200 to 600 J/m2 is considered to study phase transitions by means of a local and a nonlocal order parameter. At low laser fluences, the occurrence of layer-by-layer evaporation has been observed, which suggests a direct solid to vapor transition. The calculated amount of molten material remains very limited under the conditions studied, especially for Ni. Therefore, our results show that a kinetic equation that describes a direct solid to vapor transition might be the best approach to model laser-induced phase transitions by continuum models. Furthermore, the results provide more insight into the applicability of analytical superheating theories that were implemented in continuum models and help the understanding of nonequilibrium phase transitions.

Abstract: The energy spectrum and density of states of graphene bilayer superlattices (SLs) are evaluated. We take into account doping and/or gating of the layers as well as tunnel coupling between them. In addition, we evaluate the transmission through such SLs and through single or double barriers. The transmission exhibits a strong dependence on the direction of the incident wave vector.