Proc. SPIE 5469, Fluctuations and Noise in Materials, 337 (2004)

Noise and Dislocation Formation in Atomistic Friction Simulations

C. J. Olson Reichhardt
T-12, Los Alamos National Laboratory, Los Alamos, NM, USA

ABSTRACT

We use large-scale atomistic simulations to study the work-hardening process that occurs when two metals slide against one another. Dislocations form at the interface between the work pieces and then migrate into the bulk. We examine the relationship between the velocity noise signature at the atomistic level and the number of dislocations present. We compare these signatures to those observed in a system of a single particle dragged through a lattice, where local melting can occur.

Keywords: friction, colloids, voltage noise, melting

1. INTRODUCTION
2. MODEL
3. LOW SLIDING VELOCITY
4. DISLOCATION NOISE
5. HETEROGENEITY AND NOISE
6. CONCLUSION
REFERENCES

1  INTRODUCTION

There has been considerable recent interest in the microscopic origins of friction [1], due to a combination of the availability of new experimental probes capable of resolving nanometer length scales and new simulation models running on high speed computers. Despite its long history and its technological importance [2], nanotribology remains poorly understood. Even the seemingly simple case of dry friction between two metal surfaces, when no lubricant material is present, produces a rich array of behavior such as work hardening, which occurs when the originally crystalline sliding surface develops a complex structure of nanocrystalline grains [3].
A number of simple models have been developed to treat dry friction in a simple way [4,5,6,7,8]. We focus here on a two-dimensional system in which the two materials sliding against each other are crystals composed of identical atoms but with differing orientations. In the lower surface, the close-packed crystal direction is parallel to the interface, while in the upper surface, the close-packed direction is perpendicular to the interface. (See Fig. 1.) Consider the motion of the atoms in the top layer of the lower surface. When the surfaces move against each other at relatively low velocities, the motion of the atoms in the lower surface can be approximated by a Frenkel-Kontorova model, in which the upper surface is replaced by a sinusoidal potential and the top layer of atoms in the lower surface are connected to each other by springs. A major limitation of this model is the fact that, in the actual system, the effective potential produced by the atoms in the upper layer is not static, but responds to the motion of the atoms in the lower layer. It may be possible to represent this effect by adding noise to the model; however, although the effect of noise on friction has been studied [9], the noise resulting from friction at the microscopic level is poorly characterized.
Here we use large-scale computer simulations to measure the noise experienced by a single atom moving at the interface at relatively low velocities. We also measure the noise in the dislocation density as a function of time as the sliding velocity increases. We compare the behavior of this friction system to the noise measured in a simulation of a single colloid driven through a colloidal crystal, where local melting occurs at the interface between the crystal and the moving particle [10].

2  MODEL

To represent two metals sliding against each other at high speeds in perfect vacuum, we use the SPaSM code developed by the Los Alamos group [11] to perform large-scale two dimensional non-equilibrium molecular dynamics friction calculations. Details of the geometry of the computational cell are given in Ref. [4,12]; a schematic is shown in Fig. 1. The cell is periodic in the direction of sliding (x direction), and has open boundaries in the y direction that are connected to thermostatted reservoir regions. In non-reservoir regions, the system evolves according to Newton's equations of motion. The velocity of the upper crystal is +up ∧x and the velocity of the lower crystal is −up ∧x, so that the relative velocity difference at the interface is 2up. The interactions between atoms take the form of a Lennard-Jones spline potential, described in detail in Ref [13,14]. We consider two sample sizes, a larger sample with 29064 atoms, and a smaller sample with 3448 atoms. The initial temperature is taken to be very low, well below the melting temperature of the crystal lattice.
Fig1.png
Figure 1: Schematic of the simulation geometry. The upper crystal is oriented perpendicular to the interface, while the lower crystal is oriented parallel to it. Both crystals are otherwise identical. A normal force Fn maintains the interface in contact, while tangential forces Ft maintain constant relative velocity of the two crystals in the x direction.
To measure the noise experienced by a single particle, we identify a single atom located on the interface between the two workpieces. We measure the acceleration in the x or sliding direction, ax, experienced by this atom as a function of time. We are also interested in the fluctuations of the number of dislocations present in the system as a function of time, since this quantity shows interesting properties near a 2D melting transition [15]. It is nontrivial to identify dislocations in a system with tens of thousands of atoms, and previous methods used approximate schemes based either on the local potential energy or the local crystal orientation surrounding an atom. We have adapted an algorithm given in Ref. [16] to perform a Voronoi tessellation on the entire system and directly identify atoms which do not have six neighbors. Our method can be used on 2D systems of arbitrary size, and we have tested it with up to 975,000 atoms. The defect identification is also highly valuable for data reduction, since we can record the positions of only the defected atoms, and construct high-resolution animations of the defects that contain only a fraction of the total number of atoms in the system.

