Abstract: Edge functionalization of graphene nanoribbons with nitrogen atoms for various adatom configurations at armchair and zigzag edges are investigated. We provide comprehensive information on the electronic and magnetic properties and investigate the stability of the various systems. Two types of rippling of the nanoribbons, namely edge and bulk rippling depending on the sign of edge stress induced at the edge, are found. They are found to play the decisive role for the stability of the structures. We also propose a type of edge decoration in which every third nitrogen adatom at the zigzag edges is replaced by an oxygen atom. In this way, the electron count is compatible with a full aromatic structure, leading to additional stability and a disappearance of magnetism that is usually associated with zigzag nanoribbons.

Abstract: Massless Dirac fermions (MDFs) emerge as quasiparticles in various novel materials such as graphene and topological insulators, and they exhibit several intriguing properties, of which Veselago focusing is an outstanding example with a lot of possible applications. However, up to now Veselago focusing merely occurred in p-n junction devices based on the isotropic MDF, which lacks the tunability needed for realistic applications. Here, motivated by the emergence of novel Dirac materials, we investigate the propagation behaviors of anisotropic MDFs in such a p-n junction structure. By projecting the Hamiltonian of the anisotropic MDF to that of the isotropic MDF and deriving an exact analytical expression for the propagator, precise Veselago focusing is demonstrated without the need for mirror symmetry of the electron source and its focusing image. We show a tunable focusing position that can be used in a device to probe masked atom-scale defects. This study provides an innovative concept to realize Veselago focusing relevant for potential applications, and it paves the way for the design of novel electron optics devices by exploiting the anisotropic MDF.

Abstract: Excitonic superfluidity in double phosphorene monolayers is investigated using the BCS mean-field equations. Highly anisotropic superfluidity is predicted where we found that the maximum superfluid gap is in the Bose-Einstein condensate (BEC) regime along the armchair direction and in the BCS-BEC crossover regime along the zigzag direction. We estimate the highest Kosterlitz-Thouless transition temperature with maximum value up to similar to 90 K with onset carrier densities as high as 4 x 10(12) cm(-2). This transition temperature is significantly larger than what is found in double electron-hole few-layers graphene. Our results can guide experimental research toward the realization of anisotropic condensate states in electron-hole phosphorene monolayers.

Abstract: Water inside a nanocapillary exhibits unconventional structural and dynamical behavior due to its ordered structure. The confining walls, density, and lateral pressures control profoundly the microscopic structure of trapped water. Here we study the electrostriction of confined water subjected to pressures of the order of GPa for two different setups: (i) a graphene nanochannel containing a constant number of water molecules independent of the height of the channel, (ii) an open nanochannel where water molecules can be exchanged with those in a reservoir. For the former case, a square-rhombic structure of confined water is formed when the height of the channel is d = 6.5 angstrom having a density of rho = 1.42 g cm(-3). By increasing the height of the channel, a transition from a flat to a buckled state occurs, whereas the density rapidly decreases and reaches the bulk density for d congruent to 8.5 angstrom. When a perpendicular electric field is applied, the water structure and the lateral pressure change. For strong electric fields (similar to 1 V/angstrom), the square-rhombic structure is destroyed. For an open setup, a solid phase of confined water consisting of an imperfect square-rhombic structure is formed. By applying a perpendicular field, the density and phase of confined water change. However, the density and pressure inside the channel decrease as compared to the first setup. Our study is closely related to recent experiments on confined water, and it reveals the sensitivity of the microscopic structure of confined water to the size of the channel, the external electric field, and the experimental setup.

Abstract: The role of dielectric screening of electron-hole interaction in van der Waals heterostructures is theoretically investigated. A comparison between models available in the literature for describing these interactions is made and the limitations of these approaches are discussed. A simple numerical solution of Poisson's equation for a stack of dielectric slabs based on a transfer matrix method is developed, enabling the calculation of the electron-hole interaction potential at very low computational cost and with reasonable accuracy. Using different potential models, direct and indirect exciton binding energies in these systems are calculated within Wannier-Mott theory, and a comparison of theoretical results with recent experiments on excitons in two-dimensional materials is discussed.

Abstract: By performing density functional theory-based ab initio calculations, Raman-active phonon modes of single-layer two-dimensional (2D) materials and the effect of in-plane biaxial strain on the peak frequencies and corresponding activities of the Raman-active modes are calculated. Our findings confirm the Raman spectrum of the unstrained 2D crystals and provide expected variations in the Raman-active modes of the crystals under in-plane biaxial strain. The results are summarized as follows: (i) frequencies of the phonon modes soften (harden) under applied tensile (compressive) strains; (ii) the response of the Raman activities to applied strain for the in-plane and out-of-plane vibrational modes have opposite trends, thus, the built-in strains in the materials can be monitored by tracking the relative activities of those modes; (iii) in particular, the A peak in single-layer Si and Ge disappears under a critical tensile strain; (iv) especially in mono-and diatomic single layers, the shift of the peak frequencies is a stronger indication of the strain rather than the change in Raman activities; (v) Raman-active modes of single-layer ReX2 (X = S, Se) are almost irresponsive to the applied strain. Strain-induced modifications in the Raman spectrum of 2D materials in terms of the peak positions and the relative Raman activities of the modes could be a convenient tool for characterization.

