Abstract: We consider a classical system of two-dimensional (2D) charged particles, interacting through a repulsive Yukawa potential exp(-r/λ)/r, and confined in a parabolic channel that limits the motion of the particles in the y direction. Along the x direction, the particles are subject to a periodic potential. The ground-state configurations and the normal-mode spectra of the system are obtained as a function of the periodicity and strength of the periodic potential (V0) and density. An interesting set of tunable ground-state configurations are found, with first- or second-order structural transitions between them. A configuration with particles aligned, perpendicular to the x direction, in each minimum of the periodic potential is obtained for V0 larger than some critical value that has a power-law dependence on the density. The phonon spectrum of different configurations was also calculated. A localization of the modes into a small frequency interval is observed for sufficiently large strength of the periodic potential, and a tunable gap in the phonon spectrum is found as a function of V0.

Abstract: The properties of a new class of phosphorene nanotubes (PNT) are investigated by performing first-principles calculations. We demonstrate that it is advantageous to use blue phosphorene in order to make small nanotubes and propose a way to create low-energy PNTs by the inclusion of defect lines. Five different types of defect lines are investigated and incorporated in various combinations. The resulting defect-induced faceted PNTs have negligible bending stresses which leads to a reduction in the formation energy with respect to round PNTs. Our armchair faceted PNTs have similar formation energies than the recently proposed multiphase faceted PNTs, but they have a larger variety of possible structures. Our zigzag faceted PNTs have lower formation energies than round tubes and multiphase faceted nanotubes. The electronic properties of the defect-induced faceted PNTs are determined by the defect lines which control the band gap and the shape of the electronic states at the band edges. These band gaps increase with the radius of the nanotubes and converge to those of isolated defect lines.

Abstract: Using the tight-binding approach, we study the influence of a perpendicular magnetic field on the energy levels of hexagonal, triangular, and circular bilayer graphene (BLG) quantum dots (QDs) with zigzag and armchair edges. We obtain the energy levels for AB (Bernal)-stacked BLG QDs in both the absence and the presence of a perpendicular electric field (i.e., biased BLG QDs). We find different regions in the spectrum of biased QDs with respect to the crossing point between the lowest-electron and -hole Landau levels of a biased BLG sheet. Those different regions correspond to electron states that are localized at the center, edge, or corner of the BLG QD. Quantum Hall corner states are found to be absent in circular BLG QDs. The spatial symmetry of the carrier density distribution is related to the symmetry of the confinement potential, the position of zigzag edges, and the presence or absence of interlayer inversion symmetry.

Abstract: When subjected to the interfacially induced Dzyaloshinskii-Moriya interaction, the ground state in thin ferromagnetic films with high perpendicular anisotropy is cycloidal. The period of this cycloidal state depends on the strength of the Dzyaloshinskii-Moriya interaction. In this work, we have studied the effect of confinement on the magnetic ground state and excited states, and we determined the phase diagram of thin strips and thin square platelets by means of micromagnetic calculations. We show that multiple cycloidal states with different periods can be stable in laterally confined films, where the period of the cycloids does not depend solely on the Dzyaloshinskii-Moriya interaction strength but also on the dimensions of the film. The more complex states comprising skyrmions are also found to be stable, though with higher energy.

Abstract: We calculate the excitonic spectrum of few-layer black phosphorus by direct diagonalization of the effective mass Hamiltonian in the presence of an applied in-plane electric field. The strong attractive interaction between electrons and holes in this system allows one to investigate the Stark effect up to very high ionizing fields, including also the excited states. Our results show that the band anisotropy in black phosphorus becomes evident in the direction-dependent field-induced polarizability of the exciton.

Abstract: The structural transitions of the ground state of a system of repulsively interacting particles confined in a quasi-one-dimensional channel, and the effect of the interparticle interaction as well as the functional form of the confinement potential on those transitions are investigated. Although the nonsequential ordering of transitions (non-SOT), i.e., the 1 – 2 – 4 – 3 – 4 – 5 – 6 – ... sequence of chain configurations with increasing density, is widely robust as predicted in a number of theoretical studies, the sequential ordering of transitions (SOT), i.e., the 1 – 2 – 3 – 4 – 5 – 6 – ... chain, is found as the ground state for long-ranged interparticle interaction and hard-wall-like confinement potentials. We found an energy barrier between every two different phases around its transition point, which plays an important role in the preference of the system to follow either a SOT or a non-SOT. However, that preferential transition requires also the stability of the phases during the transition. Additionally, we analyze the effect of a small structural disorder on the transition between the two phases around its transition point. Our results show that a small deformation of the triangular structure changes dramatically the picture of the transition between two phases, removing in a considerable region the non-SOT in the system. This feature could explain the fact that the non-SOT is, up to now, not observed in experimental systems, and suggests a more advanced experimental setup to detect the non-SOT.

