Abstract: Bilayer (BL) molybdenum disulfide (MoS2) is one of the most important electronic structures not only in valleytronics but also in realizing twistronic systems on the basis of the topological mosaics in moire superlattices. In this work, BL MoS2 on sapphire substrate with 2H-stacking structure is fabricated. We apply the terahertz (THz) time-domain spectroscopy (TDS) for examining the basic optoelectronic properties of this kind of BL MoS2. The optical conductivity of BL MoS2 is obtained in temperature regime from 80 K to 280 K. Through fitting the experimental data with the theoretical formula, the key sample parameters of BL MoS2 can be determined, such as the electron density, the electronic relaxation time and the electronic localization factor. The temperature dependence of these parameters is examined and analyzed. We find that, similar to monolayer (ML) MoS2, BL MoS2 with 2H-stacking can respond strongly to THz radiation field and show semiconductor-like optoelectronic features. The theoretical calculations using density functional theory (DFT) can help us to further understand why the THz optoelectronic properties of BL MoS2 differ from those observed for ML MoS2. The results obtained from this study indicate that the THz TDS can be applied suitably to study the optoelectronic properties of BL MoS2 based twistronic systems for novel applications as optical and optoelectronic materials and devices.

Abstract: We use atomistic-based simulations to investigate the behavior of ferromagnetic skyrmions being continuously compressed against a rigid wall under dc and ac drives. The compressed skyrmions can be annihilated close to the wall and form a conformal crystal with both a size and a density gradient, making it distinct from conformal crystals observed previously for superconducting vortices and colloidal particles. For both dc and ac driving, the skyrmions can move transverse to the compression direction due to a combination of density and size gradients. Forces in the compression direction are converted by the Magnus force into transverse motion. Under ac driving, the amount of skyrmion annihilation is reduced and we find a skyrmion Magnus ratchet pump. We also observe shear banding in which skyrmions near the wall move up to twice as fast as skyrmions further from the wall. When we vary the magnitude of the applied drive, we find a critical current above which the skyrmions are completely annihilated during a time scale that depends on the magnitude of the drive. By varying the magnetic parameters, we find that the transverse motion is strongly dependent on the skyrmion size. Smaller skyrmions are more rigid, which interferes with the size gradient and destroys the transverse motion. We also confirm the role of the size gradient by comparing our atomistic simulations with a particle-based model, where we find that the transverse motion is only transient. Our results are relevant for applications where skyrmions encounter repulsive magnetic walls, domain walls, or interfaces.

Abstract: The peculiar potential transition metal dichalcogenides in regard to sensor and device applications have been exhibited by both experimental and theoretical studies. The use of these materials, thermodynamically stable even at elevated temperatures, particularly in nano- and optoelectronic technology, is about to come true. On the other hand, the distinct electronic and thermal transport properties possessing unique coherency, which may result in higher thermoelectric efficiency, have also been reported. However, exploiting this potential in terms of power generation and cooling applications requires a deeper understanding of these materials in this regard. This perspective study, concentrated with this intention, summarizes thermoelectric research based on transition metal dichalcogenides from a broad perspective and also provides a general evaluation of future theoretical investigations inevitable to shed more light on the physics of electronic and thermal transport in these materials and to lead future experimental research.

Abstract: Two-dimensional materials (2DMs) continue to attract a lot of attention, particularly for their extreme flexibility and superior thermal properties. Molecular dynamics simulations are among the most powerful methods for computing these properties, but their reliability depends on the accuracy of interatomic interactions. While first principles approaches provide the most accurate description of interatomic forces, they are computationally expensive. In contrast, classical force fields are computationally efficient, but have limited accuracy in interatomic force description. Machine learning interatomic potentials, such as Gaussian Approximation Potentials, trained on density functional theory (DFT) calculations offer a compromise by providing both accurate estimation and computational efficiency. In this work, we present a systematic procedure to develop Gaussian approximation potentials for selected 2DMs, graphene, buckled silicene, and h-XN (X = B, Al, and Ga, as binary compounds) structures. We validate our approach through calculations that require various levels of accuracy in interatomic interactions. The calculated phonon dispersion curves and lattice thermal conductivity, obtained through harmonic and anharmonic force constants (including fourth order) are in excellent agreement with DFT results. HIPHIVE calculations, in which the generated GAP potentials were used to compute higher-order force constants instead of DFT, demonstrated the first-principles level accuracy of the potentials for interatomic force description. Molecular dynamics simulations based on phonon density of states calculations, which agree closely with DFT-based calculations, also show the success of the generated potentials in high-temperature 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 detail the effects of Cl and S functionalization on the superconducting properties of layered (bulk) and monolayer niobium carbide (Nb2C) MXene crystals, based on first-principles calculations combined with Eliashberg theory. For bulk layered Nb2CCl2, the calculated superconducting transition temperature (T-c) is in very good agreement with the recently measured value of 6 K. We show that T-c is enhanced to 10 K for monolayer Nb2CCl2, due to an increase in the density of states at the Fermi level, and the corresponding electron-phonon coupling. We further demonstrate feasible gate- and strain-induced enhancements of T-c for both bulk-layered and monolayer Nb2CCl2 crystals, resulting in T-c values of around 38 K. In the S-functionalized Nb2CCl2 crystals, our calculations reveal the importance of phonon softening in understanding their superconducting properties. Finally, we predict that Nb3C2S2 in bulk-layered and monolayer forms is also superconducting, with a T-c of around 28 K. Considering that Nb2C is not superconducting in pristine form, our findings promote functionalization as a pathway towards robust superconductivity in MXenes.

