Abstract: The Green function approach to the Bardeen-Cooper-Schrieffer theory of superconductivity is used to study nanofilms. We go beyond previous models and include effects of confinement on the strength of the electron-phonon coupling as well as on the electronic spectrum and on the phonon modes. Within our approach, we find that in ultrathin films, confinement effects on the electronic screening become very important. Indeed, contrary to what has been advanced in recent years, the sudden increases of the density of states when new bands start to be occupied as the film thickness increases, tend to suppress the critical temperature rather than to enhance it. On the other hand, the increase of the number of phonon modes with increasing number of monolayers in the film leads to an increase in the critical temperature. As a consequence, the superconducting critical parameters in such nanofilms are determined by these two competing effects. Furthermore, in sufficiently thin films, the condensate consists of well-defined subcondensates associated with the occupied bands, each with a distinct coherence length. The subcondensates can interfere constructively or destructively giving rise to an interference pattern in the Cooper pair probability density.

Abstract: Using the Ginzburg-Landau (GL) theory, we calculate the stability of sample symmetry-induced vortex-antivortex molecules in a mesoscopic superconducting bilayer exposed to a homogeneous magnetic field. We demonstrate the conditions under which the two condensates cooperatively broaden the field-temperature stability range of the composite (joint) vortex-antivortex state. In cases when such broadening is not achieved, a reentrance of the vortex-antivortex state is found at lower temperatures. In a large portion of the phase diagram noncomposite states are possible, in which the antivortex is present in only one of the layers. In this case, we demonstrate that the vortex-antivortex molecule in one of the layers can be pinned and enlarged by interaction with a vortex molecule in the other. Using analogies in the respective GL formalisms, we map our findings for the bilayer onto mesoscopic two-band superconductors.

Abstract: The conductance, the transmission, and the reflection probabilities through rectangular potential barriers and p-n junctions are obtained for bilayer graphene taking into account the four bands of the energy spectrum. We have evaluated the importance of the skew hopping parameters gamma(3) and gamma(4) to these properties and show that for energies E > gamma(1)/100 their effect is negligible. For high energies two modes of propagation exist and we investigate scattering between these modes. For perpendicular incidence both propagation modes are decoupled, and scattering between them is forbidden. This extends the concept of pseudospin as defined within the two-band approximation to a four-band model and corresponds to the (anti) symmetry of the wave functions under in-plane mirroring. New transmission resonances are found that appear as sharp peaks in the conductance which are absent in the two-band approximation. The application of an interlayer bias to the system (1) breaks the pseudospin structure, (2) opens a band gap that results in a distinct feature of suppressed transmission in the conductance, and (3) breaks the angular symmetry with respect to normal incidence in the transmission and reflection.

Abstract: Density-functional tight-binding and classical molecular dynamics simulations are used to investigate the structural deformations and melting of planar carbon nanoclusters C-N with N = 2-55. The minimum-energy configurations for different clusters are used as starting configurations for the study of the temperature effects on the bond breaking and rotation in carbon lines (N < 6), carbon rings (5 < N < 19), and graphene nanoflakes. The larger the rings (graphene nanoflakes) the higher the transition temperature (melting point) with ring-to-line (perfect-to-defective) transition structures. The melting point was obtained by using the bond energy, the Lindemann criteria, and the specific heat. We found that hydrogen-passivated graphene nanoflakes (CNHM) have a larger melting temperature with a much smaller dependence on size. The edges in the graphene nanoflakes exhibit several different metastable configurations (isomers) during heating before melting occurs. DOI: 10.1103/PhysRevB.87.134103

Abstract: Motivated by recent experimental observations of Tongay et al. [Nano Lett. 12, 5576 (2012)] we show how the electronic properties and Raman characteristics of single layer MoSe2 are affected by elastic biaxial strain. We found that with increasing strain: (1) the E' and E '' Raman peaks (E-2g and E-1g in bulk) exhibit significant redshifts (up to similar to 30 cm(-1)), (2) the position of the A'(1) peak remains at similar to 180 cm(-1) (A(1g) in bulk) and does not change considerably with further strain, (3) the dispersion of low energy flexural phonons crosses over from quadratic to linear, and (4) the electronic band structure undergoes a direct to indirect band gap crossover under similar to 3% biaxial tensile strain. Thus the application of strain appears to be a promising approach for a rapid and reversible tuning of the electronic, vibrational, and optical properties of single layer MoSe2 and similar MX2 dichalcogenides. DOI:10.1103/PhysRevB.87.125415

