Section of Theory and computer simulations in materials science and condensed matter physics

The critical value of the frequency as a function of driving amplitude A. The horizontal dashed line demonstrates the upper edge of the phonon spectrum. The green-shaded region below the curve corresponds to the frequency-amplitude range in which the supratransmission effect is observed. The solid lines connecting the data points are provided as visual guide.

Simulation of supratransmission effect in nickel

The critical value of the frequency as a function of driving amplitude A. The horizontal dashed line demonstrates the upper edge of the phonon spectrum. The green-shaded region below the curve corresponds to the frequency-amplitude range in which the supratransmission effect is observed. The solid lines connecting the data points are provided as visual guide.
D.V. Bachurin, R.T. Murzaev
PDF PDF  Article
Year: 2025 Volume: 15     Issue: 4 Pages: 383-387

Electronic band structures of graphene. Solid and dashed lines correspond to DFT and non-orthogonal tight-binding model calculations, respectively.

DFT-based parameterization of a non-orthogonal tight-binding model for electronic band structure calculations of carbon and hydrocarbon materials

Electronic band structures of graphene. Solid and dashed lines correspond to DFT and non-orthogonal tight-binding model calculations, respectively.
N.V. Kryakvin, V.A. Kurakin, T.N. Kobernik, M.M. Maslov
PDF PDF  Article
Year: 2025 Volume: 15     Issue: 4 Pages: 362-368

Comparison of MTP and DeePMD models for Fe-V alloys shows similar trends in the Cowley SRO parameter across vanadium concentrations, with both methods accurately capturing short-range order behavior. Structural snapshots illustrate atomic configurations at different V contents.

A comparative analysis of DeePMD and MTP models for short-range order prediction in Fe-V solid solutions

Comparison of MTP and DeePMD models for Fe-V alloys shows similar trends in the Cowley SRO parameter across vanadium concentrations, with both methods accurately capturing short-range order behavior. Structural snapshots illustrate atomic configurations at different V contents.
P.I. Chernovol, A.V. Verkhovykh, A.A. Mirzoev
PDF PDF  Article
Year: 2025 Volume: 15     Issue: 4 Pages: 331-335

Using the MD/MC Metropolis scheme, the most probable fracture-surface structures (FS)   formed on the cleavage plane were studied in simulations of the Al–Mg and Al–Ni alloys containing symmetric Σ5 boundaries with superstructures of segregations. The resulting FS structures were analyzed depending on the type of grain boundaries and the type of impurity atoms.

Simulation of fracture along grain boundaries containing superstructures of segregations in Al-Mg and Al-Ni systems

Using the MD/MC Metropolis scheme, the most probable fracture-surface structures (FS) formed on the cleavage plane were studied in simulations of the Al–Mg and Al–Ni alloys containing symmetric Σ5 boundaries with superstructures of segregations. The resulting FS structures were analyzed depending on the type of grain boundaries and the type of impurity atoms.
L.E. Karkina, I.N. Karkin
PDF PDF  Article
Year: 2025 Volume: 15     Issue: 4 Pages: 324-330


Combining first-principles calculations with experimental results, it is demonstrated that the middle layer of multi-layered aluminum sheets serves as an effective barrier, preventing the diffusion of Si into the core layer. This significantly enhances both the mechanical strength and corrosion resistance of the material.

Enhanced mechanical properties and corrosion resistance synergy of multi-layered aluminum sheets fabricated with interlayer

Combining first-principles calculations with experimental results, it is demonstrated that the middle layer of multi-layered aluminum sheets serves as an effective barrier, preventing the diffusion of Si into the core layer. This significantly enhances both the mechanical strength and corrosion resistance of the material.
Z. Yuan, Y. Lu, Y. Jiang, X. Wang, Z. Ni, Y. Tu, J. Peng, Z. Liu, Q. Zhang, P. Wang, Y. Peng
PDF PDF  Article
Year: 2025 Volume: 15     Issue: 4 Pages: 236-248


