Interface thermal conductivities induced by van der Waals interactions.

The interface heat transfer of two layers induced by van der Waals (vdW) contacts is theoretically investigated, based on first-principles calculations at low temperatures. The results suggest that out-of-plane acoustic phonons with low frequencies dominate the interface thermal transport due to the vdW interaction. The interface thermal conductivity is proportional to the cubic of temperature at very low temperatures, but becomes linearly proportional to temperature as temperature increases. We show that manipulating the strain alters vdW coupling, leading to increased interfacial thermal conductivity at the interface. Our findings provide valuable insights into the interface heat transport in vdW heterostructures and support further design and optimization of electronic and optoelectronic nanodevices based on vdW contacts.

Electronic and valleytronic properties of crystalline boron-arsenide tuned by strain and disorder.

Ab initio density functional theory (DFT) and DFT plus coherent potential approximation (DFT + CPA) are employed to reveal, respectively, the effect of in-plane strain and site-diagonal disorder on the electronic structure of cubic boron arsenide (BAs). It is demonstrated that tensile strain and static diagonal disorder both reduce the semiconducting one-particle band gap of BAs, and a V-shaped p-band electronic state emerges - enabling advanced valleytronics based on strained and disordered semiconducting bulk crystals. At biaxial tensile strains close to 15% the valence band lineshape relevant for optoelectronics is shown to coincide with one reported for GaAs at low energies. The role played by static disorder on the As sites is to promote p-type conductivity in the unstrained BAs bulk crystal, consistent with experimental observations. These findings illuminate the intricate and interdependent changes in crystal structure and lattice disorder on the electronic degrees of freedom of semiconductors and semimetals.

Distinctive g-Factor of Moiré-Confined Excitons in van der Waals Heterostructures.

We investigated the valley Zeeman splitting of excitonic peaks in the microphotoluminescence (µPL) spectra of high-quality hBN/WS2/MoSe2/hBN heterostructures under perpendicular magnetic fields up to 20 T. We identify two neutral exciton peaks in the µPL spectra; the lower-energy peak exhibits a reduced g-factor relative to that of the higher energy peak and much lower than the recently reported values for interlayer excitons in other van der Waals (vdW) heterostructures. We provide evidence that such a discernible g-factor stems from the spatial confinement of the exciton in the potential landscape created by the moiré pattern due to lattice mismatch or interlayer twist in heterobilayers. This renders magneto-µPL an important tool to reach a deeper understanding of the effect of moiré patterns on excitonic confinement in vdW heterostructures.

Zitterbewegung of Moiré Excitons in Twisted MoS_{2}/WSe_{2} Heterobilayers.

The moiré pattern observed in stacked noncommensurate crystal lattices, such as heterobilayers of transition metal dichalcogenides, produces a periodic modulation of their band gap. Excitons subjected to this potential landscape exhibit a band structure that gives rise to a quasiparticle dubbed the moiré exciton. In the case of MoS_{2}/WSe_{2} heterobilayers, the moiré trapping potential has honeycomb symmetry and, consequently, the moiré exciton band structure is the same as that of a Dirac-Weyl fermion, whose mass can be further tuned down to zero with a perpendicularly applied field. Here we show that, analogously to other Dirac-like particles, the moiré exciton exhibits a trembling motion, also known as Zitterbewegung, whose long timescales are compatible with current experimental techniques for exciton dynamics. This promotes the study of the dynamics of moiré excitons in van der Waals heterostructures as an advantageous solid-state platform to probe Zitterbewegung, broadly tunable by gating and interlayer twist angle.

Skyrmion-(Anti)Vortex Coupling in a Chiral Magnet-Superconductor Heterostructure.

We report experimental coupling of chiral magnetism and superconductivity in [IrFeCoPt]/Nb heterostructures. The stray field of skyrmions with radius ≈50 nm is sufficient to nucleate antivortices in a 25 nm Nb film, with unique signatures in the magnetization, critical current, and flux dynamics, corroborated via simulations. We also detect a thermally tunable Rashba-Edelstein exchange coupling in the isolated skyrmion phase. This realization of a strongly interacting skyrmion-(anti)vortex system opens a path toward controllable topological hybrid materials, unattainable to date.

Transition-metal adatoms on 2D-GaAs: a route to chiral magnetic 2D materials by design.

