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 (2D) carbon materials play an important role in nanomaterials. We propose a new carbon monolayer, named hexagonal-4,4,4-graphyne (H-4,H-4,H-4-graphyne), which is a nanoporous structure composed of rectangular carbon rings and triple bonds of carbon. Using first-principles calculations, we systematically studied the structure, stability, and band structure of this new material. We found that its total energy is lower than that of experimentally synthesized beta-graphdiyne and it is stable at least up to 1500 K. In contrast to the single Dirac point band structure of other 2D carbon monolayers, the band structure of H-4,H-4,H-4-graphyne exhibits double Dirac points along the high-symmetry points and the corresponding Fermi velocities (1.04-1.27 x 10(6) m/s) are asymmetric and higher than that of graphene. The origin of these double Dirac points is traced back to the nodal line states, which can be well explained by a tight-binding model. The H-4,H-4,H-4-graphyne forms a moire superstructure when placed on top of a hexagonal boron nitride substrate. These properties make H-4,H-4,H-4-graphyne a promising semimetal material for applications in high-speed electronic devices. (C) 2018 Elsevier Ltd. All rights reserved.

Abstract: Hyperbolic polaritons in ultra thin materials such as few layers of van derWaals heterostructures provide a unique control over light-matter interaction at the nanoscale and with various applications in flat optics. Natural hyperbolic surface plasmons have been observed on thin films of WTe2 in the light wavelength range of 16-23 mu m (similar or equal to 13-18 THz) [Nat. Commun. 11, 1158 (2020)]. Using time-dependent density functional theory, it is found that carbon doped monolayer phosphorene (beta-allotrope of carbon phosphide monolayer) exhibits natural hyperbolic plasmons at frequencies above similar or equal to 5 THz which is not observed in its parent materials, i.e. monolayer of black phosphorous and graphene. Furthermore, we found that by electrostatic doping the plasmonic frequency range can be extended to the mid-infrared. (C) 2021 Elsevier Ltd. All rights reserved.

Abstract: Using the four-band continuum model we derive a general expression for the infinite-mass boundary condition in bilayer graphene. Applying this new boundary condition we analytically calculate the confined states and the corresponding wave functions in a bilayer graphene quantum dot in the absence and presence of a perpendicular magnetic field. Our results for the energy spectrum show an energy gap between the electron and hole states at small magnetic fields. Furthermore the electron (e) and hole (h) energy levels corresponding to the K and K' valleys exhibit the E-K(e(h)) (m) = E-K'(e(h)) (m) symmetry, where m is the angular momentum quantum number. (C) 2014 Elsevier Ltd. All rights reserved.

Abstract: We present the Tight-Binding Studio (TB Studio) software package that calculates the different parameters of a tight-binding Hamiltonian from a set of Bloch energy bands obtained from first principle theories such as density functional theory, Hartree-Fock calculations or semi-empirical band-structure theory. This will be helpful for scientists who are interested in studying electronic and optical properties of structures using Green's function theory within the tight-binding approximation. TB Studio is a cross-platform application written in C++ with a graphical user interface design that is user-friendly and easy to work with. This software is powered by Linear Algebra Package C interface library for solving the eigenvalue problems and the standard high performance OpenGL graphic library for real time plotting. TB Studio and its examples together with the tutorials are available for download from tight-binding.com. Program summary Program Title: Tight-Binding Studio Program Files doi:http://dx.doi.org/10.17632/j6x5mwzm2d.1 Licensing provisions: LGPL Programming language: C++ External routines: BLAS, LAPACK, LAPACKE, wxWidgets, OpenGL, MathGL Nature of problem: Obtaining Tight-Binding Hamiltonian from a set of Bloch energy bands obtained from first-principles calculations. Solution method: Starting from the simplified LCAO method, a tight-binding model in the two-center approximation is constructed. The Slater and Koster (SK) approach is used to calculate the parameters of the TB Hamiltonian. By using non-linear fitting approaches the optimal values of the SK parameters are obtained such that the TB energy eigenvalues are as close as possible to those from first-principles calculations. We obtain the expression for the Hamiltonian and the overlap matrix elements between the different orbitals of the different atoms in an orthogonal or non-orthogonal basis set. (C) 2020 Elsevier B.V. All rights reserved.

Abstract: Molecular dynamics simulations are used to investigate the effect of an AFM tip when indenting graphene nano bubbles filled by a noble gas (i.e. He, Ne and Ar) up to the breaking point. The failure points resemble those of viral shells as described by the Foppl-von Karman (FvK) dimensionless number defined in the context of elasticity theory of thin shells. At room temperature, He gas inside the bubbles is found to be in the liquid state while Ne and Ar atoms are in the solid state although the pressure inside the nano bubble is below the melting pressure of the bulk. The trapped gases are under higher hydrostatic pressure at low temperatures than at room temperature.

