openmx的非平衡格林函数(negf)范例修改:NEGF-8ZGNR-0.3.dat

在OpenMX当中有一个NEGF范例,它得出的transmission图是沿0轴对称的,事实上却是应该是这样的。而我做出的却并不是。所以需要在原来的文件上修改,以逐步达成需要。

先来看左电极:

#
# File Name
#

System.CurrrentDirectory         ./    # default=./
System.Name                     lead-l-8zgnr
level.of.stdout                   1    # default=1 (1-3)
level.of.fileout                  0    # default=1 (0-2)


#
# Definition of Atomic Species
#

Species.Number       2
<Definition.of.Atomic.Species
 C  C5.0-s2p1  C_CA13
 H  H5.0-s2    H_CA13
Definition.of.Atomic.Species>


NEGF.output_hks    true
NEGF.filename.hks  lead-l-8zgnr.hks


#
# Atoms
#

Atoms.Number       36
Atoms.SpeciesAndCoordinates.Unit   Ang # Ang|AU
<Atoms.SpeciesAndCoordinates 
   1  C     0.0000005     0.0361626     0.0012037     2.5000000     1.5000000
   2  C     0.0000005     0.7140385     1.2211879     2.5000000     1.5000000
   3  C     0.0000004     2.1457201     1.2213117     2.5000000     1.5000000
   4  C     0.0000001     2.8501725     0.0013209     2.5000000     1.5000000
   5  C    -0.0000001     4.2700744     0.0012023     2.5000000     1.5000000
   6  C    -0.0000002     4.9766306     1.2212146     2.5000000     1.5000000
   7  C     0.0000001     6.3942050     1.2213979     2.5000000     1.5000000
   8  C     0.0000001     7.1014018     0.0014044     2.5000000     1.5000000
   9  C    -0.0000001     8.5186603     0.0014042     2.5000000     1.5000000
  10  C     0.0000001     9.2258519     1.2213977     2.5000000     1.5000000
  11  C     0.0000004    10.6434033     1.2212143     2.5000000     1.5000000
  12  C     0.0000003    11.3499296     0.0012020     2.5000000     1.5000000
  13  C    -0.0000002    12.7697977     0.0013211     2.5000000     1.5000000
  14  C    -0.0000004    13.4742418     1.2213119     2.5000000     1.5000000
  15  C    -0.0000004    14.9059644     1.2211875     2.5000000     1.5000000
  16  C    -0.0000004    15.5838921     0.0012032     2.5000000     1.5000000
  17  H     0.0000003    -1.0818549     0.0012134     0.5000000     0.5000000
  18  H    -0.0000005    16.7019089     0.0012129     0.5000000     0.5000000
  19  C     0.0000005     0.0361626     2.4412037     2.5000000     1.5000000
  20  C     0.0000005     0.7140385     3.6611879     2.5000000     1.5000000
  21  C     0.0000004     2.1457201     3.6613117     2.5000000     1.5000000
  22  C     0.0000001     2.8501725     2.4413209     2.5000000     1.5000000
  23  C    -0.0000001     4.2700744     2.4412023     2.5000000     1.5000000
  24  C    -0.0000002     4.9766306     3.6612146     2.5000000     1.5000000
  25  C     0.0000001     6.3942050     3.6613979     2.5000000     1.5000000
  26  C     0.0000001     7.1014018     2.4414044     2.5000000     1.5000000
  27  C    -0.0000001     8.5186603     2.4414042     2.5000000     1.5000000
  28  C     0.0000001     9.2258519     3.6613977     2.5000000     1.5000000
  29  C     0.0000004    10.6434033     3.6612143     2.5000000     1.5000000
  30  C     0.0000003    11.3499296     2.4412020     2.5000000     1.5000000
  31  C    -0.0000002    12.7697977     2.4413211     2.5000000     1.5000000
  32  C    -0.0000004    13.4742418     3.6613119     2.5000000     1.5000000
  33  C    -0.0000004    14.9059644     3.6611875     2.5000000     1.5000000
  34  C    -0.0000004    15.5838921     2.4412032     2.5000000     1.5000000
  35  H     0.0000003    -1.0818549     2.4412134     0.5000000     0.5000000
  36  H    -0.0000005    16.7019089     2.4412129     0.5000000     0.5000000
Atoms.SpeciesAndCoordinates>
Atoms.UnitVectors.Unit             Ang # Ang|AU
<Atoms.UnitVectors 
 0.000  0.000  4.88
 10.000 0.000 0.000 
 0.000 27.000 0.000 
Atoms.UnitVectors>


#
# SCF or Electronic System
#


scf.XcType                 LSDA-CA     # LDA|LSDA-CA|LSDA-PW|GGA-PBE
scf.SpinPolarization        on         # On|Off|NC
scf.ElectronicTemperature  300.0       # default=300 (K)
scf.energycutoff           120.0       # default=150 (Ry)
#scf.Ngrid                  140 140 100  # about=200 (Ry)
scf.maxIter                1000         # default=40
scf.EigenvalueSolver       Band     # DC|GDC|Cluster|Band
scf.lapack.dste            dstevx      # dstegr|dstedc|dstevx, default=dstegr
scf.Kgrid                  201 1 1       # means n1 x n2 x n3
scf.Mixing.Type           rmm-diisk    # Simple|Rmm-Diis|Gr-Pulay|Kerker|Rmm-Diisk
scf.Kerker.factor          10.0        # default=1
scf.Init.Mixing.Weight     0.001       # default=0.30
scf.Min.Mixing.Weight      0.001       # default=0.001
scf.Max.Mixing.Weight      0.105       # default=0.40
scf.Mixing.History          40       # default=5
scf.Mixing.StartPulay       10       # default=6
scf.Mixing.EveryPulay        1          # default=6
scf.criterion             1.0e-9       # default=1.0e-6 (Hartree)

#
# MD or Geometry Optimization
#

MD.Type               Nomd    # Opt|EF|BFGS|RF|DIIS
MD.Opt.DIIS.History   18    # default=3
MD.Opt.StartDIIS      6     # default=5
MD.Opt.EveryDIIS      5     # default=10
MD.maxIter            2000   #
MD.Opt.criterion    1.0e-4  # default=1.0e-4 (a.u.)


#
# Band dispersion
#

Band.dispersion              off      # on|off, default=off
<Band.KPath.UnitCell
 0.000 0.000  4.88
 10.000 0.000 0.000 
 0.000 27.000 0.000 
Band.KPath.UnitCell>

Band.Nkpath                1
<Band.kpath                
   601  0.0 0.0 0.0   0.5 0.0 0.0   G X
Band.kpath>


#
# DOS and PDOS
#

Dos.fileout                  off      # on|off, default=off
Dos.Erange              -20.0  20.0    # default = -20 20
Dos.Kgrid                601 1 1      # default = Kgrid1 Kgrid2 Kgrid3

DosGauss.fileout      off
DosGauss.Num.Mesh    4000
DosGauss.Width      0.01

MO.fileout                   off       # on|off
num.HOMOs                     1        # default=1
num.LUMOs                     1        # default=1
MO.Nkpoint                    1        # default=1
<MO.kpoint
  0.0  0.0  0.0
MO.kpoint>

中间标红的一块区域为Central区域,左边有四个氢,右边也有四个氢。

由这个图出发,可以将现有的左右电极都化作中心,将现有的左右电极都复制分别向外侧移动,形成新的研究结构。为此需要确定坐标的轴与原子的实际关系。

通过简单分析,发现第一列全是0,这代表z轴。中心部分的第三列全是4.8812037,意味着第三列是x轴。剩下的第二列是y轴。

一、中心区域重构

经过一番研究。发现如果可以重构原子位置,进行一下剪接,可以避免原子数目急剧扩大的风险。

#
# File Name
#

System.CurrrentDirectory         ./    # default=./
System.Name                    negf-8zgnr-0.3
level.of.stdout                   1    # default=1 (1-3)
level.of.fileout                  1    # default=1 (0-2)


NEGF.filename.hks.l     lead-l-8zgnr.hks
NEGF.filename.hks.r     lead-r-8zgnr.hks

NEGF.SCF.Iter.Band          20
NEGF.Num.Poles             100       # defalut=150
NEGF.scf.Kgrid             1 1       # defalut=1 1

NEGF.bias.voltage         -0.3       # default=0.0 (eV)
NEGF.bias.neq.im.energy    0.01      # default=0.01 (eV)
NEGF.bias.neq.energy.step  0.02      # default=0.02 (eV)

Dos.fileout                 off       # on|off, default=off
NEGF.Dos.energyrange      -20.0 20.0 1.0e-4  #default=-10.0 10.0 5.0e-3 (eV)
NEGF.Dos.energy.div        4000       # default=200
NEGF.Dos.Kgrid             1 1       # default=1 1

#
# Transmission & Eigen Channel Analysis
#
 
NEGF.tran.Analysis          on       # default on
NEGF.tran.energydiv        800       # default=200
NEGF.tran.energyrange   -4 4 0.001   # default=-10.0 10.0 5.0e-3 (eV)
NEGF.tran.Kgrid            1 1       # default= 1 1

NEGF.tran.CurrentDensity    on       # default on

NEGF.tran.Channel           on       # default on
NEGF.Channel.Nkpoint         1       # default 1
<NEGF.Channel.kpoint
0.0  0.0
NEGF.Channel.kpoint>
# default 0.0  0.0
NEGF.Channel.Nenergy         1       # default 1 
<NEGF.Channel.energy
0.0
NEGF.Channel.energy>
# default 0.0
NEGF.Channel.Num            5        # defualt=5(for collinear), 10(for Non-collinear)


#
# Definition of Atomic Species
#

Species.Number       2
<Definition.of.Atomic.Species
 C  C5.0-s2p1  C_CA13
 H  H5.0-s2    H_CA13
Definition.of.Atomic.Species>


#
# Atoms
#

Atoms.SpeciesAndCoordinates.Unit   Ang # Ang|AU


Atoms.Number        10
<Atoms.SpeciesAndCoordinates
   1  C     0.0000005     0.0361626     4.8812037     2.1000000     1.9000000
   2  C     0.0000001     2.8501725     4.8813209     2.1000000     1.9000000
   3  C    -0.0000001     4.2700744     4.8812023     2.1000000     1.9000000
   4  C     0.0000001     7.1014018     4.8814044     2.1000000     1.9000000
   5  C    -0.0000001     8.5186603     4.8814042     2.1000000     1.9000000
   6  C     0.0000003    11.3499296     4.8812020     2.1000000     1.9000000
   7  C    -0.0000002    12.7697977     4.8813211     2.1000000     1.9000000
   8  C    -0.0000004    15.5838921     4.8812032     2.1000000     1.9000000
   9  H     0.0000003    -1.0818549     4.8812134     0.5000000     0.5000000
  10  H    -0.0000005    16.7019089     4.8812129     0.5000000     0.5000000
Atoms.SpeciesAndCoordinates>

Atoms.UnitVectors.Unit             Ang # Ang|AU

#<Atoms.UnitVectors
# 0.000 0.000  4.88
# 10.000 0.000 0.000
# 0.000 24.000 0.000
#Atoms.UnitVectors>


#
# Lead-Left
#

LeftLeadAtoms.Number         36
<LeftLeadAtoms.SpeciesAndCoordinates           # Unit=Ang.
   1  C     0.0000005     0.0361626     0.0012037     2.1000000     1.9000000
   2  C     0.0000005     0.7140385     1.2211879     2.1000000     1.9000000
   3  C     0.0000004     2.1457201     1.2213117     2.1000000     1.9000000
   4  C     0.0000001     2.8501725     0.0013209     2.1000000     1.9000000
   5  C    -0.0000001     4.2700744     0.0012023     2.1000000     1.9000000
   6  C    -0.0000002     4.9766306     1.2212146     2.1000000     1.9000000
   7  C     0.0000001     6.3942050     1.2213979     2.1000000     1.9000000
   8  C     0.0000001     7.1014018     0.0014044     2.1000000     1.9000000
   9  C    -0.0000001     8.5186603     0.0014042     2.1000000     1.9000000
  10  C     0.0000001     9.2258519     1.2213977     2.1000000     1.9000000
  11  C     0.0000004    10.6434033     1.2212143     2.1000000     1.9000000
  12  C     0.0000003    11.3499296     0.0012020     2.1000000     1.9000000
  13  C    -0.0000002    12.7697977     0.0013211     2.1000000     1.9000000
  14  C    -0.0000004    13.4742418     1.2213119     2.1000000     1.9000000
  15  C    -0.0000004    14.9059644     1.2211875     2.1000000     1.9000000
  16  C    -0.0000004    15.5838921     0.0012032     2.1000000     1.9000000
  17  H     0.0000003    -1.0818549     0.0012134     0.5000000     0.5000000
  18  H    -0.0000005    16.7019089     0.0012129     0.5000000     0.5000000
  19  C     0.0000005     0.0361626     2.4412037     2.1000000     1.9000000
  20  C     0.0000005     0.7140385     3.6611879     2.1000000     1.9000000
  21  C     0.0000004     2.1457201     3.6613117     2.1000000     1.9000000
  22  C     0.0000001     2.8501725     2.4413209     2.1000000     1.9000000
  23  C    -0.0000001     4.2700744     2.4412023     2.1000000     1.9000000
  24  C    -0.0000002     4.9766306     3.6612146     2.1000000     1.9000000
  25  C     0.0000001     6.3942050     3.6613979     2.1000000     1.9000000
  26  C     0.0000001     7.1014018     2.4414044     2.1000000     1.9000000
  27  C    -0.0000001     8.5186603     2.4414042     2.1000000     1.9000000
  28  C     0.0000001     9.2258519     3.6613977     2.1000000     1.9000000
  29  C     0.0000004    10.6434033     3.6612143     2.1000000     1.9000000
  30  C     0.0000003    11.3499296     2.4412020     2.1000000     1.9000000
  31  C    -0.0000002    12.7697977     2.4413211     2.1000000     1.9000000
  32  C    -0.0000004    13.4742418     3.6613119     2.1000000     1.9000000
  33  C    -0.0000004    14.9059644     3.6611875     2.1000000     1.9000000
  34  C    -0.0000004    15.5838921     2.4412032     2.1000000     1.9000000
  35  H     0.0000003    -1.0818549     2.4412134     0.5000000     0.5000000
  36  H    -0.0000005    16.7019089     2.4412129     0.5000000     0.5000000
LeftLeadAtoms.SpeciesAndCoordinates>


#
# Lead-Right
#


RightLeadAtoms.Number        36

<RightLeadAtoms.SpeciesAndCoordinates           # Unit=Ang.
   1  C     0.0000005     0.7140385     6.1011879     1.9000000     2.1000000
   2  C     0.0000004     2.1457201     6.1013117     1.9000000     2.1000000
   3  C    -0.0000002     4.9766306     6.1012146     1.9000000     2.1000000
   4  C     0.0000001     6.3942050     6.1013979     1.9000000     2.1000000
   5  C     0.0000001     9.2258519     6.1013977     1.9000000     2.1000000
   6  C     0.0000004    10.6434033     6.1012143     1.9000000     2.1000000
   7  C    -0.0000004    13.4742418     6.1013119     1.9000000     2.1000000
   8  C    -0.0000004    14.9059644     6.1011875     1.9000000     2.1000000
   9  C     0.0000005     0.0361626     7.3212037     1.9000000     2.1000000
  10  C     0.0000005     0.7140385     8.5411879     1.9000000     2.1000000
  11  C     0.0000004     2.1457201     8.5413117     1.9000000     2.1000000
  12  C     0.0000001     2.8501725     7.3213209     1.9000000     2.1000000
  13  C    -0.0000001     4.2700744     7.3212023     1.9000000     2.1000000
  14  C    -0.0000002     4.9766306     8.5412146     1.9000000     2.1000000
  15  C     0.0000001     6.3942050     8.5413979     1.9000000     2.1000000
  16  C     0.0000001     7.1014018     7.3214044     1.9000000     2.1000000
  17  C    -0.0000001     8.5186603     7.3214042     1.9000000     2.1000000
  18  C     0.0000001     9.2258519     8.5413977     1.9000000     2.1000000
  19  C     0.0000004    10.6434033     8.5412143     1.9000000     2.1000000
  20  C     0.0000003    11.3499296     7.3212020     1.9000000     2.1000000
  21  C    -0.0000002    12.7697977     7.3213211     1.9000000     2.1000000
  22  C    -0.0000004    13.4742418     8.5413119     1.9000000     2.1000000
  23  C    -0.0000004    14.9059644     8.5411875     1.9000000     2.1000000
  24  C    -0.0000004    15.5838921     7.3212032     1.9000000     2.1000000
  25  H     0.0000003    -1.0818549     7.3212134     0.5000000     0.5000000
  26  H    -0.0000005    16.7019089     7.3212129     0.5000000     0.5000000
  27  C     0.0000005     0.0361626     9.7612037     1.9000000     2.1000000
  28  C     0.0000001     2.8501725     9.7613209     1.9000000     2.1000000
  29  C    -0.0000001     4.2700744     9.7612023     1.9000000     2.1000000
  30  C     0.0000001     7.1014018     9.7614044     1.9000000     2.1000000
  31  C    -0.0000001     8.5186603     9.7614042     1.9000000     2.1000000
  32  C     0.0000003    11.3499296     9.7612020     1.9000000     2.1000000
  33  C    -0.0000002    12.7697977     9.7613211     1.9000000     2.1000000
  34  C    -0.0000004    15.5838921     9.7612032     1.9000000     2.1000000
  35  H     0.0000003    -1.0818549     9.7612134     0.5000000     0.5000000
  36  H    -0.0000005    16.7019089     9.7612129     0.5000000     0.5000000
RightLeadAtoms.SpeciesAndCoordinates>

#
# SCF or Electronic System
#


scf.XcType                LSDA-CA      # LDA|LSDA-CA|LSDA-PW|GGA-PBE
scf.SpinPolarization        on         # On|Off|NC
scf.ElectronicTemperature  300.0       # default=300 (K)
scf.energycutoff           120.0       # default=150 (Ry)
#scf.Ngrid                  140 140 100 # about=200 (Ry)
scf.maxIter                2000         # default=40
scf.EigenvalueSolver       NEGF        # DC|GDC|Cluster|Band
scf.lapack.dste            dstevx      # dstegr|dstedc|dstevx, default=dstegr
scf.Kgrid                  1 1 1       # means n1 x n2 x n3
scf.Mixing.Type           rmm-diisk    # Simple|Rmm-Diis|Gr-Pulay|Kerker|Rmm-Diisk
scf.Kerker.factor          20.0        # default=1
scf.Init.Mixing.Weight     0.1         # default=0.30
scf.Min.Mixing.Weight      0.04        # default=0.001
scf.Max.Mixing.Weight      0.4         # default=0.40
scf.Mixing.History          60         # default=5
scf.Mixing.StartPulay       20         # default=6
scf.Mixing.EveryPulay       1          # default=6
scf.criterion             1.0e-7       # default=1.0e-6 (Hartree)


#
# MD or Geometry Optimization
#

MD.Type               Nomd    # Opt|EF|BFGS|RF|DIIS
MD.Opt.DIIS.History   10    # default=3
MD.Opt.StartDIIS      6     # default=5
MD.Opt.EveryDIIS      5     # default=10
MD.maxIter            200   #
MD.Opt.criterion    1.0e-4  # default=1.0e-4 (a.u.)

方案如下图所示:

将黑色方框所标的设置为中心,并且补上三个氢原子或者直接删掉下面的三个氢原子。分别新增6个碳原子到左、右电极处,并且分别在左右电极上方补上一个氢原子。

为完成目标先进行第一步,分别将左右电极的12个原子放进中心区域。

在左电极当中,第三列坐标为3.66的8个碳原子选取第二列坐标较小的6个。这六个碳原子在左电极当中的编号为20、21、24、25、28、29。

在左电极当中,第三列坐标为2.44的8个碳原子选取第二列坐标较小的6个。这六个碳原子在左电极当中的编号为19、22、23、26、27、30。

在右电极当中,第三列坐标为6.1的8个碳原子选取第二列坐标较小的6个。这六个碳原子在左电极当中的编号为1~6。

在右电极当中,第三列坐标为7.32的8个碳原子选取第二列坐标较小的6个。这六个碳原子在左电极当中的编号为9、12、13、16、17、20。

将左右电极的碳原子放入中心区域后,还需要在中心区域增加氢原子。中心区域要保留的原有氢原子编号为距离最近的碳原子9,坐标为

9  H     0.0000003    -1.0818549     4.8812134     0.5000000     0.5000000

与之相邻的碳原子编号为1,坐标为

   1  C     0.0000005     0.0361626     4.8812037     2.1000000     1.9000000

二者第二列坐标差为0.0361626-(-1.0818549)=1.1180。这也就是C-H键长。接下来需要在原中心区的8个碳原子去掉第二列坐标最高的两个编号为7、8的这两个碳原子。再去掉编号为10的氢原子。再在6号碳原子坐标为:

   6  C     0.0000003    11.3499296     4.8812020     2.1000000     1.9000000

头上构建一个第二列坐标与之距离1.118的氢原子。计算11.3499296+1.118=12.4679

新的氢原子坐标为:

   7  H     0.0000003    12.4679        4.8812020     0.5000000     0.5000000

接下来的工作需要按照上面计算出来的新氢原子坐标,在新的中心区域中补齐剩下的4个氢原子,它们的第二列坐标分别为2.44和7.32。

4个氢原子里面有两个是现成的,分别是原左电极的35号原子:

  35  H     0.0000003    -1.0818549     2.4412134     0.5000000     0.5000000

和右电极的25号氢原子:

  25  H     0.0000003    -1.0818549     7.3212134     0.5000000     0.5000000

我们可以直接把上面两个氢原子的第二列修改为需要的12.4679,即可完成剩下两个原子的构建。

  36  H     0.0000003    12.4679        2.4412134     0.5000000     0.5000000
  26  H     0.0000003    12.4679        7.3212134     0.5000000     0.5000000

现在已经完成了,中心区域的重构,标号未修整的原子列表如下:

   1  C     0.0000005     0.0361626     4.8812037     2.1000000     1.9000000
   2  C     0.0000001     2.8501725     4.8813209     2.1000000     1.9000000
   3  C    -0.0000001     4.2700744     4.8812023     2.1000000     1.9000000
   4  C     0.0000001     7.1014018     4.8814044     2.1000000     1.9000000
   5  C    -0.0000001     8.5186603     4.8814042     2.1000000     1.9000000
   6  C     0.0000003    11.3499296     4.8812020     2.1000000     1.9000000
   7  H     0.0000003    12.4679        4.8812020     0.5000000     0.5000000
   9  H     0.0000003    -1.0818549     4.8812134     0.5000000     0.5000000
  20  C     0.0000005     0.7140385     3.6611879     2.1000000     1.9000000
  21  C     0.0000004     2.1457201     3.6613117     2.1000000     1.9000000
  24  C    -0.0000002     4.9766306     3.6612146     2.1000000     1.9000000
  25  C     0.0000001     6.3942050     3.6613979     2.1000000     1.9000000
  28  C     0.0000001     9.2258519     3.6613977     2.1000000     1.9000000
  29  C     0.0000004    10.6434033     3.6612143     2.1000000     1.9000000
  19  C     0.0000005     0.0361626     2.4412037     2.1000000     1.9000000
  22  C     0.0000001     2.8501725     2.4413209     2.1000000     1.9000000
  23  C    -0.0000001     4.2700744     2.4412023     2.1000000     1.9000000
  26  C     0.0000001     7.1014018     2.42414044     2.1000000     1.9000000
  27  C    -0.0000001     8.5186603     2.4414042     2.1000000     1.9000000
  30  C     0.0000003    11.3499296     2.4412020     2.1000000     1.9000000
  1  C     0.0000005     0.7140385     6.1011879     1.9000000     2.1000000
  2  C     0.0000004     2.1457201     6.1013117     1.9000000     2.1000000
  3  C    -0.0000002     4.9766306     6.1012146     1.9000000     2.1000000
  4  C     0.0000001     6.3942050     6.1013979     1.9000000     2.1000000
  5  C     0.0000001     9.2258519     6.1013977     1.9000000     2.1000000
  6  C     0.0000004    10.6434033     6.1012143     1.9000000     2.1000000
  9  C     0.0000005     0.0361626     7.3212037     1.9000000     2.1000000
  12  C     0.0000001     2.8501725     7.3213209     1.9000000     2.1000000
  13  C    -0.0000001     4.2700744     7.3212023     1.9000000     2.1000000
  16  C     0.0000001     7.1014018     7.3214044     1.9000000     2.1000000
  17  C    -0.0000001     8.5186603     7.3214042     1.9000000     2.1000000
  20  C     0.0000003    11.3499296     7.3212020     1.9000000     2.1000000
  35  H     0.0000003    -1.0818549     2.4412134     0.5000000     0.5000000
  25  H     0.0000003    -1.0818549     7.3212134     0.5000000     0.5000000
  36  H     0.0000003    12.4679        2.4412134     0.5000000     0.5000000
  26  H     0.0000003    12.4679        7.3212134     0.5000000     0.5000000

经过修改后的原子列表如下:

Atoms.Number        36
<Atoms.SpeciesAndCoordinates
   1  C     0.0000005     0.0361626     4.8812037     2.1000000     1.9000000
   2  C     0.0000001     2.8501725     4.8813209     2.1000000     1.9000000
   3  C    -0.0000001     4.2700744     4.8812023     2.1000000     1.9000000
   4  C     0.0000001     7.1014018     4.8814044     2.1000000     1.9000000
   5  C    -0.0000001     8.5186603     4.8814042     2.1000000     1.9000000
   6  C     0.0000003    11.3499296     4.8812020     2.1000000     1.9000000
   7  H     0.0000003    12.4679        4.8812020     0.5000000     0.5000000
   8  H     0.0000003    -1.0818549     4.8812134     0.5000000     0.5000000
   9  C     0.0000005     0.7140385     3.6611879     2.1000000     1.9000000
  10  C     0.0000004     2.1457201     3.6613117     2.1000000     1.9000000
  11  C    -0.0000002     4.9766306     3.6612146     2.1000000     1.9000000
  12  C     0.0000001     6.3942050     3.6613979     2.1000000     1.9000000
  13  C     0.0000001     9.2258519     3.6613977     2.1000000     1.9000000
  14  C     0.0000004    10.6434033     3.6612143     2.1000000     1.9000000
  15  C     0.0000005     0.0361626     2.4412037     2.1000000     1.9000000
  16  C     0.0000001     2.8501725     2.4413209     2.1000000     1.9000000
  17  C    -0.0000001     4.2700744     2.4412023     2.1000000     1.9000000
  18  C     0.0000001     7.1014018     2.42414044     2.1000000     1.9000000
  19  C    -0.0000001     8.5186603     2.4414042     2.1000000     1.9000000
  20  C     0.0000003    11.3499296     2.4412020     2.1000000     1.9000000
  21  C     0.0000005     0.7140385     6.1011879     1.9000000     2.1000000
  22  C     0.0000004     2.1457201     6.1013117     1.9000000     2.1000000
  23  C    -0.0000002     4.9766306     6.1012146     1.9000000     2.1000000
  24  C     0.0000001     6.3942050     6.1013979     1.9000000     2.1000000
  25  C     0.0000001     9.2258519     6.1013977     1.9000000     2.1000000
  26  C     0.0000004    10.6434033     6.1012143     1.9000000     2.1000000
  27  C     0.0000005     0.0361626     7.3212037     1.9000000     2.1000000
  28  C     0.0000001     2.8501725     7.3213209     1.9000000     2.1000000
  29  C    -0.0000001     4.2700744     7.3212023     1.9000000     2.1000000
  30  C     0.0000001     7.1014018     7.3214044     1.9000000     2.1000000
  31  C    -0.0000001     8.5186603     7.3214042     1.9000000     2.1000000
  32  C     0.0000003    11.3499296     7.3212020     1.9000000     2.1000000
  33  H     0.0000003    -1.0818549     2.4412134     0.5000000     0.5000000
  34  H     0.0000003    -1.0818549     7.3212134     0.5000000     0.5000000
  35  H     0.0000003    12.4679        2.4412134     0.5000000     0.5000000
  36  H     0.0000003    12.4679        7.3212134     0.5000000     0.5000000
Atoms.SpeciesAndCoordinates>

一共36个原子,与我们设计的目标一致。

二、左电极

首先需要将橙色框内的六个碳原子复制到浅蓝色的位置,再删除原有的上方的6个碳原子和一个最靠近中心区的上方氢原子(右上角),并且把最上方的氢原子向下移动到12.4679。

在左电极原子列表当中,橙色框中的碳原子的第三列坐标为1.22。我们需要保留其中第二列坐标中较小的6个,被删除的两个碳原子编号为14号和15号。

  14  C    -0.0000004    13.4742418     1.2213119     2.1000000     1.9000000
  15  C    -0.0000004    14.9059644     1.2211875     2.1000000     1.9000000

接下来需要做一道算术题目,石墨烯C-C键长为1.42埃,计算出需要平移的距离为2.4595埃。用1.2211879-2.4595=-1.2383。将编号为2、3、6、7、10、11这六个碳原子分别复制,并修改第三列为-1.2383。

将原18号氢原子的第二列坐标由16.7019089修改为12.4679。

  18  H    -0.0000005    12.4679        0.0012129     0.5000000     0.5000000

要被删除的的6个碳原子和一个氢原子编号分别为:13、16、31、32、33 、34、36。

  13  C    -0.0000002    12.7697977     0.0013211     2.1000000     1.9000000
  16  C    -0.0000004    15.5838921     0.0012032     2.1000000     1.9000000
  31  C    -0.0000002    12.7697977     2.4413211     2.1000000     1.9000000
  32  C    -0.0000004    13.4742418     3.6613119     2.1000000     1.9000000
  33  C    -0.0000004    14.9059644     3.6611875     2.1000000     1.9000000
  34  C    -0.0000004    15.5838921     2.4412032     2.1000000     1.9000000
  36  H    -0.0000005    16.7019089     2.4412129     0.5000000     0.5000000

完成后,未修改编号的左电极变为:

   1  C     0.0000005     0.0361626     0.0012037     2.1000000     1.9000000
   2  C     0.0000005     0.7140385     1.2211879     2.1000000     1.9000000
   3  C     0.0000004     2.1457201     1.2213117     2.1000000     1.9000000
   4  C     0.0000001     2.8501725     0.0013209     2.1000000     1.9000000
   5  C    -0.0000001     4.2700744     0.0012023     2.1000000     1.9000000
   6  C    -0.0000002     4.9766306     1.2212146     2.1000000     1.9000000
   7  C     0.0000001     6.3942050     1.2213979     2.1000000     1.9000000
   8  C     0.0000001     7.1014018     0.0014044     2.1000000     1.9000000
   9  C    -0.0000001     8.5186603     0.0014042     2.1000000     1.9000000
  10  C     0.0000001     9.2258519     1.2213977     2.1000000     1.9000000
  11  C     0.0000004    10.6434033     1.2212143     2.1000000     1.9000000
  12  C     0.0000003    11.3499296     0.0012020     2.1000000     1.9000000
  17  H     0.0000003    -1.0818549     0.0012134     0.5000000     0.5000000
  18  H    -0.0000005    12.4679        0.0012129     0.5000000     0.5000000
  2  C     0.0000005     0.7140385     -1.2383     2.1000000     1.9000000
  3  C     0.0000004     2.1457201     -1.2383     2.1000000     1.9000000
  6  C    -0.0000002     4.9766306     -1.2383     2.1000000     1.9000000
  7  C     0.0000001     6.3942050     -1.2383     2.1000000     1.9000000
  10  C     0.0000001     9.2258519     -1.2383     2.1000000     1.9000000
  11  C     0.0000004    10.6434033     -1.2383     2.1000000     1.9000000

进行修改编号,得到总共20个原子的列表:

LeftLeadAtoms.Number         20
<LeftLeadAtoms.SpeciesAndCoordinates           # Unit=Ang.
   1  C     0.0000005     0.0361626     0.0012037     2.1000000     1.9000000
   2  C     0.0000005     0.7140385     1.2211879     2.1000000     1.9000000
   3  C     0.0000004     2.1457201     1.2213117     2.1000000     1.9000000
   4  C     0.0000001     2.8501725     0.0013209     2.1000000     1.9000000
   5  C    -0.0000001     4.2700744     0.0012023     2.1000000     1.9000000
   6  C    -0.0000002     4.9766306     1.2212146     2.1000000     1.9000000
   7  C     0.0000001     6.3942050     1.2213979     2.1000000     1.9000000
   8  C     0.0000001     7.1014018     0.0014044     2.1000000     1.9000000
   9  C    -0.0000001     8.5186603     0.0014042     2.1000000     1.9000000
  10  C     0.0000001     9.2258519     1.2213977     2.1000000     1.9000000
  11  C     0.0000004    10.6434033     1.2212143     2.1000000     1.9000000
  12  C     0.0000003    11.3499296     0.0012020     2.1000000     1.9000000
  13  H     0.0000003    -1.0818549     0.0012134     0.5000000     0.5000000
  14  H    -0.0000005    12.4679        0.0012129     0.5000000     0.5000000
  15  C     0.0000005     0.7140385     -1.2383     2.1000000     1.9000000
  16  C     0.0000004     2.1457201     -1.2383     2.1000000     1.9000000
  17 C    -0.0000002     4.9766306     -1.2383     2.1000000     1.9000000
  18  C     0.0000001     6.3942050     -1.2383     2.1000000     1.9000000
  19  C     0.0000001     9.2258519     -1.2383     2.1000000     1.9000000
  20  C     0.0000004    10.6434033     -1.2383     2.1000000     1.9000000
LeftLeadAtoms.SpeciesAndCoordinates>

三、右电极

按照左电极的制作方法制作右电极。

应该删除的8个碳原子为:7、8、21、22、23、24、33、34。要删除的氢原子编号为26。它们的坐标如下:

  7  C    -0.0000004    13.4742418     6.1013119     1.9000000     2.1000000
  8  C    -0.0000004    14.9059644     6.1011875     1.9000000     2.1000000
  21  C    -0.0000002    12.7697977     7.3213211     1.9000000     2.1000000
  22  C    -0.0000004    13.4742418     8.5413119     1.9000000     2.1000000
  23  C    -0.0000004    14.9059644     8.5411875     1.9000000     2.1000000
  24  C    -0.0000004    15.5838921     7.3212032     1.9000000     2.1000000
  33  C    -0.0000002    12.7697977     9.7613211     1.9000000     2.1000000
  34  C    -0.0000004    15.5838921     9.7612032     1.9000000     2.1000000
  26  H    -0.0000005    16.7019089     7.3212129     0.5000000     0.5000000

将编号为26的氢原子第二列的16.7019089修改为12.4679。

将6个碳原子复制,并将第三列坐标由8.5411879修改为8.5411879+2.4595=11.0007。这六个碳原子原来的编号为10、11、14、15、18、19。

修改编号之前,原子列表如下:

  10  C     0.0000005     0.7140385     8.5411879     1.9000000     2.1000000
  11  C     0.0000004     2.1457201     8.5413117     1.9000000     2.1000000
  14  C    -0.0000002     4.9766306     8.5412146     1.9000000     2.1000000
  15  C     0.0000001     6.3942050     8.5413979     1.9000000     2.1000000
  18  C     0.0000001     9.2258519     8.5413977     1.9000000     2.1000000
  19  C     0.0000004    10.6434033     8.5412143     1.9000000     2.1000000
  27  C     0.0000005     0.0361626     9.7612037     1.9000000     2.1000000
  28  C     0.0000001     2.8501725     9.7613209     1.9000000     2.1000000
  29  C    -0.0000001     4.2700744     9.7612023     1.9000000     2.1000000
  30  C     0.0000001     7.1014018     9.7614044     1.9000000     2.1000000
  31  C    -0.0000001     8.5186603     9.7614042     1.9000000     2.1000000
  32  C     0.0000003    11.3499296     9.7612020     1.9000000     2.1000000
  35  H     0.0000003    -1.0818549     9.7612134     0.5000000     0.5000000
  36  H    -0.0000005    12.4679        9.7612129     0.5000000     0.5000000
  10  C     0.0000005     0.7140385     11.0007     1.9000000     2.1000000
  11  C     0.0000004     2.1457201     11.0007     1.9000000     2.1000000
  14  C    -0.0000002     4.9766306     11.0007     1.9000000     2.1000000
  15  C     0.0000001     6.3942050     11.0007     1.9000000     2.1000000
  18  C     0.0000001     9.2258519     11.0007     1.9000000     2.1000000
  19  C     0.0000004    10.6434033     11.0007     1.9000000     2.1000000

修改编号之后,原子列表如下,总共20个原子:

  1  C     0.0000005     0.7140385     8.5411879     1.9000000     2.1000000
  2  C     0.0000004     2.1457201     8.5413117     1.9000000     2.1000000
  3  C    -0.0000002     4.9766306     8.5412146     1.9000000     2.1000000
  4  C     0.0000001     6.3942050     8.5413979     1.9000000     2.1000000
  5  C     0.0000001     9.2258519     8.5413977     1.9000000     2.1000000
  6  C     0.0000004    10.6434033     8.5412143     1.9000000     2.1000000
  7  C     0.0000005     0.0361626     9.7612037     1.9000000     2.1000000
  8  C     0.0000001     2.8501725     9.7613209     1.9000000     2.1000000
  9  C    -0.0000001     4.2700744     9.7612023     1.9000000     2.1000000
  10  C     0.0000001     7.1014018     9.7614044     1.9000000     2.1000000
  11  C    -0.0000001     8.5186603     9.7614042     1.9000000     2.1000000
  12  C     0.0000003    11.3499296     9.7612020     1.9000000     2.1000000
  13  H     0.0000003    -1.0818549     9.7612134     0.5000000     0.5000000
  14  H    -0.0000005    12.4679        9.7612129     0.5000000     0.5000000
  15  C     0.0000005     0.7140385     11.0007     1.9000000     2.1000000
  16  C     0.0000004     2.1457201     11.0007     1.9000000     2.1000000
  17  C    -0.0000002     4.9766306     11.0007     1.9000000     2.1000000
  18  C     0.0000001     6.3942050     11.0007     1.9000000     2.1000000
  19  C     0.0000001     9.2258519     11.0007     1.9000000     2.1000000
  20  C     0.0000004    10.6434033     11.0007     1.9000000     2.1000000

Atoms.UnitVectors

在原来的左电极和右电极文件中的晶格矢量为

Atoms.UnitVectors.Unit             Ang # Ang|AU
<Atoms.UnitVectors 
 0.000  0.000  4.88
 10.000 0.000 0.000 
 0.000 27.000 0.000 
Atoms.UnitVectors>

在原NEGF-8ZGNR-0.3.dat文件中的晶格矢量为

#<Atoms.UnitVectors
# 0.000 0.000  4.88
# 10.000 0.000 0.000
# 0.000 24.000 0.000
#Atoms.UnitVectors>

看来前两个一样,最后一个影响不大。参考Openmx用户手册第139页

这段话的意思是说,不要管那么多,像上面那样写,就能把参考方向改成沿原子坐标第三列,而且能保持其他两列的意义不变。具体为什么,我也不知道。

现在沿着第三列的方向模型整体变宽了,原来的半宽恰好是4.88,而现在的半宽是6.1,所以,将Atoms.UnitVectors改成

<Atoms.UnitVectors 
 0.000 0.000  6.1
 10.000 0.000 0.000 
 0.000 24.000 0.000 
Atoms.UnitVectors>

四、左右电极文件

上面仅仅完成了NEGF-8ZGNR-0.3.dat的修改,要想开始仿真还要修改Lead-L-8ZGNR.dat和Lead-R-8ZGNR.dat文件。

修改的地方就是原子个数,原子坐标。

新的左电极文件如下。

#
# File Name
#

System.CurrrentDirectory         ./    # default=./
System.Name                     lead-l-8zgnr
level.of.stdout                   1    # default=1 (1-3)
level.of.fileout                  0    # default=1 (0-2)


#
# Definition of Atomic Species
#

Species.Number       2
<Definition.of.Atomic.Species
 C  C5.0-s2p1  C_CA13
 H  H5.0-s2    H_CA13
Definition.of.Atomic.Species>


NEGF.output_hks    true
NEGF.filename.hks  lead-l-8zgnr.hks


#
# Atoms
#

Atoms.Number       20
Atoms.SpeciesAndCoordinates.Unit   Ang # Ang|AU
<Atoms.SpeciesAndCoordinates 
   1  C     0.0000005     0.0361626     0.0012037     2.1000000     1.9000000
   2  C     0.0000005     0.7140385     1.2211879     2.1000000     1.9000000
   3  C     0.0000004     2.1457201     1.2213117     2.1000000     1.9000000
   4  C     0.0000001     2.8501725     0.0013209     2.1000000     1.9000000
   5  C    -0.0000001     4.2700744     0.0012023     2.1000000     1.9000000
   6  C    -0.0000002     4.9766306     1.2212146     2.1000000     1.9000000
   7  C     0.0000001     6.3942050     1.2213979     2.1000000     1.9000000
   8  C     0.0000001     7.1014018     0.0014044     2.1000000     1.9000000
   9  C    -0.0000001     8.5186603     0.0014042     2.1000000     1.9000000
  10  C     0.0000001     9.2258519     1.2213977     2.1000000     1.9000000
  11  C     0.0000004    10.6434033     1.2212143     2.1000000     1.9000000
  12  C     0.0000003    11.3499296     0.0012020     2.1000000     1.9000000
  13  H     0.0000003    -1.0818549     0.0012134     0.5000000     0.5000000
  14  H    -0.0000005    12.4679        0.0012129     0.5000000     0.5000000
  15  C     0.0000005     0.7140385     -1.2383     2.1000000     1.9000000
  16  C     0.0000004     2.1457201     -1.2383     2.1000000     1.9000000
  17 C    -0.0000002     4.9766306     -1.2383     2.1000000     1.9000000
  18  C     0.0000001     6.3942050     -1.2383     2.1000000     1.9000000
  19  C     0.0000001     9.2258519     -1.2383     2.1000000     1.9000000
  20  C     0.0000004    10.6434033     -1.2383     2.1000000     1.9000000
Atoms.SpeciesAndCoordinates>
Atoms.UnitVectors.Unit             Ang # Ang|AU
<Atoms.UnitVectors 
 0.000  0.000  6.1
 10.000 0.000 0.000 
 0.000 27.000 0.000 
Atoms.UnitVectors>


#
# SCF or Electronic System
#


scf.XcType                 LSDA-CA     # LDA|LSDA-CA|LSDA-PW|GGA-PBE
scf.SpinPolarization        on         # On|Off|NC
scf.ElectronicTemperature  300.0       # default=300 (K)
scf.energycutoff           120.0       # default=150 (Ry)
#scf.Ngrid                  140 140 100  # about=200 (Ry)
scf.maxIter                1000         # default=40
scf.EigenvalueSolver       Band     # DC|GDC|Cluster|Band
scf.lapack.dste            dstevx      # dstegr|dstedc|dstevx, default=dstegr
scf.Kgrid                  201 1 1       # means n1 x n2 x n3
scf.Mixing.Type           rmm-diisk    # Simple|Rmm-Diis|Gr-Pulay|Kerker|Rmm-Diisk
scf.Kerker.factor          10.0        # default=1
scf.Init.Mixing.Weight     0.001       # default=0.30
scf.Min.Mixing.Weight      0.001       # default=0.001
scf.Max.Mixing.Weight      0.105       # default=0.40
scf.Mixing.History          40       # default=5
scf.Mixing.StartPulay       10       # default=6
scf.Mixing.EveryPulay        1          # default=6
scf.criterion             1.0e-9       # default=1.0e-6 (Hartree)

#
# MD or Geometry Optimization
#

MD.Type               Nomd    # Opt|EF|BFGS|RF|DIIS
MD.Opt.DIIS.History   18    # default=3
MD.Opt.StartDIIS      6     # default=5
MD.Opt.EveryDIIS      5     # default=10
MD.maxIter            2000   #
MD.Opt.criterion    1.0e-4  # default=1.0e-4 (a.u.)


#
# Band dispersion
#

Band.dispersion              off      # on|off, default=off
<Band.KPath.UnitCell
 0.000 0.000  4.88
 10.000 0.000 0.000 
 0.000 27.000 0.000 
Band.KPath.UnitCell>

Band.Nkpath                1
<Band.kpath                
   601  0.0 0.0 0.0   0.5 0.0 0.0   G X
Band.kpath>


#
# DOS and PDOS
#

Dos.fileout                  off      # on|off, default=off
Dos.Erange              -20.0  20.0    # default = -20 20
Dos.Kgrid                601 1 1      # default = Kgrid1 Kgrid2 Kgrid3

DosGauss.fileout      off
DosGauss.Num.Mesh    4000
DosGauss.Width      0.01

MO.fileout                   off       # on|off
num.HOMOs                     1        # default=1
num.LUMOs                     1        # default=1
MO.Nkpoint                    1        # default=1
<MO.kpoint
  0.0  0.0  0.0
MO.kpoint>

五、输入效果图

这个是一个仿真结束之后的图,是我们的输入文件的原子排布的xyz文件。用xcrysden画出来

六、仿真结果

得到文件negf-8zgnr-0.3.dat,文件的第六列就是transmission。

