PCDAT: Difference between revisions

From VASP Wiki
No edit summary
No edit summary
 
(15 intermediate revisions by 3 users not shown)
Line 1: Line 1:
The {{TAG|PCDAT}} file contains the pair correlation function. For dynamic simulations ({{TAG|IBRION}}<math>\ge</math>0) an averaged pair correlation is written to the file (see also {{TAG|NBLOCK}}, {{TAG|KBLOCK}}, {{TAG|NPACO}} and {{TAG|APACO}}).
The {{TAG|PCDAT}} file contains the pair correlation function. For dynamic simulations ({{TAG|IBRION}}<math>\ge</math>0) an averaged pair correlation is written to the file (see also {{TAG|NBLOCK}}, {{TAG|KBLOCK}}, {{TAG|NPACO}} and {{TAG|APACO}}).


A sample output of the {{TAG|PCDAT}} file looks as follows:
A sample output of the {{TAG|PCDAT}} file for a system containing two element types looks as follows:
   1  8  1  0  0.8163705E+01  0.1000000E+04
   1  8  1  0  0.8163705E+01  0.1000000E+04
   CAR
   CAR
Line 55: Line 55:
*Line  3: Header of {{TAG|INCAR}} file (the tag {{TAG|SYSTEM}}).
*Line  3: Header of {{TAG|INCAR}} file (the tag {{TAG|SYSTEM}}).
*Line  4: 0, 0, 0 (all fixed output).
*Line  4: 0, 0, 0 (all fixed output).
*Line  5: 1 (fixed output), {{TAG|KBLOCK}}<math>\times</math>{{TAG|NBLOCK}}.
*Line  5: 1 (fixed output), {{TAG|NBLOCK}}.
*Line  6: {{TAG|NPACO}}, {{TAG|NPACO}}, {{TAG|NPACO}}.
*Line  6: {{TAG|NPACO}}, {{TAG|NPACO}}, {{TAG|NPACO}}.
*Line  7: {{TAG|NPACO}}.
*Line  7: {{TAG|NPACO}}.
*Line  8: 10<sup>-10</sup> (fixed output).
*Line  8: 0.1*10<sup>-9</sup> (fixed output).
*Line  9: {{TAG|APACO}}<math>\times</math>10<sup>-10</sup>/{{TAG|NPACO}}.
*Line  9: {{TAG|APACO}}<math>\times</math>10<sup>-10</sup>/{{TAG|NPACO}}.
*Line 10: {{TAG|NSW}}/{{TAG|NBLOCK}}/{{TAG|KBLOCK}}.
*Line 10: {{TAG|NSW}}/{{TAG|NBLOCK}}/{{TAG|KBLOCK}}.
Line 66: Line 66:
*Optional ({{TAG|KBLOCK}}<math>\times</math>{{TAG|NBLOCK}}/{{TAG|NSW}})<math>\times</math>{{TAG|NPACO}}+1 lines: The above is repeated {{TAG|KBLOCK}}<math>\times</math>{{TAG|NBLOCK}}/{{TAG|NSW}} times.
*Optional ({{TAG|KBLOCK}}<math>\times</math>{{TAG|NBLOCK}}/{{TAG|NSW}})<math>\times</math>{{TAG|NPACO}}+1 lines: The above is repeated {{TAG|KBLOCK}}<math>\times</math>{{TAG|NBLOCK}}/{{TAG|NSW}} times.


 
The order of species combinations (columns of the pair correlation function) follows column-wise the lower triangle of the species correlation matrix. That means for 3 species the order is the following:
The order of species combinations (columns of the pair correlation function) follow column-wise the lower triangle of the species correlation matrix. That means for 3 species the order is the following:
  total 1-1  1-2  1-3  2-2  2-3  3-3
  1-1  1-2  1-3  2-2  2-3  3-3
The numbers listed above corresponds to the species as encountered in the {{TAG|POSCAR}}/{{TAG|POTCAR}} file. The first column (total) reports the total pair correlation function.
The numbers are the species with respect to the {{TAG|POSCAR}}/{{TAG|POTCAR}} file.


The PCDAT file contains no abscissa. To obtain the pair correlation functions with the corresponding abscissa the following 'bash/awk' script can be used:
The PCDAT file contains no abscissa. To obtain the pair correlation functions with the corresponding abscissa the following 'bash/awk' script can be used:
<div class="toccolours mw-customtoggle-script">'''Click to show/''pair_correlation_xny.sh''</div>
<div class="toccolours mw-customtoggle-script">'''Click to show/''pair_correlation_xny.sh''</div>
<div class="mw-collapsible mw-collapsed" id="mw-customcollapsible-script">{{:pair_correlation_function_xny_script}}</div>
<div class="mw-collapsible mw-collapsed" id="mw-customcollapsible-script">
<pre>
file=PCDAT
awk <$file >PCDAT.xy '
NR==8 { pcskal=$1}
NR==9 { pcfein=$1}
NR==7 { npaco=$1}
NR>=13 { 
  line=line+1
  if (line==1) s=s+1
  if (line==(npaco+1))  {
    print " "
    line=0
  }
  else  {
    a1[line]=  a1[line] + $1
    a2[line]=  a2[line] + $2
    a3[line]=  a3[line] + $3
    a4[line]=  a4[line] + $4
    print (line-0.5)*pcfein/pcskal,$1,$2, $3, $4, $5
  }
}
END {
print "final sets=", s
for (line=1 ; line<=npaco ; line++)
    print (line-0.5)*pcfein/pcskal,a1[line]/s,a2[line]/s,a3[line]/s,a4[line]/s
}
'
</pre>
</div>


