Abstract: Dynamics of magnetic skyrmions in hybrid ferromagnetic films harbors interesting physical phenomena and holds promise for technological applications. In this work, we discuss the behavior of magnetic skyrmions when coupled to superconducting vortices in a ferromagnet-superconductor heterostructure. We use numerical simulations and analytic arguments within London and Thiele formalisms to reveal broader possibilities for manipulating the skyrmion-vortex dynamic correlations in the hybrid system, that are not possible in its separated constituents. We explore the thresholds of particular dynamic phases, and quantify the phase diagram as a function of the relevant material parameters, applied current, and induced magnetic torques. Finally, we demonstrate the broad and precise tunability of the skyrmion Hall angle in the presence of vortices, with respect to currents applied to either or both the superconductor and the ferromagnet within the heterostructure.

Abstract: We analyze the paradigmatic competition between intraband and crossband Cooper-pair formation in twoband superconductors, neglected in most works to date. We derive the phase-sensitive gap equations and describe the crossover between the intraband-dominated and the crossband-dominated regimes, delimited by a “gapless” state. Experimental signatures of crosspairing comprise notable gap splitting in the excitation spectrum, non-BCS behavior of gaps versus temperature, as well as changes in the pairing symmetry as a function of temperature. The consequences of these findings are illustrated on the examples of MgB2 and Ba0.6K0.4Fe2As2.

Abstract: High tunability of two-dimensional magnetic materials (by strain, gating, heterostructuring, or otherwise) provides unique conditions for studying versatile magnetic properties and controlling emergent magnetic phases. Expanding the scope of achievable magnetic phenomena in such materials is important for both fundamental and technological advances. Here we perform atomistic spin-dynamics simulations to explore the (chiral) magnetic phases of atomic monolayers in the limit of suppressed first-neighbors exchange interaction. We report the rich phase diagram of exotic magnetic configurations, obtained for both square and honeycomb lattice symmetries, comprising coexistence of ferromagnetic and antiferromagnetic spin cycloids, as well as multiple types of magnetic skyrmions. We perform a minimum-energy path analysis for the skyrmion collapse to evaluate the stability of such topological objects and reveal that magnetic monolayers could be good candidates to host the antiferromagnetic skyrmions that are experimentally evasive to date.

Abstract: The optoelectronic properties of confined water form one of the most active research areas in the past few years. Here we present the multiscale methodology to discern the out-of-plane electronic and dipolar dielectric constants (epsilon(el)(perpendicular to) and epsilon(diP)(perpendicular to)) of strongly confined water. We reveal that epsilon(perpendicular to el) and epsilon(diP)(perpendicular to) become comparable for water confined in angstrom-scale channels (with a height of less than 15 angstrom) within graphene (GE) and hexagonal boron nitride (hBN) bilayers. Channel height (h) associated with a minimum in both epsilon(e)(l)(perpendicular to) and epsilon(dip)(perpendicular to) is linked to the formation of the ordered structure of ice for h approximate to (7 -7.5) angstrom. The recently measured total dielectric constant epsilon(T)(perpendicular to) of nanoconfined water [L. Fumagalli et al., Science 360, 1339 (2018)] is corroborated by our results. Furthermore, we evaluate the contribution from the encapsulating membranes to the dielectric properties, as a function of the interlayer spacing, i.e., the height of the confining channel for water. Finally, we conduct analysis of the optical properties of both confined water and GE membranes, and show that the electron energy loss function of confined water strongly differs from that of bulk water.

Abstract: After the discovery of magnetism in monolayer CrI3, the magnetic properties of different 2D materials from the chromium-trihalide family are intuitively assumed to be similar, yielding magnetic anisotropy from the spin-orbit coupling on halide ligands. Here we reveal significant differences between the CrI3 and CrBr3 magnetic monolayers in their magnetic anisotropy, resulting Curie temperature, hysteresis in external magnetic field, and evolution of magnetism with strain, all predominantly attributed to distinctly different interplay of atomic contributions to spin-orbit coupling in two materials.

