Gate-controlled suppression of light-driven proton transport through graphene electrodes

Recent experiments demonstrated that proton transport through graphene electrodes can be accelerated by over an order of magnitude with low intensity illumination. Here we show that this photo-effect can be suppressed for a tuneable fraction of the infra-red spectrum by applying a voltage bias. Using photocurrent measurements and Raman spectroscopy, we show that such fraction can be selected by tuning the Fermi energy of electrons in graphene with a bias, a phenomenon controlled by Pauli blocking of photo-excited electrons. These findings demonstrate a dependence between graphene's electronic and proton transport properties and provide fundamental insights into molecularly thin electrode-electrolyte interfaces and their interaction with light.

Proton transport through nanoscale corrugations in two-dimensional crystals

Defect-free graphene is impermeable to all atoms1-5 and ions6,7 under ambient conditions. Experiments that can resolve gas flows of a few atoms per hour through micrometre-sized membranes found that monocrystalline graphene is completely impermeable to helium, the smallest atom2,5. Such membranes were also shown to be impermeable to all ions, including the smallest one, lithium6,7. By contrast, graphene was reported to be highly permeable to protons, nuclei of hydrogen atoms8,9. There is no consensus, however, either on the mechanism behind the unexpectedly high proton permeability10-14 or even on whether it requires defects in graphene's crystal lattice6,8,15-17. Here, using high-resolution scanning electrochemical cell microscopy, we show that, although proton permeation through mechanically exfoliated monolayers of graphene and hexagonal boron nitride cannot be attributed to any structural defects, nanoscale non-flatness of two-dimensional membranes greatly facilitates proton transport. The spatial distribution of proton currents visualized by scanning electrochemical cell microscopy reveals marked inhomogeneities that are strongly correlated with nanoscale wrinkles and other features where strain is accumulated. Our results highlight nanoscale morphology as an important parameter enabling proton transport through two-dimensional crystals, mostly considered and modelled as flat, and indicate that strain and curvature can be used as additional degrees of freedom to control the proton permeability of two-dimensional materials.

Clogging and Unclogging of Hydrocarbon-Contaminated Nanochannels

The recent advantages of the fabrication of artificial nanochannels enabled new research on the molecular transport, permeance, and selectivity of various gases and molecules. However, the physisorption/chemisorption of the unwanted molecules (usually hydrocarbons) inside nanochannels results in the alteration of the functionality of the nanochannels. We investigated contamination due to hydrocarbon molecules, nanochannels made of graphene, hexagonal boron nitride, BC₂N, and molybdenum disulfide using molecular dynamics simulations. We found that for a certain size of nanochannel (i.e., h = 0.7 nm), as a result of the anomalous hydrophilic nature of nanochannels made of graphene, the hydrocarbons are fully adsorbed in the nanochannel, giving rise to full uptake. An increasing temperature plays an important role in unclogging, while pressure does not have a significant role. The results of our pioneering work contribute to a better understanding and highlight the important factors in alleviating the contamination and unclogging of nanochannels, which are in good agreement with the results of recent experiments.

Vacancy clustering effect on the electronic and transport properties of bilayer graphene nanoribbons

Experimental realizations of two-dimensional materials are hardly free of structural defects such as e.g. vacancies, which, in turn, modify drastically its pristine physical defect-free properties. In this work, we explore effects due to point defect clustering on the electronic and transport properties of bilayer graphene nanoribbons, for AA and AB stacking and zigzag and armchair boundaries, by means of the tight-binding approach and scattering matrix formalism. Evident vacancy concentration signatures exhibiting a maximum amplitude and an universality regardless of the system size, stacking and boundary types, in the density of states around the zero-energy level are observed. Our results are explained via the coalescence analysis of the strong sizeable vacancy clustering effect in the system and the breaking of the inversion symmetry at high vacancy densities, demonstrating a similar density of states for two equivalent degrees of concentration disorder, below and above the maximum value.

Wien effect in interfacial water dissociation through proton-permeable graphene electrodes

Strong electric fields can accelerate molecular dissociation reactions. The phenomenon known as the Wien effect was previously observed using high-voltage electrolysis cells that produced fields of about 10⁷ V m^-1, sufficient to accelerate the dissociation of weakly bound molecules (e.g., organics and weak electrolytes). The observation of the Wien effect for the common case of water dissociation (H₂O [Formula: see text] H⁺ + OH^-) has remained elusive. Here we study the dissociation of interfacial water adjacent to proton-permeable graphene electrodes and observe strong acceleration of the reaction in fields reaching above 10⁸ V m^-1. The use of graphene electrodes allows measuring the proton currents arising exclusively from the dissociation of interfacial water, while the electric field driving the reaction is monitored through the carrier density induced in graphene by the same field. The observed exponential increase in proton currents is in quantitative agreement with Onsager's theory. Our results also demonstrate that graphene electrodes can be valuable for the investigation of various interfacial phenomena involving proton transport.

Interface-dependent phononic and optical properties of GeO/MoSO heterostructures

