Abstract: We numerically investigate the spin-orbit (SO) coupling effect on the magnetic flux evolution of energy and supercurrent in mesoscopic d-wave superconducting loops by solving the spin-generalized Bogoliubov-de Gennes equations self-consistently. It is found that the energy spectrum splits when the SO interaction is involved and the Majorana zero mode can be realized in the [100] edges of square systems for an appropriate SO coupling strength. Superconducting phase transitions appear when the energy gap closes, accompanied by energy jumps between different energy parabolas in the ground state, which provides a possible mechanism to support fractional flux periodicity of supercurrent. Moreover, in the case of rectangular loops with SO coupling, the jumps of the ground-state energy gradually disappear by increasing the ratio of length to height of the sample, and a paramagnetic response with opposite direction of the screening current around zero flux value can occur in such systems.

Abstract: Bi2Se3 is a three-dimensional topological insulator which has been extensively studied because it has a single Dirac cone on the surface, inside a relatively large bulk band gap. However, the effect of two-dimensional topological insulator Bi bilayers on the properties of Bi2Se3 and vice versa, has not been explored much. Bi bilayers are often present between the quintuple layers of Bi2Se3, since (Bi2)n(Bi2Se3)m form stable ground-state structures. Moreover, Bi2Se3 is a good substrate for growing ultrathin Bi bilayers. By first-principles techniques, we first show that there is no preferable surface termination by either Bi or Se. Next, we investigate the electronic structure of Bi bilayers on top of, or inside a Bi2Se3 slab. If the Bi bilayers are on top, we observe a charge transfer to the quintuple layers that increases the binding energy of the surface Dirac cones. The extra states, originating from the Bi bilayers, were declared to form a topological Dirac cone, but here we show that these are ordinary Rashba-split states. This result, together with the appearance of a new Dirac cone that is localized slightly deeper, might necessitate the reinterpretation of several experimental results. When the Bi bilayers are located inside the Bi2Se3 slab, they tend to split the slab into two topological insulators with clear surface states. Interface states can also be observed, but an energy gap persists because of strong coupling between the neighboring quintuple layers and the Bi bilayers.

Abstract: We revisit the problem of the dynamic response of a superconducting bridge after abruptly switching on the supercritical current. In contrast to previous theoretical works we take into account spatial gradients and use both the local temperature approach and the kinetic equation for the distribution function of quasiparticles. We find that the temperature dependence of the finite delay time t(d) in the voltage response is model dependent and relatively large t(d) is connected with temporary cooling of quasiparticles during decay of superconducting order parameter vertical bar Delta vertical bar in time. It turns out that the presence of even small inhomogeneities in the bridge or finite length of the homogenous bridge favors a local suppression of vertical bar Delta vertical bar during the dynamic response. It results in a decrease of the delay time, in comparison with the spatially uniform model, due to the diffusion of nonequilibrium quasiparticles from the region with locally suppressed vertical bar Delta vertical bar. In the case when the current density is maximal near the edge of a not very wide bridge the delay time is mainly connected with the time needed for the nucleation (entrance) of the first vortex and t(d) could be tuned by a weak external magnetic field. We also find that a short alternating current pulse (sinusoidlike) with zero time average may result in a nonzero time- averaged voltage response where its sign depends on the phase of the ac current.

Abstract: Using advanced three-dimensional simulations, we show that an Abrikosov vortex, trapped inside a cavity perpendicular to an artificial Josephson junction, can serve as a very efficient source for generation of Josephson vortex-antivortex pairs in the presence of the applied electric current. In such a case, the nucleation rate of the pairs can be tuned in a broad range by an out-of-plane ac magnetic field in a broad range of frequencies. This parametrically amplified vortex-antivortex nucleation can be considered as a macroscopic analog of the dynamic Casimir effect, where fluxon pairs mimic the photons and the ac magnetic field plays the role of the oscillating mirrors. The emerging vortex pairs in our system can be detected by the pronounced features in the measured voltage characteristics, or through the emitted electromagnetic radiation, and exhibit resonant dynamics with respect to the frequency of the applied magnetic field. Reported tunability of the Josephson oscillations can be useful for developing high-frequency emission devices.

Abstract: It is well-known that the Ginzburg-Landau (GL) theory is a reliable and powerful theoretical tool to investigate the magnetic response of a superconducting state. However, in its standard form, this approach is not applicable to atomically uniform nano-thin superconducting films which are effective multiband superconductors. Here we discuss a relevant generalization of the GL theory, focusing on the underlying intraband-pairing approximation. (C) 2014 Elsevier B.V. All rights reserved.

Abstract: Phase-slip centers/lines and hot spots are the main mechanisms for dissipation in current-carrying superconducting thin films. The pulsed-current method has recently been shown to be an effective tool in studying the dynamics of phase-slip centers and their evolution to hot spots. We use the time-dependent Ginzburg-Landau theory in the study of the dynamics of the superconducting condensate in superconducting strips under external current and zero external magnetic field. We show that both the flux-flow state (i.e., slow-moving vortices) and the phase-slip line state (i.e., fast-moving vortices) are dynamically stable dissipative units with temperature smaller than the critical one, whereas hot spots, which are localized normal regions where the local temperature exceeds the critical value, expand in time, resulting ultimately in a complete destruction of the condensate. The response time of the system to abrupt switching on of the overcritical current decreases with increasing both the value of the current (at all temperatures) and temperature (for a given value of the applied current). Our results are in good qualitative agreement with experiments we have conducted on Nb thin strips.