Abstract: Using the tight-binding approach, we investigate the energy spectrum of square, triangular, and hexagonal MoS2 quantum dots (QDs) in the presence of a perpendicular magnetic field. Novel edge states emerge in MoS2 QDs, which are distributed over the whole edge which we call ring states. The ring states are robust in the presence of spin-orbit coupling (SOC). The corresponding energy levels of the ring states oscillate as a function of the perpendicular magnetic field which are related to Aharonov-Bohm oscillations. Oscillations in the magnetic field dependence of the energy levels and the peaks in the magneto-optical spectrum emerge (disappear) as the ring states are formed (collapsed). The period and the amplitude of the oscillation decrease with the size of the MoS2 QDs.

Abstract: Defects are almost inevitably present in realistic materials and defective materials are expected to exhibit very different properties than their nondefective (perfect) counterparts. Here, using a combination of the tight-binding approach and the scattering matrix formalism, we investigate the electronic transport properties of defective phosphorene nanoribbons (PNRs) containing atomic vacancies. We find that for both armchair PNRs (APNRs) and zigzag PNRs (ZPNRs), single vacancies can create quasilocalized states, which can affect their conductance. With increasing vacancy concentration, three different transport regimes are identified: ballistic, diffusive, and Anderson localized ones. In particular, ZPNRs that are known to be metallic due to the presence of edge states become semiconducting: edge conductance vanishes and transport gap appears due to Anderson localization. Moreover, we find that for a fixed vacancy concentration, both APNRs and ZPNRs of narrower width and/or longer length are more sensitive to vacancy disorder than their wider and/or shorter counterparts, and that for the same ribbon length and width, ZPNRs are more sensitive to vacancy disorder than APNRs.

Abstract: <script type='text/javascript'>document.write(unpmarked('We show that delaminations in bilayer graphene (BLG) with electrostatically induced interlayer symmetry can provide one with ballistic channels for electrons with energies inside the electrostatically induced BLG gap. These channels are formed by a combination of valley-polarized evanescent states propagating along the delamination edges (which persist in the presence of a strong magnetic field) and standing waves bouncing between them inside the delaminated region (in a strong magnetic field, these transform into Landau levels in the monolayers). For inverted stackings in BLGs on the left and right of the delamination (AB-2ML-BA or BA-2ML-AB, where 2ML indicates two decoupled monolayers of graphene), the lowest-energy ballistic channels are gapless, have linear dispersion, and appear to be weakly topologically protected. When BLG stackings on both sides of the delamination are the same (AB-2ML-AB or BA-2ML-BA), the lowest-energy ballistic channels are gapped, with a gap epsilon(g) scaling as epsilon(g) alpha W-1 with delamination width and epsilon(g) alpha delta(-1) with the on-layer energy difference in the delaminated part of the structure. Depending on the width, delaminations may also support several \u0022higher-energy\u0022 waveguide modes. Our results are based on both the analytical study of the wave matching of Dirac states and tight-binding model calculations, and we analyze in detail the dependence of the delamination spectrum on the electrostatic conditions in the structure, such as the vertical displacement field.'));

Abstract: <script type='text/javascript'>document.write(unpmarked('We present an advanced method to study spin fluctuations in superconductors quantitatively and entirely from first principles. This method can be generally applied to materials where electron-phonon coupling and spin fluctuations coexist. We employ it here to examine the recently synthesized superconductor iron tetraboride (FeB4) with experimental T-c similar to 2.4 K [H. Gou et al., Phys. Rev. Lett, 111, 157002 (2013)]. We prove that FeB4 is particularly prone to ferromagnetic spin fluctuations due to the presence of iron, resulting in a large Stoner interaction strength, I = 1.5 eV, as calculated from first principles. The other important factor is its Fermi surface that consists of three separate sheets, among which two are nested ellipsoids. The resulting susceptibility has a ferromagnetic peak around q = 0, from which we calculated the repulsive interaction between Cooper pair electrons using the random phase approximation. Subsequently, we combined the electron-phonon interaction calculated from first principles with the spin fluctuation interaction in fully anisotropic Eliashberg theory calculations. We show that the resulting superconducting gap spectrum is conventional, yet very strongly depleted due to coupling to the spin fluctuations. The critical temperature decreases from T-c = 41 K, if they are not taken into account, to T-c = 1.7 K, in good agreement with the experimental value.'));