The interface-dependent electronic, vibrational, piezoelectric, and optical properties of van der Waals heterobilayers, formed by buckled GeO (b-GeO) and Janus MoSO structures, are investigated by means of first-principles calculations. The electronic band dispersions show that O/Ge and S/O interface formations result in a type-II band alignment with direct and indirect band gaps, respectively. In contrast, O/O and S/Ge interfaces give rise to the formation of a type-I band alignment with an indirect band gap. By considering the Bethe-Salpeter equation (BSE) on top of G₀W₀ approximation, it is shown that different interfaces can be distinguished from each other by means of the optical absorption spectra as a consequence of the band alignments. Additionally, the low- and high-frequency regimes of the Raman spectra are also different for each interface type. The alignment of the individual dipoles, which is interface-dependent, either weakens or strengthens the net dipole of the heterobilayers and results in tunable piezoelectric coefficients. The results indicate that the possible heterobilayers of b-GeO/MoSO asymmetric structures possess various electronic, optical, and piezoelectric properties arising from the different interface formations and can be distinguished by means of various spectroscopic techniques.

Electronic Properties of Oxidized Graphene: Effects of Strain and an Electric Field on Flat Bands and the Energy Gap

A multiscale modeling and simulation approach, including first-principles calculations, ab initio molecular dynamics simulations, and a tight binding approach, is employed to study band flattening of the electronic band structure of oxidized monolayer graphene. The width of flat bands can be tuned by strain, the external electric field, and the density of functional groups and their distribution. A transition to a conducting state is found for monolayer graphene with impurities when it is subjected to an electric field of ∼1.0 V/Å. Several parallel impurity-induced flat bands appear in the low-energy spectrum of monolayer graphene when the number of epoxy groups is changed. The width of the flat band decreases with an increase in tensile strain but is independent of the electric field strength. Here an alternative and easy route for obtaining band flattening in thermodynamically stable functionalized monolayer graphene is introduced. Our work discloses a new avenue for research on band flattening in monolayer graphene.

Stability of adsorption of Mg and Na on sulfur-functionalized MXenes

Two-dimensional materials composed of transition metal carbides and nitrides (MXenes) are poised to revolutionize energy conversion and storage. In this work, we used density functional theory (DFT) to investigate the adsorption of Mg and Na adatoms on five M₂CS₂ monolayers (where M = Mo, Nb, Ti, V, and Zr) for battery applications. We assessed the stability of the adatom (i.e. Na and Mg)-monolayer systems by calculating adsorption and formation energies, as well as voltages as a function of surface coverage. For instance, we found that Mo₂CS₂ cannot support a full layer of Na nor even a single Mg atom. Na and Mg exhibit the strongest binding on Zr₂CS₂, followed by Ti₂CS₂, Nb₂CS₂ and V₂CS₂. Using the nudged elastic band method (NEB), we computed promising diffusion barriers for both dilute and nearly full ion surface coverage cases. In the dilute ion adsorption case, a single Mg and Na atom on Ti₂CS₂ experience ∼0.47 eV and ∼0.10 eV diffusion barriers between the lowest energy sites, respectively. For a nearly full surface coverage, a Na ion moving on Ti₂CS₂ experiences a ∼0.33 eV energy barrier, implying a concentration-dependent diffusion barrier. Our molecular dynamics results indicate that the three (one) layers (layer) of the Mg (Na) ion on both surfaces of Ti₂CS₂ remain stable at T = 300 K. While, according to voltage calculations, Zr₂CS₂ can store Na up to three atomic layers, our MD simulations predict that the outermost layers detach from the Zr₂CS₂ monolayer due to the weak interaction between Na ions and the monolayer. This suggests that MD simulations are essential to confirm the stability of an ion-electrode system - an insight that is mostly absent in previous studies.

Photoluminescence and electronic transition behaviors of single-stranded DNA

Due to the potential application of DNA for biophysics and optoelectronics, the electronic energy states and transitions of this genetic material have attracted a great deal of attention recently. However, the fluorescence and corresponding physical process of DNA under optical excitation with photon energies below ultraviolet are still not fully clear. In this work, we experimentally investigate the photoluminescence (PL) properties of single-stranded DNA (ssDNA) samples under near-ultraviolet (NUV) and visible excitations (270∼440 nm). Based on the dependence of the PL peak wavelength (λ_{em}) upon the excitation wavelength (λ_{ex}), the PL behaviors of ssDNA can be approximately classified into two categories. In the relatively short excitation wavelength regime, λ_{em} is nearly constant due to exciton-like transitions associated with delocalized excitonic states and excimer states. In the relatively long excitation wavelength range, a linear relation of λ_{em}=Aλ_{ex}+B with A>0 or A<0 can be observed, which comes from electronic transitions related to coupled vibrational-electronic levels. Moreover, the transition channels in different excitation wavelength regimes and the effects of strand length and base type can be analyzed on the basis of these results. These important findings not only can give a general description of the electronic energy states and transitional behaviors of ssDNA samples under NUV and visible excitations, but also can be the basis for the application of DNA in nanoelectronics and optoelectronics.

Substrate dependent terahertz magneto-optical properties of monolayer WS2

Terahertz (THz) magneto-optical (MO) properties of monolayer (ML) tungsten disulfide (WS₂), placed on different substrates and subjected to external magnetic fields, are studied using THz time-domain spectroscopy (TDS). We find that the THz MO conductivity exhibits a nearly linear response in a weak magnetic field, while a distinctly nonlinear/oscillating behavior is found in strong magnetic fields owing to strong substrate-induced random impurity scattering and interactions. The THz MO response of ML WS₂ depends sensitively on the choice of the substrates, which we trace back to electronic localization and the impact of the substrates on the Landau level (LL) spectrum. Our results provide an in-depth understanding of the THz MO properties of ML WS₂/substrate systems, especially the effect of substrates, which can be utilized to realize atomically thin THz MO nano-devices.

