1. Introduction
Two-dimensional (2D) crystals are a relatively new class of structure that includes several families of various materials. The most popular and well-known among 2D crystals is graphene—an atomically thin layer of carbon atoms with hexagonal structure and covalent bonds, discovered in 2004 [
1]. Other popular representatives of 2D crystals are 2D boron nitride [
2], 2D molybdenum disulfide MoS
2 [
3], 2D transition metal dichalcogenides (2D TMDs), and other materials [
4]. Since their discovery, two-dimensional materials have opened new possibilities to significantly improve nanoelectronics, due to their ability to be integrated into devices with atomic-scale precision, leading to the development of smaller, faster, and more efficient electronic components. Thus, each of the 2D materials has unique physical properties that define their potential application. For instance, graphene has exceptional electron mobility, mechanical strength, and thermal conductivity. Its application in nanoelectronics includes the development of high-speed transistors, flexible transparent conductive films, and advanced sensors [
5,
6]. Graphene-based transistors, for instance, can operate at higher frequencies compared to traditional silicon-based transistors, enabling faster data processing and communication [
7,
8].
In 2011, the family of 2D materials became even bigger, as a new group of crystals—transition metals carbides and nitrides (as well as carbonitrides) with general chemical formula M
n+1X
n (M is an early transition metal, and X is carbon C and/or nitrogen N) called MXenes—were discovered [
9]. MXenes are synthesized from bulk precursors, so-called MAX-phases (nanolaminated materials with general chemical formula M
n+1AX
n, where M, as above, is an early transition metal, A is an element from group IIIA or IVA (or group 13 or 14), and X is carbon C and/or nitrogen N, and n = 1, 2, or 3) by etching a layer of element A atoms from MAX phases [
10,
11].
Together with unique electrochemical properties [
12,
13,
14], MXenes exhibit exceptional mechanical properties [
15,
16], including high strength, flexibility, and tunability, making them highly attractive for applications in nanoelectronics. Their ability to withstand mechanical deformation while maintaining electronic functionality opens up new possibilities for the development of flexible, durable, and high-performance electronic devices [
17,
18].
Additionally, the layered structure of MXenes contributes to their mechanical stability and flexibility. The layers can slide relative to each other, providing an intrinsic mechanism for accommodating strain and preventing fracture. This characteristic is particularly beneficial for creating multilayered or composite structures that combine MXenes with other materials, resulting in devices that are both mechanically robust and electrically superior. The integration of MXenes into nanoelectromechanical systems (NEMSs), leverages their mechanical properties to enhance device performance and longevity [
17,
18].
The flexibility and mechanical robustness of MXenes are further enhanced by their unique surface chemistry. MXenes can be functionalized with various surface groups, such as hydroxyl, oxygen, and fluorine, which not only modify their electronic properties but also influence their mechanical behavior. This tunability allows researchers to design MXenes with tailored mechanical properties to meet specific application requirements. For instance, by adjusting the surface chemistry, MXenes can achieve an optimal balance between flexibility and strength, making them ideal for applications in nanoelectronics where mechanical reliability is paramount.
One of the important properties of nanomaterials that defines their use in nanoelectromechanical system devices is bending or flexural rigidity [
19,
20]. It characterizes the ability of a particular sample to resist bending deformation and maintain its flexibility and elasticity. However, the experimental measurement of the mechanical properties of nanomaterials requires special facilities and represents a serious challenge due to the small sizes of the samples (see for example [
21,
22,
23]). Therefore, mechanical properties of nanomaterials are also studied by computational [
15,
16,
17,
18,
19,
20] and theoretical [
24,
25,
26] methods. Thus, one of the illustrative examples of such study is the calculation of bending rigidity of graphene nanoribbons (GNRs) by classical molecular dynamics simulations [
20]. The developed approach also allowed the authors to propose a model for an ultrasensitive pressure sensor based on GNRs, and also suggested further applications of graphene nanoribbons in NEMSs devices in general [
27].
A similar approach was recently adopted to study the bending properties of the Ti
n+1C
n MXene nanoribbons [
19]. The performed calculations confirmed the assumption that Ti
n+1C
n samples are characterized by higher bending rigidity than graphene, due to the higher thickness of MXene nanoribbons, as well as the increasing of bending rigidity with the growth of the thickness of the samples. Moreover, recent discoveries of the ordered double transition metal carbides [
28,
29] opened a new way of controlling and tuning the electronic properties of MXenes by altering the outer metal layer in M
n+1X
n sheets. Thus, after the synthesis of double transition metals, MXene Mo
2Ti
2C
3 was reported [
29], which naturally opened a question of how the changing of two outer layers of titanium in Ti
4C
3 to molybdenum affects the flexural properties of MXene nanoribbons. In this paper, we report the comparative analysis of the bending properties of Mo
2Ti
2C
3 and Ti
4C
3 by classical molecular dynamics simulations that were performed using a previously developed computational scheme of numerical experiments on bending of Ti
n+1C
n MXene nanoribbons [
19]. As no experiments on bending rigidity of pristine Mo
2Ti
2C
3 and Ti
4C
3 nanoribbons have been reported by now, our study may provide the first insights on how mechanical properties of MXenes can be affected by replacing monometallic Ti sample by (Mo, Ti) double transition metal carbide.
2. Model
We followed a previous study on the bending rigidity of Ti
n+1C
n MXene nanoribbons [
19], where preliminary results on the bending of Ti
4C
3 nanoribbon were reported. In our simulations, we consider two different MXene nanoribbons, Mo
2Ti
2C
3 and Ti
4C
3. The chemical composition of the considered samples was chosen due to the already developed MD approach for the simulation of Ti
n+1C
n samples [
15] and similar structures of Mo
2Ti
2C
3 and Ti
4C
3 MXenes [
28,
29]. As it is reported in [
29], the double transition metal Mo
2Ti
2C
3 is characterized by almost equivalent interatomic distances between Mo-C and Ti-C atoms, similar lattice parameters, and the same sample thickness as Ti
4C
3 MXene. The structures of both Mo
2Ti
2C
3 and Ti
4C
3 MXenes and the general view of the studied nanoribbons are shown in
Figure 1 and
Figure 2, respectively (all snapshots of atomistic configurations of the samples were prepared with visual molecular dynamics software version 1.9 [
30]).
Moreover, with the assumption that bending rigidity depends on the thickness
h of the sample being
[
19], simulations of the thickest (seven atomic layers) M
4X
3 samples under bending deformation may provide the approximal value of maximal bending rigidity that can be expected among (Mo,Ti)
n+1C
n and Ti
n+1C
n MXenes.
It is worth to note that several approaches for MD simulations of MXenes were reported in the literature so far (see for example [
15,
16,
31,
32]). The most precise and yet most complicated approach employs ReaxFF interatomic potential [
33]. Such an approach allowed us to obtain important data on chemical reactions and phase transitions within MXenes, and to study their mechanical properties as well [
16,
33]. However, ReaxFF can be straightforwardly used for simulations of certain MXenes only, and its application for Mo
2Ti
2C
3 requires additional parametrization, which is a rather hard challenge [
33]. Therefore, in our simulations, we adopted a previously developed approach that is based on combinations of interatomic potentials [
15]. The full description of the developed model is provided in the original paper; therefore, here we focused on the adjustments that were introduced in the model to adapt for simulations of Mo
2Ti
2C
3 MXene. Following the assumption that interactions between titanium atoms in Ti
n+1C
n MXenes can be described within embedded atom method (EAM) [
15], we suppose that EAM potential can also be used to describe interactions between molybdenum atoms in Mo
2Ti
2C
3 samples in a similar manner, as EAM potential is already parametrized for molybdenum [
34]. Moreover, with this assumption, interactions between Ti and Mo atoms in Mo
2Ti
2C
3 MXene sample can be described with an EAM alloy model for Ti-Mo [
34]. As in Mo
2Ti
2C
3, Mo-C and Ti-C bonds have almost similar lengths [
29], interatomic forces between both metals and carbon atoms are calculated in the same way as forces between titanium and carbon in Ti
4C
3 [
19]. This is a rough approximation, which makes our model rather qualitative; nevertheless, it makes it possible to obtain preliminary data on flexural properties of double metals Mo
2Ti
2C
3 MXene, without the complicated parametrization of other interatomic potentials.
In our experiments on bending, we followed an already existing methodic that was proposed for studying the flexural properties of graphene nanoribbons [
20], and which was further adopted for the calculation of the bending rigidity of Ti
n+1C
n MXenes [
19]. We consider Mo
2Ti
2C
3 and Ti
4C
3 nanoribbons with sizes of 1.7 nm × 12.0 nm located in a Cartesian coordinate box with periodic boundary conditions in X and Y direction, and free boundary condition in
Z, as shown by schematics in
Figure 3. For the same reason, as in [
19], we applied similar constraints to the boundary atoms across long edges of the nanoribbons, not allowing them to relax in the XY plane during bending.
Thus, by applying the external force
F, as is shown in
Figure 3, and measuring the corresponding central deflection of nanoribbon
, it is possible to calculate the bending rigidity
D of the sample, from the well-known equation [
19,
20]:
where
L is the length of the nanoribbon and
β is the coefficient in the range from 0.0056 to 0.00725 for a nanoribbon length-to-width ratio from 1 to
, respectively.
In our simulation, we adopted an in-house code, developed for the study of mechanical properties and thermal stability of Ti
n+1C
n MXenes [
15,
19,
35,
36]. The code was implemented on a graphics processing unit (GPU) for parallel calculations. The calculations were performed on an NVIDIA Tesla P100 GPU (NVIDIA Corporation, Santa Clara, CA, USA). Other details of the simulation setup can be found in [
19]. The obtained results are described in the following section.
3. Results
As described above, the simulations of the bending deformation of MXene nanoribbons were performed by applying external force in the –Z direction to the central part of the sample (see
Figure 3). Applications of the external load results in a deflection of the nanoribbon from the initial configuration and its bending. Similar to the simulations described in [
19,
20], we performed the simulation by gradually increasing the force applied to the sample and analyzing its behavior. During the experiments, the time dependencies of the immediate position of the nanoribbon center and other data needed for calculations of bending parameters were recorded. The typical behavior of the samples under an external bending load is shown in
Figure 4.
Presented in
Figure 4, data show that at the beginning of experiments, the application of the external force of the smallest magnitude
F0 (equivalent to 0.72 nN) initiates damping oscillations of both nanoribbons (see opening part of the dependencies in
Figure 4a), after which the sample deflects from the original position. After applying forces of larger magnitudes (with an increment
), central deflection
was measured after slight damping oscillations when samples reached the stationary configuration. In our experiments, samples were allowed to reach a stationary position for 2 × 10
6 time steps (4 × 10
6 time steps for initial indentation, as shown in
Figure 4a). The central deflection was measured by time averaging of the
Z coordinate of the central part of the sample, to which the external force is applied (see
Figure 3) over the last 5 × 10
5 time steps during each force increment.
Figure 4b shows that each force increment
resulted in approximately constant changes in each sample’s increment of central deflection
, which suggests that in this part of experiment, nanoribbons are bending in an almost linear regime.
To illustrate the behavior of the samples in elastic regimes of bending deformation, we performed additional indentation of the Ti
4C
3 sample with the doubled force increment
, as well as the gradual relaxation of Mo
2Ti
2C
3 nanoribbon. Obtained data are shown in
Figure 5.
As it follows from the data presented in
Figure 5a, the application of a two-times larger force increment
to Ti
4C
3 sample results in approximately two times larger increments of related central deflection
compared to the initial conditions of the experiments. The presented dependencies are also characterized by a larger amplitude of damping oscillations after each increase in the bending force.
Figure 5b shows the behavior of the Mo
2Ti
2C
3 sample during the first 10 steps of indentation (increasing of the applied force
with constant increment
) and following relaxation (decreasing of
with constant increment
) to the previous position. As it follows in the figure, consecutive relaxation of the applied force results in decreases in the central deflection by almost the same increment
as during the indentation. Therefore, the presented dependence illustrates the reversibility of the performed experiments within the elastic mode of the bending deformation of nanoribbons.
Further in our experiments, the application of the increasing external force continued until the complete fracture of the nanoribbon. The total time dependencies of the
Z coordinate of nanoribbon centers and the load curves
for both samples are shown in
Figure 6a and
Figure 6b, respectively.
The fracture of nanoribbon is indicated by the sharp drop in
Z(
t) dependencies to large negative values of
Z (see
Figure 6a). As it follows from the figure, both samples are characterized by a close magnitude of critical deflection
nm before the destruction of the nanoribbon. At the same time, as
Figure 6b shows, the fracturing of Mo
2Ti
2C
3 nanoribbon was observed at an almost two times larger magnitude of the critical external force
(note that there were significantly more indentation steps and correspondingly more time steps in simulations for Mo
2Ti
2C
3 nanoribbon before fracturing, as shown in
Figure 6a).
It is important to note that nonlinear modes of bending, related to plastic deformation and partial destruction of the samples, are also clearly visible from the data presented in
Figure 6. The mentioned modes of bending relate to the part of the dependencies where the growth of the external force results in increases in the corresponding increments of central deflection
or in decreases in the current magnitude of central deflection
presented in
Figure 6b. The latter case relates to the plastic deformation and partial destruction of the sample, which may lead to the local rearrangement of atoms and relaxation of mechanical stresses. Moreover, it is worth noting that the nonlinear dependence of central deflection on the applied load
may also be caused by the axial extension effect [
37]. An example of atomistic configuration of Mo
2Ti
2C
3 nanoribbon with clearly visible plastic deformation and partial destruction of the central part is shown in
Figure 7.
Figure 7 shows a typical atomistic configuration of studied samples in different modes of bending deformation. The top panels show configurations of the samples at central deflection
nm, which relate to the elastic mode of bending, corresponding to linear parts of the dependencies of central deflection on external force
, as shown in
Figure 6b. The middle panels of
Figure 7 depict the configuration of the samples at larger central deflection
nm. At these conditions, force increment
no longer leads to constant increments of central deflection
; however, no plastic deformation or fracture of the samples is observed in atomistic configurations, and samples can be relaxed to initial configuration. The bottom panels of
Figure 7 show the fracture of the samples. Notably, samples are characterized by different fracture dynamics. Thus, Mo
2Ti
2C
3 nanoribbon is characterized by severe plastic deformation in the central area where the external force is applied, following the destruction of the MXene crystal structure and the fracture of the nanoribbon. At the same time, fracture dynamics of Ti
4C
3 nanoribbon are more typical for brittle materials, with the destruction of the sample in several areas, and following the formation of fragments that preserve the two-dimensional structure of Ti
4C
3 MXene. Such behavior of Ti
4C
3 MXenes was already observed in our earlier MD simulations [
19]. The bending and fracture dynamics of both nanoribbons are also shown in
Supplementary Videos S1 and S2.
Performing simulations and recorded data allowed us to calculate elastic parameters that characterize the bending properties of MXene nanoribbons. However, the straightforward installation of obtained
dependence (shown in
Figure 6b) in Equation (1) leads to the bending rigidity
D of the samples, as presented in
Figure 8.
As it can be seen from the figure, the calculated
D strongly depends on applied force even in the region of small magnitudes of central deflections. This makes it hard to obtain the exact value of the elastic parameters that characterize the bending properties of MXene nanoribbons. However, this feature can be resolved by analyzing the dependencies shown in
Figure 4a and
Figure 6.
As can be seen in
Figure 4a, the first increment of a central deflection
in the opening part of dependencies has a larger magnitude than the following deflection increments for both samples. Moreover, the system requires more simulation time to reach the stationary configuration, and related damped oscillations of the nanoribbon last longer, compared to the next steps of indentation. This may be explained by the presence of initial deflection
of the samples and related initial stresses (pre-tension) due to the applied boundary conditions in the
XY plane. It is worth to mention that the pre-tension of the samples is often observed during the calculation of mechanical parameters in MD simulations [
16]. The pre-tension of the samples can be taken into account by introducing the correction
to measure deflection
for each sample correspondingly. The magnitude of
can be estimated as a residual part of the first increment
after subtracting the deflection increment averaged over the next few steps of indentation,
=
. Thus, magnitudes of initial deflection, calculated according to this assumption, equal
nm and 0.103 nm for Mo
2Ti
2C
3, and Ti
4C
3 nanoribbons, respectively. Dependencies of bending rigidity
on external force, calculated from Equation (1) after introducing of such correction, are shown in
Figure 9.
As it follows from the figure, dependence calculated by taking into account of the initial deflection of Mo2Ti2C3 sample almost independently of applied force in the elastic region of bending deformation nN and has an almost constant value of approximately eV. At the same time, Ti4C3 nanoribbon is characterized by a shorter linear mode compared to Mo2Ti2C3; however, it is possible to estimate the approximate value of bending rigidity eV from the initial part of related dependence.