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: We numerically study the vortex-vortex interaction in multicomponent homogeneous Bose-Einstein condensates within the realm of the Gross-Pitaevskii theory. We provide strong evidence that pairwise vortex interaction captures the underlying mechanisms which determine the geometric configuration of the vortices, such as different lattices in many-vortex states, as well as the bound vortex states with two (dimer) or three (trimer) vortices. Specifically, we discuss and apply our theoretical approach to investigate intra- and intercomponent vortex-vortex interactions in two- and three-component Bose-Einstein condensates, thereby shedding light on the formation of the exotic vortex configurations. These results correlate with current experimental efforts in multicomponent Bose-Einstein condensates and the understanding of the role of vortex interactions in multiband superconductors.

Abstract: The control of side gates on the quantum electron transport in narrow graphene ribbons of different widths and edge types (armchair and zigzag) is investigated. The conductance exhibits Fano resonances with varying side gate potential. Resonant and antiresonant peaks in the conductance can be associated with the eigenstates of a closed system, and these peaks can be accurately fitted with a Fano line shape. The local density of states (LDOS) and the electron current show a specific behavior at these resonances, which depends on the ribbon edge type. In zigzag ribbons, transport is dominated by intervalley scattering, which is reflected in the transmission functions of individual modes. The side gates induce p-n interfaces near the edges at which the LDOS exhibits peaks. Near the resonance points, the electron current flows uniformly through the constriction, while near the antiresonances it creates vortices. In the armchair ribbons the LDOS spreads in areas of high potential, with current flowing near the edges.

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: 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: We numerically investigate the mixed pairing symmetries inmesoscopic superconducting loops in the presence of spin correlations by solving the Bogoliubov-de Gennes equations self-consistently. The spatial variations of the superconducting order parameters and the spontaneous magnetization are determined by the band structure. When the threaded magnetic flux turns on, the charge and spin currents both emerge and depict periodic evolution. In the case of a mesoscopic loop with dominant triplet p(x) +/- ip(y)-wave symmetry, a slight change of the chemical potential may lead to novel flux-dependent evolution patterns of the ground-state energy and the magnetization. The spin-polarized currents show pronounced quantum oscillations with fractional periods due to the appearance of energy jumps in flux, accompanied with a steplike feature of the enhanced spin current. Particularly, at some appropriate flux, the peaks of the zero-energy local density of states clearly indicate the occurrence of the odd-frequency pairing. In the case of a superconducting loop with dominant singlet d(x2-y2)-wave symmetry, the spatial profiles of the zero-energy local density of states and the magnetization show spin-dependent features on different sample diagonals. Moreover, the evolution of the flux-induced spin current always exhibits an hc/e periodicity.

Abstract: Ca(OH)(2) crystals, well known as portlandite, are grown in layered form, and we found that they can be exfoliated on different substrates. We performed first principles calculations to investigate the structural, electronic, vibrational, and mechanical properties of bulk, bilayer, and monolayer structures of this material. Different from other lamellar structures such as graphite and transition-metal dichalcogenides, intralayer bonding in Ca(OH)(2) is mainly ionic, while the interlayer interaction remains a weak dispersion-type force. Unlike well-known transition-metal dichalcogenides that exhibit an indirect-to-direct band gap crossover when going from bulk to a single layer, Ca(OH)(2) is a direct band gap semiconductor independent of the number layers. The in-plane Young's modulus and the in-plane shear modulus of monolayer Ca(OH)(2) are predicted to be quite low while the in-plane Poisson ratio is larger in comparison to those in the monolayer of ionic crystal BN. We measured the Raman spectrum of bulk Ca(OH)(2) and identified the high-frequency OH stretching mode A(1g) at 3620 cm(-1). In this study, bilayer and monolayer portlandite [Ca(OH)(2)] are predicted to be stable and their characteristics are analyzed in detail. Our results can guide further research on ultrathin hydroxites.