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: Water confined between two graphene layers with a separation of a few A forms a layered two-dimensional ice structure. Using large scale molecular dynamics simulations with the adoptable ReaxFF interatomic potential we found that flat monolayer ice with a rhombic-square structure nucleates between the graphene layers which is nonpolar and nonferroelectric. We provide different energetic considerations and H-bonding results that explain the interlayer and intralayer properties of two-dimensional ice. The controversial AA stacking found experimentally [Algara-Siller et al., Nature (London) 519, 443 (2015)] is consistent with our minimum-energy crystal structure of bilayer ice. Furthermore, we predict that an odd number of layers of ice has the same lattice structure as monolayer ice, while an even number of ice layers exhibits the square ice AA stacking of bilayer ice.

Abstract: In order to study quantum effects in a two-dimensional crystal lattice of a free-standing monolayer graphene, we have performed both path-integral Monte Carlo (PIMC) and classical Monte Carlo (MC) simulations for temperatures up to 2000 K. The REBO potential is used for the interatomic interaction. The total energy, interatomic distance, root-mean-square displacement of the atom vibrations, and the free energy of the graphene layer are calculated. The obtained lattice vibrational energy per atom from the classical MC simulation is very close to the energy of a three-dimensional harmonic oscillator 3k(B)T. The PIMC simulation shows that quantum effects due to zero-point vibrations are significant for temperatures T < 1000 K. The quantum contribution to the lattice vibrational energy becomes larger than that of the classical lattice for T < 400 K. The lattice expansion due to the zero-point motion causes an increase of 0.53% in the lattice parameter. A minimum in the lattice parameter appears at T similar or equal to 500 K. Quantum effects on the atomic vibration amplitude of the graphene lattice and its free energy are investigated.

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: 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: Fermi's golden rule underpins the investigation of mobile carriers propagating through various solids, being a standard tool to calculate their scattering rates. As such, it provides a perturbative estimate under the implicit assumption that the effect of the interaction Hamiltonian which causes the scattering events is sufficiently small. To check the validity of this assumption, we present a general framework to derive simple validity criteria in order to assess whether the scattering rates can be trusted for the system under consideration, given its statistical properties such as average size, electron density, impurity density et cetera. We derive concrete validity criteria for metallic nanowires with conduction electrons populating a single parabolic band subjected to different elastic scattering mechanisms: impurities, grain boundaries and surface roughness.

Abstract: The effects of strain, induced by a Gaussian bump, on the magnetic field dependent transport properties of a graphene Hall bar are investigated. The numerical simulations are performed using both classical and quantum mechanical transport theory and we found that both approaches exhibit similar characteristic features. The effects of the Gaussian bump are manifested by a decrease of the bend resistance, RB, around zero-magnetic field and the occurrence of side-peaks in RB. These features are explained as a consequence of bump-assisted scattering of electrons towards different terminals of the Hall bar. Using these features we are able to give an estimate of the size of the bump. Additional oscillations in RB are found in the quantum description that are due to the population/depopulation of Landau levels. The bump has a minor influence on the Hall resistance even for very high values of the pseudo-magnetic field. When the bump is placed outside the center of the Hall bar valley polarized electrons can be collected in the leads.

Abstract: The interaction of monolayer graphene with specific substrates may break its sublattice symmetry and results in unidirectional chiral states with opposite group velocities in the different Dirac cones (Zarenia et al 2012 Phys. Rev. B 86 085451). Taking advantage of this feature, we propose a valley filter based on a transversal mass kink for low energy electrons in graphene, which is obtained by assuming a defect region in the substrate that provides a change in the sign of the substrate-induced mass and thus creates a non-biased channel, perpendicular to the kink, for electron motion. By solving the time-dependent Schrodinger equation for the tight-binding Hamiltonian, we investigate the time evolution of a Gaussian wave packet propagating through such a system and obtain the transport properties of this graphene-based substrate-induced quantum point contact. Our results demonstrate that efficient valley filtering can be obtained, provided: (i) the electron energy is sufficiently low, i.e. with electrons belonging mostly to the lowest sub-band of the channel, and (ii) the channel length (width) is sufficiently long (narrow). Moreover, even though the transmission probabilities for each valley are significantly affected by impurities and defects in the channel region, the valley polarization in this system is shown to be robust against their presence.

Abstract: Bilayer graphene can exhibit deformations such that the two graphene sheets are locally detached from each other resulting in a structure consisting of domains with different van der Waals inter-layer coupling. Here we investigate how the presence of these domains affects the transport properties of bilayer graphene. We derive analytical expressions for the transmission probability, and the corresponding conductance, across walls separating different inter-layer coupling domains. We find that the transmission can exhibit a valley-dependent layer asymmetry and that the domain walls have a considerable effect on the chiral tunnelling properties of the charge carriers. We show that transport measurements allow one to obtain the strength with which the two layers are coupled. We perform numerical calculations for systems with two domain walls and find that the availability of multiple transport channels in bilayer graphene significantly modifies the conductance dependence on inter-layer potential asymmetry.

Abstract: A theoretical and experimental study of the work function of few-layer graphene is reported. The influence of the number of layers on the work function is investigated in the presence of a substrate, a molecular dipole layer, and combinations of the two. The work function of few-layer graphene is almost independent of the number of layers with only a difference between monolayer and multilayer graphene of about 60 meV. In the presence of a charge-donating substrate the charge distribution is found to decay exponentially away from the substrate and this is directly reflected in the work function of few-layer graphene. A dipole layer changes the work function only when placed in between the substrate and few-layer graphene through a change of the charge transfer between the two.

