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 electronic structure and transport properties of monolayer MoS2 are studied using a tight-binding approach coupled with the nonequilibrium Green's function method. A zigzag nanoribbon of MoS2 is conducting due to the intersection of the edge states with the Fermi level that is located within the bulk gap. We show that applying a transverse electric field results in the disappearance of this intersection and turns the material into a semiconductor. By increasing the electric field the band gap undergoes a two stage linear increase after which it decreases and ultimately closes. It is shown that in the presence of a uniform exchange field, this electric field tuning of the gap can be exploited to open low energy domains where only one of the spin states contributes to the electronic conductance. This introduces possibilities in designing spin filters for spintronic applications.

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.

Abstract: The Aharonov-Bohm (AB) effect in square phosphorene quantum rings, with armchair and zigzag edges, is investigated using the tight-binding method. The energy spectra and wave functions of such rings, obtained as a function of the magnetic flux Phi threading the ring, are strongly influenced by the ringwidthW, an in-plane electric field E-p, and a side-gating potential V-g. Compared to a square dot, the ring shows an enhanced confinement due to its inner edges and an interedge coupling along the zigzag direction, both of which strongly affect the energy spectrum and the wave functions. The energy spectrum that is gapped consists of a regular part, of conduction (valence) band states, that shows the usual AB oscillations in the higher-(lower-) energy region, and of edge states, in the gap, that exhibit no AB oscillations. As the width W decreases, the AB oscillations become more distinct and regular and their period is close to Phi(0)/2, where the flux quantum Phi(0) = h/e is the period of an ideal circular ring (W -> 0). Both the electric field E-p and the side-gating potential V-g reduce the amplitude of the AB oscillations. The amplitude can be effectively tuned by E-p or V-g and exhibits an anisotropic behavior for different field directions or side-gating configurations.

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: 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: Electrically gating of trilayer graphene (TLG) opens a band gap offering the possibility to electrically engineer TLG quantum dots. We study the energy levels of such quantum dots and investigate their dependence on a perpendicular magnetic field B and different types of stacking of the graphene layers. The dots are modeled as circular and confined by a truncated parabolic potential which can be realized by nanostructured gates or position-dependent doping. The energy spectra exhibit the intervalley symmetry E-K(e) (m) = -E (h)(K') (m) for the electron (e) and hole (h) states, where m is the angular momentum quantum number and K and K' label the two valleys. The electron and hole spectra for B = 0 are twofold degenerate due to the intervalley symmetry E-K (m) = E-K' [-(m + 1)]. For both ABC [alpha = 1.5 (1.2) for large (small) R] and ABA (alpha = 1) stackings, the lowest-energy levels show approximately a R-alpha dependence on the dot radius R in contrast with the 1/R-3 one for ABC-stacked dots with infinite-mass boundary. As functions of the field B, the oscillator strengths for dipole-allowed transitions differ drastically for the two types of stackings.

Abstract: Two-dimensional systems of inverse patchy colloids modeled as disks with a central charge and having their surface decorated with oppositely pointlike charged patches are investigated using molecular dynamics simulations. The self-assembly of the patchy colloids leads to diverse ground state configurations ranging from crystalline arrangements of monomers to linear clusters, ramified linear clusters and to percolated configurations. Two structural phase diagrams are constructed: (1) as a function of the net charge and area fraction, and (2) as a function of the net charge and the range of the pair interaction potential. An interesting reentrant percolation transition is obtained as a function of the net charge of the colloids. We identify distinct mechanisms that lead to the percolation transition.

Abstract: We demonstrate the controlled spatial nucleation of superconductivity in a thin film deposited on periodic arrays of ferromagnetic dots with gradually increasing diameter. The perpendicular magnetization of the dots induces vortex-antivortex molecules in the sample, with the number of (anti)vortices increasing with magnet size. The resulting gradient of antivortex density between the dots predetermines local nucleation of superconductivity in the sample as a function of the applied external field and temperature. In addition, the compensation between the applied magnetic field and the antivortices results in an unprecedented enhancement of the critical temperature.

Abstract: We find localized electron and hole states in a ring-shaped potential kink in biased bilayer graphene. Within the continuum description, we show that for sharp potential steps the Dirac equation describing carrier states close to the K (or K′) point of the first Brillouin zone can be solved analytically for a circular kink/antikink dot. The solutions exhibit interfacial states which exhibit AharonovBohm oscillations as functions of the height of the potential step and/or the radius of the ring.