Abstract: Fluorinated graphene is one of the most important derivatives of graphene and has been found to have great potential in optoelectronic and photonic nanodevices. However, the stability of F atoms on fluorinated graphene under different conditions, which is essential to maintain the desired properties of fluorinated graphene, is still unclear. In this work, we investigate the diffusion of F atoms on pristine graphene, graphene with defects, and at graphene/fluorographene interfaces by using density functional theory calculations. We find that an isolated F atom diffuses easily on graphene, but those F atoms can be localized by inducing vacancies or absorbates in graphene and by creating graphene/fluorographene interfaces, which would strengthen the binding energy of F atoms on graphene and increase the diffusion energy barrier of F atoms remarkably.

Abstract: We investigate the diffusive electron-transport properties of charge-doped graphene ribbons and nanoribbons with imperfect edges. We consider different regimes of edge scattering, ranging from wide graphene ribbons with (partially) diffusive edge scattering to ribbons with large width variations and nanoribbons with atomistic edge roughness. For the latter, we introduce an approach based on pseudopotentials, allowing for an atomistic treatment of the band structure and the scattering potential, on the self-consistent solution of the Boltzmann transport equation within the relaxation-time approximation and taking into account the edge-roughness properties and statistics. The resulting resistivity depends strongly on the ribbon orientation, with zigzag (armchair) ribbons showing the smallest (largest) resistivity and intermediate ribbon orientations exhibiting intermediate resistivity values. The results also show clear resistivity peaks, corresponding to peaks in the density of states due to the confinement-induced subband quantization, except for armchair-edge ribbons that show a very strong width dependence because of their claromatic behavior. Furthermore, we identify a strong interplay between the relative position of the two valleys of graphene along the transport direction, the correlation profile of the atomistic edge roughness, and the chiral valley modes, leading to a peculiar strongly suppressed resistivity regime, most pronounced for the zigzag orientation.

Abstract: The thermal signatures of superconductivity in mesoscopic disks, rings and cylinders are calculated within the Ginzburg-Landau theory. In an applied perpendicular magnetic field H the heat capacity of mesoscopic samples shows a strong dependence on the realized vortex state; discontinuities are found at the critical field for different vorticities, as well as at the superconducting-to-normal state transition. The same applies to the intermediate state of type-I superconductors. Even the subtle changes in the fluxoid distribution inside the sample leave clear signatures on heat capacity, which is particularly useful for fully three-dimensional samples whose interior is often inaccessible by magnetometry. The heat-capacity jump ΔC(H) at the critical temperature exhibits quasiperiodic modulations as a function of magnetic field. In mesoscopic superconducting rings, these oscillations provide calorimetric verification of the Little-Parks effect.

Abstract: Nanoindentation of bilayer graphene is studied using molecular-dynamics simulations. We compared our simulation results with those from elasticity theory as based on the nonlinear Föppl-Hencky equations with rigid boundary condition. The force-deflection values of bilayer graphene are compared to those of monolayer graphene. Youngs modulus of bilayer graphene is estimated to be 0.8 TPa which is close to the value for graphite. Moreover, an almost flat bilayer membrane at low temperature under central load has a 14% smaller Youngs modulus as compared to the one at room temperature.

Abstract: Bi2Se3 is a three-dimensional topological insulator which has been extensively studied because it has a single Dirac cone on the surface, inside a relatively large bulk band gap. However, the effect of two-dimensional topological insulator Bi bilayers on the properties of Bi2Se3 and vice versa, has not been explored much. Bi bilayers are often present between the quintuple layers of Bi2Se3, since (Bi2)n(Bi2Se3)m form stable ground-state structures. Moreover, Bi2Se3 is a good substrate for growing ultrathin Bi bilayers. By first-principles techniques, we first show that there is no preferable surface termination by either Bi or Se. Next, we investigate the electronic structure of Bi bilayers on top of, or inside a Bi2Se3 slab. If the Bi bilayers are on top, we observe a charge transfer to the quintuple layers that increases the binding energy of the surface Dirac cones. The extra states, originating from the Bi bilayers, were declared to form a topological Dirac cone, but here we show that these are ordinary Rashba-split states. This result, together with the appearance of a new Dirac cone that is localized slightly deeper, might necessitate the reinterpretation of several experimental results. When the Bi bilayers are located inside the Bi2Se3 slab, they tend to split the slab into two topological insulators with clear surface states. Interface states can also be observed, but an energy gap persists because of strong coupling between the neighboring quintuple layers and the Bi bilayers.

