CC BY
23
NANOSYSTEMS: PHYSICS, CHEMISTRY, MATHEMATICS, 2011, 2 (3), P. 60-64
UDC 539.3
INFLUENCE OF SHEAR STRAIN ON STABILITY OF 2D TRIANGULAR LATTICE
E. A. Podolskaya, A. Yu. Panchenko, K. S. Bukovskaya
St. Petersburg, Russia katepodolskaya@gmail.com
PACS 46.25.Cc, 68.35.Rh
Stability of 2D triangular lattice under finite arbitrary strain is investigated. The lattice is considered infinite and consisting of particles which interact by pair force central potential. Dynamic stability criterion is used: frequency of elastic waves is required to be real for any real wave vector. Two stability regions corresponding to horizontal and vertical orientations of the lattice are obtained. It means that a structural transition, which is equal to the change of lattice orientation, is possible.
Keywords: stability, triangular lattice, finite strain, biaxial strain, pair potential, elastic wave, structural transition. 1. Introduction
In work [1] stability of plane triangular lattice under finite biaxial strain was investigated. Two stability regions, which correspond to vertical and horizontal orientations of the lattice, were obtained both analytically and using MD simulation. It was shown that taking more than one coordination sphere into account leads to a new effect: possibility of structural transition, which is equal to the change of lattice orientation. In this work shear strain is added. Modeling based on discrete atomistic methods [2] is proposed. The medium is represented by a set of particles interacting by a pair force central potential, in particular Lennard-Jones and Morse. Direct tensor calculus [3] is used.
Following [1,4,5], let us introduce the following notation to describe the geometry:
ak = rk-r0, (1)
where rk is radius vector of a particle k, r0 is radius vector of reference particle. If a lattice is simple, then any particle can be named "reference", each particle k has a pair —k and a_k = — ak. Triangular lattice is simple and close-packed: it coincides with its Bravais lattice and possesses maximum concentration of nods in elementary volume VO with the given minimum distance between the nods. Let us refer to the geometry which is described by rk and ak as reference configuration.
o
Let V and V be Hamilton's operators in reference and current configurations [3]:
A d ^ d
V = V = "eW (2)
Vectors <3j form an orthonormal basis. If vector r has projections .r,: in reference configuration, then in current configuration r will turn into R with projections X, in the same basis.
o
Suppose that the lattice is subject to strain characterized by V_R. According to long-wave approximation [2, 6]
Ak = E(r-Qk)-R (r) « ak-VR. (3)
Long-wave approximation takes into account those wave lengths that are much greater than the interatomic distance. The thermal motion is neglected.
Morse and Lennard-Jones potentials are used in this work to describe the interaction between particles
n(r) = D
r2^"1) - 2e-d^-1)
nLJ (r) = D
(4)
Here a is equilibrium distance in the system of two particles, D is the depth of potential well, d characterizes the width of the well. If B = 6, these potentials coincide in the elastic zone. Morse potential is preferable in this work, because, firstly, it decreases faster, so less particles may be taken into consideration, secondly, if r ^ 0 Morse potential remains finite.
Let Fk = F(Ak) = -n'(Ak) be interaction force and Ck = C(Ak) = n"(Ak) be the bond stiffness, both calculated in current configuration. Then we can introduce
J/, — J/. J/.'
44
AkAkAkAk,
$ = -
Fk Ak
Bk
A\
(Ck - $),
ik
$
Qk—k i
2
4
ik
lB
Qk—k—k—k
'E^ 4-
2
(5)
k
2. Stability criterion and deformation of triangular lattice
In the previous works [1,4,5] the following stability criterion was applied
Q > 0, (6)
where Q is determined from equation
det \D —
0.
Here
D=AC• • K, 4C = E<$> +AB, I< = I<K.
(7)
K is a real wave vector. This means that frequency of elastic waves is required to be real for any real wave vector.
6
1
1
? +
"> A.
Fig. 1. Reference and current configurations
Fig. 1 shows the typical part of triangular lattice before and after deformation. In reference configuration a1 = a2 = 60°. It is sufficient to take only 0 < p < 30° in account due to symmetry and infiniteness of the lattice. It was shown in [1, 5] that at least two coordination spheres should be considered.
In 2D case (7) takes the form
Q2 - QtrU + detU = 0. (8)
According to (6) roots of equation (8) are positive for stable current configurations. Thus, stability criterion is
tr22 > 0, det|2>0, 2trQ2 - (tr^)2 > 0. (9)
Inequality 2trD2 — (trD)2 > 0 is always true in 2D. Let G = he equations (9) yield
trD > 0 ^ GuKl + Gi2XiX2 + G22Kl > 0,
detD > 0 ^ AK\ + BK\Kl + CK"2 + DK\K2 + EKXK\ > 0,
where
G11 = C11 + C21, G12 = C14 + C24, G22 = G12 + C22, A = C11C21 — C41, B = 4C14C24 + C11C22 + ^12^21 — 2C41C42 — 4C24, C = C12 C22 — G42, D = 2C11C24 — 4C41C44 + 2C14 C21, E = 2C12C24 + 2C14C22 — 4C42 C44 •
(10)
(11)
Here Cij are the components of tensor4^.
The left part of tr D > 0 is a quadratic form in the components of the wave vector Kx and K2. It is positive definite, if
G11 > 0, 4G11G22 - G212 > 0. (12)
The left part of det D > 0 is a homogeneous polynomial of degree four. In this case, a general analytical criterion cannot be constructed.
Due to the fact that both K1 and K2 may be equal to zero, two necessary stability conditions are obtained, which help to narrow down the set of current configurations e1} e2, tp
A > 0, C> 0. (13)
Then, there are two ways to obtain sufficient conditions:
(1) For each e\, e2, we can construct detD, and check it for a set of K\ and K2 (Monte Carlo method). The inequality is homogeneous and even, so it is sufficient to consider only —1 < K1 < 1 and 0 < K2 < 1, which increases the efficiency.
(2) We can divide det D by and look into the problem of determining the coefficients so that a fourth-degree equation has no real roots, again for each e1, e2, ip. This method is much faster, but it causes a problem of distinguishing between complex and real roots, which leads to inaccurate results at the border.
In Fig. 2 stability regions, obtained by inequalities (12) and (13) and by the second method, are drawn. Here e1 and e2 are linear parts of Cauchy-Green tensor. There are several points, marked black, which were added by the first method. The stability regions are symmetric with respect to the plane tg p = 0. Two major areas correspond to horizontal and vertical orientations of the lattice [1]. Two additional small stability areas are connected with square lattices at p & 0o and p & 26o (see Fig. 3).
Let us draw a series of stress-strain diagrams, e.g. Fig. 4. According to [2] Cauchy stress tensor has the form
1
2F
g = (14)
k
Fig. 2. Stability regions
Fig. 3. Square lattices at p & 0° and p & 26°
Fig. 4. Pure shear (a11 = a22 = 0)
where V = v3/2(l + £i)(l + £2)- Grey zone in Fig. 4 corresponds to stability region, a 12 is diagonal component of Cauchy stress tensor.
In Fig. 4 we can see, that the loss of stability is strongly connected with the sign of the first derivative.
3. Concluding remarks
Stability analysis of 2D triangular lattice under finite arbitrary strain was carried out. In addition to [1] shear was taken into account. Two stability regions were obtained, when more than one coordination sphere were regarded, and a possibility of structural transition, which is equal to the change of lattice orientation, was noticed. Monte Carlo and analytical methods were used, and they proved to give practically equal results. Thus, Monte Carlo method can be applied to more complex cases, where it is impossible to accomplish analytical investigation, e.g. 3D lattices.
Acknowledgements
The authors are grateful to professor Anton M. Krivtsov for attention to this work and for useful discussions.
References
[1] Podolskaya E. A., Panchenko A. Yu., Krivtsov A. M. Stability of 2D triangular lattice under finite strain // Nanosystems: Physics, Chemistry, Mathematics. — 2011. — V. 2, No. 2. — P. 84-90.
[2] Krivtsov A. M. Deformation and fracture of solids with microstructure [In Russian]. — M.: Fizmatlit, 2007. — 304 p.
[3] Lurie A. I. Nonlinear theory of elasticity. — Amsterdam: North-Holland, 1990. — 617 p.
[4] Tkachev P. V., Krivtsov A. M. Stability criterion of internal structure of the material [In Russian] // XXXIII Science week, St. Petersburg State Polytechnical University, Part V, 2004. — P. 4-6.
[5] Podolskaya E. A., Tkachev P. V., Krivtsov A. M. FCC lattice stability analysis [In Russian] // XXXIX Science week, St. Petersburg State Polytechnical University, Part V, 2010. — P. 111-112.
[6] Born M., Huang K. Dynamical theory of crystal lattices. — Oxford: Clarendon Press, 1954. — 420 p.
