Abstract: By using density functional theory-based first-principles calculations, we investigate the structural, electronic, optical, and transport properties of pristine single-layer BC6N. Under different external actions and functionalization. Increasing the thickness of the structure results in a decrease of the band gap. Applying a perpendicular electric field decreases the band gap and a semiconductor-to-topological insulator transition is revealed. Uniaxial and biaxial strains of +8% result in a semiconductor-to-metal transition. Nanoribbons of BC6N having zigzag edge with even (odd) values of widths, become metal (semiconductor), while the armchair edge nanoribbons exhibit robust semiconducting behavior. In addition, we systematically investigate the effect of surface adatom and molecule, substitutional impurity and defect engineering on the electronic properties of single-layer BC6N and found transitions from metal to half-metal, to ferromagnetic metal, to dilute magnetic semiconductor, and even to spin-glass semiconductor. Furthermore we found that, topological defects including vacancies and Stone–Wales type, induce magnetism in single-layer BC6N.

Abstract: Two-dimensional (2D) materials are promising candidates for atomic-scale piezotronics and piezophototronics. Quantum edge states show fascinating fundamental physics such as nontrivial topological behavior and hold promising practical applications for low-power electronic devices. Here, using the tight-binding approach and quantum transport simulations, we investigate the piezotronic effect on the spin polarization of edge states in a zigzag-terminated monolayer MoS2 nanoribbon. We find that the strain-induced piezoelectric potential induces a phase transition of edge states from metal to semiconductor. However, in the presence of exchange field, edge states become semi-metallic with significant spin splitting and polarization that can be tuned by external strain. We show that quantum transport conductance exhibits a 100% spin polarization over a wide range of strain magnitudes. This effect is used in a propose prototype of piezotronic spin transistor. Our results provide a fundamental understanding of the piezotronic effect on edge states in zigzag monolayer MoS2 nanoribbons and are relevant for designing high-performance piezotronic spin devices.

Abstract: Water permeation between stacked layers of hBN sheets forming 2D nanochannels is investigated using large-scale ab initio-quality molecular dynamics simulations. A high-dimensional neural network potential trained on density-functional theory calculations is employed. We simulate water in van der Waals nanocapillaries and study the impact of nanometric confinement on the structure and dynamics of water using both equilibrium and nonequilibrium methods. At an interlayer distance of 10.2 A confinement induces a first-order phase transition resulting in a well-defined AA-stacked bilayer of hexagonal ice. In contrast, for h < 9 A, the 2D water monolayer consists of a mixture of different locally ordered patterns of squares, pentagons, and hexagons. We found a significant change in the transport properties of confined water, particularly for monolayer water where the water-solid friction coefficient decreases to half and the diffusion coefficient increases by a factor of 4 as compared to bulk water. Accordingly, the slip-velocity is found to increase under confinement and we found that the overall permeation is dominated by monolayer water adjacent to the hBN membranes at extreme confinements. We conclude that monolayer water in addition to bilayer ice has a major contribution to water transport through 2D nanochannels.

Abstract: Two-dimensional graphene-based van der Waals heterostructures have received considerable interest because of their intriguing characteristics compared with the constituent single-layer two-dimensional materials. Here, we investigate the interfacial characteristics, Schottky contact, and optical performance of graphene/Ga2SSe van der Waals (vdW) heterostructure using first-principles calculations. The effects of stacking patterns, electric gating, and interlayer coupling on the interfacial properties of graphene/Ga2SSe heterostructures are also examined. Our results demonstrate that the Dirac cone of graphene is well preserved at the F point in all stacking patterns due to the weak vdW interactions, which keep the heterostructures feasible such that they can be obtained in further experiments. Moreover, depending on the stacking patterns, a small band gap of about 13-17 meV opens in graphene and has a high carrier mobility, indicating that the graphene/Ga2SSe heterostructures are potential candidates for future high-speed nanoelectronic applications. In the ground state, the graphene/Ga2SSe heterostructures form an n-type Schottky contact. The transformation from an n-type to a p-type Schottky contact or to an Ohmic contact can be forced by electric gating or by varying the interlayer coupling. Our findings could provide physical guidance for designing controllable Schottky nanodevices with high electronic and optical performances.

