View source for LAMMPS Dislocation Mobility


==Objective==
In this tutorial, you will use the Modified Embedded Atom Method (MEAM) in Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) to acquire the dislocation mobility drag coefficient for your material.

You can obtain LAMMPS and find the user manual [http://lammps.sandia.gov/doc/Manual.html here].
Visit the [[LAMMPS tutorials]] page for a wide range of examples.
LAMMPS video tutorials - [https://youtu.be/GXA2PyqKYdY Intro] [https://www.youtube.com/watch?v=7Ila18g8zSY Input file Part 1] [https://www.youtube.com/watch?v=BOJPl9A7-K8 Input file Part 2]
[https://www.youtube.com/watch?v=TTDXJXJJi18/ Molecular Dynamics Theory and Examples]

==Environment Setup==

LAMMPS is an open source code and can be downloaded [http://lammps.sandia.gov here]. The source code can be easily compiled, or binary distributions are available for easy installation. 

In addition to the software, you will need an input file, an atom position file, a post-processing script, and two files for your atomistic potential.

* The input file for this simulation is [[Media:DislocationVelocity_PAD.txt | DislocationVelocity]].
* The atom position file will be created by [[Atomistic Dislocation Generation]].
* The post-processing script uses [[Ovito |Ovito's ]] scripting and analysis capabilities: [[Single Defect Velocity in Ovitos]].
* The atomistic potential files will be specific to your material.

==Calculating Dislocation Velocity==

===Walkthrough===

Video tutorial: [https://www.youtube.com/watch?v=S7w2wFHzX7o Part 1] [https://youtu.be/oZwmPV5mK0s Part 2]

====Step 1====
Generate the atomic structure using [[Atomistic Dislocation Generation]]. Once you have downloaded and compiled the fortran routine, as directed, simply run it at the command line as explained on the script page.
For this calculation, you want to make an edge dislocation in a PAD geometry. It will create the required data file as atoms.*.edge.pad, where * is either fcc or bcc, depending on your material. You will need to copy the data file to the directory where you will run the MD simulation.

====Step 2====
Edit the LAMMPS input file for your material, applied stress, and desired temperature. These can all be set by variables at the top of the input file.
<pre>
# Variable definitions
variable initTemp equal 300.          	# desired temperature
variable sigma equal 15000.    		# applied stress in bar
variable material string Ta		# material symbol
variable atom_file string atoms.bcc.edge.pad # the configuration was generated by SG with the preprocessor dislocation.f90
</pre>

<code>initTemp</code> will determine the temperature of the system. 

<code>sigma</code> will determine the applied shear stress in bar.

<code>material</code> will determine the name of your material used in your MEAM potential files. It expects to find a library file, <code>*.library.meam</code>, and a parameter file, <code>*.meam</code> where the * represents the material symbol. 

<code>atom_file</code> is the name of the data file that contains the initial atom positions.

====Step 3====

LAMMPS will create several output files from the simulation. You can load each of them in [[Ovito]] to see what each includes. You will run the post-processing script, [[Single Defect Velocity in Ovitos]], on the <code>dump.shear.unwrap</code> file, following the directions on the script page.

It will output a time versus position file which you can use to calculate the velocity of the dislocation.

== Go Back ==
*[[LAMMPS tutorials]]

== Acknowledgments ==


=== References ===

<references>
[[Category: LAMMPS]]
[[Category: Tutorial]]
[[Category: Metals]]