Source author record

Alexander V. Shapeev

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

24works
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

24 published item(s)

preprint2022arXiv

Mechanical, optical, and thermoelectric properties of semiconducting ZnIn2X4 (X= S, Se, Te) monolayers

Mechanical stability of the ZnIn2X4 monolayers. The ZnIn2S4 and ZnIn2Se4 are semiconductors with direct band gaps of 3.94 and 2.77 eV, respectively whereas the ZnIn2Te4 shows an indirect band gap of 1.84 eV at the G0W0 level. The optical properties achieved from the solution of the Bethe-Salpeter equation predict the exciton binding energy of the ZnIn2S4, ZnIn2Se4, and ZnIn2Te4 monolayers to be 0.51, 0.41, and 0.34 eV, respectively, suggesting the high stability of the excitonic states against thermal dissociation. Using the iterative solutions of the Boltzmann transport equation accelerated by machine learning interatomic potentials, the room-temperature lattice thermal conductivity of the ZnIn2S4, ZnIn2Se4, and ZnIn2Te4 monolayers is predicted to be remarkably low as 5.8, 2.0, and 0.4 W/mK, respectively. Due to the low lattice thermal conductivity, high thermopower, and large figure of merit, we propose the ZnIn2Se4 and ZnIn2Te4 monolayers as promising candidates for thermoelectric energy conversion systems. This study provides an extensive vision concerning the intrinsic physical properties of the ZnIn2X4 nanosheets and highlights their characteristics for energy conversion and optoelectronics applications.

preprint2021arXiv

Validation of moment tensor potentials for fcc and bcc metals using EXAFS spectra

Machine-learning potentials for materials, namely the moment tensor potentials (MTPs), were validated using experimental EXAFS spectra for the first time. The MTPs for four metals (bcc W and Mo, fcc Cu and Ni) were obtained by the active learning algorithm of fitting to the results of the calculations using density functional theory (DFT). The MTP accuracy was assessed by comparing metal K-edge EXAFS spectra obtained experimentally and computed from the results of molecular dynamics (MD) simulations. The sensitivity of the method to various aspects of the MD and DFT models was demonstrated using Ni as an example. Good agreement was found for W, Mo and Cu using the recommended PAW pseudopotentials, whereas a more accurate pseudopotential with 18 valence electrons was required for Ni to achieve a similar agreement. The use of EXAFS spectra allows one to estimate the MTP ability in reproducing both average and dynamic atomic structures.

preprint2020arXiv

Ab initio analysis of some Ge-based 2D nanomaterials

The structural, electronic and dynamical properties of a group of 2D germanium-based compounds, including GeC, GeN, GeO, GeSi, GeS, GeSe, and germanene, are investigated by employing first-principles calculations. The most stable structure of each of these systems is identified after considering the most probable configurations and performing accurate phonon calculations. We introduce a new phase of germanene, which we name the tile germanene, which is significantly more stable than the known hexagonal germanene. We apply the modern modified Becke-Johnson (mBJ) and DFT1/2 schemes to obtain an accurate band structure for our selected 2D materials. It is seen that GeO and GeC exhibit the highest band gaps of more than 3 eV in this group of materials. Moreover, we argue that, in contrast to the semi-metallic nature of hexagonal germanene, the tile germanene is a very good conductor. The band edges of our semiconducting 2D materials are accurately aligned to the vacuum level to address the potential photocatalytic application of this system for water splitting and carbon dioxide reduction. The optical properties, including dielectric functions, refractive index, reflectivity, and Loss function of the samples are investigated in the framework of the Bethe-Salpeter approach.

preprint2020arXiv

Accelerating first-principles estimation of thermal conductivity by machine-learning interatomic potentials: A MTP/ShengBTE solution

Accurate evaluation of the thermal conductivity of a material can be a challenging task from both experimental and theoretical points of view. In particular for the nanostructured materials, the experimental measurement of thermal conductivity is associated with diverse sources of uncertainty. As a viable alternative to experiment, the combination of density functional theory (DFT) simulations and the solution of Boltzmann transport equation is currently considered as the most trusted approach to examine thermal conductivity. The main bottleneck of the aforementioned method is to acquire the anharmonic interatomic force constants using the computationally demanding DFT calculations. In this work we propose a substantially accelerated approach for the evaluation of anharmonic interatomic force constants via employing machine-learning interatomic potentials (MLIPs) trained over short ab-initio molecular dynamics trajectories. The remarkable accuracy of the proposed accelerated method is confirmed by comparing the estimated thermal conductivities of several bulk and two-dimensional materials with those computed by the full-DFT approach. The MLIP-based method proposed in this study can be employed as a standard tool, which would substantially accelerate and facilitate the estimation of lattice thermal conductivity in comparison with the commonly used full-DFT solution.

