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: The structural transitions of the ground state of a system of repulsively interacting particles confined in a quasi-one-dimensional channel, and the effect of the interparticle interaction as well as the functional form of the confinement potential on those transitions are investigated. Although the nonsequential ordering of transitions (non-SOT), i.e., the 1 – 2 – 4 – 3 – 4 – 5 – 6 – ... sequence of chain configurations with increasing density, is widely robust as predicted in a number of theoretical studies, the sequential ordering of transitions (SOT), i.e., the 1 – 2 – 3 – 4 – 5 – 6 – ... chain, is found as the ground state for long-ranged interparticle interaction and hard-wall-like confinement potentials. We found an energy barrier between every two different phases around its transition point, which plays an important role in the preference of the system to follow either a SOT or a non-SOT. However, that preferential transition requires also the stability of the phases during the transition. Additionally, we analyze the effect of a small structural disorder on the transition between the two phases around its transition point. Our results show that a small deformation of the triangular structure changes dramatically the picture of the transition between two phases, removing in a considerable region the non-SOT in the system. This feature could explain the fact that the non-SOT is, up to now, not observed in experimental systems, and suggests a more advanced experimental setup to detect the non-SOT.

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: Within a tight binding approach we investigate the energy levels of hexagonal and triangular bilayer graphene (BLG) quantum dots (QDs) with zigzag and armchair edges. We study AA- and AB-(Bernal) stacked BLG QDs and obtain the energy levels in both the absence and the presence of a perpendicular electric field (i.e., biased BLG QDs). Our results show that the size dependence of the energy levels is different from that of monolayer graphene QDs. The energy spectrum of AB-stacked BLG QDs with zigzag edges exhibits edge states which spread out into the opened energy gap in the presence of a perpendicular electric field. We found that the behavior of these edges states is different for the hexagonal and triangular geometries. In the case of AA-stacked BLG QDs, the electron and hole energy levels cross each other in both cases of armchair and zigzag edges as the dot size or the applied bias increases.

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 structural properties of partially hydrogenated and fluorinated graphene with different percentages of H/F atoms are investigated using molecular dynamics simulations based on reactive force field (ReaxFF) potentials. We found that the roughness of graphene varies with the percentage (p) of H or F and in both cases is maximal around p = 50%. Similar results were obtained for partially oxidized graphene. The two-dimensional area size of partially fluorinated and hydrogenated graphene exhibits a local minimum around p = 35% coverage. The lattice thermal contraction in partially functionalized graphene is found to be one order of magnitude larger than that of fully covered graphene. We also show that the armchair structure for graphene oxide (similar to the structure of fully hydrogenated and fluorinated graphene) is unstable. Our results show that the structure of partially functionalized graphene changes nontrivially with the C : H and C : F ratio as well as with temperature.

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: 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: The electronic properties, carrier mobility, and strain response of TiS3 nanoribbons (TiS3 NRs) are investigated by first-principles calculations. We found that the electronic properties of TiS3 NRs strongly depend on the edge type (a or b). All a-TiS3 NRs are metallic with a magnetic ground state, while b-TiS3 NRs are direct band gap semiconductors. Interestingly, the size of the band gap and the band edge position are almost independent of the ribbon width. This feature promises a constant band gap in a b-TiS3 NR with rough edges, where the ribbon width differs in different regions. The maximum carrier mobility of b-TiS3 NRs is calculated by using the deformation potential theory combined with the effective mass approximation and is found to be of the order 10(3) cm(2) V-1 s(-1). The hole mobility of the b-TiS3 NRs is one order of magnitude lower, but it is enhanced compared to the monolayer case due to the reduction in hole effective mass. The band gap and the band edge position of b-TiS3 NRs are quite sensitive to applied strain. In addition we investigate the termination of ribbon edges by hydrogen atoms. Upon edge passivation, the metallic and magnetic features of a-TiS3 NRs remain unchanged, while the band gap of b-TiS3 NRs is increased significantly. The robust metallic and ferromagnetic nature of a-TiS3 NRs is an essential feature for spintronic device applications. The direct, width-independent, and strain-tunable band gap, as well as the high carrier mobility, of b-TiS3 NRs is of potential importance in many fields of nanoelectronics, such as field-effect devices, optoelectronic applications, and strain sensors.