Abstract: Dirac plasmons in graphene are very sensitive to the dielectric properties of the environment. We show that this can be used to probe the structure and composition of van der Waals heterostructures (vdWh) put underneath a single graphene layer. In order to do so, we assess vdWh composed of hexagonal boron nitride and different types of transition metal dichalcogenides (TMDs). By performing realistic simulations that account for the contribution of each layer of the vdWh separately and including the importance of the substrate phonons, we show that one can achieve single-layer resolution by investigating the nonlocal nature of the Dirac plasmon-polaritons. The composition of the vdWh stack can be inferred from the plasmon-phonon coupling once it is composed by more than two TMD layers. Furthermore, we show that the bulk character of TMD stacks for plasmonic screening properties in the terahertz regime is reached only beyond 100 layers.

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

Abstract: We show that a domain wall separating single-layer graphene and AA-stacked bilayer graphene (AA-BLG) can be used to generate highly collimated electron beams which can be steered by a magnetic field. Two distinct configurations are studied, namely, locally delaminated AA-BLG and terminated AA-BLG whose terminal edge types are assumed to be either zigzag or armchair. We investigate the electron scattering using semiclassical dynamics and verify the results independently with wave-packet dynamics simulations. We find that the proposed system supports two distinct types of collimated beams that correspond to the lower and upper cones in AA-BLG. Our computational results also reveal that collimation is robust against the number of layers connected to AA-BLG and terminal edges.

Abstract: Recent experimental measurements of light absorption in few-layer black phosphorus (BP) revealed a series of high and sharp peaks, interspersed by pairs of lower and broader features. Here, we propose a theoretical model for these excitonic states in few-layer BP within a continuum approach for the in-plane degrees of freedom and a tight-binding approximation that accounts for interlayer couplings. This yields excitonic transitions between different combinations of the subbands created by the coupled BP layers, which leads to a series of high and low oscillator strength excitonic states, consistent with the experimentally observed bright and dark exciton peaks, respectively. The main characteristics of such subband exciton states, as well as the possibility to control their energies and oscillator strengths via applied electric and magnetic fields, are discussed, towards a full understanding of the excitonic spectrum of few-layer BP and its tunability.

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: The McMillan-Ginzburg-Landau (MGL) model for charge density waves (CDW) is employed in a systematic phenomenological study of the different phases that have been probed in recent experiments involving transition metal dichalcogenides. We implemented an efficient imaginary time evolution method to solve the MGL equations, which enabled us to investigate the role of different coupling parameters on the CDW patterns and to perform calculations with different energy functionals that lead to several experimentally observed singularities in the CDW phase profiles. In particular, by choosing the appropriate energy functionals, we were able to obtain phases that go beyond the well-known periodic phase slips (discommensurations), exhibiting also topological defects (i.e., vortex-antivortex pairs), domain walls where the CDW order parameter is suppressed, and even CDW with broken rotational symmetry. Finally, we briefly discuss the effect of these different CDW phases on the profile and critical temperature of the competing superconducting state.

Abstract: The time evolution of a wave packet injected into a semiconductor quantum ring is investigated in order to obtain the transmission and reflection probabilities. Within the effective-mass approximation, the time-dependent Schrödinger equation is solved for a system with nonzero width of the ring and leads and finite potential-barrier heights, where we include smooth lead-ring connections. In the absence of a magnetic field, an analysis of the projection of the wave function over the different subband states shows that when the injected wave packet is within a single subband, the junction can scatter this wave packet into different subbands but remarkably at the second junction the wave packet is scattered back into the subband state of the incoming wave packet. If a magnetic field is applied perpendicularly to the ring plane, transmission and reflection probabilities exhibit Aharonov-Bohm (AB) oscillations and the outgoing electrons may end up in different subband states from those of the incoming electrons. Localized impurities, placed in the ring arms, influence the AB oscillation period and amplitude. For a single impurity or potential barrier of sufficiently strong strength, the period of the AB oscillations is halved while for two impurities localized in diametrically opposite points of the ring, the original AB period is recovered. A theoretical investigation of the confined states and time evolution of wave packets in T wires is also made, where a comparison between this system and the lead-ring junction is drawn.