Abstract: We expand on a recent study of a lattice model of interacting particles [Phys. Rev. Lett. 111, 110601 (2013)]. The adsorption isotherm and equilibrium fluctuations in particle number are discussed as a function of the interaction. Their behavior is similar to that of interacting particles in porous materials. Different expressions for the particle jump rates are derived from transition-state theory. Which expression should be used depends on the strength of the interparticle interactions. Analytical expressions for the self-and transport diffusion are derived when correlations, caused by memory effects in the environment, are neglected. The diffusive behavior is studied numerically with kinetic Monte Carlo (kMC) simulations, which reproduces the diffusion including correlations. The effect of correlations is studied by comparing the analytical expressions with the kMC simulations. It is found that the Maxwell-Stefan diffusion can exceed the self-diffusion. To our knowledge, this is the first time this is observed. The diffusive behavior in one-dimensional and higher-dimensional systems is qualitatively the same, with the effect of correlations decreasing for increasing dimension. The length dependence of both the self-and transport diffusion is studied for one-dimensional systems. For long lengths the self-diffusion shows a 1/L dependence. Finally, we discuss when agreement with experiments and simulations can be expected. The assumption that particles in different cavities do not interact is expected to hold quantitatively at low and medium particle concentrations if the particles are not strongly interacting.

Abstract: Carbon-based molecular semiconductors are explored for application in spintronics because their small spinorbit coupling promises long spin lifetimes. We calculate the electronic transport from first principles through spin valves comprising bi-and tri-layers of the fullerene molecules C-60 and C-70, sandwiched between two Fe electrodes. The spin polarization of the current, and the magnetoresistance depend sensitively on the interactions at the interfaces between the molecules and the metal surfaces. They are much less affected by the thickness of the molecular layers. A high current polarization (CP > 90%) and magnetoresistance (MR > 100%) at small bias can be attained using C-70 layers. In contrast, the current polarization and the magnetoresistance at small bias are vanishingly small for C-60 layers. Exploiting a generalized Julliere model we can trace the differences in spin-dependent transport between C-60 and C-70 layers to differences between the molecule-metal interface states. These states also allow one to interpret the current polarization and the magnetoresistance as a function of the applied bias voltage.

Abstract: The electric-field response of a one-dimensional ring of interacting fermions, where the interactions are described by the extended Hubbard model, is investigated. By using an accurate real-time propagation scheme based on the Chebyshev expansion of the evolution operator, we uncover various nonlinear regimes for a range of interaction parameters that allows modeling of metallic and insulating (either charge density wave or spin density wave insulators) rings. The metallic regime appears at the phase boundary between the two insulating phases and provides the opportunity to describe either weakly or strongly correlated metals. We find that the fidelity susceptibility of the ground state as a function of magnetic flux piercing the ring provides a very good measure of the short-time response. Even completely different interacting regimes behave in a similar manner at short time scales as long as the ground-state fidelity susceptibility is the same. Depending on the strength of the electric field we find various types of responses: persistent currents in the insulating phase, a dissipative regime, or damped Bloch-like oscillations with varying frequencies or even irregular in nature. Furthermore, we also consider the dimerization of the ring and describe the response of a correlated band insulator. In this case the distribution of the energy levels is more clustered and the Bloch-like oscillations become even more irregular.

Abstract: The influence of silicene's strong spin-orbit interaction and of an external electric field E-z on the transport coefficients are investigated in the presence of a perpendicular magnetic field B. For finite E-z the spin and valley degeneracy of the Landau levels is lifted and leads to additional plateaus in the Hall conductivity, at half-integer values of 4e(2)/h, due to spin intra-Landau-level transitions that are absent in graphene. These plateaus are more sensitive to disorder and thermal broadening than the main plateaus, occurring at integral values of 4e(2)/h, when the Fermi level passes through the Landau levels. We also evaluate the Hall and longitudinal resistivities and critically contrast the results with those for graphene on a substrate.

Abstract: Molecular-dynamic simulations were performed in order to investigate the melting processes of isotropically confined binary systems. We considered two species of particles, which differ by their amount of electric charge. A Lindemann type of criterion was used to determine the angular melting temperature. We demonstrate that the magic-to-normal cluster transition can evolve in two distinct ways, that is, through a structural phase transition of the first order or via a smooth transition where an increase of the shells' width leads to a continuous decreasing mechanical stability of the system. Moreover, for large systems, we demonstrate that the internal cluster exerts a minor effect on the mechanical stability of the external shell. Furthermore, we show that highly symmetric configurations, such as those found for multiple ring structures, have large mechanical stability, i.e., high angular melting temperature.

Abstract: Ripples in pristine freestanding graphene naturally orient themselves in an array that is alternately curved-up and curved-down; maintaining an average height of zero. Using scanning tunneling microscopy (STM) to apply a local force, the graphene sheet will reversibly rise and fall in height until the height reaches 60%-70% of its maximum at which point a sudden, permanent jump occurs. We successfully model the ripples as a spin-half Ising magnetic system, where the height of the graphene plays the role of the spin. The permanent jump in height, controlled by the tunneling current, is found to be equivalent to an antiferromagnetic-to-ferromagnetic phase transition. The thermal load underneath the STM tip alters the local tension and is identified as the responsible mechanism for the phase transition. Four universal critical exponents are measured from our STM data, and the model provides insight into the statistical role of graphene's unusual negative thermal expansion coefficient.