3  LOW SLIDING VELOCITY

In Fig. 2 we show the locations of defects in a system with 3448 atoms as the sliding velocity is increased from up=0.2 to up=1.0. Sixfold-coordinated atoms are not plotted in the figure, but fill the remainder of the boxed area. At low velocities, the interface between the workpiece on top (with close-packed crystal direction perpendicular to the interface) and the workpiece on the bottom (with close-packed direction parallel to the interface) remains well defined and the behavior of the system is elastic. For higher velocities, local heating and plastic distortion begins to occur at the interface. Dislocations proliferate at the interface and begin to move outward into the surrounding crystal [5,17]. This has the effect of broadening the interface region, and results in local mixing of the atoms from the original crystals.
Fig2.png"
Figure 2: Positions of dislocations (atoms that do not have six neighbors) in a sample containing 3448 atoms at different sliding velocities shown after a sliding time of 80t0. The remaining atoms (not shown) fill the boxed region. The sliding velocity up=(a) 0.2, (b) 0.4, (c) 0.6, (d) 0.8, and (e) 1.0. In (a) the original interface between the two workpieces can be seen clearly. At higher velocities, dislocations proliferate at the interface.
We measure the acceleration in the sliding direction, ax, of a single particle located at the interface for velocities up < 0.5 in the elastic limit for a system containing 29064 atoms. The power spectrum S(ν)=|∫ax(t)e−2πiνtdt|2 of the single atom motion has a 1/fα signature, with α ∼ 3/4, as shown in the left panel of Fig. 3. We integrate S(ν) over the first octave of frequencies to obtain the noise power S0 [18,19], which we plot in the right panel of Fig. 3 as a function of sliding velocity up. We observe that the noise power increases roughly as the square of the sliding velocity in the elastic regime at lower up. This noise dependence could be incorporated into simple models for sliding friction in the elastic limit.
Fig3.png
Figure 3: Left: Representative power spectra of the acceleration in the sliding direction, ax, at up=0.1 (bottom curve), up=0.3 (middle curve), and up=0.5 (top curve), taken from a sample containing 29064 atoms. Right: Noise power S0 of the x-acceleration noise as a function of sliding velocity up.

4  DISLOCATION NOISE

As the driving speed is further increased, dislocations proliferate at the interface, as shown in Fig. 2(b-e). We measure the percentage of sixfold-coordinated atoms, P6, as a function of time for increasing sliding velocity in a sample containing 3448 atoms. The power spectrum S(ν) of P6(t), plotted in the left panel of Fig. 4, is broad, with a 1/fα signature where α ≈ 1.5. We find a dramatic increase in noise power when the velocity is increased above the elastic limit, as illustrated by the gap between the spectra in the left panel of Fig. 4. This is more clearly shown by plotting the corresponding noise power S0 (right panel of Fig. 4). The noise in the dislocation density increases monotonically with increasing sliding velocity once the dislocations begin to proliferate at the interface, which occurs for up > 0.2 in a sample of this size.
Fig4.png
Figure 4: Left: Power spectra S(ν) of the percentage of sixfold-coordinated atoms, P6(t), at velocities up= 0.05 (bottom), 0.1, 0.15, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, and 0.9 (top), in a system with 3448 atoms. A gap in spectral power occurs between up=0.2 and up=0.3. Right: Spectral power S0 for P6 as a function of sliding velocity up.

5  HETEROGENEITY AND NOISE

We compare the strongly heterogeneous structures observed in the high-speed dry friction simulations with heterogeneity introduced in a 2D colloidal crystal by driving a single colloid. Colloids are micron-sized particles suspended in a liquid solution. Each colloidal particle carries an electrostatic charge, leading to a repulsive interaction between the colloids with a range that can be adjusted by varying the salt concentration of the solution. The colloids can be confined to two dimensions by means of laser trapping, and can be manipulated individually using optical tweezers [20,21,22,23] or magnetic inclusions [24]. An important advantage of studying the colloidal system is that the colloids can be imaged directly using video microscopy techniques due to their large size. By studying colloids moving through a periodic substrate created by other colloids, it may be possible to gain insight into related phenomena such as the role played by dislocations in friction and depinning.
To model the colloids, we perform numerical Langevin simulations of the type employed in Refs [15,25,26]. We consider overdamped particles in a two-dimensional (2D) sample with periodic boundary conditions in the x and y directions. The initial configuration of the colloids is a triangular lattice. An additional colloid is placed in the system, and a constant driving force fd is applied only to that particle. The equation of motion for colloid i is
dri
dt
 =fij + fd + fT .
Here fij=−∑j ≠ iN ∇i U(rij) is the interaction force from the other colloids. The colloids interact with a screened Coulomb interaction,
U(rij)= qi qj
r
 exp(−κr).