Abstract: We study Dirac plasmons and their damping in spatially separated N-layer graphene structures at finite doping and temperatures. The plasmon spectrum consists of one optical excitation with square-root dispersion and N – 1 acoustical excitations with linear dispersion, which are undamped at zero temperature and finite doping within a triangular energy region outside the electron-hole continuum. In the long-wavelength limit the energy and weight of the optical plasmon modes increase, respectively, as the square root and linearly with N in agreement with recent experimental findings. The energy and weight of the upper-lying acoustical branches also increase with N. This increase is strongest for the uppermost acoustical mode, and we find that its energy can exceed at some value of momentum the plasmon energy in an individual graphene sheet. Meanwhile, the energy of the low-lying acoustical branches decreases weakly with N as compared with the single acoustical mode in double-layer graphene structures. Our numerical calculations provide a detailed understanding of the overall behavior of the wave-vector dependence of the optical and acoustical multilayer plasmon modes and show how their dispersion and damping are modified as a function of temperature, interlayer spacing, and inlayer carrier density in (un)balanced graphene multilayer structures. DOI: 10.1103/PhysRevB.87.085401

Abstract: The instability of a free-standing one-sided hydrogenated/fluorinated graphene nanoribbon, i.e., graphone/fluorographene, is studied using ab initio, semiempirical, and large-scale molecular dynamics simulations. Free-standing semi-infinite armchairlike hydrogenated/fluorinated graphene (AC-GH/AC-GF) and boatlike hydrogenated/fluorinated graphene (B-GH/B-GF) (nanoribbons which are periodic along the zigzag direction) are unstable and spontaneously transform into spiral structures. We find that rolled, spiral B-GH and B-GF are energetically more favorable than spiral AC-GH and AC-GF which is opposite to the double-sided flat hydrogenated/fluorinated graphene, i.e., graphane/fluorographene. We found that the packed, spiral structures exhibit an unexpected localized highest occupied molecular orbital and lowest occupied molecular orbital at the edges with increasing energy gap during rolling. These rolled hydrocarbon structures are stable beyond room temperature up to at least T = 1000 K within our simulation time of 1 ns. DOI: 10.1103/PhysRevB.87.075448

Abstract: Nonlinear dynamics of Josephson vortices (fluxons) in artificial stacks of superconducting-normal-superconducting Josephson junctions under simultaneously applied time-periodic ac and constant biasing dc currents is studied using the time dependent Ginzburg-Landau formalism with a Lawrence-Doniach extension. At zero external magnetic field and dc biasing current the resistive state of the system is characterized by periodic nucleation and annihilation of fluxon-antifluxon pairs, relative positions of which are determined by the state of neighboring junctions. Due to the mutual repulsive interaction, fluxons in different junctions move out of phase. Their collective motion can be synchronized by adding a small ac component to the biasing dc current. Coherent motion of fluxons is observed for a broad frequency range of the applied drive. In the coherent state the maximal output voltage, which is proportional to the number of junctions in the stack, is observed near the characteristic frequency of the system determined by the crossing of the fluxons across the sample. However, in this frequency range the dynamically synchronized state has an alternative-a less ordered state with smaller amplitude of the output voltage. Collective behavior of the junctions is strongly affected by the sloped sidewalls of the stack. Synchronization is observed only for weakly trapezoidal cross sections, whereas irregular motion of fluxons is observed for larger slopes of the sample edge.