Using relativistic density-functional calculations, we examine the magneto-crystalline anisotropy and exchange properties of transition-metal atoms adsorbed on 2D-GaAs. We show that single Mn and Mo atom (Co and Os) strongly bind on 2D-GaAs, and induce local out-of-plane (in-plane) magnetic anisotropy. When a pair of TM atoms is adsorbed on 2D-GaAs in a close range from each other, magnetisation properties change (become tunable) with respect to concentrations and ordering of the adatoms. In all cases, we reveal presence of strong Dzyaloshinskii-Moriya interaction. These results indicate novel pathways towards two-dimensional chiral magnetic materials by design, tailored for desired applications in magneto-electronics.

Boosting Monte Carlo simulations of spin glasses using autoregressive neural networks.

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.

Two-dimensional hydrogenated buckled gallium arsenide: an ab initio study.

First-principles calculations have been carried out to investigate the stability, structural and electronic properties of two-dimensional (2D) hydrogenated GaAs with three possible geometries: chair, zigzag-line and boat configurations. The effect of van der Waals interactions on 2D H-GaAs systems has also been studied. These configurations were found to be energetic and dynamic stable, as well as having a semiconducting character. Although 2D GaAs adsorbed with H tends to form a zigzag-line configuration, the energy differences between chair, zigzag-line and boat are very small which implies the metastability of the system. Chair and boat configurations display a [Formula: see text]-[Formula: see text] direct bandgap nature, while pristine 2D-GaAs and zigzag-line are indirect semiconductors. The bandgap sizes of all configurations are also hydrogen dependent, and wider than that of pristine 2D-GaAs with both PBE and HSE functionals. Even though DFT-vdW interactions increase the adsorption energies and reduce the equilibrium distances of H-GaAs systems, it presents, qualitatively, the same physical results on the stability and electronic properties of our studied systems with PBE functional. According to our results, 2D buckled gallium arsenide is a good candidate to be synthesized by hydrogen surface passivation as its group III-V partners 2D buckled gallium nitride and boron nitride. The hydrogenation of 2D-GaAs tunes the bandgap of pristine 2D-GaAs, which makes it a potential candidate for optoelectronic applications in the blue and violet ranges of the visible electromagnetic spectrum.

Hydrogen-Induced High-Temperature Superconductivity in Two-Dimensional Materials: The Example of Hydrogenated Monolayer MgB_{2}.

Hydrogen-based compounds under ultrahigh pressure, such as the polyhydrides H_{3}S and LaH_{10}, superconduct through the conventional electron-phonon coupling mechanism to attain the record critical temperatures known to date. Here we exploit the intrinsic advantages of hydrogen to strongly enhance phonon-mediated superconductivity in a completely different system, namely, a two-dimensional material with hydrogen adatoms. We find that van Hove singularities in the electronic structure, originating from atomiclike hydrogen states, lead to a strong increase of the electronic density of states at the Fermi level, and thus of the electron-phonon coupling. Additionally, the emergence of high-frequency hydrogen-related phonon modes in this system boosts the electron-phonon coupling further. As a concrete example, we demonstrate the effect of hydrogen adatoms on the superconducting properties of monolayer MgB_{2}, by solving the fully anisotropic Eliashberg equations, in conjunction with a first-principles description of the electronic and vibrational states, and their coupling. We show that hydrogenation leads to a high critical temperature of 67 K, which can be boosted to over 100 K by biaxial tensile strain.

Tunable 2D-gallium arsenide and graphene bandgaps in a graphene/GaAs heterostructure: an ab initio study.

The bandgap behavior of 2D-GaAs and graphene have been investigated with van der Waals heterostructured into a yet unexplored graphene/GaAs bilayer, under both uniaxial stress along c axis and different planar strain distributions. The 2D-GaAs bandgap nature changes from [Formula: see text]-K indirect in isolated monolayer to [Formula: see text]-[Formula: see text] direct in graphene/GaAs bilayer. In the latter, graphene exhibits a bandgap of 5 meV. The uniaxial stress strongly affects the graphene electronic bandgap, while symmetric in-plane strain does not open the bandgap in graphene. Nevertheless, it induces remarkable changes on the GaAs bandgap-width around the Fermi level. However, when applying asymmetric in-plane strain to graphene/GaAs, the graphene sublattice symmetry is broken, and the graphene bandgap is open at the Fermi level to a maximum width of 814 meV. This value is much higher than that reported for just graphene under asymmetric strain. The [Formula: see text]-[Formula: see text] direct bandgap of GaAs remains unchanged in graphene/GaAs under different types of applied strain. The analyses of phonon dispersion and the elastic constants yield the dynamical and mechanical stability of the graphene/GaAs system, respectively. The calculated mechanical properties for bilayer heterostructure are better than those of their constituent monolayers. This finding, together with the tunable graphene bandgap not only by the strength but also by the direction of the strain, enhance the potential for strain engineering of ultrathin group-III-V electronic devices hybridized by graphene.