Abstract: The time evolution of a wave packet in strained graphene is studied within the tight-binding model and continuum model. The effect of an external magnetic field, as well as a strain-induced pseudomagnetic field, on the wave-packet trajectories and zitterbewegung are analyzed. Combining the effects of strain with those of an external magnetic field produces an effective magnetic field which is large in one of the Dirac cones, but can be practically zero in the other. We construct an efficient valley filter, where for a propagating incoming wave packet consisting of momenta around the K and K' Dirac points, the outgoing wave packet exhibits momenta in only one of these Dirac points while the components of the packet that belong to the other Dirac point are reflected due to the Lorentz force. We also found that the zitterbewegung is permanent in time in the presence of either external or strain-induced magnetic fields, but when both the external and strain-induced magnetic fields are present, the zitterbewegung is transient in one of the Dirac cones, whereas in the other cone the wave packet exhibits permanent spatial oscillations.

Abstract: We describe a semianalytic approach to the two-band Ginzburg-Landau theory, which predicts the behavior of vortices in two-band superconductors. We show that the character of the short-range vortex-vortex interaction is determined by the sign of the normal domain-superconductor interface energy, in analogy with the conventional differentiation between type I and type II superconductors. However, we also show that the long-range interaction is determined by a modified Ginzburg-Landau parameter κ*, different from the standard κ of a bulk superconductor. This opens the possibility for nonmonotonic vortex-vortex interaction, which is temperature dependent, and can be further tuned by alterations of the material on the microscopic scale.

Abstract: The Dirac equation is solved for triangular and hexagonal graphene quantum dots for different boundary conditions in the presence of a perpendicular magnetic field. We analyze the influence of the dot size and its geometry on their energy spectrum. A comparison between the results obtained for graphene dots with zigzag and armchair edges, as well as for infinite-mass boundary condition, is presented and our results show that the type of graphene dot edge and the choice of the appropriate boundary conditions have a very important influence on the energy spectrum. The single-particle energy levels are calculated as a function of an external perpendicular magnetic field that lifts degeneracies. Comparing the energy spectra obtained from the tight-binding approximation to those obtained from the continuum Dirac equation approach, we verify that the behavior of the energies as a function of the dot size or the applied magnetic field are qualitatively similar, but in some cases quantitative differences can exist.

Abstract: We theoretically demonstrate that the transport inefficiency recently found experimentally for branched-out mesoscopic networks can also be observed in a quantum ring of finite width with an attached central horizontal branch. This is done by investigating the time evolution of an electron wave packet in such a system. Our numerical results show that the conductivity of the ring does not necessary improve if one adds an extra channel. This ensures that there exists a quantum analog of the Braess paradox, originating from quantum scattering and interference.

Abstract: We calculate the excitonic spectrum of few-layer black phosphorus by direct diagonalization of the effective mass Hamiltonian in the presence of an applied in-plane electric field. The strong attractive interaction between electrons and holes in this system allows one to investigate the Stark effect up to very high ionizing fields, including also the excited states. Our results show that the band anisotropy in black phosphorus becomes evident in the direction-dependent field-induced polarizability of the exciton.

Abstract: Within a tight binding approach we investigate the energy levels of hexagonal and triangular bilayer graphene (BLG) quantum dots (QDs) with zigzag and armchair edges. We study AA- and AB-(Bernal) stacked BLG QDs and obtain the energy levels in both the absence and the presence of a perpendicular electric field (i.e., biased BLG QDs). Our results show that the size dependence of the energy levels is different from that of monolayer graphene QDs. The energy spectrum of AB-stacked BLG QDs with zigzag edges exhibits edge states which spread out into the opened energy gap in the presence of a perpendicular electric field. We found that the behavior of these edges states is different for the hexagonal and triangular geometries. In the case of AA-stacked BLG QDs, the electron and hole energy levels cross each other in both cases of armchair and zigzag edges as the dot size or the applied bias increases.

Abstract: Trions and biexcitons in anisotropic two-dimensional materials are investigated within an effective mass theory. Explicit results are obtained for phosphorene and arsenene, materials that share features such as a direct quasiparticle gap and anisotropic conduction and valence bands. Trions are predicted to have remarkably high binding energies and an elongated electron-hole structure with a preference for alignment along the armchair direction, where the effective masses are lower. We find that biexciton binding energies are also notably large, especially for monolayer phosphorene, where they are found to be twice as large as those for typical monolayer transition metal dichalcogenides.

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.