Abstract: <script type='text/javascript'>document.write(unpmarked('Spin-triplet chiral p-wave superconductivity is typically described by a two-component order parameter, and as such is prone to unique emergent effects when compared to the standard single-component superconductors. Here we present the equilibrium phase diagram for small mesoscopic chiral p-wave superconducting disks in the presence of magnetic field, obtained by solving the microscopic Bogoliubov-de Gennes equations self-consistently. In the ultrasmall limit, the cylindrically symmetric giant-vortex states form the ground state of the system. However, with increasing sample size, the cylindrical symmetry is broken as the two components of the order parameter segregate into domains, and the number of fragmented domain walls between them characterizes the resulting states. Such domain walls are topological defects unique for the p-wave order, and constitute a dominant phase in the mesoscopic regime. Moreover, we find two possible types of domain walls, identified by their chirality-dependent interaction with the edge states.'));

Abstract: <script type='text/javascript'>document.write(unpmarked('gamma-Fe and related alloys are model systems of the coupling between structure and magnetism in solids. Since different electronic states (with different volumes and magnetic ordering states) are closely spaced in energy, small perturbations can alter which one is the actual ground state. Here, we demonstrate that the ferromagnetic state of gamma-Fe nanoparticles is associated with a tetragonal distortion of the fcc structure. Combining a wide range of complementary experimental techniques, including low-temperature Mossbauer spectroscopy, advanced transmission electron microscopy, and synchrotron radiation techniques, we unambiguously identify the tetragonally distorted ferromagnetic ground state, with lattice parameters a = 3.76(2) angstrom and c = 3.50(2) angstrom, and a magnetic moment of 2.45(5) mu(B) per Fe atom. Our findings indicate that the ferromagnetic order in nanostructured gamma-Fe is generally associated with a tetragonal distortion. This observation motivates a theoretical reassessment of the electronic structure of gamma-Fe taking tetragonal distortion into account.'));

Abstract: The energy spectrum and eigenstates of single-layer black phosphorus nanoribbons in the presence of a perpendicular magnetic field and an in-plane transverse electric field are investigated by means of a tight-binding method, and the effect of different types of edges is examined analytically. A description based on a continuum model is proposed using an expansion of the tight-binding model in the long-wavelength limit. Thewave functions corresponding to the flatband part of the spectrum are obtained analytically and are shown to agree well with the numerical results from the tight-binding method for both narrow (10 nm) and wide (100 nm) nanoribbons. Analytical expressions for the critical magnetic field at which Landau levels are formed and the ranges of wave numbers in the dispersionless flatband segments in the energy spectra are derived. We examine the evolution of the Landau levels when an in-plane lateral electric field is applied, and we determine analytically how the edge states shift withmagnetic field. For wider nanoribbons, the conductance is shown to have a characteristic staircase shape in combined magnetic and electric fields. Some of the stairs in zigzag and skewed armchair nanoribbons originate from edge states that are found in the band gap.

Abstract: Using the tight-binding approach coupled with mean-field Hubbard model, we theoretically study the effect of mechanical deformations on the magnetic properties of bilayer graphene (BLG) quantum dots (QDs). Results are obtained for AA-and AB(Bernal)-stacked BLG QDs, considering different geometries (hexagonal, triangular and square shapes) and edge types (armchair and zigzag edges). In the absence of strain, our results show that (i) the magnetization is affected by taking different dot sizes only for hexagonal BLG QDs with zigzag edges, exhibiting different critical Hubbard interactions, and (ii) the magnetization does not depend on the interlayer hopping energies, except for the geometries with zigzag edges and AA stacking. In the presence of in-plane and uniaxial strain, for all geometries we obtain two different magnetization regimes depending on the applied strain amplitude. The appearance of such different regimes is due to the breaking of layer and sublattice symmetries in BLG QDs.

Abstract: Recent experiments show that an atomic step on the surface of atomically thin metallic films can strongly affect electronic transport. Here we reveal multiple and versatile effects that such a surface step can have on superconductivity in ultrathin films. By solving the Bogoliubov-de Gennes equations self-consistently in this regime, where quantum confinement dominates the emergent physics, we show that the electronic structure is profoundly modified on the two sides of the step, as is the spatial distribution of the superconducting order parameter and its dependence on temperature and electronic gating. Furthermore, the surface step changes nontrivially the transport properties both in the proximity-induced superconducting pair correlations and the Josephson effect, depending on the step height. These results offer a new route to tailor superconducting circuits and design atomically thin heterojunctions made of one same material.