Author Topic: Lattice rotation  (Read 350 times)

rajdeep

  • Newbie
  • *
  • Posts: 3
Lattice rotation
« on: March 14, 2019, 02:03:38 PM »
Dear all,

I have a LAMMPS data file (unit cell) of TiAl which is a face centered tetragonal (FCT) lattice. In this data file my X, Y and Z axes are oriented along  <100>, <010>, and <001>, respectively.
I would like to rotate the lattice in such a way so that my new X, Y and Z are aligned with <-211>, <111>, and <01-1>, respectively. In order to accomplish this within OVITO, I applied the affine transformation option but the new lattice become triclinic. Please see the attached snapshot.

My question is does this affine transformation option can rotate my lattice in the desired crystallographic orientation? If so, after the transformation how I can wrap the atoms in a orthogonal cell?

Thanks,
Rajdeep     

Constanze Kalcher

  • Administrator
  • Sr. Member
  • *****
  • Posts: 299
Re: Lattice rotation
« Reply #1 on: March 15, 2019, 10:03:25 AM »
Dear Rajdeep,

note that the Affine Transformation modifier takes the transformation matrix that you enter and applies it onto your system. So what you need to enter there is the rotation matrix between your two coordinate systems.

-Constanze


rajdeep

  • Newbie
  • *
  • Posts: 3
Re: Lattice rotation
« Reply #2 on: March 18, 2019, 10:44:12 AM »
Dear Constanze,

Thank you very much for your reply. I would like to transform the coordinate system from (100, 010, 001) to (111, 1-10,11-2). So my rotation matrix is (111, 1-10, 11-2). If I provide this rotation matrix to OVITO along with 'Transform simulation cell' option, it gives me a triclinic box. Later on, I minimized the triclinic box in LAMMPS, and it seems that it has a large residual pressure even with stringent energy and force criterion which indicates that the rotated structure is not stable. In that case how I can wrap the rotated structure in a orthogonal box? Please see the initial and rotated structure. Thanks for your help.

Best regards,
Rajdeep   

Alexander Stukowski

  • Administrator
  • Hero Member
  • *****
  • Posts: 637
Re: Lattice rotation
« Reply #3 on: March 18, 2019, 11:32:33 AM »
Dear Rajdeep,

Note that your transformation matrix is not a proper rotation matrix: ([111], [1-10], [11-2]), because its column vectors are not normalized. A rotation matrix must have unit column vectors, all orthogonal, and its determinant must be 1. Since your matrix doesn't fulfill these criteria, you are getting something that is not a pure rotation.

It should suffice to normalize each column vector to obtain a rotation matrix, i.e.

([111]/sqrt(3), [1-10]/sqrt(2), [11-2]/sqrt(6))

-Alex

rajdeep

  • Newbie
  • *
  • Posts: 3
Re: Lattice rotation
« Reply #4 on: March 27, 2019, 04:45:14 PM »
Dear Alexander,

Thank you for your email. Sorry for incorrect rotation matrix.
Now, I'm providing the normalized rotation matrix and the option "transform simulation cell" is selected. In the output it looks like the entire simulation box is rotated with respect to the X, Y and Z axis, but atom positions are not changed. I did minimization of the rotated structure and got exactly similar structure as the input structure was. It is expected that after rotation cell size along X, Y and Z should change, but it is not changed in the rotated system. Please find the attached figures of the initial and rotated structure. Am I missing anything or misunderstanding OVITO?

Thank you very much for your help.

Best regards,
Rajdeep
   
 

Constanze Kalcher

  • Administrator
  • Sr. Member
  • *****
  • Posts: 299
Re: Lattice rotation
« Reply #5 on: March 28, 2019, 05:22:46 PM »
Hello Rajdeep,

instead of rotating everything you need to apply the rotation only on the atoms.

-Constanze