Abstract: By using the full 3D generalized time-dependent Ginzbug-Landau equation, we study a long superconducting film of finite width and thickness under an applied transport current. We show that, for sufficiently large thickness, the vortices and the antivortices become curved before they annihilate each other. As they approach the center of the sample, their ends combine, producing a single closed vortex. We also determine the critical values of the thickness for which the closed vortex sets in for different values of the Ginzburg-Ladau parameter. Finally, we propose a model of how to detect a closed vortex experimentally.

Abstract: Van der Waals quantum wells, naturally formed in two-dimensional layered materials with nanoscale thickness, possess many inherent advantages over conventional molecular beam epitaxy grown counterparts, and could bring up intriguing physics and applications. However, optical transitions originated from the series of quantized states in these emerging quantum wells are still elusive. Here, we show that multilayer black phosphorus appears to be an excellent candidate for van der Waals quantum wells with well-defined subbands and high optical quality. Using infrared absorption spectroscopy, we probe subband structures of multilayer black phosphorus with tens of atomic layers, revealing clear signatures for optical transitions with subband index as high as 10, far from what was attainable previously. Surprisingly, in addition to allowed transitions, an unexpected series of “forbidden” transitions is also evidently observed, which enables us to determine energy spacings separately for conduction and valence subbands. Furthermore, the linear tunability of subband spacings by temperature and strain is demonstrated. Our results are expected to facilitate potential applications for infrared optoelectronics based on tunable van der Waals quantum wells.

Abstract: investigate non-equilibrium states created by a laser beam incident on a superconducting NbTiN filament subject to an electrical pulse at 4 K. In absence of the laser excitation, when the amplitude of the current pulse applied to the filament exceeds the critical current value, we monitored the delay time td that marks the collapse of the superconducting phase which is then followed by a voltage rise. We linked the delay time to the applied current using the time-dependent Ginzburg-Landau (TDGL) theory, which enabled us to deduce the cooling (or heat-removal) time from the fit to the experimental data. Subsequently, we exposed the filament biased with a current pulse close to its critical value to a focused laser beam, inducing a normal state in the impact region of the laser beam. We showed that the energy of the incident beam and the incurred delay time are related to each other by a simple expression, that enables direct measurement of incident beam energy by temporal monitoring of the transport response. This method can be extended for usage in single-photon detection regime, and be used for accurate calibration of an arbitrary light source.

Abstract: When using hexagonal boron-nitride (hBN) as a substrate for graphene, the resulting moire pattern creates secondary Dirac points. By encapsulating a multilayer graphene within aligned hBN sheets the controlled moire stacking may offer even richer benefits. Using advanced tight-binding simulations on atomistically-relaxed heterostructures, here we show that the gap at the secondary Dirac point can be opened in selected moire-stacking configurations, and is independent of any additional vertical gating of the heterostructure. On the other hand, gating can broadly tune the gap at the principal Dirac point, and may thereby strongly compress the first moire mini-band in width against the moire-induced gap at the secondary Dirac point. We reveal that in hBN-encapsulated bilayer graphene this novel mechanism can lead to isolated bands flatter than 10 meV under moderate gating, hence presenting a convenient pathway towards electronically-controlled strongly-correlated states on demand.

Abstract: Flat bands in the energy spectrum have attracted a lot of attention in recent years because of their unique properties and promising applications. Special arrangement of impurities on monolayer graphene are proposed to generate multiflat bands in the electronic band structure. In addition to the single midgap states in the spectrum of graphene with low hydrogen density, we found closely spaced bands around the Fermi level with increasing impurity density, which are similar to discrete lines in the spectrum of quantum dots, as well as the unusual Landau-level energy spectrum of graphene in the presence of a strong magnetic field. The presence of flat bands crucially depends on whether or not there are odd or even electrons of H(F) atoms bound to graphene. Interestingly, we found that a fully hydrogenated (fluoridated) of a hexagon of graphene sheet with six hydrogen (fluorine) atoms sitting on top and bottom in consecutive order exhibits Dirac cones in the electronic band structure with a 20% smaller Fermi velocity as compared to the pristine graphene. Functionalizing graphene introduces various C-C bond lengths resulting in nonuniform strains. Such a nonuniform strain may induce a giant pseudomagnetic field in the system, resulting in quantum Hall effect.

