Source author record

Alexander V. Savin

Alexander V. Savin 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

20works
7topics
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

20 published item(s)

preprint2023arXiv

Chiral organic molecular structures supported by multilayer surfaces

We study numerically the dynamics of acetanilide (ACN) molecules placed on a flat surface of a multilayer hexagonal boron nitride structure. We demonstrate that the ACN molecules, being achiral in three dimensions, become chiral after being placed on the substrate. Homochirality of the ACN molecules leads to stable secondary structures stabilized by hydrogen bonds between peptide groups of the molecules. Numerical simulations of systems of such molecules reveal that the structure of the resulting hydrogen-bond chains depends on the isomeric composition of the molecules. If all molecules are homochiral (i.e. only one isomer is present), they form secondary structures (chains of hydrogen bonds in the shapes of arcs, circles, and spirals). If the molecules at the substrate form a racemic mixture, then no regular secondary structures appear, and only curvilinear chains of hydrogen bonds of random shapes can emerge. A hydrogen-bond chain can form a straight zigzag only if it has an alternation of isomers. Such chains can create two-dimensional (2D) regular lattices, or 2D crystals. The melting scenarios of such 2D crystals depend on density of its coverage of the substrate. At 25% coverage, melting occurs continuously in a certain temperature interval. For a complete coverage, melting occurs at 415-470 K due to a shift of 11% of all molecules into the second layer of the substrate.

preprint2023arXiv

Localized nonlinear excitations of a columnar chain of coronene molecules

The nonlinear dynamics of a one-dimensional molecular crystal in the form of a chain of planar coronene molecules is analyzed. Using molecular dynamics, it is shown that a chain of coronene molecules supports acoustic solitons, rotobreathers, and discrete breathers. An increase in the size of planar molecules in a chain leads to an increase in the number of internal degrees of freedom. This results in an increase in the rate of emission of phonons from spatially localized nonlinear excitations and a decrease in their lifetime. Presented results contribute to the understanding of the effect of the rotational and internal vibrational modes of molecules on the nonlinear dynamics of molecular crystals.

preprint2022arXiv

Modeling of second sound in carbon nanostructures

The study of thermal transport in low-dimensional materials has attracted a lot of attention recently after discovery of high thermal conductivity of graphene. Here we study numerically phonon transport in low-dimensional carbon structures being interested in the hydrodynamic regime revealed through the observation of second sound. We demonstrate that correct numerical modeling of such two-dimensional systems requires semi-classical molecular dynamics simulations of temperature waves that take into account quantum statistics of thermalized phonons. We reveal that second sound can be attributed to the maximum group velocity of bending optical oscillations of carbon structures, and the hydrodynamic effects disappear for $T>200$K, being replaced by diffusive dynamics of thermal waves. Our numerical results suggest that the velocity of second sound in such low-dimensional structures is about 6 km/s, and the hydrodynamic effects are manifested stronger in carbon nanotubes rather than in carbon nanoribbons.

preprint2021arXiv

Multistability of carbon nanotube packings on flat substrate

It is shown by the method of molecular dynamics using a chain model that a multilayer packaging of identical single-walled carbon nanotubes with a diameter of D>2.5 nm located on a flat substrate is a multistable system. The system has many stationary states, which are characterized by the portion of collapsed nanotubes. The thickness of the package monotonically decreases with an increase in the portion of such nanotubes. For nanotubes with a chirality index (60,0), depending on the portion of collapsed nanotubes, the thickness of the 11-layer package can vary from 12 to 36 nm. All stationary states of the package are stable to thermal fluctuations at T=300K. The transverse compression of the package is not elastic; due to the collapse of a part of the nanotubes, it only transfers the package from one stationary state to another with a smaller thickness.

preprint2019arXiv

One-side surface modifications of graphene

