Abstract: Nonlinear dynamics of Josephson vortices (fluxons) in artificial stacks of superconducting-normal-superconducting Josephson junctions under simultaneously applied time-periodic ac and constant biasing dc currents is studied using the time dependent Ginzburg-Landau formalism with a Lawrence-Doniach extension. At zero external magnetic field and dc biasing current the resistive state of the system is characterized by periodic nucleation and annihilation of fluxon-antifluxon pairs, relative positions of which are determined by the state of neighboring junctions. Due to the mutual repulsive interaction, fluxons in different junctions move out of phase. Their collective motion can be synchronized by adding a small ac component to the biasing dc current. Coherent motion of fluxons is observed for a broad frequency range of the applied drive. In the coherent state the maximal output voltage, which is proportional to the number of junctions in the stack, is observed near the characteristic frequency of the system determined by the crossing of the fluxons across the sample. However, in this frequency range the dynamically synchronized state has an alternative-a less ordered state with smaller amplitude of the output voltage. Collective behavior of the junctions is strongly affected by the sloped sidewalls of the stack. Synchronization is observed only for weakly trapezoidal cross sections, whereas irregular motion of fluxons is observed for larger slopes of the sample edge.

Abstract: Flux entry in type-II superconductors without prominent symmetry is a first-order phase transition, where flux enters conventionally gradual in units of a flux quantum. Here we show that neither is necessarily the case in a mesoscopic superconducting sphere with a perforation. In axially applied magnetic field, vortices initially occupy the hole, and can oppose further flux entry in the sample. As a result, multiple-quanta flux entry is found at significantly higher field, and it can manifest as a second-order transition due to suppressed geometric barrier at the equatorial belt of the sample. At high fields a new state is found, with gradually destroyed condensate from the equator inwards, the exact opposite of surface superconductivity.

Abstract: Using first-principles calculations, we investigate the electronic, mechanical, and optical properties of monolayer WTe2. Atomic structure and ground state properties of monolayer WTe2 (T-d phase) are anisotropic which are in contrast to similar monolayer crystals of transition metal dichalcogenides, such as MoS2, WS2, MoSe2, WSe2, and MoTe2, which crystallize in the H-phase. We find that the Poisson ratio and the in-plane stiffness is direction dependent due to the symmetry breaking induced by the dimerization of the W atoms along one of the lattice directions of the compound. Since the semimetallic behavior of the T-d phase originates from this W-W interaction (along the a crystallographic direction), tensile strain along the dimer direction leads to a semimetal to semiconductor transition after 1% strain. By solving the Bethe-Salpeter equation on top of single shot G(0)W(0) calculations, we predict that the absorption spectrum of T-d-WTe2 monolayer is strongly direction dependent and tunable by tensile strain. (C) 2016 AIP Publishing LLC.

Abstract: Using the effective continuum model, the electron energy spectrum of gated bilayer graphene with a step-like region of decoupled graphene layers at the edge of the sample is studied. Different types of coupled-decoupled interfaces are considered, i.e., zigzag (ZZ) and armchair junctions, which result in significant different propagating states. Two non-valley-polarized conducting edge states are observed for ZZ type, which are mainly located around the ZZ-ended graphene layers. Additionally, we investigated both BA-BA and BA-AB domain walls in the gated bilayer graphene within the continuum approximation. Unlike the BA-BA domain wall, which exhibits gapped insulating behaviour, the domain walls surrounded by different stackings of bilayer regions feature valley-polarized edge states. Our findings are consistent with other theoretical calculations, such as from the tight-binding model and first-principles calculations, and agree with experimental observations. Published by AIP Publishing.

Abstract: Two-dimensional polyaniline with structural unit C3N is an indirect semiconductor with 0.4 eV band gap, which has attracted a lot of interest because of its unusual electronic, optoelectronic, thermal, and mechanical properties useful for various applications. Adsorption of adatoms is an effective method to improve and tune the properties of C3N. Using first-principles calculations, we investigated the adsorption of adatoms, including H, O, S, F, Cl, B, C, Si, N, P, Al, Li, Na, K, Be, Mg, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, and Zn, on C3N. Depending on the adatom size and the number of valence electrons, they may induce metallic, half-metallic, semiconducting, and ferromagnetic-metallic behavior. In addition, we investigate the effects of an electrical field, charging, and strain on C3N and found how the electronic and magnetic properties are modified. Semi- and full hydrogenation are studied. From the mechanical and thermal stability of C3N monolayer, we found it to be a hard material that can withstand large strain. From our calculations, we gained novel insights into the properties of C3N demonstrating its unique electronic and magnetic properties that can be useful for semiconducting, nanosensor, and catalytic applications.

