Abstract: We propose systems that allow a tuning of the phonon transmission function T(omega) in graphene nanoribbons by using C-13 isotope barriers, antidot structures, and distinct boundary conditions. Phonon modes are obtained by an interatomic fifth-nearest neighbor force-constant model (5NNFCM) and T(omega) is calculated using the non-equilibrium Green's function formalism. We show that by imposing partial fixed boundary conditions it is possible to restrict contributions of the in-plane phonon modes to T(omega) at low energy. On the contrary, the transmission functions of out-of-plane phonon modes can be diminished by proper antidot or isotope arrangements. In particular, we show that a periodic array of them leads to sharp dips in the transmission function at certain frequencies omega(nu) which can be pre-defined as desired by controlling their relative distance and size. With this, we demonstrated that by adequate engineering it is possible to govern the magnitude of the ballistic transmission functions T(omega) in graphene nanoribbons. We discuss the implications of these results in the design of controlled thermal transport at the nanoscale as well as in the enhancement of thermo-electric features of graphene-based materials.

Abstract: Superconductors, ideally diamagnetic when in the Meissner state, can also exhibit paramagnetic behavior due to trapped magnetic flux. In the absence of pinning such paramagnetic response is weak, and ceases with increasing sample thickness. Here we show that in multiband superconductors paramagnetic response can be observed even in slab geometries, and can be far larger than any previous estimate – even multiply larger than the diamagnetic Meissner response for the same applied magnetic field. We link the appearance of this giant paramagnetic response to the broad crossover between conventional Type-I and Type-II superconductors, where Abrikosov vortices interact non-monotonically and multibody effects become important, causing unique flux configurations and their locking in the presence of surfaces.

Abstract: The properties of a new class of phosphorene nanotubes (PNT) are investigated by performing first-principles calculations. We demonstrate that it is advantageous to use blue phosphorene in order to make small nanotubes and propose a way to create low-energy PNTs by the inclusion of defect lines. Five different types of defect lines are investigated and incorporated in various combinations. The resulting defect-induced faceted PNTs have negligible bending stresses which leads to a reduction in the formation energy with respect to round PNTs. Our armchair faceted PNTs have similar formation energies than the recently proposed multiphase faceted PNTs, but they have a larger variety of possible structures. Our zigzag faceted PNTs have lower formation energies than round tubes and multiphase faceted nanotubes. The electronic properties of the defect-induced faceted PNTs are determined by the defect lines which control the band gap and the shape of the electronic states at the band edges. These band gaps increase with the radius of the nanotubes and converge to those of isolated defect lines.

Abstract: A tight-binding approach based on the Chebyshev-Bogoliubov-de Gennes method is used to describe disordered single-layer graphene Josephson junctions. Scattering by vacancies, ripples, or charged impurities is included. We compute the Josephson current and investigate the nature of multiple Andreev reflections, which induce bound states appearing as peaks in the density of states for energies below the superconducting gap. In the presence of single-atom vacancies, we observe a strong suppression of the supercurrent, which is a consequence of strong intervalley scattering. Although lattice deformations should not induce intervalley scattering, we find that the supercurrent is still suppressed, which is due to the presence of pseudomagnetic barriers. For charged impurities, we consider two cases depending on whether the average doping is zero, i.e., existence of electron-hole puddles, or finite. In both cases, short-range impurities strongly affect the supercurrent, similar to the vacancies scenario.

Abstract: The feasibility of Veselago lensing in graphene with a p-n junction is investigated numerically for realistic injection leads. Two different set-ups with two narrow leads are considered with absorbing or reflecting side edges. This allows us to separately determine the influence of scattering on electron focusing for the edges and the p-n interface. Both semiclassical and tight-binding simulations show a distinctive peak in the transmission probability that is attributed to the Veselago lensing effect. We investigate the robustness of this peak on the width of the injector, the position of the p-n interface, and different gate potential profiles. Furthermore, the influence of scattering by both short- and long-range impurities is considered.

Abstract: We study the effect of mechanical deformations on the elec- tronic properties of hexagonal flakes of bilayer graphene. The behavior of electrons induced by triaxial strain can be de- scribed by an effective pseudo-magnetic field which is homo- geneous in the center of the flake. We find that in-plane strain, applied to both layers equally, can break the layer symmetry leading to different behavior in the top and bottom layers of graphene. At low energy, just one of the layers feels
the pseudo-magnetic field: the zero-energy pseudo-Landau level is missing in the second layer, thus creating a gap be- tween the lowest non-zero levels. While the layer asymmetry is most significant at zero energy, interaction with the edges of the flake extends the effect to higher pseudo-Landau lev- els. The behavior of the top and bottom layers may be re- versed by rotating the triaxial strain by 60°.

Abstract: Using the tight-binding approach, we study the influence of a perpendicular magnetic field on the energy levels of hexagonal, triangular, and circular bilayer graphene (BLG) quantum dots (QDs) with zigzag and armchair edges. We obtain the energy levels for AB (Bernal)-stacked BLG QDs in both the absence and the presence of a perpendicular electric field (i.e., biased BLG QDs). We find different regions in the spectrum of biased QDs with respect to the crossing point between the lowest-electron and -hole Landau levels of a biased BLG sheet. Those different regions correspond to electron states that are localized at the center, edge, or corner of the BLG QD. Quantum Hall corner states are found to be absent in circular BLG QDs. The spatial symmetry of the carrier density distribution is related to the symmetry of the confinement potential, the position of zigzag edges, and the presence or absence of interlayer inversion symmetry.