Zitterbewegung of Moiré Excitons in Twisted MoS_{2}/WSe_{2} Heterobilayers

The moiré pattern observed in stacked noncommensurate crystal lattices, such as heterobilayers of transition metal dichalcogenides, produces a periodic modulation of their band gap. Excitons subjected to this potential landscape exhibit a band structure that gives rise to a quasiparticle dubbed the moiré exciton. In the case of MoS_{2}/WSe_{2} heterobilayers, the moiré trapping potential has honeycomb symmetry and, consequently, the moiré exciton band structure is the same as that of a Dirac-Weyl fermion, whose mass can be further tuned down to zero with a perpendicularly applied field. Here we show that, analogously to other Dirac-like particles, the moiré exciton exhibits a trembling motion, also known as Zitterbewegung, whose long timescales are compatible with current experimental techniques for exciton dynamics. This promotes the study of the dynamics of moiré excitons in van der Waals heterostructures as an advantageous solid-state platform to probe Zitterbewegung, broadly tunable by gating and interlayer twist angle.

Strong anisotropic optical properties of 8-Pmmn borophene: a many-body perturbation study

Using first-principles many-body perturbation theory, we investigate the optical properties of 8-Pmmn borophene at two levels of approximations; the GW method considering only the electron-electron interaction and the GW in combination with the Bethe-Salpeter equation including electron-hole coupling. The band structure exhibits anisotropic Dirac cones with semimetallic character. The optical absorption spectra are obtained for different light polarizations and we predict strong optical absorbance anisotropy. The absorption peaks undergo a global redshift when the electron-hole interaction is taken into account due to the formation of bound excitons which have an anisotropic excitonic wave function.

Effect of zitterbewegung on the propagation of wave packets in ABC-stacked multilayer graphene: an analytical and computational approach

The time evolution of a low-energy two-dimensional Gaussian wave packet in ABC-stacked n-layer graphene (ABC-NLG) is investigated. Expectation values of the position (x, y) of center-of-mass and the total probability densities of the wave packet are calculated analytically using the Green's function method. These results are confirmed using an alternative numerical method based on the split-operator technique within the Dirac approach for ABC-NLG, which additionally allows to include external fields and potentials. The main features of the zitterbewegung (trembling motion) of wave packets in graphene are demonstrated and are found to depend not only on the wave packet width and initial pseudospin polarization, but also on the number of layers. Moreover, the analytical and numerical methods proposed here allow to investigate wave packet dynamics in graphene systems with an arbitrary number of layers and arbitrary potential landscapes.

Transition-metal adatoms on 2D-GaAs: a route to chiral magnetic 2D materials by design

Using relativistic density-functional calculations, we examine the magneto-crystalline anisotropy and exchange properties of transition-metal atoms adsorbed on 2D-GaAs. We show that single Mn and Mo atom (Co and Os) strongly bind on 2D-GaAs, and induce local out-of-plane (in-plane) magnetic anisotropy. When a pair of TM atoms is adsorbed on 2D-GaAs in a close range from each other, magnetisation properties change (become tunable) with respect to concentrations and ordering of the adatoms. In all cases, we reveal presence of strong Dzyaloshinskii-Moriya interaction. These results indicate novel pathways towards two-dimensional chiral magnetic materials by design, tailored for desired applications in magneto-electronics.

Out-of-plane permittivity of confined water

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Å, 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Å. Our results are in excellent agreement with a recent experiment [L. Fumagalli et al., Science 360, 1339 (2018)SCIEAS0036-807510.1126/science.aat4191].

Strain and electric field tuning of semi-metallic character WCrCO2MXenes with dual narrow band gap

Motivated by the recent successful synthesis of double-M carbides, we investigate structural and electronic properties of WCrC and WCrCO2monolayers and the effects of biaxial and out-of-plane strain and electric field using density functional theory. WCrC and WCrCO2monolayers are found to be dynamically stable. WCrC is metallic and WCrCO2display semi-metallic character with narrow band gap, which can be controlled by strain engineering and electric field. WCrCO2monolayer exhibits a dual band gap which is preserved in the presence of an electric field. The band gap of WCrCO2monolayer increases under uniaxial strain while it becomes metallic under tensile strain, resulting in an exotic 2D double semi-metallic behavior. Our results demonstrate that WCrCO2is a new platform for the study of novel physical properties in two-dimensional Dirac materials and which may provide new opportunities to realize high-speed low-dissipation devices.

Stable single-layers of calcium halides (CaX2, X = F, Cl, Br, I)

By means of density functional theory based first-principles calculations, the structural, vibrational, and electronic properties of 1H- and 1T-phases of single-layer CaX₂ (X = F, Cl, Br, or I) structures are investigated. Our results reveal that both the 1H- and 1T-phases are dynamically stable in terms of their phonon band dispersions with the latter being the energetically favorable phase for all single-layers. In both phases of single-layer CaX₂ structures, significant phonon softening occurs as the atomic radius increases. In addition, each structural phase exhibits distinctive Raman active modes that enable one to characterize either the phase or the structure via Raman spectroscopy. The electronic band dispersions of single-layer CaX₂ structures reveal that all structures are indirect bandgap insulators with a decrease in bandgaps from fluorite to iodide crystals. Furthermore, the calculated linear elastic constants, in-plane stiffness, and Poisson ratio indicate the ultra-soft nature of CaX₂ single-layers, which is quite important for their nanoelastic applications. Overall, our study reveals that with their dynamically stable 1T- and 1H-phases, single-layers of CaX₂ crystals can be alternative ultra-thin insulators.

