Author Topic: Ph.-DOS /transmission eigen values  (Read 16329 times)

0 Members and 1 Guest are viewing this topic.

Offline Dipankar Saha

  • Supreme QuantumATK Wizard
  • *****
  • Posts: 516
  • Country: in
  • Reputation: 5
    • View Profile
Ph.-DOS /transmission eigen values
« on: July 1, 2015, 22:31 »
How can I get some states, at negative energy values....in the Ph.-DOS diagram??
(Even though the "Acoustic  sum rule" is checked.....!!!!)

Best_
Dipankar
« Last Edit: July 6, 2015, 10:57 by Dipankar Saha »

Offline Dipankar Saha

  • Supreme QuantumATK Wizard
  • *****
  • Posts: 516
  • Country: in
  • Reputation: 5
    • View Profile
Re: Ph.-DOS /transmission eigen values
« Reply #1 on: July 6, 2015, 11:00 »

In ATK..... whatever "transmission eigen values"  that we get ......., how they are related to
the Trace[Gr  Sigma1  Sigma2 Ga] .....??

Thanks & Regards__
Dipankar

Offline Julian Schneider

  • QuantumATK Staff
  • QuantumATK Guru
  • *****
  • Posts: 160
  • Country: dk
  • Reputation: 25
    • View Profile
Re: Ph.-DOS /transmission eigen values
« Reply #2 on: July 6, 2015, 12:39 »
Most likely, the atomic positions in your structure are not accurately optimized before calculating the PhononDOS.
If that still does not help, and if you are using classical potentials with a long range (e.g. Coulomb interactions) then you need to increase the maximum interaction range parameter to be larger than the largest cutoff distance of the potential.

Offline Dipankar Saha

  • Supreme QuantumATK Wizard
  • *****
  • Posts: 516
  • Country: in
  • Reputation: 5
    • View Profile
Re: Ph.-DOS /transmission eigen values
« Reply #3 on: July 6, 2015, 14:15 »
Most likely, the atomic positions in your structure are not accurately optimized before calculating the PhononDOS.
If that still does not help, and if you are using classical potentials with a long range (e.g. Coulomb interactions) then you need to increase the maximum interaction range parameter to be larger than the largest cutoff distance of the potential.

Thank you Julian for your reply....../ Optimization May not be the reason in this case.... / However, I didn't get that.... while using classical potetials...... how can I  "increase the maximum interaction range parameter"  (Because, in calculator, all I have .... "potentialSet=     " and  "calculator=   ").....??

Offline Julian Schneider

  • QuantumATK Staff
  • QuantumATK Guru
  • *****
  • Posts: 160
  • Country: dk
  • Reputation: 25
    • View Profile
Re: Ph.-DOS /transmission eigen values
« Reply #4 on: July 6, 2015, 15:04 »
The maximum interaction range is an optional parameter that can be set in the DynamicalMatrix analysis object. When calculating the dynamical matrix, this parameter specifies the distance around a displaced central atom (i) in which the derivative of the forces on atoms (j)  are taken into account to calculate d^2 U / dr_i dr_j .
If this parameter is not specified, the default values are used (which are the covalent distances multiplied by a fuzz factor of 4.0).  In most cases this works well, but when e.g. a classical potential is used, which has a large cutoff radius (e.g. CoulombDSF (http://www.quantumwise.com/documents/manuals/latest/ReferenceManual/index.html/ref.coulombdsf.html) with a default of 9 Angstrom), then the maximum interaction range parameter has to be increased to be larger than this cutoff radius, i.e. larger than 9 Angstrom, to avoid discontinuities.
You can find out if large cutoff values are used by your potential, if you send the script to the Editor with the "Script details" parameter set to "Show defaults". Then all the details of the classical potential are reflected in the script.
Which exact potential are you using?

Apart from that, negative frequencies might also occur if you are sitting exactly at a local maximum (e.g. due to symmetries in your structure). In this case, small random displacements to the atoms followed by another accurate geometry optimization should lead you to a minimum.

Offline Dipankar Saha

  • Supreme QuantumATK Wizard
  • *****
  • Posts: 516
  • Country: in
  • Reputation: 5
    • View Profile
Re: Ph.-DOS /transmission eigen values
« Reply #5 on: July 6, 2015, 19:45 »
Thanks a lot Julian for this detailed answer......   :) / It's actually  the parameterization of SW.... (i.e., StillingerWeber_MoS_2013) ..../ I could not really follow the steps....as you said, like... "if you send the script to the Editor with the "Script details" parameter set to "Show defaults"...then all the details of the classical potential are reflected in the script...." ; but still looking at source paper, I find that they....fitted the phonon spectrum with SW parameters for the two-body and the three body(angle bending as well) potentials.....with a max. cut-off of  4.27 Ang.  / In such a case... should I need to change the default "atomic_displacement=0.01*Angstrom" in dynamical matrix..??

Regards_
Dipankar
« Last Edit: July 6, 2015, 19:48 by Dipankar Saha »

Offline Julian Schneider

  • QuantumATK Staff
  • QuantumATK Guru
  • *****
  • Posts: 160
  • Country: dk
  • Reputation: 25
    • View Profile
Re: Ph.-DOS /transmission eigen values
« Reply #6 on: July 7, 2015, 09:12 »
For the Stillinger-Weber and Tersoff-potentials the default values usually work well. You don't need to specify the maximum interaction range in the DynamicalMatrix parameters.

In your case, I'd rather suspect that there is something odd with the structure. Negative frequencies mean that there are soft modes present in the system, either by atoms sitting in a local saddle point.
You could try to calculate the phonon bandstructure (without checking Accoustic sum rule and Symmetrize) and then the location of the negative frequencies might give you some hints. If the negative frequencies occur at gamma, then you might need to increase the maximum interaction range, if the negative frequencies occur elsewhere, then it is likely that your structure is not in a local minimum.

Offline Dipankar Saha

  • Supreme QuantumATK Wizard
  • *****
  • Posts: 516
  • Country: in
  • Reputation: 5
    • View Profile
Re: Ph.-DOS /transmission eigen values
« Reply #7 on: July 7, 2015, 11:56 »
Thanks a lot.... Julian.....  :)