Abstract: The energy spectrum and local current patterns in graphene quantum dots (QD) are investigated for different geometries in the presence of an external perpendicular magnetic field. Our results demonstrate that, for specific geometries and edge configurations, the QD exhibits vortex and anti-vortex patterns in the local current density, in close analogy to the vortex patterns observed in the probability density current of semiconductor QD, as well as in the order parameter of mesoscopic superconductors.

Abstract: Atomic collapse can be observed in graphene because of its large “effective” fine structure constant, which enables this phenomenon to occur for an impurity charge as low as Z(c) similar to 1-2. Here we investigate the effect of the sublattice symmetry on molecular collapse in two spatially separated charge tunable vacancies, which are located on the same (A-A type) or different (A-B type) sublattices. We find that the broken sublattice symmetry: (1) does not affect the location of the main bonding and antibonding molecular collapse peaks, (2) but shifts the position of the satellite peaks, because they are a consequence of the breaking of the local sublattice symmetry, and (3) there are vacancy characteristic collapse peaks that only occur for A-B type vacancies, which can be employed to distinguish them experimentally from the A-A type. As the charge, energy, and separation distance increase, the additional collapse features merge with the main molecular collapse peaks. We show that the spatial distribution around the vacancy site of the collapse states allows us to differentiate the molecular from the frustrated collapse.

Abstract: The dielectric properties of confined water is of fundamental interest and is still controversial. For water confined in channels with height smaller than h = 8 angstrom, we found a commensurability effect and an extraordinary decrease in the out-of-plane dielectric constant down to the limit of the dielectric constant of optical water. Spatial resolved polarization density data obtained from molecular dynamics simulations are found to be antisymmetric across the channel and are used as input in a mean-field model for the dielectric constant as a function of the height of the channel for h > 15 angstrom. Our results are in excellent agreement with a recent experiment [L. Fumagalli et al., Science 360, 1339 (2018)].

Abstract: Using evolutionary algorithm for crystal structure prediction, we present a new stable two-dimensional (2D) carbon allotrope composed of polymerized as-indacenes (PAI) in a zigzag pattern, namely PAI-graphene whose energy is lower than most of the reported 2D allotropes of graphene. Crucially, the crystal structure realizes a nonsymmorphic layer group that enforces a nontrivial global topology of the band structure with two Dirac cones lying perfectly at the Fermi level. The absence of electron/hole pockets makes PAI-graphene a pristine crystalline topological semimetal having anisotropic Fermi velocities with a high value of 7.0 x 10(5) m/s. We show that while the semimetallic property of the allotrope is robust against the application of strain, the positions of the Dirac cone and the Fermi velocities can be modified significantly with strain. Moreover, by combining strain along both the x- and y-directions, two band inversions take place at G leading to the annihilation of the Dirac nodes demonstrating the possibility of strain-controlled conversion of a topological semimetal into a semiconductor. Finally we formulate the bulk-boundary correspondence of the topological nodal phase in the form of a generalized Zak-phase argument finding a perfect agreement with the topological edge states computed for different edge-terminations. (C) 2020 The Author(s). Published by Elsevier Ltd.

Abstract: Defect-free graphene is impermeable to gases and liquids but highly permeable to thermal protons. Atomic-scale defects such as vacancies, grain boundaries, and Stone-Wales defects are predicted to enhance graphene's proton permeability and may even allow small ions through, whereas larger species such as gas molecules should remain blocked. These expectations have so far remained untested in experiment. Here, we show that atomically thin carbon films with a high density of atomic-scale defects continue blocking all molecular transport, but their proton permeability becomes similar to 1000 times higher than that of defect-free graphene. Lithium ions can also permeate through such disordered graphene. The enhanced proton and ion permeability is attributed to a high density of eight-carbon-atom rings. The latter pose approximately twice lower energy barriers for incoming protons compared to that of the six-atom rings of graphene and a relatively low barrier of similar to 0.6 eV for Li ions. Our findings suggest that disordered graphene could be of interest as membranes and protective barriers in various Li-ion and hydrogen technologies.

