Abstract: Coupled electron-hole states are realized in a system consisting of a combination of an electrostatic potential barrier and ring-shaped potential well, which resembles a circular dipole. A perpendicular magnetic field induces confined states inside the Landau gaps which are mainly located at the barrier or ring. Hybridizations between the barrier and ring states are seen as anticrossings in the energy spectrum. As a consequence, the energy levels show an oscillating dependence on the electrostatic potential strength in combination with an oscillating migration of the wave functions between the barrier and ring. At the anticrossing points the quantum state consists of a mixture of electron and hole. The present system mimics closely the behavior of a relativistic dipole on gapped graphene.

Abstract: Considering a noncentrosymmetric pinning texture composed of a square array of triangular holes, the magnetic flux penetration and expulsion are investigated experimentally and theoretically. A direct visualization of the magnetic landscape obtained using a magneto-optical technique on a Nb film is complemented by a multiscale numerical modeling. This combined approach allows the magnetic flux dynamics to be identified from the single flux quantum limit up to the macroscopic electromagnetic response. Within the theoretical framework provided by time-dependent Ginzburg-Landau simulations, an estimation of the in-plane current anisotropy is obtained and its dependence with the radius of the curvature of hole vertices is addressed. These simulations show that current crowding plays an important role in channeling the flux motion, favoring hole-to-hole flux hopping rather than promoting interstitial flux displacement in between the holes. The resulting anisotropy of the critical current density gives rise to a distinct pattern of discontinuity lines for increasing and decreasing applied magnetic fields, in sharp contrast to the invariable patterns reported for centrosymmetric pinning potentials. This observation is partially accounted for by the rectification effect, as demonstrated by finite-element modeling. At low temperatures, where magnetic field penetration is dominated by thermomagnetic instabilities, highly directional magnetic flux avalanches with a fingerlike shape are observed to propagate along the easy axis of the pinning potential. This morphology is reproduced by numerical simulations. Our findings demonstrate that anisotropic pinning landscapes and, in particular, ratchet potentials produce subtle modifications to the critical state field profile that are reflected in the distribution of discontinuity lines.

Abstract: Chiral p-wave superconductors in applied magnetic field can exhibit more complex topological defects than just conventional superconducting vortices, due to the two-component order parameter (OP) and the broken time-reversal symmetry. We investigate the electronic properties of those exotic states, some of which contain clusters of one-component vortices in chiral components of the OP and/or exhibit skyrmionic character in the relative OP space, all obtained as a self-consistent solution of the microscopic Bogoliubov-de Gennes equations. We reveal the link between the local density of states (LDOS) of the novel topological states and the behavior of the chiral domain wall between the OP components, enabling direct identification of those states in scanning tunneling microscopy. For example, a skyrmion always contains a closed chiral domain wall, which is found to be mapped exactly by zero-bias peaks in LDOS. Moreover, the LDOS exhibits electron-hole asymmetry, which is different from the LDOS of conventional vortex states with same vorticity. Finally, we present the magnetic field and temperature dependence of the properties of a skyrmion, indicating that this topological defect can be surprisingly large in size, and can be pinned by an artificially indented nonsuperconducting closed path in the sample. These features are expected to facilitate the experimental observation of skyrmionic states, thereby enabling experimental verification of chirality in emerging superconducting materials.

Abstract: By using density functional theory (DFT) calculations, we predict a puckered, dynamically stable Janus single-layer black arsenic-phosphorus (b-AsP), which is composed of two different atomic sublayers, arsenic and phosphorus atoms. The calculated phonon spectrum reveals that Janus single-layer b-AsP is dynamically stable with either pure or coupled optical phonon branches arising from As and P atoms. The calculated Raman spectrum indicates that due to the relatively strong P-P bonds, As atoms have no contribution to the highfrequency optical vibrations. In addition, the orientation-dependent isovolume heat capacity reveals anisotropic contributions of LA and TA phonon branches to the low-temperature thermal properties. Unlike pristine single layers of b-As and b-P, Janus single-layer b-AsP exhibits additional out-of-plane asymmetry which leads to important consequences for its electronic, optical, and elastic properties. In contrast to single-layer b-As, Janus single-layer b-AsP is found to possess a direct band gap dominated by the P atoms. Moreover, real and imaginary parts of the dynamical dielectric function, including excitonic effects, reveal the highly anisotropic optical feature of the Janus single-layer. A tight-binding (TB) model is also presented for Janus single-layer b-AsP, and it is shown that, with up to seven nearest hoppings, the TB model reproduces well the DFT band structure in the low-energy region around the band gap. This TB model can be used in combination with the Green's function approach to study, e.g., quantum transport in finite systems based on Janus single-layer b-AsP. Furthermore, the linear-elastic properties of Janus single-layer b-AsP are investigated, and the orientation-dependent in-plane stiffness and Poisson ratio are calculated. It is found that the Janus single layer exhibits strong in-plane anisotropy in its Poisson ratio much larger than that of single-layer b-P. This Janus single layer is relevant for promising applications in optical dichroism and anisotropic nanoelasticity.