Josephson vortex loops in nanostructured Josephson junctions.

Linked and knotted vortex loops have recently received a revival of interest. Such three-dimensional topological entities have been observed in both classical- and super-fluids, as well as in optical systems. In superconductors, they remained obscure due to their instability against collapse - unless supported by inhomogeneous magnetic field. Here we reveal a new kind of vortex matter in superconductors - the Josephson vortex loops - formed and stabilized in planar junctions or layered superconductors as a result of nontrivial cutting and recombination of Josephson vortices around the barriers for their motion. Engineering latter barriers opens broad perspectives on loop manipulation and control of other possible knotted/linked/entangled vortex topologies in nanostructured superconductors. In the context of Josephson devices proposed to date, the high-frequency excitations of the Josephson loops can be utilized in future design of powerful emitters, tunable filters and waveguides of high-frequency electromagnetic radiation, thereby pushing forward the much needed Terahertz technology.

Free surfaces recast superconductivity in few-monolayer MgB2: Combined first-principles and ARPES demonstration.

Two-dimensional materials are known to harbour properties very different from those of their bulk counterparts. Recent years have seen the rise of atomically thin superconductors, with a caveat that superconductivity is strongly depleted unless enhanced by specific substrates, intercalants or adatoms. Surprisingly, the role in superconductivity of electronic states originating from simple free surfaces of two-dimensional materials has remained elusive to date. Here, based on first-principles calculations, anisotropic Eliashberg theory, and angle-resolved photoemission spectroscopy (ARPES), we show that surface states in few-monolayer MgB2 make a major contribution to the superconducting gap spectrum and density of states, clearly distinct from the widely known, bulk-like σ- and π-gaps. As a proof of principle, we predict and measure the gap opening on the magnesium-based surface band up to a critical temperature as high as ~30 K for merely six monolayers thick MgB2. These findings establish free surfaces as an unavoidable ingredient in understanding and further tailoring of superconductivity in atomically thin materials.

Microfluidic manipulation of magnetic flux domains in type-I superconductors: droplet formation, fusion and fission.

The magnetic flux domains in the intermediate state of type-I superconductors are known to resemble fluid droplets, and their dynamics in applied electric current is often cartooned as a "dripping faucet". Here we show, using the time-depended Ginzburg-Landau simulations, that microfluidic principles hold also for the determination of the size of the magnetic flux-droplet as a function of the applied current, as well as for the merger or splitting of those droplets in the presence of the nanoengineered obstacles for droplet motion. Differently from fluids, the flux-droplets in superconductors are quantized and dissipative objects, and their pinning/depinning, nucleation, and splitting occur in a discretized form, all traceable in the voltage measured across the sample. At larger applied currents, we demonstrate how obstacles can cause branching of laminar flux streams or their transformation into mobile droplets, as readily observed in experiments.

Imaging of super-fast dynamics and flow instabilities of superconducting vortices.

Quantized magnetic vortices driven by electric current determine key electromagnetic properties of superconductors. While the dynamic behavior of slow vortices has been thoroughly investigated, the physics of ultrafast vortices under strong currents remains largely unexplored. Here, we use a nanoscale scanning superconducting quantum interference device to image vortices penetrating into a superconducting Pb film at rates of tens of GHz and moving with velocities of up to tens of km/s, which are not only much larger than the speed of sound but also exceed the pair-breaking speed limit of superconducting condensate. These experiments reveal formation of mesoscopic vortex channels which undergo cascades of bifurcations as the current and magnetic field increase. Our numerical simulations predict metamorphosis of fast Abrikosov vortices into mixed Abrikosov-Josephson vortices at even higher velocities. This work offers an insight into the fundamental physics of dynamic vortex states of superconductors at high current densities, crucial for many applications.Ultrafast vortex dynamics driven by strong currents define eletromagnetic properties of superconductors, but it remains unexplored. Here, Embon et al. use a unique scanning microscopy technique to image steady-state penetration of super-fast vortices into a superconducting Pb film at rates of tens of GHz and velocities up to tens of km/s.

Velocimetry of superconducting vortices based on stroboscopic resonances.

