Abstract: We study topological phase transitions and topological quantum field effect transistor in monolayer molybdenum disulfide (MoS2) using a two-band Hamiltonian model. Without considering the quadratic (q(2)) diagonal term in the Hamiltonian, we show that the phase diagram includes quantum anomalous Hall effect, quantum spin Hall effect, and spin quantum anomalous Hall effect regions such that the topological Kirchhoff law is satisfied in the plane. By considering the q(2) diagonal term and including one valley, it is shown that MoS2 has a non-trivial topology, and the valley Chern number is non-zero for each spin. We show that the wave function is (is not) localized at the edges when the q(2) diagonal term is added (deleted) to (from) the spin-valley Dirac mass equation. We calculate the quantum conductance of zigzag MoS2 nanoribbons by using the nonequilibrium Green function method and show how this device works as a field effect topological quantum transistor.

Abstract: In van der Waals heterostructures consisting of graphene and a substrate, lattice mismatch often leads to a moire pattern with a huge supercell, preventing its treatment within first- principles calculations. Previous theoretical works considered mostly simple stacking models such as AB, AA with straining the lattice of graphene to match that of the substrate. Here, we propose a moire superlattice build from graphene and porous graphene or graphyne like monolayers, having a lower interlayer binding energy, needing little strain in order to match the lattices. In contrast to the results from the simple stacking models, the present ab initio calculations for the moire superlattices show different properties in lattice structure, energy, and band structures. For example, the Dirac cone at the K point is preserved and a linear energy dispersion near the Fermi level is obtained.

Abstract: The electronic and vibrational properties of germanane and fluorinated germanene are studied within density functional theory (DFT) and density functional perturbation theory frameworks. Different structural configurations of germanane and fluorinated germanene are investigated. The energy difference between the different configurations are consistently smaller than the energy of thermal fluctuations for all the analyzed DFT functionals LDA, GGA, and hybrid functionals, which implies that, in principle, it is possible to find these different configurations in different regions of the sample as minority phases or local defects. We calculate the Raman and infrared spectra for these configurations by using ab initio calculations and compare it with available experimental spectra for germanane. Our results show the presence of minority phases compatible with the configurations analyzed in this work. As these low energy configurations are metastable the present work shows that the synthesis of these energy competing phases is feasible by selectively changing the synthesis conditions, which is an opportunity to expand in this way the availability of new two-dimensional compounds.

Abstract: The energy spectrum and local current patterns in graphene quantum dots (QD) are investigated for different geometries in the presence of an external perpendicular magnetic field. Our results demonstrate that, for specific geometries and edge configurations, the QD exhibits vortex and anti-vortex patterns in the local current density, in close analogy to the vortex patterns observed in the probability density current of semiconductor QD, as well as in the order parameter of mesoscopic superconductors.

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: Recent advances in manufacturing arrays of artificial pinning sites, i.e., antidots, blind holes and magnetic dots, allowed an effective control of magnetic flux in superconductors. An array of magnetic bars deposited on top of a superconducting film was shown to display different pinning regimes depending on the direction of the in-plane magnetization of the bars. Changing the sign of their magnetization results in changes in the induced magnetic pinning potentials. By numerically solving the time-dependent Ginzburg-Landau equations in a superconducting film with periodic arrays of zigzag-arranged magnetic bars, we revealed various flux dynamics regimes. In particular, we demonstrate flux pinning and flux flow, depending on the direction of the magnetization of the magnetic bars. Remarkably, the revealed different flux-motion regimes are associated with different mechanisms of vortex-antivortex dynamics. For example, we found that for an 'antiparallel' configuration of magnetic bars this dynamics involves a repeating vortex-antivortex generation and annihilation. We show that the depinning transition and the onset of flux flow can be manipulated by the magnetization of the bars and the geometry of the array. This provides an effective control of the depinning critical current that can be useful for possible fluxonics applications.

Abstract: We study vortex configurations in mesoscopic superconducting squares. Our theoretical approach is based on the analytical solution of the London equation using the Green's function method. The potential energy landscape found is then used in Langevin-type molecular-dynamics simulations to obtain stable vortex configurations. We show that the filling rules for vortices in squares with increasing applied magnetic field can be formulated, although in a different manner than in disks, in terms of the formation of vortex 'shells'. We discuss metastable states and the stability of the vortex configurations found with respect to variations of the material parameters and deformations of the shape of the sample.

Abstract: The dynamics of Josephson vortices (fluxons) in artificial stacks of superconducting-normal-superconducting Josephson junctions is investigated using the anisotropic time-dependent Ginzburg-Landau theory in the presence of a square/rectangular array of pinning centers (holes). For small values of the applied drive, fluxons in different junctions move out of phase, forming a periodic triangular lattice. A rectangular lattice of moving fluxons is observed at larger currents, which is in agreement with previous theoretical predictions (Koshelev and Aranson 2000 Phys. Rev. Lett. 85 3938). This 'superradiant' flux-flow state is found to be stable in a wide region of applied current. The stability range of this ordered state is considerably larger than the one obtained for the pinning-free sample. Clear commensurability features are observed in the current-voltage characteristics of the system with pronounced peaks in the critical current at (fractional) matching fields. The effect of density and strength of the pinning centers on the stability of the rectangular fluxon lattice is discussed. Predicted synchronized motion of fluxons in the presence of ordered pinning can be detected experimentally using the rf response of the system, where enhancement of the Shapiro-like steps is expected due to the synchronization.