Abstract: We explore the flatness of conduction and valence bands of interlayer excitons in MoS2/WSe2 van der Waals heterobilayers, tuned by interlayer twist angle, pressure, and external electric field. We employ an efficient continuum model where the moire pattern from lattice mismatch and/or twisting is represented by an equivalent mesoscopic periodic potential. We demonstrate that the mismatch moire potential is too weak to produce significant flattening. Moreover, we draw attention to the fact that the quasi-particle effective masses around the Gamma-point and the band flattening are reduced with twisting. As an alternative approach, we show (i) that reducing the interlayer distance by uniform vertical pressure can significantly increase the effective mass of the moire hole, and (ii) that the moire depth and its band flattening effects are strongly enhanced by accessible electric gating fields perpendicular to the heterobilayer, with resulting electron and hole effective masses increased by more than an order of magnitude – leading to record-flat bands. These findings impose boundaries on the commonly generalized benefits of moire twistronics, while also revealing alternative feasible routes to achieve truly flat electron and hole bands to carry us to strongly correlated excitonic phenomena on demand.

Abstract: Optical and electronic properties of black phosphorus strongly depend on the number of layers and type of stacking. Using first-principles calculations within the framework of density functional theory, we investigate the electronic properties of bilayer black phosphorus with an interlayer twist angle of 90 degrees. These calculations are complemented with a simple (k) over right arrow . (p) over right arrow model which is able to capture most of the low energy features and is valid for arbitrary twist angles. The electronic spectrum of 90 degrees twisted bilayer black phosphorus is found to be x-y isotropic in contrast to the monolayer. However x-y anisotropy, and a partial return to monolayer-like behavior, particularly in the valence band, can be induced by an external out-of-plane electric field. Moreover, the preferred hole effective mass can be rotated by 90 degrees simply by changing the direction of the applied electric field. In particular, a +0.4 (-0.4) V angstrom(1) out-of-plane electric field results in a similar to 60% increase in the hole effective mass along the y (x) axis and enhances the m(y)*/m(x)* (m(x)*/m(y)*) ratio as much as by a factor of 40. Our DFT and (k) over right arrow . (p) over right arrow simulations clearly indicate that the twist angle in combination with an appropriate gate voltage is a novel way to tune the electronic and optical properties of bilayer phosphorus and it gives us a new degree of freedom to engineer the properties of black phosphorus based devices.

Abstract: The self rotation of a graphene flake over graphite is controlled by the size, initial misalignment and temperature. Using both ab initio calculations and molecular dynamics simulations, we investigate annealing effects on the self rotation of a graphene flake on a graphene substrate. The energy barriers for rotation and drift of a graphene flake over graphene is found to be smaller than 25 meV/atom which is comparable to thermal energy. We found that small flakes (of about similar to 4 nm) are more sensitive to temperature and initial misorientation angles than larger one (beyond 10 nm). The initial stacking configuration of the flake is found to be important for its dynamics and time evolution of misalignment. Large flakes, which are initially in the AA-or AB-stacking state with small misorientation angle, rotate and end up in the AB-stacking configuration. However small flakes can they stay in an incommensurate state specially when the initial misorientation angle is larger than 2 degrees. Our results are in agreement with recent experiments.

Abstract: Atomic collapse is a phenomenon inherent to relativistic quantum mechanics where electron states dive in the positron continuum for highly charged nuclei. This phenomenon was recently observed in graphene. Here we investigate a novel collapse phenomenon when multiple sub- and supercritical charges of equal strength are put close together as in a molecule. We construct a phase diagram which consists of three distinct regions: (1) subcritical, (2) frustrated atomic collapse, and (3) molecular collapse. We show that the single impurity atomic collapse resonances rearrange themselves to form molecular collapse resonances which exhibit a distinct bonding, anti-bonding and non-bonding character. Here we limit ourselves to systems consisting of two and three charges. We show that by tuning the distance between the charges and their strength a high degree of control over the molecular collapse resonances can be achieved.

Abstract: Dirac plasmons in graphene hybridize with phonons of transition metal dichalcogenides (TMDs) when the materials are combined in so-called van der Waals heterostructures (vdWh), thus forming surface plasmon-phonon polaritons (SPPPs). The extend to which these modes are coupled depends on the TMD composition and structure, but also on the plasmons' properties. By performing realistic simulations that account for the contribution of each layer of the vdWh separately, we calculate how the strength of plasmon-phonon coupling depends on the number and composition of TMD layers, on the graphene Fermi energy and the specific phonon mode. From this, we present a semiclassical theory that is capable of capturing all relevant characteristics of the SPPPs. We find that it is possible to realize both strong and ultra-strong coupling regimes by tuning graphene's Fermi energy and changing TMD layer number.

Abstract: Two-dimensional atomic crystals can radically change their properties in response to external influences, such as substrate orientation or strain, forming materials with novel electronic structure(1-5). An example is the creation of weakly dispersive, 'flat' bands in bilayer graphene for certain 'magic' angles of twist between the orientations of the two layers(6). The quenched kinetic energy in these flat bands promotes electron-electron interactions and facilitates the emergence of strongly correlated phases, such as superconductivity and correlated insulators. However, the very accurate fine-tuning required to obtain the magic angle in twisted-bilayer graphene poses challenges to fabrication and scalability. Here we present an alternative route to creating flat bands that does not involve fine-tuning. Using scanning tunnelling microscopy and spectroscopy, together with numerical simulations, we demonstrate that graphene monolayers placed on an atomically flat substrate can be forced to undergo a buckling transition(7-9), resulting in a periodically modulated pseudo-magnetic field(10-14), which in turn creates a 'post-graphene' material with flat electronic bands. When we introduce the Fermi level into these flat bands using electrostatic doping, we observe a pseudogap-like depletion in the density of states, which signals the emergence of a correlated state(15-17). This buckling of two-dimensional crystals offers a strategy for creating other superlattice systems and, in particular, for exploring interaction phenomena characteristic of flat bands. Buckled monolayer graphene superlattices are found to provide an alternative to twisted bilayer graphene for the study of flat bands and correlated states in a carbon-based material.

Abstract: The basal plane of graphene can function as a selective barrier that is permeable to protons but impermeable to all ions and gases, stimulating its use in applications such as membranes, catalysis and isotope separation. Protons can chemically adsorb on graphene and hydrogenate it, inducing a conductor–insulator transition that has been explored intensively in graphene electronic devices. However, both processes face energy barriersand various strategies have been proposed to accelerate proton transport, for example by introducing vacancies, incorporating catalytic metalsor chemically functionalizing the lattice. But these techniques can compromise other properties, such as ion selectivity or mechanical stability. Here we show that independent control of the electric field,<italic>E</italic>, at around 1 V nm⁻¹, and charge-carrier density,<italic>n</italic>, at around 1 × 10¹⁴ cm⁻², in double-gated graphene allows the decoupling of proton transport from lattice hydrogenation and can thereby accelerate proton transport such that it approaches the limiting electrolyte current for our devices. Proton transport and hydrogenation can be driven selectively with precision and robustness, enabling proton-based logic and memory graphene devices that have on–off ratios spanning orders of magnitude. Our results show that field effects can accelerate and decouple electrochemical processes in double-gated 2D crystals and demonstrate the possibility of mapping such processes as a function of<italic>E</italic>and<italic>n</italic>, which is a new technique for the study of 2D electrode–electrolyte interfaces.

Abstract: Quantum spin Hall (QSH) insulators with a large topologically nontrivial bulk gap are crucial for future applications of the QSH effect. Among these, group III-V monolayers and their halides, which have a chair structure (regular hexagonal framework), have been widely studied. Using first-principles calculations, we formulate a new structure model for the functionalized group III-V monolayers, which consist of rectangular GaBi-X-2 (X = I, Br, Cl) monolayers with a distorted hexagonal framework (DHF). These structures have a far lower energy than the GaBi-X-2 monolayers with a chair structure. Remarkably, the DHF GaBi-X-2 monolayers are all QSH insulators, which exhibit sizeable nontrivial band gaps ranging from 0.17 to 0.39 eV. The band gaps can be widely tuned by applying different spin-orbit coupling strengths, resulting in a distorted Dirac cone.

Abstract: The quantum anomalous Hall (QAH) effect is a topologically nontrivial phase, characterized by a non-zero Chern number defined in the bulk and chiral edge states in the boundary. Using first-principles calculations, we demonstrate the presence of the QAH effect in a 1T-YN2 monolayer, which was recently predicted to be a Dirac half metal without spin-orbit coupling (SOC). We show that the inclusion of SOC opens up a large nontrivial bandgap of nearly 0.1 eV in the electronic band structure. This results in the nontrivial bulk topology, which is confirmed by the calculation of Berry curvature, anomalous Hall conductance and the presence of chiral edge states. Remarkably, a QAH phase of high Chern number C = 3 is found, and there are three corresponding gapless chiral edge states emerging inside the bulk gap. Different substrates are also chosen to study the possible experimental realization of the 1T-YN2 monolayer, while retaining its nontrivial topological properties. Our results open a new avenue in searching for QAH insulators with high temperature and high Chern numbers, which can have nontrivial practical applications.