Abstract: Recently, interlayer and intralayer excitons in transition metal dichalcogenide heterobilayers have been studied both experimentally and theoretically. In spite of a growing interest, these layer-resolved excitons in the presence of external stimuli, such as strain, remain not fully understood. Here, using density-functional theory calculations with many-body effects, we explore the excitonic properties of vertically stacked MoSe2/WSe2 heterobilayer in the presence of in-plane biaxial strain of up to 5%. We calculate the strain dependence of exciton absorption spectrum, oscillator strength, wave function, and binding energy by solving the Bethe-Salpeter equation on top of the standard GW approach. We identify the interlayer and intralayer excitons by analyzing their electron-hole weights and spatial wave functions. We show that with the increase in strain magnitude, the absorption spectrum of the interlayer and intralayer excitons is red-shifted and re-ordered, and the binding energies of these layer-resolved excitons decrease monotonically and almost linearly. We derive the sensitivity of exciton binding energy to the applied strain and find that the intralayer excitons are more sensitive to strain than the interlayer excitons. For instance, a sensitivity of -7.9 meV/% is derived for the intra-MoSe2-layer excitons, which is followed by -7.4 meV/% for the intra-WSe2-layer excitons, and by -4.2 meV/% for the interlayer excitons. Our results indicate that interlayer and intralayer excitons in vertically stacked MoSe2/WSe2 heterobilayer are efficiently tunable by in-plane biaxial strain.

Abstract: Ab initio density functional theory (DFT) and DFT plus coherent potential approximation (DFT + CPA) are employed to reveal, respectively, the effect of in-plane strain and site-diagonal disorder on the electronic structure of cubic boron arsenide (BAs). It is demonstrated that tensile strain and static diagonal disorder both reduce the semiconducting one-particle band gap of BAs, and a V-shaped p-band electronic state emerges – enabling advanced valleytronics based on strained and disordered semiconducting bulk crystals. At biaxial tensile strains close to 15% the valence band lineshape relevant for optoelectronics is shown to coincide with one reported for GaAs at low energies. The role played by static disorder on the As sites is to promote p-type conductivity in the unstrained BAs bulk crystal, consistent with experimental observations. These findings illuminate the intricate and interdependent changes in crystal structure and lattice disorder on the electronic degrees of freedom of semiconductors and semimetals.

Abstract: We detail the effects of Cl and S functionalization on the superconducting properties of layered (bulk) and monolayer niobium carbide (Nb2C) MXene crystals, based on first-principles calculations combined with Eliashberg theory. For bulk layered Nb2CCl2, the calculated superconducting transition temperature (T-c) is in very good agreement with the recently measured value of 6 K. We show that T-c is enhanced to 10 K for monolayer Nb2CCl2, due to an increase in the density of states at the Fermi level, and the corresponding electron-phonon coupling. We further demonstrate feasible gate- and strain-induced enhancements of T-c for both bulk-layered and monolayer Nb2CCl2 crystals, resulting in T-c values of around 38 K. In the S-functionalized Nb2CCl2 crystals, our calculations reveal the importance of phonon softening in understanding their superconducting properties. Finally, we predict that Nb3C2S2 in bulk-layered and monolayer forms is also superconducting, with a T-c of around 28 K. Considering that Nb2C is not superconducting in pristine form, our findings promote functionalization as a pathway towards robust superconductivity in MXenes.

Abstract: A wave-packet time evolution method, based on the split-operator technique, is developed to investigate the scattering of quasiparticles at a normal-superconductor interface of arbitrary profile and shape. As a practical application, we consider a system where low-energy electrons can be described as Dirac particles, which is the case for most two-dimensional materials, such as graphene and transition-metal dichalcogenides. However, the method is easily adapted for other cases such as electrons in few-layer black phosphorus or any Schrodinger quasiparticles within the effective mass approximation in semiconductors. We employ the method to revisit Andreev reflection in mono-, bi-, and trilayer graphene, where specular-and retro-reflection cases are observed for electrons scattered by a steplike superconducting region. The effect of opening a zero-gap channel across the superconducting region on the electron and hole scattering is also addressed, as an example of the versatility of the technique proposed here.

