Abstract: We analyze the combined effect of three ingredients of an aggregation model-surface tension, particle flow and particle source-representing typical characteristics of many aggregation growth processes in nature. Through extensive numerical experiments and for different underlying lattice structures we demonstrate that the location of incoming particles and their preferential direction of flow can significantly affect the resulting general shape of the aggregate, while the surface tension controls the surface roughness. Combining all three ingredients increases the aggregate shape plasticity, yielding a wider spectrum of shapes as compared to earlier works that analyzed these ingredients separately. Our results indicate that the considered combination of effects is fundamental for modeling the polymorphic growth of a wide variety of structures in confined geometries and/or in the presence of external fields, such as rocks, crystals, corals, and biominerals.

Abstract: The autoregressive neural networks are emerging as a powerful computational tool to solve relevant problems in classical and quantum mechanics. One of their appealing functionalities is that, after they have learned a probability distribution from a dataset, they allow exact and efficient sampling of typical system configurations. Here we employ a neural autoregressive distribution estimator (NADE) to boost Markov chain Monte Carlo (MCMC) simulations of a paradigmatic classical model of spin-glass theory, namely, the two-dimensional Edwards-Anderson Hamiltonian. We show that a NADE can be trained to accurately mimic the Boltzmann distribution using unsupervised learning from system configurations generated using standard MCMC algorithms. The trained NADE is then employed as smart proposal distribution for the Metropolis-Hastings algorithm. This allows us to perform efficient MCMC simulations, which provide unbiased results even if the expectation value corresponding to the probability distribution learned by the NADE is not exact. Notably, we implement a sequential tempering procedure, whereby a NADE trained at a higher temperature is iteratively employed as proposal distribution in a MCMC simulation run at a slightly lower temperature. This allows one to efficiently simulate the spin-glass model even in the low-temperature regime, avoiding the divergent correlation times that plague MCMC simulations driven by local-update algorithms. Furthermore, we show that the NADE-driven simulations quickly sample ground-state configurations, paving the way to their future utilization to tackle binary optimization problems.

Abstract: We describe a semianalytic approach to the two-band Ginzburg-Landau theory, which predicts the behavior of vortices in two-band superconductors. We show that the character of the short-range vortex-vortex interaction is determined by the sign of the normal domain-superconductor interface energy, in analogy with the conventional differentiation between type I and type II superconductors. However, we also show that the long-range interaction is determined by a modified Ginzburg-Landau parameter κ*, different from the standard κ of a bulk superconductor. This opens the possibility for nonmonotonic vortex-vortex interaction, which is temperature dependent, and can be further tuned by alterations of the material on the microscopic scale.

Abstract: Using the Ginzburg-Landau theory extended to the next-to-leading order, we determine numerically the healing lengths of the two order parameters at the two-gap superconductor/normal metal interface. We demonstrate on several examples that those can be different even in the strict domain of applicability of the Ginzburg-Landau theory. This justifies the use of this theory to describe relevant physics of two-gap superconductors, distinguishing them from their single-gap counterparts. The calculational degree of complexity increases only slightly with respect to the conventional Ginzburg-Landau expansion, thus the extended Ginzburg-Landau model remains numerically far less demanding compared to the full microscopic approaches.

Abstract: Based on the Gor'kov formalism for a clean s-wave superconductor, we develop an extended version of the single-band Ginzburg-Landau (GL) theory by means of a systematic expansion in the deviation from the critical temperature T(c), i.e., tau = 1 – T/T(c). We calculate different contributions to the order parameter and the magnetic field: the leading contributions (proportional to tau(1/2) in the order parameter and. t in the magnetic field) are controlled by the standard GL theory, while the next-to-leading terms (proportional to tau(3/2) in the gap and proportional to tau(2) in the magnetic field) constitute the extended GL (EGL) approach. We derive the free-energy functional for the extended formalism and the corresponding expression for the current density. To illustrate the usefulness of our formalism, we calculate, in a semianalytical form, the temperature-dependent correction to the GL parameter at which the surface energy becomes zero, and analytically, the temperature dependence of the thermodynamic critical field. We demonstrate that the EGL formalism is not just a mathematical extension to the theory: variations of both the gap and the thermodynamic critical field with temperature calculated within the EGL theory are found in very good agreement with the full BCS results down to low temperatures, which dramatically improves the applicability of the formalism compared to its standard predecessor.