preprint2020arXiv

Exceptional piezoelectricity, high thermal conductivity and stiffness and promising photocatalysis in two-dimensional MoSi2N4 family confirmed by first-principles

Chemical vapor deposition has been most recently employed to fabricate centimeter-scale high-quality single-layer MoSi2N4 (Science; 2020;369; 670). Motivated by this exciting experimental advance, herein we conduct extensive first-principles based simulations to explore the stability, mechanical properties, lattice thermal conductivity, piezoelectric and flexoelectric response, and photocatalytic and electronic features of MA2Z4 (M = Cr, Mo, W; A = Si, Ge; Z = N, P) monolayers. The considered nanosheets are found to exhibit dynamical stability and remarkably high mechanical properties. Moreover, they show diverse electronic properties from antiferromagnetic metal to half metal and to semiconductors with band gaps ranging from 0.31 to 2.57 eV. Among the studied nanosheets, the MoSi2N4 and WSi2N4 monolayers yield appropriate band edge positions, high electron and hole mobilities, and strong visible light absorption, highly promising for applications in optoelectronics and photocatalytic water splitting. The MoSi2N4 and WSi2N4 monolayers are also predicted to show outstandingly high lattice thermal conductivity of 440 and 500 W/mK, respectively. For the first time we show that machine learning interatomic potentials trained over small supercells can be employed to examine the flexoelectric and piezoelectric properties of complex structures. As the most exciting finding, WSi2N4, CrSi2N4 and MoSi2N4 are found to exhibit the highest piezoelectric coefficients, outperforming all other-known 2D materials. Our results highlight that MA2Z4 nanosheets not only undoubtedly outperform the transition metal dichalcogenides family but also can compete with graphene for applications in nanoelectronics, optoelectronic, energy storage/conversion and thermal management systems.

preprint2020arXiv

Exploring Phononic Properties of Two-Dimensional Materials using Machine Learning Interatomic Potentials

Phononic properties are commonly studied by calculating force constants using the density functional theory (DFT) simulations. Although DFT simulations offer accurate estimations of phonon dispersion relations or thermal properties, but for low-symmetry and nanoporous structures the computational cost quickly becomes very demanding. Moreover, the computational setups may yield nonphysical imaginary frequencies in the phonon dispersion curves, impeding the assessment of phononic properties and the dynamical stability of the considered system. Here, we compute phonon dispersion relations and examine the dynamical stability of a large ensemble of novel materials and compositions. We propose a fast and convenient alternative to DFT simulations which derived from machine-learning interatomic potentials passively trained over computationally efficient ab-initio molecular dynamics trajectories. Our results for diverse two-dimensional (2D) nanomaterials confirm that the proposed computational strategy can reproduce fundamental thermal properties in close agreement with those obtained via the DFT approach. The presented method offers a stable, efficient, and convenient solution for the examination of dynamical stability and exploring the phononic properties of low-symmetry and porous 2D materials.

preprint2020arXiv

High thermal conductivity in semiconducting Janus and non-Janus diamanes

Most recently, F-diamane monolayer was experimentally realized by the fluorination of bilayer graphene. In this work we elaborately explore the electronic and thermal conductivity responses of diamane lattices with homo or hetero functional groups, including: non-Janus C2H, C2F and C2Cl diamane and Janus counterparts of C4HF, C4HCl and C4FCl. Noticeably, C2H, C2F, C2Cl, C4HF, C4HCl and C4FCl diamanes are found to show electronic diverse band gaps of, 3.86, 5.68, 2.42, 4.17, 0.86, and 2.05 eV, on the basis of HSE06 method estimations. The thermal conductivity of diamane nanosheets was acquired using the full iterative solutions of the Boltzmann transport equation, with substantially accelerated calculations by employing machine-learning interatomic potentials in obtaining the anharmonic force constants. According to our results, the room temperature lattice thermal conductivity of graphene and C2H, C2F, C2Cl, C4HF, C4HCl and C4FCl diamane monolayers are estimated to be 3636, 1145, 377, 146, 454, 244 and 196 W/mK, respectively. The underlying mechanisms resulting in significant effects of functional groups on the thermal conductivity of diamane nanosheets were thoroughly explored. Our results highlight the substantial role of functional groups on the electronic and thermal conduction responses of diamane nanosheets.

preprint2020arXiv

Machine-Learning Interatomic Potentials Enable First-Principles Multiscale Modeling of Lattice Thermal Conductivity in Graphene/Borophene Heterostructures

