Source author record

Irina V. Lebedeva

Irina V. Lebedeva appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

21works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

21 published item(s)

preprint2021arXiv

Transformation of a graphene nanoribbon into a hybrid 1D nanoobject with alternating double chains and polycyclic regions

Molecular dynamics simulations show that a graphene nanoribbon with alternating regions which are one and three hexagons wide can transform into a hybrid 1D nanoobject with alternating double chains and polycyclic regions under electron irradiation in HRTEM. A scheme of synthesis of such a nanoribbon using Ullmann coupling and dehydrogenation reactions is proposed. The reactive REBO-1990EVC potential is adapted for simulations of carbon-hydrogen systems and is used in combination with the CompuTEM algorithm for modeling of electron irradiation effects. The atomistic mechanism of formation of the new hybrid 1D nanoobject is found to be the following. Firstly hydrogen is removed by electron impacts. Then spontaneous breaking of bonds between carbon atoms leads to decomposition of narrow regions of the graphene nanoribbon into double chains. Simultaneously thermally activated growth of polycyclic regions occurs. Density functional theory calculations give barriers along the growth path of polycyclic regions consistent with this mechanism. Electronic properties of the new 1D nanoobject are shown to be strongly affected by the edge magnetism and make this nanostructure promising nanoelectronic and spintronic applications. The way of synthesis the 1D nanoobject proposed here can be considered as an example of the general three-stage strategy of production of nanoobjects and macromolecules. 1) precursors are synthesized by traditional chemical method, 2) precursors are placed in HRTEM with the electron energy that is sufficient only to remove hydrogen atoms, 3) as a result of hydrogen removal, the precursors become unstable or metastable and transform into new nanoobjects or macromolecules.

preprint2020arXiv

Energetics and Structure of Domain Wall Networks in Minimally Twisted Bilayer Graphene under Strain

The parameters of the triangular domain wall network in bilayer graphene with a simultaneously twisted and biaxially stretched bottom layer are studied using the two-chain Frenkel-Kontorova model. It is demonstrated that if the graphene layers are free to rotate, they prefer to stay co-aligned upon stretching the bottom layer and the regular triangular network of tensile domain walls is formed upon the commensurate-incommensurate phase transition. If the angle between the layers is fixed, the regular triangular network of shear domain walls is observed at zero elongation of the bottom layer. Upon stretching the bottom layer, however, the domain walls transform into the tensile ones and the size of the commensurate domains decreases. We also show that the parameters of the isosceles triangular domain wall network in twisted bilayer graphene under shear strain can be determined through purely geometrical considerations. Experimental analysis of the orientation of domain walls and period of the triangular network would, on the one hand, contribute to understanding of the interlayer interaction of graphene layers, and, on the other hand, serve for detection of relative strains and rotation between the layers. Vice versa external strains can be used to control the parameters of the triangular domain wall network and, therefore, electronic properties of twisted bilayer graphene.

preprint2020arXiv

The CECAM Electronic Structure Library and the modular software development paradigm

First-principles electronic structure calculations are very widely used thanks to the many successful software packages available. Their traditional coding paradigm is monolithic, i.e., regardless of how modular its internal structure may be, the code is built independently from others, from the compiler up, with the exception of linear-algebra and message-passing libraries. This model has been quite successful for decades. The rapid progress in methodology, however, has resulted in an ever increasing complexity of those programs, which implies a growing amount of replication in coding and in the recurrent re-engineering needed to adapt to evolving hardware architecture. The Electronic Structure Library (\esl) was initiated by CECAM (European Centre for Atomic and Molecular Calculations) to catalyze a paradigm shift away from the monolithic model and promote modularization, with the ambition to extract common tasks from electronic structure programs and redesign them as free, open-source libraries. They include "heavy-duty" ones with a high degree of parallelisation, and potential for adaptation to novel hardware within them, thereby separating the sophisticated computer science aspects of performance optimization and re-engineering from the computational science done by scientists when implementing new ideas. It is a community effort, undertaken by developers of various successful codes, now facing the challenges arising in the new model. This modular paradigm will improve overall coding efficiency and enable specialists (computer scientists or computational scientists) to use their skills more effectively. It will lead to a more sustainable and dynamic evolution of software as well as lower barriers to entry for new developers.

preprint2020arXiv

Two phases with different domain wall networks and a reentrant phase transition in bilayer graphene under strain

