Potential parameters for W.

## Abstract

Fundamentals of the atomic mechanisms governing the strength of nanosized metallic crystals are described. An attempt is made to explain on this basis the size and orientation effects, temperature dependence of strength and atomism of fracture of bcc crystals under triaxial uniform (hydrostatic) tension. A feature of the proposed material is that it combines the results of molecular dynamics simulation with the data of experimental research findings on failure of metallic nanosized crystals under the high-field mechanical loading. It is exhibited that local instability of the lattice is the main mechanism governing the strength of defect-free nanosized crystals (NSC). Based on the concept of local instability, an explanation is given of the nature of the size effect in NSC, as well as of the differences in its manifestation in nanocrystals with bcc and fcc lattices. The concept of the mechanism of thermal activation of local instability is outlined. This enables to explain the specific features of the temperature dependence of NSC. The results of experimental studies and molecular dynamics simulation of the failure of tungsten nanocrystals under hydrostatic tension are presented. The ideas about the atomism of the bcc-fcc transition in these conditions are articulated.

### Keywords

- molecular dynamic simulation
- nanosized crystals
- lattice instability
- strength
- nanowire
- size effect
- surface tension

## 1. Introduction

Significant advances in the fabrication of nanoscale objects with controlled size and properties (nanorods, nanowires, nanopillars and carbyne chains) are a characteristic feature of nanotechnology development in the last decade. In turn, this stimulates the development of experimental methods for characterising the properties of these nanoobjects, as well as computer methods for predicting their behaviour under various external influences. Nanosized metallic crystals occupy a special place among these objects. Due to their remarkable physical and mechanical properties, they can be directly used in various nanodevices; in addition, they are considered as building blocks for the creation of nanocrystalline materials. On the other hand, study of nanosized crystals enables to ascertain the fundamental mechanisms of deformation and fracture inherent in the nanoscale. Fabrication of defect-free nanocrystals and their compression test [1, 2, 3], as well as the development of a high-field technique for the preparation and *in-situ* tensile testing of nanoneedles [4, 5, 6], becomes a milestone event in experimental studies of nanosized crystals. Absence of defects (dislocations and twins) in these crystals made it possible to reach extremely high levels of strength, which are close to the value of “theoretical strength”. Therefore, molecular dynamics (MD) simulation is the most effective tool for studying atomic mechanisms to reach extremely high levels of strength in crystals.

Nanosized crystals of bcc metals are a classical object of MD simulation [7, 8, 9, 10, 11, 12, 13]. Currently, considerable attention is paid to metals with a bcc lattice [14, 15, 16, 17, 18, 19, 20]. From a physical point of view, bcc metals are more interesting because they provide a greater variety of atomic rearrangements, governing the strength. However, the main difficulty in MD simulation of bcc metal is lack of sufficiently reliable potentials, which limits the accuracy of such simulation. In this connection, the appearance of direct experimental data on the strength of bcc nanocrystals substantially facilitated the problem of selection and verification of these potentials. A great number of works on MD-simulation enabled to predict the key effects controlling the level of strength of nanosized nanocrystals, namely: (i) the size effect, (ii) dependence of strength on orientation, (iii) the temperature dependence of strength, and (iv) dependence of strength on the stress state mode [7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24]. At the same time, the absence of a completed theory of strength of nanosized crystals is a characteristic feature of science of the strength of these objects. Lack of sufficient experimental data is the second feature of state-of-the-art in researches of strength of nanosized crystals. Difficulties in performing mechanical tests of nanospecimens are the reason for this. In this regard, the high-field mechanical tests of nanosized specimens should be indicated. It has made it possible *in situ* mechanical loading of nanosized specimens up to level close to theoretical strength combined with direct observation of the atomic structure of these objects under well-controlled crystallographic conditions [25, 26].