(1)
Here qi(j) is the charge on colloid i(j) and 1/κ is the screening length, which can be adjusted experimentally by varying the salt concentration of the solution. We take κ = 2/a, where a is the lattice constant. All the particles have the same charge q=Q except the driven particle, which we take here to have q=3Q. The thermal noise fiT arises from random Langevin kicks with < fT(t) > =0 and < fTi(t)fTj(t′) > =2ηkB T δij δ(t−t′) Further details of the simulation appear in Ref [10].
When a single particle is driven along the zero degree angle with respect to the background lattice, we find three distinct phases for increasing temperature. For low temperatures, as in Fig. 5(a), the flow is elastic and the particle moves along a 1D path while causing small distortions in the surrounding lattice. Although a group of defects surrounds the driven particle, which is an extra or interstitial particle in the triangular lattice, there are no other dislocations induced elsewhere in the lattice. At intermediate temperatures, as in Fig. 5(b), we find a transition to plastic flow where the particle no longer moves strictly in a straight line, and defects form in the surrounding lattice. A localized molten region surrounds the moving particle, and ejects 5-7 defects into the lattice over time. Above the bulk melting temperature of the lattice, as in Fig. 5(c), defects are generated in the sample even in the absence of the driven particle. The colloid velocity drops once dislocations begin to proliferate in the lattice, as shown in Ref [10].
Fig5.png
Figure 5: Voronoi construction of the colloid configurations for a single colloid driven through a background lattice as the temperature is increased: (a) fT=0.25, (b) fT=0.45, and (c) fT=0.70. The Voronoi polygons are colored according to the number of neighbors: five (black), six (white), and seven (gray).
We measure the fraction of sixfold-coordinated particles, P6, as a function of time at different temperatures. Representative noise spectra from each of the three phases are illustrated in Fig. 6. The noise is broad with a 1/fα signature, where α ∼ 1.2. In previous work on a system with no driven particle, we observed a sharp increase in noise power of P6(t) upon crossing the melting temperature Tm [15]. When local melting is induced by a single driven particle, we find that this transition is rounded, as illustrated in the left panel of Fig. 7. The increase in noise power S0 as a function of temperature occurs at the same temperature where local melting begins, as indicated by the right panel of Fig. 7 where P6 is plotted as a function of T.
Fig6.png
Figure 6: Noise spectra S(ν) of the percentage of sixfold coordinated particles P6(t) at increasing temperatures: (a) fT=0.25 in the elastic flow regime, (b) fT=0.45 in the locally melted regime, and (c) fT=0.70 in the globally melted regime.
Fig7.png
Figure 7: Left panel: Noise power S0 of P6(t) as a function of temperature in the local melting system. Right panel: Corresponding average value of P6 showing the onset of local melting at fT=0.4.
This increase in the noise of P6 when a dislocated region forms around the driven particle is similar to the increase in noise observed in the friction system when dislocations begin to proliferate at the interface. It will be interesting to further explore possible connections between these two systems in the future.

6  CONCLUSION

We have studied the noise produced in two heterogeneous processes: the proliferation of dislocations at a friction interface for an atomic friction system at different sliding velocities, and the similar proliferation of dislocations in a colloidal crystal when a single colloid is dragged through the crystal at different temperatures. Both systems undergo a combination of effective local heating and local plastic flow. In both cases, we find broad noise signatures and an increase in the noise power as either the sliding velocity or the temperature increase.
The author thanks T.C. Germann, J.E. Hammerberg, B.L. Holian, C. Reichhardt, and J. Röder for useful discussions. This work was supported by the US Department of Energy under Contract No. W-7405-ENG-36.

References