The analytical two-chain Frenkel-Kontorova model is used to describe domain wall networks in bilayer graphene upon biaxial stretching of one of the layers. We show that the commensurate-incommensurate phase transition leading to formation of a regular triangular domain wall network at the relative biaxial elongation of $3.0\cdot 10^{-3}$ is followed by the transition to another incommensurate phase with a striped network at the elongation of $3.7\cdot 10^{-3}$. The reentrant transition to the phase with a triangular domain wall network is predicted for the elongation $\sim10^{-2}$.

preprint2020arXiv

Universal description of potential energy surface of interlayer interaction in two-dimensional materials by first spatial Fourier harmonics

We propose a hypothesis that the potential energy surface (PES) of interlayer interaction in diverse 2D materials can be universally described by the first spatial Fourier harmonics. This statement (checked previously for the interactions between graphene and hexagonal boron nitride layers in different combinations) is verified in the present paper for the case of hydrofluorinated graphene (HFG) bilayer with hydrogen bonding between fluorine and hydrogen at the interlayer interface. The PES for HFG bilayer is obtained through density functional theory calculations with van der Waals corrections. An analytical expression based on the first Fourier harmonics describing the PES which corresponds to the symmetry of HFG layers is derived. It is found that the calculated PES can be described by the first Fourier harmonics with the accuracy of 3\% relative to the PES corrugation. The shear mode frequency, shear modulus and barrier for relative rotation of the layers to incommensurate states of HFG bilayer are estimated. Additionally it is shown that HFG bilayer is stable relative to the formation of HF molecules as a result of chemical reactions between the layers.

preprint2019arXiv

Octopus, a computational framework for exploring light-driven phenomena and quantum dynamics in extended and finite systems

Over the last years extraordinary advances in experimental and theoretical tools have allowed us to monitor and control matter at short time and atomic scales with a high-degree of precision. An appealing and challenging route towards engineering materials with tailored properties is to find ways to design or selectively manipulate materials, especially at the quantum level. To this end, having a state-of-the-art ab initio computer simulation tool that enables a reliable and accurate simulation of light-induced changes in the physical and chemical properties of complex systems is of utmost importance. The first principles real-space-based Octopus project was born with that idea in mind, providing an unique framework allowing to describe non-equilibrium phenomena in molecular complexes, low dimensional materials, and extended systems by accounting for electronic, ionic, and photon quantum mechanical effects within a generalized time-dependent density functional theory framework. The present article aims to present the new features that have been implemented over the last few years, including technical developments related to performance and massive parallelism. We also describe the major theoretical developments to address ultrafast light-driven processes, like the new theoretical framework of quantum electrodynamics density-functional formalism (QEDFT) for the description of novel light-matter hybrid states. Those advances, and other being released soon as part of the Octopus package, will enable the scientific community to simulate and characterize spatial and time-resolved spectroscopies, ultrafast phenomena in molecules and materials, and new emergent states of matter (QED-materials).

preprint2016arXiv

Interlayer interaction and related properties of bilayer hexagonal boron nitride: ab initio study

Principal characteristics of interlayer interaction and relative motion of hexagonal boron nitride (h-BN) layers are investigated by the first-principles method taking into account van der Waals interactions. Dependences of the interlayer interaction energy on relative translational displacement of h-BN layers (potential energy surfaces) are calculated for two relative orientations of the layers, namely, for the layers aligned in the same direction and in the opposite directions upon relative rotation of the layers by 180 degrees. It is shown that the potential energy surfaces of bilayer h-BN can be approximated by the first Fourier components determined by symmetry. As a result, a wide set of physical qualintities describing relative motion of h-BN layers aligned in the same direction including barriers to their relative sliding and rotation, shear mode frequency and shear modulus are determined by a single parameter corresponding to the roughness of the potential energy surface, similar to bilayer graphene. The properties of h-BN layers aligned in the opposite directions are described by two such parameters. The possibility of partial and full dislocations in stacking of the layers is predicted for h-BN layers aligned in the same and opposite directions, respectively. The extended two-chain Frenkel-Kontorova model is used to estimate the width and formation energy of these dislocations on the basis of the calculated potential energy surfaces.

preprint2014arXiv

Atomistic mechanism of carbon nanotube cutting catalyzed by nickel under the electron beam