Magnetic field induced vortices in graphene quantum dots

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.

Two-dimensional hydrogenated buckled gallium arsenide: an ab initio study

First-principles calculations have been carried out to investigate the stability, structural and electronic properties of two-dimensional (2D) hydrogenated GaAs with three possible geometries: chair, zigzag-line and boat configurations. The effect of van der Waals interactions on 2D H-GaAs systems has also been studied. These configurations were found to be energetic and dynamic stable, as well as having a semiconducting character. Although 2D GaAs adsorbed with H tends to form a zigzag-line configuration, the energy differences between chair, zigzag-line and boat are very small which implies the metastability of the system. Chair and boat configurations display a [Formula: see text]-[Formula: see text] direct bandgap nature, while pristine 2D-GaAs and zigzag-line are indirect semiconductors. The bandgap sizes of all configurations are also hydrogen dependent, and wider than that of pristine 2D-GaAs with both PBE and HSE functionals. Even though DFT-vdW interactions increase the adsorption energies and reduce the equilibrium distances of H-GaAs systems, it presents, qualitatively, the same physical results on the stability and electronic properties of our studied systems with PBE functional. According to our results, 2D buckled gallium arsenide is a good candidate to be synthesized by hydrogen surface passivation as its group III-V partners 2D buckled gallium nitride and boron nitride. The hydrogenation of 2D-GaAs tunes the bandgap of pristine 2D-GaAs, which makes it a potential candidate for optoelectronic applications in the blue and violet ranges of the visible electromagnetic spectrum.

Composite super-moiré lattices in double-aligned graphene heterostructures

When two-dimensional (2D) atomic crystals are brought into close proximity to form a van der Waals heterostructure, neighbouring crystals may influence each other's properties. Of particular interest is when the two crystals closely match and a moiré pattern forms, resulting in modified electronic and excitonic spectra, crystal reconstruction, and more. Thus, moiré patterns are a viable tool for controlling the properties of 2D materials. However, the difference in periodicity of the two crystals limits the reconstruction and, thus, is a barrier to the low-energy regime. Here, we present a route to spectrum reconstruction at all energies. By using graphene which is aligned to two hexagonal boron nitride layers, one can make electrons scatter in the differential moiré pattern which results in spectral changes at arbitrarily low energies. Further, we demonstrate that the strength of this potential relies crucially on the atomic reconstruction of graphene within the differential moiré super cell.

Transport characteristics of multi-terminal pristine and defective phosphorene systems

Atomic vacancies and nanopores act as local scattering centers and modify the transport properties of charge carriers in phosphorene nanoribbons (PNRs). We investigate the influence of such atomic defects on the electronic transport of multi-terminal PNR. We use the non-equilibrium Green's function approach within the tight-binding framework to calculate the transmission coefficient and the conductance. Terminals induce band mixing resulting in oscillations in the conductance. In the presence of atomic vacancies and nanopores the conductance between non-axial terminals exhibit constructive scattering, which is in contrast to mono-axial two-terminal systems where the conductance exhibits destructive scattering. This can be understood from the spatial local density of states of the transport modes in the system. Our results provide fundamental insights into the electronic transport in PNR-based multi-terminal systems and into the ability of atomic defects and nanopores through tuning the transport properties.

Atomically thin micas as proton-conducting membranes

Monolayers of graphene and hexagonal boron nitride (hBN) are highly permeable to thermal protons^1,2. For thicker two-dimensional (2D) materials, proton conductivity diminishes exponentially, so that, for example, monolayer MoS₂ that is just three atoms thick is completely impermeable to protons¹. This seemed to suggest that only one-atom-thick crystals could be used as proton-conducting membranes. Here, we show that few-layer micas that are rather thick on the atomic scale become excellent proton conductors if native cations are ion-exchanged for protons. Their areal conductivity exceeds that of graphene and hBN by one to two orders of magnitude. Importantly, ion-exchanged 2D micas exhibit this high conductivity inside the infamous gap for proton-conducting materials³, which extends from ∼100 °C to 500 °C. Areal conductivity of proton-exchanged monolayer micas can reach above 100 S cm^-2 at 500 °C, well above the current requirements for the industry roadmap⁴. We attribute the fast proton permeation to ~5-Å-wide tubular channels that perforate micas' crystal structure, which, after ion exchange, contain only hydroxyl groups inside. Our work indicates that there could be other 2D crystals⁵ with similar nanometre-scale channels, which could help close the materials gap in proton-conducting applications.

Tunable 2D-gallium arsenide and graphene bandgaps in a graphene/GaAs heterostructure: an ab initio study