In the present work, fundamentals of the atomic mechanisms governing strength of nanosized metallic crystals are described; besides, an attempt is made to explain on this basis the size and orientation effects, temperature dependence of strength, and atomism of fracture of bcc crystals under triaxial uniform (hydrostatic) tension.

## 2. Methods of physical and numerical experiments

### 2.1. Field-induced mechanical test

The *in situ* mechanical loading of nanosized crystals in the FIM experiments was realised using the Maxwell mechanical stress induced by high electric fields. As was shown in these studies, the vast majority of dislocations was removed under the action of the electric field–induced mechanical stresses. So, most of the dislocation-free nanosized crystals were fractured at mechanical stresses corresponded to substantial parts of the ideal strength of solids. This technique was used for uniaxial tensile tests of tungsten and molybdenum nanoneedles and for uniform triaxial (hydrostatic) tensile tests of a tungsten nanocrystal.

The intrinsic hydrostatic tensile strength of tungsten was experimentally determined using a high-field testing in field-ion microscope (FIM) with needle-shaped specimens (Figure 1a) [6]. The specimens with an initial radius of curvature of about 100 nm and the taper angle of

where ^{−1} related to an evaporation rate of ~10^{−2} atomic layers per second. The electric field decreases severely under the hemispherical cap in the conical part. So, the failure of needle-shaped specimens was usually initiated at the nanosized cap in the region of the ultimate field strength. An electric field applied to a metal specimen induces the Maxwell surface stress

where

where

For nanocrystal failure under multiaxial uniform (hydrostatic) tension, a protrusion was formed on the spherical part of specimen. It was created by selection of special schedules for the high-field evaporation of atoms from the nanospecimen surface [19]. This is shown schematically in Figure 2a. Electric field strength and, accordingly, mechanical stresses increase with decrease in the curvature radius of a surface. So, maximal stresses act in this protrusion. This makes it possible to fracture the metal exactly in this volume. It is done by a steep increase in the electric field strength. Appearance of the crater on a nanoneedle surface is the direct evidence that fracture *occurs* in the local region, where hydrostatic tension acts (Figure 2b).

The Maxwell stress was determined with accuracy of ±4%. The systematic error due to uncertainties in the field calibration is ±3% (see [28, 29]), and, correspondingly, the systematic error in determination of the Maxwell stress is ±6%. The total error in determination of the Maxwell stress is ±10%.

### 2.2. Molecular dynamics simulation

To analyse atomic mechanisms of instability and failure of defect-free nanosized crystal under uniform triaxial (hydrostatic) tension, MD simulation was employed. Simulation was performed using the software (program package) XMD [

Extended Finnis–Sinclair semiempirical potential [30] (Tables 1 and 2) was utilised for MD simulation of hydrostatic tension of tungsten. This potential overcomes the “soft” behaviour of the original Finnis–Sinclair potential [31].

^{−1}) | ^{−2}) | ^{−3}) | ^{−4}) | ^{−5}) | ^{−6}) | ||
---|---|---|---|---|---|---|---|

1.885948 | 4.41 | 3.25 | 48.527 96 | −33.79621 | 5.854334 | −0.0098221 | 0.033338 |

^{−1}) | ^{−2}) | ^{−3}) | ^{−4}) | ^{−5}) | ^{−6}) | ||
---|---|---|---|---|---|---|---|

1.848648 | 4.1472 | 3.2572 | 47.98066 | −34.09924 | 5.832293 | 0.017494 | 0.020393 |

According to EAM or FS formalism, the total energy of a system is given by the expression:

where the first term is the conventional central pair-potential summation:

and the second term is the n-body term:

Uniform triaxial tension of ball-shaped nanosized tungsten crystal was also simulated. Initial radius of nanocrystal was 10.8 nm. Such a size enabled to eliminate the influence of rigid surface shell on the behaviour of atoms in the inner region of specimen. To describe the distribution of atomic velocities for given temperature, the Boltzmann distribution was used. To hold simulated temperatures, velocity of atoms was scaled by ^{8} s^{−1}. Tension of nanowires of molybdenum and tungsten in three crystallographic directions [100], [110] and [111] was modelled also. Cylinder with the axis oriented along a corresponding crystallographic direction was cut from the set of atoms placed exactly at a lattice site. Relation between the diameter and the length of cylinder is one-quarter. After that, boundary conditions are applied. The central part (gauge length) of the cylinder with length of 3 diameters is still free, but atoms of border parts are made “frozen.” Time step is 10^{−15} sec. The diameter of nanospecimens varied within the range 1÷13 nm. In addition, to interpret the experimental data on the destruction of nanopillars under uniaxial compression, MD simulation of the compression of cylindrical Mo specimens in the [100] direction was performed.

For the quantitative analysis of the driving forces of atomic rearrangements, in addition to global stresses

where *k* refers to the considered atom; *m* refers to the neighbouring atom; *k* and *m* determined as the gradient of energy functional; *k* and *m*; *n* is the number of the nearest neighbouring atoms; *i*, *j* are the stress tensor indexes.

An expression for the global stresses

where *N* is the total number of atoms in the nanowire.

For calculation of the local shear stresses, the following formula was used:

where

True strain value was employed as a measure of longitudinal strain of the nanowire as a whole:

where

When modelling hydrostatic tension of a nanosized ball-shape specimen, engineering strain value was determined as:

where

## 3. Results and discussion

As it is known, microscopic defects (dislocations and twins) are the main reason for decrease in strength of materials. Modern nanotechnologies enable us to create materials, which, in most cases, are defect-free. Therefore, initially, it was believed that their strength must tend to the ideal strength and be the intrinsic material trait much like the elastic constants. However, it was appeared that the strength of nanosized crystals changes within the wide range depending on their sizes, temperature, loading condition, etc. However, the concept of ”ideal strength” can be used as a starting point for the analysis of *atomic mechanisms* governing the strength of nanosized crystals. The ideal strength of a material is defined to be the maximal *homogenous* stress that an ideal crystal can withstand. Reaching of the limit state of such crystal is related to uniform and simultaneous break of atomic bonds or its reformation. In general case, two instability modes exist for metals, namely (i) instability of crystal under the tensile stresses (instability on the Bain path for bcc crystals) and (ii) shear instability. The first mode of instability of an ideal crystal results in its disintegration. Shear instability gives rise to change in shape and orientation of an ideal crystal [33, 34, 35].

A fundamental difference of nanosized crystals from the ideal ones is that the nanocrystals can be defect-free but not ideal since the ideal position of atoms in the crystal lattice is disturbed by both (i) action of surface tension forces and (ii) thermal vibrations of atoms. Action of the surface tension forces leads to inhomogeneous distribution of local stresses. Thermal vibration of atoms causes local stress fluctuations. These two factors are the main reasons for *localisation* of process of breaking and reformation of atomic bonds in nanosized crystals [36, 37]. As a result, the above two modes of an ideal lattice instability occur in a nanosized crystal, but this instability is realised in a limited local region of the crystal. This gives rise to the fact that defects form as a result of above instability. For instance, formation of dislocations and twins is a direct consequence of this localised shear instability. The key difference of these defects from those in macrosized single crystals lies in the fact that they are *highly* non-equilibrium, since they are formed in a defect-free crystal at stresses much exceeding the critical stress of their propagation in a crystal (Figure 3). As shown in Figure 3, the stress of formation for these defects predetermines the level of strength of nanosized crystals [37].

### 3.1. The effect of surface