The cutting of single-walled carbon nanotubes by an 80 keV electron beam catalyzed by nickel clusters is imaged in situ using aberration-corrected high-resolution transmission electron microscopy. Extensive molecular dynamics simulations within the CompuTEM approach provide insight into the mechanism of this process and demonstrate that the combination of irradiation and nickel catalyst is crucial for the cutting process to take place. The atomistic mechanism of cutting is revealed by detailed analysis of irradiation-induced reactions of bonds reorganization and atom ejection in the vicinity of the nickel cluster, showing a highly complex interplay of different chemical transformations catalysed by the metal cluster. One of the most prevalent pathways includes three consecutive stages: formation of polyyne carbon chains from carbon nanotube, dissociation of the carbon chains into single and pairs of adatoms adsorbed on the nickel cluster, and ejection of these adatoms leading to the cutting of nanotube. Significant variations in the atom ejection rate are discovered depending on the process stage and nanotube diameter. The revealed mechanism and kinetic characteristics of cutting process provide fundamental knowledge for the development of new methodologies for control and manipulation of carbon structures at the nanoscale.

preprint2014arXiv

Identification of structural motifs as tunneling two-level systems in amorphous alumina at low temperatures

One of the most accepted models that describe the anomalous thermal behavior of amorphous materials at temperatures below 1 K relies on the quantum mechanical tunneling of atoms between two nearly equivalent potential energy wells forming a two-level system (TLS). Indirect evidence for TLSs is widely available. However, the atomistic structure of these TLSs remains an unsolved topic in the physics of amorphous materials. Here, using classical molecular dynamics, we found several hitherto unknown bistable structural motifs that may be key to understanding the anomalous thermal properties of amorphous alumina at low temperatures. We show through free energy profiles that the complex potential energy surface can be reduced to canonical TLSs. The predicted tunnel splittings from instanton theory, the number density, dipole moment, and coupling to external strain of the discovered motifs are consistent with experiments.

preprint2013arXiv

AA stacking, tribological and electronic properties of double-layer graphene with krypton spacer

Structural, energetic and tribological characteristics of double-layer graphene with commensurate and incommensurate krypton spacers of nearly monolayer coverage are studied within the van der Waals-corrected density functional theory. It is shown that when the spacer is in the commensurate phase, the graphene layers have the AA stacking. For this phase, the barriers to relative in-plane translational and rotational motion and the shear mode frequency of the graphene layers are calculated. For the incommensurate phase, both of the barriers are found to be negligibly small. A considerable change of tunneling conductance between the graphene layers separated by the commensurate krypton spacer at their relative subangstrom displacement is revealed by the use of the Bardeen method. The possibility of nanoelectromechanical systems based on the studied tribological and electronic properties of the considered heterostructures is discussed.

preprint2013arXiv

Ab initio study of edge effect on relative motion of walls in carbon nanotubes

Interwall interaction energies of double-walled nanotubes (DWNTs) with long inner and short outer walls are calculated as functions of coordinates describing relative rotation and displacement of the walls using van der Waals corrected density functional theory. The magnitude of corrugation and the shape of the potential energy relief are found to be very sensitive to changes of the shorter wall length at subnanometer scale and atomic structure of the edges if at least one of the walls is chiral. Threshold forces required to start relative motion of the short walls and temperatures at which the transition between diffusive and free motion of the short walls takes place are estimated. The edges are also shown to provide a considerable contribution to the barrier to relative rotation of commensurate nonchiral walls. For such walls, temperatures of orientational melting, i.e. the crossover from rotational diffusion to free relative rotation, are estimated. The possibility to produce nanotube-based bolt/nut pairs and nanobearings is discussed.

preprint2013arXiv

Approaches to modelling irradiation-induced processes in transmission electron microscopy

The recent progress in high-resolution transmission electron microscopy (HRTEM) has given rise to the possibility of in situ observations of nanostructure transformations and chemical reactions induced by electron irradiation. In this article we briefly summarise experimental observations and discuss in detail atomistic modelling of irradiation-induced processes in HRTEM, as well as mechanisms of such processes recognised due to modelling. Accurate molecular dynamics (MD) techniques based on first principles or tight-binding models are employed in the analysis of single irradiation-induced events, and classical MD simulations are combined with a kinetic Monte Carlo algorithm to simulate continuous irradiation of nanomaterials. It has been shown that sulphur-terminated graphene nanoribbons are formed inside carbon nanotubes as a result of an irradiation-selective chemical reaction. The process of fullerene formation in HRTEM during continuous electron irradiation of a small graphene flake has been simulated, and mechanisms driving this transformation analysed.

