Abstract: Bogoliubov-de Gennes theory is used to investigate the effect of the size of a superconducting square on the vortex states in the quantum confinement regime. When the superconducting coherence length is comparable to the Fermi wavelength, the shape resonances of the superconducting order parameter have strong influence on the vortex configuration. Several unconventional vortex states, including asymmetric ones, giant-multivortex combinations, and states comprising giant antivortices, were found as ground states and their stability was found to be very sensitive on the value of k(F)xi(0), the size of the sample W, and the magnetic flux Phi. By increasing the temperature and/or enlarging the size of the sample, quantum confinement is suppressed and the conventional mesoscopic vortex states as predicted by the Ginzburg-Laudau (GL) theory are recovered. However, contrary to the GL results we found that the states containing symmetry-induced vortex-antivortex pairs are stable over the whole temperature range. It turns out that the inhomogeneous order parameter induced by quantum confinement favors vortex-antivortex molecules, as well as giant vortices with a rich structure in the vortex core-unattainable in the GL domain.

Abstract: Using an effective one-dimensional cluster expansion in combination with first-principles electronic structure calculations we have studied the energetics and electronic properties of Sb-Te layered systems. For a Te concentration between 0 and 60 at. % an almost continuous series of metastable structures is obtained consisting of consecutive Sb bilayers next to consecutive Sb2Te3 units, with the general formula (Sb-2)(n)(Sb2Te3)(m) (n, m = 1,2, ... ). Between 60 and 100 at.% no stable structures are found. We account explicitly for the weak van derWaals bonding between Sb bilayers and Sb2Te3 units by using a recently developed functional, which strongly improves the interlayer bonding distances. At T = 0 K, no evidence is found for the existence of two separate single-phase regions delta and gamma and a two-phase region delta + gamma. Metastable compounds with a Te concentration between 0 and 40 at. % are semimetallic, whereas compounds with a Te concentration between 50 and 60 at. % are semiconducting. Compounds with an odd number of Sb layers are metallic and have a much higher formation energy than those with an even number of consecutive Sb layers, thereby favoring the formation of Sb bilayers.

Abstract: Ballistic transport through nanoscale devices with time-dependent Rashba-type spin- orbit interaction (SOI) can lead to spin-polarized wavepackets that appear even for completely unpolarized input. The SOI that oscillates in a finite domain generates density and spin polarization fluctuations that leave the region as propagating waves. In particular, spin polarization has space and time dependence even in regions without SOI. Our results are based on an analytical solution of the time-dependent Schrodinger equation. The relevant Floquet quasi-energies that are obtained appear in the energy spectrum of both the transmitted and the reflected waves.

Abstract: A procedure to derive the Ginzburg-Landau (GL) theory from the multiband BCS Hamiltonian is developed in a general case with an arbitrary number of bands and arbitrary interaction matrix. It combines the standard Gor'kov truncation and a subsequent reconstruction in order to match accuracies of the obtained terms. This reconstruction recovers the phenomenological GL theory as obtained from the Landau model of phase transitions but offers explicit microscopic expressions for the relevant parameters. Detailed calculations are presented for a three-band system treated as a prototype multiband superconductor. It is demonstrated that the symmetry in the coupling matrix may lead to the chiral ground state with the phase frustration, typical for systems with broken time-reversal symmetry. DOI: 10.1103/PhysRevB.87.134510

Abstract: Using highly efficient GPU-based simulations of the tight-binding Bogoliubov-de Gennes equations we solve self-consistently for the pair correlation in rhombohedral (ABC) and Bernal (ABA) multilayer graphene by considering a finite intrinsic s-wave pairing potential. We find that the two different stacking configurations have opposite bulk/surface behavior for the order parameter. Surface superconductivity is robust for ABC stacked multilayer graphene even at very low pairing potentials for which the bulk order parameter vanishes, in agreement with a recent analytical approach. In contrast, for Bernal stacked multilayer graphene, we find that the order parameter is always suppressed at the surface and that there exists a critical value for the pairing potential below which no superconducting order is achieved. We considered different doping scenarios and find that homogeneous doping strongly suppresses surface superconductivity while nonhomogeneous field-induced doping has a much weaker effect on the superconducting order parameter. For multilayer structures with hybrid stacking (ABC and ABA) we find that when the thickness of each region is small (few layers), high-temperature surface superconductivity survives throughout the bulk due to the proximity effect between ABC/ABA interfaces where the order parameter is enhanced. DOI: 10.1103/PhysRevB.87.134509