Abstract: We study by magnetoresonant tunneling spectroscopy single-particle energy spectra of the constituent weakly coupled dots in vertical quantum dot molecules over a wide energy window. The measured energy spectra are well modeled by calculated spectra for dots with in-plane confinement potentials that are elliptical and parabolic in form. However, in the regions where two, three, or four single-particle energy levels are naively expected to cross, we observe pronounced level anticrossing behavior and strong variations in the resonant currents as a consequence of coherent mixing induced by small deviations in the nearly ideal dot confinement potentials. We present detailed analysis of the energy spectra, and focus on two examples of three-level crossings whereby the coherent mixing leads to concurrent suppression and enhancement of the resonant currents when the anticrossing levels are minimally separated. The suppression of resonant current is of particular interest since it is a signature of dark state formation due to destructive interference. We also describe in detail and compare two measurement strategies to reliably extract the resonant currents required to characterize the level mixing.

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: The fine structure of the Dirac energy spectrum in graphene induced by electron-optical phonon coupling is investigated in the portion of the spectrum near the phonon emission threshold. The derived new dispersion equation in the immediate neighborhood below the threshold corresponds to an electron-phonon bound state. We find that the singular vertex corrections beyond perturbation theory strongly increase the electron-phonon binding energy scale. The predicted enhancement of the effective electron-phonon coupling can be measured using angle-resolved spectroscopy.

Abstract: <script type='text/javascript'>document.write(unpmarked('Quantum electrodynamics predicts that heavy atoms (Z \u003E Z(c) approximate to 170) will undergo the process of atomic collapse where electrons sink into the positron continuum and a new family of so-called collapsing states emerges. The relativistic electrons in graphene exhibit the same physics but at a much lower critical charge (Z(c) approximate to 1) which has made it possible to confirm this phenomenon experimentally. However, there exist conflicting predictions on the effect of a magnetic field on atomic collapse. These theoretical predictions are based on the continuum Dirac-Weyl equation, which does not have an exact analytical solution for the interplay of a supercritical Coulomb potential and the magnetic field. Approximative solutions have been proposed, but because the two effects compete on similar energy scales, the theoretical treatment varies depending on the regime which is being considered. These limitations are overcome here by starting from a tight-binding approach and computing exact numerical results. By avoiding special limit cases, we found a smooth evolution between the different regimes. We predict that the atomic collapse effect persists even after the magnetic field is activated and that the critical charge remains unchanged. We show that the atomic collapse regime is characterized: (1) by a series of Landau level anticrossings and (2) by the absence of root B scaling of the Landau levels with regard to magnetic field strength.'));

Abstract: We have calculated the plasmon modes in graphene double layer structures at finite temperatures, taking into account the inhomogeneity of the dielectric background of the system. The effective dielectric function is obtained from the solution of the Poisson equation of a three-layer dielectric medium with graphene sheets located at the interfaces, separating the different materials. Due to the momentum dispersion of the effective dielectric function, the intra- and interlayer bare Coulomb interactions in the graphene double layer system acquires an additional momentum dependence-an effect that is of the order of the interlayer interaction itself. We show that the energies of the in-phase and out-of-phase plasmon modes are determined largely by different values of the spatially dependent effective dielectric function. The effect of the dielectric inhomogeneity increases with temperature, and even at high temperatures the energy shift induced by the dielectric inhomogeneity and temperature itself remains larger than the broadening of the plasmon energy dispersions due to the Landau damping. The obtained new features of the plasmon dispersions can be observed in frictional drag measurements and in inelastic light scattering and electron energy-loss spectroscopies.

Abstract: Using atomistic simulations we investigate the morphological properties of graphene deposited on top of a nanostructured substrate. Sinusoidally corrugated surfaces, steps, elongated trenches, one-dimensional and cubic barriers, spherical bubbles, Gaussian bumps, and Gaussian depressions are considered as support structures for graphene. The graphene-substrate interaction is governed by van der Waals forces and the profile of the graphene layer is determined by minimizing the energy using molecular dynamics simulations. Based on the obtained optimum configurations, we found that (i) for graphene placed over sinusoidally corrugated substrates with corrugation wavelengths longer than 2 nm, the graphene sheet follows the substrate pattern while for supported graphene it is always suspended across the peaks of the substrate, (ii) the conformation of graphene to the substrate topography is enhanced when increasing the energy parameter in the van der Waals model, (iii) the adhesion of graphene into the trenches depends on the width of the trench and on the graphene's orientation, i. e., in contrast to a small-width (3 nm) nanoribbon with armchair edges, the one with zigzag edges follows the substrate profile, (iv) atomic-scale graphene follows a Gaussian bump substrate but not the substrate with a Gaussian depression, and (v) the adhesion energy due to van der Waals interaction varies in the range [0.1-0.4] J/m(2).

