Abstract: Unidirectional chiral states are predicted in single layer graphene which originate from the breaking of the sublattice symmetry due to an asymmetric mass potential. The latter can be created experimentally using boron-nitride (BN) substrates with a line defect (B-B or N-N) that changes the induced mass potential in graphene. Solving the Dirac-Weyl equation, the obtained energy spectrum is compared with the one calculated using ab initio density functional calculations. We found that these one-dimensional chiral states are very robust and they can even exist in the presence of a small gap between the mass regions. In the latter case additional bound states are found that are topologically different from those chiral states.

Abstract: The interplay between quantum Hall states and Cooper pairs is usually hindered by the suppression of the superconducting state due to the strong magnetic fields needed to observe the quantum Hall effect. From this point of view, graphene is special since it allows the creation of strong pseudomagnetic fields due to strain. We show that in a Josephson junction made of strained graphene, Cooper pairs will diffuse into the strained region. The pair correlation function will be sublattice polarized due to the polarization of the local density of states in the zero pseudo-Landau level. We uncover two regimes: (1) one in which the cyclotron radius is larger than the junction length, in which case the supercurrent will be enhanced, and (2) the long junction regime where the supercurrent is strongly suppressed because the junction becomes an insulator. In the latter case quantized Hall states form and Andreev scattering at the normal/superconducting interface will induce edge states. Our numerical calculation has become possible due to an extension of the Chebyshev-Bogoliubovde Gennes method to computations on video cards (GPUs).

Abstract: To date, several experimental groups reported measurements of the thickness dependence of T-c of atomically uniform single-crystalline Pb nanofilms. The reported amplitude of the T-c oscillations varies significantly from one experiment to another. Here we propose that the reason for this unresolved issue is an interplay of the quantum-size variations in the single-electron density of states with thickness-dependent oscillations in the phonon-mediated electron-electron coupling. Such oscillations in the coupling depend on the substrate material, the quality of the interface, the protection cover, and other details of the fabrication process, changing from one experiment to another. This explains why the available data do not exhibit one-voice consistency about the amplitude of the T-c oscillations. Our analyses are based on a numerical solution of the Bogoliubov-de Gennes equations for a superconducting slab.

Abstract: The structural phase transitions with electric-quadrupole long-range order in NpO2 (Fm (3) over barm-->Pn (3) over barm) and UO2 (Fm (3) over barm-->Pa (3) over bar) are analyzed from a group theoretical point of view. In both cases, the symmetry lowering involves three quadrupolar components belonging to the irreducible representation T-2g (Gamma(5)) of O-h and condensing in a triple-q structure at the X point of the Brillouin zone. The Pa (3) over bar structure is close to Pn (3) over barm, but allows for oxygen displacements. The Pa (3) over bar ordering leads to an effective electrostatic attraction between electronic quadrupoles while the Pn (3) over barm ordering results in a repulsion between them. It is concluded that the Pn (3) over barm structure can be stabilized only through some additional process such as strengthening of the chemical bonding between Np and O. We also derive the relevant structure-factor amplitudes for Pn (3) over barm and Pa (3) over bar, and the effect of domains on resonant x-ray scattering experiments.

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: The lattice dynamics of low-frequency rigid-plane modes in metallic (graphene multilayers, GML) and in insulating (hexagonal boron-nitride multilayers, BNML) layered crystals is investigated. The frequencies of shearing and compression (stretching) modes depend on the layer number N and are presented in the form of fan diagrams. The results for GML and BNML are very similar. In both cases, only the interactions (van der Waals and Coulomb) between nearest-neighbor planes are effective, while the interactions between more distant planes are screened. A comparison with recent Raman scattering results on low-frequency shear modes in GML [Tan et al., Nat. Mater., in press, doi: 10.1038/nmat3245, (2012)] is made. Relations with the low-lying rigid-plane phonon dispersions in the bulk materials are established. Master curves, which connect the fan diagram frequencies for any given N, are derived. Static and dynamic thermal correlation functions for rigid-layer shear and compression modes are calculated. The results might be of use for the interpretation of friction force experiments on multilayer crystals.

Abstract: Large-scale atomistic simulations using the reactive force field approach are implemented to investigate the thermomechanical properties of fluorinated graphene (FG). A set of parameters for the reactive force field potential optimized to reproduce key quantum mechanical properties of relevant carbon-fluorine cluster systems are presented. Molecular dynamics simulations are used to investigate the thermal rippling behavior of FG and its mechanical properties and compare them with graphene, graphane and a sheet of boron nitride. The mean square value of the height fluctuations < h(2)> and the height-height correlation function H(q) for different system sizes and temperatures show that FG is an unrippled system in contrast to the thermal rippling behavior of graphene. The effective Young's modulus of a flake of fluorinated graphene is obtained to be 273 N/m and 250 N/m for a flake of FG under uniaxial strain along armchair and zigzag directions, respectively. DOI: 10.1103/PhysRevB.87.104114

Abstract: The thermodynamic properties of the electron gas in multilayer graphene depend strongly on the number of layers and the type of stacking. Here we analyze how those properties change when we vary the number of layers for rhombohedral stacked multilayer graphene and compare our results with those from a conventional two-dimensional electron gas. We show that the highly degenerate zero-energy Landau level which is partly filled with electrons and partly with holes has a strong influence on the values of the different thermodynamic quantities.

Abstract: Using highly efficient simulations of the tight-binding Bogoliubov-de-Gennes model, we solved self-consistently for the pair correlation and the Josephson current in a superconducting-bilayer graphene-superconducting Josephson junction. Different doping levels for the non-superconducting link are considered in the short- and long-junction regimes. Self-consistent results for the pair correlation and superconducting current resemble those reported previously for single-layer graphene except at the Dirac point, where remarkable differences in the proximity effect are found, as well as a suppression of the superconducting current in the long-junction regime. Inversion symmetry is broken by considering a potential difference between the layers and we found that the supercurrent can be switched if the junction length is larger than the Fermi length.