We study the dependence of mechanical conformations of graphene sheets located on flat substrates on the density of unilateral (one-side) attachment of hydrogen, fluorine or chlorine atoms to them. It is shown that chemically modified graphene sheet can take four main forms on a flat substrate: the form of a flat sheet located parallel to the surface of the substrate, the form of convex sheet partially detached from the substrate with bent edges adjacent to the substrate, and the forms of single and double rolls on the substrate. On the surface of crystalline graphite, the flat form of the sheet is lowest in energy for hydrogenation density p <0.21, fluorination density p<0.20 and chlorination density p<0.16. The surface of crystalline nickel has higher adsorption energy for graphene monolayer and the flat form of chemically modified sheet on such substrate is lowest in energy for hydrogenation density p<0.47, fluorination density p<0.30 and chlorination density p<0.21. The flat shape of the graphene sheet remains basic on a substrate also when molecular groups CH3, CH2-CH3 or rings C6H5 are one-side attached to its outer surface. At the attachment density p=1/6 (one group per 6 sheet atoms) the sheet becomes the nanocarpet the basis of which is formed by a sheet of graphene and the pile of which is formed by the attached molecular groups forming a tightly packed regular lattice. The addition of hydroxyl groups OH with attachment density p=1/4 leads to the formation of hexagonal lattices of hydroxyl groups on the outer surface of graphene sheet on a substrate. In this lattice, the groups can form various configurations of hydrogen bonds, which turns the chemically modified sheet into a multistable system.

preprint2015arXiv

Confinement and non-universality of anomalous heat transport and superdiffusion of energy in low-dimensional systems

We provide molecular dynamics simulation of heat transport and thermal energy diffusion in one-dimensional molecular chains with different interparticle pair potentials at zero and non-zero temperature. We model the thermal conductivity (TC) and energy diffusion in the coupled rotator chain and in the Lennard-Jones chain either without or with the confining parabolic interatomic potential. The considered chains without the confining potential have normal TC and energy diffusion, while the corresponding chains with the confining potential are characterized by anomalous (diverging with the system length) TC and superdiffusion of energy. We confirm in such a way that, surprisingly, the confinement makes both heat transport and energy diffusion anomalous in low-dimensional phononic systems. We show that the chain, which has a finite TC, is also characterized by the normal energy diffusion in the thermalized chain (at non-zero temperature), while the superdiffusion of thermal energy occurs in the thermalized chains with only anomalous TC. We present the arguments, supported by our simulations, that the scaling relation between the exponents in time dependence of the mean square displacement of thermal energy distribution and in length dependence of the anomalous TC is not universal and can be different, depending on the main mechanism of energy transport: by weakly-scattered waves or by noninteracting colliding particles performing Levy flights.

preprint2015arXiv

Heat conduction in diatomic chains with correlated disorder

The paper considers heat transport in diatomic one-dimensional lattices, containing equal amounts of particles with different masses. Ordering of the particles in the chain is governed by single correlation parameter -- the probability for two neighboring particles to have the same mass. As this parameter grows from zero to unity, the structure of the chain varies from regular staggering chain to completely random configuration, and then -- to very long clusters of particles with equal masses. Therefore, this correlation parameter allows a control of typical cluster size in the chain. In order to explore different regimes of the heat transport, two interatomic potentials are considered. The first one is an infinite potential wall, corresponding to instantaneous elastic collisions between the neighboring particles. In homogeneous chains such interaction leads to an anomalous heat transport. The other one is classical Lennard-Jones interatomic potential, which leads to a normal heat transport. The simulations demonstrate that the correlated disorder of the particle arrangement does not change the convergence properties of the heat conduction coefficient, but essentially modifies its value. For the collision potential, one observes essential growth of the coefficient for fixed chain length as the limit of large homogeneous clusters is approached. The thermal transport in these models remains superdiffusive. In the Lennard-Jones chain the effect of correlation appears to be not monotonous in the limit of low temperatures. This behavior stems from the competition between formation of long clusters mentioned above, and Anderson localization close to the staggering ordered state.

preprint2015arXiv

Scroll configurations of carbon nanoribbons