To use this script, in your folder with the {{TAG|PCDAT}} file, please copy the content to ''pair_correlation_xny.sh'' and type following:
To use this script, in your folder with the {{TAG|PCDAT}} file, please copy the content to ''pair_correlation_xny.sh'' and type the following:
  bash pair_correlation_xny.sh  
  bash pair_correlation_xny.sh  


The resulting pair correlation function is written to  
The resulting pair correlation function is written to  
  PCDAT.xy
  PCDAT.xy
== Related tags and articles ==
{{TAG|IBRION}}, {{TAG|MDALGO}}, {{TAG|NBLOCK}}, {{TAG|KBLOCK}}, {{TAG|NSW}}, {{TAG|NPACO}}, {{TAG|APACO}}
----
----


[[Category:Files]][[Category:Output Files]]
[[Category:Files]][[Category:Output files|Output files]][[Category:Ensemble properties|Ensemble properties]][[Category:Symmetry|Symmetry]]

Latest revision as of 12:42, 16 February 2024

The PCDAT file contains the pair correlation function. For dynamic simulations (IBRION0) an averaged pair correlation is written to the file (see also NBLOCK, KBLOCK, NPACO and APACO).

A sample output of the PCDAT file for a system containing two element types looks as follows:

  1   8   1   0  0.8163705E+01  0.1000000E+04
 CAR
 structure name
   0   0   0
   1   1
 350 350 350
 350
  0.1000000E-09
  0.2857143E-11
   1
  0.1000000E-14  0.4027100E-09  0.4027100E-09  0.4027100E-09
  0.2410163E+04  0.2410163E+04
  0.000  0.000  0.000  0.000
  0.000  0.000  0.000  0.000
  0.000  0.000  0.000  0.000
  0.000  0.000  0.000  0.000
  0.000  0.000  0.000  0.000
  0.000  0.000  0.000  0.000
  0.000  0.000  0.000  0.000
  0.000  0.000  0.000  0.000
  0.000  0.000  0.000  0.000
  0.000  0.000  0.000  0.000
  0.000  0.000  0.000  0.000
  0.000  0.000  0.000  0.000
  0.000  0.000  0.000  0.000
  ...
  0.000  0.000  0.000  0.000
  0.165  0.000  0.331  0.000
  0.000  0.000  0.000  0.000
  0.152  0.000  0.304  0.000
  0.293  0.000  0.585  0.000
  0.844  0.000  1.688  0.000
  1.218  0.000  2.436  0.000
  1.173  0.000  2.346  0.000
  1.256  0.000  2.512  0.000
  1.453  0.000  2.906  0.000
  1.168  0.000  2.337  0.000
  1.918  0.000  3.836  0.000
  0.981  0.000  1.962  0.000
  1.580  0.000  3.161  0.000
  0.917  0.000  1.834  0.000
  0.986  0.000  1.972  0.000
  1.528  0.000  3.056  0.000
  1.203  0.000  2.405  0.000
  0.538  0.000  1.076  0.000
  0.869  0.000  1.739  0.000
  0.759  0.000  1.518  0.000 
  ...

Here is the description of each line:

  • Line 1: 1 (fixed output), number of ions, 1 (fixed output), 0 (fixed output), unit cell volume divided by number of atoms, temperature.
  • Line 2: CAR (fixed output).
  • Line 3: Header of INCAR file (the tag SYSTEM).
  • Line 4: 0, 0, 0 (all fixed output).
  • Line 5: 1 (fixed output), NBLOCK.
  • Line 6: NPACO, NPACO, NPACO.
  • Line 7: NPACO.
  • Line 8: 0.1*10-9 (fixed output).
  • Line 9: APACO10-10/NPACO.
  • Line 10: NSW/NBLOCK/KBLOCK.
  • Line 11: POTIM10-15, norm of lattice vector 1 times 10-10, norm of lattice vector 2 times 10-10, norm of lattice vector 3 times 10-10.
  • Line 12-(12+NPACO): Input mean temperature/(NBLOCKKBLOCK), actual mean temperature.
  • Following that the next NPACO lines show the pair correlation function for each species combination.
  • Optional (KBLOCKNBLOCK/NSW)NPACO+1 lines: The above is repeated KBLOCKNBLOCK/NSW times.

The order of species combinations (columns of the pair correlation function) follows column-wise the lower triangle of the species correlation matrix. That means for 3 species the order is the following:

total 1-1  1-2  1-3  2-2  2-3  3-3

The numbers listed above corresponds to the species as encountered in the POSCAR/POTCAR file. The first column (total) reports the total pair correlation function.

The PCDAT file contains no abscissa. To obtain the pair correlation functions with the corresponding abscissa the following 'bash/awk' script can be used:

Click to show/pair_correlation_xny.sh
file=PCDAT
awk <$file >PCDAT.xy '
NR==8 { pcskal=$1}
NR==9 { pcfein=$1}
NR==7 { npaco=$1}
NR>=13 {  
  line=line+1
  if (line==1) s=s+1
  if (line==(npaco+1))  {
     print " "
     line=0
  }
  else  {
     a1[line]=  a1[line] + $1
     a2[line]=  a2[line] + $2
     a3[line]=  a3[line] + $3
     a4[line]=  a4[line] + $4
     print (line-0.5)*pcfein/pcskal,$1,$2, $3, $4, $5
  }
}
END {
 print "final sets=", s
 for (line=1 ; line<=npaco ; line++)
     print (line-0.5)*pcfein/pcskal,a1[line]/s,a2[line]/s,a3[line]/s,a4[line]/s
}
'

To use this script, in your folder with the PCDAT file, please copy the content to pair_correlation_xny.sh and type the following:

bash pair_correlation_xny.sh 

The resulting pair correlation function is written to

PCDAT.xy

Related tags and articles

IBRION, MDALGO, NBLOCK, KBLOCK, NSW, NPACO, APACO