Abstract: Rashba-type spin-orbit interaction (SOI) controlled band structure of a two-dimensional superlattice allows for the modulation of the conductance of finite size devices by changing the strength of the SOI. We consider rectangular arrays and find that the temperature dependence of the conductance disappears for high temperatures, but the strength of the SOI still affects the conductance at these temperatures. The modulation effect can be seen even in the presence of strong dephasing, which can be important for practical applications.

Abstract: Atomic collapse is a phenomenon inherent to relativistic quantum mechanics where electron states dive in the positron continuum for highly charged nuclei. This phenomenon was recently observed in graphene. Here we investigate a novel collapse phenomenon when multiple sub- and supercritical charges of equal strength are put close together as in a molecule. We construct a phase diagram which consists of three distinct regions: (1) subcritical, (2) frustrated atomic collapse, and (3) molecular collapse. We show that the single impurity atomic collapse resonances rearrange themselves to form molecular collapse resonances which exhibit a distinct bonding, anti-bonding and non-bonding character. Here we limit ourselves to systems consisting of two and three charges. We show that by tuning the distance between the charges and their strength a high degree of control over the molecular collapse resonances can be achieved.

Abstract: We study the dynamics of vortices in an asymmetric (i.e., consisting of triangular cells) ring channel driven by an external ac current I in a Corbino setup. The asymmetric potential rectifies the motion of vortices and induces a net vortex flow without any unbiased external drive, i.e., the ratchet effect. We show that the net flow of vortices strongly depends on vortex density and frequency of the driving current. Depending on the density, we distinguish a single-vortex rectification regime (for low density, when each vortex is rectified individually) determined by the potential-energy landscape inside each cell of the channel (i.e., hard and easy directions) and multi-vortex, or collective, rectification (high-density case) when the inter-vortex interaction becomes important. We analyze the average angular velocity ω of vortices as a function of I and study commensurability effects between the numbers of vortices and cells in the channel and the role of frequency of the applied ac current. We have shown that the commensurability effect results in a stepwise ω-I curve. Besides the integer steps, i.e., the large steps found in the single-vortex case, we also found fractional steps corresponding to fractional ratios between the numbers of vortices and triangular cells. We have performed preliminary measurements on a device containing a single weak-pinning circular ratchet channel in a Corbino geometry and observed a substantial asymmetric vortex response.

Abstract: The finding that triggering the redox activity of oxygen ions within the lattice of transition metal oxides can boost the performances of materials used in energy storage and conversion devices such as Li-ion batteries or oxygen evolution electrocatalysts has recently spurred intensive and innovative research in the field of energy. While experimental and theoretical efforts have been critical in understanding the role of oxygen nonbonding states in the redox activity of oxygen ions, a clear picture of the redox chemistry of the oxygen species formed upon this oxidation process is still missing. This can be, in part, explained by the complexity in stabilizing and studying these species once electrochemically formed. In this work, we alleviate this difficulty by studying the phase Ba5Ru2O11, which contains peroxide O-2(2-) groups, as oxygen evolution reaction electrocatalyst and Li-ion battery material. Combining physical characterization and electrochemical measurements, we demonstrate that peroxide groups can easily be oxidized at relatively low potential, leading to the formation of gaseous dioxygen and to the instability of the oxide. Furthermore, we demonstrate that, owing to the stabilization at high energy of peroxide, the high-lying energy of the empty sigma* antibonding O-O states limits the reversibility of the electrochemical reactions when the O-2(2-)/O2- redox couple is used as redox center for Li-ion battery materials or as OER redox active sites. Overall, this work suggests that the formation of true peroxide O-2(2-) states are detrimental for transition metal oxides used as OER catalysts and Li-ion battery materials. Rather, oxygen species with O-O bond order lower than 1 would be preferred for these applications.