Next Article in Journal
A pH-Gated Functionalized Hollow Mesoporous Silica Delivery System for Photodynamic Sterilization in Staphylococcus aureus Biofilm
Previous Article in Journal
Tuning of Magnetic Damping in Y3Fe5O12/Metal Bilayers for Spin-Wave Conduit Termination
Previous Article in Special Issue
Scaling Scientific Cellular Automata Microstructure Evolution Model of Static Recrystallization toward Practical Industrial Calculations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimization of Monolayer MoS2 with Prescribed Mechanical Properties

by
Wacław Kuś
1,*,
Mohammed Javeed Akhter
2 and
Tadeusz Burczyński
2
1
Department of Computational Mechanics and Engineering, Silesian University of Technology, 44-100 Gliwice, Poland
2
Institute of Fundamental Technological Research, Polish Academy of Sciences, 02-106 Warsaw, Poland
*
Author to whom correspondence should be addressed.
Materials 2022, 15(8), 2812; https://doi.org/10.3390/ma15082812
Submission received: 10 February 2022 / Revised: 6 April 2022 / Accepted: 8 April 2022 / Published: 12 April 2022
(This article belongs to the Special Issue Computational Modelling and Design of Novel Engineering Materials)

Abstract

:
Various technological challenges are essentially material problems in our times. New functional and functional graded nanomaterials are constructed of components with predefined properties. The design of nanostructures with predefined mechanical properties was considered in this paper. This study applies the evolutionary algorithm (EA) to the optimization problem in the design of nanomaterials. The optimal design combined EA with molecular dynamics to identify the size of the void for the prescribed elastic properties in monolayer 2D MoS2 nanostructures. The numerical results show that the proposed EA and the use of optimization method allowed accurately obtaining nanostructures with predefined mechanical material properties by introducing elliptical voids in the 2D MoS2 nanosheets.

1. Introduction

In the ever-growing nanotechnology field, materials with unique properties have continuously been presented to meet the growing demand for practical industries and engineering. Therefore, it becomes vital to investigate the material properties of these nanomaterials to provide information to understand and then design nanostructures. Studies on the numerous properties of monolayer 2D molybdenum disulfide MoS2 nanomaterial have attracted many researchers in this field. Transition metal dichalcogenide (TMD) flat MoS2 is a triple layer of molybdenum and sulfur atoms arranged in a hexagonal crystal lattice. It has excellent mechanical [1] electrical and chemical properties [2,3]. Consequently, MoS2 has been the focus of substantial research in recent years, ensuring that the material can be used in a wide range of emerging technologies and future applications, including, for example, nanoporous MoS2 membranes for efficient reverse osmosis desalination and gas separation, membrane separation, desalination of water [4], DNA sequencing [5], and power generation [6,7]. Therefore, the properties and behavior of monolayer MoS2 must be understood and accurately predicted under various conditions to introduce MoS2 to novel applications. However, fabricated MoS2 sheets typically contain a variety of defects, including nanopores [8,9]. Thus, microstructural voids can profoundly impact the properties of MoS2, which helps design the nanomaterial, ultimately influencing the performance of MoS2-based devices. Likewise, during the growth and processing of MoS2, different topological defects (such as vacancies, inclusions, dislocation, and grain boundaries) and other sizable defects (such as nano-holes and nano-cracks) are inevitable, which can compromise the expected performance in the preparation and handling of MoS2-based nanodevices [10]. In addition, such defects can also occur due to the conditions in which MoS2-based devices are used [11]. Functional materials in nanosystems are often based on well-known materials but with tuned material properties. The 2D MoS2 sheet can be functionalized by introducing voids of appropriate sizes. The sizes of the voids determine the mechanical properties of the material. In recent years, researchers have carried out investigations on the mechanical performance of monolayer MoS2 related to aspects such as defects, inclusions, strength, damage, debonding, and failure that greatly influence its properties [1,12].
This work aims to study the identification problem in order to find the parameters that describe the shape of the void as a function of the properties of the prescribed material. This approach can be used in future applications where functional or functional graded materials or nanomaterials can be designed. Materials with a priori prescribed material properties can be created by the introduction of voids. The method presented in the paper is based on two components: a developed optimization algorithm and a direct problem solver. The nanostructures were simulated using the molecular dynamics (MD) method. Monolayer MoS2 with various mechanical properties can be computed by molecular statistics or MD methods with a set of numerical tests comprising uniaxial tension, compression, and shear [1,13,14,15]. The well-known MD code LAMMPS [16] was used in this study for direct problem solving. In recent years, inverse methods have been widely applied to predict the properties of structures and materials. These methods allow solving problems regarding parameters using optimization techniques and a set of direct problem solutions. Inverse methods were used for mechanical and thermomechanical problems in which the properties of materials and their structures were searched [17]. Inverse problems can be solved using direct problem formulations computed using numerical methods such as the finite element method (FEM), boundary element method (BEM), and MD. The objective function of optimization algorithms in inverse methods is, in most cases, multimodal; thus, global optimization techniques are often used during the problem-solving process. Applications of new approaches of evolutionary algorithms coupled with BEM computation in optimization and identification for cracked structures and internal void defects under thermomechanical and dynamical loading were shown in [17,18,19]. Sigmund [20] used the inverse homogenization method to tune the elastic properties of the material for periodic truss, frame, and continuum structures, as well as design microstructures with prescribed elastic properties and negative Poisson’s ratios. An in-house implementation of EA was used to search for new stable molecular graphene-like 2D materials [21,22,23].
In our earlier work [1], we computed the mechanical properties, i.e., independent elastic properties, for monolayer MoS2 with single and multiple random defects. The computational results from our previous work showed the significant influence of defects on the elastic material properties of MoS2 nanosheets. In this work, we aimed to identify the void size for the prescribed elastic properties defined by the user. We intuitively define the prescribed elastic properties, assuming that the sheet contains a void, which should be lower than the material without a void. Thus, we use the tools discussed earlier in the implementation of EA coupled with LAMMPS. As a result, we obtained the size of the void for the prescribed elastic properties.
The remainder of this work is arranged as follows: Section 2 presents the evolutionary optimization and the evaluation of the objective function. Section 3 details the molecular dynamics modeling of 2D MoS2 with voids and the evolutionary identification of voids with prescribed properties by minimizing the objective function. Numerical identification examples proving the ability of this method in solving the intended optimization problem iteratively are also provided.