A characteristic feature of nanosized objects is the significant influence of the surface on their properties. For mechanical properties, the surface effect is realised through the surface tension. As it is known, in nanosized specimens, the level of stresses created by surface tension may reach significant values even without loading [8]. Surface tension forces induce tensile stresses in the surface layer of the crystal, which are balanced by compressive stresses in the bulk. It gives rise to *inhomogeneous* local stress distribution within interior of the nanosized specimen. Figure 4 presents the cross-sectional distribution of local shear stresses _{,} which act in shear systems, where instability of crystal is observed. According to these data, even in the case of the unloaded nanowire, the value of shear stresses in a surface layer of Mo nanowire may reach values of the order of 4 GPa. This is only 3 times less than theoretical shear stress for Mo. At loading of such specimen, these stresses cause the lattice instability in sub-surface layer (Figure 4d). Surface tension not only gives rise to localization of the instability region in surface layer, but it is also the reason for size effect, which manifests itself in growth or decrease in strength of nanosized crystals (Figure 5). Recently, two approaches exist to explain size effect. The first of them is based on a “global” criterion of limit state [8, 10, 16], and the second one employs the “local” criterion of nanocrystal instability [14, 15, 18, 19].

Because of the surface tension, an inner part of a nanosized crystal is compressed, and the thin surface region is stretched. If thickness of this stretched layer is neglected, then compressive stresses

where *d* is the diameter of specimen.

Global approach supposes that at compression of a nanosized specimen, stress *local* shear stresses

where

To derive the criterion of nanocrystals instability within the framework of continuum approach, the concept of “effective” thickness

where

Using the criterion (10) with (11) gives the first approximation for strength

According to this dependence, at _{,} the surface tension must give rise to decrease in strength

The third term in this dependence accounts for the effect of decrease in total value of shear stresses influencing the twin, which is due to the fact that part of twin is located in a compressed region. As it is exhibited in Figure 5, dependence (15) enables to describe the effect of increase in strength at tension in the direction [100]. Thus, this rather simple model enables to explain not only the anomaly of the size effect at tension of nanowires of bcc metals in the <100> directions but also the difference between the regularities of size effect in fcc and bcc metals. For all three orientations <100>, <110>, and <111>, local instability in fcc metals is due to formation of a non-equilibrium stacking fault, the critical size of which exceeds a thickness of a tensile surface layer. As a result, for these three crystallographic orientations, the effect of increase in

### 3.2. The temperature effect

As noted above, there are two main reasons for localisation of the lattice instability in nanocrystals, namely: (i) the effect of surface tension and (ii) the local stress fluctuation due to oscillations of atoms in the lattice. Statistical distribution of local shear stresses in nanosized crystal of iron at two temperatures is exhibited on Figure 6. Limit state of a nanosized crystal is reached when the local shear stresses

The fluctuation component

where

Thus, temperature dependence of the local stress fluctuations *nucleation* of defects. In the second case, the lattice vibration facilitates *mobility* of already existent defects. It is necessary to emphasise the fundamental difference in the mechanisms and structural levels of realisation of these processes. In the first case, we are talking about the thermal activation of instability of the crystal lattice (atomic-scale event). In the second case, the oscillation of the dislocation line (a microscopic object) plays a key role, leading to the formation of kinks pairs, by means of which a thermally activated overcoming by dislocation of the potential relief occurs. As a result, the laws of temperature dependence differ as well as change in the absolute strength value. For instance, over the temperature range 77К…300К, the value of yield strength of typical bcc transition metals decreases 3–5 times, while the changes of a critical stress of NSC instability do not exceed 15–25%.

### 3.3. Nanosized crystal instability under hydrostatic tension

