Showing 1–50 of 72 results for author: Fasolino, A

Alternatives for Testing of Context-Aware Contemporary Software Systems in industrial settings: Results from a Rapid review

Authors: Santiago Matalonga, Domenico Amalfitano, Andrea Doreste, Ana Rita Fasolino, Guilherme H. Travassos

Abstract: Context: Context-aware contemporary software systems (CACSS) are mainstream. Furthermore, they present challenges for current engineering practices. These challenges are distinctively present when testing CACSS, as the variation of context deepens the limitations of available software testing practices and technologies. Objective: To understand how the industry deals with the variation of context… ▽ More Context: Context-aware contemporary software systems (CACSS) are mainstream. Furthermore, they present challenges for current engineering practices. These challenges are distinctively present when testing CACSS, as the variation of context deepens the limitations of available software testing practices and technologies. Objective: To understand how the industry deals with the variation of context when testing CACSS. Method: A Rapid Review was commissioned to uncover the necessary evidence to achieve the objectives. Results: Our results show that current research initiatives aim to generate or improve Test Suites that can deal with the variation of context and the sheer volume of test input possibilities. To achieve this, they mostly rely on modelling the systems' dynamic behavior and increasing computing resources to generate test inputs. We found no evidence of research results aiming at managing context variation through the testing lifecycle process. Conclusions: We discuss how the identified solutions are not ready for mainstream adoption. They are all domain-specific, and while the ideas and approaches can be reproduced in different settings, the technologies noon to be re-engineered and tailor to the specific CACSS. △ Less

Submitted 14 April, 2021; v1 submitted 3 April, 2021; originally announced April 2021.

Comments: Submitted to JSS for peer review

Graphene on h-BN: to align or not to align?

Authors: Roberto Guerra, Merel van Wijk, Andrea Vanossi, Annalisa Fasolino, Erio Tosatti

Abstract: The contact strength, adhesion and friction, between graphene and an incommensurate crystalline substrate such as {\it h}-BN depends on their relative alignment angle $θ$. The well established Novaco-McTague (NM) theory predicts for a monolayer graphene on a hard bulk {\it h}-BN crystal face a small spontaneous misalignment, here $θ_{NM}$\,$\simeq$\,0.45 degrees which if realized would be relevant… ▽ More The contact strength, adhesion and friction, between graphene and an incommensurate crystalline substrate such as {\it h}-BN depends on their relative alignment angle $θ$. The well established Novaco-McTague (NM) theory predicts for a monolayer graphene on a hard bulk {\it h}-BN crystal face a small spontaneous misalignment, here $θ_{NM}$\,$\simeq$\,0.45 degrees which if realized would be relevant to a host of electronic properties besides the mechanical ones. Because experimental equilibrium is hard to achieve, we inquire theoretically about alignment or misalignment by simulations based on dependable state-of-the-art interatomic force fields. Surprisingly at first, we find compelling evidence for $θ= 0$, i.e., full energy-driven alignment in the equilibrium state of graphene on {\it h}-BN. Two factors drive this deviation from NM theory. First, graphene is not flat, developing on {\it h}-BN a long-wavelength out-of-plane corrugation. Second, {\it h}-BN is not hard, releasing its contact stress by planar contractions/expansions that accompany the interface moiré structure. Repeated simulations by artificially forcing graphene to keep flat, and {\it h}-BN to keep rigid, indeed yield an equilibrium misalignment similar to $θ_{NM}$ as expected. Subsequent sliding simulations show that friction of graphene on {\it h}-BN, small and essentially independent of misalignments in the artificial frozen state, strongly increases in the more realistic corrugated, strain-modulated, aligned state. △ Less

Submitted 26 May, 2017; originally announced May 2017.

Theoretical models of Rashba spin splitting in asymmetric SrTiO3-based heterostructures

Authors: L. W. van Heeringen, A. McCollam, G. A. de Wijs, A. Fasolino

Abstract: Rashba spin splitting in two-dimensional (2D) semiconductor systems is generally calculated in a ${\bf k} \cdot {\bf p}$ Luttinger-Kohn approach where the spin splitting due to asymmetry emerges naturally from the bulk band structure. In recent years, several new classes of 2D systems have been discovered where electronic correlations are believed to have an important role. In these correlated sys… ▽ More Rashba spin splitting in two-dimensional (2D) semiconductor systems is generally calculated in a ${\bf k} \cdot {\bf p}$ Luttinger-Kohn approach where the spin splitting due to asymmetry emerges naturally from the bulk band structure. In recent years, several new classes of 2D systems have been discovered where electronic correlations are believed to have an important role. In these correlated systems, the effects of asymmetry leading to Rashba splitting have typically been treated phenomenologically. We compare these two approaches for the case of 2D electron systems in SrTiO$_3$-based heterostructures, and find that the two models produce fundamentally different behavior in regions of the Brillouin zone that are particularly relevant for magnetotransport. Our results demonstrate the importance of identifying the correct approach in the quantitative interpretation of experimental data, and are likely to be relevant to a range of 2D systems in correlated materials. △ Less

Submitted 3 April, 2017; originally announced April 2017.

Comments: Paper (4 figures) and supplementary material (2 figures)

Mechanics of thermally fluctuating membranes

Authors: J. H. Los, A. Fasolino, M. I. Katsnelson

