Author Topic: how to calculate Na diffusion coefficient in Si bulk  (Read 4018 times)

0 Members and 1 Guest are viewing this topic.

Offline nadeemqasr

  • New QuantumATK user
  • *
  • Posts: 3
  • Country: kr
  • Reputation: 0
    • View Profile
how to calculate Na diffusion coefficient in Si bulk
« on: December 18, 2015, 15:01 »
I am very new user to ATK and especially beginner to DFT studies. Could anyone please suggest me tutorial or example to calculate the diffusion coefficient for a case when Na ion diffuses into Sn bulk crystal structure? I have heard about NEB method and MD simulation but not know how to do that. Any response will be appreciated.

Offline nadeemqasr

  • New QuantumATK user
  • *
  • Posts: 3
  • Country: kr
  • Reputation: 0
    • View Profile
Re: how to calculate Na diffusion coefficient in Si bulk
« Reply #1 on: December 20, 2015, 02:10 »
Actually i can see a tutorial that tells how to calculate the barrier that barrier is required for diffusion calculations. Here is the problem that i do not know how to do that. I have seen some papers who use some kind of mean square displacement (r) and simulation time (t) etc. I can share these text or papers.
please do guide me details if anyone knows.

Offline zh

  • QuantumATK Support
  • Supreme QuantumATK Wizard
  • *****
  • Posts: 1141
  • Reputation: 24
    • View Profile
Re: how to calculate Na diffusion coefficient in Si bulk
« Reply #2 on: December 21, 2015, 00:58 »
When you do the MD simulation of your system, you save the trajectory (i.e., the  atomic coordinates) of every MD simulation step (or, every several MD stimulation steps) and then calculate the mean square displacement from these trajectory.  In the last step, you have to do a bit but easy coding. If you can't do it, you may find available tools from internet:
http://www.pages.drexel.edu/~cfa22/msim/node29.html
http://www.multiscalelab.org/utilities/DiffusionCoefficientTool
https://nanohub.org/resources/22854
http://www.mdanalysis.org/

Offline nadeemqasr

  • New QuantumATK user
  • *
  • Posts: 3
  • Country: kr
  • Reputation: 0
    • View Profile
Re: how to calculate Na diffusion coefficient in Si bulk
« Reply #3 on: December 21, 2015, 09:16 »
I am very thankful for your reply. I am now practicing the KMC/MD simulations for trajectory save. and meanwhile studying the links you provided. Just in case,could you please tell me how to save the trajactory (.nc) file to vasp DATACAR file? or other formats to see the details of positions during the diffusion?

Offline zh

  • QuantumATK Support
  • Supreme QuantumATK Wizard
  • *****
  • Posts: 1141
  • Reputation: 24
    • View Profile
Re: how to calculate Na diffusion coefficient in Si bulk
« Reply #4 on: December 22, 2015, 00:45 »
The ATK already provides a class for calculating the mean-square-displacement from an MD trajectory and an example for the diffusion cofficient:
http://www.quantumwise.com/documents/manuals/latest/ReferenceManual/index.html/ref.meansquaredisplacement.html