The bandgap behavior of 2D-GaAs and graphene have been investigated with van der Waals heterostructured into a yet unexplored graphene/GaAs bilayer, under both uniaxial stress along c axis and different planar strain distributions. The 2D-GaAs bandgap nature changes from [Formula: see text]-K indirect in isolated monolayer to [Formula: see text]-[Formula: see text] direct in graphene/GaAs bilayer. In the latter, graphene exhibits a bandgap of 5 meV. The uniaxial stress strongly affects the graphene electronic bandgap, while symmetric in-plane strain does not open the bandgap in graphene. Nevertheless, it induces remarkable changes on the GaAs bandgap-width around the Fermi level. However, when applying asymmetric in-plane strain to graphene/GaAs, the graphene sublattice symmetry is broken, and the graphene bandgap is open at the Fermi level to a maximum width of 814 meV. This value is much higher than that reported for just graphene under asymmetric strain. The [Formula: see text]-[Formula: see text] direct bandgap of GaAs remains unchanged in graphene/GaAs under different types of applied strain. The analyses of phonon dispersion and the elastic constants yield the dynamical and mechanical stability of the graphene/GaAs system, respectively. The calculated mechanical properties for bilayer heterostructure are better than those of their constituent monolayers. This finding, together with the tunable graphene bandgap not only by the strength but also by the direction of the strain, enhance the potential for strain engineering of ultrathin group-III-V electronic devices hybridized by graphene.

Vibrational properties of germanane and fluorinated germanene in the chair, boat, and zigzag-line configurations

The electronic and vibrational properties of germanane and fluorinated germanene are studied within density functional theory (DFT) and density functional perturbation theory frameworks. Different structural configurations of germanane and fluorinated germanene are investigated. The energy difference between the different configurations are consistently smaller than the energy of thermal fluctuations for all the analyzed DFT functionals LDA, GGA, and hybrid functionals, which implies that, in principle, it is possible to find these different configurations in different regions of the sample as minority phases or local defects. We calculate the Raman and infrared spectra for these configurations by using ab initio calculations and compare it with available experimental spectra for germanane. Our results show the presence of minority phases compatible with the configurations analyzed in this work. As these low energy configurations are metastable the present work shows that the synthesis of these energy competing phases is feasible by selectively changing the synthesis conditions, which is an opportunity to expand in this way the availability of new two-dimensional compounds.

Single-layer structures of a100- and b010-Gallenene: a tight-binding approach

Using the simplified linear combination of atomic orbitals (LCAO) method in combination with ab initio calculations, we construct a tight-binding (TB) model for two different crystal structures of monolayer gallium: a₁₀₀- and b₀₁₀-Gallenene.

Ionized water confined in graphene nanochannels

When confined between graphene layers, water in the presence of additional hydronium and hydroxide ions exhibits distinct properties such as ion layering structure determined by the channel size, disruption of the ion solvation shell, and slower ion recombination rate as compared to bulk water.

Confined states in graphene quantum blisters

Bilayer graphene samples may exhibit regions where the two layers are locally delaminated forming a so-called quantum blister in the graphene sheet. Electron and hole states can be confined in this graphene quantum blisters (GQB) by applying a global electrostatic bias. We scrutinize the electronic properties of these confined states under the variation of interlayer bias, coupling, and blister's size. The spectra display strong anti-crossings due to the coupling of the confined states on upper and lower layers inside the blister. These spectra are layer localized where the respective confined states reside on either layer or equally distributed. For finite angular momentum, this layer localization can be at the edge of the blister and corresponds to degenerate modes of opposite momenta. Furthermore, the energy levels in GQB exhibit electron-hole symmetry that is sensitive to the electrostatic bias. Finally, we demonstrate that confinement in GQB persists even in the presence of a variation in the inter-layer coupling.

Electrical generation of terahertz blackbody radiation from graphene

Recent experimental work on the application of graphene for novel illumination motivated us to present a theoretical study of the blackbody radiation emission from a freely suspended graphene driven by a dc electric field. Strong terahertz (THz) emission, with intensity up to mW/cm², can be generated with increasing electric field strength due to the heating of electrons in graphene. We show that the intensity of the THz emission generated electrically from graphene depends rather sensitively on the lattice temperature in relatively weak electric fields, whereas it is less sensitive to the lattice temperature in relative strong electric fields. Our study highlights the practical application of graphene as intense THz source where the radiation is generated electrically.

Multiband Mechanism for the Sign Reversal of Coulomb Drag Observed in Double Bilayer Graphene Heterostructures

Coupled 2D sheets of electrons and holes are predicted to support novel quantum phases. Two experiments of Coulomb drag in electron-hole (e-h) double bilayer graphene (DBLG) have reported an unexplained and puzzling sign reversal of the drag signal. However, we show that this effect is due to the multiband character of DBLG. Our multiband Fermi liquid theory produces excellent agreement and captures the key features of the experimental drag resistance for all temperatures. This demonstrates the importance of multiband effects in DBLG: they have a strong effect not only on superfluidity, but also on the drag.

Phase transition and field effect topological quantum transistor made of monolayer MoS2

We study topological phase transitions and topological quantum field effect transistor in monolayer molybdenum disulfide (MoS₂) using a two-band Hamiltonian model. Without considering the quadratic (q ²) diagonal term in the Hamiltonian, we show that the phase diagram includes quantum anomalous Hall effect, quantum spin Hall effect, and spin quantum anomalous Hall effect regions such that the topological Kirchhoff law is satisfied in the plane. By considering the q ² diagonal term and including one valley, it is shown that MoS₂ has a non-trivial topology, and the valley Chern number is non-zero for each spin. We show that the wave function is (is not) localized at the edges when the q ² diagonal term is added (deleted) to (from) the spin-valley Dirac mass equation. We calculate the quantum conductance of zigzag MoS₂ nanoribbons by using the nonequilibrium Green function method and show how this device works as a field effect topological quantum transistor.