Abstract: Besides having unique electronic properties, graphene is claimed to be the strongest material in nature. In the press release of the Nobel committee it is claimed that a hammock made of a squared meter of one-atom thick graphene could sustain the wight of a 4 kg cat. More practically important are applications of graphene like scaffolds and sensors which are crucially dependent on the mechanical s… ▽ More Besides having unique electronic properties, graphene is claimed to be the strongest material in nature. In the press release of the Nobel committee it is claimed that a hammock made of a squared meter of one-atom thick graphene could sustain the wight of a 4 kg cat. More practically important are applications of graphene like scaffolds and sensors which are crucially dependent on the mechanical strength. Meter-sized graphene is even being considered for the lightsails in the starshot project to reach the star alpha centaury. The predicted strength of graphene is based on its very large Young modulus which is, per atomic layer, much larger than that of steel. This reasoning however would apply to conventional thin plates but does not take into account the peculiar properties of graphene as a thermally fluctuating crystalline membrane. It was shown recently both experimentally and theoretically that thermal fluctuations lead to a dramatic reduction of the Young modulus and increase of the bending rigidity for micron-sized graphene samples in comparison with atomic scale values. This makes the use of the standard Föppl-von Karman elasticity (FvK) theory for thin plates not directly applicable to graphene and other single atomic layer membranes. This fact is important because the current interpretation of experimental results is based on the FvK theory. In particular, we show that the FvK-derived Schwerin equation, routinely used to derive the Young modulus from indentation experiments has to be essentially modified for graphene at room temperature and for micron sized samples. Based on scaling analysis and atomistic simulation we investigate the mechanics of graphene under transverse load up to breaking. We determine the limits of applicability of the FvK theory and provide quantitative estimates for the different regimes. △ Less

Submitted 24 March, 2017; originally announced March 2017.

Comments: to appear in npj 2D Materials and Applications

Density Functional Based Simulations of Proton Permeation of Graphene and Hexagonal Boron Nitride

Authors: J. M. H. Kroes, A. Fasolino, M. I. Katsnelson

Abstract: Using density functional theory, we study proton permeation through graphene and hexagonal boron nitride. We consider several factors influencing the barriers for permeation, including structural optimization, the role of the solvent, surface curvature and proton transport through hydrogenated samples. Furthermore, we discuss the ground state charge transfer from the membrane to the proton and the… ▽ More Using density functional theory, we study proton permeation through graphene and hexagonal boron nitride. We consider several factors influencing the barriers for permeation, including structural optimization, the role of the solvent, surface curvature and proton transport through hydrogenated samples. Furthermore, we discuss the ground state charge transfer from the membrane to the proton and the strong tendency for bond formation. If the process is assumed to be slow we find that none of these effects lead to a satisfactory answer to the observed discrepancies between theory and experiment. △ Less

Submitted 13 February, 2017; originally announced February 2017.

Comments: Physical Chemistry Chemical Physics, 2017

Emergent friction in two-dimensional Frenkel-Kontorova models

Authors: Jesper Norell, Annalisa Fasolino, Astrid S. de Wijn

Abstract: Simple models for friction are typically one-dimensional, but real interfaces are two-dimensional. We investigate the effects of the second dimension on static and dynamic friction by using the Frenkel-Kontorova (FK) model. We study the two most straightforward extensions of the FK model to two dimensions and simulate both the static and dynamic properties. We show that the behavior of the static… ▽ More Simple models for friction are typically one-dimensional, but real interfaces are two-dimensional. We investigate the effects of the second dimension on static and dynamic friction by using the Frenkel-Kontorova (FK) model. We study the two most straightforward extensions of the FK model to two dimensions and simulate both the static and dynamic properties. We show that the behavior of the static friction is robust and remains similar in two dimensions for physically reasonable parameter values. The dynamic friction, however, is strongly influenced by the second dimension and the accompanying additional dynamics and parameters introduced into the models. We discuss our results in terms of the thermal equilibration and phonon dispersion relations of the lattices, establishing a physically realistic and suitable two-dimensional extension of the FK model. We find that the presence of additional dissipation channels can increase the friction and produces significantly different temperature-dependence when compared to the one-dimensional case. We also briefly study the anisotropy of the dynamic friction and show highly nontrivial effects, including that the friction anisotropy can lead to motion in different directions depending on the value of the initial velocity. △ Less

Submitted 29 September, 2016; v1 submitted 6 April, 2016; originally announced April 2016.

Comments: 14 pages, 13 figures

Journal ref: Phys. Rev. E 94, 023001 (2016)

Macroscopic self-reorientation of interacting two-dimensional crystals

Authors: C. R. Woods, F. Withers, M. J. Zhu, Y. Cao, G. Yu, A. Kozikov, M. Ben Shalom, S. V. Morozov, M. M. van Wijk, A. Fasolino, M. I. Katsnelson, K. Watanabe, T. Taniguchi, A. K. Geim, A. Mishchenko, K. S. Novoselov

Abstract: Microelectromechanical systems, which can be moved or rotated with nanometre precision, already find applications in such fields like radio-frequency electronics, micro-attenuators, sensors and many others. Especially interesting are those which allow fine control over the motion on atomic scale due to self-alignment mechanisms and forces acting on the atomic level. Such machines can produce well-… ▽ More Microelectromechanical systems, which can be moved or rotated with nanometre precision, already find applications in such fields like radio-frequency electronics, micro-attenuators, sensors and many others. Especially interesting are those which allow fine control over the motion on atomic scale due to self-alignment mechanisms and forces acting on the atomic level. Such machines can produce well-controlled movements as a reaction to small changes of the external parameters. Here we demonstrate that, for the system of graphene on hexagonal boron nitride, the interplay between the van der Waals and elastic energies results in graphene mechanically self-rotating towards the hexagonal boron nitride crystallographic directions. Such rotation is macroscopic (for graphene flakes of tens of micrometres the tangential movement can be on hundreds of nanometres) and can be used for reproducible manufacturing of aligned van der Waals heterostructures. △ Less

Submitted 15 March, 2016; originally announced March 2016.

Comments: 32 pages, 17 figures, includes supplementary information, published in Nature Communications (2016)

Journal ref: Nature Communications 7, Article number:10800, 2016

Two-dimensional dispersion of magnetostatic volume spin waves

Authors: F. J. Buijnsters, L. J. A. van Tilburg, A. Fasolino, M. I. Katsnelson