Abstract: An effective tight-binding (TB) Hamiltonian for monolayer GeSnH2 is constructed which has an inversion-asymmetric honeycomb structure. The low-energy band structure of our TB model agrees very well with previous ab initio calculations even under biaxial tensile strain. Our model predicts a phase transition at 7.5% biaxial tensile strain in agreement with DFT calculations. Upon 8.5% strain the system exhibits a band gap of 134 meV, suitable for room temperature applications. It is shown that an external applied magnetic field produces a special phase which is a combination of the quantum Hall (QH) and quantum spin Hall (QSH) phases; and at a critical magnetic field strength the QSH phase completely disappears. The topological nature of the phase transition is confirmed from: (1) the calculation of the Z(2) topological invariant, and (2) quantum transport properties of disordered GeSnH2 nanoribbons which allows us to determine the universality class of the conductance fluctuations. The application of an external applied magnetic field reduces the conductance fluctuations by a factor of root 2.

Abstract: The electron states in graphene-based magnetic dot and magnetic ring structures and combinations of both are investigated. The corresponding spectra are studied as a function of the radii, the strengths of the inhomogeneous magnetic field and of a uniform background field, the strength of an electrostatic barrier and the angular momentum quantum number. In the absence of an external magnetic field we have only long-lived quasi-bound and scattering states and we assess their influence on the density of states. In addition, we consider elastic electron scattering by a magnetic dot, whose average B vanishes, and show that the Hall and longitudinal resistivities, as a function of the Fermi energy, exhibit a pronounced oscillatory structure due to the presence of quasi-bound states. Depending on the dot parameters this oscillatory structure differs substantially for energies below and above the first Landau level.

Abstract: The scattering of two-dimensional (2D) massless Dirac electrons is investigated in the presence of a random array of circular mass barriers. The inverse momentum relaxation time and the Hall factor are calculated and used to obtain parallel and perpendicular resistivity components within linear transport theory. We found a nonzero perpendicular resistivity component which has opposite sign for electrons in the different K and K′ valleys. This property can be used for valley filter purposes. The total cross section for scattering on penetrable barriers exhibits resonances due to the presence of quasibound states in the barriers that show up as sharp gaps in the cross section while for Schrödinger electrons they appear as peaks.

Abstract: The conductance, the transmission, and the reflection probabilities through rectangular potential barriers and p-n junctions are obtained for bilayer graphene taking into account the four bands of the energy spectrum. We have evaluated the importance of the skew hopping parameters gamma(3) and gamma(4) to these properties and show that for energies E > gamma(1)/100 their effect is negligible. For high energies two modes of propagation exist and we investigate scattering between these modes. For perpendicular incidence both propagation modes are decoupled, and scattering between them is forbidden. This extends the concept of pseudospin as defined within the two-band approximation to a four-band model and corresponds to the (anti) symmetry of the wave functions under in-plane mirroring. New transmission resonances are found that appear as sharp peaks in the conductance which are absent in the two-band approximation. The application of an interlayer bias to the system (1) breaks the pseudospin structure, (2) opens a band gap that results in a distinct feature of suppressed transmission in the conductance, and (3) breaks the angular symmetry with respect to normal incidence in the transmission and reflection.

Abstract: The scattering of two-dimensional (2D) massless electrons as presented in graphene in the presence of a random array of circular magnetic flux tubes is investigated. The momentum relaxation time and the Hall factor are obtained using optical theorem techniques for scattering. Electrons with energy close to those of the Landau levels of the flux tubes exhibit resonant scattering and have a long life-time to reside inside the magnetic flux tube. These resonances appear as sharp structures in the Hall factor and the magneto-resistance.