preprint2013arXiv

Structure, energetic and tribological properties, and possible applications in NEMS of argon-separated double-layer graphene

The possibility to control the commensurability and distance between graphene layers separated by a dielectric spacer is considered by the example of a heterostructure consisting of double-layer graphene separated by atomic layers of argon. Van der Waals corrected density functional theory (DFT-D) is applied to study structural, energetic and tribological characteristics of this heterostructure. It is found that in the ground state, monolayer and bilayer argon spacers are incommensurate with the graphene layers, whereas submonolayer argon spacers which are commensurate with the graphene layers can exist only as metastable states. The calculations show that this incommensurability provides negligibly small static friction for relative motion of argon-separated graphene layers and, therefore, such a heterostructure holds great promise for the use in nanoelectromechanical systems. A scheme and operational principles of the nanorelay based on the revealed tribological properties of the heterostructure are proposed.

preprint2012arXiv

Barriers to motion and rotation of graphene layers based on measurements of shear mode frequencies

Both van der Waals corrected density functional theory and classical calculations show that the potential relief of interaction energy between layers of graphite and few-layer graphene can be described by a simple expression containing only the first Fourier components. Thus a set of physical quantities and phenomena associated with in-plane relative vibration, translational motion and rotation of graphene layers are interrelated and are determined by a single parameter characterizing the roughness of the potential energy relief. This relationship is used to estimate the barriers to relative motion and rotation of graphene layers based on experimental measurements of shear mode frequencies.

preprint2012arXiv

Modeling of graphene-based NEMS

The possibility of designing nanoelectromechanical systems (NEMS) based on relative motion or vibrations of graphene layers is analyzed. Ab initio and empirical calculations of the potential relief of interlayer interaction energy in bilayer graphene are performed. A new potential based on the density functional theory calculations with the dispersion correction is developed to reliably reproduce the potential relief of interlayer interaction energy in bilayer graphene. Telescopic oscillations and small relative vibrations of graphene layers are investigated using molecular dynamics simulations. It is shown that these vibrations are characterized with small Q-factor values. The perspectives of nanoelectromechanical systems based on relative motion or vibrations of graphene layers are discussed.

preprint2012arXiv

Nanotube-based scanning rotational microscope

A scheme of the scanning rotational microscope is designed. This scheme is based on using carbon nanotubes simultaneously as a probe tip and as a bolt/nut pair which converts translational displacements of two piezo actuators into pure rotation of the probe tip. First-principles calculations of the interaction energy between movable and rotational parts of the microscope confirms the capability for its operation. The scanning rotational microscope with a chemically functionalized nanotube-based tip can be used to study how the interaction between individual molecules or a molecule and a surface depends on their relative orientation.

preprint2012arXiv

Ni-assisted transformation of graphene flakes to fullerenes

Transformation of graphene flakes to fullerenes assisted by Ni clusters is investigated using molecular dynamics simulations. The bond-order potential for Ni-C systems is developed. The potential reproduces the experimental and first-principles data on the physical properties of pure Ni as well as on relative energies of carbon species on Ni surfaces and in Ni bulk. The potential is applied for molecular dynamics simulations of the transformation of graphene flakes consisting of 50 - 400 atoms with and without Ni clusters attached. Free fullerenes, fullerenes with Ni clusters attached from outside and fullerenes encapsulating Ni clusters (Ni endofullerenes) are observed to form in the presence of Ni clusters consisting of 5 - 80 atoms. Moreover, a new type of heterofullerenes with a patch made of a Ni cluster is found to form as an intermediate structure during the transformation. The Ni clusters are shown to reduce the activation energy for the graphene-fullerene transformation from 4.0 eV to 1.5 - 1.9 eV, providing the decrease of the minimal temperature at which such a transformation can be observed experimentally from about 1400 K for free graphene flakes to about 700 - 800 K. While the transformation of free graphene flakes is found to occur through formation of chains of two-coordinated carbon atoms at the flake edges, the mechanism of the Ni-assisted graphene-fullerene transformation is revealed to be based on the transfer of carbon atoms from the graphene flake to the Ni cluster and back. The way of controlled synthesis of endofullerenes with a transition metal cluster inside and heterofullerenes with a transition metal patch is also proposed.