Abstract: We study topological phase transitions and topological quantum field effect transistor in monolayer molybdenum disulfide (MoS2) using a two-band Hamiltonian model. Without considering the quadratic (q(2)) diagonal term in the Hamiltonian, we show that the phase diagram includes quantum anomalous Hall effect, quantum spin Hall effect, and spin quantum anomalous Hall effect regions such that the topological Kirchhoff law is satisfied in the plane. By considering the q(2) diagonal term and including one valley, it is shown that MoS2 has a non-trivial topology, and the valley Chern number is non-zero for each spin. We show that the wave function is (is not) localized at the edges when the q(2) diagonal term is added (deleted) to (from) the spin-valley Dirac mass equation. We calculate the quantum conductance of zigzag MoS2 nanoribbons by using the nonequilibrium Green function method and show how this device works as a field effect topological quantum transistor.

Abstract: Bilayer graphene samples may exhibit regions where the two layers are locally delaminated forming a so-called quanttun blister in the graphene sheet. Electron and hole states can be confined in this graphene quantum blisters (GQB) by applying a global electrostatic bias. We scrutinize the electronic properties of these confined states under the variation of interlayer bias, coupling, and blister's size. The spectra display strong anti-crossings due to the coupling of the confined states on upper and lower layers inside the blister. These spectra are layer localized where the respective confined states reside on either layer or equally distributed. For finite angular momentum, this layer localization can be at the edge of the blister and corresponds to degenerate modes of opposite momenta. Furthermore, the energy levels in GQB exhibit electron-hole symmetry that is sensitive to the electrostatic bias. Finally, we demonstrate that confinement in GQB persists even in the presence of a variation in the interlayer coupling.

Abstract: In van der Waals heterostructures consisting of graphene and a substrate, lattice mismatch often leads to a moire pattern with a huge supercell, preventing its treatment within first- principles calculations. Previous theoretical works considered mostly simple stacking models such as AB, AA with straining the lattice of graphene to match that of the substrate. Here, we propose a moire superlattice build from graphene and porous graphene or graphyne like monolayers, having a lower interlayer binding energy, needing little strain in order to match the lattices. In contrast to the results from the simple stacking models, the present ab initio calculations for the moire superlattices show different properties in lattice structure, energy, and band structures. For example, the Dirac cone at the K point is preserved and a linear energy dispersion near the Fermi level is obtained.

Abstract: The bandgap behavior of 2D-GaAs and graphene have been investigated with van der Waals heterostructured into a yet unexplored graphene/GaAs bilayer, under both uniaxial stress along c axis and different planar strain distributions. The 2D-GaAs bandgap nature changes from Gamma-K indirect in isolated monolayer to Gamma-Gamma direct in graphene/GaAs bilayer. In the latter, graphene exhibits a bandgap of 5 meV. The uniaxial stress strongly affects the graphene electronic bandgap, while symmetric in-plane strain does not open the bandgap in graphene. Nevertheless, it induces remarkable changes on the GaAs bandgap-width around the Fermi level. However, when applying asymmetric in-plane strain to graphene/GaAs, the graphene sublattice symmetry is broken, and the graphene bandgap is open at the Fermi level to a maximum width of 814 meV. This value is much higher than that reported for just graphene under asymmetric strain. The Gamma-Gamma direct bandgap of GaAs remains unchanged in graphene/ GaAs under different types of applied strain. The analyses of phonon dispersion and the elastic constants yield the dynamical and mechanical stability of the graphene/GaAs system, respectively. The calculated mechanical properties for bilayer heterostructure are better than those of their constituent monolayers. This finding, together with the tunable graphene bandgap not only by the strength but also by the direction of the strain, enhance the potential for strain engineering of ultrathin group-III-V electronic devices hybridized by graphene.

Abstract: First-principles calculations have been carried out to investigate the stability, structural and electronic properties of two-dimensional (2D) hydrogenated GaAs with three possible geometries: chair, zigzag-line and boat configurations. The effect of van der Waals interactions on 2D H-GaAs systems has also been studied. These configurations were found to be energetic and dynamic stable, as well as having a semiconducting character. Although 2D GaAs adsorbed with H tends to form a zigzag-line configuration, the energy differences between chair, zigzag-line and boat are very small which implies the metastability of the system. Chair and boat configurations display a – direct bandgap nature, while pristine 2D-GaAs and zigzag-line are indirect semiconductors. The bandgap sizes of all configurations are also hydrogen dependent, and wider than that of pristine 2D-GaAs with both PBE and HSE functionals. Even though DFT-vdW interactions increase the adsorption energies and reduce the equilibrium distances of H-GaAs systems, it presents, qualitatively, the same physical results on the stability and electronic properties of our studied systems with PBE functional. According to our results, 2D buckled gallium arsenide is a good candidate to be synthesized by hydrogen surface passivation as its group III-V partners 2D buckled gallium nitride and boron nitride. The hydrogenation of 2D-GaAs tunes the bandgap of pristine 2D-GaAs, which makes it a potential candidate for optoelectronic applications in the blue and violet ranges of the visible electromagnetic spectrum.