As noted above, there are two modes of crystal instability, namely shear instability and instability under the action of tensile stresses. In ideal bcc crystals, this is manifested in the existence of instability on the orthorhombic deformation path and instability on the Bain path. As mentioned above, a shear instability in nanowires is usually realised at tension in directions <110> and <111>. In the general case, when the bcc crystals are stretched in the <100> direction, *ab-initio* computation predicts the possibility of instability, both on the orthorhombic and Bain paths. As a result, the mode of instability is realised, for which a lower value of strains is required. So, for example, for an ideal niobium crystal, the value of critical strain of instability on the orthorhombic path is almost 1.5 times less than the corresponding strain on the Bain path [34]. This means that the strength of Nb nanowires at tension in the direction <100> should be governed by local shear instability. The situation with Mo, W, and Fe is more complicated. According to the results of *ab-initio* calculations, for ideal crystals of these metals, the values of critical strains and stresses on both paths are very close [33, 34, 35]. However, the results of the MD simulation on tension of nanowires of these metals indicate the realisation of local shear instability. One of the results of such instability is the reorientation of the crystal lattice from [100] to [110]. According to the *ab-initio* calculations, such a reorientation should be a consequence of the instability on the orthorhombic path. The only difference is that in an ideal crystal, this reorientation must occur simultaneously in the entire crystal, and in a nanosized specimen, this occurs locally with subsequent propagation to the entire volume of nanospecimen (Figure 8). These regularities were considered in the works [14, 40]. Later, similar results were obtained on nanowires with a square cross-section [16]. As a result of realisation of the shear instability mechanism at tension of Mo, W, and α-Fe nanowires in the direction [100], brittle fracture (Bain instability) does not occur, but on the contrary, the plastic strain is larger as compared with that for orientations [110] and [111] [14, 40]. In this connection, the question arises whether it is possible to realise instability on the Bain path in *nanosized* specimens made of the bcc metals. Ваіn instability was found in [14] at MD-simulation of hydrostatic tension of ball-shaped specimen of iron. Under hydrostatic tension, the value of the average shear stresses is zero. This contributes to the realisation of the Bain instability, since there is no driving force for the shear instability of the crystal. However, the bcc→fcc transition requires the realisation of a certain relationship between the strains of lattice along the directions OX, OY, and OZ [34]. This is impossible to be realised under conditions of *uniform* triaxial tension. To analyse the details of this mechanism of bcc→fcc transition in nanosized crystals, hydrostatic tension of a molybdenum nanoball with a diameter of 10.8 nm was simulated [19]. MD simulations were performed for T = 77 K (the temperature for which the experimental data were obtained) and also for a lower temperature of 30 K. Simulation at lower temperatures makes it possible to obtain a clearer picture of atomic rearrangements, since thermal “smearing” of atom locations is reduced. Based on the results of this simulation, the diagram of deformation of the ball was built (Figure 9). To analyse local bcc-fcc rearrangements, for each atom, the number of neighbouring atoms in the first coordination sphere was calculated. This permits to determine the regions in crystal with both initial bcc lattice and formed unstable fcc lattice, as well as heavily deformed regions (Figure 9).

According to the MD-simulation findings, when critical stress of instability of a specimen under uniform triaxial tension,

To analyse quantitatively the instability initiation for the temperature 77К, distribution of local tensile stresses *uniaxial* tension *ab-initio* calculations for iron presented in [42] show that when applying a tensile stresses in the orthogonal direction, the Bain instability stress increases as compared to purely uniaxial tension. In our case, the Bain transition occurs in one of the directions <100> as a result of *fluctuation* of *local* tensile stress in this direction. This is the feature of the atomic mechanism of instability initiation in nanosized specimen under global uniform triaxial tension. It should be emphasised that in this case, the appearance of fluctuations gives rise not only to excess of the local stresses over the average ones but also to deviation from uniform triaxial tension within the *local region* where the Bain transition is realised.

For quantitative estimates, the expression for value of local stress of Bain transition,

where

The probability of *simultaneous* fluctuations of local stresses in three mutually orthogonal directions is many orders of magnitude less than probability of fluctuation in *only one* direction, so, in the first approximation, at the moment of local instability:

where

Respectively:

Dependence (20) enables to estimate the value of local stress of initiation of local instability,