Abstract: In the Comment, the authors state that the separation of the relative and center of mass variables in our work is not correct. Here we point out that there is a typographical error, i.e., qi instead of -e, in two of our equations which, when corrected, makes the Comment redundant. Within the ansatzes mentioned in our paper all our results are correct, in contrast to the claims of the Comment.

Abstract: By using density functional theory (DFT) calculations, we predict a puckered, dynamically stable Janus single-layer black arsenic-phosphorus (b-AsP), which is composed of two different atomic sublayers, arsenic and phosphorus atoms. The calculated phonon spectrum reveals that Janus single-layer b-AsP is dynamically stable with either pure or coupled optical phonon branches arising from As and P atoms. The calculated Raman spectrum indicates that due to the relatively strong P-P bonds, As atoms have no contribution to the highfrequency optical vibrations. In addition, the orientation-dependent isovolume heat capacity reveals anisotropic contributions of LA and TA phonon branches to the low-temperature thermal properties. Unlike pristine single layers of b-As and b-P, Janus single-layer b-AsP exhibits additional out-of-plane asymmetry which leads to important consequences for its electronic, optical, and elastic properties. In contrast to single-layer b-As, Janus single-layer b-AsP is found to possess a direct band gap dominated by the P atoms. Moreover, real and imaginary parts of the dynamical dielectric function, including excitonic effects, reveal the highly anisotropic optical feature of the Janus single-layer. A tight-binding (TB) model is also presented for Janus single-layer b-AsP, and it is shown that, with up to seven nearest hoppings, the TB model reproduces well the DFT band structure in the low-energy region around the band gap. This TB model can be used in combination with the Green's function approach to study, e.g., quantum transport in finite systems based on Janus single-layer b-AsP. Furthermore, the linear-elastic properties of Janus single-layer b-AsP are investigated, and the orientation-dependent in-plane stiffness and Poisson ratio are calculated. It is found that the Janus single layer exhibits strong in-plane anisotropy in its Poisson ratio much larger than that of single-layer b-P. This Janus single layer is relevant for promising applications in optical dichroism and anisotropic nanoelasticity.

Abstract: We report characteristic vortex configurations in s + id superconductors with time-reversal symmetry breaking, exposed to magnetic field. A vortex in the s + id state tends to have an opposite phase winding between s- and d-wave condensates. We find that this peculiar feature together with the competition between s- and d-wave symmetry results in three distinct classes of vortical configurations. When either s or d condensate absolutely dominates, vortices form a conventional lattice. However, when one condensate is relatively dominant, vortices organize in chains that exhibit skyrmionic character, separating the chiral components of the s +/- id order parameter into domains within and outside the chain. Such skyrmionic chains are found stable even at high magnetic field. When s and d condensates have comparable strength, vortices split cores in two chiral components to form full-fledged skyrmions, i.e., coreless topological structures with an integer topological charge, organized in a lattice. We provide characteristic magnetic field distributions of all states, enabling their identification in, e.g., scanning Hall probe and scanning SQUID experiments. These unique vortex states are relevant for high-T-c cuprate and iron-based superconductors, where the relative strength of competing pairing symmetries is expected to be tuned by temperature and/or doping level, and can help distinguish s + is and s + id superconducting phases.

Abstract: High tunability of two-dimensional magnetic materials (by strain, gating, heterostructuring, or otherwise) provides unique conditions for studying versatile magnetic properties and controlling emergent magnetic phases. Expanding the scope of achievable magnetic phenomena in such materials is important for both fundamental and technological advances. Here we perform atomistic spin-dynamics simulations to explore the (chiral) magnetic phases of atomic monolayers in the limit of suppressed first-neighbors exchange interaction. We report the rich phase diagram of exotic magnetic configurations, obtained for both square and honeycomb lattice symmetries, comprising coexistence of ferromagnetic and antiferromagnetic spin cycloids, as well as multiple types of magnetic skyrmions. We perform a minimum-energy path analysis for the skyrmion collapse to evaluate the stability of such topological objects and reveal that magnetic monolayers could be good candidates to host the antiferromagnetic skyrmions that are experimentally evasive to date.