That's a nice way.....!!!  / Let me....try that....   I will let you know...

Best _
Dipankar

Offline Dipankar Saha

  • Supreme QuantumATK Wizard
  • *****
  • Posts: 516
  • Country: in
  • Reputation: 5
    • View Profile
Re: Ph.-DOS /transmission eigen values
« Reply #8 on: July 7, 2015, 15:56 »
Interestingly.... making the "acoustic_sum_rule=   ",  symmetrize=  " True or False... produces no change in  dispersion...!!!  :o / On the other hand .... "max_interaction_range= " set to default.... or, a very specific range of values ..... could only provide the reqd. Ph-DOS... (otherwise it's arbitrary)!!!

Instead...I went for further opti. of the strcut. ...which somewhat soothes out the density states at negative energy values(see, the attachment)...../ So your guess, that it is due to the struct.  ....seems to be correct....   :)


Offline Dipankar Saha

  • Supreme QuantumATK Wizard
  • *****
  • Posts: 516
  • Country: in
  • Reputation: 5
    • View Profile
Re: Ph.-DOS /transmission eigen values
« Reply #9 on: July 7, 2015, 16:26 »

In ATK..... whatever "transmission eigen values"  that we get ......., how they are related to
the Trace[Gr  Sigma1  Sigma2 Ga] .....??

Besides, it will be of great help ....if you can please comment on this.....!!!  :)

Offline Dipankar Saha

  • Supreme QuantumATK Wizard
  • *****
  • Posts: 516
  • Country: in
  • Reputation: 5
    • View Profile
Re: Ph.-DOS /transmission eigen values
« Reply #10 on: July 10, 2015, 20:03 »
 Dr. Julian Schneider,

For any particular energy 'E', one way to find the Transmission co-eff. is......

T(E) = Summation of [tk tkdagger delta(E-Ek)] for all 'k' states at energy 'E'

Which is again same as.....  Trace[Gr  Sigma1  Sigma2 Ga]   .... :)

On the other hand, we have transmission eigen values of the transmission matrix...which can also be summed up...
for getting the transmission co-efficient...

1) How can one get these tk s (or, as you say transmission amplitudes).....??

Actually, why I'm looking for this...is simply to find an analogy...with the ideal transmisson....either 0 or, 1.

2) Is it incorporated through the delta function....??

Offline Anders Blom

  • QuantumATK Staff
  • Supreme QuantumATK Wizard
  • *****
  • Posts: 5394
  • Country: dk
  • Reputation: 89
    • View Profile
    • QuantumATK at Synopsys

Offline Dipankar Saha

  • Supreme QuantumATK Wizard
  • *****
  • Posts: 516
  • Country: in
  • Reputation: 5
    • View Profile
Re: Ph.-DOS /transmission eigen values
« Reply #12 on: July 12, 2015, 08:37 »
Dr. Anders Blom,

To be more specific,
can I extract the individual channel contributions form the total transmission....??
If yes.... then, which is the parameter...I should look into ??

Thanks in advance...... :)

Best_
Dipankar

Offline Anders Blom

  • QuantumATK Staff
  • Supreme QuantumATK Wizard
  • *****
  • Posts: 5394
  • Country: dk
  • Reputation: 89
    • View Profile
    • QuantumATK at Synopsys
Re: Ph.-DOS /transmission eigen values
« Reply #13 on: July 12, 2015, 16:13 »
That is precisely what the transmission eigenvalues are

Offline Dipankar Saha

  • Supreme QuantumATK Wizard
  • *****
  • Posts: 516
  • Country: in
  • Reputation: 5
    • View Profile
Re: Ph.-DOS /transmission eigen values
« Reply #14 on: July 12, 2015, 19:14 »
Okay....  :) / Thanks a lot...!!!
« Last Edit: July 12, 2015, 19:20 by Dipankar Saha »