Abstract: Molecular-dynamic simulations were performed in order to investigate the melting processes of isotropically confined binary systems. We considered two species of particles, which differ by their amount of electric charge. A Lindemann type of criterion was used to determine the angular melting temperature. We demonstrate that the magic-to-normal cluster transition can evolve in two distinct ways, that is, through a structural phase transition of the first order or via a smooth transition where an increase of the shells' width leads to a continuous decreasing mechanical stability of the system. Moreover, for large systems, we demonstrate that the internal cluster exerts a minor effect on the mechanical stability of the external shell. Furthermore, we show that highly symmetric configurations, such as those found for multiple ring structures, have large mechanical stability, i.e., high angular melting temperature.

Abstract: We expand on a recent study of a lattice model of interacting particles [Phys. Rev. Lett. 111, 110601 (2013)]. The adsorption isotherm and equilibrium fluctuations in particle number are discussed as a function of the interaction. Their behavior is similar to that of interacting particles in porous materials. Different expressions for the particle jump rates are derived from transition-state theory. Which expression should be used depends on the strength of the interparticle interactions. Analytical expressions for the self-and transport diffusion are derived when correlations, caused by memory effects in the environment, are neglected. The diffusive behavior is studied numerically with kinetic Monte Carlo (kMC) simulations, which reproduces the diffusion including correlations. The effect of correlations is studied by comparing the analytical expressions with the kMC simulations. It is found that the Maxwell-Stefan diffusion can exceed the self-diffusion. To our knowledge, this is the first time this is observed. The diffusive behavior in one-dimensional and higher-dimensional systems is qualitatively the same, with the effect of correlations decreasing for increasing dimension. The length dependence of both the self-and transport diffusion is studied for one-dimensional systems. For long lengths the self-diffusion shows a 1/L dependence. Finally, we discuss when agreement with experiments and simulations can be expected. The assumption that particles in different cavities do not interact is expected to hold quantitatively at low and medium particle concentrations if the particles are not strongly interacting.

Abstract: Phase-slip centers/lines and hot spots are the main mechanisms for dissipation in current-carrying superconducting thin films. The pulsed-current method has recently been shown to be an effective tool in studying the dynamics of phase-slip centers and their evolution to hot spots. We use the time-dependent Ginzburg-Landau theory in the study of the dynamics of the superconducting condensate in superconducting strips under external current and zero external magnetic field. We show that both the flux-flow state (i.e., slow-moving vortices) and the phase-slip line state (i.e., fast-moving vortices) are dynamically stable dissipative units with temperature smaller than the critical one, whereas hot spots, which are localized normal regions where the local temperature exceeds the critical value, expand in time, resulting ultimately in a complete destruction of the condensate. The response time of the system to abrupt switching on of the overcritical current decreases with increasing both the value of the current (at all temperatures) and temperature (for a given value of the applied current). Our results are in good qualitative agreement with experiments we have conducted on Nb thin strips.

Abstract: We investigate the dependence of the structural phase transitions in an infinite quasi-one-dimensional system of repulsively interacting particles on the profile of the confining channel. Three different functional expressions for the confinement potential related to real experimental systems are used that can be tuned continuously from a parabolic to a hard-wall potential in order to find a thorough understanding of the ordering of the chainlike structure transitions. We resolve the long-standing issue why the most theories predicted a 1-2-4-3-4 sequence of chain configurations with increasing density, while some experiments found the 1-2-3-4 sequence.

Abstract: The influence of silicene's strong spin-orbit interaction and of an external electric field E-z on the transport coefficients are investigated in the presence of a perpendicular magnetic field B. For finite E-z the spin and valley degeneracy of the Landau levels is lifted and leads to additional plateaus in the Hall conductivity, at half-integer values of 4e(2)/h, due to spin intra-Landau-level transitions that are absent in graphene. These plateaus are more sensitive to disorder and thermal broadening than the main plateaus, occurring at integral values of 4e(2)/h, when the Fermi level passes through the Landau levels. We also evaluate the Hall and longitudinal resistivities and critically contrast the results with those for graphene on a substrate.

Abstract: Carbon-based molecular semiconductors are explored for application in spintronics because their small spinorbit coupling promises long spin lifetimes. We calculate the electronic transport from first principles through spin valves comprising bi-and tri-layers of the fullerene molecules C-60 and C-70, sandwiched between two Fe electrodes. The spin polarization of the current, and the magnetoresistance depend sensitively on the interactions at the interfaces between the molecules and the metal surfaces. They are much less affected by the thickness of the molecular layers. A high current polarization (CP > 90%) and magnetoresistance (MR > 100%) at small bias can be attained using C-70 layers. In contrast, the current polarization and the magnetoresistance at small bias are vanishingly small for C-60 layers. Exploiting a generalized Julliere model we can trace the differences in spin-dependent transport between C-60 and C-70 layers to differences between the molecule-metal interface states. These states also allow one to interpret the current polarization and the magnetoresistance as a function of the applied bias voltage.