Carbon nanoscroll is a unique topologically open configuration of graphene nanoribbon possessing outstanding properties and application perspectives due to its morphology. However molecular dynamics study of nanoscrolls with more than a few coils is limited by computational power. Here, we propose a simple model of the molecular chain moving in the plane, allowing to describe the folded and rolled packaging of long graphene nanoribbons. The model is used to describe a set of possible stationary states and the low-frequency oscillation modes of isolated single-layer nanoribbon scrolls as the function of the nanoribbon length. Possible conformational changes of scrolls due to thermal fluctuations are analyzed and their thermal stability is examined. Using the full-atomic model, frequency spectrum of thermal vibrations is calculated for the scroll and compared to that of the flat nanoribbon. It is shown that the density of phonon states of the scroll differs from the one of the flat nanoribbon only in the low (omega<100 cm-1) and high (omega>1450 cm-1) frequency ranges. Finally, the linear thermal expansion coefficient for the scroll outer radius is calculated from the long-term dynamics with the help of the developed planar chain model. The scrolls demonstrate anomalously high coefficient of thermal expansion and this property can find new applications.

preprint2013arXiv

Localized defect modes in graphene

We study the properties of localized vibrational modes associated with structural defects in a sheet of graphene. For the example of the Stone-Wales defects, one- and two-atom vacancies, many-atom linear vacancies, and adatoms in a honeycomb lattice, we demonstrate that the local defect modes are characterized by stable oscillations with the frequencies lying outside the linear frequency bands of an ideal graphene. In the frequency spectral density of thermal oscillations, such localized defect modes lead to the additional peaks from the right side of the frequency band of the ideal sheet of graphene. Thus, the general structure of the frequency spectral density can provide a fingerprint of its quality and the type of quantity of the structural defect a graphene sheet may contain.

preprint2013arXiv

Mechanical control of heat conductivity in microscopic models of dielectrics

We discuss a possibility to control a heat conductivity in simple one-dimensional models of dielectrics by means of external mechanical loads. To illustrate such possibilities we consider first a well-studied chain with degenerate double-well potential of the interparticle interaction. Contrary to previous studies, we consider varying length of the chain with fixed number of particles. Number of possible energetically degenerate ground states strongly depends on the overall length of the chain, or, in other terms, on average length of the link between neighboring particles. These degenerate states correspond to mechanical equilibrium, therefore one can say that the transition between them mimics to some extent a process of plastic deformation. We demonstrate that such modification of the chain length can lead to quite profound (almost five-fold) reduction of the heat conduction coefficient. Even more profound effect is revealed for a model with single-well non-convex potential. It is demonstrated that in certain range of constant external forcing this model becomes "effectively"\ double-well, and has a multitude of possible states of equilibrium for the same value of the external load. Thus, the heat conduction coefficient can be reduced by two orders of magnitude. We suggest a mechanical model of a chain with periodic double-well potential, which allows control over the heat conduction. The models considered may be useful for description of heat transport in biological macromolecules and for control of the heat transport in microsystems.

preprint2013arXiv

Two-phase stretching of molecular chains

While stretching of most polymer chains leads to rather featureless force-extension diagrams, some, notably DNA, exhibit non-trivial behavior with a distinct plateau region. Here we propose a unified theory that connects force-extension characteristics of the polymer chain with the convexity properties of the extension energy profile of its individual monomer subunits. Namely, if the effective monomer deformation energy as a function of its extension has a non-convex (concave up) region, the stretched polymer chain separates into two phases: the weakly and strongly stretched monomers. Simplified planar and 3D polymer models are used to illustrate the basic principles of the proposed model. Specifically, we show rigorously that when the secondary structure of a polymer is mostly due to weak non-covalent interactions, the stretching is two-phase, and the force-stretching diagram has the characteristic plateau. We then use realistic coarse-grained models to confirm the main findings and make direct connection to the microscopic structure of the monomers. We demostrate in detail how the two-phase scenario is realized in the α-helix, and in DNA double helix. The predicted plateau parameters are consistent with single molecules experiments. Detailed analysis of DNA stretching demonstrates that breaking of Watson-Crick bonds is not necessary for the existence of the plateau, although some of the bonds do break as the double-helix extends at room temperature. The main strengths of the proposed theory are its generality and direct microscopic connection.

preprint2012arXiv

Semi-quantum approach to molecular dynamics simulation of thermal properties of low-dimensional nanostructures