#/************************************************************/
#
#  Current:
#
#    convert the unit from a.u to ampere
#
#    The unit of current is given by eEh/bar{h} in
#    atomic unit, where e is the elementary charge,
#    Eh Hartree, and bar{h} h/{2pi} given by
#
#    e = 1.60217653 * 10^{-19} C
#    Eh = 4.35974417 * 10^{-18} J
#    bar{h} = 1.05457168 * 10^{-34} J s
#
#    Therefore,
#
#    1 a.u.
#    = 1.60217653 * 10^{-19} C * 4.35974417 * 10^{-18} J
#    / (1.05457168 * 10^{-34} J s )
#    = 6.6236178 * 10^{-3} [Cs^{-1}=ampere]
#
#  Electric potential:
#
#    convert the unit from a.u to volt
#
#    The unit of electric potential is given by Eh/e in
#    atomic unit.
#
#    Therefore,
#
#    1 a.u.
#    = 4.35974417 * 10^{-18} J/ (1.60217653 * 10^{-19} C)
#    = 27.21138 [JC^{-1}=volt]
#
#    This allows us to consider that the difference of 1 eV in
#    the chemical potential of leads can be regarded as 1 volt.
#
#***********************************************************/
#
# k2=   0.0000  k3=  -0.0000
# SpinP_switch= 1
# Chemical potential (eV)      Left, Right Leads=-7.092117e+00 -7.382248e+00
# Chemical potential (Hartree) Left, Right Leads=-2.606305e-01 -2.712926e-01
# diff Chemical potential (eV)     = 2.901306e-01
# diff Chemical potential (Hartree)= 1.066210e-02
# tran_transmission_energydiv= 800
#
# iw w.real(au) w.imag(au) w.real(eV) w.imag(eV) trans.real(up) trans.imag(up) trans.real(down) trans.imag(down)
#
  0 -1.469973e-01  3.674932e-05 -4.000000e+00  1.000000e-03  1.504455e-10  1.521637e-24  1.071494e-10  1.816968e-25
  1 -1.466293e-01  3.674932e-05 -3.989987e+00  1.000000e-03  1.158592e-10  1.305401e-24  1.056194e-10  5.332574e-25
  2 -1.462614e-01  3.674932e-05 -3.979975e+00  1.000000e-03  9.618697e-11  1.190603e-24  1.508025e-10  1.605111e-24
  3 -1.458934e-01  3.674932e-05 -3.969962e+00  1.000000e-03  8.818859e-11  1.187069e-24  4.791283e-10  9.338231e-24
  4 -1.455255e-01  3.674932e-05 -3.959950e+00  1.000000e-03  1.072549e-10  1.517041e-24  1.124427e-08  3.594742e-22
  5 -1.451575e-01  3.674932e-05 -3.949937e+00  1.000000e-03  1.522907e-09  2.749582e-23  4.817887e-08  3.401657e-22
  6 -1.447896e-01  3.674932e-05 -3.939925e+00  1.000000e-03  1.626574e-10  8.560909e-24  2.177308e-08 -8.016507e-22
  7 -1.444216e-01  3.674932e-05 -3.929912e+00  1.000000e-03  9.518922e-10  1.539321e-22  6.147583e-09  2.502252e-22
  8 -1.440537e-01  3.674932e-05 -3.919900e+00  1.000000e-03  5.208821e-11  2.141974e-24  5.415061e-11  5.973337e-24
  9 -1.436857e-01  3.674932e-05 -3.909887e+00  1.000000e-03  4.604393e-11  1.121868e-24  1.741657e-11  1.166194e-24
 10 -1.433178e-01  3.674932e-05 -3.899875e+00  1.000000e-03  4.730882e-11  1.079651e-24  1.339429e-11  4.096322e-25
 11 -1.429498e-01  3.674932e-05 -3.889862e+00  1.000000e-03  5.243957e-11  1.253502e-24  1.152121e-11  2.036732e-25
 12 -1.425819e-01  3.674932e-05 -3.879850e+00  1.000000e-03  6.327638e-11  1.629360e-24  1.014272e-11  1.311523e-25
 13 -1.422139e-01  3.674932e-05 -3.869837e+00  1.000000e-03  8.614049e-11  2.396678e-24  9.064420e-12  1.020612e-25
 14 -1.418460e-01  3.674932e-05 -3.859825e+00  1.000000e-03  1.434529e-10  4.338940e-24  8.225351e-12  9.068025e-26
 15 -1.414780e-01  3.674932e-05 -3.849812e+00  1.000000e-03  3.665718e-10  1.168167e-23  7.588897e-12  8.751189e-26
 16 -1.411100e-01  3.674932e-05 -3.839800e+00  1.000000e-03  5.329663e-09  1.670105e-22  7.132836e-12  8.934003e-26
 17 -1.407421e-01  3.674932e-05 -3.829787e+00  1.000000e-03  8.216993e-10  3.142963e-23  6.849852e-12  9.587737e-26
 18 -1.403741e-01  3.674932e-05 -3.819775e+00  1.000000e-03  1.207264e-10  4.881948e-24  6.751390e-12  1.066073e-25
 19 -1.400062e-01  3.674932e-05 -3.809762e+00  1.000000e-03  4.093870e-11  1.742862e-24  6.877140e-12  1.233159e-25
 20 -1.396382e-01  3.674932e-05 -3.799750e+00  1.000000e-03  1.880501e-11  8.345210e-25  7.317351e-12  1.478377e-25
 21 -1.392703e-01  3.674932e-05 -3.789737e+00  1.000000e-03  1.015922e-11  4.621147e-25  8.267899e-12  1.862367e-25
 22 -1.389023e-01  3.674932e-05 -3.779725e+00  1.000000e-03  6.091086e-12  2.821568e-25  1.018285e-11  2.503448e-25
 23 -1.385344e-01  3.674932e-05 -3.769712e+00  1.000000e-03  3.950824e-12  1.810851e-25  1.428460e-11  3.701793e-25
 24 -1.381664e-01  3.674932e-05 -3.759700e+00  1.000000e-03  2.772059e-12  1.202969e-25  2.487813e-11  6.470174e-25
 25 -1.377985e-01  3.674932e-05 -3.749687e+00  1.000000e-03  2.187164e-12  8.038813e-26  6.802053e-11  1.560588e-24
 26 -1.374305e-01  3.674932e-05 -3.739675e+00  1.000000e-03  2.134428e-12  5.042547e-26  1.267607e-09  1.074098e-23
 27 -1.370626e-01  3.674932e-05 -3.729662e+00  1.000000e-03  2.869026e-12  2.564927e-26  1.426098e-10 -1.862343e-26
 28 -1.366946e-01  3.674932e-05 -3.719650e+00  1.000000e-03  5.360067e-12 -2.664920e-27  3.526859e-11 -8.543947e-25
 29 -1.363267e-01  3.674932e-05 -3.709637e+00  1.000000e-03  1.314221e-11 -4.911759e-26  4.678516e-11 -1.078344e-24
 30 -1.359587e-01  3.674932e-05 -3.699625e+00  1.000000e-03  4.386952e-11 -1.805363e-25  1.545314e-10 -1.704623e-24
 31 -1.355907e-01  3.674932e-05 -3.689612e+00  1.000000e-03  2.966769e-10 -1.017261e-24  1.126718e-09 -5.046092e-24
 32 -1.352228e-01  3.674932e-05 -3.679599e+00  1.000000e-03  7.100023e-09 -1.789084e-23  7.247211e-08 -4.207550e-23
 33 -1.348548e-01  3.674932e-05 -3.669587e+00  1.000000e-03  2.822794e-10 -6.667357e-25  2.275713e-09 -1.096274e-24
 34 -1.344869e-01  3.674932e-05 -3.659574e+00  1.000000e-03  1.127237e-10 -2.112687e-25  1.265338e-09  5.022547e-26
 35 -1.341189e-01  3.674932e-05 -3.649562e+00  1.000000e-03  7.179669e-11 -1.148708e-25  1.484029e-09  4.069964e-25
 36 -1.337510e-01  3.674932e-05 -3.639549e+00  1.000000e-03  5.625037e-11 -7.810395e-26  9.781181e-09  3.237784e-24
 37 -1.333830e-01  3.674932e-05 -3.629537e+00  1.000000e-03  4.969053e-11 -6.188473e-26  3.467560e-08  1.070045e-23
 38 -1.330151e-01  3.674932e-05 -3.619524e+00  1.000000e-03  4.761920e-11 -5.431095e-26  3.342954e-08  9.359108e-24
 39 -1.326471e-01  3.674932e-05 -3.609512e+00  1.000000e-03  4.855281e-11 -5.155850e-26  2.605274e-08  5.942553e-24
 40 -1.322792e-01  3.674932e-05 -3.599499e+00  1.000000e-03  5.214527e-11 -5.274451e-26  1.263828e-08  2.259038e-24
 41 -1.319112e-01  3.674932e-05 -3.589487e+00  1.000000e-03  5.872580e-11 -5.792881e-26  5.780934e-10  7.460541e-26
 42 -1.315433e-01  3.674932e-05 -3.579474e+00  1.000000e-03  6.931780e-11 -6.555243e-26  2.056356e-10  1.672435e-26
 43 -1.311753e-01  3.674932e-05 -3.569462e+00  1.000000e-03  8.603865e-11 -7.967716e-26  1.074563e-10  5.042484e-27
 44 -1.308074e-01  3.674932e-05 -3.559449e+00  1.000000e-03  1.132055e-10 -1.041915e-25  6.606227e-11  2.491676e-28
 45 -1.304394e-01  3.674932e-05 -3.549437e+00  1.000000e-03  1.604048e-10 -1.452273e-25  4.472853e-11 -1.769312e-27
 46 -1.300715e-01  3.674932e-05 -3.539424e+00  1.000000e-03  2.526890e-10 -2.273466e-25  3.235554e-11 -2.575535e-27
 47 -1.297035e-01  3.674932e-05 -3.529412e+00  1.000000e-03  4.781689e-10 -4.323546e-25  2.458560e-11 -2.854848e-27
 48 -1.293355e-01  3.674932e-05 -3.519399e+00  1.000000e-03  1.495405e-09 -1.353496e-24  1.941658e-11 -3.083979e-27
 49 -1.289676e-01  3.674932e-05 -3.509387e+00  1.000000e-03  2.768624e-08 -2.509132e-23  1.582421e-11 -3.166232e-27
 50 -1.285996e-01  3.674932e-05 -3.499374e+00  1.000000e-03  4.312680e-08 -3.951316e-23  1.324137e-11 -3.288564e-27
 51 -1.282317e-01  3.674932e-05 -3.489362e+00  1.000000e-03  4.713764e-08 -4.443972e-23  1.133433e-11 -3.360617e-27
 52 -1.278637e-01  3.674932e-05 -3.479349e+00  1.000000e-03  3.923903e-08 -3.736261e-23  9.896896e-12 -3.458078e-27
 53 -1.274958e-01  3.674932e-05 -3.469337e+00  1.000000e-03  3.755944e-09 -3.694013e-24  8.796633e-12 -3.594925e-27
 54 -1.271278e-01  3.674932e-05 -3.459324e+00  1.000000e-03  9.047723e-10 -9.105710e-25  7.945706e-12 -3.730936e-27
 55 -1.267599e-01  3.674932e-05 -3.449312e+00  1.000000e-03  5.036687e-10 -5.240214e-25  7.284385e-12 -3.893545e-27
 56 -1.263919e-01  3.674932e-05 -3.439299e+00  1.000000e-03  3.523417e-10 -3.782280e-25  6.771292e-12 -4.130408e-27
 57 -1.260240e-01  3.674932e-05 -3.429287e+00  1.000000e-03  2.796797e-10 -3.136106e-25  6.377455e-12 -4.444875e-27
 58 -1.256560e-01  3.674932e-05 -3.419274e+00  1.000000e-03  2.428110e-10 -2.859129e-25  6.082588e-12 -4.794995e-27
 59 -1.252881e-01  3.674932e-05 -3.409262e+00  1.000000e-03  2.270653e-10 -2.804191e-25  5.872748e-12 -5.201756e-27
 60 -1.249201e-01  3.674932e-05 -3.399249e+00  1.000000e-03  2.277420e-10 -2.973182e-25  5.738861e-12 -5.740843e-27
 61 -1.245522e-01  3.674932e-05 -3.389237e+00  1.000000e-03  2.460753e-10 -3.432497e-25  5.675850e-12 -6.460071e-27
 62 -1.241842e-01  3.674932e-05 -3.379224e+00  1.000000e-03  2.909510e-10 -4.374288e-25  5.682215e-12 -7.366291e-27
 63 -1.238162e-01  3.674932e-05 -3.369212e+00  1.000000e-03  3.910178e-10 -6.343367e-25  5.759997e-12 -8.630228e-27
 64 -1.234483e-01  3.674932e-05 -3.359199e+00  1.000000e-03  6.650584e-10 -1.177298e-24  5.915154e-12 -1.014495e-26
 65 -1.230803e-01  3.674932e-05 -3.349186e+00  1.000000e-03  2.868030e-09 -5.495228e-24  6.158423e-12 -1.223142e-26
 66 -1.227124e-01  3.674932e-05 -3.339174e+00  1.000000e-03  2.463514e-08 -5.278808e-23  6.506912e-12 -1.531251e-26
 67 -1.223444e-01  3.674932e-05 -3.329161e+00  1.000000e-03  2.681316e-08 -6.574997e-23  6.986817e-12 -1.942004e-26
 68 -1.219765e-01  3.674932e-05 -3.319149e+00  1.000000e-03  2.233045e-08 -6.445702e-23  7.638102e-12 -2.568374e-26
 69 -1.216085e-01  3.674932e-05 -3.309136e+00  1.000000e-03  1.274037e-08 -4.378581e-23  8.522741e-12 -3.541486e-26
 70 -1.212406e-01  3.674932e-05 -3.299124e+00  1.000000e-03  6.128597e-10 -2.509816e-24  9.739936e-12 -5.105997e-26
 71 -1.208726e-01  3.674932e-05 -3.289111e+00  1.000000e-03  1.885553e-10 -9.419296e-25  1.145595e-11 -7.810440e-26
 72 -1.205047e-01  3.674932e-05 -3.279099e+00  1.000000e-03  9.293218e-11 -5.931038e-25  1.396749e-11 -1.279084e-25
 73 -1.201367e-01  3.674932e-05 -3.269086e+00  1.000000e-03  5.520950e-11 -4.625778e-25  1.785136e-11 -2.349648e-25
 74 -1.197688e-01  3.674932e-05 -3.259074e+00  1.000000e-03  3.740481e-11 -4.308453e-25  2.437325e-11 -4.930322e-25
 75 -1.194008e-01  3.674932e-05 -3.249061e+00  1.000000e-03  2.940056e-11 -4.755947e-25  3.687277e-11 -1.295415e-24
 76 -1.190329e-01  3.674932e-05 -3.239049e+00  1.000000e-03  3.004818e-11 -6.766494e-25  6.749555e-11 -4.972459e-24
 77 -1.186649e-01  3.674932e-05 -3.229036e+00  1.000000e-03  5.643561e-11 -1.543457e-24  2.103918e-10 -5.030076e-23
 78 -1.182969e-01  3.674932e-05 -3.219024e+00  1.000000e-03  9.152355e-10 -2.238614e-23  7.132123e-08 -9.537630e-20
 79 -1.179290e-01  3.674932e-05 -3.209011e+00  1.000000e-03  1.370247e-10 -2.491376e-24  5.429508e-09 -1.552977e-21
 80 -1.175610e-01  3.674932e-05 -3.198999e+00  1.000000e-03  2.944216e-11 -3.468291e-25  4.428135e-09 -4.090001e-22
 81 -1.171931e-01  3.674932e-05 -3.188986e+00  1.000000e-03  1.420050e-11 -1.098330e-25  3.288562e-09 -1.455189e-22
 82 -1.168251e-01  3.674932e-05 -3.178974e+00  1.000000e-03  8.956238e-12 -4.407557e-26  5.405029e-10 -1.381539e-23
 83 -1.164572e-01  3.674932e-05 -3.168961e+00  1.000000e-03  6.410662e-12 -1.632631e-26  6.443738e-11 -9.762806e-25
 84 -1.160892e-01  3.674932e-05 -3.158949e+00  1.000000e-03  4.948325e-12  2.720330e-27  2.981577e-11 -2.478668e-25
 85 -1.157213e-01  3.674932e-05 -3.148936e+00  1.000000e-03  4.088220e-12  4.227336e-26  1.774907e-11  1.151085e-26
 86 -1.153533e-01  3.674932e-05 -3.138924e+00  1.000000e-03  5.129022e-12  8.824905e-25  1.265900e-11  7.545932e-25
 87 -1.149854e-01  3.674932e-05 -3.128911e+00  1.000000e-03  2.619178e-12  1.468834e-25  9.149691e-12  1.137743e-24
 88 -1.146174e-01  3.674932e-05 -3.118899e+00  1.000000e-03  2.220264e-12  2.341624e-26  6.276022e-12  6.136421e-26
 89 -1.142495e-01  3.674932e-05 -3.108886e+00  1.000000e-03  1.951705e-12  8.050330e-27  4.901522e-12  7.820913e-27
 90 -1.138815e-01  3.674932e-05 -3.098874e+00  1.000000e-03  1.739869e-12  3.818837e-27  3.970001e-12 -6.485789e-28
 91 -1.135136e-01  3.674932e-05 -3.088861e+00  1.000000e-03  1.572246e-12  2.091528e-27  3.308960e-12 -2.332732e-27
 92 -1.131456e-01  3.674932e-05 -3.078849e+00  1.000000e-03  1.440483e-12  1.347384e-27  2.828956e-12 -2.436993e-27
 93 -1.127777e-01  3.674932e-05 -3.068836e+00  1.000000e-03  1.338598e-12  9.173054e-28  2.476767e-12 -2.034295e-27
 94 -1.124097e-01  3.674932e-05 -3.058824e+00  1.000000e-03  1.262439e-12  7.394275e-28  2.219198e-12 -1.577361e-27
 95 -1.120417e-01  3.674932e-05 -3.048811e+00  1.000000e-03  1.209344e-12  6.328542e-28  2.035228e-12 -1.121657e-27
 96 -1.116738e-01  3.674932e-05 -3.038798e+00  1.000000e-03  1.177922e-12  6.117639e-28  1.911802e-12 -7.217091e-28
 97 -1.113058e-01  3.674932e-05 -3.028786e+00  1.000000e-03  1.167973e-12  6.240918e-28  1.841554e-12 -3.239118e-28
 98 -1.109379e-01  3.674932e-05 -3.018773e+00  1.000000e-03  1.180532e-12  6.850789e-28  1.821660e-12  5.953935e-29
 99 -1.105699e-01  3.674932e-05 -3.008761e+00  1.000000e-03  1.218055e-12  8.048018e-28  1.853478e-12  4.889071e-28
