Abstract: Using the tight-binding model, we investigate the influence of vacancy disorder on electrical transport in graphene Hall bars in the presence of quantizing magnetic fields. Disorder, induced by a random distribution of monovacancies, breaks the graphene sublattice symmetry and creates states localized on the vacancies. These states are observable in the bend resistance, as well as in the total DOS. Their energy is proportional to the square root of the magnetic field, while their localization length is proportional to the cyclotron radius. At the energies of these localized states, the electron current flows around the monovacancies and, as we show, it can follow unexpected paths depending on the particular arrangement of vacancies. We study how these localized states change with the vacancy concentration, and what are the effects of including the next-nearest-neighbor hopping term. Our results are also compared with the situation when double vacancies are present in the system. Double vacancies also induce localized states, but their energy and magnetic field dependencies are different. Their localization energy scales linearly with the magnetic field, and their localization length appears not to depend on the field strength.

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: The report on square ice sandwiched between two graphene layers by Algara-Siller et al. [Nature (London) 519, 443 (2015)] has generated a large interest in this system. By applying high lateral pressure on nanoconfined water, we found that monolayer ice is transformed to bilayer ice when the two graphene layers are separated by H = 6,7 angstrom. It was also found that three layers of a denser phase of ice with smaller lattice constant are formed if we start from bilayer ice and apply a lateral pressure of about 0.7 GPa with H = 8,9 angstrom. The lattice constant (2.5-2.6 angstrom) in both transitions is found to be smaller than those typical for the known phases of ice and water, i.e., 2.8 angstrom. We validate these results using ab initio calculations and find good agreement between ab initio O-O distance and those obtained from classical molecular dynamics simulations. The reversibility of the mentioned transitions is confirmed by decompressing the systems.

Abstract: We introduce a new approach for the correlation energy of one- and two-valley two-dimensional electron gas (2DEG) systems. Our approach is based on an interpolation between two limits, a random phase approximation at high densities and a classical approach at low densities which gives excellent agreement with available Quantum Monte Carlo (QMC) calculations. The two-valley 2DEG model is introduced to describe the electron correlations in monolayer transition metal dichalcogenides (TMDs). We study the zero-temperature transition from a Fermi liquid to a quantum Wigner crystal phase in monolayer TMDs. Consistent with QMC, we find that electrons crystallize at r(s) = 31 in one-valley 2DEG. For two valleys, we predict Wigner crystallization at r(s) = 30, implying that valley degeneracy has little effect on the critical r(s), in contrast to an earlier claim.

Abstract: Two-dimensional crystals with strong interactions between layers has attracted increasing attention in recent years in a variety of fields. In particular, the growth of a single layer of oxide materials (e.g., MgO, CaO, NiO, and MnO) over metallic substrates were found to display different physical properties than their bulk. In this study, we report on the physical properties of a single layer of metallic oxide materials and compare their properties with their bulk and other two-dimensional (2D) crystals. We found that the planar structure of metallic monoxides are unstable whereas the buckled structures are thermodynamically stable. Also, the 2D-MnO and NiO exhibit different magnetic (ferromagnetic) and optical properties than their bulk, whereas band-gap energy and linear stiffness are found to be decreasing from NiO to MgO. Our findings provide insight into oxide thin-film technology applications.

Abstract: The effects of external electric and magnetic fields on the energy spectrum of quantum rings made out of a bidimensional semiconductor material with anisotropic band structures are investigated within the effective-mass model. The interplay between the effective-mass anisotropy and the radial confinement leads to wave functions that are strongly localized at two diametrically opposite regions where the kinetic energy is lowest due to the highest effective mass. We show that this quantum phenomenon has clear consequences on the behavior of the energy states in the presence of applied in-plane electric fields and out-of-plane magnetic fields. In the former, the quantum confined Stark effect is observed with either linear or quadratic shifts, depending on the direction of the applied field. As for the latter, the usual Aharonov-Bohm oscillations are not observed for a circularly symmetric confining potential, however they can be reinstated if an elliptic ring with an appropriate aspect ratio is chosen.

