Search

Heterostructures made of transition metal oxides are new tailor-made materials which are attracting much attention. We have constructed a 6-band k.p Hamiltonian and used it within the envelope function method to calculate the subband structure of a variety of LaAlO3/SrTiO3 heterostructures. By use of density functional calculations, we determine the k.p parameters describing the conduction band edge of SrTiO3: the three effective mass parameters, L=0.6104 eV AA^2, M=9.73 eV AA^2, N=-1.616 eV AA^2, the spin orbit splitting Delta_SO=28.5 meV and the low temperature tetragonal distortion energy splitting Delta_T=2.1 meV. For confined systems we find strongly anisotropic non-parabolic subbands. As an application we calculate bands, density of states and magnetic energy levels and compare the results to Shubnikov-de Haas quantum oscillations observed in high magnetic fields. For typical heterostructures we find that electric field strength at the interface of F = 0.1 meV/AA for a carrier density of 7.2 10^{12} cm^-2 results in a subband structure that is similar to experimental results., Comment: 8 pages, 7 figures

Motivated by the observation of ferromagnetism in carbon foams, a massive search for (meta)stable disorder structures of elemental carbon is performed by a generate and test approach. We use the Density Functional based program SIESTA to optimize the structures and calculate the electronic spectra and spin densities. About 1% of the 24000 optimized structures presents magnetic moments, a necessary but not sufficient condition for intrinsic magnetic order. We analyze the results using elements of graph theory. Although the relation between structure and the occurrence of magnetic moments is not yet fully clarified, we give some minimal requirements for this possibility, such as the existence of three-fold coordinated atoms surrounded by four-fold coordinated atoms. We discuss in detail the most promising structures., Comment: 14 pages, 11 figures

We study by molecular dynamics the structural properties of single layer h-BN in comparison to graphene. We show that the Tersoff bond order potential developed for BN by Albe, Moller and Heinig gives a thermally stable hexagonal single layer with a bending constant kappa= 0.54 eV at T = 0. We find that the non-monotonic behaviour of the lattice parameter, the expansion of the interatomic distance and the growth of the bending rigidity with temperature are qualitatively similar to those of graphene. Conversely, the energetics of point defects is extremely different: instead of Stone Wales defects, the two lowest energy defects in h-BN involve either a broken bond or an out of plane displacement of a N atom to form a tetrahedron with three B atoms in the plane. We provide the formation energies and an estimate of the energy barriers., Comment: 13 pages, 7 figures

The understanding of the structural and thermal properties of membranes, low-dimensional flexible systems in a space of higher dimension, is pursued in many fields from string theory to chemistry and biology. The case of a two-dimensional (2D) membrane in three dimensions is the relevant one for dealing with real materials. Traditionally, membranes are primarily discussed in the context of biological membranes and soft matter in general. The complexity of these systems hindered a realistic description of their interatomic structures based on a truly microscopic approach. Therefore, theories of membranes were developed mostly within phenomenological models. From the point of view of statistical mechanics, membranes at finite temperature are systems governed by interacting long-range fluctuations. Graphene, the first truly two-dimensional system consisting of just one layer of carbon atoms, provides a model system for the development of a microscopic description of membranes. In this Account, we review key results in the microscopic theory of structural and thermal properties of graphene and compare them with the predictions of phenomenological theories. The two approaches are in good agreement for the various scaling properties of correlation functions of atomic displacements. However, some other properties, such as the temperature dependence of the bending rigidity, cannot be understood based on phenomenological approaches. We also consider graphene at very high temperature and compare the results with existing models for two-dimensional melting. The melting of graphene presents a different scenario, and we describe that process as the decomposition of the graphene layer into entangled carbon chains., Comment: 20 pages, 9 figures. Misprints are fixed and text style is changed

We study by density functional and large scale tight-binding transport calculations the electronic structure, magnetism and transport properties of the recently proposed graphene ribbons with edges rolled to form nanotubes. Edges with armchair nanotubes present magnetic moments localized either in the tube or the ribbon and metallic or half-metallic character, depending on the symmetry of the junction. These properties have potential for spin valve and spin filter devices with advantages over other proposed systems. Edges with zigzag nanotubes are either metallic or semiconducting without affecting the intrinsic mobility of the ribbon. By varying the type and size of the nanotubes and ribbons offers the possibility to tailor the magnetic and transport properties, making these systems very promising for applications., Comment: 9 pages, 5 figures

We study the effects of temperature and sliding velocity on superlubricity in numerical simulations of the Frenkel-Kontorova model. We show that resonant excitations of the phonons in an incommensurate sliding body lead to an effective friction and to thermal equilibrium with energy distributed over the internal degrees of freedom. For finite temperature, the effective friction can be described well by a viscous damping force, with a damping coefficient that emerges naturally from the microscopic dynamics. This damping coefficient is a non-monotonic function of the sliding velocity which peaks around resonant velocities and increases with temperature. At low velocities, it remains finite and nonzero, indicating the preservation of superlubricity in the zero-velocity limit. Finally, we propose experimental systems in which our results could be verified.