Abstract: Superconducting correlations in an isolated metallic grain are governed by the interplay between two energy scales: the mean level spacing delta and the bulk pairing gap Delta(0), which are strongly influenced by the position of the chemical potential with respect to the closest single-electron level. In turn superconducting correlations affect the position of the chemical potential. Within the parity projected BCS model we investigate the probability distribution of the chemical potential in a superconducting grain with randomly distributed single-electron levels. Taking into account statistical fluctuations of the chemical potential due to the pairing interaction, we find that such fluctuations have a significant impact on the critical level spacing delta(c) at which the superconducting correlations cease: the critical ratio delta(c)/Delta(0) at which superconductivity disappears is found to be increased.

Abstract: Themean-field Hubbard model is used to investigate the formation of the antiferromagnetic phase in hexagonal graphene rings with inner zigzag edges. The outer edge of the ring was taken to be either zigzag or armchair, and we found that both types of structures can have a larger antiferromagnetic interaction as compared with hexagonal dots. This difference could be partially ascribed to the larger number of zigzag edges per unit area in rings than in dots. Furthermore, edge states localized on the inner ring edge are found to hybridize differently than the edge states of dots, which results in important differences in the magnetism of graphene rings and dots. The largest staggered magnetization is found when the outer edge has a zigzag shape. However, narrow rings with armchair outer edge are found to have larger staggered magnetization than zigzag hexagons. The edge defects are shown to have the least effect on magnetization when the outer ring edge is armchair shaped. DOI: 10.1103/PhysRevB.87.085434

Abstract: We investigate the energy levels and optical properties of a circular graphene quantum dot in the presence of an external magnetic field perpendicular to the dot. Based on the Dirac-Weyl equation and assuming zero outward current at the edge of the dot we present the results for two different types of boundary conditions, i.e. infinite-mass (IMBC) and zigzag boundary conditions. We found that the dot with zigzag edges displays a zero-energy state in the energy spectra while this is not the case for the IMBCs. For both boundary conditions, the confinement becomes dominated by the magnetic field, where the energy levels converge to the Landau levels as the magnetic field increases. The effect of boundary conditions on the electron-and hole-energy states is found to affect the interband absorption spectra, where we found larger absorption in the case of IMBCs. The selection rules for interband optical transitions are determined and discussed for both boundary conditions.

Abstract: The electronic properties of a triaxially strained hexagonal graphene flake with either armchair or zigzag edges are investigated using molecular dynamics simulations and tight-binding calculations. We found that (i) the pseudomagnetic field in strained graphene flakes is not uniform neither in the center nor at the edge of zigzag terminated flakes, (ii) the pseudomagnetic field is almost zero in the center of armchair terminated flakes but increases dramatically near the edges, (iii) the pseudomagnetic field increases linearly with strain, for strains lower than 15% but increases nonlinearly beyond it, (iv) the local density of states in the center of the zigzag hexagon exhibits pseudo-Landau levels with broken sublattice symmetry in the zeroth pseudo-Landau level, and in addition there is a shift in the Dirac cone due to strain induced scalar potentials, and (v) there is size effect in pseudomagnetic field. This study provides a realistic model of the electronic properties of inhomogeneously strained graphene where the relaxation of the atomic positions is correctly included together with strain induced modifications of the hopping terms up to next-nearest neighbors.

Abstract: Density-functional tight-binding and classical molecular dynamics simulations are used to investigate the structural deformations and melting of planar carbon nanoclusters C-N with N = 2-55. The minimum-energy configurations for different clusters are used as starting configurations for the study of the temperature effects on the bond breaking and rotation in carbon lines (N < 6), carbon rings (5 < N < 19), and graphene nanoflakes. The larger the rings (graphene nanoflakes) the higher the transition temperature (melting point) with ring-to-line (perfect-to-defective) transition structures. The melting point was obtained by using the bond energy, the Lindemann criteria, and the specific heat. We found that hydrogen-passivated graphene nanoflakes (CNHM) have a larger melting temperature with a much smaller dependence on size. The edges in the graphene nanoflakes exhibit several different metastable configurations (isomers) during heating before melting occurs. DOI: 10.1103/PhysRevB.87.134103