Abstract: Twisted graphene layers produce a moire pattern (MP) structure with a predetermined wavelength for a given twist angle. However, predicting the membrane corrugation amplitude for any angle other than pure AB-stacked or AA-stacked graphene is impossible using first-principles density functional theory (DFT) due to the large supercell. Here, within elasticity theory, we define the MP structure as the minimum-energy configuration, thereby leaving the height amplitude as the only unknown parameter. The latter is determined from DFT calculations for AB-and AA-stacked bilayer graphene in order to eliminate all fitting parameters. Excellent agreement with scanning tunneling microscopy results across multiple substrates is reported as a function of twist angle.

Abstract: In honeycomb Dirac systems with broken inversion symmetry, orbital magnetic moments coupled to the valley degree of freedom arise due to the topology of the band structure, leading to valley-selective optical dichroism. On the other hand, in Dirac systems with prominent spin-orbit coupling, similar orbital magnetic moments emerge as well. These moments are coupled to spin, but otherwise have the same functional form as the moments stemming from spatial inversion breaking. After reviewing the basic properties of these moments, which are relevant for a whole set of newly discovered materials, such as silicene and germanene, we study the particular impact that these moments have on graphene nanoengineered barriers with artificially enhanced spin-orbit coupling. We examine transmission properties of such barriers in the presence of a magnetic field. The orbital moments are found to manifest in transport characteristics through spin-dependent transmission and conductance, making them directly accessible in experiments. Moreover, the Zeeman-type effects appear without explicitly incorporating the Zeeman term in the models, i.e., by using minimal coupling and Peierls substitution in continuum and the tight-binding methods, respectively. We find that a quasiclassical view is able to explain all the observed phenomena.

Abstract: Quantum thermal transport in armchair and zigzag graphene nanoribbons is investigated in the presence of single atomic vacancies and subject to different boundary conditions. We start with a full comparison of the phonon polarizations and energy dispersions as given by a fifth-nearest-neighbor force-constant model (5NNFCM) and by elasticity theory of continuum membranes (ETCM). For free-edge ribbons, we discuss the behavior of an additional acoustic edge-localized flexural mode, known as fourth acoustic branch (4ZA), which has a small gap when it is obtained by the 5NNFCM. Then, we show that ribbons with supported edges have a sample-size dependent energy gap in the phonon spectrum which is particularly large for in-plane modes. Irrespective to the calculation method and the boundary condition, the dependence of the energy gap for the low-energy optical phonon modes against the ribbon width W is found to be proportional to 1/W for in-plane, and 1/W-2 for out-of-plane phonon modes. Using the 5NNFCM, the ballistic thermal conductance and its contributions from every single phonon mode are then obtained by the nonequilibrium Green's function technique. We found that, while edge and central localized single atomic vacancies do not affect the low-energy transmission function of in-plane phonon modes, they reduce considerably the contributions of the flexural modes. On the other hand, in-plane modes contributions are strongly dependent on the boundary conditions and at low temperatures can be highly reduced in supported-edge samples. These findings could open a route to engineer graphene based devices where it is possible to discriminate the relative contribution of polarized phonons and to tune the thermal transport on the nanoscale.

Abstract: The electronic properties of the two-dimensional material silicene are strongly influenced by the application of a perpendicular electric field E-z and of an exchange field M due to adatoms positioned on the surface or a ferromagnetic substrate. Within the random phase approximation, we investigate how electron-electron interactions are affected by these fields and present analytical and numerical results for the dispersion of plasmons, their lifetime, and their oscillator strength. We find that the combination of the fields E-z and M brings a spin and valley texture to the particle-hole excitation spectrum and allows the formation of spin-and valley-polarized plasmons. When the Fermi level lies in the gap of one spin in one valley, the intraband region of the corresponding spectrum disappears. For zero E-z and finite M the spin symmetry is broken and spin polarization is possible. The lifetime and oscillator strength of the plasmons are shown to depend strongly on the number of spin and valley type electrons that form the electron-hole pairs.