Abstract: Using atomistic simulations we investigate the thermodynamical properties of a single atomic layer of hexagonal boron nitride (h-BN). The thermal induced ripples, heat capacity, and thermal lattice expansion of large scale h-BN sheets are determined and compared to those found for graphene (GE) for temperatures up to 1000 K. By analyzing the mean-square height fluctuations < h(2)> and the height-height correlation function H(q) we found that the h-BN sheet is a less stiff material as compared to graphene. The bending rigidity of h-BN (i) is about 16% smaller than the one of GE at room temperature (300 K), and (ii) increases with temperature as in GE. The difference in stiffness between h-BN and GE results in unequal responses to external uniaxial and shear stress and different buckling transitions. In contrast to a GE sheet, the buckling transition of a h-BN sheet depends strongly on the direction of the applied compression. The molar heat capacity, thermal-expansion coefficient, and Gruneisen parameter are estimated to be 25.2 J mol(-1) K-1, 7.2 x 10(-6) K-1, and 0.89, respectively.

Abstract: The diffusion of a system of ferromagnetic dipoles confined in a quasi-one-dimensional parabolic trap is studied using Brownian dynamics simulations. We show that the dynamics of the system is tunable by an in-plane external homogeneous magnetic field. For a strong applied magnetic field, we find that the mobility of the system, the exponent of diffusion, and the crossover time among different diffusion regimes can be tuned by the orientation of the magnetic field. For weak magnetic fields, the exponent of diffusion in the subdiffusive regime is independent of the orientation of the external field. DOI: 10.1103/PhysRevE.87.012307

Abstract: We modeled the electron and hole states in Si/SiO2 quantum wells within a basis of standing waves using the 30-band k . p theory. The hard-wall confinement potential is assumed, and the influence of the peculiar band structure of bulk silicon on the quantum-well sub-bands is explored. Numerous spurious solutions in the conduction-band and valence-band energy spectra are found and are identified to be of two types: (1) spurious states which have large contributions of the bulk solutions with large wave vectors (the high-k spurious solutions) and (2) states which originate mainly from the spurious valley outside the Brillouin zone (the extravalley spurious solutions). An algorithm to remove all those nonphysical solutions from the electron and hole energy spectra is proposed. Furthermore, slow and oscillatory convergence of the hole energy levels with the number of basis functions is found and is explained by the peculiar band mixing and the confinement in the considered quantum well. We discovered that assuming the hard-wall potential leads to numerical instability of the hole states computation. Nonetheless, allowing the envelope functions to exponentially decay in a barrier of finite height is found to improve the accuracy of the computed hole states.

Abstract: We theoretically demonstrate that the transport inefficiency recently found experimentally for branched-out mesoscopic networks can also be observed in a quantum ring of finite width with an attached central horizontal branch. This is done by investigating the time evolution of an electron wave packet in such a system. Our numerical results show that the conductivity of the ring does not necessary improve if one adds an extra channel. This ensures that there exists a quantum analog of the Braess paradox, originating from quantum scattering and interference.

Abstract: We investigate the effect of Rashba and intrinsic spin-orbit couplings on the electronic properties and spin configurations of Dirac fermions confined in: (i) a flat graphene sheet, (ii) a graphene wire with p-n-p structure, and (iii) a superlattice of graphene wires. The interplay between the spin-orbit interaction mechanisms breaks the electron-hole symmetry and the spin configuration induced by Rashba spin-orbit coupling lacks inversion symmetry in k space. We show that the Rashba spin-orbit interaction doubles the Fabry-Perot resonant modes in the transmission spectrum of a graphene wire and opens new channels for the electron transmission. Moreover, it leads to the appearance of spin split extra Dirac cones in the energy spectrum of a graphene superlattice. It is shown that the spin of the electrons and holes confined in a flat graphene sheet is always perpendicular to their motion while this is not the case for the other nanostructures.

Abstract: The effect of strain in graphene is usually modeled by a pseudomagnetic vector potential which is, however, derived in the limit of small strain. In realistic cases deviations are expected in view of graphene's very high strain tolerance, which can be up to 25%. Here we investigate the pseudomagnetic field generated by a Gaussian bump and we show that it exhibits significant differences with numerical tight-binding results. Furthermore, we calculate the electronic states in the strained region for a hexagon shaped flake with armchair edges. We find that the sixfold symmetry of the wave functions inside the Gaussian bump is directly related to the different effects of strain along the fundamental directions of graphene: zigzag and armchair. Low energy electrons are strongly confined in the armchair directions and are localized on the carbon atoms of a single sublattice.