Abstract: Structural, mechanical and electronic properties of two-dimensional single-layer hexagonal structures in the (111) crystal plane of IIIAs-ZnS systems (III = B, Ga and In) are studied by first-principles calculations based on density functional theory (DFT). Elastic and phonon dispersion relation display that 2D h-IIIAs systems (III = B, Ga and In) are both mechanical and dynamically stable. Electronic structures analysis show that the semiconducting nature of the 3D-IIIAs compounds is retained by their 2D single layer counterpart. Furthermore, density of states reveals the influence of sigma and pi bonding in the most stable geometry (planar or buckled) for 2D h-IIIAs systems. Calculations of elastic constants show that the Young's modulus, bulk modulus and shear modulus decrease for 2D h-IIIAs binary compounds as we move down on the group of elements of the periodic table. In addition, as the bond length between the neighboring cation-anion atoms increases, the 2D h-IIIAs binary compounds display less stiffness and more plasticity. Our findings can be used to understand the contribution of the r and p bonding in the most stable geometry (planar or buckled) for 2D h-IIIAs systems. Structural and electronic properties of h-IIIAs systems as a function of the number of layers have been also studied. It is shown that h-BAs keeps its planar geometry while both h-GAs and h-InAs retained their buckled ones obtained by their single layers. Bilayer h-IIIAs present the same bandgap nature of their counterpart in 3D. As the number of layers increase from 2 to 4, the bandgap width for layered h-IIIAs decreases until they become semimetal or metal. Interestingly, these results are different to those found for layered h-GaN. The results presented in this study for single and few-layer h-IIIAs structures could give some physical insights for further theoretical and experimental studies of 2D h-IIIV-like systems. (C) 2017 Elsevier B.V. All rights reserved.

Abstract: Oxocarbons, known for over two centuries, have recently revealed a long-awaited facet: two-dimensional crystalline structures. Employing an intelligent global optimization algorithm (IGOA) alongside densityfunctional calculations, we unearthed a quasi -flat oxocarbon (C 6 0 6 ), featuring an oxygen -decorated hole, and a novel 3D-borocarbon. Comparative analyses with recently synthesized isostructures, such as 2D -porous carbon nitride (C 6 N 6 ) and 2D -porous boroxine (B 6 0 6 ), highlight the unique attributes of these compounds. All structures share a common stoichiometry of X 6 Y 6 (which we call COF-66), where X = B, C, and Y = B, N, O (with X not equal Y), exhibiting a 2D -crystalline structure, except for borocarbon C 6 B 6 , which forms a 3D crystal. In our comprehensive study, we conducted a detailed exploration of the electronic structure of X 6 Y 6 compounds, scrutinizing their thermodynamic properties and systematically evaluating phonon stability criteria. With expansive surface areas, diverse pore sizes, biocompatibility, pi-conjugation, and distinctive photoelectric properties, these structures, belonging to the covalent organic framework (COF) family, present enticing prospects for fundamental research and hold potential for biosensing applications.

Abstract: Using density-functional theory in combination with the nonequilibrium Green's function formalism, we study the effect of iodide/chloride and iodide/bromide mixing on the electronic transport in lead based organometallic perovskite CH3NH3PbI3, which is known to be an effective tool to tune the electronic and optical properties of such materials. We found that depending on the level and position of the halide mixing, the electronic transport can be increased by more than a factor of 4 for a given voltage biasing. The largest current is observed for small concentration of bromide substitutions located at the equatorial sites. However, full halide substitution has a negative effect on the transport properties of this material: the current drops by an order of magnitude for both CH3NH3PbCl3 and CH3NH3PbBr3 samples. (C) 2015 Elsevier B.V. All rights reserved.

Abstract: On the basis of the time-dependent Ginzburg-Landau equation we studied the dynamics of the superconducting condensate in a wide two-dimensional sample in the presence of a perpendicular magnetic field and applied current. We could identify two critical currents: the current at which the pure superconducting state becomes unstable (J(c2)(1)) and the current at which the system transits from the resistive state to the superconducting state (J(c1) < J(c2)). The current J(c2) decreases monotonically with external magnetic field, while J(c1) exhibits a maximum at H*. For sufficient large magnetic fields the hysteresis disappears and J(c1) = J(c2) = Jc. In this high magnetic field region and for currents close to Jc the voltage appears as a result of the motion of separate vortices. With increasing current the moving vortices form,channels' with suppressed order parameter along which the vortices can move very fast. This leads to a sharp increase of the voltage. These 'channels' resemble in some respect the phase slip lines which occur at zero magnetic field. (C) 2004 Elsevier B.V. All rights reserved.