Abstract: The envelope function method traditionally employs a single basis set which, in practice, relates to a single material because the k.p matrix elements are generally only known in a particular basis. In this work, we defined a basis function transformation to alleviate this restriction. The transformation is completely described by the known inter-band momentum matrix elements. The resulting envelope function equation can solve the electronic structure in lattice matched heterostructures without resorting to boundary conditions at the interface between materials, while all unit-cell averaged observables can be calculated as with the standard envelope function formalism. In the case of two coupled bands, this heterostructure formalism is equivalent to the standard formalism while taking position dependent matrix elements. (C) 2015 Elsevier B.V. All rights reserved

Abstract: When floating on a two-dimensional surface of superfluid He-4, electrons arrange themselves in two-dimensional crystalline structure known as Wigner crystal. In channels, the boundaries interfere the crystalline order and in case of very narrow channels one observes a quasi-one-dimensional (quasi-1D) Wigner crystal formed by just a few rows of electrons and, ultimately, one row in the “quantum wire” regime. Recently, the “quantum wire” regime was accessed experimentally [D.G. Rees, H. Totsuji, K. Kono, Phys. Rev. Lett. 108, 176801 (2012)] resulting in unusual transport phenomena such as, e.g., oscillations in the electron conductance. Using molecular dynamics simulations, we study the nonlinear transport of electrons in channels with various types of constrictions: single and multiple symmetric and asymmetric geometrical constrictions with varying width and length, and saddle-point-type potentials with varying gate voltage. In particular, we analyze the average particle velocity of the particles and the corresponding electron current versus the driving force or the gate voltage. We have revealed a significant difference in the dynamics for long and short constrictions: The oscillations of the average velocity of the particles for the systems with short constrictions exhibit a clear correlation with the transitions between the states with different numbers of rows of particles; on the other hand, for the systems with longer constrictions these oscillations are suppressed. The obtained results qualitatively agree with the experimental observations. Next, we propose a FET-like structure that consists of a channel with asymmetric constrictions. We show that applying a transverse bias results either in increase of the average particle velocity or in its suppression thus allowing a flexible control tool over the electron transport. The advantage of the asymmetric FET is that it does not have a gate and it allows an easy control of relatively large electron flow. Furthermore, the asymmetric device can be used for rectification of an ac-driven electron flow. Our results bring important insights into the dynamics of electrons floating on the surface of superfluid He-4 in channels with constrictions and allow the effective control over the electron transport.

Abstract: A thorough investigation of three-dimensional kagome lattices of group IV elements is performed with first-principles calculations. The investigated kagome lattices of silicon and germanium are found to be of similar stability as the recently proposed carbon kagome lattice. Carbon and silicon kagome lattices are both direct-gap semiconductors but they have qualitatively different electronic band structures. While direct optical transitions between the valence and conduction bands are allowed in the carbon case, no such transitions can be observed for silicon. The kagome lattice of germanium exhibits semimetallic behavior but can be transformed into a semiconductor after compression.

Abstract: Since the 1960's it has been well known that the basic superconductive quantities can exhibit oscillations as functions of the thickness (diameter) in superconducting nanofilms (nanowires) due to the size quantization of the electronic spectrum. However, very little is known about the effects of quantum confinement on the microscopic properties of vortices. Based on a numerical solution to the Bogoliubov-de Gennes equations, we study the quantum-size oscillations of the vortex core resulting from the sequential interchange of the Kramer-Pesch and anti-Kramer-Pesch regimes with changing nanocylinder radius. The physics behind the anti-Kramer-Pesch anomaly is displayed by utilizing a semi-analytical Anderson approximate solution. We also demonstrate that the anti-Kramer-Pesch vortex core is robust against thermal smearing and results in a distinctive two-maxima structure in the local density of states, which can be used to identify the existence of the anti-Kramer-Pesch vortex.

Abstract: Motivated by a recent experiment that reported the successful synthesis of hexagonal (h) AlN [Tsipas et al., Appl. Phys. Lett. 103, 251605 (2013)], we investigate structural, electronic, and vibrational properties of bulk, bilayer, and monolayer structures of h-AlN by using first-principles calculations. We show that the hexagonal phase of the bulk h-AlN is a stable direct-band-gap semiconductor. The calculated phonon spectrum displays a rigid-layer shear mode at 274 cm(-1) and an E-g mode at 703 cm(-1), which are observable by Raman measurements. In addition, single-layer h-AlN is an indirect-band-gap semiconductor with a nonmagnetic ground state. For the bilayer structure, AA'-type stacking is found to be the most favorable one, and interlayer interaction is strong. While N-layered h-AlN is an indirect-band-gap semiconductor for N = 1 – 9, we predict that thicker structures (N >= 10) have a direct band gap at the Gamma point. The number-of-layer-dependent band-gap transitions in h-AlN is interesting in that it is significantly different from the indirect-to-direct crossover obtained in the transition-metal dichalcogenides.

Abstract: Ripples in pristine freestanding graphene naturally orient themselves in an array that is alternately curved-up and curved-down; maintaining an average height of zero. Using scanning tunneling microscopy (STM) to apply a local force, the graphene sheet will reversibly rise and fall in height until the height reaches 60%-70% of its maximum at which point a sudden, permanent jump occurs. We successfully model the ripples as a spin-half Ising magnetic system, where the height of the graphene plays the role of the spin. The permanent jump in height, controlled by the tunneling current, is found to be equivalent to an antiferromagnetic-to-ferromagnetic phase transition. The thermal load underneath the STM tip alters the local tension and is identified as the responsible mechanism for the phase transition. Four universal critical exponents are measured from our STM data, and the model provides insight into the statistical role of graphene's unusual negative thermal expansion coefficient.

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: 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: The electric-field response of a one-dimensional ring of interacting fermions, where the interactions are described by the extended Hubbard model, is investigated. By using an accurate real-time propagation scheme based on the Chebyshev expansion of the evolution operator, we uncover various nonlinear regimes for a range of interaction parameters that allows modeling of metallic and insulating (either charge density wave or spin density wave insulators) rings. The metallic regime appears at the phase boundary between the two insulating phases and provides the opportunity to describe either weakly or strongly correlated metals. We find that the fidelity susceptibility of the ground state as a function of magnetic flux piercing the ring provides a very good measure of the short-time response. Even completely different interacting regimes behave in a similar manner at short time scales as long as the ground-state fidelity susceptibility is the same. Depending on the strength of the electric field we find various types of responses: persistent currents in the insulating phase, a dissipative regime, or damped Bloch-like oscillations with varying frequencies or even irregular in nature. Furthermore, we also consider the dimerization of the ring and describe the response of a correlated band insulator. In this case the distribution of the energy levels is more clustered and the Bloch-like oscillations become even more irregular.

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: 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.