Abstract: The dipolar (magnetostatic) interaction dominates the behavior of spin waves in magnetic films in the long-wavelength regime. In an in-plane magnetized film, volume modes exist with a negative group velocity (backward volume magnetostatic spin waves), in addition to the forward surface-localized mode (Damon-Eshbach). Inside the film of finite thickness $L$, the volume modes have a nontrivial spati… ▽ More The dipolar (magnetostatic) interaction dominates the behavior of spin waves in magnetic films in the long-wavelength regime. In an in-plane magnetized film, volume modes exist with a negative group velocity (backward volume magnetostatic spin waves), in addition to the forward surface-localized mode (Damon-Eshbach). Inside the film of finite thickness $L$, the volume modes have a nontrivial spatial dependence, and their two-dimensional dispersion relations $ω(\mathbf{k})$ can be calculated only numerically. We present explicit perturbative expressions for the profiles and frequencies of the volume modes, taking into account an in-plane applied field and uniaxial anisotropy, for the regimes $\lVert \mathbf{k}L \rVert \gg 1$ and $\lVert \mathbf{k}L \rVert \ll 1$, which together provide a good indication of the behavior of the modes for arbitrary wavevector $\mathbf{k}$. Moreover, we derive a very accurate semianalytical expression for the dispersion relation $ω(\mathbf{k})$ of the lowest-frequency mode that is straightforward to evaluate using standard numerical routines. Our results are useful to quickly interpret and control the excitation and propagation of spin waves in (opto-)magnetic experiments. △ Less

Submitted 3 February, 2016; originally announced February 2016.

Comments: 16 pages, 8 figures

Nonlinear effects in the propagation of optically generated magnetostatic volume mode spin waves

Authors: L. J. A. van Tilburg, F. J. Buijnsters, A. Fasolino, T. Rasing, M. I. Katsnelson

Abstract: Recent experimental work has demonstrated optical control of spin wave emission by tuning the shape of the optical pulse (Satoh et al.\ Nature Photonics, 6, 662 (2012)). We reproduce these results and extend the scope of the control by investigating nonlinear effects for large amplitude excitations. We observe an accumulation of spin wave power at the center of the initial excitation combined with… ▽ More Recent experimental work has demonstrated optical control of spin wave emission by tuning the shape of the optical pulse (Satoh et al.\ Nature Photonics, 6, 662 (2012)). We reproduce these results and extend the scope of the control by investigating nonlinear effects for large amplitude excitations. We observe an accumulation of spin wave power at the center of the initial excitation combined with short-wavelength spin waves. These kind of nonlinear effects have not been observed in earlier work on nonlinearities of spin waves. Our observations pave the way for the manipulation of magnetic structures at a smaller scale than the beam focus, for instance in devices with all-optical control of magnetism. △ Less

Submitted 31 August, 2017; v1 submitted 11 January, 2016; originally announced January 2016.

Comments: Added new figures to further illustrate the nonlinear effects to show time evolution and spectral flow. Added references. Changed perspective on nonlinear effects w.r.t. applicability of NSE. Added acknowledgement

Journal ref: Phys. Rev. B 96, 054437 (2017)

Scaling behavior and strain dependence of in-plane elastic properties of graphene

Authors: J. H. Los, A. Fasolino, M. I. Katsnelson

Abstract: We show by atomistic simulations that, in the thermodynamic limit, the in-plane elastic moduli of graphene at finite temperature vanish with system size $ L $ as a power law $ ~ L^{-η_u} $ with $ η_u \simeq 0.325 $, in agreement with the membrane theory. Our simulations clearly reveal the size and strain dependence of graphene's elastic moduli, allowing comparison to experimental data. Although th… ▽ More We show by atomistic simulations that, in the thermodynamic limit, the in-plane elastic moduli of graphene at finite temperature vanish with system size $ L $ as a power law $ ~ L^{-η_u} $ with $ η_u \simeq 0.325 $, in agreement with the membrane theory. Our simulations clearly reveal the size and strain dependence of graphene's elastic moduli, allowing comparison to experimental data. Although the recently measured difference of a factor 2 between the asymptotic value of the Young modulus for tensilely strained systems and the value from {\it ab initio} calculations remains unsolved, our results do explain the experimentally observed increase of more than a factor 2 for a tensile strain of only a few permille. We also discuss the scaling of the Poisson ratio, for which our simulations disagree with the predictions of the self-consistent screening approximation. △ Less

Submitted 21 August, 2015; originally announced August 2015.

Comments: 5 figures

Journal ref: Phys. Rev. Lett. 116, 015901 (2016)

Effect of structural relaxation on the electronic structure of graphene on hexagonal boron nitride

Authors: G. J. Slotman, M. M. van Wijk, Pei-Liang Zhao, A. Fasolino, M. I. Katsnelson, Shengjun Yuan

Abstract: We performed calculations of electronic, optical and transport properties of graphene on hBN with realistic moiré patterns. The latter are produced by structural relaxation using a fully atomistic model. This relaxation turns out to be crucially important for electronic properties. We describe experimentally observed features such as additional Dirac points and the "Hofstadter butterfly" structure… ▽ More We performed calculations of electronic, optical and transport properties of graphene on hBN with realistic moiré patterns. The latter are produced by structural relaxation using a fully atomistic model. This relaxation turns out to be crucially important for electronic properties. We describe experimentally observed features such as additional Dirac points and the "Hofstadter butterfly" structure of energy levels in a magnetic field. We find that the electronic structure is sensitive to many-body renormalization of the local energy gap. △ Less

Submitted 11 July, 2015; originally announced July 2015.

Comments: 5 pages, 6 figures. Supplementary material is available at http://www.theorphys.science.ru.nl/people/yuan/attachments/sm_hbn.pdf

Journal ref: Phys. Rev. Lett. 115, 186801 (2015)

Chirality-dependent transmission of spin waves through domain walls

Authors: F. J. Buijnsters, Y. Ferreiros, A. Fasolino, M. I. Katsnelson

Abstract: Spin-wave technology (magnonics) has the potential to further reduce the size and energy consumption of information processing devices. In the submicrometer regime (exchange spin waves), topological defects such as domain walls may constitute active elements to manipulate spin waves and perform logic operations. We predict that spin waves that pass through a domain wall in an ultrathin perpendicul… ▽ More Spin-wave technology (magnonics) has the potential to further reduce the size and energy consumption of information processing devices. In the submicrometer regime (exchange spin waves), topological defects such as domain walls may constitute active elements to manipulate spin waves and perform logic operations. We predict that spin waves that pass through a domain wall in an ultrathin perpendicular-anisotropy film experience a phase shift that depends on the orientation of the domain wall (chirality). The effect, which is absent in bulk materials, originates from the interfacial Dzyaloshinskii-Moriya interaction and can be interpreted as a geometric phase. We demonstrate analytically and by means of micromagnetic simulations that the phase shift is strong enough to switch between constructive and destructive interference. The two chirality states of the domain wall may serve as a memory bit or spin-wave switch in magnonic devices. △ Less

Submitted 7 March, 2016; v1 submitted 23 June, 2015; originally announced June 2015.

Comments: 11 pages, 10 figures (incl. supp. mat.); Phys. Rev. Lett. (accepted)

Journal ref: Phys. Rev. Lett. 116, 147204 (2016)

Superlubric to stick-slip sliding of incommensurate graphene flakes on graphite

Authors: M. M. van Wijk, M. Dienwiebel, J. W. M. Frenken, A. Fasolino

Abstract: We calculate the friction of fully mobile graphene flakes sliding on graphite. For incommensurately stacked flakes, we find a sudden and reversible increase in friction with load, in agreement with experimental observations. The transition from smooth sliding to stick-slip and the corresponding increase in friction is neither due to rotations to commensurate contact nor to dislocations but to a pi… ▽ More We calculate the friction of fully mobile graphene flakes sliding on graphite. For incommensurately stacked flakes, we find a sudden and reversible increase in friction with load, in agreement with experimental observations. The transition from smooth sliding to stick-slip and the corresponding increase in friction is neither due to rotations to commensurate contact nor to dislocations but to a pinning caused by vertical distortions of edge atoms also when they are saturated by Hydrogen. This behavior should apply to all layered materials with strong in-plane bonding. △ Less

Submitted 13 March, 2015; originally announced March 2015.

Journal ref: Phys. Rev. B 88, 235423 (2013)

Melting temperature of graphene

Authors: J. H. Los, K. V. Zakharchenko, M. I. Katsnelson, Annalisa Fasolino

Abstract: We present an approach to the melting of graphene based on nucleation theory for a first order phase transition from the 2D solid to the 3D liquid via an intermediate quasi-2D liquid. The applicability of nucleation theory, supported by the results of systematic atomistic Monte Carlo simulations, provides an intrinsic definition of the melting temperature of graphene, $ T_m $, and allows us to d… ▽ More We present an approach to the melting of graphene based on nucleation theory for a first order phase transition from the 2D solid to the 3D liquid via an intermediate quasi-2D liquid. The applicability of nucleation theory, supported by the results of systematic atomistic Monte Carlo simulations, provides an intrinsic definition of the melting temperature of graphene, $ T_m $, and allows us to determine it. We find $T_m \simeq 4510$ K, about 250 K higher than that of graphite using the same interatomic interaction model. The found melting temperature is shown to be in good agreement with the asymptotic results of melting simulations for finite disks and ribbons of graphene. Our results strongly suggest that graphene is the most refractory of all known materials. △ Less

Submitted 12 March, 2015; originally announced March 2015.

Journal ref: Phys. Rev. B 91, 045415 (2015)

Toward Reverse Engineering of VBA Based Excel Spreadsheet Applications

Authors: Domenico Amalfitano, Nicola Amatucci, Vincenzo De Simone, Anna Rita Fasolino, Porfirio Tramontana

Abstract: Modern spreadsheet systems can be used to implement complex spreadsheet applications including data sheets, customized user forms and executable procedures written in a scripting language. These applications are often developed by practitioners that do not follow any software engineering practice and do not produce any design documentation. Thus, spreadsheet applications may be very difficult to b… ▽ More Modern spreadsheet systems can be used to implement complex spreadsheet applications including data sheets, customized user forms and executable procedures written in a scripting language. These applications are often developed by practitioners that do not follow any software engineering practice and do not produce any design documentation. Thus, spreadsheet applications may be very difficult to be maintained or restructured. In this position paper we present in a nutshell two reverse engineering techniques and a tool that we are currently realizing for the abstraction of conceptual data models and business logic models. △ Less

Submitted 11 March, 2015; originally announced March 2015.

