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: External fields are a powerful tool to probe optical excitations in a material. The linear energy shift of an excitation in a magnetic field is quantified by its effective g factor. Here we show how exciton g factors and their sign can be determined by converged first-principles calculations. We apply the method to monolayer excitons in semiconducting transition metal dichalcogenides and to interlayer excitons in MoSe2/WSe2 heterobilayers and obtain good agreement with recent experimental data. The precision of our method allows us to assign measured g factors of optical peaks to specific transitions in the band structure and also to specific regions of the samples. This revealed the nature of various, previously measured interlayer exciton peaks. We further show that, due to specific optical selection rules, g factors in van der Waals heterostructures are strongly spin and stacking-dependent. The calculation of orbital angular momenta requires the summation over hundreds of bands, indicating that for the considered two-dimensional materials the basis set size is a critical numerical issue. The presented approach can potentially be applied to a wide variety of semiconductors.

Abstract: Two-dimensional crystals with strong interactions between layers has attracted increasing attention in recent years in a variety of fields. In particular, the growth of a single layer of oxide materials (e.g., MgO, CaO, NiO, and MnO) over metallic substrates were found to display different physical properties than their bulk. In this study, we report on the physical properties of a single layer of metallic oxide materials and compare their properties with their bulk and other two-dimensional (2D) crystals. We found that the planar structure of metallic monoxides are unstable whereas the buckled structures are thermodynamically stable. Also, the 2D-MnO and NiO exhibit different magnetic (ferromagnetic) and optical properties than their bulk, whereas band-gap energy and linear stiffness are found to be decreasing from NiO to MgO. Our findings provide insight into oxide thin-film technology applications.

Abstract: Predicting the mechanical and thermal properties of quasi-two-dimensional (2D) transition metal dichalco-genides (TMDs) is an essential task necessary for their implementation in device applications. Although rigorous density-functional-theory-based calculations are able to predict mechanical and electronic properties, mostly they are limited to zero temperature. Classical molecular dynamics facilitates the investigation of temperature-dependent properties, but its performance highly depends on the potential used for defining interactions between the atoms. In this study, we calculated temperature-dependent phonon properties of single-layer TMDs, namely, MoS2, MoSe2, WS2, and WSe2, by utilizing Stillinger-Weber-type potentials with optimized sets of parameters with respect to the first-principles results. The phonon lifetimes and contribution of each phonon mode in thermal conductivities in these monolayer crystals are systematically investigated by means of the spectralenergy-density method based on molecular dynamics simulations. The obtained results from this approach are in good agreement with previously available results from the Green-Kubo method. Moreover, detailed analysis of lattice thermal conductivity, including temperature-dependent mode decomposition through the entire Brillouin zone, shed more light on the thermal properties of these 2D crystals. The LA and TA acoustic branches contribute most to the lattice thermal conductivity, while ZA mode contribution is less because of the quadratic dispersion around the Brillouin zone center, particularly in MoSe2 due to the phonon anharmonicity, evident from the redshift, especially in optical modes, by increasing temperature. For all the considered 2D crystals, the phonon lifetime values are compelled by transition metal atoms, whereas the group velocity spectrum is dictated by chalcogen atoms. Overall, the lattice thermal conductivity is linearly proportional with inverse temperature.

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: 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: The superconducting properties of two recently synthesized ternary silicides with unit formula X₄CoSi (X = Nb, Ta) are investigated through ab initio calculations combined with Eliashberg theory. Interestingly, their crystal structure comprises interlocking honeycomb networks of Nb/Ta atoms. Nb₄CoSi is found to harbor better conditions for phonon-mediated superconductivity, as it possesses a higher density of states at the Fermi level, fostering stronger electron-phonon coupling. The superconducting critical temperatures (T_c) follow the same trend, with Nb₄CoSi having a twice higher value than Ta₄CoSi. Furthermore, the calculated T_c values (5.9 K vs 3.1 K) agree excellently with the experimentally obtained ones, establishing superconductivity in this new materials class as mediated by the electron-phonon coupling. Furthermore, my calculations show that the superconducting properties of these compounds do not simply correlate with the parameters of their honeycomb networks, contrary to proposals raised in the literature. Rather, their complete fermiology and phonon spectrum should be taken into account in order to explain their respective superconducting properties.

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: We investigate the intralayer correlations acting within the layers in a superfluid system of electron -hole spatially separated layers. In this system, superfluidity is predicted to be almost exclusively confined to the Bose-Einstein condensate (BEC) and crossover regimes where the electron -hole pairs are well localized. In this case, Hartree-Fock is an excellent approximation for the intralayer correlations. We find in the BEC regime that the effect of the intralayer correlations on superfluid properties is negligible but in the BCS-BEC crossover regime the superfluid gap is significantly weakened by the intralayer correlations. This is caused by the intralayer correlations boosting the number of low -energy particle -hole excitations that drive the screening. We further find that the intralayer correlations suppress the predicted phenomenon in which the average pair size passes through a minimum as the crossover regime is traversed.