100 -1.102020e-01  3.674932e-05 -2.998748e+00  1.000000e-03  1.284790e-12  9.582122e-28  1.942883e-12  9.832756e-28
101 -1.098340e-01  3.674932e-05 -2.988736e+00  1.000000e-03  1.387418e-12  1.198643e-27  2.101378e-12  1.563683e-27
102 -1.094661e-01  3.674932e-05 -2.978723e+00  1.000000e-03  1.536129e-12  1.518037e-27  2.348314e-12  2.320628e-27
103 -1.090981e-01  3.674932e-05 -2.968711e+00  1.000000e-03  1.746403e-12  1.980186e-27  2.714908e-12  3.382600e-27
104 -1.087302e-01  3.674932e-05 -2.958698e+00  1.000000e-03  2.042055e-12  2.634865e-27  3.251463e-12  4.882699e-27
105 -1.083622e-01  3.674932e-05 -2.948686e+00  1.000000e-03  2.460609e-12  3.599825e-27  4.040705e-12  7.078544e-27
106 -1.079943e-01  3.674932e-05 -2.938673e+00  1.000000e-03  3.063184e-12  5.094057e-27  5.223495e-12  1.048589e-26
107 -1.076263e-01  3.674932e-05 -2.928661e+00  1.000000e-03  3.953735e-12  7.389908e-27  7.051481e-12  1.604270e-26
108 -1.072584e-01  3.674932e-05 -2.918648e+00  1.000000e-03  5.319081e-12  1.112131e-26  1.000337e-11  2.549600e-26
109 -1.068904e-01  3.674932e-05 -2.908636e+00  1.000000e-03  7.519534e-12  1.777441e-26  1.506784e-11  4.300300e-26
110 -1.065224e-01  3.674932e-05 -2.898623e+00  1.000000e-03  1.131769e-11  3.025201e-26  2.452607e-11  7.844788e-26
111 -1.061545e-01  3.674932e-05 -2.888611e+00  1.000000e-03  1.854678e-11  5.636480e-26  4.453636e-11  1.604613e-25
112 -1.057865e-01  3.674932e-05 -2.878598e+00  1.000000e-03  3.451546e-11  1.213547e-25  9.628040e-11  3.962858e-25
113 -1.054186e-01  3.674932e-05 -2.868586e+00  1.000000e-03  8.002424e-11  3.348970e-25  2.923406e-10  1.398144e-24
114 -1.050506e-01  3.674932e-05 -2.858573e+00  1.000000e-03  3.046347e-10  1.618561e-24  2.528210e-09  1.463738e-23
115 -1.046827e-01  3.674932e-05 -2.848561e+00  1.000000e-03  1.519410e-08  1.234596e-22  4.555257e-09  3.535453e-23
116 -1.043147e-01  3.674932e-05 -2.838548e+00  1.000000e-03  6.175709e-10  1.180077e-23  5.203039e-10  7.485546e-24
117 -1.039468e-01  3.674932e-05 -2.828536e+00  1.000000e-03  2.625862e-10  1.796733e-22  2.664489e-10  7.078965e-23
118 -1.035788e-01  3.674932e-05 -2.818523e+00  1.000000e-03  7.165168e-11  1.278857e-24  1.361795e-10  2.718514e-24
119 -1.032109e-01  3.674932e-05 -2.808511e+00  1.000000e-03  4.281935e-11  4.563581e-25  9.676243e-11  1.152334e-24
120 -1.028429e-01  3.674932e-05 -2.798498e+00  1.000000e-03  2.971567e-11  3.226415e-25  7.902689e-11  1.005880e-24
121 -1.024750e-01  3.674932e-05 -2.788486e+00  1.000000e-03  2.249908e-11  2.894296e-25  7.068562e-11  1.100359e-24
122 -1.021070e-01  3.674932e-05 -2.778473e+00  1.000000e-03  1.798993e-11  2.995107e-25  6.791805e-11  1.394816e-24
123 -1.017391e-01  3.674932e-05 -2.768461e+00  1.000000e-03  1.482840e-11  3.462459e-25  6.937795e-11  2.031653e-24
124 -1.013711e-01  3.674932e-05 -2.758448e+00  1.000000e-03  1.225728e-11  4.545375e-25  7.472338e-11  3.505963e-24
125 -1.010031e-01  3.674932e-05 -2.748436e+00  1.000000e-03  9.628229e-12  7.150483e-25  8.351072e-11  7.893066e-24
126 -1.006352e-01  3.674932e-05 -2.738423e+00  1.000000e-03  5.976104e-12  1.623188e-24  9.425322e-11  2.940393e-23
127 -1.002672e-01  3.674932e-05 -2.728411e+00  1.000000e-03  6.003821e-12  1.151748e-23  8.325935e-10  6.264401e-22
128 -9.989929e-02  3.674932e-05 -2.718398e+00  1.000000e-03  1.110016e-10  2.131048e-23  4.769916e-07  8.501360e-20
129 -9.953134e-02  3.674932e-05 -2.708385e+00  1.000000e-03  3.813614e-11  3.894623e-24  2.059684e-08  2.323462e-21
130 -9.916338e-02  3.674932e-05 -2.698373e+00  1.000000e-03  2.945818e-11  3.612176e-24  8.391792e-10  1.327036e-22
131 -9.879543e-02  3.674932e-05 -2.688360e+00  1.000000e-03  2.617559e-11  6.217046e-24  1.659393e-10  6.244065e-23
132 -9.842748e-02  3.674932e-05 -2.678348e+00  1.000000e-03  2.265276e-11  1.823029e-23  3.823377e-11  7.129448e-23
133 -9.805952e-02  3.674932e-05 -2.668335e+00  1.000000e-03  2.925152e-11  1.883479e-22  1.369402e-10  3.016190e-22
134 -9.769157e-02  3.674932e-05 -2.658323e+00  1.000000e-03  2.059304e-10  2.140251e-22  6.306279e-10  2.732462e-22
135 -9.732362e-02  3.674932e-05 -2.648310e+00  1.000000e-03  1.148746e-10  3.299027e-23  1.128033e-10  1.900576e-23
136 -9.695566e-02  3.674932e-05 -2.638298e+00  1.000000e-03  1.249063e-10  1.660202e-23  5.488178e-11  4.903417e-24
137 -9.658771e-02  3.674932e-05 -2.628285e+00  1.000000e-03  1.639158e-10  1.270710e-23  3.428359e-11  1.904495e-24
138 -9.621976e-02  3.674932e-05 -2.618273e+00  1.000000e-03  2.482167e-10  1.274912e-23  2.407099e-11  9.175051e-25
139 -9.585181e-02  3.674932e-05 -2.608260e+00  1.000000e-03  4.498758e-10  1.659239e-23  1.814443e-11  5.047947e-25
140 -9.548385e-02  3.674932e-05 -2.598248e+00  1.000000e-03  1.128821e-09  3.163843e-23  1.437483e-11  3.057294e-25
141 -9.511590e-02  3.674932e-05 -2.588235e+00  1.000000e-03  1.055105e-08  2.355175e-22  1.182961e-11  1.991198e-25
142 -9.474795e-02  3.674932e-05 -2.578223e+00  1.000000e-03  8.496405e-08  1.519663e-21  1.003885e-11  1.371681e-25
143 -9.437999e-02  3.674932e-05 -2.568210e+00  1.000000e-03  8.026925e-09  1.186398e-22  8.742198e-12  9.909886e-26
144 -9.401204e-02  3.674932e-05 -2.558198e+00  1.000000e-03  1.283713e-09  1.623884e-23  7.785248e-12  7.438082e-26
145 -9.364409e-02  3.674932e-05 -2.548185e+00  1.000000e-03  5.921542e-10  6.472751e-24  7.071701e-12  5.790998e-26
146 -9.327613e-02  3.674932e-05 -2.538173e+00  1.000000e-03  3.689535e-10  3.531696e-24  6.539214e-12  4.635924e-26
147 -9.290818e-02  3.674932e-05 -2.528160e+00  1.000000e-03  2.705662e-10  2.291293e-24  6.146371e-12  3.820148e-26
148 -9.254023e-02  3.674932e-05 -2.518148e+00  1.000000e-03  2.217758e-10  1.679925e-24  5.865284e-12  3.226059e-26
149 -9.217227e-02  3.674932e-05 -2.508135e+00  1.000000e-03  1.985246e-10  1.353566e-24  5.677225e-12  2.781922e-26
150 -9.180432e-02  3.674932e-05 -2.498123e+00  1.000000e-03  1.922184e-10  1.188105e-24  5.569995e-12  2.455697e-26
151 -9.143637e-02  3.674932e-05 -2.488110e+00  1.000000e-03  2.011917e-10  1.134804e-24  5.536352e-12  2.212126e-26
152 -9.106841e-02  3.674932e-05 -2.478098e+00  1.000000e-03  2.295478e-10  1.191149e-24  5.573122e-12  2.027515e-26
153 -9.070046e-02  3.674932e-05 -2.468085e+00  1.000000e-03  2.916113e-10  1.394436e-24  5.680811e-12  1.897511e-26
154 -9.033251e-02  3.674932e-05 -2.458073e+00  1.000000e-03  4.325267e-10  1.917078e-24  5.863639e-12  1.806550e-26
155 -8.996455e-02  3.674932e-05 -2.448060e+00  1.000000e-03  8.475837e-10  3.496506e-24  6.129991e-12  1.751960e-26
156 -8.959660e-02  3.674932e-05 -2.438048e+00  1.000000e-03  3.978584e-09  1.544326e-23  6.493390e-12  1.729318e-26
157 -8.922865e-02  3.674932e-05 -2.428035e+00  1.000000e-03  5.387061e-08  1.931055e-22  6.974192e-12  1.738111e-26
158 -8.886069e-02  3.674932e-05 -2.418023e+00  1.000000e-03  9.577381e-09  3.174242e-23  7.602387e-12  1.789647e-26
159 -8.849274e-02  3.674932e-05 -2.408010e+00  1.000000e-03  6.772481e-10  2.128822e-24  8.422238e-12  1.873287e-26
160 -8.812479e-02  3.674932e-05 -2.397997e+00  1.000000e-03  2.383704e-10  7.045014e-25  9.500105e-12  2.000060e-26
161 -8.775684e-02  3.674932e-05 -2.387985e+00  1.000000e-03  1.173251e-10  3.270302e-25  1.093812e-11  2.192108e-26
162 -8.738888e-02  3.674932e-05 -2.377972e+00  1.000000e-03  6.819772e-11  1.778604e-25  1.289917e-11  2.469016e-26
163 -8.702093e-02  3.674932e-05 -2.367960e+00  1.000000e-03  4.399727e-11  1.074692e-25  1.565531e-11  2.857411e-26
164 -8.665298e-02  3.674932e-05 -2.357947e+00  1.000000e-03  3.056516e-11  6.969610e-26  1.968832e-11  3.432995e-26
165 -8.628502e-02  3.674932e-05 -2.347935e+00  1.000000e-03  2.247962e-11  4.795111e-26  2.591866e-11  4.344044e-26
166 -8.591707e-02  3.674932e-05 -2.337922e+00  1.000000e-03  1.732133e-11  3.407185e-26  3.629215e-11  5.789789e-26
167 -8.554912e-02  3.674932e-05 -2.327910e+00  1.000000e-03  1.388836e-11  2.510856e-26  5.555069e-11  8.463581e-26
168 -8.518116e-02  3.674932e-05 -2.317897e+00  1.000000e-03  1.153407e-11  1.904695e-26  9.806116e-11  1.422776e-25
169 -8.481321e-02  3.674932e-05 -2.307885e+00  1.000000e-03  9.888682e-12  1.481006e-26  2.282954e-10  3.140932e-25
170 -8.444526e-02  3.674932e-05 -2.297872e+00  1.000000e-03  8.730557e-12  1.171720e-26  1.388937e-09  1.775064e-24
171 -8.407730e-02  3.674932e-05 -2.287860e+00  1.000000e-03  7.922162e-12  9.515314e-27  1.819738e-08  2.082512e-23
172 -8.370935e-02  3.674932e-05 -2.277847e+00  1.000000e-03  7.376229e-12  7.764705e-27  3.318357e-09  3.402086e-24
173 -8.334140e-02  3.674932e-05 -2.267835e+00  1.000000e-03  7.037085e-12  6.464989e-27  2.930134e-10  2.678657e-25
174 -8.297344e-02  3.674932e-05 -2.257822e+00  1.000000e-03  6.870014e-12  5.435121e-27  1.195917e-10  9.071509e-26
175 -8.260549e-02  3.674932e-05 -2.247810e+00  1.000000e-03  6.855160e-12  4.671194e-27  6.805458e-11  4.069716e-26
176 -8.223754e-02  3.674932e-05 -2.237797e+00  1.000000e-03  6.984155e-12  4.080118e-27  4.589452e-11  1.898068e-26
177 -8.186958e-02  3.674932e-05 -2.227785e+00  1.000000e-03  7.258572e-12  3.592488e-27  3.457531e-11  7.282051e-27
178 -8.150163e-02  3.674932e-05 -2.217772e+00  1.000000e-03  7.689798e-12  3.144739e-27  2.828659e-11 -6.004491e-29
179 -8.113368e-02  3.674932e-05 -2.207760e+00  1.000000e-03  8.300258e-12  2.959560e-27  2.473776e-11 -5.435966e-27
180 -8.076572e-02  3.674932e-05 -2.197747e+00  1.000000e-03  9.126261e-12  2.794066e-27  2.290248e-11 -1.029665e-26
181 -8.039777e-02  3.674932e-05 -2.187735e+00  1.000000e-03  1.022315e-11  2.729758e-27  2.230192e-11 -1.491137e-26
182 -8.002982e-02  3.674932e-05 -2.177722e+00  1.000000e-03  1.167407e-11  2.812580e-27  2.274100e-11 -1.997260e-26
183 -7.966186e-02  3.674932e-05 -2.167710e+00  1.000000e-03  1.360508e-11  3.105540e-27  2.421027e-11 -2.644213e-26
184 -7.929391e-02  3.674932e-05 -2.157697e+00  1.000000e-03  1.621158e-11  3.673366e-27  2.686705e-11 -3.464589e-26
185 -7.892596e-02  3.674932e-05 -2.147685e+00  1.000000e-03  1.980688e-11  4.794791e-27  3.107622e-11 -4.605576e-26
186 -7.855801e-02  3.674932e-05 -2.137672e+00  1.000000e-03  2.491620e-11  7.066465e-27  3.753138e-11 -6.289122e-26
187 -7.819005e-02  3.674932e-05 -2.127660e+00  1.000000e-03  3.247171e-11  1.026466e-26  4.753413e-11 -8.911743e-26
188 -7.782210e-02  3.674932e-05 -2.117647e+00  1.000000e-03  4.425387e-11  1.688146e-26  6.365368e-11 -1.317723e-25
189 -7.745415e-02  3.674932e-05 -2.107635e+00  1.000000e-03  6.400910e-11  2.873077e-26  9.145030e-11 -2.100196e-25
190 -7.708619e-02  3.674932e-05 -2.097622e+00  1.000000e-03  1.007444e-10  5.657672e-26  1.447419e-10 -3.694054e-25
191 -7.671824e-02  3.674932e-05 -2.087610e+00  1.000000e-03  1.808068e-10  1.209383e-25  2.659388e-10 -7.574535e-25
192 -7.635029e-02  3.674932e-05 -2.077597e+00  1.000000e-03  4.106810e-10  3.252432e-25  6.415674e-10 -2.064930e-24
193 -7.598233e-02  3.674932e-05 -2.067584e+00  1.000000e-03  1.634727e-09  1.506390e-24  3.109938e-09 -1.147759e-23
194 -7.561438e-02  3.674932e-05 -2.057572e+00  1.000000e-03  1.700695e-07  1.984503e-22  6.851535e-08 -2.763164e-22
195 -7.524643e-02  3.674932e-05 -2.047559e+00  1.000000e-03  1.895991e-09  2.883017e-24  1.847612e-09 -8.921721e-24
196 -7.487847e-02  3.674932e-05 -2.037547e+00  1.000000e-03  4.901338e-10  8.290215e-25  5.827117e-10 -3.262446e-24
197 -7.451052e-02  3.674932e-05 -2.027534e+00  1.000000e-03  2.300157e-10  4.724304e-25  3.102663e-10 -2.025558e-24
198 -7.414257e-02  3.674932e-05 -2.017522e+00  1.000000e-03  1.393284e-10  3.362060e-25  2.182343e-10 -1.646650e-24
199 -7.377461e-02  3.674932e-05 -2.007509e+00  1.000000e-03  9.826279e-11  2.754421e-25  1.919769e-10 -1.644784e-24
200 -7.340666e-02  3.674932e-05 -1.997497e+00  1.000000e-03  7.734587e-11  2.469836e-25  2.136467e-10 -1.959273e-24
201 -7.303871e-02  3.674932e-05 -1.987484e+00  1.000000e-03  6.664714e-11  2.323881e-25  3.186788e-10 -2.959670e-24
202 -7.267075e-02  3.674932e-05 -1.977472e+00  1.000000e-03  6.230327e-11  2.347724e-25  7.389773e-10 -6.509050e-24
203 -7.230280e-02  3.674932e-05 -1.967459e+00  1.000000e-03  6.292694e-11  2.459057e-25  5.020654e-09 -4.044941e-23
204 -7.193485e-02  3.674932e-05 -1.957447e+00  1.000000e-03  6.854227e-11  2.697226e-25  1.124344e-08 -7.826845e-23
205 -7.156689e-02  3.674932e-05 -1.947434e+00  1.000000e-03  8.047425e-11  3.067714e-25  1.085324e-09 -6.448350e-24
206 -7.119894e-02  3.674932e-05 -1.937422e+00  1.000000e-03  1.019698e-10  3.652061e-25  4.290936e-10 -2.190461e-24
207 -7.083099e-02  3.674932e-05 -1.927409e+00  1.000000e-03  1.401000e-10  4.595479e-25  2.567122e-10 -1.110054e-24
208 -7.046304e-02  3.674932e-05 -1.917397e+00  1.000000e-03  2.111263e-10  6.141428e-25  1.911190e-10 -7.008479e-25
209 -7.009508e-02  3.674932e-05 -1.907384e+00  1.000000e-03  3.578235e-10  9.024498e-25  1.659155e-10 -5.077893e-25
210 -6.972713e-02  3.674932e-05 -1.897372e+00  1.000000e-03  7.209766e-10  1.540398e-24  1.648237e-10 -4.130191e-25
211 -6.935918e-02  3.674932e-05 -1.887359e+00  1.000000e-03  1.992029e-09  3.470415e-24  1.888455e-10 -3.732299e-25
212 -6.899122e-02  3.674932e-05 -1.877347e+00  1.000000e-03  1.326678e-08  1.898848e-23  2.582188e-10 -4.088399e-25
213 -6.862327e-02  3.674932e-05 -1.867334e+00  1.000000e-03  5.925457e-08  5.932667e-23  4.568955e-10 -6.202743e-25
214 -6.825532e-02  3.674932e-05 -1.857322e+00  1.000000e-03  4.381058e-09  3.425656e-24  1.313171e-09 -1.592769e-24
215 -6.788736e-02  3.674932e-05 -1.847309e+00  1.000000e-03  1.629064e-09  9.852609e-25  2.334978e-08 -3.919192e-23
216 -6.751941e-02  3.674932e-05 -1.837297e+00  1.000000e-03  9.207737e-10  4.113180e-25  4.069637e-09 -8.266554e-24
217 -6.715146e-02  3.674932e-05 -1.827284e+00  1.000000e-03  6.378600e-10  2.689170e-25  7.654046e-10 -2.285235e-24
218 -6.678350e-02  3.674932e-05 -1.817272e+00  1.000000e-03  5.016436e-10  2.818548e-25  3.366600e-10 -1.417482e-24
219 -6.641555e-02  3.674932e-05 -1.807259e+00  1.000000e-03  4.328842e-10  4.511876e-25  2.038374e-10 -1.172898e-24
220 -6.604760e-02  3.674932e-05 -1.797247e+00  1.000000e-03  4.034805e-10  8.558355e-25  1.477749e-10 -1.144466e-24
221 -6.567964e-02  3.674932e-05 -1.787234e+00  1.000000e-03  4.039278e-10  1.776993e-24  1.226970e-10 -1.249769e-24
222 -6.531169e-02  3.674932e-05 -1.777222e+00  1.000000e-03  4.348931e-10  4.173669e-24  1.216853e-10 -1.498325e-24
223 -6.494374e-02  3.674932e-05 -1.767209e+00  1.000000e-03  5.069559e-10  1.211430e-23  5.420624e-10 -1.324845e-24
224 -6.457578e-02  3.674932e-05 -1.757196e+00  1.000000e-03  6.448509e-10  5.331499e-23  1.294629e-09 -8.653993e-25
225 -6.420783e-02  3.674932e-05 -1.747184e+00  1.000000e-03  8.412555e-10  8.540600e-22  1.815162e-09 -2.029278e-24
226 -6.383988e-02  3.674932e-05 -1.737171e+00  1.000000e-03  2.274791e-09  1.570686e-21  2.298141e-09 -1.139150e-23
227 -6.347192e-02  3.674932e-05 -1.727159e+00  1.000000e-03  5.189058e-09  5.626120e-22  3.368100e-09 -2.865749e-23
228 -6.310397e-02  3.674932e-05 -1.717146e+00  1.000000e-03  3.073655e-08  1.450626e-21  3.302295e-09  5.136355e-24
229 -6.273602e-02  3.674932e-05 -1.707134e+00  1.000000e-03  2.233511e-07  6.360471e-21  3.743012e-09  6.352640e-24
230 -6.236807e-02  3.674932e-05 -1.697121e+00  1.000000e-03  1.436467e-08  2.801368e-22  4.228243e-09  6.222131e-24
231 -6.200011e-02  3.674932e-05 -1.687109e+00  1.000000e-03  6.066772e-09  8.383812e-23  4.746905e-09  8.733485e-24
232 -6.163216e-02  3.674932e-05 -1.677096e+00  1.000000e-03  4.350081e-09  3.982069e-23  5.302439e-09  1.386926e-23
233 -6.126421e-02  3.674932e-05 -1.667084e+00  1.000000e-03  4.419309e-09  1.791988e-23  5.904768e-09  1.480101e-23
234 -6.089625e-02  3.674932e-05 -1.657071e+00  1.000000e-03  6.666442e-09 -1.199875e-23  6.711931e-09  5.048255e-24
235 -6.052830e-02  3.674932e-05 -1.647059e+00  1.000000e-03  2.259998e-08 -1.851127e-22  7.603099e-09  1.063152e-23
236 -6.016035e-02  3.674932e-05 -1.637046e+00  1.000000e-03  4.785712e-07 -5.364440e-21  8.587864e-09  2.001682e-23
237 -5.979239e-02  3.674932e-05 -1.627034e+00  1.000000e-03  1.586955e-08 -4.733979e-23  9.713972e-09  2.044248e-23
238 -5.942444e-02  3.674932e-05 -1.617021e+00  1.000000e-03  1.335280e-08  1.744469e-22  1.103755e-08  2.342096e-23
239 -5.905649e-02  3.674932e-05 -1.607009e+00  1.000000e-03  9.611082e-08  2.287650e-21  1.265109e-08  4.681123e-23
240 -5.868853e-02  3.674932e-05 -1.596996e+00  1.000000e-03  3.313356e-08  8.765979e-22  1.473694e-08  7.111269e-23
241 -5.832058e-02  3.674932e-05 -1.586984e+00  1.000000e-03  4.762835e-09  1.140928e-22  1.773066e-08  1.124978e-22
242 -5.795263e-02  3.674932e-05 -1.576971e+00  1.000000e-03  2.026994e-09  3.817717e-23  2.296956e-08  2.163570e-22
243 -5.758467e-02  3.674932e-05 -1.566959e+00  1.000000e-03  1.367278e-09  1.738767e-23  3.643275e-08  4.187392e-22
244 -5.721672e-02  3.674932e-05 -1.556946e+00  1.000000e-03  1.361259e-09  9.267159e-24  1.126001e-07  8.557603e-22
245 -5.684877e-02  3.674932e-05 -1.546934e+00  1.000000e-03  7.394166e-09  7.847311e-24  1.440245e-05 -1.289366e-19
246 -5.648081e-02  3.674932e-05 -1.536921e+00  1.000000e-03  2.411678e-08  1.271548e-23  3.162604e-07  8.646370e-21
247 -5.611286e-02  3.674932e-05 -1.526909e+00  1.000000e-03  3.882820e-08  2.483077e-23  2.940631e-07 -4.670133e-21
248 -5.574491e-02  3.674932e-05 -1.516896e+00  1.000000e-03  5.593093e-08  4.410230e-23  3.712698e-06 -3.773859e-20
249 -5.537695e-02  3.674932e-05 -1.506884e+00  1.000000e-03  7.751932e-08  6.994200e-23  3.868446e-06 -5.055653e-20
250 -5.500900e-02  3.674932e-05 -1.496871e+00  1.000000e-03  1.062384e-07  1.086030e-22  2.856884e-07 -1.102141e-21
251 -5.464105e-02  3.674932e-05 -1.486859e+00  1.000000e-03  1.463945e-07  1.126193e-22  1.067298e-07 -4.841281e-21
252 -5.427310e-02  3.674932e-05 -1.476846e+00  1.000000e-03  2.060668e-07  2.231966e-22  1.068655e-07  8.145700e-21
253 -5.390514e-02  3.674932e-05 -1.466834e+00  1.000000e-03  3.030748e-07  4.461541e-22  7.839063e-07  4.020802e-19
254 -5.353719e-02  3.674932e-05 -1.456821e+00  1.000000e-03  4.885692e-07  6.765045e-22  1.185961e-07  9.768189e-21
255 -5.316924e-02  3.674932e-05 -1.446809e+00  1.000000e-03  1.029261e-06  1.383589e-21  1.112471e-07  1.205063e-21
256 -5.280128e-02  3.674932e-05 -1.436796e+00  1.000000e-03  1.462866e-05  4.434350e-20  8.947914e-08 -1.704590e-22
257 -5.243333e-02  3.674932e-05 -1.426783e+00  1.000000e-03  2.722523e-05  1.477014e-19  1.009865e-07 -6.018886e-22
258 -5.206538e-02  3.674932e-05 -1.416771e+00  1.000000e-03  8.466523e-06  6.378687e-19  1.160763e-07 -8.707700e-22
259 -5.169742e-02  3.674932e-05 -1.406758e+00  1.000000e-03  1.030819e-05  8.022990e-20  1.339705e-07 -1.041505e-21
260 -5.132947e-02  3.674932e-05 -1.396746e+00  1.000000e-03  4.885840e-05  1.914294e-19  1.544220e-07 -1.353053e-21
261 -5.096152e-02  3.674932e-05 -1.386733e+00  1.000000e-03  1.305059e-03  3.147574e-18  1.765037e-07 -1.746863e-21
262 -5.059356e-02  3.674932e-05 -1.376721e+00  1.000000e-03  3.380551e-05  5.939819e-20  1.921613e-07 -2.657484e-21
263 -5.022561e-02  3.674932e-05 -1.366708e+00  1.000000e-03  1.099959e-05  1.461132e-20  2.382506e-07 -2.875577e-21
264 -4.985766e-02  3.674932e-05 -1.356696e+00  1.000000e-03  7.066768e-06 -2.316106e-21  5.273366e-07 -2.032434e-21
265 -4.948970e-02  3.674932e-05 -1.346683e+00  1.000000e-03  5.136976e-06 -7.980639e-21  5.556519e-06  3.016392e-20
266 -4.912175e-02  3.674932e-05 -1.336671e+00  1.000000e-03  4.151915e-06 -1.109084e-20  1.514058e-05  9.193115e-20
267 -4.875380e-02  3.674932e-05 -1.326658e+00  1.000000e-03  3.845294e-06 -7.206397e-21  4.846527e-05  2.995010e-19
268 -4.838584e-02  3.674932e-05 -1.316646e+00  1.000000e-03  5.732226e-06  5.215705e-20  3.078481e-04  1.681073e-18
269 -4.801789e-02  3.674932e-05 -1.306633e+00  1.000000e-03  4.241849e-04 -9.620746e-18  5.816108e-04  2.386046e-18
270 -4.764994e-02  3.674932e-05 -1.296621e+00  1.000000e-03  1.330062e-05 -5.388718e-19  1.090192e-04 -7.115904e-19
271 -4.728198e-02  3.674932e-05 -1.286608e+00  1.000000e-03  1.135091e-05 -5.400894e-19  4.638622e-05 -1.758757e-18
272 -4.691403e-02  3.674932e-05 -1.276596e+00  1.000000e-03  1.322529e-05 -6.625915e-19  2.627466e-05 -1.220249e-18
273 -4.654608e-02  3.674932e-05 -1.266583e+00  1.000000e-03  1.738622e-05 -8.567474e-19  1.510109e-05 -7.186793e-19
274 -4.617812e-02  3.674932e-05 -1.256571e+00  1.000000e-03  3.211993e-05 -1.207657e-18  2.321346e-06 -7.778145e-20
275 -4.581017e-02  3.674932e-05 -1.246558e+00  1.000000e-03  3.475332e-04 -1.827685e-18  1.416788e-06  1.829227e-20
276 -4.544222e-02  3.674932e-05 -1.236546e+00  1.000000e-03  1.449861e-03  1.930679e-18  1.456570e-06 -2.101039e-20
277 -4.507427e-02  3.674932e-05 -1.226533e+00  1.000000e-03  3.548485e-02  7.341044e-17  1.623100e-06 -1.287755e-20
278 -4.470631e-02  3.674932e-05 -1.216521e+00  1.000000e-03  4.541104e-02  3.261416e-16  1.877045e-06 -1.393965e-20
279 -4.433836e-02  3.674932e-05 -1.206508e+00  1.000000e-03  4.198670e-03 -4.500794e-17  2.233158e-06 -1.498851e-20
280 -4.397041e-02  3.674932e-05 -1.196496e+00  1.000000e-03  2.750126e-03  3.515797e-16  2.737808e-06 -1.007837e-20
281 -4.360245e-02  3.674932e-05 -1.186483e+00  1.000000e-03  8.877257e-04  2.868392e-16  3.514383e-06  1.766858e-21
282 -4.323450e-02  3.674932e-05 -1.176471e+00  1.000000e-03  3.014298e-03 -3.919933e-16  5.313196e-06  3.388132e-20
283 -4.286655e-02  3.674932e-05 -1.166458e+00  1.000000e-03  2.145991e-03 -2.659006e-16  3.303077e-05  1.346782e-19
284 -4.249859e-02  3.674932e-05 -1.156446e+00  1.000000e-03  6.399446e-04 -1.386695e-16  5.148712e-05  1.210410e-18
285 -4.213064e-02  3.674932e-05 -1.146433e+00  1.000000e-03  7.634939e-03  3.226586e-16  1.227250e-04  2.092171e-19
286 -4.176269e-02  3.674932e-05 -1.136421e+00  1.000000e-03  9.287372e-05 -2.564816e-18  5.387708e-04 -3.867202e-18
287 -4.139473e-02  3.674932e-05 -1.126408e+00  1.000000e-03  3.181948e-05 -5.717472e-19  3.360611e-03 -1.119257e-16
288 -4.102678e-02  3.674932e-05 -1.116395e+00  1.000000e-03  2.235658e-05 -4.053053e-19  2.973755e-04  1.806847e-18
289 -4.065883e-02  3.674932e-05 -1.106383e+00  1.000000e-03  2.280180e-05 -4.280693e-19  1.762713e-04  1.300484e-18
290 -4.029087e-02  3.674932e-05 -1.096370e+00  1.000000e-03  3.171286e-05 -8.262277e-19  1.460784e-04  1.947034e-18
291 -3.992292e-02  3.674932e-05 -1.086358e+00  1.000000e-03  7.526752e-05 -3.513493e-18  1.856144e-04  8.672688e-18
292 -3.955497e-02  3.674932e-05 -1.076345e+00  1.000000e-03  2.402236e-03 -1.086452e-15  2.721075e-04  2.011318e-17
293 -3.918701e-02  3.674932e-05 -1.066333e+00  1.000000e-03  1.933767e-03 -7.018176e-17  1.624278e-04  7.108426e-19
294 -3.881906e-02  3.674932e-05 -1.056320e+00  1.000000e-03  4.305343e-03  4.846384e-17  2.181555e-04 -9.244033e-20
295 -3.845111e-02  3.674932e-05 -1.046308e+00  1.000000e-03  8.289160e-03  2.092612e-16  3.407428e-04 -2.916755e-19
296 -3.808315e-02  3.674932e-05 -1.036295e+00  1.000000e-03  1.746463e-02  5.696127e-16  4.824902e-04 -3.533997e-19
297 -3.771520e-02  3.674932e-05 -1.026283e+00  1.000000e-03  3.772649e-02  1.489694e-15  5.800771e-04 -1.131285e-19
298 -3.734725e-02  3.674932e-05 -1.016270e+00  1.000000e-03  3.178349e-02  1.299654e-15  8.520689e-05  1.554736e-19
299 -3.697930e-02  3.674932e-05 -1.006258e+00  1.000000e-03  9.923406e-03  4.098885e-16  1.330876e-05  6.751696e-20
300 -3.661134e-02  3.674932e-05 -9.962453e-01  1.000000e-03  1.197692e-03  1.097439e-17  1.321045e-05  1.703578e-19
301 -3.624339e-02  3.674932e-05 -9.862328e-01  1.000000e-03  1.126955e-03 -1.564214e-18  5.450618e-05  1.596287e-18
302 -3.587544e-02  3.674932e-05 -9.762203e-01  1.000000e-03  7.124784e-04 -3.064969e-18  8.558087e-05  6.393828e-18
303 -3.550748e-02  3.674932e-05 -9.662078e-01  1.000000e-03  8.669562e-05 -5.392754e-19  2.423433e-04 -4.145167e-18
304 -3.513953e-02  3.674932e-05 -9.561952e-01  1.000000e-03  3.638341e-05 -2.560157e-19  2.807318e-04 -4.576731e-18
305 -3.477158e-02  3.674932e-05 -9.461827e-01  1.000000e-03  2.915623e-05 -2.054055e-19  3.042129e-04 -2.372539e-18
306 -3.440362e-02  3.674932e-05 -9.361702e-01  1.000000e-03  1.905839e-05 -1.295551e-19  3.280227e-04 -1.778769e-19
307 -3.403567e-02  3.674932e-05 -9.261577e-01  1.000000e-03  1.449627e-05 -8.094541e-20  3.543543e-04  9.618059e-19
308 -3.366772e-02  3.674932e-05 -9.161452e-01  1.000000e-03  9.372773e-06 -4.399567e-20  3.863144e-04  1.497427e-17
309 -3.329976e-02  3.674932e-05 -9.061327e-01  1.000000e-03  5.527415e-06 -2.497755e-20  5.205911e-04  3.462671e-18
310 -3.293181e-02  3.674932e-05 -8.961202e-01  1.000000e-03  5.343158e-06 -5.637484e-20  8.401358e-04 -1.948049e-17
311 -3.256386e-02  3.674932e-05 -8.861076e-01  1.000000e-03  8.592297e-06 -1.346981e-19  1.421498e-03 -2.838407e-18
312 -3.219590e-02  3.674932e-05 -8.760951e-01  1.000000e-03  1.929719e-05 -1.977160e-19  1.809330e-03 -2.627411e-17
313 -3.182795e-02  3.674932e-05 -8.660826e-01  1.000000e-03  2.650483e-04  1.100465e-17  2.837051e-03 -3.897283e-17
314 -3.146000e-02  3.674932e-05 -8.560701e-01  1.000000e-03  5.096868e-03 -1.244664e-15  4.503881e-04 -4.055806e-18
315 -3.109204e-02  3.674932e-05 -8.460576e-01  1.000000e-03  1.132572e-03 -5.215012e-17  3.504890e-05 -4.290751e-20
316 -3.072409e-02  3.674932e-05 -8.360451e-01  1.000000e-03  9.234237e-04 -2.172809e-17  7.596542e-06  1.424074e-20
317 -3.035614e-02  3.674932e-05 -8.260325e-01  1.000000e-03  8.985320e-04 -1.361351e-17  6.380895e-06  9.714409e-21
318 -2.998818e-02  3.674932e-05 -8.160200e-01  1.000000e-03  8.932542e-04 -1.768605e-17  1.174404e-05  8.851494e-20
319 -2.962023e-02  3.674932e-05 -8.060075e-01  1.000000e-03  8.861326e-04 -1.127570e-17  2.839158e-05 -4.480804e-19
320 -2.925228e-02  3.674932e-05 -7.959950e-01  1.000000e-03  8.723716e-04 -1.666961e-17  5.539631e-05  2.196992e-18
321 -2.888433e-02  3.674932e-05 -7.859825e-01  1.000000e-03  8.548207e-04 -3.268870e-17  9.457444e-04 -3.299447e-17
322 -2.851637e-02  3.674932e-05 -7.759700e-01  1.000000e-03  9.645385e-04  2.181957e-18  1.285687e-03 -2.018649e-17
323 -2.814842e-02  3.674932e-05 -7.659574e-01  1.000000e-03  1.051183e-03  3.722066e-16  1.041417e-03 -6.698337e-18
324 -2.778047e-02  3.674932e-05 -7.559449e-01  1.000000e-03  2.507074e-02 -1.032507e-14  9.629145e-04 -3.693064e-18
325 -2.741251e-02  3.674932e-05 -7.459324e-01  1.000000e-03  2.769938e-03 -7.285839e-16  9.207945e-04 -2.141299e-18
326 -2.704456e-02  3.674932e-05 -7.359199e-01  1.000000e-03  3.446403e-03 -5.590146e-16  8.896283e-04  9.656176e-19
327 -2.667661e-02  3.674932e-05 -7.259074e-01  1.000000e-03  8.515860e-04  1.222438e-16  8.609471e-04 -4.179261e-18
328 -2.630865e-02  3.674932e-05 -7.158949e-01  1.000000e-03  9.095329e-04  2.797242e-17  8.322558e-04 -4.690868e-18
329 -2.594070e-02  3.674932e-05 -7.058824e-01  1.000000e-03  9.035550e-04  2.873136e-18  7.548159e-04 -2.295375e-17
330 -2.557275e-02  3.674932e-05 -6.958698e-01  1.000000e-03  9.079707e-04  3.458605e-17  4.726920e-04  8.504211e-19
331 -2.520479e-02  3.674932e-05 -6.858573e-01  1.000000e-03  9.594777e-04  5.003593e-17  4.573380e-04  1.253609e-18
332 -2.483684e-02  3.674932e-05 -6.758448e-01  1.000000e-03  2.071232e-03  1.024029e-16  5.457591e-04  1.457574e-17
333 -2.446889e-02  3.674932e-05 -6.658323e-01  1.000000e-03  7.360329e-03  7.779151e-17  8.107770e-05  5.895349e-19
334 -2.410093e-02  3.674932e-05 -6.558198e-01  1.000000e-03  3.601717e-03  8.239937e-18  2.003133e-05  2.632155e-19
335 -2.373298e-02  3.674932e-05 -6.458073e-01  1.000000e-03  2.174620e-03 -1.036768e-17  1.391887e-05  4.279634e-19
336 -2.336503e-02  3.674932e-05 -6.357947e-01  1.000000e-03  1.471888e-03 -1.959018e-17  7.335863e-05  5.133020e-19
337 -2.299707e-02  3.674932e-05 -6.257822e-01  1.000000e-03  1.209639e-03 -3.226857e-17  3.139517e-05 -7.411538e-22
338 -2.262912e-02  3.674932e-05 -6.157697e-01  1.000000e-03  1.807273e-03 -8.700722e-17  1.748556e-05 -2.486836e-20
339 -2.226117e-02  3.674932e-05 -6.057572e-01  1.000000e-03  6.574085e-03 -2.586364e-16  2.426399e-05 -1.795710e-19
340 -2.189321e-02  3.674932e-05 -5.957447e-01  1.000000e-03  6.351087e-03  9.627715e-17  4.860648e-05 -9.740879e-21
341 -2.152526e-02  3.674932e-05 -5.857322e-01  1.000000e-03  3.705721e-03 -9.787635e-17  2.258132e-04  1.680513e-18
342 -2.115731e-02  3.674932e-05 -5.757196e-01  1.000000e-03  1.311954e-03 -2.092341e-17  8.855188e-04  6.783040e-18
343 -2.078936e-02  3.674932e-05 -5.657071e-01  1.000000e-03  7.989029e-04  1.125029e-17  1.550644e-02  5.303917e-16
344 -2.042140e-02  3.674932e-05 -5.556946e-01  1.000000e-03  8.950340e-05 -3.635889e-19  4.480590e-04  7.413232e-18
345 -2.005345e-02  3.674932e-05 -5.456821e-01  1.000000e-03  1.048280e-04 -7.595768e-19  1.227017e-04  2.686789e-18
346 -1.968550e-02  3.674932e-05 -5.356696e-01  1.000000e-03  9.074847e-05 -6.005464e-19  5.847040e-05  1.524659e-18
347 -1.931754e-02  3.674932e-05 -5.256571e-01  1.000000e-03  1.782526e-03 -5.645983e-17  3.171321e-05  1.160766e-18
348 -1.894959e-02  3.674932e-05 -5.156446e-01  1.000000e-03  9.087688e-04 -1.155353e-17  1.408246e-05  9.762055e-19
349 -1.858164e-02  3.674932e-05 -5.056320e-01  1.000000e-03  7.378039e-04 -7.704612e-18  1.206454e-05  7.517417e-19
350 -1.821368e-02  3.674932e-05 -4.956195e-01  1.000000e-03  6.580432e-04 -4.553649e-18  1.185217e-05  6.024522e-19
351 -1.784573e-02  3.674932e-05 -4.856070e-01  1.000000e-03  6.056448e-04 -4.438453e-18  1.237695e-05  4.893733e-19
352 -1.747778e-02  3.674932e-05 -4.755945e-01  1.000000e-03  5.646280e-04 -3.361874e-18  1.397044e-05  5.175371e-19
353 -1.710982e-02  3.674932e-05 -4.655820e-01  1.000000e-03  5.290651e-04 -3.455894e-18  5.988842e-05  1.143494e-20
354 -1.674187e-02  3.674932e-05 -4.555695e-01  1.000000e-03  4.962271e-04 -2.619026e-18  9.166604e-05  4.518921e-19
355 -1.637392e-02  3.674932e-05 -4.455569e-01  1.000000e-03  4.646392e-04 -2.425479e-18  8.627911e-05  5.734413e-18
356 -1.600596e-02  3.674932e-05 -4.355444e-01  1.000000e-03  4.333635e-04 -1.721171e-18  1.634675e-04  1.152600e-18
357 -1.563801e-02  3.674932e-05 -4.255319e-01  1.000000e-03  4.016634e-04 -1.799098e-18  2.690201e-04  2.820620e-18
358 -1.527006e-02  3.674932e-05 -4.155194e-01  1.000000e-03  3.687729e-04 -1.346782e-18  1.255373e-04  9.930614e-18
359 -1.490210e-02  3.674932e-05 -4.055069e-01  1.000000e-03  3.335816e-04 -1.280714e-18  1.218616e-04  7.312541e-19
360 -1.453415e-02  3.674932e-05 -3.954944e-01  1.000000e-03  2.937160e-04 -1.026604e-18  1.321394e-04  5.620064e-19
361 -1.416620e-02  3.674932e-05 -3.854819e-01  1.000000e-03  2.390290e-04 -5.302426e-19  1.298648e-04  4.275399e-19
362 -1.379824e-02  3.674932e-05 -3.754693e-01  1.000000e-03  1.988870e-04  6.437450e-20  1.132423e-04  1.355888e-18
363 -1.343029e-02  3.674932e-05 -3.654568e-01  1.000000e-03  8.295743e-04  3.486388e-18  7.560175e-05  2.508806e-18
364 -1.306234e-02  3.674932e-05 -3.554443e-01  1.000000e-03  2.183521e-05  1.917736e-20  1.052568e-05  2.123406e-19
365 -1.269438e-02  3.674932e-05 -3.454318e-01  1.000000e-03  1.013296e-05  3.838118e-22  1.250647e-05 -1.053497e-20
366 -1.232643e-02  3.674932e-05 -3.354193e-01  1.000000e-03  1.309316e-05 -1.755840e-19  3.866878e-05 -1.365841e-19
367 -1.195848e-02  3.674932e-05 -3.254068e-01  1.000000e-03  1.557656e-04 -7.355422e-19  3.273647e-03  1.463673e-18
368 -1.159053e-02  3.674932e-05 -3.153942e-01  1.000000e-03  1.733833e-04 -2.259884e-18  3.092304e-05  2.159934e-19
369 -1.122257e-02  3.674932e-05 -3.053817e-01  1.000000e-03  1.676647e-04 -5.573053e-18  5.968612e-06  3.633771e-19
370 -1.085462e-02  3.674932e-05 -2.953692e-01  1.000000e-03  1.597986e-04 -4.751855e-19  4.880924e-05  1.557185e-17
371 -1.048667e-02  3.674932e-05 -2.853567e-01  1.000000e-03  1.525399e-04 -3.523657e-19  1.347130e-06  1.153024e-19
372 -1.011871e-02  3.674932e-05 -2.753442e-01  1.000000e-03  1.404022e-04 -5.183842e-19  7.701066e-07  5.340278e-20
373 -9.750759e-03  3.674932e-05 -2.653317e-01  1.000000e-03  1.362740e-04 -2.100642e-19  7.135652e-06  1.543718e-19
374 -9.382806e-03  3.674932e-05 -2.553191e-01  1.000000e-03  1.360538e-04 -1.702536e-19  5.095760e-05 -4.179578e-20
375 -9.014853e-03  3.674932e-05 -2.453066e-01  1.000000e-03  1.386807e-04  1.024910e-19  7.408960e-05  4.383396e-20
376 -8.646900e-03  3.674932e-05 -2.352941e-01  1.000000e-03  1.539967e-04  1.627997e-18  8.403537e-05 -1.312901e-19
377 -8.278947e-03  3.674932e-05 -2.252816e-01  1.000000e-03  4.286601e-04  4.090830e-17  8.744512e-05 -3.144610e-20
378 -7.910994e-03  3.674932e-05 -2.152691e-01  1.000000e-03  1.018177e-04  1.129942e-18  8.690984e-05 -1.852885e-20
379 -7.543040e-03  3.674932e-05 -2.052566e-01  1.000000e-03  1.084911e-04  2.837560e-20  8.357675e-05 -4.716915e-20
380 -7.175087e-03  3.674932e-05 -1.952441e-01  1.000000e-03  1.097296e-04 -1.814900e-19  7.786327e-05  6.935082e-20
381 -6.807134e-03  3.674932e-05 -1.852315e-01  1.000000e-03  1.093740e-04 -1.757593e-20  6.958735e-05  6.797439e-20
382 -6.439181e-03  3.674932e-05 -1.752190e-01  1.000000e-03  1.083177e-04 -6.289220e-20  5.736959e-05  9.613824e-20
383 -6.071228e-03  3.674932e-05 -1.652065e-01  1.000000e-03  1.068702e-04 -1.876178e-19  3.206382e-05  8.364450e-20
384 -5.703274e-03  3.674932e-05 -1.551940e-01  1.000000e-03  1.051336e-04 -1.539482e-19  2.590839e-05  6.712736e-20
385 -5.335321e-03  3.674932e-05 -1.451815e-01  1.000000e-03  1.022785e-04  7.835055e-21  2.341005e-05  3.832824e-20
386 -4.967368e-03  3.674932e-05 -1.351690e-01  1.000000e-03  1.004414e-04 -1.060909e-19  2.159492e-05  6.839791e-20
387 -4.599415e-03  3.674932e-05 -1.251564e-01  1.000000e-03  9.700375e-05 -5.442187e-20  2.009474e-05  7.729176e-20
388 -4.231462e-03  3.674932e-05 -1.151439e-01  1.000000e-03  8.274609e-05  2.488159e-20  1.877032e-05  7.856231e-20
389 -3.863508e-03  3.674932e-05 -1.051314e-01  1.000000e-03  4.707047e-05  6.301396e-19  1.754941e-05  9.698527e-20
390 -3.495555e-03  3.674932e-05 -9.511890e-02  1.000000e-03  2.551290e-05  5.669827e-20  1.638695e-05  9.719703e-20
391 -3.127602e-03  3.674932e-05 -8.510638e-02  1.000000e-03  1.964263e-04  3.858235e-19  1.525059e-05  1.725830e-19
392 -2.759649e-03  3.674932e-05 -7.509387e-02  1.000000e-03  1.769201e-06  5.075986e-19  1.411380e-05  4.279634e-19
393 -2.391696e-03  3.674932e-05 -6.508135e-02  1.000000e-03  3.670097e-07  1.123579e-19  1.295461e-05  3.902704e-18
394 -2.023743e-03  3.674932e-05 -5.506884e-02  1.000000e-03  2.748050e-06  3.814554e-20  1.173768e-05  1.873849e-18
395 -1.655789e-03  3.674932e-05 -4.505632e-02  1.000000e-03  7.791056e-06  8.420037e-20  1.045839e-05  2.433102e-19
396 -1.287836e-03  3.674932e-05 -3.504380e-02  1.000000e-03  1.037775e-04  9.641882e-19  9.660836e-06  7.348011e-20
397 -9.198830e-04  3.674932e-05 -2.503129e-02  1.000000e-03  2.722809e-05  2.262670e-19  7.446556e-06  4.817500e-20
398 -5.519298e-04  3.674932e-05 -1.501877e-02  1.000000e-03  4.820818e-06  4.235020e-20  5.017174e-06  6.140989e-21
399 -1.839766e-04  3.674932e-05 -5.006258e-03  1.000000e-03  2.365693e-06  2.204230e-20  7.748268e-07 -5.029258e-22
400  1.839766e-04  3.674932e-05  5.006258e-03  1.000000e-03  1.941406e-06  1.986226e-20  2.666209e-07  5.459392e-22
401  5.519298e-04  3.674932e-05  1.501877e-02  1.000000e-03  3.140833e-06 -1.248112e-20  2.091704e-07  6.815968e-22
402  9.198830e-04  3.674932e-05  2.503129e-02  1.000000e-03  3.549866e-05  2.239873e-19  8.358861e-07  8.134494e-21
403  1.287836e-03  3.674932e-05  3.504380e-02  1.000000e-03  5.126060e-05  1.157259e-19  1.560464e-07  5.260869e-22
404  1.655789e-03  3.674932e-05  4.505632e-02  1.000000e-03  5.636964e-05  9.751467e-20  1.367331e-07  5.293956e-23
405  2.023743e-03  3.674932e-05  5.506884e-02  1.000000e-03  5.743693e-05  4.198107e-20  1.768546e-07 -4.235165e-22
406  2.391696e-03  3.674932e-05  6.508135e-02  1.000000e-03  5.648948e-05 -7.940934e-21  5.961466e-07 -4.542876e-21
407  2.759649e-03  3.674932e-05  7.509387e-02  1.000000e-03  5.442986e-05 -4.023406e-21  3.656520e-06 -3.552244e-20
408  3.127602e-03  3.674932e-05  8.510638e-02  1.000000e-03  5.176081e-05 -1.895236e-20  3.911606e-07 -1.052174e-21
409  3.495555e-03  3.674932e-05  9.511890e-02  1.000000e-03  4.884972e-05 -4.161049e-20  1.938870e-07 -7.974021e-22
410  3.863508e-03  3.674932e-05  1.051314e-01  1.000000e-03  4.611184e-05 -4.510450e-20  1.394030e-07 -5.624828e-23
411  4.231462e-03  3.674932e-05  1.151439e-01  1.000000e-03  4.438478e-05 -7.337423e-20  1.149692e-07  7.113753e-23
412  4.599415e-03  3.674932e-05  1.251564e-01  1.000000e-03  4.702762e-05 -1.129730e-19  1.009993e-07  3.374897e-22
413  4.967368e-03  3.674932e-05  1.351690e-01  1.000000e-03  2.485241e-05  2.196992e-20  9.150108e-08  2.762783e-22
414  5.335321e-03  3.674932e-05  1.451815e-01  1.000000e-03  2.061059e-06  7.179928e-22  8.358920e-08  4.797648e-22
415  5.703274e-03  3.674932e-05  1.551940e-01  1.000000e-03  8.060318e-07  1.675868e-21  1.075292e-07 -4.329463e-21
416  6.071228e-03  3.674932e-05  1.652065e-01  1.000000e-03  4.335940e-07  4.605742e-21  2.560097e-07 -1.311081e-21
417  6.439181e-03  3.674932e-05  1.752190e-01  1.000000e-03  4.688398e-07  1.476228e-20  4.455494e-07  9.365752e-21
418  6.807134e-03  3.674932e-05  1.852315e-01  1.000000e-03  4.585498e-06  4.449570e-19  7.181822e-07 -7.808585e-21
419  7.175087e-03  3.674932e-05  1.952441e-01  1.000000e-03  8.467568e-06  1.208081e-19  1.363060e-05 -2.055378e-19
420  7.543040e-03  3.674932e-05  2.052566e-01  1.000000e-03  9.632561e-06  7.840349e-20  4.935093e-07 -1.013410e-20
421  7.910994e-03  3.674932e-05  2.152691e-01  1.000000e-03  9.979720e-06  5.926584e-20  2.243556e-07 -1.326565e-20
422  8.278947e-03  3.674932e-05  2.252816e-01  1.000000e-03  9.680897e-06  4.327809e-20  1.280631e-06 -2.435782e-19
423  8.646900e-03  3.674932e-05  2.352941e-01  1.000000e-03  8.704418e-06  3.030790e-20  2.484217e-06 -9.383537e-20
424  9.014853e-03  3.674932e-05  2.453066e-01  1.000000e-03  6.789016e-06  2.096407e-20  3.092935e-06 -5.750560e-20
425  9.382806e-03  3.674932e-05  2.553191e-01  1.000000e-03  2.415950e-06  6.471861e-21  3.564189e-06 -5.360130e-20
426  9.750759e-03  3.674932e-05  2.653317e-01  1.000000e-03  3.182887e-07  6.865793e-22  3.991621e-06 -4.224577e-20
427  1.011871e-02  3.674932e-05  2.753442e-01  1.000000e-03  2.087335e-07  3.518620e-22  4.415085e-06 -4.841323e-20
428  1.048667e-02  3.674932e-05  2.853567e-01  1.000000e-03  1.638903e-07  2.686834e-22  4.862826e-06 -5.062345e-20
429  1.085462e-02  3.674932e-05  2.953692e-01  1.000000e-03  1.382955e-07  2.090428e-22  5.365060e-06 -5.526890e-20
430  1.122257e-02  3.674932e-05  3.053817e-01  1.000000e-03  1.243560e-07  1.706675e-22  5.963683e-06 -7.276542e-20
431  1.159053e-02  3.674932e-05  3.153942e-01  1.000000e-03  1.145006e-07  1.351194e-22  6.726171e-06 -8.462389e-20
432  1.195848e-02  3.674932e-05  3.254068e-01  1.000000e-03  1.071616e-07  1.312721e-22  7.773553e-06 -7.705353e-20
433  1.232643e-02  3.674932e-05  3.354193e-01  1.000000e-03  1.017469e-07  1.096927e-22  9.347804e-06 -1.426787e-19
434  1.269438e-02  3.674932e-05  3.454318e-01  1.000000e-03  9.782470e-08  1.126181e-22  1.200410e-05 -1.528349e-19
435  1.306234e-02  3.674932e-05  3.554443e-01  1.000000e-03  9.512742e-08  8.840816e-23  1.729987e-05 -3.433263e-19
436  1.343029e-02  3.674932e-05  3.654568e-01  1.000000e-03  9.350232e-08  1.317132e-22  3.136377e-05 -5.991699e-19
437  1.379824e-02  3.674932e-05  3.754693e-01  1.000000e-03  9.288265e-08  1.148012e-22  1.050256e-04 -2.654601e-18
438  1.416620e-02  3.674932e-05  3.854819e-01  1.000000e-03  9.327696e-08  1.509968e-22  2.617562e-03 -7.287194e-17
439  1.453415e-02  3.674932e-05  3.954944e-01  1.000000e-03  9.477444e-08  1.619319e-22  2.494767e-05 -7.345893e-19
440  1.490210e-02  3.674932e-05  4.055069e-01  1.000000e-03  9.757104e-08  1.973593e-22  1.381869e-06 -1.125495e-19
441  1.527006e-02  3.674932e-05  4.155194e-01  1.000000e-03  1.020324e-07  2.391715e-22  5.939308e-07 -6.275985e-20
442  1.563801e-02  3.674932e-05  4.255319e-01  1.000000e-03  1.088412e-07  2.716077e-22  3.897848e-07 -4.290751e-20
443  1.600596e-02  3.674932e-05  4.355444e-01  1.000000e-03  1.193801e-07  3.696617e-22  1.246282e-06 -1.279020e-19
444  1.637392e-02  3.674932e-05  4.455569e-01  1.000000e-03  1.369489e-07  4.958802e-22  5.024770e-08 -1.818221e-21
445  1.674187e-02  3.674932e-05  4.555695e-01  1.000000e-03  1.723485e-07  8.285022e-22  3.369676e-08 -4.456953e-22
446  1.710982e-02  3.674932e-05  4.655820e-01  1.000000e-03  3.095912e-07  1.890640e-21  2.971382e-08 -1.757694e-22
447  1.747778e-02  3.674932e-05  4.755945e-01  1.000000e-03  2.908826e-06  2.630683e-20  8.871651e-08  1.265690e-21
448  1.784573e-02  3.674932e-05  4.856070e-01  1.000000e-03  4.546801e-06  4.837683e-20  1.252778e-07  3.852594e-22
449  1.821368e-02  3.674932e-05  4.956195e-01  1.000000e-03  5.582189e-06  6.756246e-20  1.300977e-07  9.388500e-22
450  1.858164e-02  3.674932e-05  5.056320e-01  1.000000e-03  6.297736e-06  8.670921e-20  1.286782e-07  6.392038e-22
451  1.894959e-02  3.674932e-05  5.156446e-01  1.000000e-03  6.780581e-06  1.070889e-19  1.158131e-07  3.704735e-22
452  1.931754e-02  3.674932e-05  5.256571e-01  1.000000e-03  7.067896e-06  1.354293e-19  8.570481e-08  1.137373e-22
453  1.968550e-02  3.674932e-05  5.356696e-01  1.000000e-03  7.175727e-06  1.701511e-19  9.029781e-08  6.203855e-24
454  2.005345e-02  3.674932e-05  5.456821e-01  1.000000e-03  7.104649e-06  2.114241e-19  9.683089e-08  4.342698e-23
455  2.042140e-02  3.674932e-05  5.556946e-01  1.000000e-03  6.809323e-06  2.768276e-19  1.043237e-07  3.226004e-23
456  2.078936e-02  3.674932e-05  5.657071e-01  1.000000e-03  6.464214e-06  3.282517e-19  1.127917e-07  2.646978e-23
457  2.115731e-02  3.674932e-05  5.757196e-01  1.000000e-03  5.995214e-06  4.157079e-19  1.224100e-07 -6.203855e-24
458  2.152526e-02  3.674932e-05  5.857322e-01  1.000000e-03  5.401609e-06  5.036140e-19  1.337550e-07 -1.364848e-23
459  2.189321e-02  3.674932e-05  5.957447e-01  1.000000e-03  4.705606e-06  6.213251e-19  1.620033e-07 -1.220091e-23
460  2.226117e-02  3.674932e-05  6.057572e-01  1.000000e-03  3.956781e-06  7.962639e-19  1.863229e-07  1.571643e-23
461  2.262912e-02  3.674932e-05  6.157697e-01  1.000000e-03  3.268877e-06  1.101143e-18  2.076894e-07  2.605619e-23
462  2.299707e-02  3.674932e-05  6.257822e-01  1.000000e-03  2.924139e-06  1.502954e-18  2.294218e-07 -1.675041e-23
463  2.336503e-02  3.674932e-05  6.357947e-01  1.000000e-03  3.723083e-06  2.342867e-18  2.529374e-07 -9.636654e-23
464  2.373298e-02  3.674932e-05  6.458073e-01  1.000000e-03  8.521300e-06  3.758815e-18  2.793027e-07 -3.639595e-23
465  2.410093e-02  3.674932e-05  6.558198e-01  1.000000e-03  3.320900e-05  9.031065e-18  3.095823e-07 -2.812414e-23
466  2.446889e-02  3.674932e-05  6.658323e-01  1.000000e-03  4.082916e-04  5.755081e-17  3.450110e-07 -8.437242e-23
467  2.483684e-02  3.674932e-05  6.758448e-01  1.000000e-03  2.002804e-04  2.224309e-18  3.871433e-07 -1.819797e-22
468  2.520479e-02  3.674932e-05  6.858573e-01  1.000000e-03  4.836627e-06 -1.027557e-19  4.380386e-07 -3.805031e-22
469  2.557275e-02  3.674932e-05  6.958698e-01  1.000000e-03  1.575154e-06 -9.600589e-20  5.005275e-07 -2.316106e-22
470  2.594070e-02  3.674932e-05  7.058824e-01  1.000000e-03  8.255423e-07 -7.182575e-20  5.786318e-07 -2.729696e-22
471  2.630865e-02  3.674932e-05  7.158949e-01  1.000000e-03  5.437580e-07 -5.829307e-20  6.782611e-07 -4.466775e-22
472  2.667661e-02  3.674932e-05  7.259074e-01  1.000000e-03  4.200419e-07 -5.100727e-20  8.084283e-07 -8.007108e-22
473  2.704456e-02  3.674932e-05  7.359199e-01  1.000000e-03  3.748178e-07 -4.908903e-20  9.834926e-07 -9.859993e-22
474  2.741251e-02  3.674932e-05  7.459324e-01  1.000000e-03  3.919601e-07 -5.301793e-20  1.227599e-06 -4.566037e-22
475  2.778047e-02  3.674932e-05  7.559449e-01  1.000000e-03  4.994084e-07 -6.475832e-20  1.588983e-06  5.415717e-20
476  2.814842e-02  3.674932e-05  7.659574e-01  1.000000e-03  8.376688e-07 -9.851390e-20  2.136349e-06 -3.838118e-22
477  2.851637e-02  3.674932e-05  7.759700e-01  1.000000e-03  2.329558e-06 -2.190110e-19  3.069000e-06 -3.824883e-21
478  2.888433e-02  3.674932e-05  7.859825e-01  1.000000e-03  4.397661e-05 -3.300464e-18  4.854929e-06 -9.926167e-21
479  2.925228e-02  3.674932e-05  7.959950e-01  1.000000e-03  5.729439e-06 -2.971233e-19  9.059211e-06 -1.694066e-20
480  2.962023e-02  3.674932e-05  8.060075e-01  1.000000e-03  9.903732e-07 -3.339824e-20  2.393249e-05 -4.838676e-20
481  2.998818e-02  3.674932e-05  8.160200e-01  1.000000e-03  3.609830e-07 -6.926810e-21  2.148235e-04 -5.022905e-19
482  3.035614e-02  3.674932e-05  8.260325e-01  1.000000e-03  1.383442e-07 -7.079632e-22  1.310695e-04 -3.633771e-19
483  3.072409e-02  3.674932e-05  8.360451e-01  1.000000e-03  6.952405e-08 -8.840493e-23  1.389657e-05 -4.166343e-20
484  3.109204e-02  3.674932e-05  8.460576e-01  1.000000e-03  4.339399e-08 -3.703895e-23  3.494397e-06 -1.308269e-20
485  3.146000e-02  3.674932e-05  8.560701e-01  1.000000e-03  2.925063e-08 -1.731829e-23  2.662292e-07 -3.589343e-21
486  3.182795e-02  3.674932e-05  8.660826e-01  1.000000e-03  2.085674e-08 -9.254891e-24  1.327910e-07 -2.003018e-21
487  3.219590e-02  3.674932e-05  8.760951e-01  1.000000e-03  1.555155e-08 -5.177553e-24  1.126191e-07 -1.401244e-21
488  3.256386e-02  3.674932e-05  8.861076e-01  1.000000e-03  1.204028e-08 -3.165945e-24  1.093097e-07 -1.026531e-21
489  3.293181e-02  3.674932e-05  8.961202e-01  1.000000e-03  9.636613e-09 -1.860147e-24  1.124342e-07 -7.494256e-22
490  3.329976e-02  3.674932e-05  9.061327e-01  1.000000e-03  7.953233e-09 -1.156962e-24  1.192561e-07 -6.088049e-22
491  3.366772e-02  3.674932e-05  9.161452e-01  1.000000e-03  6.761431e-09 -5.283980e-25  1.288536e-07 -4.417144e-22
492  3.403567e-02  3.674932e-05  9.261577e-01  1.000000e-03  5.922352e-09 -2.401166e-25  1.407726e-07 -3.441071e-22
493  3.440362e-02  3.674932e-05  9.361702e-01  1.000000e-03  5.351585e-09 -1.645879e-26  1.543227e-07 -2.862045e-22
494  3.477158e-02  3.674932e-05  9.461827e-01  1.000000e-03  5.000613e-09  2.119448e-25  1.674718e-07 -2.051408e-22
495  3.513953e-02  3.674932e-05  9.561952e-01  1.000000e-03  4.847775e-09  4.268179e-25  1.735655e-07 -1.356576e-22
496  3.550748e-02  3.674932e-05  9.662078e-01  1.000000e-03  4.895723e-09  5.450587e-25  1.458690e-07 -9.926167e-24
497  3.587544e-02  3.674932e-05  9.762203e-01  1.000000e-03  5.175028e-09  8.246563e-25  3.138411e-08  1.408792e-24
498  3.624339e-02  3.674932e-05  9.862328e-01  1.000000e-03  5.756465e-09  1.133031e-24  3.455416e-08  7.666769e-24
499  3.661134e-02  3.674932e-05  9.962453e-01  1.000000e-03  6.779547e-09  1.557527e-24  5.058404e-08  2.130717e-23
500  3.697930e-02  3.674932e-05  1.006258e+00  1.000000e-03  8.517627e-09  2.200430e-24  8.827871e-08  5.591709e-23
501  3.734725e-02  3.674932e-05  1.016270e+00  1.000000e-03  1.153709e-08  3.281257e-24  1.941822e-07  1.533806e-22
502  3.771520e-02  3.674932e-05  1.026283e+00  1.000000e-03  1.713711e-08  4.856051e-24  6.610657e-07 -6.261628e-22
503  3.808315e-02  3.674932e-05  1.036295e+00  1.000000e-03  2.880298e-08  8.447501e-24  9.882918e-06  4.308701e-20
504  3.845111e-02  3.674932e-05  1.046308e+00  1.000000e-03  5.852485e-08  1.885713e-23  6.896614e-06  3.881607e-20
505  3.881906e-02  3.674932e-05  1.056320e+00  1.000000e-03  1.713636e-07  6.127599e-23  9.819928e-06  1.258741e-20
506  3.918701e-02  3.674932e-05  1.066333e+00  1.000000e-03  1.623682e-06 -1.080094e-20  7.535612e-06  1.019252e-20
507  3.955497e-02  3.674932e-05  1.076345e+00  1.000000e-03  1.641654e-06 -2.405674e-21  5.168258e-06  6.889019e-21
508  3.992292e-02  3.674932e-05  1.086358e+00  1.000000e-03  2.101608e-07  7.440829e-23  3.754223e-06  4.487662e-21
509  4.029087e-02  3.674932e-05  1.096370e+00  1.000000e-03  8.375466e-08  2.150750e-23  2.867108e-06  3.497966e-21
510  4.065883e-02  3.674932e-05  1.106383e+00  1.000000e-03  4.745193e-08  1.224221e-23  2.275213e-06  2.185127e-21
511  4.102678e-02  3.674932e-05  1.116395e+00  1.000000e-03  3.206314e-08  8.513286e-24  1.859702e-06  2.176855e-21
512  4.139473e-02  3.674932e-05  1.126408e+00  1.000000e-03  2.412795e-08  6.643875e-24  1.555961e-06  1.632350e-21
513  4.176269e-02  3.674932e-05  1.136421e+00  1.000000e-03  1.954964e-08  5.013874e-24  1.326589e-06  1.342630e-21
514  4.213064e-02  3.674932e-05  1.146433e+00  1.000000e-03  1.673537e-08  4.401162e-24  1.148727e-06  1.259673e-21
515  4.249859e-02  3.674932e-05  1.156446e+00  1.000000e-03  1.496327e-08  3.999790e-24  1.007746e-06  1.235633e-21
516  4.286655e-02  3.674932e-05  1.166458e+00  1.000000e-03  1.387295e-08  3.299483e-24  8.939144e-07  9.463011e-22
517  4.323450e-02  3.674932e-05  1.176471e+00  1.000000e-03  1.327610e-08  3.210121e-24  8.005445e-07  1.077697e-21
518  4.360245e-02  3.674932e-05  1.186483e+00  1.000000e-03  1.307772e-08  3.045988e-24  7.229081e-07  6.064591e-22
519  4.397041e-02  3.674932e-05  1.196496e+00  1.000000e-03  1.324204e-08  2.832124e-24  6.575800e-07  8.046884e-22
520  4.433836e-02  3.674932e-05  1.206508e+00  1.000000e-03  1.378018e-08  3.031346e-24  6.020245e-07  8.168893e-22
521  4.470631e-02  3.674932e-05  1.216521e+00  1.000000e-03  1.475202e-08  3.032558e-24  5.543275e-07  6.278511e-22
522  4.507427e-02  3.674932e-05  1.226533e+00  1.000000e-03  1.628267e-08  3.189775e-24  5.130182e-07  8.360082e-22
523  4.544222e-02  3.674932e-05  1.236546e+00  1.000000e-03  1.860246e-08  3.464021e-24  4.769462e-07  9.656946e-22
524  4.581017e-02  3.674932e-05  1.246558e+00  1.000000e-03  2.213530e-08  4.043411e-24  4.451963e-07  5.564870e-22
525  4.617812e-02  3.674932e-05  1.256571e+00  1.000000e-03  2.770350e-08  4.390964e-24  4.170257e-07  6.461500e-22
526  4.654608e-02  3.674932e-05  1.266583e+00  1.000000e-03  3.705585e-08  4.180736e-24  3.918187e-07  7.324175e-22
527  4.691403e-02  3.674932e-05  1.276596e+00  1.000000e-03  5.446238e-08  5.440894e-24  3.690514e-07  1.416812e-21
528  4.728198e-02  3.674932e-05  1.286608e+00  1.000000e-03  9.280885e-08  4.736901e-24  3.482689e-07  1.745918e-21
529  4.764994e-02  3.674932e-05  1.296621e+00  1.000000e-03  2.079843e-07 -1.518813e-23  3.290856e-07  2.622857e-21
530  4.801789e-02  3.674932e-05  1.306633e+00  1.000000e-03  9.511794e-07 -1.934052e-22  3.112803e-07  4.039284e-21
531  4.838584e-02  3.674932e-05  1.316646e+00  1.000000e-03  1.533898e-05 -5.352686e-21  2.953893e-07  5.500053e-21
532  4.875380e-02  3.674932e-05  1.326658e+00  1.000000e-03  4.373779e-07 -1.493384e-22  2.869830e-07  1.334297e-20
533  4.912175e-02  3.674932e-05  1.336671e+00  1.000000e-03  1.292708e-07 -1.156914e-22  3.561074e-07  5.108233e-20
534  4.948970e-02  3.674932e-05  1.346683e+00  1.000000e-03  7.251232e-08 -2.088890e-22  2.055567e-05  4.724194e-18
535  4.985766e-02  3.674932e-05  1.356696e+00  1.000000e-03  1.800555e-07 -1.404966e-21  6.456081e-07  4.583935e-20
536  5.022561e-02  3.674932e-05  1.366708e+00  1.000000e-03  2.549123e-07 -1.671732e-21  4.063832e-07  1.265067e-20
537  5.059356e-02  3.674932e-05  1.376721e+00  1.000000e-03  3.565499e-07 -2.306180e-21  3.450694e-07  4.725075e-21
538  5.096152e-02  3.674932e-05  1.386733e+00  1.000000e-03  6.060130e-07 -3.249993e-21  3.163312e-07  2.741949e-21
539  5.132947e-02  3.674932e-05  1.396746e+00  1.000000e-03  6.047441e-06 -2.829619e-20  2.990417e-07  1.502518e-21
540  5.169742e-02  3.674932e-05  1.406758e+00  1.000000e-03  1.194393e-05 -6.220398e-20  2.872948e-07  1.076567e-21
541  5.206538e-02  3.674932e-05  1.416771e+00  1.000000e-03  1.591692e-05 -7.866818e-20  2.788584e-07  9.411215e-22
542  5.243333e-02  3.674932e-05  1.426783e+00  1.000000e-03  1.892019e-05 -1.308666e-19  2.727543e-07  8.419309e-22
543  5.280128e-02  3.674932e-05  1.436796e+00  1.000000e-03  2.114044e-05 -1.151965e-19  2.685740e-07  7.160411e-22
544  5.316924e-02  3.674932e-05  1.446809e+00  1.000000e-03  2.261010e-05 -1.452662e-19  2.663383e-07  5.937541e-22
545  5.353719e-02  3.674932e-05  1.456821e+00  1.000000e-03  2.329322e-05 -1.653832e-19  2.670614e-07  5.189524e-22
546  5.390514e-02  3.674932e-05  1.466834e+00  1.000000e-03  2.310297e-05 -1.880413e-19  3.017236e-07  3.510477e-22
547  5.427310e-02  3.674932e-05  1.476846e+00  1.000000e-03  2.190196e-05 -1.776652e-19  4.792700e-07  2.949480e-22
548  5.464105e-02  3.674932e-05  1.486859e+00  1.000000e-03  1.950342e-05 -2.346281e-19  6.185053e-07  3.215018e-22
549  5.500900e-02  3.674932e-05  1.496871e+00  1.000000e-03  1.571279e-05 -2.820620e-19  7.749586e-07  2.515081e-22
550  5.537695e-02  3.674932e-05  1.506884e+00  1.000000e-03  1.059172e-05 -4.040347e-19  9.691321e-07  2.352101e-22
551  5.574491e-02  3.674932e-05  1.516896e+00  1.000000e-03  6.010544e-06 -6.111343e-19  1.224557e-06  3.515098e-22
552  5.611286e-02  3.674932e-05  1.526909e+00  1.000000e-03  1.736201e-05 -4.789971e-19  1.577208e-06  4.774254e-22
553  5.648081e-02  3.674932e-05  1.536921e+00  1.000000e-03  2.541447e-04 -1.294266e-18  2.088882e-06  3.280159e-22
554  5.684877e-02  3.674932e-05  1.546934e+00  1.000000e-03  9.792577e-03  8.435093e-17  2.875643e-06  9.382813e-22
555  5.721672e-02  3.674932e-05  1.556946e+00  1.000000e-03  7.701886e-04  6.329030e-18  4.178346e-06  5.882805e-22
556  5.758467e-02  3.674932e-05  1.566959e+00  1.000000e-03  3.594099e-04  1.097755e-18  6.566743e-06  2.386209e-21
557  5.795263e-02  3.674932e-05  1.576971e+00  1.000000e-03  2.505040e-04  2.486889e-18  1.166946e-05  3.646212e-21
558  5.832058e-02  3.674932e-05  1.586984e+00  1.000000e-03  2.036512e-04  1.934623e-18  2.581594e-05  8.994348e-21
559  5.868853e-02  3.674932e-05  1.596996e+00  1.000000e-03  1.790081e-04  6.691560e-19  9.474774e-05 -6.382956e-19
560  5.905649e-02  3.674932e-05  1.607009e+00  1.000000e-03  1.652126e-04  1.727947e-18  5.302247e-03  4.904533e-18
561  5.942444e-02  3.674932e-05  1.617021e+00  1.000000e-03  1.587912e-04  1.069464e-17  1.539132e-04  8.006711e-19
562  5.979239e-02  3.674932e-05  1.627034e+00  1.000000e-03  1.616832e-04  3.864164e-18  3.565563e-05  2.267964e-20
563  6.016035e-02  3.674932e-05  1.637046e+00  1.000000e-03  2.035219e-04 -1.322727e-17  1.592167e-05  1.402071e-20
564  6.052830e-02  3.674932e-05  1.647059e+00  1.000000e-03  9.384758e-05  1.653002e-16  9.372461e-06  1.014785e-20
565  6.089625e-02  3.674932e-05  1.657071e+00  1.000000e-03  8.906825e-05  8.354286e-18  9.243151e-06  4.077670e-20
566  6.126421e-02  3.674932e-05  1.667084e+00  1.000000e-03  9.635458e-05  6.723748e-18  1.605623e-05  1.217610e-19
567  6.163216e-02  3.674932e-05  1.677096e+00  1.000000e-03  9.804809e-05 -3.274418e-18  1.847237e-05  1.545570e-19
568  6.200011e-02  3.674932e-05  1.687109e+00  1.000000e-03  9.789612e-05  8.894375e-19  1.953769e-05  1.608171e-19
569  6.236807e-02  3.674932e-05  1.697121e+00  1.000000e-03  9.696858e-05  4.528979e-19  1.972127e-05  1.829724e-19
570  6.273602e-02  3.674932e-05  1.707134e+00  1.000000e-03  9.565932e-05  2.397103e-19  1.917862e-05  1.598245e-19
571  6.310397e-02  3.674932e-05  1.717146e+00  1.000000e-03  9.413880e-05  1.809474e-19  1.795964e-05  1.462455e-19
572  6.347192e-02  3.674932e-05  1.727159e+00  1.000000e-03  9.249982e-05  1.341488e-19  1.604279e-05  1.275314e-19
573  6.383988e-02  3.674932e-05  1.737171e+00  1.000000e-03  9.143819e-05  3.721651e-20  1.329670e-05  9.560884e-20
574  6.420783e-02  3.674932e-05  1.747184e+00  1.000000e-03  9.022538e-05 -3.102258e-20  9.190958e-06  4.182225e-20
575  6.457578e-02  3.674932e-05  1.757196e+00  1.000000e-03  8.873274e-05  1.164670e-19  3.309640e-06 -1.147730e-19
576  6.494374e-02  3.674932e-05  1.767209e+00  1.000000e-03  8.708961e-05  1.867708e-19  3.065690e-06 -1.823238e-19
577  6.531169e-02  3.674932e-05  1.777222e+00  1.000000e-03  8.533401e-05  9.200895e-20  3.053091e-06 -3.034496e-19
578  6.567964e-02  3.674932e-05  1.787234e+00  1.000000e-03  8.348318e-05  6.373923e-20  3.098682e-06 -5.353248e-19
579  6.604760e-02  3.674932e-05  1.797247e+00  1.000000e-03  8.154635e-05  6.352747e-22  3.230745e-06 -1.024910e-18
580  6.641555e-02  3.674932e-05  1.807259e+00  1.000000e-03  7.952901e-05 -7.824467e-20  3.726161e-06 -2.515688e-18
581  6.678350e-02  3.674932e-05  1.817272e+00  1.000000e-03  7.743465e-05 -5.219841e-20  7.925263e-06 -1.315569e-17
582  6.715146e-02  3.674932e-05  1.827284e+00  1.000000e-03  7.526564e-05  3.483423e-20  5.926969e-05 -2.270048e-19
583  6.751941e-02  3.674932e-05  1.837297e+00  1.000000e-03  7.302374e-05  4.002231e-20  7.873497e-06 -5.246946e-18
584  6.788736e-02  3.674932e-05  1.847309e+00  1.000000e-03  7.071043e-05  7.093901e-21  5.195479e-06 -1.673737e-18
585  6.825532e-02  3.674932e-05  1.857322e+00  1.000000e-03  6.832714e-05 -1.164670e-21  4.489228e-06 -8.300923e-19
586  6.862327e-02  3.674932e-05  1.867334e+00  1.000000e-03  6.587548e-05 -2.027585e-20  4.188536e-06 -4.451158e-19
587  6.899122e-02  3.674932e-05  1.877347e+00  1.000000e-03  6.335771e-05  7.040961e-21  4.024163e-06 -2.860854e-19
588  6.935918e-02  3.674932e-05  1.887359e+00  1.000000e-03  6.077742e-05  2.879912e-20  3.916844e-06 -1.848649e-19
589  6.972713e-02  3.674932e-05  1.897372e+00  1.000000e-03  5.814141e-05 -3.070494e-21  3.835707e-06 -1.306548e-19
590  7.009508e-02  3.674932e-05  1.907384e+00  1.000000e-03  5.546510e-05 -3.758709e-21  3.766444e-06 -9.296187e-20
591  7.046304e-02  3.674932e-05  1.917397e+00  1.000000e-03  5.279738e-05 -2.260519e-20  3.702840e-06 -1.022792e-19
592  7.083099e-02  3.674932e-05  1.927409e+00  1.000000e-03  5.039858e-05  1.815827e-20  3.636496e-06 -7.263308e-20
593  7.119894e-02  3.674932e-05  1.937422e+00  1.000000e-03  5.933642e-05  1.604069e-20  3.568219e-06 -6.013934e-20
594  7.156689e-02  3.674932e-05  1.947434e+00  1.000000e-03  6.682305e-05 -1.005852e-20  3.495526e-06 -3.494011e-20
595  7.193485e-02  3.674932e-05  1.957447e+00  1.000000e-03  7.049800e-05  2.869324e-20  3.417267e-06 -6.839791e-20
596  7.230280e-02  3.674932e-05  1.967459e+00  1.000000e-03  7.250999e-05  5.259545e-20  3.331004e-06 -2.435220e-20
597  7.267075e-02  3.674932e-05  1.977472e+00  1.000000e-03  7.345750e-05  1.228198e-20  3.237268e-06 -2.424632e-20
598  7.303871e-02  3.674932e-05  1.987484e+00  1.000000e-03  7.360552e-05 -3.520481e-21  3.134877e-06 -2.043467e-20
599  7.340666e-02  3.674932e-05  1.997497e+00  1.000000e-03  7.309123e-05  3.017555e-21  3.023306e-06 -1.736418e-20
600  7.377461e-02  3.674932e-05  2.007509e+00  1.000000e-03  7.198398e-05 -1.418780e-19  2.902883e-06 -1.905824e-21
601  7.414257e-02  3.674932e-05  2.017522e+00  1.000000e-03  7.029834e-05 -1.926735e-19  2.784149e-06  8.559268e-19
602  7.451052e-02  3.674932e-05  2.027534e+00  1.000000e-03  6.796753e-05 -2.038173e-20  2.945488e-06  1.852885e-20
603  7.487847e-02  3.674932e-05  2.037547e+00  1.000000e-03  6.467574e-05  1.778769e-20  3.031070e-06 -3.054613e-20
604  7.524643e-02  3.674932e-05  2.047559e+00  1.000000e-03  5.854100e-05  2.853442e-20  3.070353e-06 -4.350308e-20
605  7.561438e-02  3.674932e-05  2.057572e+00  1.000000e-03  5.874266e-05  3.168433e-20  3.076052e-06 -4.548170e-20
606  7.598233e-02  3.674932e-05  2.067584e+00  1.000000e-03  5.929798e-05  4.420453e-21  3.043790e-06 -4.923379e-20
607  7.635029e-02  3.674932e-05  2.077597e+00  1.000000e-03  5.964297e-05 -5.241016e-21  2.959964e-06 -4.293398e-20
608  7.671824e-02  3.674932e-05  2.087610e+00  1.000000e-03  5.967760e-05 -6.379217e-21  2.792726e-06 -3.456953e-20
609  7.708619e-02  3.674932e-05  2.097622e+00  1.000000e-03  5.932727e-05  1.053497e-20  2.421654e-06 -3.959879e-20
610  7.745415e-02  3.674932e-05  2.107635e+00  1.000000e-03  5.946613e-05  3.891058e-21  2.069521e-06 -4.266928e-20
611  7.782210e-02  3.674932e-05  2.117647e+00  1.000000e-03  5.946591e-05  1.950823e-20  2.211801e-06 -3.091670e-20
612  7.819005e-02  3.674932e-05  2.127660e+00  1.000000e-03  5.931214e-05 -2.011703e-21  2.355952e-06 -3.229313e-20
613  7.855801e-02  3.674932e-05  2.137672e+00  1.000000e-03  5.900801e-05  2.398162e-20  2.490552e-06 -3.165786e-20
614  7.892596e-02  3.674932e-05  2.147685e+00  1.000000e-03  5.855783e-05  9.899698e-21  2.612941e-06 -2.731681e-20
615  7.929391e-02  3.674932e-05  2.157697e+00  1.000000e-03  5.796566e-05  1.307607e-20  2.721858e-06 -2.339929e-20
616  7.966186e-02  3.674932e-05  2.167710e+00  1.000000e-03  5.723505e-05  2.797856e-20  2.816482e-06 -2.636390e-20
617  8.002982e-02  3.674932e-05  2.177722e+00  1.000000e-03  5.636898e-05  1.000558e-20  2.896194e-06 -2.339929e-20
618  8.039777e-02  3.674932e-05  2.187735e+00  1.000000e-03  5.536980e-05  3.705769e-22  2.960474e-06 -2.488159e-20
619  8.076572e-02  3.674932e-05  2.197747e+00  1.000000e-03  5.423921e-05  2.901088e-20  3.008863e-06 -2.223461e-20
620  8.113368e-02  3.674932e-05  2.207760e+00  1.000000e-03  5.297823e-05  2.456396e-20  3.040933e-06 -1.545835e-20
621  8.150163e-02  3.674932e-05  2.217772e+00  1.000000e-03  5.158713e-05  2.263166e-20  3.056292e-06 -1.328783e-20
622  8.186958e-02  3.674932e-05  2.227785e+00  1.000000e-03  5.006540e-05  1.815827e-20  3.054671e-06 -1.249374e-20
623  8.223754e-02  3.674932e-05  2.237797e+00  1.000000e-03  4.841198e-05  3.847382e-20  3.036698e-06 -1.434662e-20
624  8.260549e-02  3.674932e-05  2.247810e+00  1.000000e-03  4.665398e-05  1.502160e-20  3.060907e-06 -1.355253e-20
625  8.297344e-02  3.674932e-05  2.257822e+00  1.000000e-03  4.476541e-05  2.564922e-20  3.089258e-06 -1.275843e-20
626  8.334140e-02  3.674932e-05  2.267835e+00  1.000000e-03  4.271792e-05  2.578157e-20  3.061634e-06 -1.707301e-20
627  8.370935e-02  3.674932e-05  2.277847e+00  1.000000e-03  4.051227e-05  1.050056e-19  2.999316e-06 -1.965381e-19
628  8.407730e-02  3.674932e-05  2.287860e+00  1.000000e-03  3.813957e-05  3.456622e-18  2.904738e-06 -5.148372e-20
629  8.444526e-02  3.674932e-05  2.297872e+00  1.000000e-03  3.558320e-05  1.803254e-19  2.777704e-06 -1.468411e-20
630  8.481321e-02  3.674932e-05  2.307885e+00  1.000000e-03  3.282527e-05  3.822236e-20  2.614724e-06 -1.011973e-20
631  8.518116e-02  3.674932e-05  2.317897e+00  1.000000e-03  2.984820e-05  2.144714e-20  2.409056e-06 -6.074814e-21
632  8.554912e-02  3.674932e-05  2.327910e+00  1.000000e-03  2.675564e-05  4.916762e-21  2.147712e-06 -4.830735e-21
633  8.591707e-02  3.674932e-05  2.337922e+00  1.000000e-03  2.945357e-05  8.486873e-21  1.801584e-06 -2.355810e-21
634  8.628502e-02  3.674932e-05  2.347935e+00  1.000000e-03  3.005177e-05  2.421985e-21  1.272680e-06 -1.535247e-21
635  8.665298e-02  3.674932e-05  2.357947e+00  1.000000e-03  2.769001e-05 -3.030790e-21  5.671924e-07  3.441071e-22
636  8.702093e-02  3.674932e-05  2.367960e+00  1.000000e-03  2.037436e-05 -3.096964e-21  5.442254e-07 -3.970467e-22
637  8.738888e-02  3.674932e-05  2.377972e+00  1.000000e-03  2.225097e-05  3.427836e-21  5.428534e-07  5.823352e-22
638  8.775684e-02  3.674932e-05  2.387985e+00  1.000000e-03  2.443307e-05  1.482308e-21  5.445999e-07  1.852885e-21
639  8.812479e-02  3.674932e-05  2.397997e+00  1.000000e-03  2.643006e-05 -3.705769e-21  5.467553e-07  1.641126e-21
640  8.849274e-02  3.674932e-05  2.408010e+00  1.000000e-03  2.822887e-05 -7.993873e-21  5.485159e-07  3.970467e-22
641  8.886069e-02  3.674932e-05  2.418023e+00  1.000000e-03  2.984076e-05  9.582060e-21  5.495887e-07  3.176374e-22
642  8.922865e-02  3.674932e-05  2.428035e+00  1.000000e-03  3.127873e-05 -4.870439e-21  5.498768e-07  1.429368e-21
643  8.959660e-02  3.674932e-05  2.438048e+00  1.000000e-03  3.255449e-05 -7.596827e-21  5.493852e-07  1.058791e-22
644  8.996455e-02  3.674932e-05  2.448060e+00  1.000000e-03  3.367821e-05 -5.082198e-21  5.481934e-07 -9.529121e-22
645  9.033251e-02  3.674932e-05  2.458073e+00  1.000000e-03  3.465871e-05 -1.469073e-20  5.464619e-07  6.882143e-22
646  9.070046e-02  3.674932e-05  2.468085e+00  1.000000e-03  3.550372e-05 -9.529121e-21  5.444667e-07  5.293956e-22
647  9.106841e-02  3.674932e-05  2.478098e+00  1.000000e-03  3.622004e-05 -7.146840e-22  5.426818e-07  2.911676e-22
648  9.143637e-02  3.674932e-05  2.488110e+00  1.000000e-03  3.681375e-05  0.000000e+00  5.419635e-07 -1.058791e-22
649  9.180432e-02  3.674932e-05  2.498123e+00  1.000000e-03  3.729031e-05  4.235165e-22  5.440077e-07  3.441071e-22
650  9.217227e-02  3.674932e-05  2.508135e+00  1.000000e-03  3.765465e-05  1.037615e-20  5.526879e-07 -1.098496e-21
651  9.254023e-02  3.674932e-05  2.518148e+00  1.000000e-03  3.791126e-05  4.499863e-21  5.793060e-07  5.293956e-22
652  9.290818e-02  3.674932e-05  2.528160e+00  1.000000e-03  3.806419e-05 -1.389663e-20  6.804233e-07  1.085261e-21
653  9.327613e-02  3.674932e-05  2.538173e+00  1.000000e-03  3.811712e-05 -6.617445e-21  3.543882e-06 -1.091878e-21
654  9.364409e-02  3.674932e-05  2.548185e+00  1.000000e-03  3.807336e-05  7.252720e-21  6.600562e-06 -1.315217e-21
655  9.401204e-02  3.674932e-05  2.558198e+00  1.000000e-03  3.793578e-05  1.402898e-21  7.642944e-06 -1.452529e-21
656  9.437999e-02  3.674932e-05  2.568210e+00  1.000000e-03  3.770670e-05 -2.752857e-21  7.641166e-06 -1.058791e-21
657  9.474795e-02  3.674932e-05  2.578223e+00  1.000000e-03  3.738948e-05 -1.471720e-20  6.679659e-06 -2.607273e-21
658  9.511590e-02  3.674932e-05  2.588235e+00  1.000000e-03  3.697920e-05 -4.685151e-21  4.277776e-06 -1.985233e-22
659  9.548385e-02  3.674932e-05  2.598248e+00  1.000000e-03  3.645883e-05 -5.082198e-21  3.726654e-07  4.181398e-22
660  9.585181e-02  3.674932e-05  2.608260e+00  1.000000e-03  3.591986e-05  1.950823e-20  1.898171e-07  8.952162e-22
661  9.621976e-02  3.674932e-05  2.618273e+00  1.000000e-03  3.530989e-05  5.066316e-20  1.364732e-07  3.060775e-21
662  9.658771e-02  3.674932e-05  2.628285e+00  1.000000e-03  3.464540e-05  1.314516e-18  1.103926e-07  1.210496e-20
663  9.695566e-02  3.674932e-05  2.638298e+00  1.000000e-03  3.387004e-05  1.677125e-19  8.790689e-08 -1.716400e-23
664  9.732362e-02  3.674932e-05  2.648310e+00  1.000000e-03  3.304058e-05  3.245195e-20  7.542827e-08 -1.695720e-23
665  9.769157e-02  3.674932e-05  2.658323e+00  1.000000e-03  3.214308e-05  1.275843e-20  6.603449e-08 -3.515518e-24
666  9.805952e-02  3.674932e-05  2.668335e+00  1.000000e-03  3.118171e-05 -4.685151e-21  5.861947e-08 -4.601192e-24
667  9.842748e-02  3.674932e-05  2.678348e+00  1.000000e-03  3.015483e-05 -7.782115e-21  5.257340e-08  1.519944e-23
668  9.879543e-02  3.674932e-05  2.688360e+00  1.000000e-03  2.906355e-05  3.335192e-21  4.752346e-08  7.392927e-24
669  9.916338e-02  3.674932e-05  2.698373e+00  1.000000e-03  2.790885e-05 -2.938146e-21  4.322641e-08  2.326445e-24
670  9.953134e-02  3.674932e-05  2.708385e+00  1.000000e-03  2.669156e-05 -7.371834e-21  3.951563e-08 -5.014782e-24
671  9.989929e-02  3.674932e-05  2.718398e+00  1.000000e-03  2.541227e-05 -3.573420e-21  3.627266e-08  1.550964e-25
672  1.002672e-01  3.674932e-05  2.728411e+00  1.000000e-03  2.407127e-05 -5.293956e-22  3.341090e-08  8.271806e-25
673  1.006352e-01  3.674932e-05  2.738423e+00  1.000000e-03  2.266850e-05  2.091113e-21  3.086683e-08  1.364848e-23
674  1.010031e-01  3.674932e-05  2.748436e+00  1.000000e-03  2.120343e-05 -3.321957e-21  2.860033e-08  2.680582e-23
675  1.013711e-01  3.674932e-05  2.758448e+00  1.000000e-03  1.967607e-05 -8.417390e-21  2.668657e-08  8.002972e-23
676  1.017391e-01  3.674932e-05  2.768461e+00  1.000000e-03  1.808006e-05 -4.301339e-21  2.588024e-08  5.307139e-22
677  1.021070e-01  3.674932e-05  2.778473e+00  1.000000e-03  1.641535e-05 -4.698386e-21  2.294740e-08  2.414333e-23
678  1.024750e-01  3.674932e-05  2.788486e+00  1.000000e-03  1.467427e-05  2.276401e-21  2.134073e-08  9.098987e-24
679  1.028429e-01  3.674932e-05  2.798498e+00  1.000000e-03  1.284521e-05 -5.618211e-21  1.988681e-08 -8.271806e-25
680  1.032109e-01  3.674932e-05  2.808511e+00  1.000000e-03  1.090714e-05  4.963084e-22  1.854673e-08  1.589738e-24
681  1.035788e-01  3.674932e-05  2.818523e+00  1.000000e-03  8.817276e-06 -1.522012e-21  1.730391e-08  4.161752e-24
682  1.039468e-01  3.674932e-05  2.828536e+00  1.000000e-03  6.469987e-06 -7.841672e-22  1.614612e-08 -1.757759e-24
683  1.043147e-01  3.674932e-05  2.838548e+00  1.000000e-03  3.456269e-06 -9.065900e-22  1.506327e-08  1.305394e-24
684  1.046827e-01  3.674932e-05  2.848561e+00  1.000000e-03  2.546859e-07 -7.413606e-23  1.404706e-08  4.911385e-24
685  1.050506e-01  3.674932e-05  2.858573e+00  1.000000e-03  1.151403e-07 -4.347868e-23  1.309396e-08  2.714186e-25
686  1.054186e-01  3.674932e-05  2.868586e+00  1.000000e-03  7.104586e-08 -3.567216e-24  1.239953e-08  2.869283e-24
687  1.057865e-01  3.674932e-05  2.878598e+00  1.000000e-03  4.720981e-08  1.274375e-23  1.215647e-08  1.242063e-23
688  1.061545e-01  3.674932e-05  2.888611e+00  1.000000e-03  3.113387e-08  2.102202e-23  1.148966e-08  3.083833e-23
689  1.065224e-01  3.674932e-05  2.898623e+00  1.000000e-03  1.811602e-08  2.338724e-23  1.066136e-08  4.023458e-23
690  1.068904e-01  3.674932e-05  2.908636e+00  1.000000e-03  2.519707e-09  1.283947e-23  9.558444e-09  4.867441e-23
691  1.072584e-01  3.674932e-05  2.918648e+00  1.000000e-03  4.695473e-10  8.753611e-25  8.824775e-09  4.919140e-23
692  1.076263e-01  3.674932e-05  2.928661e+00  1.000000e-03  2.537167e-10  4.301911e-25  8.166028e-09  4.882304e-23
693  1.079943e-01  3.674932e-05  2.938673e+00  1.000000e-03  1.640624e-10  2.818545e-25  7.461891e-09  3.895504e-23
694  1.083622e-01  3.674932e-05  2.948686e+00  1.000000e-03  1.153528e-10  2.056900e-25  6.545025e-09  9.234696e-24
695  1.087302e-01  3.674932e-05  2.958698e+00  1.000000e-03  8.534227e-11  1.608298e-25  6.034599e-09  8.853417e-25
696  1.090981e-01  3.674932e-05  2.968711e+00  1.000000e-03  6.540795e-11  1.291434e-25  5.637427e-09 -3.554292e-25
697  1.094661e-01  3.674932e-05  2.978723e+00  1.000000e-03  5.147532e-11  1.093030e-25  5.264142e-09 -2.488004e-25
698  1.098340e-01  3.674932e-05  2.988736e+00  1.000000e-03  4.136931e-11  9.435329e-26  4.909795e-09  1.938705e-25
699  1.102020e-01  3.674932e-05  2.998748e+00  1.000000e-03  3.382768e-11  8.280722e-26  4.571948e-09  1.486340e-25
700  1.105699e-01  3.674932e-05  3.008761e+00  1.000000e-03  2.807171e-11  7.439787e-26  4.248731e-09  5.040632e-25
701  1.109379e-01  3.674932e-05  3.018773e+00  1.000000e-03  2.359792e-11  6.794203e-26  3.938475e-09  8.368741e-25
702  1.113058e-01  3.674932e-05  3.028786e+00  1.000000e-03  2.006844e-11  6.331343e-26  3.639592e-09  5.299126e-25
703  1.116738e-01  3.674932e-05  3.038798e+00  1.000000e-03  1.724967e-11  5.932988e-26  3.350495e-09  3.554292e-26
704  1.120417e-01  3.674932e-05  3.048811e+00  1.000000e-03  1.497601e-11  5.680730e-26  3.069529e-09  1.260158e-25
705  1.124097e-01  3.674932e-05  3.058824e+00  1.000000e-03  1.312771e-11  5.551249e-26  2.794888e-09  5.961517e-25
706  1.127777e-01  3.674932e-05  3.068836e+00  1.000000e-03  1.161675e-11  5.452067e-26  2.524488e-09  3.812786e-25
707  1.131456e-01  3.674932e-05  3.078849e+00  1.000000e-03  1.037769e-11  5.444743e-26  2.255760e-09  2.689953e-25
708  1.135136e-01  3.674932e-05  3.088861e+00  1.000000e-03  9.361512e-12  5.561442e-26  1.985277e-09  4.660969e-25
709  1.138815e-01  3.674932e-05  3.098874e+00  1.000000e-03  8.531609e-12  5.709107e-26  1.707992e-09  4.362085e-26
710  1.142495e-01  3.674932e-05  3.108886e+00  1.000000e-03  7.861056e-12  5.997118e-26  1.415448e-09  1.712522e-25
711  1.146174e-01  3.674932e-05  3.118899e+00  1.000000e-03  7.330989e-12  6.419762e-26  1.090431e-09  1.657996e-25
712  1.149854e-01  3.674932e-05  3.128911e+00  1.000000e-03  6.929868e-12  6.993153e-26  6.821414e-10  2.244656e-25
713  1.153533e-01  3.674932e-05  3.138924e+00  1.000000e-03  6.653764e-12  7.814292e-26  6.461532e-11  2.048049e-25
714  1.157213e-01  3.674932e-05  3.148936e+00  1.000000e-03  6.508067e-12  8.923570e-26  2.715370e-11  2.490785e-25
715  1.160892e-01  3.674932e-05  3.158949e+00  1.000000e-03  6.511792e-12  1.050394e-25  1.972204e-11  3.198980e-25
716  1.164572e-01  3.674932e-05  3.168961e+00  1.000000e-03  6.707380e-12  1.288015e-25  1.728053e-11  4.520914e-25
717  1.168251e-01  3.674932e-05  3.178974e+00  1.000000e-03  7.184015e-12  1.653669e-25  1.796458e-11  7.157694e-25
718  1.171931e-01  3.674932e-05  3.188986e+00  1.000000e-03  8.139425e-12  2.228429e-25  2.426308e-11  1.376499e-24
719  1.175610e-01  3.674932e-05  3.198999e+00  1.000000e-03  1.007433e-11  3.295577e-25  5.779362e-11  4.176669e-24
720  1.179290e-01  3.674932e-05  3.209011e+00  1.000000e-03  1.458779e-11  5.545447e-25  1.440692e-09  1.071257e-22
721  1.182969e-01  3.674932e-05  3.219024e+00  1.000000e-03  2.947165e-11  1.204336e-24  1.315603e-10  7.209739e-24
722  1.186649e-01  3.674932e-05  3.229036e+00  1.000000e-03  1.643038e-10  5.829298e-24  3.814975e-11  1.088881e-24
723  1.190329e-01  3.674932e-05  3.239049e+00  1.000000e-03  4.713177e-10  6.419514e-24  2.729174e-11  2.113689e-25
724  1.194008e-01  3.674932e-05  3.249061e+00  1.000000e-03  5.762966e-11 -9.834811e-25  3.027215e-11 -1.552672e-25
725  1.197688e-01  3.674932e-05  3.259074e+00  1.000000e-03  4.676878e-11 -1.542572e-24  4.993282e-11 -5.540681e-25
726  1.201367e-01  3.674932e-05  3.269086e+00  1.000000e-03  9.261972e-11 -3.017177e-24  1.562091e-10 -1.974774e-24
727  1.205047e-01  3.674932e-05  3.279099e+00  1.000000e-03  7.962439e-10 -1.993506e-23  8.679693e-09 -9.194913e-23
728  1.208726e-01  3.674932e-05  3.289111e+00  1.000000e-03  6.228309e-10 -1.017160e-23  2.569537e-10 -2.001835e-24
729  1.212406e-01  3.674932e-05  3.299124e+00  1.000000e-03  8.916435e-11 -8.526522e-25  6.809849e-11 -3.379646e-25
730  1.216085e-01  3.674932e-05  3.309136e+00  1.000000e-03  4.001840e-11 -1.827192e-25  3.664703e-11 -1.067328e-25
731  1.219765e-01  3.674932e-05  3.319149e+00  1.000000e-03  2.637131e-11 -4.347240e-26  2.646977e-11 -4.335402e-26
732  1.223444e-01  3.674932e-05  3.329161e+00  1.000000e-03  2.123156e-11  1.671606e-27  2.250137e-11 -2.398808e-26
733  1.227124e-01  3.674932e-05  3.339174e+00  1.000000e-03  1.936121e-11  1.880574e-26  2.121719e-11 -2.039464e-26
734  1.230803e-01  3.674932e-05  3.349186e+00  1.000000e-03  1.926076e-11  2.653984e-26  2.152869e-11 -2.653611e-26
735  1.234483e-01  3.674932e-05  3.359199e+00  1.000000e-03  2.050554e-11  3.099939e-26  2.314394e-11 -4.027397e-26
736  1.238162e-01  3.674932e-05  3.369212e+00  1.000000e-03  2.317630e-11  3.568787e-26  2.619401e-11 -6.500804e-26
737  1.241842e-01  3.674932e-05  3.379224e+00  1.000000e-03  2.783274e-11  4.282013e-26  3.124924e-11 -1.084401e-25
738  1.245522e-01  3.674932e-05  3.389237e+00  1.000000e-03  3.589293e-11  5.658937e-26  3.968410e-11 -1.892964e-25
739  1.249201e-01  3.674932e-05  3.399249e+00  1.000000e-03  5.099101e-11  8.454550e-26  5.491391e-11 -3.536981e-25
740  1.252881e-01  3.674932e-05  3.409262e+00  1.000000e-03  8.436903e-11  1.485787e-25  8.711658e-11 -7.654339e-25
741  1.256560e-01  3.674932e-05  3.419274e+00  1.000000e-03  1.863021e-10  3.568302e-25  1.789907e-10 -2.113600e-24
742  1.260240e-01  3.674932e-05  3.429287e+00  1.000000e-03  8.909078e-10  1.649429e-24  7.114729e-10 -1.100265e-23
743  1.263919e-01  3.674932e-05  3.439299e+00  1.000000e-03  9.291679e-09  7.627878e-24  3.067986e-08 -6.049551e-22
744  1.267599e-01  3.674932e-05  3.449312e+00  1.000000e-03  3.904914e-10 -8.994207e-25  5.054969e-10 -1.158615e-23
745  1.271278e-01  3.674932e-05  3.459324e+00  1.000000e-03  1.520057e-10 -1.321419e-24  1.712195e-10 -4.054437e-24
746  1.274958e-01  3.674932e-05  3.469337e+00  1.000000e-03  1.234138e-10 -2.257872e-24  1.128313e-10 -2.479110e-24
747  1.278637e-01  3.674932e-05  3.479349e+00  1.000000e-03  1.982078e-10 -5.479317e-24  1.055850e-10 -1.959762e-24
748  1.282317e-01  3.674932e-05  3.489362e+00  1.000000e-03  9.620125e-10 -2.927442e-23  1.260677e-10 -1.901416e-24
749  1.285996e-01  3.674932e-05  3.499374e+00  1.000000e-03  1.303906e-08 -3.506335e-22  1.921939e-10 -2.229819e-24
750  1.289676e-01  3.674932e-05  3.509387e+00  1.000000e-03  7.676050e-10 -1.602489e-23  4.317026e-10 -2.583485e-24
751  1.293355e-01  3.674932e-05  3.519399e+00  1.000000e-03  3.915481e-10 -6.222730e-24  2.497480e-09  5.059511e-25
752  1.297035e-01  3.674932e-05  3.529412e+00  1.000000e-03  3.328642e-10 -4.373515e-24  2.875034e-08  1.796640e-22
753  1.300715e-01  3.674932e-05  3.539424e+00  1.000000e-03  3.674758e-10 -4.366952e-24  2.286519e-09  3.066425e-23
754  1.304394e-01  3.674932e-05  3.549437e+00  1.000000e-03  5.013331e-10 -6.000626e-24  1.518282e-09  2.719243e-23
755  1.308074e-01  3.674932e-05  3.559449e+00  1.000000e-03  9.051672e-10 -1.275567e-23  1.933642e-09  4.246474e-23
756  1.311753e-01  3.674932e-05  3.569462e+00  1.000000e-03  3.041848e-09 -5.342153e-23  4.818846e-09  1.371507e-22
757  1.315433e-01  3.674932e-05  3.579474e+00  1.000000e-03  2.624562e-07 -5.782613e-21  1.938011e-07  6.617801e-21
758  1.319112e-01  3.674932e-05  3.589487e+00  1.000000e-03  2.791507e-09 -6.865899e-23  6.775196e-09  2.350874e-22
759  1.322792e-01  3.674932e-05  3.599499e+00  1.000000e-03  9.953545e-10 -2.244314e-23  2.271780e-09  5.741077e-23
760  1.326471e-01  3.674932e-05  3.609512e+00  1.000000e-03  7.655837e-10 -1.341885e-23  1.943860e-09  2.900775e-23
761  1.330151e-01  3.674932e-05  3.619524e+00  1.000000e-03  8.275581e-10 -1.050661e-23  2.396218e-09  1.918292e-23
762  1.333830e-01  3.674932e-05  3.629537e+00  1.000000e-03  1.060165e-09 -9.787099e-24  3.499396e-09  1.502298e-23
763  1.337510e-01  3.674932e-05  3.639549e+00  1.000000e-03  1.503902e-09 -1.032946e-23  5.690400e-09  1.297651e-23
764  1.341189e-01  3.674932e-05  3.649562e+00  1.000000e-03  2.314114e-09 -1.234072e-23  1.020410e-08  1.199336e-23
765  1.344869e-01  3.674932e-05  3.659574e+00  1.000000e-03  3.881040e-09 -1.648472e-23  2.054497e-08  1.052554e-23
766  1.348548e-01  3.674932e-05  3.669587e+00  1.000000e-03  7.289639e-09 -2.536924e-23  4.868694e-08  6.814458e-24
767  1.352228e-01  3.674932e-05  3.679599e+00  1.000000e-03  1.636717e-08 -4.831601e-23  1.509845e-07 -2.121709e-23
768  1.355907e-01  3.674932e-05  3.689612e+00  1.000000e-03  5.227449e-08 -1.343003e-22  8.269503e-07 -1.868824e-22
769  1.359587e-01  3.674932e-05  3.699625e+00  1.000000e-03  5.204287e-07 -1.197995e-21  5.688304e-05 -2.697602e-20
770  1.363267e-01  3.674932e-05  3.709637e+00  1.000000e-03  6.458411e-07 -1.324443e-21  3.101231e-05 -1.331659e-20
771  1.366946e-01  3.674932e-05  3.719650e+00  1.000000e-03  8.632297e-08 -1.545684e-22  3.045666e-05 -1.307342e-20
772  1.370626e-01  3.674932e-05  3.729662e+00  1.000000e-03  3.702972e-08 -6.182842e-23  2.104830e-05 -9.999926e-21
773  1.374305e-01  3.674932e-05  3.739675e+00  1.000000e-03  2.271001e-08 -3.554679e-23  1.489760e-05 -8.909320e-21
774  1.377985e-01  3.674932e-05  3.749687e+00  1.000000e-03  1.667577e-08 -2.455133e-23  1.063399e-05 -7.035196e-21
775  1.381664e-01  3.674932e-05  3.759700e+00  1.000000e-03  1.369605e-08 -1.935310e-23  7.335843e-06 -4.436391e-21
776  1.385344e-01  3.674932e-05  3.769712e+00  1.000000e-03  1.217578e-08 -1.654883e-23  4.250848e-06 -2.606491e-21
777  1.389023e-01  3.674932e-05  3.779725e+00  1.000000e-03  1.151028e-08 -1.518213e-23  2.571681e-07 -1.577984e-22
778  1.392703e-01  3.674932e-05  3.789737e+00  1.000000e-03  1.145558e-08 -1.470505e-23  5.571531e-08 -3.053818e-23
779  1.396382e-01  3.674932e-05  3.799750e+00  1.000000e-03  1.193857e-08 -1.499423e-23  2.502531e-08 -1.145753e-23
780  1.400062e-01  3.674932e-05  3.809762e+00  1.000000e-03  1.300077e-08 -1.620618e-23  1.400739e-08 -4.138332e-24
781  1.403741e-01  3.674932e-05  3.819775e+00  1.000000e-03  1.480360e-08 -1.830546e-23  8.879886e-09  4.788995e-25
782  1.407421e-01  3.674932e-05  3.829787e+00  1.000000e-03  1.769251e-08 -2.201223e-23  6.191191e-09  6.064529e-24
783  1.411100e-01  3.674932e-05  3.839800e+00  1.000000e-03  2.237492e-08 -2.814109e-23  5.119885e-09  3.330776e-23
784  1.414780e-01  3.674932e-05  3.849812e+00  1.000000e-03  3.041820e-08 -3.838269e-23  1.156713e-08  3.823049e-22
785  1.418460e-01  3.674932e-05  3.859825e+00  1.000000e-03  4.594176e-08 -5.912364e-23  1.262449e-08  4.363928e-22
786  1.422139e-01  3.674932e-05  3.869837e+00  1.000000e-03  8.409808e-08 -1.111886e-22  1.231210e-08  4.257914e-22
787  1.425819e-01  3.674932e-05  3.879850e+00  1.000000e-03  3.031398e-07 -4.110469e-22  1.155114e-08  3.893799e-22
788  1.429498e-01  3.674932e-05  3.889862e+00  1.000000e-03  3.807021e-06 -5.213947e-21  1.055913e-08  3.546546e-22
789  1.433178e-01  3.674932e-05  3.899875e+00  1.000000e-03  5.809968e-06 -8.171151e-21  9.307613e-09  3.054008e-22
790  1.436857e-01  3.674932e-05  3.909887e+00  1.000000e-03  7.176932e-06 -1.019868e-20  7.001269e-09  1.975998e-22
791  1.440537e-01  3.674932e-05  3.919900e+00  1.000000e-03  8.112722e-06 -1.218669e-20  3.472513e-09 -7.225933e-25
792  1.444216e-01  3.674932e-05  3.929912e+00  1.000000e-03  8.455338e-06 -1.363785e-20  6.168728e-09 -2.162496e-23
793  1.447896e-01  3.674932e-05  3.939925e+00  1.000000e-03  7.241001e-06 -1.377718e-20  8.057683e-08 -4.233631e-22
794  1.451575e-01  3.674932e-05  3.949937e+00  1.000000e-03  6.777021e-06 -4.819164e-21  2.203196e-07 -1.054760e-21
795  1.455255e-01  3.674932e-05  3.959950e+00  1.000000e-03  1.988081e-06  3.694937e-21  2.909059e-07 -1.341085e-21
796  1.458934e-01  3.674932e-05  3.969962e+00  1.000000e-03  1.370780e-06  4.299847e-21  3.319833e-07 -1.515685e-21
797  1.462614e-01  3.674932e-05  3.979975e+00  1.000000e-03  1.323759e-06  5.660296e-21  3.464341e-07 -1.591605e-21
798  1.466293e-01  3.674932e-05  3.989987e+00  1.000000e-03  1.686891e-06  9.339370e-21  3.260666e-07 -1.529987e-21
799  1.469973e-01  3.674932e-05  4.000000e+00  1.000000e-03  6.734467e-06  4.725014e-20  2.411280e-07 -1.169342e-21
横坐标不要当真,没有按照本正能量画。