Abstract: Using first-principles calculations and reactive force field molecular dynamics simulations, we study the structural properties and dynamics of a fluorine (F) atom, either adsorbed on the surface of single layer graphene (F/GE) or between the layers of AB stacked bilayer graphene (F@ bilayer graphene). It is found that the diffusion of the F atom is very different in those cases, and that the mobility of the F atom increases by about an order of magnitude when inserted between two graphene layers. The obtained diffusion constant for F/GE is twice larger than that experimentally found for gold adatom and theoretically found for C-60 molecule on graphene. Our study provides important physical insights into the dynamics of fluorine atoms between and on graphene layers and explains the mechanism behind the separation of graphite layers due to intercalation of F atoms.

Abstract: Simulations are carried out based on the dynamical mean-field theory (DMFT) in order to investigate the properties of correlated thin films for various values of the chemical potential, temperature, interaction strength, and applied transverse electric field. Application of a sufficiently strong field to a thin film at half filling leads to the appearance of conducting regions near the surfaces of the film, whereas in doped slabs the application of a field leads to a conductivity enhancement on one side of the film and a gradual transition to the insulating state on the opposite side. In addition to the inhomogeneous DMFT, a local density approximation (LDA) is considered in which the particle density n, quasiparticle residue Z, and spectral weight at the Fermi level A(ω=0) of each layer are approximated by a homogeneous bulk environment. A systematic comparison between the two approaches reveals that the less expensive LDA results are in good agreement with the DMFT approach, except close to the metal-to-insulator transition points and in the layers immediately at the film surfaces. LDA values for n are overall more reliable than those for Z and A(ω=0). The hysteretic behavior (memory effect) characteristic of the bulk doping driven Mott transition persists in the slab.

Abstract: We reveal the emergence of metallicKondo clouds in an atomistic carbon quantum dot, realized as a single-atom junction in a suitably patterned graphene nanoflake. Using density functional dynamical mean-field theory (DFDMFT) we show how correlation effects lead to striking features in the electronic structure of our device, and how those are enhanced by the electron-electron interactions when graphene is patterned at the atomistic scale. Our setup provides a well-controlled environment to understand the principles behind the orbital-selective Kondo physics and the interplay between orbital and spin degrees of freedom in carbon-based nanomaterials, which indicate new pathways for spintronics in atomically patterned graphene.

Abstract: Often real samples of graphene consist of islands of both monolayer and bilayer graphene. Bound states in such hybrid quantum dots are investigated for (i) a circular single-layer graphene quantum dot surrounded by an infinite bilayer graphene sheet and (ii) a circular bilayer graphene quantum dot surrounded by an infinite single-layer graphene. Using the continuum model and applying zigzag boundary conditions at the single-layer-bilayer graphene interface, we obtain analytical results for the energy levels and the corresponding wave spinors. Their dependence on perpendicular magnetic and electric fields are studied for both types of quantum dots. The energy levels exhibit characteristics of interface states, and we find anticrossings and closing of the energy gap in the presence of a bias potential.

Abstract: Nitrogenated holey graphene (NHG) is a recently synthesized two-dimensional material. In this paper the structural and electronic properties of heterostructures of graphene and NHG are investigated using first-principles and tight-binding calculations. Due to the lattice mismatch between NHG and graphene, the formation of a moire pattern is preferred in the graphene/NHG heterostructure, instead of a lattice-coherent structure. In moire-patterned graphene/NHG, the band gap opening at the K point is negligible, and the linear band dispersion of graphene survives. Applying an electric field modifies the coupling strength between the two atomic layers. The Fermi velocity upsilon(F) is reduced as compared to the one of pristine graphene, and its magnitude depends on the twist angle theta between graphene and NHG: For theta = 0 degrees, upsilon(F) is 30% of that of graphene, and it increases rapidly to a value of 80% with increasing theta. The heterostructure exhibits electron-hole asymmetry in upsilon(F), which is large for small theta. In NHG encapsulated between two graphene layers, a “Dirac ring” appears around the K point. Its presence is robust with respect to the relative stacking of the two graphene layers. These findings can be useful for future applications of graphene/NHG heterostructures.