Abstract: We perform the wavepacket dynamics simulation of a graphene-based device where propagating electron trajectories are tamed by an applied magnetic field toward a normal/superconductor interface. The magnetic field controls the incidence angle of the incoming electronic wavepacket at the interface, which results in the tunable electron-hole ratio in the reflected wave function due to the angular dependence of the Andreev reflection. Here, mapped control of the quasiparticle trajectories by the external magnetic field not only defines an experimental probe for fundamental studies of the Andreev reflection in graphene but also lays the foundation for further development of magnetic focusing devices based on nanoengineered superconducting two-dimensional materials.

Abstract: We examine from first principles the flexoelectric properties of phosphorene nanoribbons under mechanical bending along armchair and zigzag directions. In both cases we find that the radial polarization depends linearly on the strain gradient. The flexoelectricity along the armchair direction is over 40% larger than along the zigzag direction. The obtained flexoelectric coefficients of phosphorene are four orders of magnitude larger than those of graphene and comparable to transition metal dichalcogenides. Analysis of charge density shows that the flexoelectricity mainly arises from the pz orbitals of phosphorus atoms. The electron mobilities in bent phosphorene can be enhanced by over 60% along the armchair direction, which is significantly higher than previous reports of mobility tuned by uniaxial strain. Our results indicate phosphorene is a candidate for a two-dimensional material applicable in flexible-electronic devices.

Abstract: Dirac plasmons in a two-dimensional (2D) crystal are strongly affected by the dielectric properties of the environment, due to interaction of their electric field lines with the surrounding medium. Using graphene as a 2D reservoir of free carriers, one can engineer a material configuration that provides an anisotropic environment to the plasmons. In this work, we discuss the physical properties of Dirac plasmons in graphene surrounded by an arbitrary anisotropic dielectric and exemplify how h-BN-based heterostructures can be designed to bear the required anisotropic characteristics. We calculate how dielec-tric anisotropy impacts the spatial propagation of the plasmons and find that an anisotropy-induced plasmon mode emerges, together with a damping pathway, that stem from the out-of-plane off-diagonal elements in the dielectric tensor. Furthermore, we find that one can create hyperbolic plasmons by inher-iting the dielectric hyperbolicity of the designed material environment. Strong control over plasmon propagation patterns can be realized in a similar manner. Finally, we show that in this way one can also control the polarization of the light-matter excitations that constitute the plasmon. Taken together, our results promote the design of the dielectric environment as an effective path to tailor the plasmonic response of graphene on the nanoscopic level.

Abstract: In a thin-film topological insulator (TI), the edge states on two surfaces may couple by quantum tunneling, opening a gap known as the hybridization gap. Controlling the hybridization gap and transport has a variety of potential uses in photodetection and energy-harvesting applications. In this paper, we report the effect of strain, and electric and magnetic field, on the hybridization gap and transport in a thin Bi2Se3 film, investigated within the tight-binding theoretical framework. We demonstrate that vertical compression decreases the hybridization gap, as does tensile in-plane strain. Applying an electric field breaks the inversion symmetry and leads to a Rashba-like spin splitting proportional to the electric field, hence closing and reopening the gap. The influence of a magnetic field on thin-film TI is also discussed, starting from the role of an out-of-plane magnetic field on quantum Hall states. We further demonstrate that the hybridization gap can be controlled by an in-plane magnetic field, and that by applying a sufficiently strong field a quantum phase transition from an insulator to a semimetal can be achieved.

Abstract: A supersolid, a counterintuitive quantum state in which a rigid lattice of particles flows without resistance, has to date not been unambiguously realized. Here we reveal a supersolid ground state of excitons in a double-layer semiconductor heterostructure over a wide range of layer separations outside the focus of recent experiments. This supersolid conforms to the original Chester supersolid with one exciton per supersolid site, as distinct from the alternative version reported in cold-atom systems of a periodic density modulation or clustering of the superfluid. We provide the phase diagram augmented by the supersolid. This new phase appears at layer separations much smaller than the predicted exciton normal solid, and it persists up to a solid-solid transition where the quantum phase coherence collapses. The ranges of layer separations and exciton densities in our phase diagram are well within reach of the current experimental capabilities.

