Abstract: We numerically study the vortex-vortex interaction in multicomponent homogeneous Bose-Einstein condensates within the realm of the Gross-Pitaevskii theory. We provide strong evidence that pairwise vortex interaction captures the underlying mechanisms which determine the geometric configuration of the vortices, such as different lattices in many-vortex states, as well as the bound vortex states with two (dimer) or three (trimer) vortices. Specifically, we discuss and apply our theoretical approach to investigate intra- and intercomponent vortex-vortex interactions in two- and three-component Bose-Einstein condensates, thereby shedding light on the formation of the exotic vortex configurations. These results correlate with current experimental efforts in multicomponent Bose-Einstein condensates and the understanding of the role of vortex interactions in multiband superconductors.

Abstract: We present a systematic study of the energy spectra of graphene quantum rings having different geometries and edge types in the presence of a perpendicular magnetic field. Results are obtained within the tight-binding (TB) and Dirac models and we discuss which features of the former can be recovered by using the approximations imposed by the latter. Energy levels of graphene quantum rings obtained by diagonalizing the TB Hamiltonian are demonstrated to be strongly dependent on the rings geometry and the microscopical structure of the edges. This makes it difficult to recover those spectra by the existing theories that are based on the continuum (Dirac) model. Nevertheless, our results show that both approaches (i.e., TB and Dirac model) may provide similar results, but only for very specific combinations of ring geometry and edge types. The results obtained by a simplified model describing an infinitely thin circular Dirac ring show good agreement with those obtained for hexagonal and rhombus armchair graphene rings within the TB model. Moreover, we show that the energy levels of a circular quantum ring with an infinite mass boundary condition obtained within the Dirac model agree with those for a ring defined by a ring-shaped staggered potential obtained within the TB model.

Abstract: The scattering of a Gaussian wave packet in armchair and zigzag graphene edges is theoretically investigated by numerically solving the time-dependent Schrodinger equation for the tight-binding model Hamiltonian. Our theory allows us to investigate scattering in reciprocal space, and depending on the type of graphene edge we observe scattering within the same valley, or between different valleys. In the presence of an external magnetic field, the well-known skipping orbits are observed. However, our results demonstrate that in the case of a pseudomagnetic field, induced by nonuniform strain, the scattering by an armchair edge results in a nonpropagating edge state.

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: The vortex-vortex interaction potential in bulk superconductors is calculated within the Ginzburg-Landau (GL) theory and is obtained from a numerical solution of a set of two coupled nonlinear GL differential equations for the vector potential and the superconducting order parameter, where the merger of vortices into a giant vortex is allowed. Further, the interaction potentials between a vortex and a giant vortex and between a vortex and an antivortex are obtained for both type-I and type-II superconductors. Our numerical results agree asymptotically with the analytical expressions for large intervortex separations that are available in the literature. We propose empirical expressions valid over the full interaction range, which are fitted to our numerical data for different values of the GL parameter.

Abstract: Within a minimal model, we present analytical expressions for the eigenstates and eigenvalues of carriers confined in quantum rings in monolayer and bilayer graphene. The calculations were performed in the context of the continuum model by solving the Dirac equation for a zero width ring geometry, i.e., by freezing out the carrier radial motion. We include the effect of an external magnetic field and show the appearance of Aharonov-Bohm oscillations and of a nonzero gap in the spectrum. Our minimal model gives insight on the energy spectrum of graphene-based quantum rings and models different aspects of finite width rings.

Abstract: Using the four-band continuum model we derive a general expression for the infinite-mass boundary condition in bilayer graphene. Applying this new boundary condition we analytically calculate the confined states and the corresponding wave functions in a bilayer graphene quantum dot in the absence and presence of a perpendicular magnetic field. Our results for the energy spectrum show an energy gap between the electron and hole states at small magnetic fields. Furthermore the electron (e) and hole (h) energy levels corresponding to the K and K' valleys exhibit the E-K(e(h)) (m) = E-K'(e(h)) (m) symmetry, where m is the angular momentum quantum number. (C) 2014 Elsevier Ltd. All rights reserved.

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

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

