Abstract: A pseudomagnetic field kink can be realized along a graphene nanoribbon using strain engineering. Electron transport along this kink is governed by snake states that are characterized by a single propagation direction. Those pseudomagnetic fields point towards opposite directions in the K and K' valleys, leading to valley polarized snake states. In a graphene nanoribbon with armchair edges this effect results in a valley filter that is based only on strain engineering. We discuss how to maximize this valley filtering by adjusting the parameters that define the stress distribution along the graphene ribbon.

Abstract: Similar to the spontaneous formation of colonies of bacteria, flocks of birds, or schools of fish, “living crystals” can be formed by artificial self-propelled particles such as Janus colloids. Unlike usual solids, these “crystals” are far from thermodynamic equilibrium. They fluctuate in time forming a crystalline structure, breaking apart and re-forming again. We propose a method to stabilize living crystals by applying a weak confinement potential that does not suppress the ability of the particles to perform self-propelled motion, but it stabilizes the structure and shape of the dynamical clusters. This gives rise to such configurations of living crystals as “living shells” formed by Janus colloids. Moreover, the shape of the stable living clusters can be controlled by tuning the potential strength. Our proposal can be verified experimentally with either artificial microswimmers such as Janus colloids, or with living active matter.

Abstract: 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.

Abstract: We study the influence of confinement on the atomic collapse due to a Coulomb impurity placed at the center of a graphene quantum dot of radius R. We apply the zigzag or infinite-mass boundary condition and consider both a point-size and a finite-size impurity. As a function of the impurity strength Za, the energy spectra are discrete. In the case of the zigzag boundary condition, the degenerate (with respect to the angular momentum m) zero-energy levels are pulled down in energy as Z alpha increases, and they remain below epsilon = – Z alpha. Our results show that the energy levels exhibit a 1/R dependence in the subcritical regime [Z alpha < |km + 1/2|, k = 1 (-1) for the K (K') valley]. In the supercritical regime (Z alpha > |km + 1/2|) we find a qualitatively very different behavior where the levels decrease as a function of R in a nonmonotonic manner. While the valley symmetry is preserved in the presence of the impurity, we find that the impurity breaks electron-hole symmetry. We further study the energy spectrum of zigzag quantum dots in gapped graphene. Our results show that as the gap increases, the lowest electron states are pushed into the gap by the impurity.

Abstract: To unravel the effects of phonon confinement, the influence of size and morphology on the atomic vibrations is investigated in Sn nanoislands and cluster-assembled films. Nuclear resonant inelastic x-ray scattering is used to probe the phonon densities of states of the Sn nanostructures which show significant broadening of the features compared to bulk phonon behavior. Supported by ab initio calculations, the broadening is attributed to phonon scattering and can be described within the damped harmonic oscillator model. Contrary to the expectations based on previous research, the appearance of high-energy modes above the cutoff energy is not observed. From the thermodynamic properties extracted from the phonon densities of states, it was found that grain boundary Sn atoms are bound by weaker forces than bulk Sn atoms.

Abstract: We report on work function measurements on graphene, which is exfoliated over a predefined array of wells in silicon oxide, by Kelvin probe force microscopy operating in a vacuum. The obtained graphene sealed microchambers can support large pressure differences, providing controllable stretching of the nearly impermeable graphene membranes. These measurements allow detecting variations of the work function induced by the mechanical stresses in the suspended graphene where the work function varies linearly with the strain and changes by 62 +/- 2 meV for 1 percent of strain. Our related ab initio calculations result in a work function variation that is a factor of 1.4 larger than the experimental value. The limited discrepancy between the theory and the experiment can be accounted for by a charge transfer from the unstrained to the strained graphene regions. Published by AIP Publishing.

Abstract: Graphene nanosheets have recently received a revival of interest as a new class of ultrathin, high-flux, and energy-efficient sieving membranes because of their unique two-dimensional and atomically thin structure, good flexibility, and outstanding mechanical properties. However, for practical applications of graphene for advanced water purification and desalination technologies, the creation of well controlled, high-density, and subnanometer diameter pores becomes a key factor. Here, we conduct reactive force-field molecular dynamics simulations to study the effect of external strain on nanopore creation in the suspended graphene by bombardment with Si clusters. Depending on the size and energy of the clusters, different kinds of topography were observed in the graphene sheet. In all the considered conditions, tensile strain results in the creation of nanopores with regular shape and smooth edges. On the contrary, compressive strain increases the elastic response of graphene to irradiation that leads to the formation of net-like defective structures with predominantly carbon atom chains. Our findings show the possibility of creating controlled nanopores in strained graphene by bombardment with Si clusters. Published by AIP Publishing.

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

Abstract: Two-dimensional group-V elemental materials have attracted widespread attention due to their nonzero band gap while displaying high electron mobility. Using first-principles calculations, we propose a series of new elemental bilayers with group-V elements (Bi, Sb, As). Our study reveals the dynamical stability of four-, six-, and eight-atom ring structures, demonstrating their possible coexistence in such bilayer systems. The proposed structures for Sb and As are large-gap semiconductors that are potentially interesting for applications in future nanodevices. The Bi structures have nontrivial topological properties with a direct nontrivial band gap. The nontrivial gap is shown to arise from a band inversion at the Brillouin zone center due to the strong intrinsic spin-orbit coupling in Bi atoms. Moreover, we demonstrate the possibility of tuning the properties of these materials by enhancing the ratio of six-atom rings to four-and eight-atom rings, which results in wider nontrivial band gaps and lower formation energies.

Abstract: Recently proposed accurate correlation energies are used to determine the phase diagram of strongly coupled electron-hole graphene bilayers. The control parameters of the phase diagram are the charge carrier density and the insulating barrier thickness separating the bilayers. In addition to the electron-hole superfluid phase we find two new inhomogeneous ground states, a one dimensional charge density wave phase and a coupled electron-hole Wigner crystal. The elementary crystal structure of bilayer graphene plays no role in generating these new quantum phases, which are completely determined by the electrons and holes interacting through the Coulomb interaction. The experimental parameters for the new phases lie within attainable ranges and therefore coupled electron-hole bilayer graphene presents itself as an experimental system where novel emergent many-body phases can be realized.

Abstract: Starting from a Hamiltonian with anharmonic coupling between in-plane acoustic displacements and outof-plane (flexural) modes, we derived coupled equations of motion for in-plane displacements correlations and flexural mode density fluctuations. Linear response theory and time-dependent thermal Green's functions techniques are applied in order to obtain different response functions. As external perturbations we allow for stresses and thermal heat sources. The displacement correlations are described by a Dyson equation where the flexural density distribution enters as an additional perturbation. The flexural density distribution satisfies a kinetic equation where the in-plane lattice displacements act as a perturbation. In the hydrodynamic limit this system of coupled equations is at the basis of a unified description of elastic and thermal phenomena, such as isothermal versus adiabatic sound motion and thermal conductivity versus second sound. The general theory is formulated in view of application to graphene, two-dimensional h-BN, and 2H-transition metal dichalcogenides and oxides.

Abstract: We propose a tight-binding (TB) model, that includes spin-orbit coupling (SOC), to describe the electronic properties of methyl-substituted germanane (GeCH3). This model gives an electronic spectrum in agreement with first principle results close to the Fermi level. Using the Z(2) formalism, we show that a topological phase transition from a normal insulator (NI) to a quantum spin Hall (QSH) phase occurs at 11.6% biaxial tensile strain. The sensitivity of the electronic properties of this system on strain, in particular its transition to the topological insulating phase, makes it very attractive for applications in strain sensors and other microelectronic applications.

Abstract: The electron energy spectrum in monolayer-bilayer-monolayer and in bilayer-monolayer-bilayer graphene structures is investigated and the effects of a perpendicular magnetic field and electric bias are studied. Different types of monolayer-bilayer interfaces are considered as zigzag (ZZ) or armchair (AC) junctions which modify considerably the bulk Landau levels (LLs) when the spectra are plotted as a function of the center coordinate of the cyclotron orbit. Far away from the two interfaces, one obtains the well-known LLs for extended monolayer or bilayer graphene. The LL structure changes significantly at the two interfaces or junctions where the valley degeneracy is lifted for both types of junctions, especially when the distance between them is approximately equal to the magnetic length. Varying the nonuniform bias and the width of this junction-to-junction region in either structure strongly influence the resulting spectra. Significant differences exist between ZZ and AC junctions in both structures. The densities of states (DOSs) for unbiased structures are symmetric in energy whereas those for biased structures are asymmetric. An external bias creates interface LLs in the gaps between the LLs of the unbiased system in which the DOS can be quite small. Such a pattern of LLs can be probed by scanning tunneling microscopy.

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.