Abstract: We revisit the problem of the dynamic response of a superconducting bridge after abruptly switching on the supercritical current. In contrast to previous theoretical works we take into account spatial gradients and use both the local temperature approach and the kinetic equation for the distribution function of quasiparticles. We find that the temperature dependence of the finite delay time t(d) in the voltage response is model dependent and relatively large t(d) is connected with temporary cooling of quasiparticles during decay of superconducting order parameter vertical bar Delta vertical bar in time. It turns out that the presence of even small inhomogeneities in the bridge or finite length of the homogenous bridge favors a local suppression of vertical bar Delta vertical bar during the dynamic response. It results in a decrease of the delay time, in comparison with the spatially uniform model, due to the diffusion of nonequilibrium quasiparticles from the region with locally suppressed vertical bar Delta vertical bar. In the case when the current density is maximal near the edge of a not very wide bridge the delay time is mainly connected with the time needed for the nucleation (entrance) of the first vortex and t(d) could be tuned by a weak external magnetic field. We also find that a short alternating current pulse (sinusoidlike) with zero time average may result in a nonzero time- averaged voltage response where its sign depends on the phase of the ac current.

Abstract: Transition metal dichalcogenides (TMDs) are considered as promising candidates for next generation of electronic and optoelectronic devices. To make use of these materials, for instance in field effect transistor applications, it is mandatory to know the detailed properties of contacts of such TMDs with metal electrodes. Here, we investigate the role of the contact structure on the electronic and transport properties of metal-MoSe2 interfaces. Two different contact types, namely face and edge contacts, are studied. We consider both low (Sc) and high (Au) work function metals in order to thoroughly elucidate the role of the metal work function and the type of metal. First principles plane wave calculations and transport calculations based on nonequilibrium Green's function formalism reveal that the contact type has a large impact on the electronic and transport properties of metal-MoSe2 interfaces. For the Sc electrode, the Schottky barrier heights are around 0.25 eV for face contact and bigger than 0.6 eV for edge contact. For the Au case, we calculate very similar barrier heights for both contact types with an average value of 0.5 eV. Furthermore, while the face contact is found to be highly advantageous as compared to the edge contact for the Sc electrode, the latter contact becomes a better choice for the Au electrode. Our findings provide guidelines for the fabrication of TMD-based devices.

Abstract: Recently, rhenium disulfide (ReS2) monolayers were experimentally extracted by conventional mechanical exfoliation technique from as-grown ReS2 crystals. Unlike the well-known members of transition metal dichalcogenides (TMDs), ReS2 crystallizes in a stable distorted-1T structure and lacks an indirect to direct gap crossover. Here we present an experimental and theoretical study of the formation, energetics, and stability of the most prominent lattice defects in monolayer ReS2. Experimentally, irradiation with 3-MeV He+2 ions was used to break the strong covalent bonds in ReS2 flakes. Photoluminescence measurements showed that the luminescence from monolayers is mostly unchanged after highly energetic a particle irradiation. In order to understand the energetics of possible vacancies in ReS2 we performed systematic first-principles calculations. Our calculations revealed that the formation of a single sulfur vacancy has the lowest formation energy in both Re and S rich conditions and a random distribution of such defects are energetically more preferable. Sulfur point defects do not result in any spin polarization whereas the creation of Re-containing point defects induce magnetization with a net magnetic moment of 1-3 mu B. Experimentally observed easy formation of sulfur vacancies is in good agreement with first-principles calculations.

Abstract: In order to account explicitly for the existence of long-periodic layered structures and the strong structural relaxations in the most common binary and ternary alloys of the Bi-Sb-Te-Se system, we have developed a one-dimensional cluster expansion (CE) based on first-principles electronic structure calculations, which accounts for the Bi and Sb bilayer formation. Excellent interlayer distances are obtained with a van der Waals density functional. It is shown that a CE solely based on pair interactions is sufficient to provide an accurate description of the ground-state energies of Bi-Sb-Te-Se binary and ternary systems without making the data set of ab initio calculated structures unreasonably large. For the binary alloys A1−xQx (A=Sb, Bi; Q=Te, Se), a ternary CE yields an almost continuous series of (meta)stable structures consisting of consecutive A bilayers next to consecutive A2Q3 for 0<x<0.6. For x>0.6, the binary alloy segregates into pure Q and A2Q3. The Bi-Sb system is described by a quaternary CE and is found to be an ideal solid solution stabilized by entropic effects at T≠0 K but with an ordered structure of alternating Bi and Sb layers for x=0.5 at T=0 K. A quintuple CE is used for the ternary Bi-Sb-Te system, where stable ternary layered compounds with an arbitrary stacking of Sb2Te3, Bi2Te3, and Te-Bi-Te-Sb-Te quintuple units are found, optionally separated by mixed Bi/Sb bilayers. Electronic properties of the stable compounds were studied taking spin-orbit coupling into account.