Abstract: Fluorinated graphene is one of the most important derivatives of graphene and has been found to have great potential in optoelectronic and photonic nanodevices. However, the stability of F atoms on fluorinated graphene under different conditions, which is essential to maintain the desired properties of fluorinated graphene, is still unclear. In this work, we investigate the diffusion of F atoms on pristine graphene, graphene with defects, and at graphene/fluorographene interfaces by using density functional theory calculations. We find that an isolated F atom diffuses easily on graphene, but those F atoms can be localized by inducing vacancies or absorbates in graphene and by creating graphene/fluorographene interfaces, which would strengthen the binding energy of F atoms on graphene and increase the diffusion energy barrier of F atoms remarkably.

Abstract: Stable monolayer materials based on existing, well known and stable two-dimensional crystal fluorographene and molybdenum disulfide are predicted to exhibit a huge magnetocrystalline anisotropy when functionalized with adsorbed transition metal atoms at vacant sides. Ab initio calculations within the density-functional theory formalism were performed to investigate the adsorption of the transitional metals in a single S (or F) vacancy of monolayer molybdenum disulfide (or fluorographene). We found strong bonding of the transitional metal atoms to the vacant sites with binding energies ranging from 2.5 to 5.2 eV. Our calculations revealed that these systems with adsorbed metal atoms exhibit a magnetic anisotropy, specifically the structures including Os and Ir show a giant magnetocrystalline anisotropy energy of 31-101 meV. Our results demonstrate the possibility of obtaining stable monolayer materials with huge magnetocrystalline anisotropy based on preexisting, well known and stable two-dimensional crystals: fluorographene and molybdenum disulfide. We believe that the results obtained here are useful not only for deeper understanding of the origin of magnetocrystalline anisotropy but also for the design of monolayer optoelectronic devices with novel functionalities.

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 structural and mechanical properties of N-doped graphene (NG) are investigated using reactive force field (ReaxFF) potentials in large-scale molecular dynamics simulations. We found that ripples, which are induced by the dopants, change the roughness of NG, which depends on the number of dopants and their local arrangement. For any doping ratio N/C, the NG becomes ferroelectric with a net dipole moment. The formation energy increases nonlinearly with N/C ratio, while the Young's modulus, tensile strength, and intrinsic strain decrease with the number of dopants. Our results for the structural deformation and the thermoelectricity of the NG sheet are in good agreement with recent experiments and ab initio calculations.

Abstract: Motivated by a recent experiment that reported the synthesis of a new 2D material nitrogenated holey graphene (C2N) [Mahmood et al., Nat. Commun., 2015, 6, 6486], the electronic, magnetic, and mechanical properties of nitrogenated (C2N), phosphorated (C2P) and arsenicated (C2As) monolayer holey graphene structures are investigated using first-principles calculations. Our total energy calculations indicate that, similar to the C2N monolayer, the formation of the other two holey structures are also energetically feasible. Calculated cohesive energies for each monolayer show a decreasing trend going from the C2N to C2As structure. Remarkably, all the holey monolayers considered are direct band gap semiconductors. Regarding the mechanical properties (in-plane stiffness and Poisson ratio), we find that C2N has the highest in-plane stiffness and the largest Poisson ratio among the three monolayers. In addition, our calculations reveal that for the C2N, C2P and C2As monolayers, creation of N and P defects changes the semiconducting behavior to a metallic ground state while the inclusion of double H impurities in all holey structures results in magnetic ground states. As an alternative to the experimentally synthesized C2N, C2P and C2As are mechanically stable and flexible semiconductors which are important for potential applications in optoelectronics.