Abstract: We study the influence of confinement on the atomic collapse due to a Coulomb impurity placed at the center of a graphene quantum dot of radius R. We apply the zigzag or infinite-mass boundary condition and consider both a point-size and a finite-size impurity. As a function of the impurity strength Za, the energy spectra are discrete. In the case of the zigzag boundary condition, the degenerate (with respect to the angular momentum m) zero-energy levels are pulled down in energy as Z alpha increases, and they remain below epsilon = – Z alpha. Our results show that the energy levels exhibit a 1/R dependence in the subcritical regime [Z alpha < |km + 1/2|, k = 1 (-1) for the K (K') valley]. In the supercritical regime (Z alpha > |km + 1/2|) we find a qualitatively very different behavior where the levels decrease as a function of R in a nonmonotonic manner. While the valley symmetry is preserved in the presence of the impurity, we find that the impurity breaks electron-hole symmetry. We further study the energy spectrum of zigzag quantum dots in gapped graphene. Our results show that as the gap increases, the lowest electron states are pushed into the gap by the impurity.

Abstract: The electronic and structural properties of excitons and trions in monolayer transition metal dichalcogenides are investigated using both a multiband and a single- band model. In the multiband model we construct the excitonic Hamiltonian in the product base of the single-particle states at the conduction and valence band edges. We decouple the corresponding energy eigenvalue equation and solve the resulting differential equation self-consistently, using the finite element method (FEM), to determine the energy eigenvalues and the wave functions. As a comparison, we also consider the simple single-band model which is often used in numerical studies. We solve the energy eigenvalue equation using the FEM as well as with the stochastic variational method (SVM) in which a variational wave function is expanded in a basis of a large number of correlated Gaussians. We find good agreement between the results of both methods, as well as with other theoretical works for excitons, and we also compare with available experimental data. For trions the agreement between both methods is not as good due to our neglect of angular correlations when using the FEM. Finally, when comparing the two models, we see that the presence of the valence bands in the mutiband model leads to differences with the single- band model when (interband) interactions are strong.

Abstract: By making use of circularly polarized light and electrostatic gating, monolayer molybdenum disulfide (ML – MoS2) can form a platform supporting multiple types of charge carriers. They can be discriminated by their spin, valley index, or whether they are electrons or holes. We investigate the collective properties of those charge carriers and are able to identify distinct plasmon modes. We analyze the corresponding dispersion relation, lifetime, and oscillator strength, and calculate the phase relation between the oscillations in the different components of the plasmon modes. All platforms in ML-MoS2 support a long-wavelength root q plasmon branch at zero kelvins. In addition to this, for an n-component system, n-1 distinct plasmon modes appear as acoustic modes with linear dispersion in the long-wavelength limit. These modes correspond to out-of-phase oscillations in the different fermion liquids and have, although being damped, a relatively long lifetime. Additionally, we also find distinct modes at large wave vectors that are more strongly damped by intraband processes.

Abstract: Starting from a Hamiltonian with anharmonic coupling between in-plane acoustic displacements and outof-plane (flexural) modes, we derived coupled equations of motion for in-plane displacements correlations and flexural mode density fluctuations. Linear response theory and time-dependent thermal Green's functions techniques are applied in order to obtain different response functions. As external perturbations we allow for stresses and thermal heat sources. The displacement correlations are described by a Dyson equation where the flexural density distribution enters as an additional perturbation. The flexural density distribution satisfies a kinetic equation where the in-plane lattice displacements act as a perturbation. In the hydrodynamic limit this system of coupled equations is at the basis of a unified description of elastic and thermal phenomena, such as isothermal versus adiabatic sound motion and thermal conductivity versus second sound. The general theory is formulated in view of application to graphene, two-dimensional h-BN, and 2H-transition metal dichalcogenides and oxides.

Abstract: The energy spectrum and eigenstates of single-layer black phosphorus nanoribbons in the presence of a perpendicular magnetic field and an in-plane transverse electric field are investigated by means of a tight-binding method, and the effect of different types of edges is examined analytically. A description based on a continuum model is proposed using an expansion of the tight-binding model in the long-wavelength limit. Thewave functions corresponding to the flatband part of the spectrum are obtained analytically and are shown to agree well with the numerical results from the tight-binding method for both narrow (10 nm) and wide (100 nm) nanoribbons. Analytical expressions for the critical magnetic field at which Landau levels are formed and the ranges of wave numbers in the dispersionless flatband segments in the energy spectra are derived. We examine the evolution of the Landau levels when an in-plane lateral electric field is applied, and we determine analytically how the edge states shift withmagnetic field. For wider nanoribbons, the conductance is shown to have a characteristic staircase shape in combined magnetic and electric fields. Some of the stairs in zigzag and skewed armchair nanoribbons originate from edge states that are found in the band gap.

Abstract: <script type='text/javascript'>document.write(unpmarked('We show that delaminations in bilayer graphene (BLG) with electrostatically induced interlayer symmetry can provide one with ballistic channels for electrons with energies inside the electrostatically induced BLG gap. These channels are formed by a combination of valley-polarized evanescent states propagating along the delamination edges (which persist in the presence of a strong magnetic field) and standing waves bouncing between them inside the delaminated region (in a strong magnetic field, these transform into Landau levels in the monolayers). For inverted stackings in BLGs on the left and right of the delamination (AB-2ML-BA or BA-2ML-AB, where 2ML indicates two decoupled monolayers of graphene), the lowest-energy ballistic channels are gapless, have linear dispersion, and appear to be weakly topologically protected. When BLG stackings on both sides of the delamination are the same (AB-2ML-AB or BA-2ML-BA), the lowest-energy ballistic channels are gapped, with a gap epsilon(g) scaling as epsilon(g) alpha W-1 with delamination width and epsilon(g) alpha delta(-1) with the on-layer energy difference in the delaminated part of the structure. Depending on the width, delaminations may also support several \u0022higher-energy\u0022 waveguide modes. Our results are based on both the analytical study of the wave matching of Dirac states and tight-binding model calculations, and we analyze in detail the dependence of the delamination spectrum on the electrostatic conditions in the structure, such as the vertical displacement field.'));

Abstract: Defects are almost inevitably present in realistic materials and defective materials are expected to exhibit very different properties than their nondefective (perfect) counterparts. Here, using a combination of the tight-binding approach and the scattering matrix formalism, we investigate the electronic transport properties of defective phosphorene nanoribbons (PNRs) containing atomic vacancies. We find that for both armchair PNRs (APNRs) and zigzag PNRs (ZPNRs), single vacancies can create quasilocalized states, which can affect their conductance. With increasing vacancy concentration, three different transport regimes are identified: ballistic, diffusive, and Anderson localized ones. In particular, ZPNRs that are known to be metallic due to the presence of edge states become semiconducting: edge conductance vanishes and transport gap appears due to Anderson localization. Moreover, we find that for a fixed vacancy concentration, both APNRs and ZPNRs of narrower width and/or longer length are more sensitive to vacancy disorder than their wider and/or shorter counterparts, and that for the same ribbon length and width, ZPNRs are more sensitive to vacancy disorder than APNRs.

Abstract: Using the tight-binding approach, we investigate the energy spectrum of square, triangular, and hexagonal MoS2 quantum dots (QDs) in the presence of a perpendicular magnetic field. Novel edge states emerge in MoS2 QDs, which are distributed over the whole edge which we call ring states. The ring states are robust in the presence of spin-orbit coupling (SOC). The corresponding energy levels of the ring states oscillate as a function of the perpendicular magnetic field which are related to Aharonov-Bohm oscillations. Oscillations in the magnetic field dependence of the energy levels and the peaks in the magneto-optical spectrum emerge (disappear) as the ring states are formed (collapsed). The period and the amplitude of the oscillation decrease with the size of the MoS2 QDs.