Author Topic: mobility calculation  (Read 4460 times)

0 Members and 1 Guest are viewing this topic.

Offline ramkrishna

  • Supreme QuantumATK Wizard
  • *****
  • Posts: 253
  • Country: us
  • Reputation: 5
    • View Profile
mobility calculation
« on: May 8, 2016, 01:38 »
Dear Sir,
Is it possible to calculate electron-phonon coupling and mobility of a system by combing DFT (or SE) along with the classical potential, specially for a large system (i.e. large number of atoms in a unit cell)? We know that the classical potential does not support to calculate electrical properties, on the other hand it is nearly impossible to calculate dynamical matrix and Hamiltonian derivatives by DFT for a large system. So, if that can be possible,  then, I think, it will be easy to handle a large system.

Thanks
Ramkrishna

Offline Jess Wellendorff

  • QuantumATK Staff
  • Supreme QuantumATK Wizard
  • *****
  • Posts: 933
  • Country: dk
  • Reputation: 29
    • View Profile
Re: mobility calculation
« Reply #1 on: May 9, 2016, 10:10 »
Yes, provided that you have a good set of classical potentials, you can use DFT for the electronic parts of the calculations and ATK-Classical for the phonon parts.

Offline ramkrishna

  • Supreme QuantumATK Wizard
  • *****
  • Posts: 253
  • Country: us
  • Reputation: 5
    • View Profile
Re: mobility calculation
« Reply #2 on: May 9, 2016, 17:11 »
Hi Jess,
Thanks for your reply. Suppose I am looking for Silicon nanowire, can you please let me know the way to combine DFT (or Slater-Koster) electronic part and ATK-classical phonon part to calculate electron phonon coupling and mobility in the tutorial http://docs.quantumwise.com/tutorials/mobility.html

Thanks
Ramkrishna

Offline Petr Khomyakov

  • QuantumATK Staff
  • Supreme QuantumATK Wizard
  • *****
  • Posts: 1290
  • Country: dk
  • Reputation: 25
    • View Profile
Re: mobility calculation
« Reply #3 on: May 10, 2016, 00:59 »
For mobility calculation, you may do a classical molecular dynamics (MD) simulation for 'DeviceConfiguration', see http://docs.quantumwise.com/tutorials/md_basics.html, combined with conductance calculations.  In your case, Device Configuration is a silicon nanowire with a given length attached to the left and right electrodes, which can be designed as semi-infinite silicon nanowires to avoid a contact resistance issue as it would be the case for metal electrodes.

The actual purpose of the classical MD simulation is to use its trajectory to calculate the conductance for a set of frames along the trajectory.  Each frame corresponds to a certain atom arrangement in the finite-length nanowire, which is a central region in the Device Configuration described in the previous paragraph.

The mobility or conductivity of an infinite silicon nanowire can then be related to the finite-length nanowire conductance averaged over the trajectory frames selected. Assuming that the nanowire length in the device central region is more than the electron mean free path, the conductivity can be defined as the trajectory-averaged conductance multiplied by the nanowire length. If the nanowire diameter is large enough to neglect a surface effect, the conductance can be divided by the nanowire diameter to get a bulk-like conductivity.  To reach the actual bulk limit, the nanowire diameter should exceed the mean free path indeed.

The mobility can be related to the trajectory-averaged conductivity as discussed in the theory section of the tutorial on 'Mobility', http://docs.quantumwise.com/tutorials/mobility.html.

The theoretical grounds for the described procedure are based on the fact that electron dynamics is much faster than that of ions, i.e., the electrons go through the nanowire faster than ions change their positions, so that the electrons "see" a static lattice of ions at the time scale of electron dynamics. The described approach does not require explicit calculation of the electron-phonon coupling. Note that It accounts for the elastic electron-phonon scattering only. This is perhaps good enough for silicon since it is not a polar semiconductor, i.e., it does not have polar optical phonon modes like in III-V semiconductors.

To avoid direct MD simulations, one may also follow the procedure described in http://journals.aps.org/prb/abstract/10.1103/PhysRevB.91.220405 and applied to metals, but this approach is limited to harmonic approximation, and would also require an additional effort to be interfaced with ATK, unlike the MD simulator that is already a part of ATK.


Offline ramkrishna

  • Supreme QuantumATK Wizard
  • *****
  • Posts: 253
  • Country: us
  • Reputation: 5
    • View Profile
Re: mobility calculation
« Reply #4 on: May 10, 2016, 20:02 »
Hi Petr,
Thanks for a detailed procedure to calculate mobility through MD simulation (elastic limit). I can certainly try this method for my calculations.

Thanks
Ramkrishna