Abstract: We review the transmission properties of carriers interacting with potential barriers in graphene. The tunneling of electrons and holes in quantum structures in graphene is found to display features that are in marked contrast with those of other systems. In particular, the interaction between the carriers with electrostatic potential barriers can be related to the propagation of electromagnetic waves in media with negative refraction indices, also known as metamaterials. This behavior becomes evident as one calculates the time evolution of wavepackets propagating across the barrier interface. In addition, we discuss the effect of trigonal warping on the tunneling through potential barriers.

Abstract: The highly tunable band structure of the zero-energy Landau level (zLL) of bilayer graphene makes it an ideal platform for engineering novel quantum states. However, the zero-energy Landau level at high electric fields has remained largely unexplored. Here we present magnetotransport measurements of bilayer graphene in high transverse electric fields. We observe previously undetected Landau level crossings at filling factors nu = -2, 1, and 3 at high electric fields. These crossings provide constraints for theoretical models of the zero-energy Landau level and show that the orbital, valley, and spin character of the quantum Hall states at high electric fields is very different from low electric fields. At high E, new transitions between states at nu = -2 with different orbital and spin polarization can be controlled by the gate bias, while the transitions between nu = 0 -> 1 and nu = 2 -> 3 show anomalous behavior.

Abstract: The efficiency of light emission is a critical performance factor for monolayer transition metal dichalcogenides (1L-TMDs) for photonic applications. While various methods have been studied to compensate for lattice defects to improve the quantum yield (QY) of 1L-TMDs, exciton-exciton annihilation (EEA) is still a major nonradiative decay channel for excitons at high exciton densities. Here, we demonstrate that the combined use of a proximal Au plate and a negative electric gate bias (NEGB) for 1L-WS2 provides a dramatic enhancement of the exciton lifetime at high exciton densities with the corresponding QY enhanced by 30 times and the EEA rate constant decreased by 80 times. The suppression of EEA by NEGB is attributed to the reduction of the defect-assisted EEA process, which we also explain with our theoretical model. Our results provide a synergetic solution to cope with EEA to realize high-intensity 2D light emitters using TMDs.

Abstract: Van der Waals quantum wells, naturally formed in two-dimensional layered materials with nanoscale thickness, possess many inherent advantages over conventional molecular beam epitaxy grown counterparts, and could bring up intriguing physics and applications. However, optical transitions originated from the series of quantized states in these emerging quantum wells are still elusive. Here, we show that multilayer black phosphorus appears to be an excellent candidate for van der Waals quantum wells with well-defined subbands and high optical quality. Using infrared absorption spectroscopy, we probe subband structures of multilayer black phosphorus with tens of atomic layers, revealing clear signatures for optical transitions with subband index as high as 10, far from what was attainable previously. Surprisingly, in addition to allowed transitions, an unexpected series of “forbidden” transitions is also evidently observed, which enables us to determine energy spacings separately for conduction and valence subbands. Furthermore, the linear tunability of subband spacings by temperature and strain is demonstrated. Our results are expected to facilitate potential applications for infrared optoelectronics based on tunable van der Waals quantum wells.

Abstract: The short exciton lifetime and strong exciton-exciton interaction in transition metal dichalcogenides limit the efficiency of exciton emission. Here, the authors show that exciton-exciton interaction in monolayer WS2 can be screened using proximal metal plates, leading to an improved quantum yield. Monolayer transition metal dichalcogenides (1L-TMDs) have tremendous potential as atomically thin, direct bandgap semiconductors that can be used as convenient building blocks for quantum photonic devices. However, the short exciton lifetime due to the defect traps and the strong exciton-exciton interaction in TMDs has significantly limited the efficiency of exciton emission from this class of materials. Here, we show that exciton-exciton interaction in 1L-WS2 can be effectively screened using an ultra-flat Au film substrate separated by multilayers of hexagonal boron nitride. Under this geometry, induced dipolar exciton-exciton interaction becomes quadrupole-quadrupole interaction because of effective image dipoles formed within the metal. The suppressed exciton-exciton interaction leads to a significantly improved quantum yield by an order of magnitude, which is also accompanied by a reduction in the exciton-exciton annihilation (EEA) rate, as confirmed by time-resolved optical measurements. A theoretical model accounting for the screening of the dipole-dipole interaction is in a good agreement with the dependence of EEA on exciton densities. Our results suggest that fundamental EEA processes in the TMD can be engineered through proximal metallic screening, which represents a practical approach towards high-efficiency 2D light emitters.

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.