We present a detailed description of semi-quantum molecular dynamics simulation of stochastic dynamics of a system of interacting particles. Within this approach, the dynamics of the system is described with the use of classical Newtonian equations of motion in which the effects of phonon quantum statistics are introduced through random Langevin-like forces with a specific power spectral density (the color noise). The color noise describes the interaction of the molecular system with the thermostat. We apply this technique to the simulation of thermal properties and heat transport in different low-dimensional nanostructures. We describe the determination of temperature in quantum lattice systems, to which the equipartition limit is not applied. We show that one can determine the temperature of such system from the measured power spectrum and temperature- and relaxation-rate-independent density of vibrational (phonon) states. We simulate the specific heat and heat transport in carbon nanotubes, as well as the heat transport in molecular nanoribbons with perfect (atomically smooth) and rough (porous) edges, and in nanoribbons with strongly anharmonic periodic interatomic potentials. We show that the effects of quantum statistics of phonons are essential for the carbon nanotube in the whole temperature range T<500K, in which the values of the specific heat and thermal conductivity of the nanotube are considerably less than that obtained within the description based on classical statistics of phonons.

preprint2012arXiv

Transport of fullerene molecules along graphene nanoribbons

We study the motion of C60 fullerene molecules (buckyballs) and short-length carbon nanotubes on graphene nanoribbons. We demonstrate that the nanoribbon edge creates an effective potential that keeps the carbon structures on the surface. We reveal that the character of the motion of C60 molecules depends on temperature: for low temperatures (T<150K) the main type of motion is sliding along the surface, but for higher temperatures the sliding is replaced by rocking and rolling. Modeling of the buckyball with an included metal ion, such as Fe@C60, demonstrates that this molecular complex undergoes a rolling motion along the nanoribbon with the constant velocity under the action of a constant electric field. The similar effect is observed in the presence of the heat gradient applied to the nanoribbon, but mobility of carbon structures in this case depends largely on their size and symmetry, such that larger and more asymmetric structures demonstrate much lower mobility. Our results suggest that both electorphoresis and thermophoresis can be employed to control the motion of carbon molecules and fullerenes and, for example, sort them by their size, shape, and possible inclusions.

preprint2011arXiv

Heat conductivity of DNA double helix

Thermal conductivity of isolated single molecule DNA fragments is of importance for nanotechnology, but has not yet been measured experimentally. Theoretical estimates based on simplified (1D) models predict anomalously high thermal conductivity. To investigate thermal properties of single molecule DNA we have developed a 3D coarse-grained (CG) model that retains the realism of the full all-atom description, but is significantly more efficient. Within the proposed model each nucleotide is represented by 6 particles or grains; the grains interact via effective potentials inferred from classical molecular dynamics (MD) trajectories based on a well-established all-atom potential function. Comparisons of 10 ns long MD trajectories between the CG and the corresponding all-atom model show similar root-mean-square deviations from the canonical B-form DNA, and similar structural fluctuations. At the same time, the CG model is 10 to 100 times faster depending on the length of the DNA fragment in the simulation. Analysis of dispersion curves derived from the CG model yields longitudinal sound velocity and torsional stiffness in close agreement with existing experiments. The computational efficiency of the CG model makes it possible to calculate thermal conductivity of a single DNA molecule not yet available experimentally. For a uniform (polyG-polyC) DNA, the estimated conductivity coefficient is 0.3 W/mK which is half the value of thermal conductivity for water. This result is in stark contrast with estimates of thermal conductivity for simplified, effectively 1D chains ("beads on a spring") that predict anomalous (infinite) thermal conductivity. Thus, full 3D character of DNA double-helix retained in the proposed model appears to be essential for describing its thermal properties at a single molecule level.

preprint2011arXiv

Nonlinear Breathing-like Localized Modes in C60 Nanocrystals

We study the dynamics of nanocrystals composed of C60 fullerene molecules. We demonstrate that such structures can support long-lived strongly localized nonlinear oscillatory modes, which resemble discrete breathers in simple lattices. We reveal that at room temperatures the lifetime of such nonlinear localized modes may exceed tens of picoseconds; this suggests that C60 nanoclusters should demonstrate anomalously slow thermal relaxation when the temperature gradient decays in accord to a power law, thus violating the Cattaneo-Vernotte law of thermal conductivity.