Molecular dynamics simulations demonstrate that silicene experiences a smooth decline in mechanical properties accompanied by the formation of inhomogeneous wrinklons as defect size increases, while graphene undergoes a sharp drop in strength at a critical defect size and exhibits uniform ripples.  These distinct deformation stages and wrinkle morphologies are directly linked to the atomic structures of silicene and graphene, influencing their mechanical responses under tension.

Deformation behaviour of silicene with dislocation dipoles under uniaxial tension

Molecular dynamics simulations demonstrate that silicene experiences a smooth decline in mechanical properties accompanied by the formation of inhomogeneous wrinklons as defect size increases, while graphene undergoes a sharp drop in strength at a critical defect size and exhibits uniform ripples. These distinct deformation stages and wrinkle morphologies are directly linked to the atomic structures of silicene and graphene, influencing their mechanical responses under tension.
A.K. Akhunova, R.T. Murzaev
PDF PDF  Article
Year: 2025 Volume: 15     Issue: 3 Pages: 207-212





The thermal conductivity properties of the graphene-nanotube film/Al2O3 interface are investigated using the molecular dynamics method. Both films and nanoparticles of Al2O3 are considered.

Thermal conductivity properties of the graphene-carbon nanotube hybrid / Al2O3 interface

The thermal conductivity properties of the graphene-nanotube film/Al2O3 interface are investigated using the molecular dynamics method. Both films and nanoparticles of Al2O3 are considered.
M.M. Slepchenkov, O.E. Glukhova, A.A. Petrunin
PDF PDF  Article
Year: 2025 Volume: 15     Issue: 2 Pages: 127-133


Using the molecular dynamics method, a study was conducted on the influence of the mean grain size and excess energy on the melting temperature of nanocrystalline aluminum.

Influence of mean grain size and excess energy on the melting temperature of nanocrystalline aluminum

Using the molecular dynamics method, a study was conducted on the influence of the mean grain size and excess energy on the melting temperature of nanocrystalline aluminum.
G.M. Poletaev, A.A. Sitnikov, Y.V. Bebikhov, A.S. Semenov
PDF PDF  Article
Year: 2025 Volume: 15     Issue: 2 Pages: 77-83

Density profile of the simulation box along the x-axis, illustrating the shift of the solid-liquid interface toward the solid region at temperatures above the melting point. The top image shows the original simulation box of the V-5Cr-5Ti alloy, while the bottom image displays the same simulation box with structure types identified using Common Neighbor Analysis.

Thermo-mechanical properties of random ternary alloy V-5Cr-5Ti

Density profile of the simulation box along the x-axis, illustrating the shift of the solid-liquid interface toward the solid region at temperatures above the melting point. The top image shows the original simulation box of the V-5Cr-5Ti alloy, while the bottom image displays the same simulation box with structure types identified using Common Neighbor Analysis.
H.S. M. Phuong, N.T. H. Trung, M.D. Starostenkov
PDF PDF  Article
Year: 2025 Volume: 15     Issue: 1 Pages: 36-42

The ultrasonic vibration of substrate with deposited layers was simulated. The dependencies of the beam deflection on the length and height of the wall
were studied

The simulation of the variation in the ultrasonic wave characteristics during the steel sample building by wire arc additive manufacturing

The ultrasonic vibration of substrate with deposited layers was simulated. The dependencies of the beam deflection on the length and height of the wall were studied
A.M. Ivanov, S.P. Belyaev, V.V. Rubanik, N.N. Resnina, V.V. Rubanik jr., R.M. Bikbaev, I.V. Ponikarova, A.V. Sibirev
PDF PDF  Article
Year: 2024 Volume: 14     Issue: 4 Pages: 399-405

Change in sp2-hybridization to sp3 results in an increase of Young's modulus and a decrease of Poisson's ratio for bornitrane. The substitution of carbon atoms by B and N atoms in the lattice of graphyne and graphene, leads to a decrease in Young's modulus.

