Abstract: Using first-principles calculations and reactive force field molecular dynamics simulations, we study the structural properties and dynamics of a fluorine (F) atom, either adsorbed on the surface of single layer graphene (F/GE) or between the layers of AB stacked bilayer graphene (F@ bilayer graphene). It is found that the diffusion of the F atom is very different in those cases, and that the mobility of the F atom increases by about an order of magnitude when inserted between two graphene layers. The obtained diffusion constant for F/GE is twice larger than that experimentally found for gold adatom and theoretically found for C-60 molecule on graphene. Our study provides important physical insights into the dynamics of fluorine atoms between and on graphene layers and explains the mechanism behind the separation of graphite layers due to intercalation of F atoms.

Abstract: Anharmonic effects in an atomic monolayer thin crystal with honeycomb lattice structure are investigated by analytical and numerical lattice dynamical methods. Starting from a semiempirical model for anharmonic couplings of third and fourth orders, we study the in-plane and out-of-plane (flexural) mode components of the generalized wave vector dependent Gruneisen parameters, the thermal tension and the thermal expansion coefficients as a function of temperature and crystal size. From the resonances of the displacement-displacement correlation functions, we obtain the renormalization and decay rate of in-plane and flexural phonons as a function of temperature, wave vector, and crystal size in the classical and in the quantum regime. Quantitative results are presented for graphene. There, we find that the transition temperature T-alpha from negative to positive thermal expansion is lowered with smaller system size. Renormalization of the flexural mode has the opposite effect and leads to values of T-alpha approximate to 300 K for systems of macroscopic size. Extensive numerical analysis throughout the Brillouin zone explores various decay and scattering channels. The relative importance of normal and umklapp processes is investigated. The work is complementary to crystalline membrane theory and computational studies of anharmonic effects in two-dimensional crystals.

Abstract: A two-dimensional (2D) binary colloidal system consisting of interacting dipoles is investigated using an analytical approach. Within the harmonic approximation we obtain the phonon spectrum of the system as a function of the composition, dipole-moment ratio, and mass ratio between the small and big particles. Through a systematic analysis of the phonon spectra we are able to determine the stability region of the different lattice structures of the colloidal alloys. The gaps in the phonon frequency spectrum, the optical frequencies in the long-wavelength limit, and the sound velocity are discussed as well. Using the modified Lindemann criterion and within the harmonic approximation we estimate the melting temperature of the sublattice generated by the big particles.

Abstract: We investigate the dependence of the structural phase transitions in an infinite quasi-one-dimensional system of repulsively interacting particles on the profile of the confining channel. Three different functional expressions for the confinement potential related to real experimental systems are used that can be tuned continuously from a parabolic to a hard-wall potential in order to find a thorough understanding of the ordering of the chainlike structure transitions. We resolve the long-standing issue why the most theories predicted a 1-2-4-3-4 sequence of chain configurations with increasing density, while some experiments found the 1-2-3-4 sequence.

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: 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: Using first principles calculations we showed that the electronic and optical properties of single-layer black phosphorus (BP) depend strongly on the applied strain. Due to the strong anisotropic atomic structure of BP, its electronic conductivity and optical response are sensitive to the magnitude and the orientation of the applied strain. We found that the inclusion of many body effects is essential for the correct description of the electronic properties of monolayer BP; for example, while the electronic gap of strainless BP is found to be 0.90 eV by using semilocal functionals, it becomes 2.31 eV when many-body effects are taken into account within the G(0)W(0) scheme. Applied tensile strain was shown to significantly enhance electron transport along zigzag direction of BP. Furthermore, biaxial strain is able to tune the optical band gap of monolayer BP from 0.38 eV (at -8% strain) to 2.07 eV (at 5.5%). The exciton binding energy is also sensitive to the magnitude of the applied strain. It is found to be 0.40 eV for compressive biaxial strain of -8%, and it becomes 0.83 eV for tensile strain of 4%. Our calculations demonstrate that the optical response of BP can be significantly tuned using strain engineering which appears as a promising way to design novel photovoltaic devices that capture a broad range of solar spectrum.

Abstract: We study the energy and temperature relaxation of electrons in graphene on a piezoelectric substrate. Scattering from the combined potential of extrinsic piezoelectric surface acoustical (PA) phonons of the substrate and intrinsic deformation acoustical phonons of graphene is considered for a (non) degenerate gas of Dirac fermions. It is shown that in the regime of low energies or temperatures the PA phonons dominate the relaxation and change qualitatively its character. This prediction is relevant for quantum metrology and electronic applications using graphene devices and suggests an experimental setup for probing electron-phonon coupling in graphene.

Abstract: We present a systematic study of the energy spectra of graphene quantum rings having different geometries and edge types in the presence of a perpendicular magnetic field. Results are obtained within the tight-binding (TB) and Dirac models and we discuss which features of the former can be recovered by using the approximations imposed by the latter. Energy levels of graphene quantum rings obtained by diagonalizing the TB Hamiltonian are demonstrated to be strongly dependent on the rings geometry and the microscopical structure of the edges. This makes it difficult to recover those spectra by the existing theories that are based on the continuum (Dirac) model. Nevertheless, our results show that both approaches (i.e., TB and Dirac model) may provide similar results, but only for very specific combinations of ring geometry and edge types. The results obtained by a simplified model describing an infinitely thin circular Dirac ring show good agreement with those obtained for hexagonal and rhombus armchair graphene rings within the TB model. Moreover, we show that the energy levels of a circular quantum ring with an infinite mass boundary condition obtained within the Dirac model agree with those for a ring defined by a ring-shaped staggered potential obtained within the TB model.

Abstract: We present a theoretical study of the plasmon and coupled plasmon-phonon modes induced by intraband electron-electron interaction in graphene in the presence of driving dc electric field. We find that the electric field dependence of these collective excitation modes in graphene differs significantly from that in a conventional two-dimensional electron gas with a parabolic energy spectrum. This is due mainly to the fact that graphene has a linear energy spectrum and the Fermi velocity of electrons in graphene is much larger than the drift velocity of electrons. The obtained results demonstrate that the plasmon and coupled plasmon-phonon modes in graphene can be tuned by applying not only the gate voltage but also the source-to-drain field. The manipulation of plasmon and coupled plasmon-phonon modes by source-to-drain voltage can let graphene be more conveniently applied as an advanced plasmonic material.