Abstract: The Majorana zero modes and persistent spin current in mesoscopic d-wave-superconducting loops with spin-orbit (SO) interaction are investigated by numerically solving the spin-generalized Bogoliubov-de Gennes equations self-consistently. For some appropriate strength of the SO coupling, Majorana zero-energy states and sharp jumps of the spin-polarized currents can be observed when the highest energy levels cross the Fermi energy in the spectrum, leading to spin currents with opposite chirality flowing near the inner and outer edges of the sample. When the threaded magnetic flux turns on, four flux-dependent patterns of the persistent spin current with step-like features show up, accompanied by Majorana edge modes at flux values where the energy gap closes. Moreover, the Majorana zero mode is highly influenced by the direction of the Zeeman field. A finite in-plane field can lead to the gap opening since the inversion symmetry is broken. Remarkably, multiple Majorana zero-energy states occur in the presence of an out-of-plane field h(z), and the number of steps in the spin current evolution can be effectively tuned by the field strength due to the shift of Majorana zero modes. Finally, when the loop sample contains surface indentation defects, zero-energy modes can always show up in the presence of an appropriate h(z). Interestingly, multiple Majorana states may be present in the system with a corner defect even if h(z) = 0.

Abstract: The application of graphene to electronic and optoelectronic devices is limited by the absence of reliable semiconducting variants of this material. A promising candidate in this respect is graphene oxide, with a band gap on the order of similar to 5 eV, however, this has a finite density of states at the Fermi level. Here, we examine the electronic structure of three variants of half-fluorinated carbon on Sic(0001), i.e., the (6 root 3 x 6 root 3) R30 degrees C/SiC “buffer layer,” graphene on this (6 root 3 x 6 root 3) R30 degrees C/SiC buffer layer, and graphene decoupled from the SiC substrate by hydrogen intercalation. Using angle-resolved photoemission, core level photoemission, and x-ray absorption, we show that the electronic, chemical, and physical structure of all three variants is remarkably similar, exhibiting a large band gap and a vanishing density of states at the Fermi level. These results are explained in terms of first-principles calculations. This material thus appears very suitable for applications, even more so since it is prepared on a processing-friendly substrate. We also investigate two separate UV photon-induced modifications of the electronic structure that transform the insulating samples (6.2-eV band gap) into semiconducting (similar to 2.5-eV band gap) and metallic regions, respectively.

Abstract: The structural properties of partially hydrogenated and fluorinated graphene with different percentages of H/F atoms are investigated using molecular dynamics simulations based on reactive force field (ReaxFF) potentials. We found that the roughness of graphene varies with the percentage (p) of H or F and in both cases is maximal around p = 50%. Similar results were obtained for partially oxidized graphene. The two-dimensional area size of partially fluorinated and hydrogenated graphene exhibits a local minimum around p = 35% coverage. The lattice thermal contraction in partially functionalized graphene is found to be one order of magnitude larger than that of fully covered graphene. We also show that the armchair structure for graphene oxide (similar to the structure of fully hydrogenated and fluorinated graphene) is unstable. Our results show that the structure of partially functionalized graphene changes nontrivially with the C : H and C : F ratio as well as with temperature.

Abstract: We study transport properties of a MoS2 monolayer in the presence of a perpendicular magnetic field B. We derive and discuss its band structure and take into account spin and valley Zeeman effects. Compared to a conventional two-dimensional electron gas, these effects lead to new quantum Hall plateaus and new peaks in the longitudinal resistivity as functions of the magnetic field. The field B leads to a significant enhancement of the spin splitting in the conduction band, to a beating of the Shubnikov-de Haas (SdH) oscillations in the low-field regime, and to their splitting in the high-field regime. The Zeeman fields suppress significantly the beating of the SdH oscillations in the low-field regime and strongly enhance their splitting at high fields. The spin and valley polarizations show a similar beating pattern at low fields and are clearly separated at high fields in which they attain a value higher than 90%.