Abstract: By taking account of the electric-field-induced charge screening, a self-consistent calculation within the framework of the tight-binding approach is employed to obtain the electronic band structure of gated multilayer phosphorene and the charge densities on the different phosphorene layers. We find charge density and screening anomalies in single-gated multilayer phosphorene and electron-hole bilayers in dual-gated multilayer phosphorene. Due to the unique puckered lattice structure, both intralayer and interlayer charge screenings are important in gated multilayer phosphorene. We find that the electric-field tuning of the band structure of multilayer phosphorene is distinctively different in the presence and absence of charge screening. For instance, it is shown that the unscreened band gap of multilayer phosphorene decreases dramatically with increasing electric-field strength. However, in the presence of charge screening, the magnitude of this band-gap decrease is significantly reduced and the reduction depends strongly on the number of phosphorene layers. Our theoretical results of the band-gap tuning are compared with recent experiments and good agreement is found.

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: By making use of circularly polarized light and electrostatic gating, monolayer molybdenum disulfide (ML – MoS2) can form a platform supporting multiple types of charge carriers. They can be discriminated by their spin, valley index, or whether they are electrons or holes. We investigate the collective properties of those charge carriers and are able to identify distinct plasmon modes. We analyze the corresponding dispersion relation, lifetime, and oscillator strength, and calculate the phase relation between the oscillations in the different components of the plasmon modes. All platforms in ML-MoS2 support a long-wavelength root q plasmon branch at zero kelvins. In addition to this, for an n-component system, n-1 distinct plasmon modes appear as acoustic modes with linear dispersion in the long-wavelength limit. These modes correspond to out-of-phase oscillations in the different fermion liquids and have, although being damped, a relatively long lifetime. Additionally, we also find distinct modes at large wave vectors that are more strongly damped by intraband processes.

Abstract: Atomic scale defects in semiconductors enable their technological applications and realization of different quantum states. Using scanning tunneling microscopy and spectroscopy complemented by ab initio calculations we determine the nature of defects in the anisotropic van der Waals layered semiconductor ReS2. We demonstrate the in-plane anisotropy of the lattice by directly visualizing chains of rhenium atoms forming diamond-shaped clusters. Using scanning tunneling spectroscopy we measure the semiconducting gap in the density of states. We reveal the presence of lattice defects and by comparison of their topographic and spectroscopic signatures with ab initio calculations we determine their origin as oxygen atoms absorbed at lattice point defect sites. These results provide an atomic-scale view into the semiconducting transition metal dichalcogenides, paving the way toward understanding and engineering their properties.

Abstract: Atomic collapse can be observed in graphene because of its large “effective” fine structure constant, which enables this phenomenon to occur for an impurity charge as low as Z(c) similar to 1-2. Here we investigate the effect of the sublattice symmetry on molecular collapse in two spatially separated charge tunable vacancies, which are located on the same (A-A type) or different (A-B type) sublattices. We find that the broken sublattice symmetry: (1) does not affect the location of the main bonding and antibonding molecular collapse peaks, (2) but shifts the position of the satellite peaks, because they are a consequence of the breaking of the local sublattice symmetry, and (3) there are vacancy characteristic collapse peaks that only occur for A-B type vacancies, which can be employed to distinguish them experimentally from the A-A type. As the charge, energy, and separation distance increase, the additional collapse features merge with the main molecular collapse peaks. We show that the spatial distribution around the vacancy site of the collapse states allows us to differentiate the molecular from the frustrated collapse.

Abstract: An external electric field changes the physical properties of polar liquids due to the reorientation of their permanent dipoles. Using molecular dynamics simulations, we predict that an in-plane electric field applied parallel to the channel polarizes water molecules which are confined between two graphene layers, resulting in distinct ferroelectricity and electrical hysteresis. We found that electric fields alter the in-plane order of the hydrogen bonds: Reversing the electric field does not restore the system to the nonpolar initial state, instead a residual dipole moment remains in the system. The square-rhombic structure of 2D ice is transformed into two rhombic-rhombic structures. Our study provides insights into the ferroelectric state of water when confined in nanochannels and shows how this can be tuned by an electric field.