2. Materials and Methods

2.1. Optimization Problem Formulation

The goal of the optimization problem is to design in an automatic way MoS2 structures with predefined material properties. The mechanical stiffness is taken into account in this paper; however, the optimization problem can also be solved for thermal, optical, or other properties of the microstructure. The objective function depends on the prescribed material properties and the actual properties computed for each design of the microstructure. The design vector ch may define the size, shape, and topology of the microstructure. The MoS2 structure considered in this paper was modified by introducing voids with the properties described using design variables   g i . The optimization problem can be formulated as
{ f i n d                             ch = ( g 1 , g 2 , g N ) m i n i m i z e               f ( ch ) = P P r e f s . t               g i L   g i   g i U ,
where g i L and g i U are the lower and upper constraints of design variables. Tensor or vector-containing tensor element P denotes the nanostructure properties obtained for design vector ch, and P r e f describes the prescribed reference properties of the nanostructure. The goal of minimizing the difference between current and reference material properties is set as a function of design variables leading to an objective function equal to zero (i.e., identical reference and obtained properties). Optimization with a small difference between the reference and obtained properties is also acceptable. The paper is devoted to the optimization of the nanostructure taking into account mechanical properties (in our case, P and P r e f ) depending on the stiffness of the nanostructure with introduced voids. The stress–strain relationship for small strains can be expressed with Voigh notation as follows:
σ i j = P i j ε i j ,
where ε i j denotes the strain components, and P i j denotes the elastic constants to be used during objective function evaluation. The shape of the nanostructure is modified according to ch by introducing an elliptic void, as shown in Figure 1.
The 2D infinite nanostructure was modeled using periodic boundary conditions, and the maximum size of the void was defined by the unit cell size used in simulations. The mechanical properties (stiffness in two directions) were computed using the MD approach. The algorithm for determining the stiffness of the nanostructure is shown in Figure 2.
First, the void was introduced into the pristine 2D nanostructure. Next, relaxation was performed, and the stress–strain curve was obtained in two directions during uniaxial tensile load. In general, other loads can be applied to obtain, for example, shear stresses or thermal and optical properties of the nanostructure.
This paper contains examples of optimization based on the computation of an objective function using MD results for stiffness computed in two directions for monolayer MoS2. The atomic model consisted of about 10,000 atoms in the nanosheet for a domain size of 175 Å2. The Stillinger–Weber (SW) [24,25] potential was applied to describe interatomic interactions between atoms. Before tensile deformation, the model was relaxed at 300 K and 0 bar pressure through an isothermal–isobaric ensemble (NPT) for 30 ps (picoseconds). All MoS2 tensile deformations were carried out at a constant temperature of 300 K. The Nose–Hoover thermostat was used to maintain the temperature of the simulation system. The position and velocity of all atoms were updated by the Verlet integration algorithm. The LAMMPS software package and the Open Visualization Tool (OVITO) [26] were used for molecular dynamics simulation, visualization, and output data analysis. Uniaxial tensile deformation at a constant strain rate of 0.0001 ps−1 was applied to estimate the stiffness of the structure. The stress tensor components [27] were calculated as follows:
σ a b = 1 V [ 1 2 i N * j ( 1 ) N f i j a r i j b + m i u i a u i b ] ,
where a and b denote the Cartesian components, fij is the force acting on atom i due to another atom j, V, mi, and ui are the volume, mass, and velocity of atom i, and N is the number of atoms. Figure 3 shows examples of two simulation results with a stiffness of 163 GPa and 158 GPa.
The objective function value could be computed for such a case as the sum of the absolute difference between computed stiffness and prescribed material properties.

2.2. Evolutionary Optimization

Optimization can be defined as selecting the best choice (according to some criteria) from a set of available variants. There are mainly two groups of optimization methods. The first group includes local optimization methods, mostly leading to the local optimum, where the results depend on the starting vector values. The second group includes global optimization techniques, mostly based on heuristics. Heuristic methods are best suited for multimodal functions, where most algorithms fail to determine the global optimum. Many optimization techniques are based on biologically inspired methods such as natural selection, learning procedures, and probabilistic rules. Evolutionary algorithms [28] and particle swarm optimization (PSO) [29] are examples of these techniques. These techniques can solve multimodal optimization problems in mathematics and engineering. It is impossible to identify the best optimization techniques among those mentioned above, since each algorithm has its strengths and weaknesses, and the performance depends on the optimization problems, constraints, and algorithms parameters. The evolutionary algorithm (EA) searches the space of possible solutions on the basis of mechanisms taken from the evolution of species.
The flowchart of our implementation of EA used in this work is shown in Figure 4. The algorithm starts by generating an initial population of individuals generated in a random or ordered way according to the requirements of the problem. The individual, containing chromosomes, represents a single solution. Usually, in applications of EA, individuals contain only one chromosome ch, with vectors of genes representing design variables gi. Genes may contain coded design variables; however, in our approach, we used floating point genes. Hence, additional coding was not needed.
For each individual, MD simulations and individual fitness functions were calculated in the next step of the algorithm on the basis of the atomic structure defined by genes. The selection procedure was performed to choose individuals for the next iteration considering their fitness values. The probability of survival of an individual depends on the value of the fitness function. An individual with good fitness has a better chance of survival during the selection process. The ranking selection is performed in a few steps. First, individuals are classified according to the values of the fitness function; then, a rank value is assigned to each individual. This depends on the individual’s number and the rank function. The best individuals obtain the highest rank value; the worst individuals obtain the lowest ones. In the final step, individuals for the offspring generation are drawn, but the probability of drawing particular individuals is closely related to their rank value. The process is repeated iteratively until the termination condition is satisfied. The termination condition may be formulated as a maximum number of iterations. In cases when the stop condition is not fulfilled, the genes are modified using evolutionary operators. The algorithm in this paper used evolutionary operators such as uniform and Gaussian mutation, as well as simple and arithmetic crossover. The uniform and Gaussian mutations modify individuals randomly. Simple and arithmetic crossovers create new chromosomes on the basis of two randomly chosen chromosomes from the population. The modified individuals are introduced into the population.
The objective function (i.e., the fitness function for each chromosome) is computed on the basis of a direct problem solution for each individual. The number of direct problems solved in EA is quite large due to the number of chromosomes and the number of iterations of the algorithms. To minimize the total optimization time (wall time), a parallel approach is crucial, and it was used to compute the numerical examples presented in this paper.

3. Results

Nanostructure optimization by introducing a void according to prescribed mechanical properties is illustrated using a few cases in the section. The tests were performed for a few sets of prescribed material properties. The evolutionary algorithm featured two subpopulations, a total number of individuals of 32, tournament selection with a tournament size of 5, uniform selection with a probability of 0.3, Gaussian mutation with a probability of 0.5, simple crossover with a probability of 0.1, and arithmetic crossover with a probability of 0.1. The parameters of the EA were chosen on the basis of previous experience gained during the solving of structural optimization problems [19,21,22,23]. The number of EA iterations was set as 50.
An optimized nanostructure with an approximate size of 170 Å × 170 Å containing an elliptical void was used. The size of the void could range from 1 × 1 up to 50 × 50 Å (radius of the elliptic void). The stop condition was formulated as the maximum number of iterations.
The evaluation of the fitness function for each individual was performed as described in Section 2. The LAMMPS software was used to solve two problems to obtain the stiffness in two directions, and the objective function was calculated on the basis of these results. The computing of the fitness functions in each iteration can be performed in parallel way with very good efficiency when the number of processing units is equal to the number of MD problems (number of individuals times number of mechanical properties; in our case, 64). Additional processing units may also be used for parallelizing each MD simulation. Such an approach can be used to perform computations with a high number of processing units. The results presented in the section were obtained with the use of supercomputers Okeanos and Karolina. Hundreds to thousands of processor cores were used during computations.
The changes in the best objective function over few iterations for adequate nanostructures with prescribed material properties Pref11 = Pref22 = 160 GPa are shown in Figure 5.
A total of three distinct void identification analyses were performed for the MoS2 nanosheet with different prescribed material properties. The void (elliptical void) was induced at the center of the sheet. Table 1 contains the values of the best obtained solutions for the numerical tests. The obtained material properties (P11, P22) for the ellipse radius (g1, g2) were not identical to those prescribed (Pref), as denoted by the errors eP11 and eP22. The obtained properties were very close to the prescribed ones.
The resulting structures for the above cases are shown in Figure 6. The nanostructure in the first case should have an ellipse with a bigger radius in the y-direction; the radii in in the second case should be similar, whereas those in the third case should be similar to the first case when considering the prescribed material properties. The MoS2 nanostructure was not symmetric in the x (g1) and y (g2) directions; thus, we did not expect the same ellipse radii for the second case in the x- and y-directions. The results agreed with the intuitive approach and gave exact values of void size. The method can be used for any prescribed stiffness values; however, it is of course limited by the size of the void and maximum stiffness of the MoS2 sheet.

4. Conclusions

The appearance of defects in MoS2 can weaken its mechanical properties, such as fracture strength and Young’s modulus. However, such defects have potential in novel applications of MoS2, such as graded materials or nanosystems. We used the optimization method to tune the material properties for a periodic monolayer MoS2 with a void as a defect. The tailoring problem was formulated as an optimization problem to identify the simplest possible nanostructure with constraints given by the prescribed elastic constants. Numerical examples proved that the proposed method could find the void size with prescribed mechanical properties. The present study shows the potential of molecular simulations for 2D nanostructures. It is revealed as an efficient method to design nanostructures with prescribed properties. The method is generally applicable and can be used to take into account the thermal or optical properties of the nanostructure by modifying the components of the objective function and using an adequate direct problem method.

Author Contributions

Conceptualization, W.K. and T.B.; methodology, W.K., M.J.A. and T.B.; software, W.K.; writing—original draft preparation, W.K. and M.J.A.; writing—review and editing, W.K., M.J.A. and T.B. All authors read and agreed to the published version of the manuscript.

Funding

The research was funded by the National Science Center, Poland, project no. 2016/21/B/ST8/02450.

Acknowledgments

Optimizations were performed in part at the Interdisciplinary Center for Mathematical and Computational Modeling, University of Warsaw, under grant GB80-16. We acknowledge EuroHPC JU for awarding us access to supercomputer Karolina hosted by IT4Innovations, National Supercomputing Center, Ostrava, Czechia, under project 2010PA6018.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Akhter, M.J.; Kuś, W.; Mrozek, A.; Burczyński, T. Mechanical Properties of Monolayer MoS2 with Randomly Distributed Defects. Materials 2020, 13, 1307. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Khan, M.F.; Nazir, G.; Lermolenko, V.M.; Eom, J. Electrical and photo-electrical properties of MoS2 nanosheets with and without an Al2O3 capping layer under various environmental conditions. Sci. Technol. Adv. Mater. 2016, 17, 166–176. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Li, X.; Zhu, H. Two-dimensional MoS2: Properties, preparation, and applications. J. Mater. 2015, 1, 33–44. [Google Scholar] [CrossRef] [Green Version]
  4. Heiranian, M.; Farimani, A.B.; Aluru, N.R. Water desalination with a single-layer MoS2 nanopore. Nat. Commun. 2015, 6, 8616. [Google Scholar] [CrossRef] [Green Version]
  5. Feng, J.; Liu, K.; Bulushev, R.D.; Khlybov, S.; Dumcenco, D.; Kis, A.; Radenovic, A. Identification of single nucleotides in MoS2 nanopores. Nat. Nanotechnol. 2015, 10, 1070–1076. [Google Scholar] [CrossRef]
  6. Yang, T.; Lin, H.; Zheng, X.; Loh, K.P.; Jia, B. Tailoring pores in graphene-based materials: From generations to applications. J. Mater. Chem. A 2017, 5, 16537–16558. [Google Scholar] [CrossRef] [Green Version]
  7. Macha, M.; Marion, S.; Nandigana, V.V.; Radenovic, A. 2D materials as an emerging platform for nanopore-based power generation. Nat. Rev. Mater. 2019, 4, 588–605. [Google Scholar] [CrossRef]
  8. Su, S.; Wang, X.; Xue, J. Nanopores in two-dimensional materials: Accurate fabrication. Mater. Horiz. 2021, 8, 1390–1408. [Google Scholar] [CrossRef]
  9. Graf, M.; Lihter, M.; Thakur, M.; Georgiou, V.; Topolancik, J.; Ilic, B.R.; Liu, K.; Feng, J.; Astier, Y.; Radenovic, A. Fabrication and practical applications of molybdenum disulfide nanopores. Nat. Protoc. 2019, 14, 1130–1168. [Google Scholar] [CrossRef]
  10. Hong, J.; Hu, Z.; Probert, M.; Li, K.; Lv, D.; Yang, X.; Gu, L.; Mao, N.; Feng, Q.; Xie, L.; et al. Exploring atomic defects in molybdenum disulphide monolayers. Nat. Commun. 2015, 6, 6293. [Google Scholar] [CrossRef] [Green Version]
  11. Vancsó, P.; Magda, G.Z.; Pető, J.; Noh, J.; Kim, Y.; Hwang, C.; Biró, L.P.; Tapasztó, L. The intrinsic defect structure of exfoliated MoS2 single layers revealed by Scanning Tunneling Microscopy. Sci. Rep. 2016, 6, 29726. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Wang, S.; Robertson, A.; Warner, J.H. Atomic structure of defects and dopants in 2D layered transition metal dichalcogenides. Chem. Soc. Rev. 2018, 47, 6764. [Google Scholar] [CrossRef] [PubMed]
  13. Mahata, A.; Jiang, J.; Mahapatra, D.R.; Rabczuk, T. Effect of intrinsic structural defects on mechanical properties of single layer MoS2. Nano-Struct. Nano-Objects 2019, 18, 100247. [Google Scholar] [CrossRef]
  14. Li, Y.; Chen, P.; Zhang, C.; Peng, J.; Gao, F.; Liu, H. Molecular dynamics simulation on the buckling of single-layer MoS2 sheet with defects under uniaxial compression. Comput. Mater. Sci. 2019, 162, 116–123. [Google Scholar] [CrossRef]
  15. Li, M.; Wan, Y.; Tu, L.; Yang, Y.; Lou, J. The Effect of VMoS3 Point Defect on the Elastic Properties of Monolayer MoS2 with REBO Potentials. Nanoscale Res. Lett. 2016, 11, 155. [Google Scholar] [CrossRef] [Green Version]
  16. Thompson, A.P.; Aktulga, H.M.; Berger, R.; Bolintineanu, D.S.; Brown, W.M.; Crozier, P.S.; in’t Veld, P.J.; Kohlmeyer, A.; Moore, S.G.; Nguyen, T.D.; et al. LAMMPS-a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. Comp. Phys. Comm. 2022, 271, 10817. [Google Scholar] [CrossRef]
  17. Burczyński, T.; Kuś, W.; Beluch, W.; Długosz, A.; Poteralski, A.; Szczepanik, M. Intelligent Computing in Optimal Design. In Solid Mechanics and Its Applications; Springer: Berlin/Heidelberg, Germany, 2020. [Google Scholar]
  18. Kokot, G.; Burczyński, T. Evolutionary Algorithms and Boundary Element Method in Generalized Shape Optimization. J. Theor. Appl. Mech. 2003, 41, 341–364. [Google Scholar]
  19. Kokot, G.; Orantek, P. The Topology Optimization Using Evolutionary Algorithms. In IUTAM Symposium on Evolutionary Methods in Mechanics. Solid Mechanics and Its Applications 117; Burczyński, T., Osyczka, A., Eds.; Springer: Berlin/Heidelberg, Germany, 2004; pp. 173–186. [Google Scholar] [CrossRef]
  20. Sigmund, O. Tailoring materials with prescribed elastic properties. Mech. Mater. 1995, 20, 351–368. [Google Scholar] [CrossRef]
  21. Mrozek, A.; Kuś, W.; Burczyński, T. Searching of stable configurations of nanostructures using computational intelligence methods. Tech. Sci. 2010, 20, 85–97. [Google Scholar]
  22. Mrozek, A.; Kuś, W.; Burczyński, T. Nano level optimization of graphene allotropes by means of a hybrid parallel evolutionary algorithm. Comput. Mater. Sci. 2015, 106, 161–169. [Google Scholar] [CrossRef]
  23. Mrozek, A.; Kuś, W.; Burczyński, T. Method For Determining Structures of New Carbon-based 2D Materials with Predefined Mechanical Properties. Int. J. Multiscale Comput. Eng. 2017, 15, 379–394. [Google Scholar] [CrossRef]
  24. Kandemir, A.; Yapicioglu, H.; Kinaci, A.; CaÇin, T.; Sevik, C. Thermal transport properties of MoS2 and MoSe2 monolayers. Nanotechnology 2016, 27, 55703. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Wen, M.; Shirodkar, S.N.; Plech, P.; Kaxiras, E.; Elliott, R.S.; Tadmor, E.B. A force-matching Stillinger-Weber potential for MoS2: Parameterization and Fisher information theory based sensitivity analysis. J. Appl. Phys. 2017, 122, 244301. [Google Scholar] [CrossRef]
  26. Stukowski, A. Visualization and analysis of atomistic simulation data with OVITO—the Open Visualization Tool. Simul. Mater. Sci. Eng. 2010, 18, 015012. [Google Scholar] [CrossRef]
  27. Thompson, A.P.; Plimpton, S.J.; Mattson, W. General formulation of pressure and stress tensor for arbitrary many-body interaction potentials under periodic boundary conditions. J. Chem. Phys. 2009, 131, 154107. [Google Scholar] [CrossRef] [Green Version]
  28. Michalewicz, Z. Genetic Algorithms + Data Structures = Evolutionary Algorithms; Springer: Berlin, Germany, 1996. [Google Scholar]
  29. Kennedy, J.; Eberhart, R.C.; Shi, Y. Swarm Intelligence; Morgan Kaufmann Publishers: San Francisco, CA, USA, 2001. [Google Scholar]
Figure 1. The nanostructure with an elliptical void described by two design variables.
Figure 1. The nanostructure with an elliptical void described by two design variables.
Materials 15 02812 g001
Figure 2. Determination of nanostructure stiffness in two directions. An ellipse is introduced into the pristine nanostructure by removing atoms in black area. Next, the structure is relaxed, and two analyses of microstructure stretching are performed in two different directions. Then, the stiffness is computed on the basis of the MD results.
Figure 2. Determination of nanostructure stiffness in two directions. An ellipse is introduced into the pristine nanostructure by removing atoms in black area. Next, the structure is relaxed, and two analyses of microstructure stretching are performed in two different directions. Then, the stiffness is computed on the basis of the MD results.
Materials 15 02812 g002
Figure 3. Examples of stress-strain relationship during uniaxial tension in two directions.
Figure 3. Examples of stress-strain relationship during uniaxial tension in two directions.
Materials 15 02812 g003
Figure 4. Evolutionary algorithm flowchart.
Figure 4. Evolutionary algorithm flowchart.
Materials 15 02812 g004
Figure 5. Progress of the convergence of the objective function indicating the evolution of the void during iterative generations for prescribed elastic properties Pref11 = 160 GPa, Pref22 = 160 GPa. The corresponding void dimensions obtained during generations were (i) g1 = 46.06 Å and g2 = 2.95 Å, (ii) g1 = 32.18 Å and g2 = 6.07 Å, (iii) g1 = 33.92 Å and g2 = 8.86 Å, (iv) g1 = 21.36 Å and g2 = 14.88 Å, (v) g1 = 18.41 Å and g2 = 19.38 Å, and (vi) g1 = 28.12 Å and g2 = 23.15 Å.
Figure 5. Progress of the convergence of the objective function indicating the evolution of the void during iterative generations for prescribed elastic properties Pref11 = 160 GPa, Pref22 = 160 GPa. The corresponding void dimensions obtained during generations were (i) g1 = 46.06 Å and g2 = 2.95 Å, (ii) g1 = 32.18 Å and g2 = 6.07 Å, (iii) g1 = 33.92 Å and g2 = 8.86 Å, (iv) g1 = 21.36 Å and g2 = 14.88 Å, (v) g1 = 18.41 Å and g2 = 19.38 Å, and (vi) g1 = 28.12 Å and g2 = 23.15 Å.
Materials 15 02812 g005
Figure 6. Monolayer MoS2 nanosheet with void identified by optimization: (a) Pref11 = 180, Pref22 = 150 GPa (g1 = 15.35 Å and g2 = 33.44 Å); (b) Pref11 = 160, Pref22 = 160 GPa (g1 = 28.12 Å and g2 = 23.15 Å); (c) Pref11 = 150, Pref22 = 180 GPa (g1 = 36.35 Å and g2 = 12.19 Å).
Figure 6. Monolayer MoS2 nanosheet with void identified by optimization: (a) Pref11 = 180, Pref22 = 150 GPa (g1 = 15.35 Å and g2 = 33.44 Å); (b) Pref11 = 160, Pref22 = 160 GPa (g1 = 28.12 Å and g2 = 23.15 Å); (c) Pref11 = 150, Pref22 = 180 GPa (g1 = 36.35 Å and g2 = 12.19 Å).
Materials 15 02812 g006
Table 1. The prescribed and resulting stiffness, ellipse radius, and error of obtained stiffness.
Table 1. The prescribed and resulting stiffness, ellipse radius, and error of obtained stiffness.
CasePref11
(GPa)
P11
(GPa)
Pref22
(GPa)
P22
(GPa)
g1
(Å)
g2
(Å)
eP11 (%)eP22 (%)
1150.0149.2180.0179.536.3512.190.50.3
2160.0162.6160.0158.028.1223.151.61.3
3180.0179.5150.0148.015.3533.44 0.31.3
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kuś, W.; Akhter, M.J.; Burczyński, T. Optimization of Monolayer MoS2 with Prescribed Mechanical Properties. Materials 2022, 15, 2812. https://doi.org/10.3390/ma15082812

AMA Style

Kuś W, Akhter MJ, Burczyński T. Optimization of Monolayer MoS2 with Prescribed Mechanical Properties. Materials. 2022; 15(8):2812. https://doi.org/10.3390/ma15082812

Chicago/Turabian Style

Kuś, Wacław, Mohammed Javeed Akhter, and Tadeusz Burczyński. 2022. "Optimization of Monolayer MoS2 with Prescribed Mechanical Properties" Materials 15, no. 8: 2812. https://doi.org/10.3390/ma15082812

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop