Abstract: The thermodynamic properties of the electron gas in multilayer graphene depend strongly on the number of layers and the type of stacking. Here we analyze how those properties change when we vary the number of layers for rhombohedral stacked multilayer graphene and compare our results with those from a conventional two-dimensional electron gas. We show that the highly degenerate zero-energy Landau level which is partly filled with electrons and partly with holes has a strong influence on the values of the different thermodynamic quantities.

Abstract: Using highly efficient GPU-based simulations of the tight-binding Bogoliubov-de Gennes equations we solve self-consistently for the pair correlation in rhombohedral (ABC) and Bernal (ABA) multilayer graphene by considering a finite intrinsic s-wave pairing potential. We find that the two different stacking configurations have opposite bulk/surface behavior for the order parameter. Surface superconductivity is robust for ABC stacked multilayer graphene even at very low pairing potentials for which the bulk order parameter vanishes, in agreement with a recent analytical approach. In contrast, for Bernal stacked multilayer graphene, we find that the order parameter is always suppressed at the surface and that there exists a critical value for the pairing potential below which no superconducting order is achieved. We considered different doping scenarios and find that homogeneous doping strongly suppresses surface superconductivity while nonhomogeneous field-induced doping has a much weaker effect on the superconducting order parameter. For multilayer structures with hybrid stacking (ABC and ABA) we find that when the thickness of each region is small (few layers), high-temperature surface superconductivity survives throughout the bulk due to the proximity effect between ABC/ABA interfaces where the order parameter is enhanced. DOI: 10.1103/PhysRevB.87.134509

Abstract: Monolayer alpha-graphyne was recently proposed as a new all-carbon material having an electronic spectrum consisting of Dirac cones. Based on a first-principles investigation of bilayer alpha-graphyne, we show that the electronic band structure is qualitatively different from its monolayer form and depends crucially on the stacking mode of the two layers. Two stable stacking modes are found: a configuration with a gapless parabolic band structure, similar to AB stacked bilayer graphene, and another one which exhibits a doubled Dirac-cone spectrum. The latter can be tuned by an electric field with a gap opening rate of 0.3 eA. (C) 2013 AIP Publishing LLC.

Abstract: The interaction between electrons and plasmons in trilayer graphene is investigated within the Overhauser approach resulting in the “plasmaron” quasiparticle. This interaction is cast into a field theoretical problem, and its effect on the energy spectrum is calculated using improved Wigner-Brillouin perturbation theory. The plasmaron spectrum is shifted with respect to the bare electron spectrum by ΔE(k)∼150−200meV for ABC stacked trilayer graphene and for ABA trilayer graphene by ΔE(k)∼30−150 meV[ ΔE(k) ∼1 −5meV] for the hyperbolic (linear) part of the spectrum. The shift in general increases with the electron concentration and electron momentum. The dispersion of plasmarons is more pronounced in ABC stacked than in ABA stacked trilayer graphene, because of the different energy band structure and their different plasmon dispersion.

Abstract: The diffusive motion of metal nanoparticles Au and Ag on monolayer and between bilayer heterostructures of transition metal dichalcogenides and graphene are investigated in the framework of density functional theory. We found that the minimum energy barriers for diffusion and the possibility of cluster formation depend strongly on both the type of nanoparticle and the type of monolayers and bilayers. Moreover, the tendency to form clusters of Ag and Au can be tuned by creating various bilayers. Tunability of the diffusion characteristics of adatoms in van der Waals heterostructures holds promise for controllable growth of nanostructures. (C) 2014 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution 3.0 Unported License.

Abstract: Molecular-dynamic simulations were performed in order to investigate the melting processes of isotropically confined binary systems. We considered two species of particles, which differ by their amount of electric charge. A Lindemann type of criterion was used to determine the angular melting temperature. We demonstrate that the magic-to-normal cluster transition can evolve in two distinct ways, that is, through a structural phase transition of the first order or via a smooth transition where an increase of the shells' width leads to a continuous decreasing mechanical stability of the system. Moreover, for large systems, we demonstrate that the internal cluster exerts a minor effect on the mechanical stability of the external shell. Furthermore, we show that highly symmetric configurations, such as those found for multiple ring structures, have large mechanical stability, i.e., high angular melting temperature.