[1]
B.N.J. Persson, Sliding Friction: Physical Principles and Applications, Springer, Heidelberg, 1998.
[2]
J. Krim, "Surface science and the atomic-scale origins of friction: what once was old is new again," Surf. Sci. 500, pp. 741-758, 2002.
[3]
T.L. Brown, S. Swaminathan, S. Chandrasekar, W.D. Compton, A.H. King, and K.P. Trumble, "Low-cost manufacturing process for nanostructured metals and alloys," J. Mater. Res. 17, pp. 2484-2488, 2002.
[4]
J. Röder, J.E. Hammerberg, B.L. Holian, and A.R. Bishop, "Multichain Frenkel-Kontorova model for interfacial slip," Phys. Rev. B 57, pp. 2759-2766, 1998.
[5]
J.E. Hammerberg, B.L. Holian, J. Röder, A.R. Bishop, and S.J. Zhou, "Nonlinear dynamics and the problem of slip at material interfaces," Physica D 123, pp. 330-340, 1998.
[6]
J. Röder, A.R. Bishop, B.L. Holian, J.E. Hammerberg, and R.P. Mikulla, "Dry friction: modeling and energy flow," Physica D 142, pp. 306-316, 2000.
[7]
B.L. Holian and J.E. Hammerberg, "Onset of incommensurate interfacial instability in a minimal model of dry friction," Phys. Rev. E 68, pp. 036101/1-4, 2003.
[8]
J. Röder, T. Germann, J.E. Hammerberg, and B.L. Holian, "Sliding between two incommensurate workpieces. II. Reduced dimensionality modeling," preprint.
[9]
Y. Braiman, H.G.E. Hentschel, F. Family, C. Mak, and J. Krim, "Tuning friction with noise and disorder," Phys. Rev. E 59, pp. R4737-4740, 1999.
[10]
C. Reichhardt and C.J. Olson Reichhardt, "Local melting and drag for a particle driven through a colloidal crystal," Phys. Rev. Lett. 92, p. 108301/1-4, 2004.
[11]
The SPaSM ("Scalable Parallel Short-Range Molecular Dynamics") code is described in: P.S. Lomdahl, P. Tamayo, N. Grønbech-Jensen, and D.M. Beazley, "50 Gflops molecular dynamics on the Connection Machine 5," Proc. Supercomputing '93, IEEE Computer Soc. Press, Los Alamitos, Calif., pp. 520-527, 1993; D.M. Beazley and P.S. Lomdahl, "Message-passing multi-cell molecular dynamics on the Connection Machine 5," Parallel Comp. 20, pp. 173-195, 1994; D.M. Beazley and P.S. Lomdahl, "Controlling the data glut in large-scale molecular-dynamics simulations," Comp. Phys. 11, pp. 230-238, 1997.
[12]
B.L. Holian, J.E. Hammerberg, and P.S. Lomdahl, "The birth of dislocations in shock waves and high-speed friction," J. Comp. Aid. Mater. Design 5, pp. 207-224, 1998.
[13]
B.L. Holian, A.F. Voter, N.J. Wagner, R.J. Ravelo, S.P. Chen, W.G. Hoover, C.G. Hoover, J.E. Hammerberg, and T.D. Dontje, "Effect of pairwise versus many-body forces on high-stress plastic deformation," Phys. Rev. A 43, pp. 2655-2661, 1991.
[14]
B.L. Holian and R. Ravelo, "Fracture simulations using large-scale molecular dynamics," Phys. Rev. B 51, pp. 11275-11288, 1995.
[15]
C. Reichhardt and C.J. Olson Reichhardt, "Fluctuating topological defects in 2D liquids: Heterogeneous motion and noise," Phys. Rev. Lett. 90, pp. 095504/1-4, 2003.
[16]
S. Fortune, "A sweepline algorithm for Voronoi diagrams," Algorithmica 2, pp. 153-174, 1987.
[17]
D.A. Rigney and J.E. Hammerberg, "Mechanical mixing and the development of nanocrystalline material during the sliding of metals," Advanced Materials for the 21st Century, The Minerals, Metals & Materials Society, 1999.
[18]
R.D. Merithew, M.W. Rabin, M.B. Weissman, M.J. Higgins, and S. Bhattacharya, "Persistent metastable states in vortex flow at the peak effect in NbSe2," Phys. Rev. Lett. 77, pp. 3197-3199, 1996.
[19]
M.W. Rabin, R.D. Merithew, M.B. Weissman, M.J. Higgins, and S. Bhattacharya, "Noise probes of underlying static correlation lengths in the superconducting peak effect," Phys. Rev. B 57, pp. R720-R723, 1998.
[20]
For a review, see D.G. Grier, "A revolution in optical manipulation," Nature (London) 424, pp. 810-816, 2003.
[21]
B.A. Koss and D.G. Grier, "Optical peristalsis," Appl. Phys. Lett. 82, pp. 3985-3987, 2003.
[22]
P.T. Korda, M.B. Taylor, and D.G. Grier, "Kinetically locked-in colloidal transport in an array of optical tweezers," Phys. Rev. Lett. 89, pp. 128301/1-4, 2002.
[23]
P.T. Korda and D.G. Grier, "Annealing thin colloidal crystals with optical gradient forces," J. Chem. Phys. 114, pp. 7570-7573, 2001.
[24]
P. Habdas, D. Schaar, A.C. Levitt, and E.R. Weeks, "Forced motion of a probe particle near the colloidal glass transition," preprint (cond-mat/0308622).
[25]
C. Reichhardt and C.J. Olson, "Colloidal dynamics on disordered substrates," Phys. Rev. Lett. 89, pp. 078301/1-4, 2002.
[26]
M.B. Hastings, C.J. Olson Reichhardt, and C. Reichhardt, "Depinning by fracture in a glassy background," Phys. Rev. Lett. 90, pp. 098302/1-4, 2003.


File translated from TEX by TTHgold, version 4.00.
 Back to Home