Abstract: Atomic collapse in graphene nanoribbons behaves in a fundamentally different way as compared to monolayer graphene due to the presence of multiple energy bands and the effect of edges. For armchair nanoribbons we find that bound states gradually transform into atomic collapse states with increasing impurity charge. This is very different in zigzag nanoribbons where multiple quasi-one-dimensional bound states are found that originates from the zero-energy zigzag edge states. They are a consequence of the flat band and the electron distribution of these bound states exhibits two peaks. The lowest-energy edge state transforms from a bound state into an atomic collapse resonance and shows a distinct relocalization from the edge to the impurity position with increasing impurity charge.

Abstract: We examine from first principles the flexoelectric properties of phosphorene nanoribbons under mechanical bending along armchair and zigzag directions. In both cases we find that the radial polarization depends linearly on the strain gradient. The flexoelectricity along the armchair direction is over 40% larger than along the zigzag direction. The obtained flexoelectric coefficients of phosphorene are four orders of magnitude larger than those of graphene and comparable to transition metal dichalcogenides. Analysis of charge density shows that the flexoelectricity mainly arises from the pz orbitals of phosphorus atoms. The electron mobilities in bent phosphorene can be enhanced by over 60% along the armchair direction, which is significantly higher than previous reports of mobility tuned by uniaxial strain. Our results indicate phosphorene is a candidate for a two-dimensional material applicable in flexible-electronic devices.

Abstract: Structural symmetry breaking in two-dimensional materials can lead to superior physical properties and introduce an additional degree of piezoelectricity. In the present paper, we propose three structural phases (1H, 1T, and 1T') of Janus WXO (X = S, Se, and Te) monolayers and investigate their vibrational, thermal, elastic, piezoelectric, and electronic properties by using first-principles methods. Phonon spectra analysis reveals that while the 1H phase is dynamically stable, the 1T phase exhibits imaginary frequencies and transforms to the distorted 1T' phase. Ab initio molecular dynamics simulations confirm that 1H- and 1T'-WXO monolayers are thermally stable even at high temperatures without any significant structural deformations. Different from binary systems, additional Raman active modes appear upon the formation of Janus monolayers. Although the mechanical properties of 1H-WXO are found to be isotropic, they are orientation dependent for 1T'-WXO. It is also shown that 1H-WXO monolayers are indirect band-gap semiconductors and the band gap narrows down the chalcogen group. Except 1T'-WSO, 1T'-WXO monolayers have a narrow band gap correlated with the Peierls distortion. The effect of spin-orbit coupling on the band structure is also examined for both phases and the alteration in the band gap is estimated. The versatile mechanical and electronic properties of Janus WXO monolayers together with their large piezoelectric response imply that these systems are interesting for several nanoelectronic applications.

Abstract: The discovery of the hydrodynamic electron liquid (HEL) in graphene [D. Bandurin et al., Science 351, 1055 (2016) and J. Crossno et al., Science 351, 1058 (2016)] has marked the birth of the solid-state HEL which can be probed near room temperature in a table-top setup. Here we examine the terahertz (THz) magneto-optical (MO) properties of a graphene HEL. Considering the case where the magnetic length l(B) = root h/eB is comparable to the mean-free path l(ee) for electron-electron interaction in graphene, the MO conductivities are obtained by taking a momentum balance equation approach on the basis of the Boltzmann equation. We find that when l(B) similar to l(ee), the viscous effect in a HEL can weaken significantly the THz MO effects such as cyclotron resonance and Faraday rotation. The upper hybrid and cyclotron resonance magnetoplasmon modes omega(+/-) are also obtained through the RPA dielectric function. The magnetoplasmons of graphene HEL at large wave-vector regime are affected by the viscous effect, and results in red-shifts of the magnetoplasmon frequencies. We predict that the viscosity in graphene HEL can affect strongly the magneto-optical and magnetoplasmonic properties, which can be verified experimentally.

Abstract: The interplay between a magnetic field and the Coulomb potential from a charged vacancy on the electron states in graphene is investigated within the tight-binding model. The Coulomb potential removes locally Landau level degeneracy, while the vacancy introduces a satellite level next to the normal Landau level. These satellite levels are found throughout the positive-energy region, but in the negative-energy region, they turn into atomic collapse resonances. Crossings between Landau levels with different angular quantum number m are found. Unlike the point impurity system in which an anticrossing occurs between Landau levels of the same m, in this work anticrossing is found between the normal Landau level and the vacancy-induced level. The atomic collapse resonance hybridizes with the Landau levels. The charge at which the lowest Landau level m = -1, N = 1 crosses E = 0 increases with enhancing magnetic field. A Landau level scaling anomaly occurs when the charge is larger than the critical charge beta 0.6 and this critical charge is independent of the magnetic field.

Abstract: Transport measurements are readily used to probe different phases in disordered topological insulators (TIs), where determining topological invariants explicitly is challenging. On that note, universal conductance fluctuations (UCF) theory asserts the conductance G for an ensemble has a Gaussian distribution, and that standard deviation 8G depends solely on the symmetries and dimensions of the system. Using a real-space tight -binding Hamiltonian on a system with Anderson disorder, we explore conductance fluctuations in a thin Bi2Se3 film and demonstrate the agreement of their behavior with UCF hypotheses. We further show that magnetic field applied out-of-plane breaks the time -reversal symmetry and transforms the system's Wigner-Dyson class from root symplectic to unitary, increasing 8G by 2. Finally, we reveal that while Bi2Se3 is a strong TI, weak TI and metallic phases can be stabilized in presence of strain and disorder, and detected by monitoring the conductance fluctuations.

Abstract: Flat bands in the energy spectrum have attracted a lot of attention in recent years because of their unique properties and promising applications. Special arrangement of impurities on monolayer graphene are proposed to generate multiflat bands in the electronic band structure. In addition to the single midgap states in the spectrum of graphene with low hydrogen density, we found closely spaced bands around the Fermi level with increasing impurity density, which are similar to discrete lines in the spectrum of quantum dots, as well as the unusual Landau-level energy spectrum of graphene in the presence of a strong magnetic field. The presence of flat bands crucially depends on whether or not there are odd or even electrons of H(F) atoms bound to graphene. Interestingly, we found that a fully hydrogenated (fluoridated) of a hexagon of graphene sheet with six hydrogen (fluorine) atoms sitting on top and bottom in consecutive order exhibits Dirac cones in the electronic band structure with a 20% smaller Fermi velocity as compared to the pristine graphene. Functionalizing graphene introduces various C-C bond lengths resulting in nonuniform strains. Such a nonuniform strain may induce a giant pseudomagnetic field in the system, resulting in quantum Hall effect.

Abstract: The surface termination of MXenes greatly determines the electrochemical properties and ion kinetics on their surfaces. So far, hydroxyl-, oxygen-, and fluorine-terminated MXenes have been widely studied for energy storage applications. Recently, sulfur-functionalized MXene structures, which possess low diffusion barriers, have been proposed as candidate materials to enhance battery performance. We performed first-principles calculations on the structural, stability, electrochemical, and ion dynamic properties of Li-adsorbed sulfur-functionalized groups 3B, 4B, 5B, and 6B transition-metal (M)-based MXenes (i.e., M2CS2 with M = Sc, Ti, Zr, Hf, V, Nb, Ta, Cr, Mo, and W). We performed phonon calculations, which indicated that all of the above M2CS2 MXenes, except for Sc, are dynamically stable at T = 0 K. The ground-state structure of each M2CS2 monolayer depends on the type of M atom. For instance, while sulfur prefers to sit at the FCC site on Ti2CS2, it occupies the HCP site of Cr-based MXene. We determined the Li adsorption configurations at different concentrations using the cluster expansion method. The highest maximum open-circuit voltages were computed for the group 4B element (i.e., Ti, Zr, and Hf)-based M2CS2, which are larger than 2.1 V, while their average voltages are approximately 1 V. The maximum voltage for the group 6B element (i.e., Cr, Mo, W)-based M2CS2 is less than 1 V, and the average voltage is less than 0.71 V. We found that S functionalization is helpful for capacity improvements over the O-terminated MXenes. In this respect, the computed storage gravimetric capacity may reach up to 417.4 mAh/g for Ti2CS2 and 404.5 mAh/g for V2CS2. Ta-, Cr-, Mo-, and W-based M2CS2 MXenes show very low capacities, which are less than 100 mAh/g. The Li surface diffusion energy barriers for all of the considered MXenes are less than 0.22 eV, which is favorable for high charging and discharging rates. Finally, ab initio molecular dynamic simulations performed at 400 K and bond-length analysis with respect to Li concentration verify that selected promising systems are robust against thermally induced perturbations that may induce structural transformations or distortions and undesirable Li release.

Abstract: The hunt for high-temperature superfluidity has received new impetus from the discovery of atomically thin stable materials. Electron-hole superfluidity in coupled MoSe2-WSe2 monolayers is investigated using a mean-field multiband model that includes band splitting caused by strong spin-orbit coupling. This splitting leads to a large energy misalignment of the electron and hole bands which is strongly modified by interchanging the doping of the monolayers. The choice of doping determines if the superfluidity is tunable from one to two components.

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: 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: Using density functional theory based first-principles calculations, we predict the dynamically stable 1H phase of a Janus single layer composed of S-Mo-O atomic layers. It is an indirect band gap semiconductor exhibiting strong polarization arising from the charge difference on the two surfaces. In contrast to 1H phases of MoS2 and MoO2, Janus MoSO is found to possess four Raman active phonon modes and a large out-of-plane piezoelectric coefficient which is absent in fully symmetric single layers of MoS2 and MoO2. We investigated the electronic and phononic properties under applied biaxial strain and found an electronic phase transition with tensile strain while the conduction band edge displays a shift when under compressive strain. Furthermore, single-layer MoSO exhibits phononic stability up to 5% of compressive and 11% of tensile strain with significant phonon shifts. The phonon instability is shown to arise from the soft in-plane and out-of-plane acoustic modes at finite wave vector. The large strain tolerance of Janus MoSO is important for nanoelastic applications. In view of the dynamical stability even under moderate strain, we expect that Janus MoSO can be fabricated in the common 1H phase with a strong out-of-plane piezoelectric coefficient.

Abstract: Although there is strong theoretical and experimental evidence for electron-hole superfluidity in separated sheets of electrons and holes at low T, extending superfluidity to high T is limited by strong two-dimensional fluctuations and Kosterlitz-Thouless effects. We show this limitation can be overcome using a superlattice of alternating electron- and hole-doped semiconductor monolayers. The superfluid transition in a three-dimensional superlattice is not topological, and for strong electron-hole pair coupling, the transition temperature T-c can be at room temperature. As a quantitative illustration, we show T-c can reach 270 K for a superfluid in a realistic superlattice of transition metal dichalcogenide monolayers.

Abstract: The strain fields of periodically buckled graphene induce a periodic pseudomagnetic field (PMF) that modifies the electronic band structure. From the geometry, amplitude, and period of the periodic pseudomagnetic field, we determine the necessary conditions to access the regime of correlated phases by examining the band flattening. As compared to twisted bilayer graphene the proposed system has the advantages that (1) only a single layer of graphene is needed, (2) one is not limited to hexagonal superlattices, and (3) narrower flat bandwidth and larger separation between flat bands can be induced. We, therefore, propose that periodically strained graphene single layers can become a platform for the exploration of exotic many-body phases.

Abstract: We analyzed pattern formation and identified various morphologies in a system of particles interacting through a non-monotonic potential with a competing range interaction characterized by a repulsive core (r < r(c)) and an attractive tail (r > r(c)), using molecular-dynamics simulations. Depending on parameters, the interaction potential models the inter-particle interaction in various physical systems ranging from atoms, molecules and colloids to vortices in low kappa type-II superconductors and in recently discovered 'type-1.5' superconductors. We constructed a 'morphology diagram' in the plane 'critical radius r(c)-density n' and proposed a new approach to characterizing the different types of patterns. Namely, we elaborated a set of quantitative criteria in order to identify the different pattern types, using the radial distribution function (RDF), the local density function and the occupation factor.