Transport of hydrogen isotopes through interlayer spacing in van der Waals crystals

Atoms start behaving as waves rather than classical particles if confined in spaces commensurate with their de Broglie wavelength. At room temperature this length is only about one ångström even for the lightest atom, hydrogen. This restricts quantum-confinement phenomena for atomic species to the realm of very low temperatures^1-5. Here, we show that van der Waals gaps between atomic planes of layered crystals provide ångström-size channels that make quantum confinement of protons apparent even at room temperature. Our transport measurements show that thermal protons experience a notably higher barrier than deuterons when entering van der Waals gaps in hexagonal boron nitride and molybdenum disulfide. This is attributed to the difference in the de Broglie wavelengths of the isotopes. Once inside the crystals, transport of both isotopes can be described by classical diffusion, albeit with unexpectedly fast rates comparable to that of protons in water. The demonstrated ångström-size channels can be exploited for further studies of atomistic quantum confinement and, if the technology can be scaled up, for sieving hydrogen isotopes.

Multicomponent Electron-Hole Superfluidity and the BCS-BEC Crossover in Double Bilayer Graphene

Superfluidity in coupled electron-hole sheets of bilayer graphene is predicted here to be multicomponent because of the conduction and valence bands. We investigate the superfluid crossover properties as functions of the tunable carrier densities and the tunable energy band gap E_{g}. For small band gaps there is a significant boost in the two superfluid gaps, but the interaction-driven excitations from the valence to the conduction band can weaken the superfluidity, even blocking the system from entering the Bose-Einstein condensate (BEC) regime at low densities. At a given larger density, a band gap E_{g}∼80-120  meV can carry the system into the strong-pairing multiband BCS-BEC crossover regime, the optimal range for realization of high-T_{c} superfluidity.

Free surfaces recast superconductivity in few-monolayer MgB2: Combined first-principles and ARPES demonstration

Two-dimensional materials are known to harbour properties very different from those of their bulk counterparts. Recent years have seen the rise of atomically thin superconductors, with a caveat that superconductivity is strongly depleted unless enhanced by specific substrates, intercalants or adatoms. Surprisingly, the role in superconductivity of electronic states originating from simple free surfaces of two-dimensional materials has remained elusive to date. Here, based on first-principles calculations, anisotropic Eliashberg theory, and angle-resolved photoemission spectroscopy (ARPES), we show that surface states in few-monolayer MgB2 make a major contribution to the superconducting gap spectrum and density of states, clearly distinct from the widely known, bulk-like σ- and π-gaps. As a proof of principle, we predict and measure the gap opening on the magnesium-based surface band up to a critical temperature as high as ~30 K for merely six monolayers thick MgB2. These findings establish free surfaces as an unavoidable ingredient in understanding and further tailoring of superconductivity in atomically thin materials.

Quantum transport across van der Waals domain walls in bilayer graphene

Bilayer graphene can exhibit deformations such that the two graphene sheets are locally detached from each other resulting in a structure consisting of domains with different van der Waals inter-layer coupling. Here we investigate how the presence of these domains affects the transport properties of bilayer graphene. We derive analytical expressions for the transmission probability, and the corresponding conductance, across walls separating different inter-layer coupling domains. We find that the transmission can exhibit a valley-dependent layer asymmetry and that the domain walls have a considerable effect on the chiral tunnelling properties of the charge carriers. We show that transport measurements allow one to obtain the strength with which the two layers are coupled. We perform numerical calculations for systems with two domain walls and find that the availability of multiple transport channels in bilayer graphene significantly modifies the conductance dependence on inter-layer potential asymmetry.

Microfluidic manipulation of magnetic flux domains in type-I superconductors: droplet formation, fusion and fission

The magnetic flux domains in the intermediate state of type-I superconductors are known to resemble fluid droplets, and their dynamics in applied electric current is often cartooned as a "dripping faucet". Here we show, using the time-depended Ginzburg-Landau simulations, that microfluidic principles hold also for the determination of the size of the magnetic flux-droplet as a function of the applied current, as well as for the merger or splitting of those droplets in the presence of the nanoengineered obstacles for droplet motion. Differently from fluids, the flux-droplets in superconductors are quantized and dissipative objects, and their pinning/depinning, nucleation, and splitting occur in a discretized form, all traceable in the voltage measured across the sample. At larger applied currents, we demonstrate how obstacles can cause branching of laminar flux streams or their transformation into mobile droplets, as readily observed in experiments.

Structure and reentrant percolation in an inverse patchy colloidal system

Two-dimensional systems of inverse patchy colloids modeled as disks with a central charge and having their surface decorated with oppositely pointlike charged patches are investigated using molecular dynamics simulations. The self-assembly of the patchy colloids leads to diverse ground state configurations ranging from crystalline arrangements of monomers to linear clusters, ramified linear clusters and to percolated configurations. Two structural phase diagrams are constructed: (1) as a function of the net charge and area fraction, and (2) as a function of the net charge and the range of the pair interaction potential. An interesting reentrant percolation transition is obtained as a function of the net charge of the colloids. We identify distinct mechanisms that lead to the percolation transition.

Valley filtering in graphene due to substrate-induced mass potential

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 Schrödinger 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.