Abstract: Transition metal dichalcogenides (TMDs) are considered as promising candidates for next generation of electronic and optoelectronic devices. To make use of these materials, for instance in field effect transistor applications, it is mandatory to know the detailed properties of contacts of such TMDs with metal electrodes. Here, we investigate the role of the contact structure on the electronic and transport properties of metal-MoSe2 interfaces. Two different contact types, namely face and edge contacts, are studied. We consider both low (Sc) and high (Au) work function metals in order to thoroughly elucidate the role of the metal work function and the type of metal. First principles plane wave calculations and transport calculations based on nonequilibrium Green's function formalism reveal that the contact type has a large impact on the electronic and transport properties of metal-MoSe2 interfaces. For the Sc electrode, the Schottky barrier heights are around 0.25 eV for face contact and bigger than 0.6 eV for edge contact. For the Au case, we calculate very similar barrier heights for both contact types with an average value of 0.5 eV. Furthermore, while the face contact is found to be highly advantageous as compared to the edge contact for the Sc electrode, the latter contact becomes a better choice for the Au electrode. Our findings provide guidelines for the fabrication of TMD-based devices.

Abstract: The absence of a direct-to-indirect band gap transition in ReS2 when going from the monolayer to bulk makes it special among the other semiconducting transition metal dichalcogenides. The functionalization of this promising layered material emerges as a necessity for the next generation technological applications. Here, the structural, electronic, and magnetic properties of substitutionally doped ReS2 monolayers at either the S or Re site were systematically studied by using first principles density functional calculations. We found that substitutional doping of ReS2 depends sensitively on the growth conditions of ReS2. Among the large number of non-metallic atoms, namely H, B, C, Se, Te, F, Br, Cl, As, P. and N, we identified the most promising candidates for n-type and p-type doping of ReS2. While Cl is an ideal candidate for n-type doping, P appears to be the most promising candidate for p-type doping of the ReS2 monolayer. We also investigated the doping of ReS2 with metal atoms, namely Mo, W, Ti, V. Cr, Co, Fe, Mn, Ni, Cu, Nb, Zn, Ru, Os and Pt. Mo, Nb, Ti, and V atoms are found to be easily incorporated in a single layer of ReS2 as substitutional impurities at the Re site for all growth conditions considered in this work. Tuning chemical potentials of dopant atoms energetically makes it possible to dope ReS2 with Fe, Co, Cr, Mn, W, Ru, and Os at the Re site. We observe a robust trend for the magnetic moments when substituting a Re atom with metal atoms such that depending on the electronic configuration of dopant atoms, the net magnetic moment of the doped ReS2 becomes either 0 or 1 mu(B). Among the metallic dopants, Mo is the best candidate for p-type doping of ReS2 owing to its favorable energetics and promising electronic properties.

Abstract: Phase-slip centers/lines and hot spots are the main mechanisms for dissipation in current-carrying superconducting thin films. The pulsed-current method has recently been shown to be an effective tool in studying the dynamics of phase-slip centers and their evolution to hot spots. We use the time-dependent Ginzburg-Landau theory in the study of the dynamics of the superconducting condensate in superconducting strips under external current and zero external magnetic field. We show that both the flux-flow state (i.e., slow-moving vortices) and the phase-slip line state (i.e., fast-moving vortices) are dynamically stable dissipative units with temperature smaller than the critical one, whereas hot spots, which are localized normal regions where the local temperature exceeds the critical value, expand in time, resulting ultimately in a complete destruction of the condensate. The response time of the system to abrupt switching on of the overcritical current decreases with increasing both the value of the current (at all temperatures) and temperature (for a given value of the applied current). Our results are in good qualitative agreement with experiments we have conducted on Nb thin strips.

Abstract: The electronic and optical properties of zincblende nanowires are investigated in the presence of a uniform magnetic field directed along the [001] growth direction within the k . p method. We focus our numerical study on core-shell nanowires consisting of the III-V materials GaAs, AlxGa1-xAs and AlyGa1-y/0.51In0.49P. Nanowires with electrons confined in the core exhibit a Fock-Darwin-like spectrum, whereas nanowires with electrons confined in the shell show Aharonov-Bohm oscillations. Thus, by properly choosing the core and the shell materials of the nanowire, the optical properties in a magnetic field can be tuned in very different ways.

Abstract: Metallic contacts are critical components of electronic devices and the presence of a large Schottky barrier is detrimental for an optimal device operation. Here, we show by using first-principles calculations that a self-assembled monolayer (SAM) of polar molecules between the metal electrode and MoSe2 monolayer is able to convert the Schottky contact into an almost Ohmic contact. We choose -CH3 and -CF3 terminated short-chain alkylthiolate (i.e. SCH3 and fluorinated alkylthiolates (SCF3)) based SAMs to test our approach. We consider both high (Au) and low (Sc) work function metals in order to thoroughly elucidate the role of the metal work function. In the case of Sc, the Fermi level even moves into the conduction band of the MoSe2 monolayer upon SAM insertion between the metal surface and the MoSe2 monolayer, and hence possibly switches the contact type from Schottky to Ohmic. The usual Fermi level pinning at the metal-transition metal dichalcogenide (TMD) contact is shown to be completely removed upon the deposition of a SAM. Systematic analysis indicates that the work function of the metal surface and the energy level alignment between the metal electrode and the TMD monolayer can be tuned significantly by using SAMs as a buffer layer. These results clearly indicate the vast potential of the proposed interface engineering to modify the physical and chemical properties of MoSe2.