One of the ultimate goals of computational modeling in condensed matter is to be able to accurately compute materials properties with minimal empirical information. First-principles approaches such as the density functional theory (DFT) provide the best possible accuracy on electronic properties but they are limited to systems up to a few hundreds, or at most thousands of atoms. On the other hand, classical molecular dynamics (CMD) simulations and finite element method (FEM) are extensively employed to study larger and more realistic systems, but conversely depend on empirical information. Here, we show that machine-learning interatomic potentials (MLIPs) trained over short ab-initio molecular dynamics trajectories enable first-principles multiscale modeling, in which DFT simulations can be hierarchically bridged to efficiently simulate macroscopic structures. As a case study, we analyze the lattice thermal conductivity of coplanar graphene/borophene heterostructures, recently synthesized experimentally (Sci. Adv. 2019; 5: eaax6444), for which no viable classical modeling alternative is presently available. Our MLIP-based approach can efficiently predict the lattice thermal conductivity of graphene and borophene pristine phases, the thermal conductance of complex graphene/borophene interfaces and subsequently enable the study of effective thermal transport along the heterostructures at continuum level. This work highlights that MLIPs can be effectively and conveniently employed to enable first-principles multiscale modeling via hierarchical employment of DFT/CMD/FEM simulations, thus expanding the capability for computational design of novel nanostructures.

preprint2020arXiv

Nanoporous C3N4, C3N5 and C3N6 nanosheets; Novel strong semiconductors with low thermal conductivities and appealing optical/electronic properties

Carbon nitride two-dimensional (2D) materials are among the most attractive class of nanomaterials, with wide range of application prospects. As a continuous progress, most recently, two novel carbon nitride 2D lattices of C3N5 and C3N4 have been successfully experimentally realized. Motivated by these latest accomplishments and also by taking into account the well-known C3N4 triazine-based graphitic carbon nitride structures, we predicted two novel C3N6 and C3N4 counterparts. We then conducted extensive density functional theory simulations to explore the thermal stability, mechanical, electronic and optical properties of these novel nanoporous carbon-nitride nanosheets. According to our results all studied nanosheets were found to exhibit desirable thermal stability and mechanical properties. Non-equilibrium molecular dynamics simulations on the basis of machine learning interatomic potentials predict ultralow thermal conductivities for these novel nanosheets. Electronic structure analyses confirm direct band gap semiconducting electronic character and optical calculations reveal the ability of these novel 2D systems to adsorb visible range of light. Extensive first-principles based results by this study provide a comprehensive vision on the stability, mechanical, electronic and optical responses of C3N4, C3N5 and C3N6 as novel 2D semiconductors and suggest them as promising candidates for the design of advanced nanoelectronics and energy storage/conversion systems.

preprint2020arXiv

The MLIP package: Moment Tensor Potentials with MPI and Active Learning

The subject of this paper is the technology (the "how") of constructing machine-learning interatomic potentials, rather than science (the "what" and "why") of atomistic simulations using machine-learning potentials. Namely, we illustrate how to construct moment tensor potentials using active learning as implemented in the MLIP package, focusing on the efficient ways to sample configurations for the training set, how expanding the training set changes the error of predictions, how to set up ab initio calculations in a cost-effective manner, etc. The MLIP package (short for Machine-Learning Interatomic Potentials) is available at https://mlip.skoltech.ru/download/.

preprint2019arXiv

A Performance and Cost Assessment of Machine Learning Interatomic Potentials

Machine learning of the quantitative relationship between local environment descriptors and the potential energy surface of a system of atoms has emerged as a new frontier in the development of interatomic potentials (IAPs). Here, we present a comprehensive evaluation of ML-IAPs based on four local environment descriptors --- Behler-Parrinello symmetry functions, smooth overlap of atomic positions (SOAP), the Spectral Neighbor Analysis Potential (SNAP) bispectrum components, and moment tensors --- using a diverse data set generated using high-throughput density functional theory (DFT) calculations. The data set comprising bcc (Li, Mo) and fcc (Cu, Ni) metals and diamond group IV semiconductors (Si, Ge) is chosen to span a range of crystal structures and bonding. All descriptors studied show excellent performance in predicting energies and forces far surpassing that of classical IAPs, as well as predicting properties such as elastic constants and phonon dispersion curves. We observe a general trade-off between accuracy and the degrees of freedom of each model, and consequently computational cost. We will discuss these trade-offs in the context of model selection for molecular dynamics and other applications.

preprint2019arXiv

Ring Polymer Molecular Dynamics and Active Learning of Moment Tensor Potential for Gas-Phase Barrierless Reactions: Application to S + H2