Abstract: The rate of water flow through hydrophobic nanocapillaries is greatly enhanced as compared to that expected from macroscopic hydrodynamics. This phenomenon is usually described in terms of a relatively large slip length, which is in turn defined by such microscopic properties as the friction between water and capillary surfaces and the viscosity of water. We show that the viscosity of water and, therefore, its flow rate are profoundly affected by the layered structure of confined water if the capillary size becomes less than 2 nm. To this end, we study the structure and dynamics of water confined between two parallel graphene layers using equilibrium molecular dynamics simulations. We find that the shear viscosity is not only greatly enhanced for subnanometer capillaries, but also exhibits large oscillations that originate from commensurability between the capillary size and the size of water molecules. Such oscillating behavior of viscosity and, consequently, the slip length should be taken into account in designing and studying graphene-based and similar membranes for desalination and filtration.

Abstract: Depending on the Ginzburg-Landau parameter kappa, superconductors can either be fully diamagnetic if kappa < 1/root 2 (type I superconductors) or allow magnetic flux to penetrate through Abrikosov vortices if kappa > 1/root 2 (type II superconductors; refs 1,2). At the Bogomolny critical point, kappa = kappa(c) = 1/root 2, a state that is infinitely degenerate with respect to vortex spatial configurations arises(3,4). Despite in-depth investigations of conventional type I and type II superconductors, a thorough understanding of the magnetic behaviour in the near-Bogomolny critical regime at kappa similar to kappa(c) remains lacking. Here we report that in confined systems the critical regime expands over a finite interval of kappa forming a critical superconducting state. We show that in this state, in a sample with dimensions comparable to the vortex core size, vortices merge into a multi-quanta droplet, which undergoes Rayleigh instability(5) on increasing kappa and decays by emitting single vortices. Superconducting vortices realize Nielsen-Olesen singular solutions of the Abelian Higgs model, which is pervasive in phenomena ranging from quantum electrodynamics to cosmology(6-9). Our study of the transient dynamics of Abrikosov-Nielsen-Olesen vortices in systems with boundaries promises access to non-trivial effects in quantum field theory by means of bench-top laboratory experiments.

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: Using a numerical tight-binding approach based on the Chebyshev–Bogoliubov–de Gennes method we describe Josephson junctions made of multilayer graphene contacted by top superconducting gates. Both Bernal (ABA) and rhombohedral (ABC) stacking are considered and we find that the type of stacking has a strong effect on the proximity effect and the supercurrent flow. For both cases the pair amplitude shows a polarization between dimer and nondimer atoms, being more pronounced for rhombohedral stacking. Even though the proximity effect in nondimer sites is enhanced when compared to single-layer graphene, we find that the supercurrent is suppressed. The spatial distribution of the supercurrent shows that for Bernal stacking the current flows only in the topmost layers while for rhombohedral stacking the current flows throughout the whole structure.

Abstract: We present a detailed theoretical study of the electronic and transport properties of monolayer black phosphorus (BP). This study is motivated by recent experimental activities in investigating n-type few-layer BP systems. The electron density of states, the screening length, and the low-temperature electron mobility are calculated for monolayer BP (MLBP). In particular, the electron transport mobilities along the armchair and zigzag directions are examined on the basis of the momentum-balance equation derived from a semiclassical Boltzmann equation. The anisotropic electron mobilities in MLBP along different directions are demonstrated where the electron-impurity scattering is considered. Furthermore, we compare the results obtained from two electronic band structures of MLBP and find that the simplified model can describe quite rightly the electronic and transport properties of MLBP. This study is relevant to the application of few-layer BP based electronic systems as advanced electronic devices.

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: A twin boundary in one of the layers of a twisted van der Waals heterostructure separates regions with near opposite inter-layer twist angles. In a MoS₂/WSe₂bilayer, the regions with<inline-formula><tex-math><?CDATA $Rh^h$?></tex-math><math overflow=“scroll”><msubsup><mi>R</mi><mi>h</mi><mi>h</mi></msubsup></math><inline-graphic href=“tdmac529dieqn1.gif” type=“simple” /></inline-formula>and<inline-formula><tex-math><?CDATA $Rh^X$?></tex-math><math overflow=“scroll”><msubsup><mi>R</mi><mi>h</mi><mi>X</mi></msubsup></math><inline-graphic href=“tdmac529dieqn2.gif” type=“simple” /></inline-formula>stacking registry that defined the sub-lattices of the moiré honeycomb pattern would be mirror-reflected across such a twist boundary. In that case, we demonstrate that topologically protected chiral moiré exciton states are confined at the twist boundary. These are one-dimensional and uni-directional excitons with opposite velocities for excitons composed by electronic states with opposite valley/spin character, enabling intrinsic, guided, and far reaching valley-polarized exciton currents.

Abstract: The interaction of monolayer graphene with specific substrates may break its sublattice symmetry and results in unidirectional chiral states with opposite group velocities in the different Dirac cones (Zarenia et al 2012 Phys. Rev. B 86 085451). Taking advantage of this feature, we propose a valley filter based on a transversal mass kink for low energy electrons in graphene, which is obtained by assuming a defect region in the substrate that provides a change in the sign of the substrate-induced mass and thus creates a non-biased channel, perpendicular to the kink, for electron motion. By solving the time-dependent Schrodinger equation for the tight-binding Hamiltonian, we investigate the time evolution of a Gaussian wave packet propagating through such a system and obtain the transport properties of this graphene-based substrate-induced quantum point contact. Our results demonstrate that efficient valley filtering can be obtained, provided: (i) the electron energy is sufficiently low, i.e. with electrons belonging mostly to the lowest sub-band of the channel, and (ii) the channel length (width) is sufficiently long (narrow). Moreover, even though the transmission probabilities for each valley are significantly affected by impurities and defects in the channel region, the valley polarization in this system is shown to be robust against their presence.