Abstract: Considering a noncentrosymmetric pinning texture composed of a square array of triangular holes, the magnetic flux penetration and expulsion are investigated experimentally and theoretically. A direct visualization of the magnetic landscape obtained using a magneto-optical technique on a Nb film is complemented by a multiscale numerical modeling. This combined approach allows the magnetic flux dynamics to be identified from the single flux quantum limit up to the macroscopic electromagnetic response. Within the theoretical framework provided by time-dependent Ginzburg-Landau simulations, an estimation of the in-plane current anisotropy is obtained and its dependence with the radius of the curvature of hole vertices is addressed. These simulations show that current crowding plays an important role in channeling the flux motion, favoring hole-to-hole flux hopping rather than promoting interstitial flux displacement in between the holes. The resulting anisotropy of the critical current density gives rise to a distinct pattern of discontinuity lines for increasing and decreasing applied magnetic fields, in sharp contrast to the invariable patterns reported for centrosymmetric pinning potentials. This observation is partially accounted for by the rectification effect, as demonstrated by finite-element modeling. At low temperatures, where magnetic field penetration is dominated by thermomagnetic instabilities, highly directional magnetic flux avalanches with a fingerlike shape are observed to propagate along the easy axis of the pinning potential. This morphology is reproduced by numerical simulations. Our findings demonstrate that anisotropic pinning landscapes and, in particular, ratchet potentials produce subtle modifications to the critical state field profile that are reflected in the distribution of discontinuity lines.

Abstract: In a thin-film topological insulator (TI), the edge states on two surfaces may couple by quantum tunneling, opening a gap known as the hybridization gap. Controlling the hybridization gap and transport has a variety of potential uses in photodetection and energy-harvesting applications. In this paper, we report the effect of strain, and electric and magnetic field, on the hybridization gap and transport in a thin Bi2Se3 film, investigated within the tight-binding theoretical framework. We demonstrate that vertical compression decreases the hybridization gap, as does tensile in-plane strain. Applying an electric field breaks the inversion symmetry and leads to a Rashba-like spin splitting proportional to the electric field, hence closing and reopening the gap. The influence of a magnetic field on thin-film TI is also discussed, starting from the role of an out-of-plane magnetic field on quantum Hall states. We further demonstrate that the hybridization gap can be controlled by an in-plane magnetic field, and that by applying a sufficiently strong field a quantum phase transition from an insulator to a semimetal can be achieved.

Abstract: Using a particle-based model we examine the depinning motion of solitons in skyrmion chains in quasi -onedimensional (1D) and two-dimensional (2D) systems containing embedded 1D interfaces. The solitons take the form of a particle or hole in a commensurate chain of skyrmions. Under an applied drive, just above a critical depinning threshold, the soliton moves with a skyrmion Hall angle of zero. For higher drives, the entire chain depins, and in a 2D system we observe that both the solitons and chain move at zero skyrmion Hall angle and then transition to a finite skyrmion Hall angle as the drive increases. In a 2D system with a 1D interface that is at an angle to the driving direction, there can be a reversal of the sign of the skyrmion Hall angle from positive to negative. Our results suggest that solitons in skyrmion systems could be used as information carriers in racetrack geometries that would avoid the drawbacks of finite skyrmion Hall angles. The soliton states become mobile at significantly lower drives than the depinning transition of the skyrmion chains themselves.

Abstract: We fabricated an asymmetric nanoscale SQUID consisting of one nanobridge weak link and one Dayem bridge weak link. The current phase relation of these particular weak links is characterized by multivaluedness and linearity. While the latter is responsible for a particular magnetic field dependence of the critical current (so-called vorticity diamonds), the former enables the possibility of different vorticity states (phase winding numbers) existing at one magnetic field value. In experiments the observed critical current value is stochastic in nature, does not necessarily coincide with the current associated with the lowest energy state and critically depends on the measurement conditions. In this paper, we unravel the origin of the observed metastability as a result of the phase dynamics happening during the freezing process and while sweeping the current. Moreover, we employ special measurement protocols to prepare the desired vorticity state and identify the (hidden) phase slip dynamics ruling the detected state of these nanodevices. In order to gain insights into the dynamics of the condensate and, more specifically the hidden phase slips, we performed time-dependent Ginzburg-Landau simulations.

