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: Trions and biexcitons in anisotropic two-dimensional materials are investigated within an effective mass theory. Explicit results are obtained for phosphorene and arsenene, materials that share features such as a direct quasiparticle gap and anisotropic conduction and valence bands. Trions are predicted to have remarkably high binding energies and an elongated electron-hole structure with a preference for alignment along the armchair direction, where the effective masses are lower. We find that biexciton binding energies are also notably large, especially for monolayer phosphorene, where they are found to be twice as large as those for typical monolayer transition metal dichalcogenides.

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.

Abstract: The effect of an applied electric field and the effect of charging are investigated on themagnetic anisotropy (MA) of various stable two-dimensional (2D) crystals such as graphene, FeCl2, graphone, fluorographene, and MoTe2 using first-principles calculations. We found that themagnetocrystalline anisotropy energy of Co-on-graphene and Os-doped-MoTe2 systems change linearly with electric field, opening the possibility of electric field tuningMAof these compounds. In addition, charging can rotate the easy-axis direction ofCo-on-graphene andOs-doped-MoTe2 systems from the out-of-plane (in-plane) to in-plane (out-of-plane) direction. The tunable MA of the studied materials is crucial for nanoscale electronic technologies such as data storage and spintronics devices. Our results show that controlling the MA of the mentioned 2D crystal structures can be realized in various ways, and this can lead to the emergence of a wide range of potential applications where the tuning and switching of magnetic functionalities are important.

Abstract: We investigate the superconducting states that arise as a consequence of mesoscopic confinement and a multicomponent order parameter in the Ginzburg-Landau model for p-wave superconductivity. Conventional vortices, but also half-quantum vortices and skyrmions, are found as the applied magnetic field and the anisotropy parameters of the Fermi surface are varied. The solutions are well differentiated by a topological charge that for skyrmions is given by the Hopf invariant and for vortices by the circulation of the superconducting velocity. We revealed several unique states combining vortices and skyrmions, their possible reconfiguration with varied magnetic field, as well as temporal and field-induced transitions between vortical and skyrmionic states.

Abstract: Using the continuum model, we investigate the confined states and the corresponding wave functions of ABC-stacked trilayer graphene (TLG) quantum dots (QDs). First, a general infinite-mass boundary condition is derived and applied to calculate the electron and hole energy levels of a circular QD in both the absence and presence of a perpendicular magnetic field. Our analytical results for the energy spectra agree with those obtained by using the tight-binding model, where a TLG QD is surrounded by a staggered potential. Our findings show that (i) the energy spectrum exhibits intervalley symmetry E-K(e)(m) = -E-K'(h)(m) for the electron (e) and hole (h) states, where m is the angular momentum quantum number, (ii) the zero-energy Landau level (LL) is formed by the magnetic states with m <= 0 for both Dirac valleys, that is different from monolayer and bilayer graphene QD with infinite-mass potential in which only one of the cones contributes, and (iii) groups of three quantum Hall edge states in the tight-binding magnetic spectrum approach the zero LL, which results from the layer symmetry in TLG QDs.

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: Magnesium hydroxide [Mg(OH)(2)] has a layered brucitelike structure in its bulk form and was recently isolated as a new member of two-dimensional monolayer materials. We investigated the electronic and optical properties of monolayer crystals of Mg(OH)(2) and WS2 and their possible heterobilayer structure by means of first-principles calculations. It was found that both monolayers of Mg(OH)(2) and WS2 are direct-gap semiconductors and these two monolayers form a typical van der Waals heterostructure with a weak interlayer interaction and a type-II band alignment with a staggered gap that spatially separates electrons and holes. We also showed that an out-of-plane electric field induces a transition from a staggered to a straddling-type heterojunction. Moreover, by solving the Bethe-Salpeter equation on top of single-shot G(0)W(0) calculations, we show that the low-energy spectrum of the heterobilayer is dominated by the intralyer excitons of the WS2 monolayer. Because of the staggered interfacial gap and the field-tunable energy-band structure, the Mg(OH)(2)-WS2 heterobilayer can become an important candidate for various optoelectronic device applications in nanoscale.

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 magnetic-field dependence of the energy spectrum, wave function, binding energy, and oscillator strength of exciton states confined in a circular graphene quantum dot (CGQD) is obtained within the configuration interaction method. We predict that (i) excitonic effects are very significant in the CGQD as a consequence of a combination of geometric confinement, magnetic confinement, and reduced screening; (ii) two types of excitons (intravalley and intervalley) are present in the CGQD because of the valley degree of freedom in graphene; (iii) the intravalley and intervalley exciton states display different magnetic-field dependencies due to the different electron-hole symmetries of the single-particle energy spectra; (iv) with increasing magnetic field, the exciton ground state in the CGQD undergoes an intravalley to intervalley transition accompanied by a change of angular momentum; (v) the exciton binding energy does not increase monotonically with the magnetic field due to the competition between geometric and magnetic confinements; and (vi) the optical transitions of the intervalley and intravalley excitons can be tuned by the magnetic field, and valley-dependent excitonic transitions can be realized in a CGQD.

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: The elastic constant C-11 and piezoelectric stress constant e(1),(11) of two-dimensional (2D) dielectric materials comprising h-BN, 2H-MoS2, and other transition-metal dichalcogenides and dioxides are calculated using lattice dynamical theory. The results are compared with corresponding quantities obtained with ab initio calculations. We identify the difference between clamped-ion and relaxed-ion contributions with the dependence on inner strains which are due to the relative displacements of the ions in the unit cell. Lattice dynamics allows us to express the inner-strain contributions in terms of microscopic quantities such as effective ionic charges and optoacoustical couplings, which allows us to clarify differences in the piezoelectric behavior between h-BN and MoS2. Trends in the different microscopic quantities as functions of atomic composition are discussed.

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: An external electric field changes the physical properties of polar liquids due to the reorientation of their permanent dipoles. Using molecular dynamics simulations, we predict that an in-plane electric field applied parallel to the channel polarizes water molecules which are confined between two graphene layers, resulting in distinct ferroelectricity and electrical hysteresis. We found that electric fields alter the in-plane order of the hydrogen bonds: Reversing the electric field does not restore the system to the nonpolar initial state, instead a residual dipole moment remains in the system. The square-rhombic structure of 2D ice is transformed into two rhombic-rhombic structures. Our study provides insights into the ferroelectric state of water when confined in nanochannels and shows how this can be tuned by an electric field.