A new Lagrangian formulation is introduced. It can be used to make molecular dynamics (MD) calculations on systems under the most general, externally applied, conditions of stress. In this formulation the MD cell shape and size can change according to dynamical equations given by this Lagrangian. This new MD technique is well suited to the study of structural transformations in solids under external stress and at finite temperature. As an example of the use of this technique we show how a single crystal of Ni behaves under uniform uniaxial compressive and tensile loads. This work confirms some of the results of static (i.e., zero temperature) calculations reported in the literature. We also show that some results regarding the stress‐strain relation obtained by static calculations are invalid at finite temperature. We find that, under compressive loading, our model of Ni shows a bifurcation in its stress‐strain relation; this bifurcation provides a link in configuration space between cubic and hexagonal close packing. It is suggested that such a transformation could perhaps be observed experimentally under extreme conditions of shock.

1.
E. R.
Cowley
,
G.
Jacucci
,
M. L.
Klein
, and
I. R.
McDonald
,
Phys. Rev. B
14
,
1758
(
1976
).
2.
M.
Parrinello
and
A.
Rahman
,
Phys. Rev. Lett.
45
,
1196
(
1980
).
3.
F.
Milstein
,
J. Appl. Phys.
44
,
3825
(
1973
).
4.
F.
Milstein
,
J. Appl. Phys.
44
,
3833
(
1973
).
5.
F.
Milstein
and
B.
Farber
,
Phys. Rev. Lett.
44
,
277
(
1980
).
6.
J. P. Hansen and I. R. McDonald, Theory of Simple Liquids (Academic, London, 1976).
7.
H. C.
Andersen
,
J. Chem. Phys.
72
,
2384
(
1980
).
8.
W. G.
Hoover
,
D. J.
Evans
,
R. B.
Hickman
,
A. J.
Ladd
,
W. T.
Ashurst
, and
B.
Moran
,
Phys. Rev. A
22
,
1690
(
1980
).
9.
L. D. Landau and E. M. Lifshitz, Theory of Elasticity (Pergamon, Oxford, 1959).
10.
Note that in the definition of S through Eq. (2.20) we have adopted a certain sign convention: S = +p corresponds to a system under hydrostatic pressure p.
11.
A. Rahman, in NATO Advanced Study Institute Series: Correlation Functions and Quasiparticle Interactions in Condensed Matter, edited by J. W. Halley (Plenum, New York, 1978).
12.
This procedure, used systematically in Ref. 2, is extremely useful for the purpose of resolving peaks in g(r) which otherwise would merge in neighboring peaks especially if the corresponding coordination numbers are very different. For example, in an hep structure the third and the fouth shells are very close to each other with a population, respectively, of 2 and 18 neighbors.
13.
In an MD calculation the “failure” of a system is manifested by a monotonic change in the hii and in an obvious disruption of the shell‐like structure of g(r). A monotonic drop in the value of the density is simply a consequence of the behavior of the hii.
14.
C. Kittel, Introduction to Solid State Physics (Wiley, New York, 1968).
This content is only available via PDF.
You do not currently have access to this content.