Comments: In Proceedings of the 2nd Workshop on Software Engineering Methods in Spreadsheets (http://spreadsheetlab.org/sems15/)

Relaxation of moiré patterns for slightly misaligned identical lattices: graphene on graphite

Authors: M. M. van Wijk, A. Schuring, M. I. Katsnelson, A. Fasolino

Abstract: We study the effect of atomic relaxation on the structure of moiré patterns in twisted graphene on graphite and double layer graphene by large scale atomistic simulations. The reconstructed structure can be described as a superlattice of `hot spots' with vortex-like behaviour of in-plane atomic displacements and increasing out-of-plane displacements with decreasing angle. These lattice distortions… ▽ More We study the effect of atomic relaxation on the structure of moiré patterns in twisted graphene on graphite and double layer graphene by large scale atomistic simulations. The reconstructed structure can be described as a superlattice of `hot spots' with vortex-like behaviour of in-plane atomic displacements and increasing out-of-plane displacements with decreasing angle. These lattice distortions affect both scalar and vector potential and the resulting electronic properties. At low misorientation angles (<$\sim$1$^\circ$) the optimized structures deviate drastically from the sinusoidal modulation which is often assumed in calculations of the electronic properties. The proposed structure might be verified by scanning probe microscopy measurements. △ Less

Submitted 9 March, 2015; originally announced March 2015.

Journal ref: 2D materials 2, 034010 (2015)

Reply to 'Comment on "Thermodynamics of quantum crystalline membranes"'

Authors: B. Amorim, R. Roldán, E. Cappelluti, F. Guinea, A. Fasolino, M. I. Katsnelson

Abstract: In this note, we reply to the comment made by E.I.Kats and V.V.Lebedev [arXiv:1407.4298] on our recent work "Thermodynamics of quantum crystalline membranes" [Phys. Rev. B 89, 224307 (2014)]. Kats and Lebedev question the validity of the calculation presented in our work, in particular on the use of a Debye momentum as a ultra-violet regulator for the theory. We address and counter argue the criti… ▽ More In this note, we reply to the comment made by E.I.Kats and V.V.Lebedev [arXiv:1407.4298] on our recent work "Thermodynamics of quantum crystalline membranes" [Phys. Rev. B 89, 224307 (2014)]. Kats and Lebedev question the validity of the calculation presented in our work, in particular on the use of a Debye momentum as a ultra-violet regulator for the theory. We address and counter argue the criticisms made by Kats and Lebedev to our work. △ Less

Submitted 9 October, 2014; originally announced October 2014.

Comments: 5 pages, 4 figures

Motion of domain walls and the dynamics of kinks in the magnetic Peierls potential

Authors: F. J. Buijnsters, A. Fasolino, M. I. Katsnelson

Abstract: We study the dynamics of magnetic domain walls in the Peierls potential due to the discreteness of the crystal lattice. The propagation of a narrow domain wall (comparable to the lattice parameter) under the effect of a magnetic field proceeds through the formation of kinks in its profile. We predict that, despite the discreteness of the system, such kinks can behave like sine-Gordon solitons in t… ▽ More We study the dynamics of magnetic domain walls in the Peierls potential due to the discreteness of the crystal lattice. The propagation of a narrow domain wall (comparable to the lattice parameter) under the effect of a magnetic field proceeds through the formation of kinks in its profile. We predict that, despite the discreteness of the system, such kinks can behave like sine-Gordon solitons in thin films of materials such as yttrium iron garnets, and we derive general conditions for other materials. In our simulations we also observe long-lived breathers. We provide analytical expressions for the effective mass and limiting velocity of the kink in excellent agreement with our numerical results. △ Less

Submitted 29 July, 2014; originally announced July 2014.

Comments: 12 pages, 9 figures (incl. supp. mat.)

Journal ref: Phys. Rev. Lett. 113, 217202 (2014)

Phonons and electron-phonon coupling in graphene-h-BN heterostructures

Authors: Guus J. Slotman, Gilles A. de Wijs, Annalisa Fasolino, Mikhail I. Katsnelson

Abstract: First principle calculations of the phonons of graphene-BN heterostructures are presented and compared to those of the constituents. We show that AA and AB' stacking are not only energetically less favoured than AB but also dynamically unstable. We have identified low energy flat phonon branches of BN character with out of plane displacement and evaluated their coupling to electrons in graphene. First principle calculations of the phonons of graphene-BN heterostructures are presented and compared to those of the constituents. We show that AA and AB' stacking are not only energetically less favoured than AB but also dynamically unstable. We have identified low energy flat phonon branches of BN character with out of plane displacement and evaluated their coupling to electrons in graphene. △ Less

Submitted 24 July, 2014; originally announced July 2014.

Comments: 5 figures

Moir{é} patterns as a probe of interplanar interactions: graphene on h-BN

Authors: M. M. van Wijk, A. Schuring, M. I. Katsnelson, A. Fasolino

Abstract: By atomistic modeling of moir{é} patterns of graphene on a substrate with a small lattice mismatch, we find qualitatively different strain distributions for small and large misorientation angles, corresponding to the commensurate-incommensurate transition recently observed in graphene on hexagonal BN. We find that the ratio of C-N and C-B interactions is the main parameter determining the differen… ▽ More By atomistic modeling of moir{é} patterns of graphene on a substrate with a small lattice mismatch, we find qualitatively different strain distributions for small and large misorientation angles, corresponding to the commensurate-incommensurate transition recently observed in graphene on hexagonal BN. We find that the ratio of C-N and C-B interactions is the main parameter determining the different bond lengths in the center and edges of the moir{é} pattern. Agreement with experimental data is obtained only by assuming that the C-B interactions are at least twice weaker than the C-N interactions. The correspondence between the strain distribution in the nanoscale moir{é} pattern and the potential energy surface at the atomic scale found in our calculations, makes the moir{é} pattern a tool to study details of dispersive forces in van der Waals heterostructures. △ Less

Submitted 13 March, 2015; v1 submitted 10 July, 2014; originally announced July 2014.

Comments: 5 pages, 3 figures

Journal ref: Phys. Rev. Lett. 113, 135504 (2014)

Minimal graphene thickness for wear protection of diamond

Authors: M. M. van Wijk, A. Fasolino

Abstract: We show by means of molecular dynamics simulations that graphene is an excellent coating for diamond. The transformation of diamond to amorphous carbon while sliding under pressure can be prevented by having at least two graphene layers between the diamond slabs, making this combination of materials suitable for new coatings and micro- and nanoelectromechanical devices. Grain boundaries, vacancies… ▽ More We show by means of molecular dynamics simulations that graphene is an excellent coating for diamond. The transformation of diamond to amorphous carbon while sliding under pressure can be prevented by having at least two graphene layers between the diamond slabs, making this combination of materials suitable for new coatings and micro- and nanoelectromechanical devices. Grain boundaries, vacancies and adatoms on the diamond surface do not change this picture whereas reactive adsorbates between the graphene layers may have detrimental effects. Our findings can be explained by the properties of layered materials where the weak interlayer bonding evolves to a strong interlayer repulsion under pressure. △ Less

Submitted 14 April, 2014; originally announced April 2014.

Journal ref: AIP Advances 5, 017117 (2015)

Thermodynamics of quantum crystalline membranes

Authors: B. Amorim, R. Roldán, E. Cappelluti, A. Fasolino, F. Guinea, M. I. Katsnelson

Abstract: We investigate the thermodynamic properties and the lattice stability of two-dimensional crystalline membranes, such as graphene and related compounds, in the low temperature quantum regime $T\rightarrow0$. A key role is played by the anharmonic coupling between in-plane and out-of plane lattice modes that, in the quantum limit, has very different consequences than in the classical regime. The rol… ▽ More We investigate the thermodynamic properties and the lattice stability of two-dimensional crystalline membranes, such as graphene and related compounds, in the low temperature quantum regime $T\rightarrow0$. A key role is played by the anharmonic coupling between in-plane and out-of plane lattice modes that, in the quantum limit, has very different consequences than in the classical regime. The role of retardation, namely of the frequency dependence, in the effective anharmonic interactions turns out to be crucial in the quantum regime. We identify a crossover temperature, $T^{*}$, between classical and quantum regimes, which is $\sim 70 - 90$ K for graphene. Below $T^{*}$, the heat capacity and thermal expansion coefficient decrease as power laws with decreasing temperature, tending to zero for $T\rightarrow0$ as required by the third law of thermodynamics. △ Less

Submitted 30 June, 2014; v1 submitted 11 March, 2014; originally announced March 2014.

Comments: 13 pages, 1 figure

Journal ref: Phys. Rev. B 89, 224307 (2014)

Zero modes in magnetic systems: general theory and an efficient computational scheme

Authors: F. J. Buijnsters, A. Fasolino, M. I. Katsnelson

Abstract: The presence of topological defects in magnetic media often leads to normal modes with zero frequency (zero modes). Such modes are crucial for long-time behavior, describing, for example, the motion of a domain wall as a whole. Conventional numerical methods to calculate the spin-wave spectrum in magnetic media are either inefficient or they fail for systems with zero modes. We present a new effic… ▽ More The presence of topological defects in magnetic media often leads to normal modes with zero frequency (zero modes). Such modes are crucial for long-time behavior, describing, for example, the motion of a domain wall as a whole. Conventional numerical methods to calculate the spin-wave spectrum in magnetic media are either inefficient or they fail for systems with zero modes. We present a new efficient computational scheme that reduces the magnetic normal-mode problem to a generalized Hermitian eigenvalue problem also in the presence of zero modes. We apply our scheme to several examples, including two-dimensional domain walls and Skyrmions, and show how the effective masses that determine the dynamics can be calculated directly. These systems highlight the fundamental distinction between the two types of zero modes that can occur in spin systems, which we call special and inertial zero modes. Our method is suitable for both conservative and dissipative systems. For the latter case, we present a perturbative scheme to take into account damping, which can also be used to calculate dynamical susceptibilities. △ Less

Submitted 5 December, 2013; originally announced December 2013.

Comments: 64 pages, 15 figures

Journal ref: Phys. Rev. B 89, 174433 (2014)

Nanoscopic friction under electrochemical control

Authors: A. S. de Wijn, A. Fasolino, A. Filippov, M. Urbakh

Abstract: We propose a theoretical model of friction under electrochemical conditions focusing on the interaction of a force microscope tip with adsorbed polar molecules of which the orientation depends on the applied electric field. We demonstrate that the dependence of friction force on the electric field is determined by the interplay of two channels of energy dissipation: (i) the rotation of dipoles and… ▽ More We propose a theoretical model of friction under electrochemical conditions focusing on the interaction of a force microscope tip with adsorbed polar molecules of which the orientation depends on the applied electric field. We demonstrate that the dependence of friction force on the electric field is determined by the interplay of two channels of energy dissipation: (i) the rotation of dipoles and (ii) slips of the tip over potential barriers. We suggest a promising strategy to achieve a strong dependence of nanoscopic friction on the external field based on the competition between long range electrostatic interactions and short range chemical interactions between tip and adsorbed polar molecules. △ Less

Submitted 25 September, 2013; originally announced September 2013.

Comments: 5 pages, 5 figures

Journal ref: Phys. Rev. Lett. 112, 055502 (2014)

Concavity Effects on the Optical Properties of Aromatic Hydrocarbons

Authors: Caterina Cocchi, Deborah Prezzi, Alice Ruini, Marilia J. Caldas, Annalisa Fasolino, Elisa Molinari

Abstract: We study the modifications on the ground and excited state properties of polycyclic aromatic hydrocarbons (PAHs), induced by the variation of concavity and $π$-connectivity. Inspired by experimentally feasible systems, we study three series of PAHs, from H-saturated graphene flakes to geodesic buckybowls, related to the formation of fullerene C60 and C50-carbon nanotube caps. Working within the fr… ▽ More We study the modifications on the ground and excited state properties of polycyclic aromatic hydrocarbons (PAHs), induced by the variation of concavity and $π$-connectivity. Inspired by experimentally feasible systems, we study three series of PAHs, from H-saturated graphene flakes to geodesic buckybowls, related to the formation of fullerene C60 and C50-carbon nanotube caps. Working within the framework of quantum chemistry semi-empirical methods AM1 and ZINDO/S, we find that the interplay between concavity and π-connectivity shifts the bright optical lines to higher energies, and introduces symmetry-forbidden dark excitations at low energy. A generally good agreement with the available experimental data supports our results, which can be viewed as the basis for designing optical properties of novel curved aromatic molecules. △ Less

Submitted 12 June, 2013; originally announced June 2013.

Journal ref: J. Phys. Chem. C, 117,12909-12915 (2013)

k.p subband structure of the LaAlO3/SrTiO3 interface

Authors: L. W. van Heeringen, G. A. de Wijs, A. McCollam, J. C. Maan, A. Fasolino

Abstract: 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 ed… ▽ More 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. △ Less

Submitted 8 January, 2014; v1 submitted 28 March, 2013; originally announced March 2013.

Comments: 8 pages, 7 figures

Journal ref: Phys. Rev. B 88,205140 (2013)

Structure and magnetism of disordered carbon

Authors: M. A. Akhukov, M. I. Katsnelson, A. Fasolino

Abstract: 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… ▽ More 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. △ Less

Submitted 1 March, 2013; originally announced March 2013.

Comments: 14 pages, 11 figures

Structure, stability and defects of single layer h-BN in comparison to graphene

Authors: GJ Slotman, A Fasolino

Abstract: 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… ▽ More 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. △ Less

Submitted 21 February, 2013; originally announced February 2013.

Comments: 13 pages, 7 figures

Journal ref: J Phys Condens Matter. 2013 Jan 30;25(4):045009. Epub 2012 Dec 18

Graphene as a Prototype Crystalline Membrane

Authors: Mikhail I. Katsnelson, Annalisa Fasolino

Abstract: 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 biologi… ▽ More 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. △ Less

Submitted 7 February, 2013; v1 submitted 6 February, 2013; originally announced February 2013.

Comments: 20 pages, 9 figures. Misprints are fixed and text style is changed

Journal ref: Accounts of Chemical Research 46, 97 (2013)

Electronic, Magnetic and Transport Properties of Graphene Ribbons Terminated by Nanotubes

Authors: M. A. Akhukov, Shengjun Yuan, A. Fasolino, M. I. Katsnelson

Abstract: 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 junc… ▽ More 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. △ Less

Submitted 12 October, 2012; originally announced October 2012.

Comments: 9 pages, 5 figures

Effect of temperature and velocity on superlubricity

Authors: Joost A. van den Ende, Astrid S. de Wijn, Annalisa Fasolino

Abstract: 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 we… ▽ More 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. △ Less

Submitted 3 October, 2012; originally announced October 2012.

Journal ref: Journal of Physics: Condensed Matter 24, 445009 (2012)

Dangling bonds and magnetism of grain boundaries in graphene

Authors: M. A. Akhukov, A. Fasolino, Y. N. Gornostyrev, M. I. Katsnelson

Abstract: 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… ▽ More 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. △ Less

Submitted 24 February, 2012; v1 submitted 1 December, 2011; originally announced December 2011.

Comments: 10 pages, 11 figures, 4 tables, the final version accepted in prb

Journal ref: Phys. Rev. B 85, 115407 (2012)

Low friction and rotational dynamics of crystalline flakes in solid lubrication

Authors: A. S. de Wijn, A. Fasolino, A. E. Filippov, M. Urbakh

Abstract: 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 simulatio… ▽ More 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. △ Less

Submitted 11 August, 2011; originally announced August 2011.

Comments: 5 pages, 5 figures

Journal ref: EPL 95 66002 (2011)

Relating chaos to deterministic diffusion of a molecule adsorbed on a surface

Authors: Astrid S. de Wijn, Annalisa Fasolino

Abstract: 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 wi… ▽ More 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. △ Less

Submitted 13 July, 2011; originally announced July 2011.

Journal ref: J. Phys.: Condens. Matter 21 264002 (2009)

Stability of superlubric sliding on graphite

Authors: Astrid S. de Wijn, Claudio Fusco, Annalisa Fasolino

Abstract: 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: particula… ▽ More 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. △ Less

Submitted 7 July, 2011; originally announced July 2011.

Comments: 11 pages, 15 figures

Journal ref: Phys. Rev. E 81, 046105 (2010)

Melting of graphene: from two to one dimension

Authors: K. V. Zakharchenko, Annalisa Fasolino, J. H. Los, M. I. Katsnelson

Abstract: 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 liqu… ▽ More 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. △ Less

Submitted 6 April, 2011; originally announced April 2011.

Comments: 12 pages, 6 figures

Suppression of anharmonicities in crystalline membranes by external strain

Authors: Rafael Roldán, Annalisa Fasolino, Kostyantyn V. Zakharchenko, Mikhail I. Katsnelson

Abstract: 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 b… ▽ More 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. △ Less

Submitted 11 May, 2011; v1 submitted 31 January, 2011; originally announced January 2011.

Comments: 7 pages, 6 figures

Journal ref: Phys. Rev. B 83, 174104 (2011)

Phonons of graphene and graphitic materials derived from the empirical potential LCBOPII

Authors: L. J. Karssemeijer, Annalisa Fasolino

Abstract: 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 constan… ▽ More 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. △ Less

Submitted 27 October, 2010; originally announced October 2010.

Journal ref: Surf. Sci. 605, 1611 (2011)

Bending modes, elastic constants and mechanical stability of graphitic systems

Authors: G. Savini, Y. J. Dappe, S. Öberg, J. -C. Charlier, M. I. Katsnelson, A. Fasolino

Abstract: 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 exfolia… ▽ More 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. △ Less

Submitted 21 August, 2010; originally announced August 2010.

Comments: 7 pages, 6 figures

Self-Consistent Screening Approximation for Flexible Membranes: Application to Graphene

Authors: K. V. Zakharchenko, R. Roldan, A. Fasolino, M. I. Katsnelson

Abstract: 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 renorma… ▽ More 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}$ Å$^{-1}$ till 10 Å$^{-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}$ Å$^{-1}$ till 1 Å$^{-1}$. In the limit $q\rightarrow 0$ our data for the exponent $η$ of the renormalized bending rigidity $κ_R(q)\propto q^{-η}$ is compatible with the previously known analytical results for the SCSA $η\simeq 0.82$. However, this limit appears to be reached only for $q<10^{-5}$ Å$^{-1}$ whereas at intermediate $q$ the behavior of $G(q)$ cannot be described by a single exponent. △ Less