Elastic constants of two-dimensional BN nanostructures

Change in sp2-hybridization to sp3 results in an increase of Young's modulus and a decrease of Poisson's ratio for bornitrane. The substitution of carbon atoms by B and N atoms in the lattice of graphyne and graphene, leads to a decrease in Young's modulus.
P.V. Polyakova
PDF PDF  Article
Year: 2024 Volume: 14     Issue: 4 Pages: 379-385



Using the molecular dynamics method, a study was conducted on the influence of structural defects, carbon and hydrogen impurities, as well as the presence of a shell of another metal, on the melting temperature of metallic nanoparticles.

Molecular dynamics study of the effect of structural defects, impurities, and the presence of a shell on the melting temperature of metallic nanoparticles

Using the molecular dynamics method, a study was conducted on the influence of structural defects, carbon and hydrogen impurities, as well as the presence of a shell of another metal, on the melting temperature of metallic nanoparticles.
G.M. Poletaev, A.V. Sannikov, Y.Y. Gafner, S.L. Gafner
PDF PDF  Article
Year: 2024 Volume: 14     Issue: 4 Pages: 332-339




The present work aims to provide the differences in elasticity of pentagonal wires and rods within the framework of distributed and single disclination models, employing both the analytical technique and the finite element method.

Elasticity of pentagonal wires: single disclination model versus distributed disclination model

The present work aims to provide the differences in elasticity of pentagonal wires and rods within the framework of distributed and single disclination models, employing both the analytical technique and the finite element method.
R.E. Shevchuk, S.A. Krasnitckii, A.E. Romanov, A.M. Smirnov
PDF PDF  Article
Year: 2024 Volume: 14     Issue: 3 Pages: 175-182

The influence of COOH and C=O groups on the I-V characteristics of graphene nanomesh. Depending on the type of group, one can either strengthen or weaken the effect of negative differential resistance on the I-V characteristic.

Current flow patterns in graphene nanomesh films functionalized with carbonyl and carboxyl groups

The influence of COOH and C=O groups on the I-V characteristics of graphene nanomesh. Depending on the type of group, one can either strengthen or weaken the effect of negative differential resistance on the I-V characteristic.
P.V. Barkov, M.M. Slepchenkov, O.E. Glukhova
PDF PDF  Article
Year: 2024 Volume: 14     Issue: 3 Pages: 210-215


Chemically controlled size distribution of point defect clusters generated by single displacement cascades determines amount and sizes of radiation defects under prolonged irradiation of materials

Chemically controlled radiation resistance of single-phase fcc Ni-Fe-Cr concentrated solid solutions

Chemically controlled size distribution of point defect clusters generated by single displacement cascades determines amount and sizes of radiation defects under prolonged irradiation of materials
A.V. Korchuganov, O.A. Berezikov
PDF PDF  Article
Year: 2024 Volume: 14     Issue: 1 Pages: 72-78

The effect of laser thermal action on metal grains in 316L alloy obtained by supersonic laser deposition have been investigated. The treat by laser provide growth metal grains by diffusion. The results of experimental studies correlate with the results of process modeling.

Phase-field modeling of changes in the grain structure of 316L steel obtained by cold spraying followed by laser treatment

The effect of laser thermal action on metal grains in 316L alloy obtained by supersonic laser deposition have been investigated. The treat by laser provide growth metal grains by diffusion. The results of experimental studies correlate with the results of process modeling.
M.S. Orlova, A.I. Gorunov
PDF PDF  Article
Year: 2024 Volume: 14     Issue: 1 Pages: 79-84


The numerical investigation of the behaviour of a Cu-graphene composite under harmonic impact revealed that if the graphene component is subjected to harmonic driving, it becomes the predominant channel for energy transport due to its wider phonon spectrum and greater number of phonons

Nonlinear supratransmission effect in copper-graphene composite