Ring polymer molecular dynamics (RPMD) has proven to be an accurate approach for calculating thermal rate coefficients of various chemical reactions. For wider application of this methodology, efficient ways to generate the underlying full-dimensional potential energy surfaces (PESs) and the corresponding energy gradients are required. Recently, we have proposed a fully automated procedure based on combining the original RPMDrate code with active learning for PES on-the-fly using moment tensor potential and successfully applied it to two representative thermally activated chemical reactions [I. S. Novikov, Y. V. Suleimanov, A. V. Shapeev, Phys. Chem. Chem. Phys., 29503-29512 (2018)]. In this work, using a prototype insertion chemical reaction S + H$_2$, we show that this procedure works equally well for another class of chemical reactions. We find that the corresponding PES can be generated by fitting to less than 1500 automatically generated structures while the RPMD rate coefficients show deviation from the reference values within the typical convergence error of RPMDrate. We note that more structures are accumulated during the real-time propagation of the dynamic factor (the recrossing factor) as opposed to the previous study. We also observe that relatively flat free energy profile of the along the reaction coordinate before entering the complex-formation well can cause issues with locating the maximum of the free energy surface for less converged PESs. However, the final RPMD rate coefficient is independent of the position of the dividing surface that makes it invulnerable to this problem, keeping the total number of necessary structures within a few thousand. Our work concludes that, in future, the proposed methodology can be applied to realistic complex chemical reactions with various energy profiles.

preprint2014arXiv

Accuracy of computation of crystalline defects at finite temperature

The present paper aims at developing a theory of computation of crystalline defects at finite temperature. In a one-dimensional setting we introduce Gibbs distributions corresponding to such defects and rigorously establish their asymptotic expansion. We then give an example of using such asymptotic expansion to compare the accuracy of computations using the free boundary conditions and using an atomistic-to-continuum coupling method.

preprint2014arXiv

Analysis of an optimization-based atomistic-to-continuum coupling method for point defects

We formulate and analyze an optimization-based Atomistic-to-Continuum (AtC) coupling method for problems with point defects. Near the defect core the method employs a potential-based atomistic model, which enables accurate simulation of the defect. Away from the core, where site energies become nearly independent of the lattice position, the method switches to a more efficient continuum model. The two models are merged by minimizing the mismatch of their states on an overlap region, subject to the atomistic and continuum force balance equations acting independently in their domains. We prove that the optimization problem is well-posed and establish error estimates.

preprint2014arXiv

Analysis of Blended Atomistic/Continuum Hybrid Methods

We present a comprehensive error analysis of two prototypical atomistic-to-continuum coupling methods of blending type: the energy-based and the force-based quasicontinuum methods. Our results are valid in two and three dimensions, for finite range many-body interactions (e.g., EAM type), and in the presence of lattice defects (we consider point defects and dislocations). The two key ingredients in the analysis are (i) new force and energy consistency error estimates; and (ii) a new technique for proving energy norm stability of a/c couplings that requires only the assumption that the exact atomistic solution is a stable equilibrium.

preprint2013arXiv

An Optimization-Based Atomistic-to-Continuum Coupling Method

We present a new optimization-based method for atomistic-to-continuum (AtC) coupling. The main idea is to cast the coupling of the atomistic and continuum models as a constrained optimization problem with virtual Dirichlet controls on the interfaces between the atomistic and continuum subdomains. The optimization objective is to minimize the error between the atomistic and continuum solutions on the overlap between the two subdomains, while the atomistic and continuum force balance equations provide the constraints. Splitting of the atomistic and continuum problems instead of blending them and their subsequent use as constraints in the optimization problem distinguishes our approach from the existing AtC formulations. We present and analyze the method in the context of a one-dimensional chain of atoms modeled using a linearized two-body next-nearest neighbor interactions.

preprint2013arXiv

Development of an Optimization-Based Atomistic-to-Continuum Coupling Method

Atomistic-to-Continuum (AtC) coupling methods are a novel means of computing the properties of a discrete crystal structure, such as those containing defects, that combine the accuracy of an atomistic (fully discrete) model with the efficiency of a continuum model. In this note we extend the optimization-based AtC, formulated in arXiv:1304.4976 for linear, one-dimensional problems to multi-dimensional settings and arbitrary interatomic potentials. We conjecture optimal error estimates for the multidimensional AtC, outline an implementation procedure, and provide numerical results to corroborate the conjecture for a 1D Lennard-Jones system with next-nearest neighbor interactions.

preprint2013arXiv

Theory-based Benchmarking of the Blended Force-Based Quasicontinuum Method