Submitted 20 September, 2010; v1 submitted 8 June, 2010; originally announced June 2010.

Comments: 5 pages, 4 figures

Journal ref: Phys. Rev. B 82, 125435 (2010)

Non-spherical shapes of capsules within a fourth-order curvature model

Authors: O. V. Manyuhina, J. J. Hetzel, M. I. Katsnelson, A. Fasolino

Abstract: 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-crys… ▽ More 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. △ Less

Submitted 7 June, 2010; originally announced June 2010.

Comments: 6 pages, 4 figures (LaTeX macros EPJ), accepted for publication in EPJ E

Journal ref: Eur. Phys. J. E032: 223-228,2010

Atomistic simulations of structural and thermodynamic properties of bilayer graphene

Authors: K. V. Zakharchenko, J. H. Los, M. I. Katsnelson, A. Fasolino

Abstract: 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 expan… ▽ More 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}$ △ Less

Submitted 1 April, 2010; originally announced April 2010.

Comments: 6 pages, 7 figures.

Stability of low-friction surface sliding of nanocrystals with rectangular symmetry and application to W on NaF(001)

Authors: Astrid S. de Wijn, Annalisa Fasolino

Abstract: 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 incomme… ▽ More 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. △ Less

Submitted 31 March, 2010; originally announced April 2010.