Abstract: We study nonlinear flux dynamics in a current-carrying type-I superconductor. The stray magnetic field of the current induces the intermediate state, where nucleation of flux domains is discretized to a single fluxoid at a time, while their final shape (tubular or laminar), size, and nucleation rate depend on applied current and edge conditions. The current induces opposite flux domains on opposite sides of the sample, and subsequently drives them to annihilation-which is also discretized, as a sequence of vortex-antivortex pairs. The discretization of both nucleation and annihilation leaves measurable traces in the voltage across the sample and in locally probed magnetization. The reported dynamic phenomena thus provide an unambiguous proof of a flux quantum being the smallest building block of the intermediate state in type-I superconductors.

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: Magnetic skyrmions are topological spin configurations in materials with chiral Dzyaloshinskii-Moriya interaction (DMI), that are potentially useful for storing or processing information. To date, DMI has been found in few bulk materials, but can also be induced in atomically thin magnetic films in contact with surfaces with large spin-orbit interactions. Recent experiments have reported that isolated magnetic skyrmions can be stabilized even near room temperature in few-atom-thick magnetic layers sandwiched between materials that provide asymmetric spin-orbit coupling. Here we present the minimum-energy path analysis of three distinct mechanisms for the skyrmion collapse, based on ab initio input and the performed atomic-spin simulations. We focus on the stability of a skyrmion in three atomic layers of Co, either epitaxial on the Pt(111) surface or within a hybrid multilayer where DMI nontrivially varies per monolayer due to competition between different symmetry breaking from two sides of the Co film. In laterally finite systems, their constrained geometry causes poor thermal stability of the skyrmion toward collapse at the boundary, which we show to be resolved by designing the high-DMI structure within an extended film with lower or no DMI.

Abstract: The magnetic flux dependence of energy and supercurrent in mesoscopic d-wave symmetric and asymmetric superconducting loops is investigated by numerically solving the Bogoliubov-de Gennes equations self-consistently. For square loops, we find an hc/e-flux periodicity in energy and supercurrent and demonstrate that the flux periodicity is sensitive to the hole size and the superconducting pairing strength as well as temperature. The hc/2e-periodic behavior can be restored almost entirely when we displace the central hole sufficiently out of the center of the sample. In rectangular loops, the discrete current-carrying low-energy spectrum can exist for an odd winding number of the order parameter.

Abstract: We have used differential magneto-optical (MO) imaging to investigate the mixed state of superconducting Bi2Sr2CaCu2O8+ (BSCCO) microdisks fabricated on a single-crystal sample. MO difference images of the stray field distribution over a range of out-of-plane fields allow us to distinguish between flux that is penetrating the disks and that entering the underlying BSCCO platelet. We find that flux preferentially flows along linear defects into the interstitial platelet regions up to a characteristic field Hp, above which flux enters the disks. We identify this as the field of first penetration of pancake vortices over the Bean-Livingston barrier around the disks, where Hp(T) at intermediate temperatures is well described by an exponentially decaying function with a characteristic temperature T0=19 K. At a given temperature, a minority of the disks exhibit a lower penetration field and we correlate the location of these disks with the linear defects in the BSCCO crystal.

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: Our calculations, within known limitations of Ginzburg-Landau theory, are fully correct and valid for transport phenomena in asymmetric mesoscopic superconductors, deep in the superconducting state. We deemed the experiments of Burlakov et al. [JETP Lett. 86, 517 (2007)] relevant and important to mention in the general context of our paper since the observed shifts in the oscillations of different quantities are qualitatively similar, even though those measurements are performed close to the superconducting-normal state transition in the so-called Little-Parks regime.

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.