To estimate *ab-initio* calculations, the value of critical strain of the lattice instability increases at transition from uniaxial to triaxial tension, and at uniform triaxial tension, it reaches maximum value of 0.15 [42]. It means that obtained value

According *ab-initio* calculation findings, the instability stress of an ideal Mo crystal under uniform triaxial tension is 50.0–52.0 GPa [42]. However, such a value of strength can be reached only in an ideal crystal but not in nanosized specimens. Fluctuations in local tensile stresses are the reason for this. They do not enable the uniform triaxial tension in *local* regions of the crystal to be realised at a *global* (on average) *uniform* triaxial tension of the entire specimen. Probability of the fluctuation of local stresses in three directions simultaneously is much lower than the probability of fluctuation along one of the directions <100>; so, initiation of instability of nanosized tungsten crystal occurs at local stresses

Thus, the strength of nanosized crystals under hydrostatic tension will always be below the strength of an ideal crystal under the same conditions. In this case, the strength of tungsten nanosized crystals under hydrostatic tension is 28 GPa. This is approximately 1.90 times less than the strength of an ideal tungsten crystal. The reason for this is fluctuations of local tensile stresses, which leads to: (i) exceeding of the value of local stresses over global (fluctuations themselves) and (ii) to deviation from triaxial uniform tension within *the local region* where this transition realises. For nanosized Mo specimens, the first factor results in a decrease in strength by 1.43 times, while the second one—by 1.25–1.30 times.

### 3.4. Strength of nanopillars and nanoneedles

Nowadays, nanoneedles and nanopillars are the most common nanosized specimens, which strength can be determined experimentally. Two main kinds of nanopillars exist, namely: (i) nanopillars obtained by focused ion beam (FIB) technology [43, 44], and (ii) nanopillars obtained by etching of nanocomposites [3]. Surface layer of specimens obtained by FIB technique contains a great number of defects in sub-surface layer, and so their strength does not exceed 1.0–1.5 GPa [43]. The second kind of specimens has no such shortcomings, and so, their strength is much higher. For instance, strength of Mo nanopillars of diameter 300–1000 nm is approximately 9 GPa. This is 6 times greater than strength of nanopillars obtained by FIB technique. However, this strength value is more than 2 times less than the evidence on MD-simulation of compression of defect-free Mo nanopillars of this crystallographic orientation [100]; according to our data, this value is

A high-field technique of tensile testing *in-situ* of nanoneedles is free of such disadvantage. In addition, it makes it possible to test the specimens of significantly smaller diameters—from 20 till 125 nm. Data presented in Figure 12 enable to estimate experimentally strength of Mo and W nanocrystals by the value of upper scatter limit of experimental evidence on failure of nanoneedle specimens. This value for W is 23.2 GPа, and for Mo, it is 20.0 GPа. These figures also present the results of MD simulation of tension of cylindrical specimens in the direction [110]. Calculated strength values are approximately 1.4 times greater than the maximum experimental values of the nanopillars strength. At high-field treatment of the tip of nanoneedle specimen, its working volume is “cleaned” from dislocations [4, 19], and so, both decrease the values of strength of nanoneedle specimens, and their scatter may be because of stress raising related to rough lateral surface, which is formed by electrochemical polishing. The results obtained make it possible to identify the main factors leading to a decrease in the strength of nanosized crystals. According to these data, tensile strength of defect-free nanowires is always less than ideal strength. This is due to both surface tension effect (effect of a physical surface) and thermally induced local stress fluctuation (temperature effect). At transition to nanoneedles, we have additional effect of a surface roughness. This effect may give rise to a reduction in strength by at least 30–40%, and it is the cause for considerable scatter of experimental data. The strength of nanopillars is significantly influenced by the buckling effect.

## 4. Conclusions

