Abstract: Ab initio calculations within the density-functional theory formalism are performed to investigate the chemical functionalization of a graphene-like monolayer of siliconsilicenewith B, N, Al, or P atoms. The structural, electronic, magnetic, and vibrational properties are reported. The most preferable adsorption sites are found to be valley, bridge, valley and hill sites for B, N, Al, and P adatoms, respectively. All the relaxed systems with adsorbed/substituted atoms exhibit metallic behavior with strongly bonded B, N, Al, and P atoms accompanied by an appreciable electron transfer from silicene to the B, N, and P adatom/substituent. The Al atoms exhibit opposite charge transfer, with n-type doping of silicene and weaker bonding. The adatoms/substituents induce characteristic branches in the phonon spectrum of silicene, which can be probed by Raman measurements. Using molecular dynamics, we found that the systems under study are stable up to at least T=500 K. Our results demonstrate that silicene has a very reactive and functionalizable surface.

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: Neutral excitons in strained axially symmetric In(Ga)As/GaAs quantum dots with a ringlike shape are investigated. Similar to experimental self-assembled quantum rings, the analyzed quantum dots have volcano-like shapes. The continuum mechanical model is employed to determine the strain distribution, and the single-band envelope function approach is adopted to compute the electron states. The hole states are determined by the axially symmetric multiband Luttinger-Kohn Hamiltonian, and the exciton states are obtained from an exact diagonalization. We found that the presence of the inner layer covering the ring opening enhances the excitonic Aharonov-Bohm (AB) oscillations. The reason is that the hole becomes mainly localized in the inner part of the quantum dot due to strain, whereas the electron resides mainly inside the ring-shaped rim. Interestingly, larger AB oscillations are found in the analyzed quantum dot than in a fully opened quantum ring of the same width. Comparison with the unstrained ringlike quantum dot shows that the amplitude of the excitonic Aharonov-Bohm oscillations are almost doubled in the presence of strain. The computed oscillations of the exciton energy levels are comparable in magnitude to the oscillations measured in recent experiments. DOI: 10.1103/PhysRevB.87.085314

Abstract: The adsorption characteristics of alkali, alkaline-earth, and transition metal adatoms on silicene, a graphene-like monolayer structure of silicon are analyzed by means of first-principles calculations. In contrast to graphene, interaction between the metal atoms and the silicene surface is quite strong due to its highly reactive buckled hexagonal structure. In addition to structural properties, we also calculate the electronic band dispersion, net magnetic moment, charge transfer, work function, and dipole moment of the metal adsorbed silicene sheets. Alkali metals, Li, Na, and K, adsorb to hollow sites without any lattice distortion. As a consequence of the significant charge transfer from alkalis to silicene, metalization of silicene takes place. Trends directly related to atomic size, adsorption height, work function, and dipole moment of the silicene/alkali adatom system are also revealed. We found that the adsorption of alkaline-earth metals on silicene is entirely different from their adsorption on graphene. The adsorption of Be, Mg, and Ca turns silicene into a narrow gap semiconductor. Adsorption characteristics of eight transition metals Ti, V, Cr, Mn, Fe, Co, Mo, and W are also investigated. As a result of their partially occupied d orbital, transition metals show diverse structural, electronic, and magnetic properties. Upon the adsorption of transition metals, depending on the adatom type and atomic radius, the system can exhibit metal, half-metal, and semiconducting behavior. For all metal adsorbates, the direction of the charge transfer is from adsorbate to silicene, because of its high surface reactivity. Our results indicate that the reactive crystal structure of silicene provides a rich playground for functionalization at nanoscale. DOI: 10.1103/PhysRevB.87.085423