Abstract: In 1976, Larkin and Ovchinnikov [Zh. Eksp. Teor. Fiz. 68, 1915 (1975) [Sov. Phys.–JETP 41, 960 (1976)]] predicted that vortex matter in superconductors driven by an electrical current can undergo an abrupt dynamic transition from a flux-flow regime to a more dissipative state at sufficiently high vortex velocities. Typically, this transition manifests itself as a large voltage jump at a particular current density, so-called instability current density J∗, which is smaller than the depairing current. By tuning the effective pinning strength in Al films, using an artificial periodic pinning array of triangular holes, we show that a unique and well-defined instability current density exists if the pinning is strong, whereas a series of multiple voltage transitions appear in the relatively weaker pinning regime. This behavior is consistent with time-dependent Ginzburg-Landau simulations, where the multiple-step transition can be unambiguously attributed to the progressive development of vortex chains and subsequently phase-slip lines. In addition, we explore experimentally the magnetic braking effects, caused by a thick Cu layer deposited on top of the superconductor, on the instabilities and the vortex ratchet effect.

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: The coupling of angstrom-sized electron probes with spin-polarized electronic transitions shows that the inelastically scattered probe electron is in a mixed state containing electron vortices with nonzero orbital angular momentum. These electrons create an asymmetric intensity distribution in energy filtered diffraction patterns, giving access to maps of the magnetic moments with atomic resolution. A feasibility experiment shows evidence of the predicted effect. Potential applications are column-by-column maps of magnetic ordering, and the creation of angstrom-sized free electrons with orbital angular momentum by inelastic scattering in a thin ferromagnetic foil.

Abstract: Monte Carlo (MC) methods have a long-standing history as partners of molecular dynamics (MD) to simulate the evolution of materials at the atomic scale. Among these techniques, the uniform-acceptance force-bias Monte Carlo (UFMC) method [ G. Dereli Mol. Simul. 8 351 (1992)] has recently attracted attention [ M. Timonova et al. Phys. Rev. B 81 144107 (2010)] thanks to its apparent capacity of being able to simulate physical processes in a reduced number of iterations compared to classical MD methods. The origin of this efficiency remains, however, unclear. In this work we derive a UFMC method starting from basic thermodynamic principles, which leads to an intuitive and unambiguous formalism. The approach includes a statistically relevant time step per Monte Carlo iteration, showing a significant speed-up compared to MD simulations. This time-stamped force-bias Monte Carlo (tfMC) formalism is tested on both simple one-dimensional and three-dimensional systems. Both test-cases give excellent results in agreement with analytical solutions and literature reports. The inclusion of a time scale, the simplicity of the method, and the enhancement of the time step compared to classical MD methods make this method very appealing for studying the dynamics of many-particle systems.

Abstract: In this work we present an alternative way to look at electron diffraction in a transmission electron microscope.
Instead of writing the scattering amplitude in Fourier space as a set of plane waves,we use the cylindrical Fourier transform to describe the scattering amplitude in a basis of orbital angular momentum (OAM) eigenstates. We show how working in this framework can be very convenient when investigating, e.g., rotation and screw-axis symmetries. For the latter we find selection rules on the OAM coefficients that unambiguously reveal the handedness of the screw axis. Detecting the OAM coefficients of the scattering amplitude thus offers the possibility to detect the handedness of crystals without the need for dynamical simulations, the thickness of the sample, nor the exact crystal structure. We propose an experimental setup to measure the OAM components where an image of the crystal is taken after inserting a spiral phase plate in the diffraction plane and perform multislice simulations on α quartz to demonstrate how the method indeed reveals the chirality. The experimental feasibility of the technique is discussed together with its main advantages with respect to chirality determination of screw axes. The method shows how the use of a spiral phase plate can be extended from a simple phase imaging technique to a tool to measure the local OAM decomposition of an electron wave, widening the field of interest well beyond chiral space group determination.

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