On the atomic scale, the strength of nanosized crystals (NSC) is governed by two modes of lattice instability. This is instability with respect to the action of tensile stresses and shear instability. In nanosized crystals of fcc metals, the second mode of instability predominates. BCC metals are characterised by a competition between these two modes of instability, depending on the loading conditions of the crystal. Instability on the Bain path is realised at uniform triaxial (hydrostatic) tension. At uniaxial tension, even under cryogenic conditions, the strength of the bcc nanocrystals is controlled by shear instability (instability on the orthorhombic path).

Localisation of instability in a limited volume of a nanospecimen is the distinguishing feature of the mechanism of nanosized crystal instability. It occurs due to two main factors: (i) local stresses fluctuation caused by thermal vibrations of atoms and (ii) the effect of surface tension in NSC. The first factor governs regularities of strength dependence on temperature and the feature of the bcc→fcc transition under hydrostatic tension. The second one influences the absolute value of nanocrystal strength and determines the main regularities in manifestation of the size effect in nanospecimens of bcc and fcc metals, as well as the orientation dependence of the size effect in nanospecimens of bcc metals.

Dependence of the magnitude of local shear stresses in the surface layer of nanospecimen, induced by the surface tension, on its diameter (cross-sectional size) is the reason for existence of a size effect for the strength of nanospecimens. The “sign” of a size effect (increase or decrease in strength with decreasing diameter) is determined by the kind of non-equilibrium defect that is formed as a result of local instability of the nanocrystal.

Strength of nanospecimen (nanowire) increases with decrease in its diameter (cross-sectional size) if the value of this strength is controlled by the formation of a stacking fault (fcc lattice, orientations <100>, <110>, <111>) or a twin (bcc lattice, <100>). Formation of non-equilibrium dislocations (bcc lattice, <110>, <111>) gives rise to an anomaly of the scale effect, i.e., to fall in strength with decreasing diameter. This is due to the fact that non-equilibrium dislocations form in a thin surface layer, where the action of surface tension gives rise to an increase in the level of local shear stresses, i.e., promotes realisation of the local instability. Size of the region required for the formation of a non-equilibrium stacking fault or twin exceeds thickness of the stretched surface layer, and this region falls within the range of compressive stresses. The latter reduces the magnitude of the resulting local shear stresses in the region where a twin or a stacking fault is formed, which leads to an increase in the strength of nanospecimen.

Physical mechanisms of the thermally activated reduction in strength of nanosized and macrosized crystals are fundamentally different. In first case, the thermal vibration of atoms causes local instability in defect-free structure and

*formation*of non-equilibrium defects. In the second case, the atom vibrations raise*mobility*of already existent defects. As a result, the laws of temperature dependence differ as well as change in the absolute value of strength. For instance, over the temperature range 77К…300К, the value of yield strength of typical bcc transition metals decreases 3–5 times, while the changes in a strength of nanosized crystals do not exceed 15–25%.Instability in bcc nanocrystals on the Bain path is possible under uniform triaxial tension of nanosized specimens. The peculiarity of the manifestation of “localised instability” under these conditions is that the fluctuations of local tensile stresses not only lead to a decrease in the global (average) stress required for the bcc→fcc transition but also cause deviation from triaxial uniform tension within

*the local region*where this transition realises. It means that at a global uniform triaxial tension of nanospecimen, locally, bcc→fcc transition occurs under non-uniform triaxial tension. This leads to a decrease in the magnitude of stress required for such a transition. For W at T = 77 K, this decrease is 25–30%.

The maximum attainable experimental values of the strength of Mo and W nanoneedle specimens in the crystallographic direction [110] are approximately 40% less than the results of MD simulation. This may be caused by stress raising due to rough lateral surface of nanoneedle specimens forming as a result of electrochemical polishing. This means that formation of an atomically smooth surface of nanoneedles is one of the factors to reach the ultimate strength levels.

## Conflict of interest

The authors declare that they have no competing interests.

## Funding

This work was supported by the National Academy of Sciences of Ukraine [grants number #0117 U002131; #0117 U006351].