preprint2011arXiv

Commensurate-incommensurate phase transition in bilayer graphene

A commensurate-incommensurate phase transition in bilayer graphene is investigated in the framework of the Frenkel-Kontorova model extended to the case of two interacting chains of particles. Analytic expressions are derived to estimate the critical unit elongation of one of the graphene layers at which the transition to the incommensurate phase takes place, the length and formation energy of incommensurability defects (IDs) and the threshold force required to start relative motion of the layers on the basis of dispersion-corrected density functional theory calculations of the interlayer interaction energy as a function of the relative position of the layers. These estimates are confirmed by atomistic calculations using the DFT-D based classical potential. The possibility to measure the barriers for relative motion of graphene layers by the study of formation of IDs in bilayer graphene is discussed.

preprint2011arXiv

Diffusion and drift of graphene flake on graphite surface

Diffusion and drift of a graphene flake on a graphite surface are analyzed. A potential energy relief of the graphene flake is computed using ab initio and empirical calculations. Based on the analysis of this relief, different mechanisms of diffusion and drift of the graphene flake on the graphite surface are considered. A new mechanism of diffusion and drift of the flake is proposed. According to the proposed mechanism, rotational transition of the flake from commensurate to incommensurate state takes place with subsequent simultaneous rotation and translational motion until a commensurate state is reached again, and so on. Analytic expressions for the diffusion coefficient and mobility of the flake corresponding to different mechanisms are derived in wide ranges of temperatures and sizes of the flake. The molecular dynamics simulations and estimates based on ab initio and empirical calculations demonstrate that the proposed mechanism can be dominant under certain conditions. The influence of structural defects on the diffusion of the flake is examined on the basis of calculations of the potential energy relief and molecular dynamics simulations. The methods of control over the diffusion and drift of graphene components in nanoelectromechanical systems are discussed. The possibility to experimentally determine the barriers to relative motion of graphene layers based on the study of diffusion of a graphene flake is considered. The results obtained can be also applied to polycyclic aromatic molecules on graphene and should be qualitatively valid for a set of commensurate adsorbate-adsorbent systems.

preprint2011arXiv

Fast diffusion of graphene flake on graphene layer

Diffusion of a graphene flake on a graphene layer is analyzed and a new diffusion mechanism is proposed for the system under consideration. According to this mechanism, rotational transition of the flake from commensurate to incommensurate states takes place with subsequent simultaneous rotation and translational motion until the commensurate state is reached again, and so on. The molecular dynamics simulations and analytic estimates based on ab initio and semi-empirical calculations demonstrate that the proposed diffusion mechanism is dominant at temperatures T ~ Tcom, where Tcom corresponds to the barrier for transitions of the flake between adjacent energy minima in the commensurate states. For example, for the flake consisting of ~ 40, 200 and 700 atoms the contribution of the proposed diffusion mechanism through rotation of the flake to the incommensurate states exceeds that for diffusion of the flake in the commensurate states by one-two orders of magnitude at temperatures 50 - 150 K, 200 - 600 K and 800 - 2400 K, respectively. The possibility to experimentally measure the barriers to relative motion of graphene layers based on the study of diffusion of a graphene flake is considered. The results obtained are also relevant for understanding of dynamic behavior of polycyclic aromatic molecules on graphene and should be qualitatively valid for a set of commensurate adsorbate-adsorbent systems.

preprint2011arXiv

Molecular dynamics simulation of the self-retracting motion of a graphene flake

The self-retracting motion of a graphene flake on a stack of graphene flakes is studied using molecular dynamics simulations. It is shown that in the case when the extended flake is initially rotated to an incommensurate state, there is no barrier to the self-retracting motion of the flake and the flake retracts as fast as possible. If the extended flake is initially commensurate with the other flakes, the self-retracting motion is hindered by potential energy barriers. However, in this case, the rotation of the flake to incommensurate states is often observed. Such a rotation is found to be induced by the torque acting on the flake on hills of the potential relief of the interaction energy between the flakes. Contrary to carbon nanotubes, telescopic oscillations of the graphene flake are suppressed because of the high dynamic friction related to the excitation of flexural vibrations of the flake. This makes graphene promising for the use in fast-responding electromechanical memory cells.