Grain boundaries with dangling bonds (DBGB) in graphene are studied by atomistic Monte Carlo and molecular dynamics simulations in combination with density functional (SIESTA) calculations. The most stable configurations are selected and their structure is analyzed in terms of grain boundary dislocations. It is shown that the grain boundary dislocation with the core consisting of pentagon, octagon and heptagon (5-8-7 defect) is a typical structural element of DBGB with relatively low energies. Electron energy spectrum and magnetic properties of the obtained DBGB are studied by density functional calculations. It is shown that the 5-8-7 defect is magnetic and that its magnetic moment survives after hydrogenation. The effects of hydrogenation and of out of plane deformations on the magnetic properties of DBGB are studied., Comment: 10 pages, 11 figures, 4 tables, the final version accepted in prb

Solids at incommensurate contact display low-friction, 'superlubric', sliding. For graphene flakes on a graphite surface, superlubric sliding is only temporary due to rotation of the flakes from incommensurate to commensurate contact with the substrate. We examine this rotational channel of friction in a prototype geometry of meso- and macroscopic solid lubrication. By molecular dynamics simulations and theoretical arguments we find that two surfaces lubricated by mobile, rotating graphene flakes exhibit stable superlubric sliding as for ideally incommensurate contacts also when they are covered by randomly oriented pinned graphene patches. For commensurate surfaces, we find a low friction state at low temperature where incommensurate states are not destroyed by thermal fluctuations., Comment: 5 pages, 5 figures

Chaotic internal degrees of freedom of a molecule can act as noise and affect the diffusion of the molecule on a substrate. A separation of time scales between the fast internal dynamics and the slow motion of the centre of mass on the substrate makes it possible to directly link chaos to diffusion. We discuss the conditions under which this is possible, and show that in simple atomistic models with pair-wise harmonic potentials, strong chaos can arise through the geometry. Using molecular-dynamics simulations, we demonstrate that a realistic model of benzene is indeed chaotic, and that the internal chaos affects the diffusion on a graphite substrate.

Recent AFM experiments have shown that the low-friction sliding of incommensurate graphite flakes on graphite can be destroyed by torque-induced rotations. Here we theoretically investigate the stability of superlubric sliding against rotations of the flake. We find that the occurrence of superlubric motion critically depends on the physical parameters and on the experimental conditions: particular scan lines, thermal fluctuations and high loading forces can destroy the stability of superlubric orbits. We find that the optimal conditions to achieve superlubric sliding are given by large flakes, low temperature, and low loads, as well as scanning velocities higher than those used in AFM experiments., Comment: 11 pages, 15 figures

The high temperature behaviour of graphene is studied by atomistic simulations based on an accurate interatomic potential for carbon. We find that clustering of Stone-Wales defects and formation of octagons are the first steps in the process of melting which proceeds via the formation of carbon chains. The molten state forms a three-dimensional network of entangled chains rather than a simple liquid. The melting temperature estimated from the two-dimensional Lindemann criterion and from extrapolation of our simulation for different heating rates is about 4900 K., Comment: 12 pages, 6 figures

In practice, physical membranes are exposed to a certain amount of external strain (tension or compression), due to the environment where they are placed. As a result, the behavior of the phonon modes of the membrane is modified. We show that anharmonic effects in stiff two-dimensional membranes are highly suppressed under the application of tension. For this, we consider the anharmonic coupling between bending and stretching modes in the self-consistent screening approximation (SCSA), and compare the obtained height-height correlation function in the SCSA to the corresponding harmonic propagator. The elasticity theory results are compared to atomistic Monte Carlo simulations for a graphene membrane under tension. We find that, while rather high values of strain are needed to avoid anharmonicity in soft membranes, strain fields less than 1% are enough to suppress all the anharmonic effects in stiff membranes, as graphene., Comment: 7 pages, 6 figures

We present the interatomic force constants and phonon dispersions of graphite and graphene from the LCBOPII empirical bond order potential. We find a good agreement with experimental results, particularly in comparison to other bond order potentials. From the flexural mode we determine the bending rigidity of graphene to be 0.69 eV at zero temperature. We discuss the large increase of this constant with temperature and argue that derivation of force constants from experimental values should take this feature into account. We examine also other graphitic systems, including multilayer graphene for which we show that the splitting of the flexural mode can provide a tool for characterization.

The thermodynamic and mechanical properties of graphitic systems are strongly dependent on the shear elastic constant C44. Using state-of-the-art density functional calculations, we provide the first complete determination of their elastic constants and exfoliation energies. We show that stacking misorientations lead to a severe lowering of C44 of at least one order of magnitude. The lower exfoliation energy and the lower C44 (more bending modes) suggest that flakes with random stacking should be easier to exfoliate than the ones with perfect or rhombohedral stacking. We also predict ultralow friction behaviour in turbostratic graphitic systems., Comment: 7 pages, 6 figures