Abstract: We provide experimental evidence for the existence of ferromagnetism in bilayers of Pd/C-60 which is supported by theoretical calculations based on density functional theory (DFT). The observed ferromagnetism is surprising as C-60 and Pd films are both non-ferromagnetic in the non-interacting limit. Magnetization (M) versus applied field (H) data acquired at different temperatures (T) show magnetic hysteresis with typical coercive fields (H-c) on the order of 50 Oe. From the temperature-dependent magnetization M(T) we extract a Curie temperature (T-c >= 550 K) using Bloch-like power law extrapolations to high temperatures. Using DFT calculations we investigated all plausible scenarios for the interaction between the C-60 molecules and the Pd slabs, Pd single atoms and Pd clusters. DFT shows that while the C-60 molecules are nonmagnetic, Pd films have a degenerate ground state that subject to a weak perturbation, can become ferromagnetic. Calculations also show that the interaction of C-60 molecules with excess Pd atoms and with sharp edges of a Pd slab is the most likely configuration that render the system ferromagnetic Interestingly, the calculated charge transfer (0.016 e per surface Pd atom, 0.064 e per Pd for intimate contact region) between C-60 and Pd does not appear to play an important role. (C) 2013 Elsevier B.V. All rights reserved.

Abstract: Recently, rhenium disulfide (ReS2) monolayers were experimentally extracted by conventional mechanical exfoliation technique from as-grown ReS2 crystals. Unlike the well-known members of transition metal dichalcogenides (TMDs), ReS2 crystallizes in a stable distorted-1T structure and lacks an indirect to direct gap crossover. Here we present an experimental and theoretical study of the formation, energetics, and stability of the most prominent lattice defects in monolayer ReS2. Experimentally, irradiation with 3-MeV He+2 ions was used to break the strong covalent bonds in ReS2 flakes. Photoluminescence measurements showed that the luminescence from monolayers is mostly unchanged after highly energetic a particle irradiation. In order to understand the energetics of possible vacancies in ReS2 we performed systematic first-principles calculations. Our calculations revealed that the formation of a single sulfur vacancy has the lowest formation energy in both Re and S rich conditions and a random distribution of such defects are energetically more preferable. Sulfur point defects do not result in any spin polarization whereas the creation of Re-containing point defects induce magnetization with a net magnetic moment of 1-3 mu B. Experimentally observed easy formation of sulfur vacancies is in good agreement with first-principles calculations.

Abstract: We present a systematic study of the energy spectra of graphene quantum rings having different geometries and edge types in the presence of a perpendicular magnetic field. Results are obtained within the tight-binding (TB) and Dirac models and we discuss which features of the former can be recovered by using the approximations imposed by the latter. Energy levels of graphene quantum rings obtained by diagonalizing the TB Hamiltonian are demonstrated to be strongly dependent on the rings geometry and the microscopical structure of the edges. This makes it difficult to recover those spectra by the existing theories that are based on the continuum (Dirac) model. Nevertheless, our results show that both approaches (i.e., TB and Dirac model) may provide similar results, but only for very specific combinations of ring geometry and edge types. The results obtained by a simplified model describing an infinitely thin circular Dirac ring show good agreement with those obtained for hexagonal and rhombus armchair graphene rings within the TB model. Moreover, we show that the energy levels of a circular quantum ring with an infinite mass boundary condition obtained within the Dirac model agree with those for a ring defined by a ring-shaped staggered potential obtained within the TB model.

Abstract: Moire patterns in the pseudo-magnetic field and in the strain profile of graphene (GE) when put on top of a hexagonal lattice substrate are predicted from elasticity theory. The van der Waals interaction between GE and the substrate induces out-of-plane deformations in graphene which results in a strain field, and consequently in a pseudo-magnetic field. When the misorientation angle is about 0.5 degrees, a three-fold symmetric strain field is realized that results in a pseudo-magnetic field very similar to the one proposed by F. Guinea, M. I. Katsnelson, and A. K. Geim [Nature Phys. 6, 30 (2010)]. Our results show that the periodicity and length of the pseudo-magnetic field can be tuned in GE by changing the misorientation angle and substrate adhesion parameters and a considerable energy gap (23 meV) can be obtained due to out-of-plane deformation of graphene which is in the range of recent experimental measurements (20-30 meV). (C) 2014 AIP Publishing LLC.