Abstract: We realize a superconducting nanodevice in which vortex traps in the form of an aluminum square are integrated with a Dayem nanobridge. We perform field cooling of the traps arriving to different vortex configurations, dependent on the applied magnetic field, to demonstrate that the switching current of the bridge is highly sensitive to the presence and location of vortices in the trap. Our measurements exhibit unprecedented precision and ability to detect the first and successive vortex entries into all fabricated traps, from few hundred nm to 2 mu m in size. The experimental results are corroborated by Ginzburg-Landau simulations, which reveal the subtle yet crucial changes in the density of the superconducting condensate in the vicinity of the bridge with every additional vortex entry and relocation inside the trap. An ease of integration and simplicity make our design a convenient platform for studying dynamics of vortices in strongly confining geometries, involving a promise to manipulate vortex states electronically with simultaneous in situ control and monitoring.

Abstract: We investigate the coupling factor phi( mu) that quantifies the magnetic flux phi per magnetic moment mu of a point-like magnetic dipole that couples to a superconducting quantum interference device (SQUID). Representing the dipole by a tiny current-carrying (Amperian) loop, the reciprocity of mutual inductances of SQUID and Amperian loop provides an elegant way of calculating phi(mu)(r,e(mu)) vs. position r and orientation e(mu) of the dipole anywhere in space from the magnetic field B-J(r) produced by a supercurrent circulating in the SQUID loop. We use numerical simulations based on London and Ginzburg-Landau theory to calculate phi (mu) from the supercurrent density distributions in various superconducting loop geometries. We treat the far-field regime ( r greater than or similar to a= inner size of the SQUID loop) with the dipole placed on (oriented along) the symmetry axis of circular or square shaped loops. We compare expressions for phi (mu) from simple filamentary loop models with simulation results for loops with finite width w (outer size A > alpha), thickness d and London penetration depth lambda(L )and show that for thin ( d << alpha ) and narrow (w < alpha) loops the introduction of an effective loop size a(eff) in the filamentary loop-model expressions results in good agreement with simulations. For a dipole placed right in the center of the loop, simulations provide an expression phi(mu)(a,A,d,lambda(L)) that covers a wide parameter range. In the near-field regime (dipole centered at small distance z above one SQUID arm) only coupling to a single strip representing the SQUID arm has to be considered. For this case, we compare simulations with an analytical expression derived for a homogeneous current density distribution, which yields excellent agreement for lambda(L)>w,d . Moreover, we analyze the improvement of phi(mu) provided by the introduction of a narrow constriction in the SQUID arm below the magnetic dipole.

Abstract: Magnetic skyrmions are topological spin configurations in materials with chiral Dzyaloshinskii-Moriya interaction (DMI), that are potentially useful for storing or processing information. To date, DMI has been found in few bulk materials, but can also be induced in atomically thin magnetic films in contact with surfaces with large spin-orbit interactions. Recent experiments have reported that isolated magnetic skyrmions can be stabilized even near room temperature in few-atom-thick magnetic layers sandwiched between materials that provide asymmetric spin-orbit coupling. Here we present the minimum-energy path analysis of three distinct mechanisms for the skyrmion collapse, based on ab initio input and the performed atomic-spin simulations. We focus on the stability of a skyrmion in three atomic layers of Co, either epitaxial on the Pt(111) surface or within a hybrid multilayer where DMI nontrivially varies per monolayer due to competition between different symmetry breaking from two sides of the Co film. In laterally finite systems, their constrained geometry causes poor thermal stability of the skyrmion toward collapse at the boundary, which we show to be resolved by designing the high-DMI structure within an extended film with lower or no DMI.

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: We solve analytically the linearized Ginzburg-Landau (GL) equation in the presence of an uniform magnetic field with cylindrical boundary conditions. The solution of the non-linear GL equation is provided as an expansion in the basis of linearized solutions. We present examples of the resulting vortex structure for a solid and perforated superconducting cylinder. (C) 2012 IMACS. Published by Elsevier B.V. All rights reserved.