Journal ref: Tribology Letters 39, 91-99 (2010)

Stable and fast semi-implicit integration of the stochastic Landau-Lifshitz equation

Authors: J. H. Mentink, M. V. Tretyakov, A. Fasolino, M. I. Katsnelson, Th. Rasing

Abstract: 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… ▽ More 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. △ Less

Submitted 9 February, 2010; originally announced February 2010.

Comments: 24 pages, 5 figures

Journal ref: J. Phys.: Condens. Matter 22 176001 (2010)

Electron density distribution and screening in rippled graphene sheets

Authors: Marco Gibertini, Andrea Tomadin, Marco Polini, A. Fasolino, M. I. Katsnelson

Abstract: 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 ext… ▽ More 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. △ Less

Submitted 5 February, 2010; originally announced February 2010.

Comments: 13 pages, 13 figures, submitted. High-quality figures can be requested to the authors

Journal ref: Phys. Rev. B 81, 125437 (2010)

Topological defects and shape of aromatic self-assembled vesicles

Authors: O. V. Manyuhina, A. Fasolino, M. I. Katsnelson

Abstract: 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. 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. △ Less

Submitted 10 July, 2009; originally announced July 2009.

Comments: accepted as a Letter in the J. Phys. Chem. B

Journal ref: J. Phys. Chem. B 113 (31), 10549 (2009)