Abstract: To design fast memory devices, we need material combinations that can facilitate fast read and write operations. We present a heterostructure comprising a two-dimensional (2D) magnet and a 2D topological insulator (TI) as a viable option for designing fast memory devices. We theoretically model the spin-charge dynamics between 2D magnets and 2D TIs. Using the adiabatic approximation, we combine the nonequi-librium Green's function method for spin-dependent electron transport and a time-quantified Monte Carlo method for simulating magnetization dynamics. We show that it is possible to switch a magnetic domain of a ferromagnet using the spin torque from spin-polarized edge states of a 2D TI. We show further that the switching of 2D magnets by TIs is strongly dependent on the interface exchange (Jint), and an opti-mal interface exchange, is required for efficient switching. Finally, we compare experimentally grown Cr compounds and show that Cr compounds with higher anisotropy (such as CrI3) result in a lower switching speed but a more stable magnetic order.

Abstract: High-performance rechargeable batteries are becoming very important for high-end technologies with their ever increasing application areas. Hence, improving the performance of such batteries has become the main bottleneck to transferring high-end technologies to end users. In this study, we propose an argon intercalation strategy to enhance battery performance via engineering the interlayer spacing of honeycomb structures such as graphite, a common electrode material in lithium-ion batteries (LIBs). Herein, we systematically investigated the LIB performance of graphite and hexagonal boron nitride (h-BN) when argon atoms were sent into between their layers by using first-principles density-functional-theory calculations. Our results showed enhanced lithium binding for graphite and h-BN structures when argon atoms were intercalated. The increased interlayer space doubles the gravimetric lithium capacity for graphite, while the volumetric capacity also increased by around 20% even though the volume was also increased. The ab initio molecular dynamics simulations indicate the thermal stability of such graphite structures against any structural transformation and Li release. The nudged-elastic-band calculations showed that the migration energy barriers were drastically lowered, which promises fast charging capability for batteries containing graphite electrodes. Although a similar level of battery promise was not achieved for h-BN material, its enhanced battery capabilities by argon intercalation also support that the argon intercalation strategy can be a viable route to enhance such honeycomb battery electrodes.

Abstract: Using a particle based model, we investigate the skyrmion dynamical behavior in a channel where the upper wall contains divots of one depth and the lower wall contains divots of a different depth. Under an applied driving force, skyrmions in the channels move with a finite skyrmion Hall angle that deflects them toward the upper wall for -x direction driving and the lower wall for +x direction driving. When the upper divots have zero height, the skyrmions are deflected against the flat upper wall for -x direction driving and the skyrmion velocity depends linearly on the drive. For +x direction driving, the skyrmions are pushed against the lower divots and become trapped, giving reduced velocities and a nonlinear velocity-force response. When there are shallow divots on the upper wall and deep divots on the lower wall, skyrmions get trapped for both driving directions; however, due to the divot depth difference, skyrmions move more easily under -x direction driving, and become strongly trapped for +x direction driving. The preferred -x direction motion produces what we call a Magnus diode effect since it vanishes in the limit of zero Magnus force, unlike the diode effects observed for asymmetric sawtooth potentials. We show that the transport curves can exhibit a series of jumps or dips, negative differential conductivity, and reentrant pinning due to collective trapping events. We also discuss how our results relate to recent continuum modeling on a similar skyrmion diode system.

Abstract: The miniaturization of transistors for application in new processors and logic devices poses a significant challenge in the field of materials. Spintronics, which relies on controlled movement of magnetic nanostructures, offers a promising solution. Among the candidates, magnetic skyrmions are considered one of the most promising. These chiral spin structures, characterized by topological protection and enhanced stability compared to vortices or magnetic bubbles, have been extensively studied. To advance in the control of skyrmion motion, essential for practical applications, we investigated their dynamic behavior in a two-dimensional chiral magnet at zero temperature. Our study focused on the influence of periodic arrays of pinning centers. The simulations considered skyrmions as point-like particles considering the following interactions: skyrmion-skyrmion interactions, interactions with pinning center arrays, a current of polarized spins, and the Magnus force. We conducted calculations for scenarios involving a single skyrmion as well as different skyrmion density values in the material. The aim was to explore possibilities for controlled skyrmion motion, investigate different dynamic regimes, and examine collective effects. The results demonstrate that by adjusting the size, strength, and density of the pinning centers, we can effectively control the motion of individual skyrmions and manage the flow of multiple skyrmions. Furthermore, we discovered that periodic arrays of pinning centers can facilitate topological selection when different species of skyrmions with distinct Magnus components are present. Employing alternating currents, we observed the significant role of the ratchet effect in the skyrmion dynamics. By fine-tuning the amplitudes of the alternating currents, we achieved direct and controlled motion of skyrmions in specific directions. These findings hold potential for advancing our understanding of skyrmion dynamics and can inspire future technological applications involving these quasi-particles. Overall, we anticipate that our results will be valuable to the scientific community, contributing to a deeper comprehension of skyrmion dynamics and paving the way for future technological applications.