We formulate an atomistic-to-continuum coupling method based on blending atomistic and continuum forces. Our precise choice of blending mechanism is informed by theoretical predictions. We present a range of numerical experiments studying the accuracy of the scheme, focusing in particular on its stability. These experiments confirm and extend the theoretical predictions, and demonstrate a superior accuracy of B-QCF over energy-based blending schemes.

preprint2012arXiv

A priori and a posteriori $W^{1,\infty}$ error analysis of a QC method for complex lattices

In this paper we prove a priori and a posteriori error estimates for a multiscale numerical method for computing equilibria of multilattices under an external force. The error estimates are derived in a $W^{1,\infty}$ norm in one space dimension. One of the features of our analysis is that we establish an equivalent way of formulating the coarse-grained problem which greatly simplifies derivation of the error bounds (both, a priori and a posteriori). We illustrate our error estimates with numerical experiments.

preprint2012arXiv

Consistent Energy-based Atomistic/Continuum Coupling for Two-body Potentials in Three Dimensions

Very few works exist to date on development of a consistent energy-based coupling of atomistic and continuum models of materials in more than one dimension. The difficulty in constructing such a coupling consists in defining a coupled energy whose minimizers are free from uncontrollable errors on the atomistic/continuum interface. In this paper a consistent coupling in three dimensions is proposed. The main achievement of this work is to identify and efficiently treat a modified Cauchy-Born continuum model which can be coupled to the exact atomistic model. The convergence and stability of the method is confirmed with numerical tests.

preprint2012arXiv

Numerical Methods for Multilattices

Among the efficient numerical methods based on atomistic models, the quasicontinuum (QC) method has attracted growing interest in recent years. The QC method was first developed for crystalline materials with Bravais lattice and was later extended to multilattices (Tadmor et al, 1999). Another existing numerical approach to modeling multilattices is homogenization. In the present paper we review the existing numerical methods for multilattices and propose another concurrent macro-to-micro method in the numerical homogenization framework. We give a unified mathematical formulation of the new and the existing methods and show their equivalence. We then consider extensions of the proposed method to time-dependent problems and to random materials.

preprint2011arXiv

Analysis of an Energy-based Atomistic/Continuum Coupling Approximation of a Vacancy in the 2D Triangular Lattice

We present a comprehensive a priori error analysis of a practical energy based atomistic/continuum coupling method (Shapeev, arXiv:1010.0512) in two dimensions, for finite-range pair-potential interactions, in the presence of vacancy defects. The majority of the work is devoted to the analysis of consistency and stability of the method. These yield a priori error estimates in the H1-norm and the energy, which depend on the mesh size and the "smoothness" of the atomistic solution in the continuum region. Based on these error estimates, we present heuristics for an optimal choice of the atomistic region and the finite element mesh, which yields convergence rates in terms of the number of degrees of freedom. The analytical predictions are supported by extensive numerical tests.

preprint2011arXiv

Consistent Energy-Based Atomistic/Continuum Coupling for Two-Body Potentials in One and Two Dimensions

This paper addresses the problem of consistent energy-based coupling of atomistic and continuum models of materials, limited to zero-temperature statics of simple crystals. It has been widely recognized that the most practical coupled methods exhibit large errors on the atomistic/continuum interface (which are often attributed to spurious forces called "ghost forces"). There are only few existing works that propose a coupling which is sufficiently accurate near the interface under certain limitations. In this paper a novel coupling that is free from "ghost forces" is proposed for a two-body interaction potential under the assumptions of either (i) one spatial dimension, or (ii) two spatial dimensions and piecewise affine finite elements for describing the continuum deformation. The performance of the proposed coupling is demonstrated with numerical experiments. The coupling strategy is based on judiciously defining the contributions of the atomistic bonds to the discrete and the continuum potential energy. The same method in one dimension has been independently developed and analyzed in Li and Luskin (arXiv:1007.2336).

preprint2010arXiv

The Spectrum of the Force-Based Quasicontinuum Operator for a Homogeneous Periodic Chain

We show under general conditions that the linearized force-based quasicontinuum (QCF) operator has a positive spectrum, which is identical to the spectrum of the quasinonlocal quasicontinuum (QNL) operator in the case of second-neighbour interactions. Moreover, we establish a bound on the condition number of a matrix of eigenvectors that is uniform in the number of atoms and the size of the atomistic region. These results establish the validity of and improve upon recent conjectures ([arXiv:0907.3861, Conjecture 2] and [arXiv:0910.2013, Conjecture 8]) which were based on numerical experiments. As immediate consequences of our results we obtain rigorous estimates for convergence rates of (preconditioned) GMRES algorithms, as well as a new stability estimate for the QCF method.