这就是修改之后的transmission,跟修改之前的变化还是挺大滴,下面把原来的transmission图粘贴出来,进行对比一下。

为了验证是否是晶格矢量不匹配的问题,必须要对原来的NEGF-8ZGNR-0.3这个例子的晶格矢量修改,修改的内容是把4.88都改为6.1。仅改此一处,发现仿真结果如下:

从结果来看,确实变化非常大,而起没什么规律,这也就是说,必须将晶格匹配做到很好才行。

如果把原子的Atoms.UnitVectors修改为4.88的两倍,即9.76,transmission就全部成为0。可见,这个单位矢量对于这个非平衡格林函数的透射谱的影响之大。4.88就正好,6.1就很小,而到了9.76就到了全0.

下面是conductance文件:

#/************************************************************/
#
#  Conductance:
#
#    G0 = e^2/h
#       = (1.60217653 * 10^{-19})^2/(6.626076 * 10^{-34})
#       = 3.87404194169 * 10^{-4} [C^2 J^{-1} s^{-1}]
#    note that
#    e = 1.60217653 * 10^{-19} C
#    h = 6.626076 * 10^{-34} J s
#***********************************************************/
#
# SpinP_switch= 1
# Chemical potential (eV)      Left, Right Leads=-7.092117e+00 -7.382248e+00
# Chemical potential (Hartree) Left, Right Leads=-2.606305e-01 -2.712926e-01
# diff Chemical potential (eV)     = 2.901306e-01
# diff Chemical potential (Hartree)= 1.066210e-02
# Corresponding bias voltage (Volt)= 2.901305e-01
# average conductance for up and down spins (G0)= 1.022572e-04  2.338899e-05
#
# i2 i3 k2 k3  conductance (G0, up spin)  conductance (G0, down spin)
#
  0   0   0.0000   0.0000  1.022572e-04  2.338899e-05