An experimental determination of the mean vortex velocity in superconductors mostly relies on the measurement of flux-flow resistance with magnetic field, temperature, or driving current. In the present work we introduce a method combining conventional transport measurements and a frequency-tuned flashing pinning potential to obtain reliable estimates of the vortex velocity. The proposed device is characterized using the time-dependent Ginzburg-Landau formalism, where the velocimetry method exploits the resonances in mean vortex dissipation when temporal commensuration occurs between the vortex crossings and the flashing potential. We discuss the sensitivity of the proposed technique on applied current, temperature and heat diffusion, as well as the vortex core deformations during fast motion.

Phonon limited superconducting correlations in metallic nanograins.

Conventional superconductivity is inevitably suppressed in ultra-small metallic grains for characteristic sizes smaller than the Anderson limit. Experiments have shown that above the Anderson limit the critical temperature may be either enhanced or reduced when decreasing the particle size, depending on the superconducting material. In addition, there is experimental evidence that whether an enhancement or a reduction is found depends on the strength of the electron-phonon interaction in the bulk. We reveal how the strength of the e-ph interaction interplays with the quantum-size effect and theoretically obtain the critical temperature of the superconducting nanograins in excellent agreement with experimental data. We demonstrate that strong e-ph scattering smears the peak structure in the electronic density-of-states of a metallic grain and enhances the electron mass, and thereby limits the highest Tc achievable by quantum confinement.

Stroboscopic phenomena in superconductors with dynamic pinning landscape.

Introducing artificial pinning centers is a well established strategy to trap quantum vortices and increase the maximal magnetic field and applied electric current that a superconductor can sustain without dissipation. In case of spatially periodic pinning, a clear enhancement of the superconducting critical current arises when commensurability between the vortex configurations and the pinning landscape occurs. With recent achievements in (ultrafast) optics and nanoengineered plasmonics it has become possible to exploit the interaction of light with superconductivity, and create not only spatially periodic imprints on the superconducting condensate, but also temporally periodic ones. Here we show that in the latter case, temporal matching phenomena develop, caused by stroboscopic commensurability between the characteristic frequency of the vortex motion under applied current and the frequency of the dynamic pinning. The matching resonances persist in a broad parameter space, including magnetic field, driving current, or material purity, giving rise to unusual features such as externally variable resistance/impedance and Shapiro steps in current-voltage characteristics. All features are tunable by the frequency of the dynamic pinning landscape. These findings open further exploration avenues for using flashing, spatially engineered, and/or mobile excitations on superconductors, permitting us to achieve advanced functionalities.

Spontaneous symmetry breaking in vortex systems with two repulsive lengthscales.

Scanning Hall probe microscopy (SHPM) has been used to study vortex structures in thin epitaxial films of the superconductor MgB2. Unusual vortex patterns observed in MgB2 single crystals have previously been attributed to a competition between short-range repulsive and long-range attractive vortex-vortex interactions in this two band superconductor; the type 1.5 superconductivity scenario. Our films have much higher levels of disorder than bulk single crystals and therefore both superconducting condensates are expected to be pushed deep into the type 2 regime with purely repulsive vortex interactions. We observe broken symmetry vortex patterns at low fields in all samples after field-cooling from above Tc. These are consistent with those seen in systems with competing repulsions on disparate length scales, and remarkably similar structures are reproduced in dirty two band Ginzburg-Landau calculations, where the simulation parameters have been defined by experimental observations. This suggests that in our dirty MgB2 films, the symmetry of the vortex structures is broken by the presence of vortex repulsions with two different lengthscales, originating from the two distinct superconducting condensates. This represents an entirely new mechanism for spontaneous symmetry breaking in systems of superconducting vortices, with important implications for pinning phenomena and high current density applications.

Giant paramagnetic Meissner effect in multiband superconductors.

Superconductors, ideally diamagnetic when in the Meissner state, can also exhibit paramagnetic behavior due to trapped magnetic flux. In the absence of pinning such paramagnetic response is weak, and ceases with increasing sample thickness. Here we show that in multiband superconductors paramagnetic response can be observed even in slab geometries, and can be far larger than any previous estimate - even multiply larger than the diamagnetic Meissner response for the same applied magnetic field. We link the appearance of this giant paramagnetic response to the broad crossover between conventional Type-I and Type-II superconductors, where Abrikosov vortices interact non-monotonically and multibody effects become important, causing unique flux configurations and their locking in the presence of surfaces.

Interplay of model ingredients affecting aggregate shape plasticity in diffusion-limited aggregation.

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.