Scaling Properties of Flexible Membranes from Atomistic Simulations: Application to Graphene

Authors: J. H. Los, M. I. Katsnelson, O. V. Yazyev, K. V. Zakharchenko, A. Fasolino

Abstract: 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 Mo… ▽ More 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+η}$ with the same exponent $η\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. △ Less

Submitted 23 March, 2009; originally announced March 2009.

Finite temperature lattice properties of graphene beyond the quasiharmonic approximation

Authors: K. V. Zakharchenko, M. I. Katsnelson, A. Fasolino

Abstract: The thermal and mechanical stability of graphene is important for many potential applications in nanotechnology. We calculate the temperature dependence of lattice parameter, elastic properties and heat capacity by means of atomistic Monte Carlo simulations that allow to go beyond the quasiharmonic approximation. We predict an unusual, non-monotonic, behavior of the lattice parameter with minimu… ▽ More The thermal and mechanical stability of graphene is important for many potential applications in nanotechnology. We calculate the temperature dependence of lattice parameter, elastic properties and heat capacity by means of atomistic Monte Carlo simulations that allow to go beyond the quasiharmonic approximation. We predict an unusual, non-monotonic, behavior of the lattice parameter with minimum at temperature about 900 K and of the shear modulus with maximum at the same temperature. The Poisson ratio in graphene is found to be small ~0.1 in a broad temperature interval. △ Less

Submitted 22 December, 2008; originally announced December 2008.

Comments: 4 pages, 5 figures

State-of-the-art models for the phase diagram of carbon and diamond nucleation

Authors: Luca M. Ghiringhelli, C. Valeriani, J. H. Los, E. J. Meijer, A. Fasolino, D. Frenkel

Abstract: We review recent developments in the modelling of the phase diagram and the kinetics of crystallization of carbon. In particular, we show that a particular class of bond-order potentials (the so-called LCBOP models) account well for many of the known structural and thermodynamic properties of carbon at high pressures and temperatures. We discuss the LCBOP models in some detail. In addition, we b… ▽ More We review recent developments in the modelling of the phase diagram and the kinetics of crystallization of carbon. In particular, we show that a particular class of bond-order potentials (the so-called LCBOP models) account well for many of the known structural and thermodynamic properties of carbon at high pressures and temperatures. We discuss the LCBOP models in some detail. In addition, we briefly review the ``history'' of experimental and theoretical studies of the phase behaviour of carbon. Using a well-tested version of the LCBOP model (viz. LCBOPI+) we address some of the more controversial hypotheses concerning the phase behaviour of carbon, in particular: the suggestion that liquid carbon can exist in two phases separated by a first-order phase transition and the conjecture that diamonds could have formed by homogeneous nucleation in Uranus and Neptune. △ Less

Submitted 15 July, 2009; v1 submitted 10 April, 2008; originally announced April 2008.

Comments: 28 pages, 15 figures. To appear in Mol. Phys. in 2008

Topological instabilities of spherical vesicles

Authors: O. V. Manyuhina, A. Fasolino, P. C. M. Christianen, M. I. Katsnelson

Abstract: Within the framework of the Helfrich elastic theory of membranes and of differential geometry we study the possible instabilities of spherical vesicles towards double bubbles. We find that not only temperature, but also magnetic fields can induce topological transformations between spherical vesicles and double bubbles and provide a phase diagram for the equilibrium shapes. Within the framework of the Helfrich elastic theory of membranes and of differential geometry we study the possible instabilities of spherical vesicles towards double bubbles. We find that not only temperature, but also magnetic fields can induce topological transformations between spherical vesicles and double bubbles and provide a phase diagram for the equilibrium shapes. △ Less

Submitted 10 July, 2009; v1 submitted 20 February, 2008; originally announced February 2008.

Comments: 5 pages, 4 figures

Journal ref: Phys. Rev. E 80, 010403(R) (2009)