下面是未修改之前的结果conductance文件:

#/************************************************************/
#
#  Conductance:
#
#    G0 = e^2/h
#       = (1.60217653 * 10^{-19})^2/(6.626076 * 10^{-34})
#       = 3.87404194169 * 10^{-4} [C^2 J^{-1} s^{-1}]
#    note that
#    e = 1.60217653 * 10^{-19} C
#    h = 6.626076 * 10^{-34} J s
#***********************************************************/
#
# SpinP_switch= 1
# Chemical potential (eV)      Left, Right Leads=-4.864470e+00 -5.164470e+00
# Chemical potential (Hartree) Left, Right Leads=-1.787660e-01 -1.897908e-01
# diff Chemical potential (eV)     = 2.999998e-01
# diff Chemical potential (Hartree)= 1.102479e-02
# Corresponding bias voltage (Volt)= 2.999998e-01
# average conductance for up and down spins (G0)= 3.095083e-03  4.359661e-02
#
# i2 i3 k2 k3  conductance (G0, up spin)  conductance (G0, down spin)
#
  0   0   0.0000   0.0000  3.095083e-03  4.359661e-02

上面是直接修改NEGF-8ZGNR-0.3这个例子。接下来需要修改的是我们自己的例子。也就是前三节刚做的那个例子。我们把左右电极的晶格矢量改写成下面这样:

<Atoms.UnitVectors 
 0.000  0.000  4.92
 10.000 0.000 0.000 
 0.000 27.000 0.000 
Atoms.UnitVectors>

因为两个电极的最大和最小的第三列坐标之差为2.4595,又因为2.4595*2=4.92,所以需要试2.4595和4.92这两个数。另外强调一下,0.3.dat文件里不需要Atoms.UnitVectors这个文件,只有左右电极的文件才需要。4.92仿真结果如下图所示:

2.4595这一组因为提示太近而跑不下去。3.5、3.6893也不行。4虽然能跑完电极,但跑不了NEGF主文件。4.5的结果如下

unit vector 4.5
unit vector 4.2
unit vector 4.1
unit vector 4.1 并且把scf.Kgrid 7 7 7 # means n1 x n2 x n3修改完成(原来是200 1 1).文件夹是work_negf_simplified_4.1_777
unit vector 4.1 并且把scf.Kgrid 1 1 200

附录:

为方便查阅,现将原文档8ZGNR的三个文件粘贴出来:

#
# File Name
#

System.CurrrentDirectory         ./    # default=./
System.Name                     lead-l-8zgnr
level.of.stdout                   1    # default=1 (1-3)
level.of.fileout                  0    # default=1 (0-2)


#
# Definition of Atomic Species
#