Scanning gate microscopy of magnetic focusing in graphene devices: quantum versus classical simulation

We compare classical versus quantum electron transport in recently investigated magnetic focusing devices (Bhandari et al 2016 Nano Lett. 16 1690) exposed to the perturbing potential of a scanning gate microscope (SGM). Using the Landauer-Büttiker formalism for a multi-terminal device, we calculate resistance maps that are obtained as the SGM tip is scanned over the sample. There are three unique regimes in which the scanning tip can operate (focusing, repelling, and mixed regime) which are investigated. Tip interacts mostly with electrons with cyclotron trajectories passing directly underneath it, leaving a trail of modified current density behind it. Other (indirect) trajectories become relevant when the tip is placed near the edges of the sample, and current is scattered between the tip and the edge. We point out that, in contrast to SGM experiments on gapped semiconductors, the STM tip can induce a pn junction in graphene, which improves contrast and resolution in SGM. We also discuss possible explanations for spatial asymmetry of experimentally measured resistance maps, and connect it with specific configurations of the measuring probes.

Ab initio study of hydrogenic effective mass impurities in Si nanowires

The effect of B and P dopants on the band structure of Si nanowires is studied using electronic structure calculations based on density functional theory. At low concentrations a dispersionless band is formed, clearly distinguishable from the valence and conduction bands. Although this band is evidently induced by the dopant impurity, it turns out to have purely Si character. These results can be rigorously analyzed in the framework of effective mass theory. In the process we resolve some common misconceptions about the physics of hydrogenic shallow impurities, which can be more clearly elucidated in the case of nanowires than would be possible for bulk Si. We also show the importance of correctly describing the effect of dielectric confinement, which is not included in traditional electronic structure calculations, by comparing the obtained results with those of G₀W₀ calculations.

Electric- and magnetic-field dependence of the electronic and optical properties of phosphorene quantum dots

Recently, black phosphorus quantum dots were fabricated experimentally. Motivated by these experiments, we theoretically investigate the electronic and optical properties of rectangular phosphorene quantum dots (RPQDs) in the presence of an in-plane electric field and a perpendicular magnetic field. The energy spectra and wave functions of RPQDs are obtained numerically using the tight-binding approach. We find edge states within the band gap of the RPQD which are well separated from the bulk states. In an undoped RPQD and for in-plane polarized light, due to the presence of well-defined edge states, we find three types of optical transitions which are between the bulk states, between the edge and bulk states, and between the edge states. The electric and magnetic fields influence the bulk-to-bulk, edge-to-bulk, and edge-to-edge transitions differently due to the different responses of bulk and edge states to these fields.

Strained graphene Hall bar

The effects of strain, induced by a Gaussian bump, on the magnetic field dependent transport properties of a graphene Hall bar are investigated. The numerical simulations are performed using both classical and quantum mechanical transport theory and we found that both approaches exhibit similar characteristic features. The effects of the Gaussian bump are manifested by a decrease of the bend resistance, R _B, around zero-magnetic field and the occurrence of side-peaks in R _B. These features are explained as a consequence of bump-assisted scattering of electrons towards different terminals of the Hall bar. Using these features we are able to give an estimate of the size of the bump. Additional oscillations in R _B are found in the quantum description that are due to the population/depopulation of Landau levels. The bump has a minor influence on the Hall resistance even for very high values of the pseudo-magnetic field. When the bump is placed outside the center of the Hall bar valley polarized electrons can be collected in the leads.

The work function of few-layer graphene

A theoretical and experimental study of the work function of few-layer graphene is reported. The influence of the number of layers on the work function is investigated in the presence of a substrate, a molecular dipole layer, and combinations of the two. The work function of few-layer graphene is almost independent of the number of layers with only a difference between monolayer and multilayer graphene of about 60 meV. In the presence of a charge-donating substrate the charge distribution is found to decay exponentially away from the substrate and this is directly reflected in the work function of few-layer graphene. A dipole layer changes the work function only when placed in between the substrate and few-layer graphene through a change of the charge transfer between the two.

Controlled growth mechanism of poly (3-hexylthiophene) nanowires

Synthesis of 1D-polymer nanowires by a self-assembly method using marginal solvents is an attractive technique. While the formation mechanism is poorly understood, this method is essential in order to control the growth of nanowires. Here we visualized the time-dependent assembly of poly (3-hexyl-thiophene-2,5-diyl) (P3HT) nanowires by atomic force microscopy and scanning tunneling microscopy. The assembly of P3HT nanowires was carried out at room temperature by mixing cyclohexanone (CHN), as a poor solvent, with polymer solution in 1,2-dichlorobenzene (DCB). Both π-π stacking and planarization, obtained at the mix volume ratio of P3HT (in DCB):CHN (10:7), were considered during the investigation. We find that the length of nanowires was determined by the ordering of polymers in the polymer repetition direction. Additionally, our density functional theory calculations revealed that the presence of DCB and CHN molecules that stabilize the structural distortions due to tail group of polymers was essential for the core-wire formation.

Anomalous Dynamical Behavior of Freestanding Graphene Membranes