preprint2011arXiv

The simplest model of polymer crystal exhibiting polymorphism

Almost all the polymer crystals have several polymorphic modifications. Their structure and existence conditions, as well as transitions between them are not understood even in the case of the 'model' polymer polyethylene (PE). For analysis of polymorphism in polymer crystals, we consider the simplest possible model of polymer chain: an extended flat zigzag made of 'united' atoms (replacing CH2-groups in PE chain); the united atoms belonging to different zigzags interact via Lennard-Jones potential. Analysis of potential of interaction between such zigzags allowed to predict the structure of five possible equilibrium lattices in polymer crystal built out of such zigzags. Molecular dynamics simulation of the crystal built out of flexible zigzags showed that, depending on model parameters (dimensions of the zigzag and equilibrium distance of Lennard-Jones potential), one to three of these lattices are stable in bulk at low temperatures. We have determined the model parameters at which the existing stable lattices are analogous to the ones observed in real PE and linear alkanes. The triclinic lattice has the lowest potential energy, then follows the monoclinic lattice, and the orthorhombic lattice has the highest energy, exactly as in full atomic (with hydrogens) molecular dynamics models.

preprint2010arXiv

Suppression of thermal conductivity in graphene nanoribbons with rough edges

We analyze numerically the thermal conductivity of carbon nanoribbons with ideal and rough edges. We demonstrate that edge disorder can lead to a suppression of thermal conductivity by several orders of magnitude. This effect is associated with the edge-induced Anderson localization and suppression of the phonon transport, and it becomes more pronounced for longer nanoribbons and low temperatures.

preprint2010arXiv

Vibrational Tamm states at the edges of graphene nanoribbons

We study vibrational states localized at the edges of graphene nanoribbons. Such surface oscillations can be considered as a phonon analog of Tamm states well known in the electronic theory. We consider both armchair and zigzag graphene stripes and demonstrate that surface modes correspond to phonons localized at the edges of the graphene nanoribbon, and they can be classified as in-plane and out-of-plane modes. In addition, in armchair nanoribbons anharmonic edge modes can experience longitudinal localization in the form of self-localized nonlinear modes, or surface breather solitons.

preprint2009arXiv

Non-stationary heat conduction in one-dimensional chains with conserved momentum

The Letter addresses the relationship between hyperbolic equations of heat conduction and microscopic models of dielectrics. Effects of the non-stationary heat conduction are investigated in two one-dimensional models with conserved momentum: Fermi-Pasta-Ulam (FPU) chain and chain of rotators (CR). These models belong to different universality classes with respect to stationary heat conduction. Direct numeric simulations reveal in both models a crossover from oscillatory decay of short-wave perturbations of the temperature field to smooth diffusive decay of the long-wave perturbations. Such behavior is inconsistent with parabolic Fourier equation of the heat conduction. The crossover wavelength decreases with increase of average temperature in both models. For the FPU model the lowest order hyperbolic Cattaneo-Vernotte equation for the non-stationary heat conduction is not applicable, since no unique relaxation time can be determined.

preprint2009arXiv

Reduced Thermal Conductivity of Nanowires and Nanoribbons with Dynamically Rough Surfaces and the "Problem of One-Dimensional Heat Conductors"

We present analytical model and molecular dynamics simulations of phonon heat transport in nanowires and nanoribbons with anharmonic lattices and dynamically rough surfaces and edges. In agreement with recent experiments on heat transport in single-crystalline silicon nanowires with rough surfaces, our model and simulations predict finite and length-independent phonon thermal conductivity in such quasi-one-dimensional systems, in contrast to anomalous phonon thermal conductivity of corresponding momentum-conserving systems with atomically smooth surfaces, divergent with the system length. Within our model, the main cause of thermal conductivity reduction is momentum-nonconserving scattering of longitudinal acoustic phonons by anharmonic side phonon leads in quasi-one-dimensional phonon waveguide with dynamically rough surface or edge layers.