Crystalline membranes at finite temperatures have an anomalous behavior of the bending rigidity that makes them more rigid in the long wavelength limit. This issue is particularly relevant for applications of graphene in nano- and micro-electromechanical systems. We calculate numerically the height-height correlation function $G(q)$ of crystalline two-dimensional membranes, determining the renormalized bending rigidity, in the range of wavevectors $q$ from $10^{-7}$ \AA$^{-1}$ till 10 \AA$^{-1}$ in the self-consistent screening approximation (SCSA). For parameters appropriate to graphene, the calculated correlation function agrees reasonably with the results of atomistic Monte Carlo simulations for this material within the range of $q$ from $10^{-2}$ \AA$^{-1}$ till 1 \AA$^{-1}$. In the limit $q\rightarrow 0$ our data for the exponent $\eta$ of the renormalized bending rigidity $\kappa_R(q)\propto q^{-\eta}$ is compatible with the previously known analytical results for the SCSA $\eta\simeq 0.82$. However, this limit appears to be reached only for $q<10^{-5}$ \AA$^{-1}$ whereas at intermediate $q$ the behavior of $G(q)$ cannot be described by a single exponent., Comment: 5 pages, 4 figures

We minimize a discrete version of the fourth-order curvature based Landau free energy by extending Brakke's Surface Evolver. This model predicts spherical as well as non-spherical shapes with dimples, bumps and ridges to be the energy minimizers. Our results suggest that the buckling and faceting transitions, usually associated with crystalline matter, can also be an intrinsic property of non-crystalline membranes., Comment: 6 pages, 4 figures (LaTeX macros EPJ), accepted for publication in EPJ E

We study the structural and thermodynamic properties of bilayer graphene, a prototype two-layer membrane, by means of Monte Carlo simulations based on the empirical bond order potential LCBOPII. We present the temperature dependence of lattice parameter, bending rigidity and high temperature heat capacity as well as the correlation function of out-of-plane atomic displacements. The thermal expansion coefficient changes sign from negative to positive above $\approx 400$ K, which is lower than previously found for single layer graphene and close to the experimental value of bulk graphite. The bending rigidity is twice as large than for single layer graphene, making the out-of-plane fluctuations smaller. The crossover from correlated to uncorrelated out-of-plane fluctuations of the two carbon planes occurs for wavevectors shorter than $\approx 3$ nm$^{-1}$, Comment: 6 pages, 7 figures.

We investigate the stability of low-friction sliding of nanocrystal with rectangular atomic arrangement on rectangular lattices, for which analytical results can be obtained. We find that several incommensurate periodic orbits exist and are stable against thermal fluctuations and other perturbations. As incommensurate orientations lead to low corrugation, and therefore low friction, such incommensurate periodic orbits are interesting for the study of nanotribology. The analytical results compare very well with simulations of W nanocrystals on NaF(001). The geometry and high typical corrugation of substrates with square lattices increase the robustness compared to typical hexagonal lattices, such as graphite.

We propose new semi-implicit numerical methods for the integration of the stochastic Landau-Lifshitz equation with built-in angular momentum conservation. The performance of the proposed integrators is tested on the 1D Heisenberg chain. For this system, our schemes show better stability properties and allow us to use considerably larger time steps than standard explicit methods. At the same time, these semi-implicit schemes are also of comparable accuracy to and computationally much cheaper than the standard midpoint implicit method. The results are of key importance for atomistic spin dynamics simulations and the study of spin dynamics beyond the macro spin approximation., Comment: 24 pages, 5 figures

Single-layer graphene sheets are typically characterized by long-wavelength corrugations (ripples) which can be shown to be at the origin of rather strong potentials with both scalar and vector components. We present an extensive microscopic study, based on a self-consistent Kohn-Sham-Dirac density-functional method, of the carrier density distribution in the presence of these ripple-induced external fields. We find that spatial density fluctuations are essentially controlled by the scalar component, especially in nearly-neutral graphene sheets, and that in-plane atomic displacements are as important as out-of-plane ones. The latter fact is at the origin of a complicated spatial distribution of electron-hole puddles which has no evident correlation with the out-of-plane topographic corrugations. In the range of parameters we have explored, exchange and correlation contributions to the Kohn-Sham potential seem to play a minor role., Comment: 13 pages, 13 figures, submitted. High-quality figures can be requested to the authors

We show that the stacking of flat aromatic molecules on a curved surface results in topological defects. We consider, as an example, spherical vesicles, self-assembled from molecules with 5- and 6-thiophene cores. We predict that the symmetry of the molecules influences the number of topological defects and the resulting equilibrium shape., Comment: accepted as a Letter in the J. Phys. Chem. B

Structure and thermodynamics of crystalline membranes are characterized by the long wavelength behavior of the normal-normal correlation function G(q). We calculate G(q) by Monte Carlo and Molecular Dynamics simulations for a quasi-harmonic model potential and for a realistic potential for graphene. To access the long wavelength limit for finite-size systems (up to 40000 atoms) we introduce a Monte Carlo sampling based on collective atomic moves (wave moves). We find a power-law behaviour $G(q)\propto q^{-2+\eta}$ with the same exponent $\eta \approx 0.85$ for both potentials. This finding supports, from the microscopic side, the adequacy of the scaling theory of membranes in the continuum medium approach, even for an extremely rigid material like graphene.