The numerical investigation of the behaviour of a Cu-graphene composite under harmonic impact revealed that if the graphene component is subjected to harmonic driving, it becomes the predominant channel for energy transport due to its wider phonon spectrum and greater number of phonons
A.M. Kazakov, P.V. Zakharov, G.F. Korznikova, Y.V. Bebikhov, E.A. Korznikova
PDF PDF  Article
Year: 2023 Volume: 13     Issue: 4s Pages: 499-504


The influence of the rate of change in temperature and size of nickel nanoparticles on the features of the nanocrystalline structure formed during crystallization and devitrification was studied by the method of molecular dynamics

Molecular dynamics study of the formation of a nanocrystalline structure during crystallization and devitrification of nickel nanoparticles

The influence of the rate of change in temperature and size of nickel nanoparticles on the features of the nanocrystalline structure formed during crystallization and devitrification was studied by the method of molecular dynamics
G.M. Poletaev, Y.Y. Gafner, S.L. Gafner, A.V. Sannikov, I.V. Zorya
PDF PDF  Article
Year: 2023 Volume: 13     Issue: 4s Pages: 420-425

1. A theoretical study of new autointercalated carbon compounds formed from hexagonal graphyne layers has been performed.
2. The most stable nanostructured materials should be composed of autointercalated alpha-graphyne-1.
3. A porous orthorombic three-dimensional phase can be obtained on the basis of alpha-graphyne-1.

Calculations of the structure and properties of autointercalated graphyne layers

1. A theoretical study of new autointercalated carbon compounds formed from hexagonal graphyne layers has been performed. 2. The most stable nanostructured materials should be composed of autointercalated alpha-graphyne-1. 3. A porous orthorombic three-dimensional phase can be obtained on the basis of alpha-graphyne-1.
V.A. Greshnyakov, V.V. Pavlik
PDF PDF  Article
Year: 2023 Volume: 13     Issue: 4 Pages: 323-328

3D view of topological charge phase diagram. Axes are parameters of exchange and magnetoelectric interaction, height is absolute value of topological charge in ground state.

Skyrmion phases in ground state of magnetoelectric bilayer induced by planar Dzyaloshinskii-Moriya interaction

3D view of topological charge phase diagram. Axes are parameters of exchange and magnetoelectric interaction, height is absolute value of topological charge in ground state.
A.G. Nugumanov, I.F. Sharafullin, M.K. Kharrasov
PDF PDF  Article
Year: 2023 Volume: 13     Issue: 4 Pages: 317-322

It has been established that the localization of crystalline LaB6 nanoparticles determines the work function of the hybrid carbon nanomaterial single-walled nanotubes/graphene. As the mass fraction of LaB6 nanoparticles increases, the work function decreases.

Effect of LaB6 nanoparticles on the electronic and emission properties of single-walled carbon nanotubes / graphene hybrid 1D nanomaterial

It has been established that the localization of crystalline LaB6 nanoparticles determines the work function of the hybrid carbon nanomaterial single-walled nanotubes/graphene. As the mass fraction of LaB6 nanoparticles increases, the work function decreases.
O.E. Glukhova, M.M. Slepchenkov
PDF PDF  Article
Year: 2023 Volume: 13     Issue: 4 Pages: 312-316



The impulse load wave pattern of the planar target with the simple shaped long-length impactor as a two-dimensional bar is presented. Needle-like spall fracture occurs under shock loading of the target surface with spherical particles.

Specific features of target pulse destruction under impact by an elongated solid

The impulse load wave pattern of the planar target with the simple shaped long-length impactor as a two-dimensional bar is presented. Needle-like spall fracture occurs under shock loading of the target surface with spherical particles.
S.N. Buravova, E.V. Petrov, V.O. Kopytskiy
PDF PDF  Article
Year: 2023 Volume: 13     Issue: 3 Pages: 232-236

New two-dimensional structure - diamane have been studied by molecular dynamics simulation. The stiffness constants are calculated from molecular dynamics and then applied to analize Young's modulus, shear modulus, and Poisson's ratio.