Abstract: A wave-packet time evolution method, based on the split-operator technique, is developed to investigate the scattering of quasiparticles at a normal-superconductor interface of arbitrary profile and shape. As a practical application, we consider a system where low-energy electrons can be described as Dirac particles, which is the case for most two-dimensional materials, such as graphene and transition-metal dichalcogenides. However, the method is easily adapted for other cases such as electrons in few-layer black phosphorus or any Schrodinger quasiparticles within the effective mass approximation in semiconductors. We employ the method to revisit Andreev reflection in mono-, bi-, and trilayer graphene, where specular-and retro-reflection cases are observed for electrons scattered by a steplike superconducting region. The effect of opening a zero-gap channel across the superconducting region on the electron and hole scattering is also addressed, as an example of the versatility of the technique proposed here.

Abstract: We theoretically investigate the binding energy and electron-hole (e-h) overlap of excitonic states confined at the interface between two-dimensional materials with type-II band alignment, i.e., with lowest conduction and highest valence band edges placed in different materials, arranged in a side-by-side planar heterostructure. We propose a variational procedure within the effective mass approximation to calculate the exciton ground state and apply our model to a monolayer MoS2/WS2 heterostructure. The role of nonabrupt interfaces between the materials is accounted for in our model by assuming a WxMo1-xS2 alloy around the interfacial region. Our results demonstrate that (i) interface-bound excitons are energetically favorable only for small interface thickness and/or for systems under high dielectric screening by the materials surrounding the monolayer, and that (ii) the interface exciton binding energy and its e-h overlap are controllable by the interface width and dielectric environment.

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: Because of strong flux confinement in mesoscopic superconductors, a giant vortex may appear in the ground state of the system in an applied magnetic field. This multiquanta vortex can then split into individual vortices (and vice versa) as a function of, e.g., applied current, magnetic field, or temperature. Here we show that such transitions can be identified by calorimetry, as the formation or splitting of a giant vortex results in a clear jump in measured heat capacity versus external drive. We attribute this phenomenon to an abrupt change in the density of states of the quasiparticle excitations in the vortex core(s), and further link it to a sharp change of the magnetic susceptibility at the transitionproving that the formation of a giant vortex can also be detected by magnetometry.

Abstract: We show in the case of a superconducting Nb ladder that a mesoscopic superconductor typically exhibits magnetoresistance oscillations whose amplitude and temperature dependence are different from those stemming from the Little-Parks effect. We demonstrate that these large resistance oscillations (as well as the monotonic background on which they are superimposed) are due to current-excited moving vortices, where the applied current in competition with the oscillating Meissner currents imposes or removes the barriers for vortex motion in an increasing magnetic field. Because of the ever present current in transport measurements, this effect should be considered in parallel with the Little-Parks effect in low-critical temperature (T-c) samples, as well as with recently proposed thermal activation of dissipative vortex-antivortex pairs in high-T-c samples.

Abstract: We describe investigations of the largely unexplored field of mesoscopic type-I superconductors. Micromagnetometry and 3D Ginzburg-Landau simulations of our single crystal β-tin samples in this regime reveal size- and temperature-dependent supercritical fields whose behavior is radically different from the bulk critical field HcB. We find that complete suppression of the intermediate state in medium-size samples can result in a surprising reduction of the critical field significantly below HcB. We also reveal an evolution of the superconducting-to-normal phase transition from the expected irreversible first order at low temperatures through the previously unobserved reversible first-order to a second-order transition close to Tc, where the critical field can be many times larger than HcB. Finally, we have identified striking correlations between the mesoscopic Hc3 for nucleation of surface superconductivity and the thermodynamic Hc near Tc. All these observations are entirely unexpected in the conventional type-I picture.