Abstract: Finding novel atomically thin heterostructures and understanding their characteristic properties are critical for developing better nanoscale optoelectronic devices. In this study, we investigate the electronic and optical properties of a GaS-Ca(OH)(2) heterostructure using first-principle calculations. The band gap of the GaS-Ca(OH)(2) heterostructure is significantly reduced when compared to those of the isolated constituent layers. Our calculations showthat the GaS-Ca(OH)(2) heterostructure is a type-II heterojunction which can be used to separate photoinduced charge carriers where electrons are localized in GaS and holes in the Ca(OH)(2) layer. This leads to spatially indirect excitons which are important for solar energy and optoelectronic applications due to their long lifetime. By solving the Bethe-Salpeter equation on top of a single shot GW calculation (G(0)W(0)), the dielectric function and optical oscillator strength of the constituent monolayers and the heterostructure are obtained. The oscillator strength of the optical transition for the GaS monolayer is an order of magnitude larger than the Ca(OH)(2) monolayer. We also found that the calculated optical spectra of different stacking types of the heterostructure show dissimilarities, although their electronic structures are rather similar. This prediction can be used to determine the stacking type of ultrathin heterostructures.

Abstract: Anderson's theorem states that non-magnetic impurities do not change the bulk properties of conventional superconductors. However, as the dimensionality is reduced, the effect of impurities becomes more significant. Here we investigate superconducting nanowires with diameter comparable to the Fermi wavelength $\lambda_F$ (which is less than the superconducting coherence length) by using a microscopic description based on the Bogoliubov-de Gennes method. We find that: 1) impurities strongly affect the superconducting properties, 2) the effect is impurity position dependent, and 3) it exhibits opposite behavior for resonant and off-resonant wire widths. We show that this is due to the interplay between the shape resonances of the order parameter and the subband energy spectrum induced by the lateral quantum confinement. These effects can be used to manipulate the Josephson current, filter electrons by subband and investigate the symmetries of the superconducting subband gaps.

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%.

Abstract: The effect of the number of stacking layers and the type of stacking on the electronic and optical properties of bilayer and trilayer black phosphorus are investigated by using first-principles calculations within the framework of density functional theory. We find that inclusion of many-body effects (i.e., electron-electron and electron-hole interactions) modifies strongly both the electronic and optical properties of black phosphorus. While trilayer black phosphorus with a particular stacking type is found to be a metal by using semilocal functionals, it is predicted to have an electronic band gap of 0.82 eV when many-body effects are taken into account within the G(0)W(0) scheme. Though different stacking types result in similar energetics, the size of the band gap and the optical response of bilayer and trilayer phosphorene are very sensitive to the number of layers and the stacking type. Regardless of the number of layers and the type of stacking, bilayer and trilayer black phosphorus are direct band gap semiconductors whose band gaps vary within a range of 0.3 eV. Stacking arrangements that are different from the ground state structure in both bilayer and trilayer black phosphorus exhibit significant modified valence bands along the zigzag direction and result in larger hole effective masses. The optical gap of bilayer (trilayer) black phosphorus varies by 0.4 (0.6) eV when changing the stacking type. The calculated binding energy of the bound exciton hardly changes with the type of stacking and is found to be 0.44 (0.30) eV for bilayer (trilayer) phosphorous.

Abstract: We study ballistic transport of Dirac fermions in silicene through arrays of barriers, of width d, in the presence of an exchange field M and a tunable potential of height U or depth-U. The spin-and valley-resolved conductances as functions of U or M, exhibit resonances away from the Dirac point (DP) and close to it a pronounced dip that becomes a gap when a critical electric field E-z is applied. This gap widens by increasing the number of barriers and can be used to realize electric field-controlled switching of the current. The spin p(s) and valley p(v) polarizations of the current near the DP increase with Ez or M and can reach 100% for certain of their values. These field ranges widen significantly by increasing the number of barriers. Also, ps and pv oscillate nearly periodically with the separation between barriers or wells and can be inverted by reversing M.

Abstract: Experiments are supplemented with ab initio density functional theory (DFT) calculations in order to investigate how the structural, electronic and optical properties of zinc oxide (ZnO) thin films are modified upon Cu doping. Changes in characteristic properties of doped thin films, that are deposited on a glass substrate by sol-gel dip coating technique, are monitored using X-ray diffraction (XRD) and UV measurements. Our ab initio calculations show that the electronic structure of ZnO can be well described by DFT+U/G(0)W(0) method and we find that Cu atom substitutional doping in ZnO is the most favourable case. Our XRD measurements reveal that the crystallite size of the films decrease with increasing Cu doping. Moreover, we determine the optical constants such as refractive index, extinction coefficient, optical dielectric function and optical energy band gap values of the films by means of UV-Vis transmittance spectra. The optical band gap of ZnO the thin film linearly decreases from 3.25 to 3.20 eV at 5% doping. In addition, our calculations reveal that the electronic defect states that stem from Cu atoms are not optically active and the optical band gap is determined by the ZnO band edges. Experimentally observed structural and optical results are in good agreement with our theoretical results.