We report subnanometer, high-bandwidth measurements of the out-of-plane (vertical) motion of atoms in freestanding graphene using scanning tunneling microscopy. By tracking the vertical position over a long time period, a 1000-fold increase in the ability to measure space-time dynamics of atomically thin membranes is achieved over the current state-of-the-art imaging technologies. We observe that the vertical motion of a graphene membrane exhibits rare long-scale excursions characterized by both anomalous mean-squared displacements and Cauchy-Lorentz power law jump distributions.

Structural changes in a Schiff base molecular assembly initiated by scanning tunneling microscopy tip

We report the controlled self-organization and switching of newly designed Schiff base (E)-4-((4-(phenylethynyl) benzylidene) amino) benzenethiol (EPBB) molecules on a Au (111) surface at room temperature. Scanning tunneling microscopy and spectroscopy (STM/STS) were used to image and analyze the conformational changes of the EPBB molecules. The conformational change of the molecules was induced by using the STM tip while increasing the tunneling current. The switching of a domain or island of molecules was shown to be induced by the STM tip during scanning. Unambiguous fingerprints of the switching mechanism were observed via STM/STS measurements. Surface-enhanced Raman scattering was employed, to control and identify quantitatively the switching mechanism of molecules in a monolayer. Density functional theory calculations were also performed in order to understand the microscopic details of the switching mechanism. These calculations revealed that the molecular switching behavior stemmed from the strong interaction of the EPBB molecules with the STM tip. Our approach to controlling intermolecular mechanics provides a path towards the bottom-up assembly of more sophisticated molecular machines.

Giant magnetic anisotropy in doped single layer molybdenum disulfide and fluorographene

Stable monolayer materials based on existing, well known and stable two-dimensional crystal fluorographene and molybdenum disulfide are predicted to exhibit a huge magnetocrystalline anisotropy when functionalized with adsorbed transition metal atoms at vacant sides. Ab initio calculations within the density-functional theory formalism were performed to investigate the adsorption of the transitional metals in a single S (or F) vacancy of monolayer molybdenum disulfide (or fluorographene). We found strong bonding of the transitional metal atoms to the vacant sites with binding energies ranging from 2.5 to 5.2 eV. Our calculations revealed that these systems with adsorbed metal atoms exhibit a magnetic anisotropy, specifically the structures including Os and Ir show a giant magnetocrystalline anisotropy energy of 31-101 meV. Our results demonstrate the possibility of obtaining stable monolayer materials with huge magnetocrystalline anisotropy based on preexisting, well known and stable two-dimensional crystals: fluorographene and molybdenum disulfide. We believe that the results obtained here are useful not only for deeper understanding of the origin of magnetocrystalline anisotropy but also for the design of monolayer optoelectronic devices with novel functionalities.

The 30-band k ⋅ p theory of valley splitting in silicon thin layers

The valley splitting of the conduction-band states in a thin silicon-on-insulator layer is investigated using the 30-band k ⋅ p theory. The system composed of a few nm thick [Formula: see text] layer embedded within thick SiO2 layers is analyzed. The valley split states are found to cross periodically with increasing quantum well width, and therefore the energy splitting is an oscillatory function of the quantum well width, with period determined by the wave vector K 0 of the conduction band minimum. Because the valley split states are classified by parity, the optical transition between the ground hole state and one of those valley split conduction band states is forbidden. The oscillations in the valley splitting energy decrease with electric field and with smoothing of the composition profile between the well and the barrier by diffusion of oxygen from the SiO2 layers to the Si quantum well. Such a smoothing also leads to a decrease of the interband transition matrix elements. The obtained results are well parametrized by the effective two-valley model, but are found to disagree from previous 30-band calculations. This discrepancy could be traced back to the fact that the basis for the numerical solution of the eigenproblem must be restricted to the first Brillouin zone in order to obtain quantitatively correct results for the valley splitting.

Large gap electron-hole superfluidity and shape resonances in coupled graphene nanoribbons

We predict enhanced electron-hole superfluidity in two coupled electron-hole armchair-edge terminated graphene nanoribbons separated by a thin insulating barrier. In contrast to graphene monolayers, the multiple subbands of the nanoribbons are parabolic at low energy with a gap between the conduction and valence bands, and with lifted valley degeneracy. These properties make screening of the electron-hole interaction much weaker than for coupled electron-hole monolayers, thus boosting the pairing strength and enhancing the superfluid properties. The pairing strength is further boosted by the quasi one-dimensional quantum confinement of the carriers, as well as by the large density of states near the bottom of each subband. The latter magnifies superfluid shape resonances caused by the quantum confinement. Several superfluid partial condensates are present for finite-width nanoribbons with multiple subbands. We find that superfluidity is predominately in the strongly-coupled BEC and BCS-BEC crossover regimes, with large superfluid gaps up to 100 meV and beyond. When the gaps exceed the subband spacing, there is significant mixing of the subbands, a rounding of the shape resonances, and a resulting reduction in the one-dimensional nature of the system.

Characterization of the size and position of electron-hole puddles at a graphene p-n junction

The effect of an electron-hole puddle on the electrical transport when governed by snake states in a bipolar graphene structure is investigated. Using numerical simulations we show that information on the size and position of the electron-hole puddle can be obtained using the dependence of the conductance on magnetic field and electron density of the gated region. The presence of the scatterer disrupts snake state transport which alters the conduction pattern. We obtain a simple analytical formula that connects the position of the electron-hole puddle with features observed in the conductance. The size of the electron-hole puddle is estimated from the magnetic field and gate potential that maximizes the effect of the puddle on the electrical transport.