Abstract: We propose a simple way of probing the number of modes contributing to the channeling in graphene waveguides which are formed by a gauge potential produced by mechanical strain. The energy mode structure for both homogeneous and non-homogeneous strain regimes is carefully studied using the continuum description of the Dirac equation. We found that high strain values privilege negative (instead of positive) group velocities throughout the guidance, sorting the types of modes flowing through it. We also show how the effect of a substrate-induced gap competes against the strain.

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: 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: Typical Raman spectra of transition-metal dichalcogenides (TMDs) display two prominent peaks, E-2g and A(1g), that are well separated from each other. We find that these modes are degenerate in bulk WSe2 yielding one single Raman peak in contrast to other TMDs. As the dimensionality is lowered, the observed peak splits in two. In contrast, our ab initio calculations predict that the degeneracy is retained even for WSe2 monolayers. Interestingly, for minuscule biaxial strain, the degeneracy is preserved, but once the crystal symmetry is broken by a small uniaxial strain, the degeneracy is lifted. Our calculated phonon dispersion for uniaxially strained WSe2 shows a good match to the measured Raman spectrum, which suggests that uniaxial strain exists in WSe2 flakes, possibly induced during the sample preparation and/or as a result of the interaction between WSe2 and the substrate. Furthermore, we find that WSe2 undergoes an indirect-to-direct band-gap transition from bulk to monolayers, which is ubiquitous for semiconducting TMDs. These results not only allow us to understand the vibrational and electronic properties of WSe2, but also point to effects of the interaction between the monolayer TMDs and the substrate on the vibrational and electronic properties. DOI: 10.1103/PhysRevB.87.165409

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: We investigated the magnetic field dependence of the Hall and the bend resistances in the ballistic regime for a single layer graphene Hall bar structure containing a pn-junction. When both regions are n-type the Hall resistance dominates and Hall type of plateaus are formed. These plateaus occur as a consequence of the restriction on the angle imposed by Snell's law allowing only electrons with a certain initial angles to transmit though the potential step. The size of the plateau and its position is determined by the position of the potential interface as well as the value of the applied potential. When the second region is p-type, the bend resistance dominates, which is asymmetric in field due to the presence of snake states. Changing the position of the pn-interface in the Hall bar strongly affects these states and therefore the bend resistance is also changed. Changing the applied potential, we observe that the bend resistance exhibits a peak around the charge-neutrality point (CNP), which is independent of the position of the pn-interface, while the Hall resistance shows a sign reversal when the CNP is crossed, which is in very good agreement with a recent experiment [J. R. Williams and C. M. Marcus, Phys. Rev. Lett. 107, 046602 (2011)].

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 scattering of two-dimensional (2D) massless electrons as presented in graphene in the presence of a random array of circular magnetic flux tubes is investigated. The momentum relaxation time and the Hall factor are obtained using optical theorem techniques for scattering. Electrons with energy close to those of the Landau levels of the flux tubes exhibit resonant scattering and have a long life-time to reside inside the magnetic flux tube. These resonances appear as sharp structures in the Hall factor and the magneto-resistance.

Abstract: Vortex loops are generated by the inhomogeneous stray field of a magnetic dipole on top of a current-carrying mesoscopic superconductor. Cutting and recombination processes unfold under the applied drive, resulting in periodic voltage oscillations across the sample. We show that a direct and detectable consequence of the cutting and recombination of these vortex loops in the present setup is the onset of vortices at surfaces where they were absent prior to the application of the external current. The nonlinear dynamics of vortex loops is studied within the time-dependent Ginzburg-Landau theory to describe the profound three-dimensional features of their time evolution.

Abstract: Charge carriers in single and multilayered graphene systems behave as chiral particles due to the particular lattice symmetry of the crystal. We show that the interplay between the meta-material properties of graphene multilayers and the pseudospinorial properties of the charge carriers result in the occurrence of Klein and anti-Klein tunneling for rhombohedral stacked multilayers. We derive an algebraic formula predicting the angles at which these phenomena occur and support this with numerical calculations for systems up to four layers. We present a decomposition of an arbitrarily stacked multilayer into pseudospin doublets that have the same properties as rhombohedral systems with a lower number of layers. Copyright (C) EPLA, 2013

Abstract: We investigate the link between superconductivity and atomic vibrations in Nb3Sn films with a thickness ranging from 10 to 50 nm. The challenge of measuring the phonon density of states (PDOS) of these films has been tackled by employing the technique of nuclear inelastic scattering by Sn-119 isotopes to reveal the Sn-partial phonon density of states. With the support of ab initio calculations, we evaluate the effect of reduced film thickness on the PDOS. This approach allows us to estimate the changes in superconducting critical temperature T-c induced by phonon confinement, which turned out to be limited to a few tenths of K. The presented method is successful for the Nb3Sn system and paves the way for more systematic studies of the role of phonon confinement in Sn-containing superconductors.