Abstract: Amorphous semiconductors are usually characterized by a low charge carrier mobility, essentially related to their lack of long-range order. The development of such material with higher charge carrier mobility is hence challenging. Part of the issue comes from the difficulty encountered by first-principles simulations to evaluate concepts such as the electron effective mass for disordered systems since the absence of periodicity induced by the disorder precludes the use of common concepts derived from condensed matter physics. In this paper, we propose a methodology based on first-principles simulations that partially solves this problem, by quantifying the degree of delocalization of a wave function and of the connectivity between the atomic sites within this electronic state. We validate the robustness of the proposed formalism on crystalline and molecular systems and extend the insights gained to disordered/amorphous InGaZnO4 and Si. We also explore the properties of p-type oxide semiconductor candidates recently reported to have a low effective mass in their crystalline phases [G. Hautier et al., Nat. Commun. 4, 2292 (2013)]. Although in their amorphous phase none of the candidates present a valence band with delocalization properties matching those found in the conduction band of amorphous InGaZnO4, three of the seven analyzed materials show some potential. The most promising candidate, K2Sn2O3, is expected to possess in its amorphous phase a slightly higher hole mobility than the electron mobility in amorphous silicon.

Abstract: After the discovery of magnetism in monolayer CrI3, the magnetic properties of different 2D materials from the chromium-trihalide family are intuitively assumed to be similar, yielding magnetic anisotropy from the spin-orbit coupling on halide ligands. Here we reveal significant differences between the CrI3 and CrBr3 magnetic monolayers in their magnetic anisotropy, resulting Curie temperature, hysteresis in external magnetic field, and evolution of magnetism with strain, all predominantly attributed to distinctly different interplay of atomic contributions to spin-orbit coupling in two materials.

Abstract: A strong in-plane magnetic field drastically alters the low-energy spectrum of bilayer graphene by separating the parabolic energy dispersion into two linear Dirac cones. The effect of this dramatic change on the transport properties strongly depends on the orientation of the in-plane magnetic field with respect to the propagation direction of the charge carriers and the angle at which they impinge on the electrostatic potentials. For magnetic fields oriented parallel to the potential boundaries an additional propagating mode that results from the splitting into Dirac cones enhances the transmission probability for charge carriers tunneling through the potentials and increases the corresponding conductance. Our results show that the chiral suppression of transmission at normal incidence, reminiscent of bilayer graphene's 2 pi Berry phase, is turned into a chiral enhancement when the magnetic field increases, thus indicating a transition from a bilayer to a monolayer-like system at normal incidence. Further, we find that the typical transmission resonances stemming from confinement in a potential barrier are shifted to higher energy and are eventually transformed into antiresonances with increasing magnetic field. For magnetic fields oriented perpendicular to the potential boundaries we find a very pronounced transition from a bilayer system to two separated monolayer-like systems with Klein tunneling emerging at certain incident angles symmetric around 0, which also leaves a signature in the conductance. For both orientations of the magnetic field, the transmission probability is still correctly described by pseudospin conservation. Finally, to motivate the large in-plane magnetic field, we show that its energy spectrum can be mimicked by specific lattice deformations such as a relative shift of one of the layers. With this equivalence we introduce the notion of an in-plane pseudomagnetic field.

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: A four-band exciton Hamiltonian is constructed starting from the single-particle Dirac Hamiltonian for charge carriers in monolayer transition metal dichalcogenides (TMDs). The angular part of the exciton wave function can be separated from the radial part, in the case of zero center of mass momentum excitons, by exploiting the eigenstates of the total exciton angular momentum operator with which the Hamiltonian commutes. We explain why this approach fails for excitons with finite center of mass momentum or in the presence of a perpendicular magnetic field and present an approximation to resolve this issue. We calculate the (binding) energy and average interparticle distance of different excited exciton states in different TMDs and compare these with results available in the literature. Remarkably, we find that the intervalley exciton ground state in the -/+ K valley has angular momentum j = +/- 1, which is due to the pseudospin of the separate particles. The exciton mass and the exciton Landau levels are calculated and we find that the degeneracy of exciton states with opposite relative angular momentum is altered by a magnetic field.