Abstract: The valley splitting of the conduction-band states in a thin silicon-on-insulator layer is investigated using the 30-band k . p theory. The system composed of a few nm thick Si layer embedded within thick SiO2 layers is analyzed. The valley split states are found to cross periodically with increasing quantum well width, and therefore the energy splitting is an oscillatory function of the quantum well width, with period determined by the wave vector K-0 of the conduction band minimum. Because the valley split states are classified by parity, the optical transition between the ground hole state and one of those valley split conduction band states is forbidden. The oscillations in the valley splitting energy decrease with electric field and with smoothing of the composition profile between the well and the barrier by diffusion of oxygen from the SiO2 layers to the Si quantum well. Such a smoothing also leads to a decrease of the interband transition matrix elements. The obtained results are well parametrized by the effective two-valley model, but are found to disagree from previous 30-band calculations. This discrepancy could be traced back to the fact that the basis for the numerical solution of the eigenproblem must be restricted to the first Brillouin zone in order to obtain quantitatively correct results for the valley splitting.

Abstract: Anharmonic effects in an atomic monolayer thin crystal with honeycomb lattice structure are investigated by analytical and numerical lattice dynamical methods. Starting from a semiempirical model for anharmonic couplings of third and fourth orders, we study the in-plane and out-of-plane (flexural) mode components of the generalized wave vector dependent Gruneisen parameters, the thermal tension and the thermal expansion coefficients as a function of temperature and crystal size. From the resonances of the displacement-displacement correlation functions, we obtain the renormalization and decay rate of in-plane and flexural phonons as a function of temperature, wave vector, and crystal size in the classical and in the quantum regime. Quantitative results are presented for graphene. There, we find that the transition temperature T-alpha from negative to positive thermal expansion is lowered with smaller system size. Renormalization of the flexural mode has the opposite effect and leads to values of T-alpha approximate to 300 K for systems of macroscopic size. Extensive numerical analysis throughout the Brillouin zone explores various decay and scattering channels. The relative importance of normal and umklapp processes is investigated. The work is complementary to crystalline membrane theory and computational studies of anharmonic effects in two-dimensional crystals.

Abstract: A strong in-plane magnetic field drastically alters the low-energy spectrum of bilayer graphene by separating the parabolic energy dispersion into two linear Dirac cones. The effect of this dramatic change on the transport properties strongly depends on the orientation of the in-plane magnetic field with respect to the propagation direction of the charge carriers and the angle at which they impinge on the electrostatic potentials. For magnetic fields oriented parallel to the potential boundaries an additional propagating mode that results from the splitting into Dirac cones enhances the transmission probability for charge carriers tunneling through the potentials and increases the corresponding conductance. Our results show that the chiral suppression of transmission at normal incidence, reminiscent of bilayer graphene's 2 pi Berry phase, is turned into a chiral enhancement when the magnetic field increases, thus indicating a transition from a bilayer to a monolayer-like system at normal incidence. Further, we find that the typical transmission resonances stemming from confinement in a potential barrier are shifted to higher energy and are eventually transformed into antiresonances with increasing magnetic field. For magnetic fields oriented perpendicular to the potential boundaries we find a very pronounced transition from a bilayer system to two separated monolayer-like systems with Klein tunneling emerging at certain incident angles symmetric around 0, which also leaves a signature in the conductance. For both orientations of the magnetic field, the transmission probability is still correctly described by pseudospin conservation. Finally, to motivate the large in-plane magnetic field, we show that its energy spectrum can be mimicked by specific lattice deformations such as a relative shift of one of the layers. With this equivalence we introduce the notion of an in-plane pseudomagnetic field.