Species.Number       2
<Definition.of.Atomic.Species
 C  C5.0-s2p1  C_CA13
 H  H5.0-s2    H_CA13
Definition.of.Atomic.Species>


NEGF.output_hks    true
NEGF.filename.hks  lead-l-8zgnr.hks


#
# Atoms
#

Atoms.Number       36
Atoms.SpeciesAndCoordinates.Unit   Ang # Ang|AU
<Atoms.SpeciesAndCoordinates 
   1  C     0.0000005     0.0361626     0.0012037     2.5000000     1.5000000
   2  C     0.0000005     0.7140385     1.2211879     2.5000000     1.5000000
   3  C     0.0000004     2.1457201     1.2213117     2.5000000     1.5000000
   4  C     0.0000001     2.8501725     0.0013209     2.5000000     1.5000000
   5  C    -0.0000001     4.2700744     0.0012023     2.5000000     1.5000000
   6  C    -0.0000002     4.9766306     1.2212146     2.5000000     1.5000000
   7  C     0.0000001     6.3942050     1.2213979     2.5000000     1.5000000
   8  C     0.0000001     7.1014018     0.0014044     2.5000000     1.5000000
   9  C    -0.0000001     8.5186603     0.0014042     2.5000000     1.5000000
  10  C     0.0000001     9.2258519     1.2213977     2.5000000     1.5000000
  11  C     0.0000004    10.6434033     1.2212143     2.5000000     1.5000000
  12  C     0.0000003    11.3499296     0.0012020     2.5000000     1.5000000
  13  C    -0.0000002    12.7697977     0.0013211     2.5000000     1.5000000
  14  C    -0.0000004    13.4742418     1.2213119     2.5000000     1.5000000
  15  C    -0.0000004    14.9059644     1.2211875     2.5000000     1.5000000
  16  C    -0.0000004    15.5838921     0.0012032     2.5000000     1.5000000
  17  H     0.0000003    -1.0818549     0.0012134     0.5000000     0.5000000
  18  H    -0.0000005    16.7019089     0.0012129     0.5000000     0.5000000
  19  C     0.0000005     0.0361626     2.4412037     2.5000000     1.5000000
  20  C     0.0000005     0.7140385     3.6611879     2.5000000     1.5000000
  21  C     0.0000004     2.1457201     3.6613117     2.5000000     1.5000000
  22  C     0.0000001     2.8501725     2.4413209     2.5000000     1.5000000
  23  C    -0.0000001     4.2700744     2.4412023     2.5000000     1.5000000
  24  C    -0.0000002     4.9766306     3.6612146     2.5000000     1.5000000
  25  C     0.0000001     6.3942050     3.6613979     2.5000000     1.5000000
  26  C     0.0000001     7.1014018     2.4414044     2.5000000     1.5000000
  27  C    -0.0000001     8.5186603     2.4414042     2.5000000     1.5000000
  28  C     0.0000001     9.2258519     3.6613977     2.5000000     1.5000000
  29  C     0.0000004    10.6434033     3.6612143     2.5000000     1.5000000
  30  C     0.0000003    11.3499296     2.4412020     2.5000000     1.5000000
  31  C    -0.0000002    12.7697977     2.4413211     2.5000000     1.5000000
  32  C    -0.0000004    13.4742418     3.6613119     2.5000000     1.5000000
  33  C    -0.0000004    14.9059644     3.6611875     2.5000000     1.5000000
  34  C    -0.0000004    15.5838921     2.4412032     2.5000000     1.5000000
  35  H     0.0000003    -1.0818549     2.4412134     0.5000000     0.5000000
  36  H    -0.0000005    16.7019089     2.4412129     0.5000000     0.5000000
Atoms.SpeciesAndCoordinates>
Atoms.UnitVectors.Unit             Ang # Ang|AU
<Atoms.UnitVectors 
 0.000  0.000  4.88
 10.000 0.000 0.000 
 0.000 27.000 0.000 
Atoms.UnitVectors>


#
# SCF or Electronic System
#


scf.XcType                 LSDA-CA     # LDA|LSDA-CA|LSDA-PW|GGA-PBE
scf.SpinPolarization        on         # On|Off|NC
scf.ElectronicTemperature  300.0       # default=300 (K)
scf.energycutoff           120.0       # default=150 (Ry)
#scf.Ngrid                  140 140 100  # about=200 (Ry)
scf.maxIter                1000         # default=40
scf.EigenvalueSolver       Band     # DC|GDC|Cluster|Band
scf.lapack.dste            dstevx      # dstegr|dstedc|dstevx, default=dstegr
scf.Kgrid                  201 1 1       # means n1 x n2 x n3
scf.Mixing.Type           rmm-diisk    # Simple|Rmm-Diis|Gr-Pulay|Kerker|Rmm-Diisk
scf.Kerker.factor          10.0        # default=1
scf.Init.Mixing.Weight     0.001       # default=0.30
scf.Min.Mixing.Weight      0.001       # default=0.001
scf.Max.Mixing.Weight      0.105       # default=0.40
scf.Mixing.History          40       # default=5
scf.Mixing.StartPulay       10       # default=6
scf.Mixing.EveryPulay        1          # default=6
scf.criterion             1.0e-9       # default=1.0e-6 (Hartree)

#
# MD or Geometry Optimization
#

MD.Type               Nomd    # Opt|EF|BFGS|RF|DIIS
MD.Opt.DIIS.History   18    # default=3
MD.Opt.StartDIIS      6     # default=5
MD.Opt.EveryDIIS      5     # default=10
MD.maxIter            2000   #
MD.Opt.criterion    1.0e-4  # default=1.0e-4 (a.u.)


#
# Band dispersion
#

Band.dispersion              off      # on|off, default=off
<Band.KPath.UnitCell
 0.000 0.000  4.88
 10.000 0.000 0.000 
 0.000 27.000 0.000 
Band.KPath.UnitCell>

Band.Nkpath                1
<Band.kpath                
   601  0.0 0.0 0.0   0.5 0.0 0.0   G X
Band.kpath>


#
# DOS and PDOS
#

Dos.fileout                  off      # on|off, default=off
Dos.Erange              -20.0  20.0    # default = -20 20
Dos.Kgrid                601 1 1      # default = Kgrid1 Kgrid2 Kgrid3

DosGauss.fileout      off
DosGauss.Num.Mesh    4000
DosGauss.Width      0.01

MO.fileout                   off       # on|off
num.HOMOs                     1        # default=1
num.LUMOs                     1        # default=1
MO.Nkpoint                    1        # default=1
<MO.kpoint
  0.0  0.0  0.0
MO.kpoint>

#
# File Name
#

System.CurrrentDirectory         ./    # default=./
System.Name                      lead-r-8zgnr
level.of.stdout                   1    # default=1 (1-3)
level.of.fileout                  0    # default=1 (0-2)

#
# Definition of Atomic Species
#

Species.Number       2
<Definition.of.Atomic.Species
 C  C5.0-s2p1  C_CA13
 H  H5.0-s2    H_CA13
Definition.of.Atomic.Species>


NEGF.output_hks    true
NEGF.filename.hks  lead-r-8zgnr.hks


#
# Atoms
#

Atoms.Number      36
Atoms.SpeciesAndCoordinates.Unit   Ang # Ang|AU
<Atoms.SpeciesAndCoordinates 
   1  C     0.0000005     0.7140385     6.1011879     1.5000000     2.5000000
   2  C     0.0000004     2.1457201     6.1013117     1.5000000     2.5000000
   3  C    -0.0000002     4.9766306     6.1012146     1.5000000     2.5000000
   4  C     0.0000001     6.3942050     6.1013979     1.5000000     2.5000000
   5  C     0.0000001     9.2258519     6.1013977     1.5000000     2.5000000
   6  C     0.0000004    10.6434033     6.1012143     1.5000000     2.5000000
   7  C    -0.0000004    13.4742418     6.1013119     1.5000000     2.5000000
   8  C    -0.0000004    14.9059644     6.1011875     1.5000000     2.5000000
   9  C     0.0000005     0.0361626     7.3212037     1.5000000     2.5000000
  10  C     0.0000005     0.7140385     8.5411879     1.5000000     2.5000000
  11  C     0.0000004     2.1457201     8.5413117     1.5000000     2.5000000
  12  C     0.0000001     2.8501725     7.3213209     1.5000000     2.5000000
  13  C    -0.0000001     4.2700744     7.3212023     1.5000000     2.5000000
  14  C    -0.0000002     4.9766306     8.5412146     1.5000000     2.5000000
  15  C     0.0000001     6.3942050     8.5413979     1.5000000     2.5000000
  16  C     0.0000001     7.1014018     7.3214044     1.5000000     2.5000000
  17  C    -0.0000001     8.5186603     7.3214042     1.5000000     2.5000000
  18  C     0.0000001     9.2258519     8.5413977     1.5000000     2.5000000
  19  C     0.0000004    10.6434033     8.5412143     1.5000000     2.5000000
  20  C     0.0000003    11.3499296     7.3212020     1.5000000     2.5000000
  21  C    -0.0000002    12.7697977     7.3213211     1.5000000     2.5000000
  22  C    -0.0000004    13.4742418     8.5413119     1.5000000     2.5000000
  23  C    -0.0000004    14.9059644     8.5411875     1.5000000     2.5000000
  24  C    -0.0000004    15.5838921     7.3212032     1.5000000     2.5000000
  25  H     0.0000003    -1.0818549     7.3212134     0.5000000     0.5000000
  26  H    -0.0000005    16.7019089     7.3212129     0.5000000     0.5000000
  27  C     0.0000005     0.0361626     9.7612037     1.5000000     2.5000000
  28  C     0.0000001     2.8501725     9.7613209     1.5000000     2.5000000
  29  C    -0.0000001     4.2700744     9.7612023     1.5000000     2.5000000
  30  C     0.0000001     7.1014018     9.7614044     1.5000000     2.5000000
  31  C    -0.0000001     8.5186603     9.7614042     1.5000000     2.5000000
  32  C     0.0000003    11.3499296     9.7612020     1.5000000     2.5000000
  33  C    -0.0000002    12.7697977     9.7613211     1.5000000     2.5000000
  34  C    -0.0000004    15.5838921     9.7612032     1.5000000     2.5000000
  35  H     0.0000003    -1.0818549     9.7612134     0.5000000     0.5000000
  36  H    -0.0000005    16.7019089     9.7612129     0.5000000     0.5000000
Atoms.SpeciesAndCoordinates>
Atoms.UnitVectors.Unit             Ang # Ang|AU
<Atoms.UnitVectors 
 0.000 0.000  4.88
 10.000 0.000 0.000 
 0.000 27.000 0.000 
Atoms.UnitVectors>


#
# SCF or Electronic System
#


scf.XcType                 LSDA-CA     # LDA|LSDA-CA|LSDA-PW|GGA-PBE
scf.SpinPolarization        on         # On|Off|NC
scf.ElectronicTemperature  300.0       # default=300 (K)
scf.energycutoff           120.0       # default=150 (Ry)
#scf.Ngrid                  140 140 100  # about=200 (Ry)
scf.maxIter                1000         # default=40
scf.EigenvalueSolver       Band     # DC|GDC|Cluster|Band
scf.lapack.dste            dstevx      # dstegr|dstedc|dstevx, default=dstegr
scf.Kgrid                  201 1 1       # means n1 x n2 x n3
scf.Mixing.Type           rmm-diisk    # Simple|Rmm-Diis|Gr-Pulay|Kerker|Rmm-Diisk
scf.Kerker.factor          10.0        # default=1
scf.Init.Mixing.Weight     0.100       # default=0.30
scf.Min.Mixing.Weight      0.01        # default=0.001
scf.Max.Mixing.Weight      0.105       # default=0.40
scf.Mixing.History          40       # default=5
scf.Mixing.StartPulay       10       # default=6
scf.Mixing.EveryPulay        1          # default=6
scf.criterion             1.0e-9       # default=1.0e-6 (Hartree)

#
# MD or Geometry Optimization
#

MD.Type               Nomd    # Opt|EF|BFGS|RF|DIIS
MD.Opt.DIIS.History   18    # default=3
MD.Opt.StartDIIS      6     # default=5
MD.Opt.EveryDIIS      5     # default=10
MD.maxIter            2000   #
MD.Opt.criterion    1.0e-4  # default=1.0e-4 (a.u.)


#
# Band dispersion
#

Band.dispersion              off       # on|off, default=off
<Band.KPath.UnitCell
 0.000 0.000  4.88
 10.000 0.000 0.000 
 0.000 27.000 0.000 
Band.KPath.UnitCell>

Band.Nkpath                1
<Band.kpath                
   601  0.0 0.0 0.0   0.5 0.0 0.0   G X
Band.kpath>


#
# DOS and PDOS
#

Dos.fileout                  off      # on|off, default=off
Dos.Erange              -20.0  20.0    # default = -20 20
Dos.Kgrid                601 1 1      # default = Kgrid1 Kgrid2 Kgrid3

DosGauss.fileout      off
DosGauss.Num.Mesh    4000
DosGauss.Width      0.01

#
# File Name
#

System.CurrrentDirectory         ./    # default=./
System.Name                    negf-8zgnr-0.3
level.of.stdout                   1    # default=1 (1-3)
level.of.fileout                  1    # default=1 (0-2)


NEGF.filename.hks.l     lead-l-8zgnr.hks
NEGF.filename.hks.r     lead-r-8zgnr.hks

NEGF.SCF.Iter.Band          20
NEGF.Num.Poles             100       # defalut=150
NEGF.scf.Kgrid             1 1       # defalut=1 1

NEGF.bias.voltage         -0.3       # default=0.0 (eV)
NEGF.bias.neq.im.energy    0.01      # default=0.01 (eV)
NEGF.bias.neq.energy.step  0.02      # default=0.02 (eV)

Dos.fileout                 off       # on|off, default=off
NEGF.Dos.energyrange      -20.0 20.0 1.0e-4  #default=-10.0 10.0 5.0e-3 (eV)
NEGF.Dos.energy.div        4000       # default=200
NEGF.Dos.Kgrid             1 1       # default=1 1

#
# Transmission & Eigen Channel Analysis
#
 
NEGF.tran.Analysis          on       # default on
NEGF.tran.energydiv        800       # default=200
NEGF.tran.energyrange   -4 4 0.001   # default=-10.0 10.0 5.0e-3 (eV)
NEGF.tran.Kgrid            1 1       # default= 1 1

NEGF.tran.CurrentDensity    on       # default on

NEGF.tran.Channel           on       # default on
NEGF.Channel.Nkpoint         1       # default 1
<NEGF.Channel.kpoint
0.0  0.0
NEGF.Channel.kpoint>
# default 0.0  0.0
NEGF.Channel.Nenergy         1       # default 1 
<NEGF.Channel.energy
0.0
NEGF.Channel.energy>
# default 0.0
NEGF.Channel.Num            5        # defualt=5(for collinear), 10(for Non-collinear)


#
# Definition of Atomic Species
#

Species.Number       2
<Definition.of.Atomic.Species
 C  C5.0-s2p1  C_CA13
 H  H5.0-s2    H_CA13
Definition.of.Atomic.Species>


#
# Atoms
#

Atoms.SpeciesAndCoordinates.Unit   Ang # Ang|AU


Atoms.Number        10
<Atoms.SpeciesAndCoordinates
   1  C     0.0000005     0.0361626     4.8812037     2.1000000     1.9000000
   2  C     0.0000001     2.8501725     4.8813209     2.1000000     1.9000000
   3  C    -0.0000001     4.2700744     4.8812023     2.1000000     1.9000000
   4  C     0.0000001     7.1014018     4.8814044     2.1000000     1.9000000
   5  C    -0.0000001     8.5186603     4.8814042     2.1000000     1.9000000
   6  C     0.0000003    11.3499296     4.8812020     2.1000000     1.9000000
   7  C    -0.0000002    12.7697977     4.8813211     2.1000000     1.9000000
   8  C    -0.0000004    15.5838921     4.8812032     2.1000000     1.9000000
   9  H     0.0000003    -1.0818549     4.8812134     0.5000000     0.5000000
  10  H    -0.0000005    16.7019089     4.8812129     0.5000000     0.5000000
Atoms.SpeciesAndCoordinates>

Atoms.UnitVectors.Unit             Ang # Ang|AU

#<Atoms.UnitVectors
# 0.000 0.000  4.88
# 10.000 0.000 0.000
# 0.000 24.000 0.000
#Atoms.UnitVectors>


#
# Lead-Left
#

LeftLeadAtoms.Number         36
<LeftLeadAtoms.SpeciesAndCoordinates           # Unit=Ang.
   1  C     0.0000005     0.0361626     0.0012037     2.1000000     1.9000000
   2  C     0.0000005     0.7140385     1.2211879     2.1000000     1.9000000
   3  C     0.0000004     2.1457201     1.2213117     2.1000000     1.9000000
   4  C     0.0000001     2.8501725     0.0013209     2.1000000     1.9000000
   5  C    -0.0000001     4.2700744     0.0012023     2.1000000     1.9000000
   6  C    -0.0000002     4.9766306     1.2212146     2.1000000     1.9000000
   7  C     0.0000001     6.3942050     1.2213979     2.1000000     1.9000000
   8  C     0.0000001     7.1014018     0.0014044     2.1000000     1.9000000
   9  C    -0.0000001     8.5186603     0.0014042     2.1000000     1.9000000
  10  C     0.0000001     9.2258519     1.2213977     2.1000000     1.9000000
  11  C     0.0000004    10.6434033     1.2212143     2.1000000     1.9000000
  12  C     0.0000003    11.3499296     0.0012020     2.1000000     1.9000000
  13  C    -0.0000002    12.7697977     0.0013211     2.1000000     1.9000000
  14  C    -0.0000004    13.4742418     1.2213119     2.1000000     1.9000000
  15  C    -0.0000004    14.9059644     1.2211875     2.1000000     1.9000000
  16  C    -0.0000004    15.5838921     0.0012032     2.1000000     1.9000000
  17  H     0.0000003    -1.0818549     0.0012134     0.5000000     0.5000000
  18  H    -0.0000005    16.7019089     0.0012129     0.5000000     0.5000000
  19  C     0.0000005     0.0361626     2.4412037     2.1000000     1.9000000
  20  C     0.0000005     0.7140385     3.6611879     2.1000000     1.9000000
  21  C     0.0000004     2.1457201     3.6613117     2.1000000     1.9000000
  22  C     0.0000001     2.8501725     2.4413209     2.1000000     1.9000000
  23  C    -0.0000001     4.2700744     2.4412023     2.1000000     1.9000000
  24  C    -0.0000002     4.9766306     3.6612146     2.1000000     1.9000000
  25  C     0.0000001     6.3942050     3.6613979     2.1000000     1.9000000
  26  C     0.0000001     7.1014018     2.4414044     2.1000000     1.9000000
  27  C    -0.0000001     8.5186603     2.4414042     2.1000000     1.9000000
  28  C     0.0000001     9.2258519     3.6613977     2.1000000     1.9000000
  29  C     0.0000004    10.6434033     3.6612143     2.1000000     1.9000000
  30  C     0.0000003    11.3499296     2.4412020     2.1000000     1.9000000
  31  C    -0.0000002    12.7697977     2.4413211     2.1000000     1.9000000
  32  C    -0.0000004    13.4742418     3.6613119     2.1000000     1.9000000
  33  C    -0.0000004    14.9059644     3.6611875     2.1000000     1.9000000
  34  C    -0.0000004    15.5838921     2.4412032     2.1000000     1.9000000
  35  H     0.0000003    -1.0818549     2.4412134     0.5000000     0.5000000
  36  H    -0.0000005    16.7019089     2.4412129     0.5000000     0.5000000
LeftLeadAtoms.SpeciesAndCoordinates>


#
# Lead-Right
#


RightLeadAtoms.Number        36

<RightLeadAtoms.SpeciesAndCoordinates           # Unit=Ang.
   1  C     0.0000005     0.7140385     6.1011879     1.9000000     2.1000000
   2  C     0.0000004     2.1457201     6.1013117     1.9000000     2.1000000
   3  C    -0.0000002     4.9766306     6.1012146     1.9000000     2.1000000
   4  C     0.0000001     6.3942050     6.1013979     1.9000000     2.1000000
   5  C     0.0000001     9.2258519     6.1013977     1.9000000     2.1000000
   6  C     0.0000004    10.6434033     6.1012143     1.9000000     2.1000000
   7  C    -0.0000004    13.4742418     6.1013119     1.9000000     2.1000000
   8  C    -0.0000004    14.9059644     6.1011875     1.9000000     2.1000000
   9  C     0.0000005     0.0361626     7.3212037     1.9000000     2.1000000
  10  C     0.0000005     0.7140385     8.5411879     1.9000000     2.1000000
  11  C     0.0000004     2.1457201     8.5413117     1.9000000     2.1000000
  12  C     0.0000001     2.8501725     7.3213209     1.9000000     2.1000000
  13  C    -0.0000001     4.2700744     7.3212023     1.9000000     2.1000000
  14  C    -0.0000002     4.9766306     8.5412146     1.9000000     2.1000000
  15  C     0.0000001     6.3942050     8.5413979     1.9000000     2.1000000
  16  C     0.0000001     7.1014018     7.3214044     1.9000000     2.1000000
  17  C    -0.0000001     8.5186603     7.3214042     1.9000000     2.1000000
  18  C     0.0000001     9.2258519     8.5413977     1.9000000     2.1000000
  19  C     0.0000004    10.6434033     8.5412143     1.9000000     2.1000000
  20  C     0.0000003    11.3499296     7.3212020     1.9000000     2.1000000
  21  C    -0.0000002    12.7697977     7.3213211     1.9000000     2.1000000
  22  C    -0.0000004    13.4742418     8.5413119     1.9000000     2.1000000
  23  C    -0.0000004    14.9059644     8.5411875     1.9000000     2.1000000
  24  C    -0.0000004    15.5838921     7.3212032     1.9000000     2.1000000
  25  H     0.0000003    -1.0818549     7.3212134     0.5000000     0.5000000
  26  H    -0.0000005    16.7019089     7.3212129     0.5000000     0.5000000
  27  C     0.0000005     0.0361626     9.7612037     1.9000000     2.1000000
  28  C     0.0000001     2.8501725     9.7613209     1.9000000     2.1000000
  29  C    -0.0000001     4.2700744     9.7612023     1.9000000     2.1000000
  30  C     0.0000001     7.1014018     9.7614044     1.9000000     2.1000000
  31  C    -0.0000001     8.5186603     9.7614042     1.9000000     2.1000000
  32  C     0.0000003    11.3499296     9.7612020     1.9000000     2.1000000
  33  C    -0.0000002    12.7697977     9.7613211     1.9000000     2.1000000
  34  C    -0.0000004    15.5838921     9.7612032     1.9000000     2.1000000
  35  H     0.0000003    -1.0818549     9.7612134     0.5000000     0.5000000
  36  H    -0.0000005    16.7019089     9.7612129     0.5000000     0.5000000
RightLeadAtoms.SpeciesAndCoordinates>

#
# SCF or Electronic System
#


scf.XcType                LSDA-CA      # LDA|LSDA-CA|LSDA-PW|GGA-PBE
scf.SpinPolarization        on         # On|Off|NC
scf.ElectronicTemperature  300.0       # default=300 (K)
scf.energycutoff           120.0       # default=150 (Ry)
#scf.Ngrid                  140 140 100 # about=200 (Ry)
scf.maxIter                2000         # default=40
scf.EigenvalueSolver       NEGF        # DC|GDC|Cluster|Band
scf.lapack.dste            dstevx      # dstegr|dstedc|dstevx, default=dstegr
scf.Kgrid                  1 1 1       # means n1 x n2 x n3
scf.Mixing.Type           rmm-diisk    # Simple|Rmm-Diis|Gr-Pulay|Kerker|Rmm-Diisk
scf.Kerker.factor          20.0        # default=1
scf.Init.Mixing.Weight     0.1         # default=0.30
scf.Min.Mixing.Weight      0.04        # default=0.001
scf.Max.Mixing.Weight      0.4         # default=0.40
scf.Mixing.History          60         # default=5
scf.Mixing.StartPulay       20         # default=6
scf.Mixing.EveryPulay       1          # default=6
scf.criterion             1.0e-7       # default=1.0e-6 (Hartree)


#
# MD or Geometry Optimization
#

MD.Type               Nomd    # Opt|EF|BFGS|RF|DIIS
MD.Opt.DIIS.History   10    # default=3
MD.Opt.StartDIIS      6     # default=5
MD.Opt.EveryDIIS      5     # default=10
MD.maxIter            200   #
MD.Opt.criterion    1.0e-4  # default=1.0e-4 (a.u.)

发表评论