Plasticity

For large loads, solids can respond with irreversible deformation. One form of irreversibility is plasticity, that is carried by extended defects, the dislocations, in crystals. The module matscipy.dislocation implements tools for studying structure and movement of dislocations. Construction and analysis of model atomic systems is implemented for compact and dissociated screw, as well as edge dislocations in cubic crystals. The implementation supports ideal straight as well as kinked dislocations.

Creating an atomistic system containing a dislocation requires a rather deep knowledge of crystallographic system, elasticity theory as well as hands on experience with atomistic simulations. Thus, it can be challenging for people not familiar to the field. Within this module we attempt to create a flexible, friendly and pythonic tool to enable atomistic simulations of dislocations with ease. The base of the model is matscipy.dislocation.CubicCrystalDislocation class that contains all the necessary information to create an atomistic dislocation. To start experimenting, a user only has to choose the dislocation of interest and provide a lattice parameter and elastic constants.

Some dislocation systems feature additional complexity, as they can dissociate into two partial dislocations connected by a stacking fault. These kinds of dislocations are implemented within the subclass matscipy.dislocation.CubicCrystalDissociatedDislocation and follow a near-identical interface. matscipy.gamma_surface implements classes which assist in the modeling of stacking faults, as well as more general gamma surfaces. The classes can be initialised based on a dislocation system, or by known axes.

Installation and tests

matscipy.dislocation module relies on the anisotropic elasticity theory solution within Stroh formalism implemented in atomman package. atomman is not a part of default dependency of matscipy, however it can be easily installed with pip install atomman. To test your installation you can run pytest test_dislocation.py from tests directory of the repository. If atomman is not installed, most of the tests will be skipped resulting in the output similar to the following: ssssssssssssss..ssssssss..s.ssssssss.s where s stands for skipped test and . for passed tests. You can see which test are skipped together with the reasoning behind by running pytest in verbose mode with pytest -v test_dislocation.py. Once you have atomman installed you should have corresponding tests passing and should be able to use the code. If you plan to use LAMMPS and lammps python module there are couple of test that will test your installation when lammps python module is detected.

The majority of the tests (~70 %) are used during the development and adding new dislocations to the module. The idea is to check if Dislocation analysis (DXA) algorithm detects the same dislocation as we intend to create. These tests require OVITO python interface to be installed and are skipped if it is not detected. If you do not plan to add new structures to the module, you can safely ignore these tests.

Tutorials:

• Building cylindrical configurations with dislocations
  • Accessing Key Dislocation Properties
  • Available Dislocation Systems
  • Body Centered Cubic
    • 1/2<111>{110} screw dislocation
    • 1/2<111>{110} edge dislocation
    • 1/2<111>{110} mixed dislocation
    • <100>{110} edge ‘junction’ dislocation
    • <100>{001} edge dislocation
  • Face Centered Cubic
    • 1/2<110>{111} screw dislocation (perfect and dissociated)
    • 1/2<110>{111} edge dislocation (perfect and dissociated)
  • Diamond Cubic
    • 1/2<110>{111} screw dislocation (perfect and dissociated)
    • 1/2<110>{111} 60 degree screw dislocation (perfect and dissociated)
• Multi-species Dislocation Systems
  • Additional Dislocations in Multispecies Systems
  • Dislocations in Disordered Systems
• Dislocation Quadrupoles
  • Dislocation Glide in quadrupoles
  • Quadrupole Kink
  • Quadrupoles in other systems
    • BCC Screw in W
    • Dissociation of Screw dislocations in FCC Ni Quadrupoles
• Generating Gamma Surfaces & Stacking Faults
  • Method 1 - “Atom Move” Method
  • Method 2 - “Cell Move” Method
  • Gamma Surface Example using Method 1
  • Stacking Fault Example using Method 2
• Advanced Features of gamma_surface.py
  • Stacking Faults of Specific Dislocations
  • Accessing Different Stacking Fault Planes
  • Stacking Faults in more complex systems