Abstract: Transport measurements are readily used to probe different phases in disordered topological insulators (TIs), where determining topological invariants explicitly is challenging. On that note, universal conductance fluctuations (UCF) theory asserts the conductance G for an ensemble has a Gaussian distribution, and that standard deviation 8G depends solely on the symmetries and dimensions of the system. Using a real-space tight -binding Hamiltonian on a system with Anderson disorder, we explore conductance fluctuations in a thin Bi2Se3 film and demonstrate the agreement of their behavior with UCF hypotheses. We further show that magnetic field applied out-of-plane breaks the time -reversal symmetry and transforms the system's Wigner-Dyson class from root symplectic to unitary, increasing 8G by 2. Finally, we reveal that while Bi2Se3 is a strong TI, weak TI and metallic phases can be stabilized in presence of strain and disorder, and detected by monitoring the conductance fluctuations.

Abstract: In a recent paper, Arita et al. [Phys. Rev. E 90, 052108 (2014)] consider the transport properties of a class of generalized exclusion processes. Analytical expressions for the transport-diffusion coefficient are derived by ignoring correlations. It is claimed that these expressions become exact in the hydrodynamic limit. In this Comment,we point out that (i) the influence of correlations upon the diffusion does not vanish in the hydrodynamic limit, and (ii) the expressions for the self- and transport diffusion derived by Arita et al. are special cases of results derived in Becker et al. [Phys. Rev. Lett. 111, 110601 (2013)].

Abstract: Molecular dynamics simulations are used to investigate the structural properties of a two-dimensional ensemble of magnetic rods, which are modeled as aligned single dipolar beads. The obtained self-assembled configurations can be characterized as (1) clusters, (2) percolated, and (3) ordered structures, and their structural properties are investigated in detail. By increasing the aspect ratio of the magnetic rods, we show that the percolation transition is suppressed due to the reduced mobility of the rods in two dimensions. Such a behavior is opposite to the one observed in three dimensions. A magnetic bulk phase is found with local ferromagnetic order and an unusual nonmonotonic behavior of the nematic order is observed.

Abstract: The dielectric properties of confined water is of fundamental interest and is still controversial. For water confined in channels with height smaller than h = 8 angstrom, we found a commensurability effect and an extraordinary decrease in the out-of-plane dielectric constant down to the limit of the dielectric constant of optical water. Spatial resolved polarization density data obtained from molecular dynamics simulations are found to be antisymmetric across the channel and are used as input in a mean-field model for the dielectric constant as a function of the height of the channel for h > 15 angstrom. Our results are in excellent agreement with a recent experiment [L. Fumagalli et al., Science 360, 1339 (2018)].

Abstract: Molecular dynamics simulations are used to study the effects of different cations on the permeation of charged polymers through flat capillaries with heights below 2 nm. Interestingly, we found that, despite being monovalent, Li+ , Na+ , and K+ cations have different effects on polymer permeation, which consequently affects their transmission speed throughout those capillaries. We attribute this phenomenon to the interplay of the cations' hydration free energies and the hydrodynamic drag in front of the polymer when it enters the capillary. Different alkali cations exhibit different surface versus bulk preferences in small clusters of water under the influence of an external electric field. This paper presents a tool to control the speed of charged polymers in confined spaces using cations.

Abstract: The confinement of a C-60 molecule encapsulated in a cylindrical nanotube depends on the tube radius. In small tubes with radius R-T less than or similar to 7 A, a fivefold axis of the molecule coincides with the tube axis. The interaction between C-60 molecules in the nanotube is then described by a O-2-rotor model on a 1D liquid chain with coupling between orientational and displacive correlations. This coupling leads to chain contraction. The structure factor of the 1D liquid is derived. In tubes with a larger radius the molecular centers of mass are displaced off the tube axis. The distinction of two groups of peapods with on- and off-axis molecules suggests an explanation of the apparent splitting of A(g) modes of C-60 in nanotubes measured by resonant Raman scattering.