Elastic properties of diamane

New two-dimensional structure - diamane have been studied by molecular dynamics simulation. The stiffness constants are calculated from molecular dynamics and then applied to analize Young's modulus, shear modulus, and Poisson's ratio.
P.V. Polyakova, L.K. Galiakhmetova, R.T. Murzaev, D.S. Lisovenko, J.A. Baimova
PDF PDF  Article
Year: 2023 Volume: 13     Issue: 2 Pages: 171-176



Simulation of hydrides’ morphology in zirconium in 3D representation assuming classical heterogenous nucleation and diffusion growth with taking into account the mutual influence of neighbors.

Evolution of an approach to the modeling of zirconium hydrides morphology based on Monte-Carlo method in 3D representation

Simulation of hydrides’ morphology in zirconium in 3D representation assuming classical heterogenous nucleation and diffusion growth with taking into account the mutual influence of neighbors.
T.N. Aliev, M.Y. Kolesnik
PDF PDF  Review
Year: 2023 Volume: 13     Issue: 2 Pages: 143-148


1. Modified Likhachev-Volkov  model  may be successfully used for the simulation of the strain variation during cooling and isothermal holding under stress.
2. To describe the strain variation on holding after active deformation, elastic energy stored during the transformation should be decreasd.

Simulation of recoverable strain variation during isothermal holding of the Ni51Ti49 alloy under various regimes

1. Modified Likhachev-Volkov model may be successfully used for the simulation of the strain variation during cooling and isothermal holding under stress. 2. To describe the strain variation on holding after active deformation, elastic energy stored during the transformation should be decreasd.
N.N. Resnina, A.M. Ivanov, F.S. Belyaev, A.E. Volkov, S.P. Belyaev
PDF PDF  Article
Year: 2023 Volume: 13     Issue: 1 Pages: 33-38

Gravimetric density as the function of exposure time at two temperatures and pressure. The inset shows the initial structure of a crumpled graphene flake (gray atoms) in a hydrogen atmosphere (blue atoms).

Effect of external pressure on the hydrogen storage capacity of a graphene flake: molecular dynamics

Gravimetric density as the function of exposure time at two temperatures and pressure. The inset shows the initial structure of a crumpled graphene flake (gray atoms) in a hydrogen atmosphere (blue atoms).
N.G. Apkadirova, K.A. Krylova, J.A. Baimova
PDF PDF  Article
Year: 2022 Volume: 12     Issue: 4s Pages: 445-450

Modeling of the structure and electronic properties was carried out from the first principles for the polymorphic varieties of graphene with adsorbed hydroxyl groups (hydroxygraphene), where basic graphene layer structure consisted exclusively of 4-8 paired topological defects.

Modeling the structure and electronic properties of new polymorphic L4-8 hydroxygraphene varieties

Modeling of the structure and electronic properties was carried out from the first principles for the polymorphic varieties of graphene with adsorbed hydroxyl groups (hydroxygraphene), where basic graphene layer structure consisted exclusively of 4-8 paired topological defects.
M.E. Belenkov, V.A. Greshnyakov, V.M. Chernov
PDF PDF  Article
Year: 2022 Volume: 12     Issue: 4s Pages: 404-408

The synthesized material mainly consists of amorphous carbon with graphite crystallites, the maximum amount of which is formed during an arc discharge in nitrogen.

Experimental and theoretical study of the conditions for the formation of carbon nanostructures in an arc discharge in helium, argon and nitrogen

The synthesized material mainly consists of amorphous carbon with graphite crystallites, the maximum amount of which is formed during an arc discharge in nitrogen.
S.Z. Sakhapov, V.A. Andryushchenko, E.V. Boyko, M.S. Skirda, D.V. Smovzh
PDF PDF  Article
Year: 2022 Volume: 12     Issue: 4 Pages: 321-326


1 2 3 4 5  ...