diff --git a/Test/Au/i_lsms b/Test/Au/i_lsms new file mode 100644 index 000000000..92c6d815c --- /dev/null +++ b/Test/Au/i_lsms @@ -0,0 +1,75 @@ +systemid="au" +system_title = "Gold (FCC) Test for LSMS 3" +pot_in_type=1 +-- write potential out. Use this only when xRepeat, yRepeat and zRepeat = 1 +pot_out_type=1 +num_atoms=1 +nspin=3 +mtasa=0 +xcFunctional={0,1} -- von Barth Hedin (LSMS_1) +-- xcFunctional={1,1,17} -- von Barth Hedin (libxc) +-- xcFunctional={1,1,12} -- Perdew & Wang +-- xcFunctional={1,20} -- Teter 93 +-- xcFunctional={1,101,130} -- GGA PBE + +iprint=0 +default_iprint=-1 +print_node=0 +-- istop="buildKKRMatrix" +istop="main" + +nscf=50 +rmsTolerance=1.e-12 + +xRepeat=1 +yRepeat=1 +zRepeat=1 +makeTypesUnique=1 + +-- energyContour = {npts=1,grid=0,ebot=-0.2,etop=0.0,eitop=0.825,eibot=0.1} +-- energyContour = {npts=11,grid=2,ebot=-0.3,etop=0.0,eitop=0.825,eibot=0.025,maxGroupSize=5} +energyContour = {npts=25,grid=2,ebot=-0.5,etop=0.0,eitop=0.825,eibot=0.0025} + +a = 7.7067 + +bravais = {} +bravais[1]={0.5*a,0.5*a,0} +bravais[2]={0,0.5*a,0.5*a} +bravais[3]={0.5*a,0,0.5*a} + +site_default={lmax=3,rLIZ=12.5,rsteps={89.5,91.5,93.2,99.9},atom="Au",Z=79,Zc=46,Zs=22,Zv=11,rad=2} + +-- Possible choices of mixing parameters +-- quantity : charge, potential, moment_magnitude, moment_direction, noQuantity +-- algorithm : simple, broyden, noAlgorithm +-- mixing = { {quantity = "charge", algorithm = "simple", mixing_parameter = 0.05} } +-- mixing = { {quantity = "potential", algorithm = "simple", mixing_parameter = 0.05} } +-- mixing = { {quantity = "potential", algorithm = "broyden", mixing_parameter = 0.05}, {quantity = "moment_magnitude", algorithm = "noAlgorithm", mixing_parameter = 0.01} } +mixing = { {quantity = "potential", algorithm = "broyden", mixing_parameter = 0.05} } + +numberOfMixQuantities = 0 +for k,v in pairs(mixing) do + numberOfMixQuantities = numberOfMixQuantities + 1 +end + +--[[ +print("Site defaults:") +for k,v in pairs(site_default) do + print(k," : ",v) +end +--]] + +site = {} +for i =1,num_atoms do site[i]={} end + +site[1].pos={0,0,0} +site[1].evec={0,0,1} + +-- set site defaults +for i =1,num_atoms do + for k,v in pairs(site_default) do + if(site[i][k]==nil) then site[i][k]=v end + end +end + + diff --git a/Test/Au/v_au.0 b/Test/Au/v_au.0 new file mode 100644 index 000000000..0b5b5c2dc --- /dev/null +++ b/Test/Au/v_au.0 @@ -0,0 +1,407 @@ + nspin vdif :::::: Potentials for gold 1.00 + 1 0.0000000000000D+00 + NEWSS: gold z = 79. no spin + 79. 7.71010 46. 0.6772475579292E+00 + -0.1113096714020E+02 0.1002810386820E+01 501 501 +-0.1579862475472E+03-0.1579859096887E+03-0.1579855635291E+03-0.1579852088643E+03 +-0.1579848454853E+03-0.1579844731780E+03-0.1579840917229E+03-0.1579837008951E+03 +-0.1579833004645E+03-0.1579828901948E+03-0.1579824698444E+03-0.1579820391654E+03 +-0.1579815979040E+03-0.1579811458001E+03-0.1579806825872E+03-0.1579802079922E+03 +-0.1579797217354E+03-0.1579792235301E+03-0.1579787130826E+03-0.1579781900921E+03 +-0.1579776542500E+03-0.1579771052407E+03-0.1579765427403E+03-0.1579759664172E+03 +-0.1579753759317E+03-0.1579747709355E+03-0.1579741510720E+03-0.1579735159756E+03 +-0.1579728652718E+03-0.1579721985770E+03-0.1579715154979E+03-0.1579708156317E+03 +-0.1579700985657E+03-0.1579693638770E+03-0.1579686111323E+03-0.1579678398875E+03 +-0.1579670496879E+03-0.1579662400672E+03-0.1579654105480E+03-0.1579645606408E+03 +-0.1579636898444E+03-0.1579627976449E+03-0.1579618835161E+03-0.1579609469186E+03 +-0.1579599872997E+03-0.1579590040933E+03-0.1579579967192E+03-0.1579569645830E+03 +-0.1579559070755E+03-0.1579548235727E+03-0.1579537134350E+03-0.1579525760073E+03 +-0.1579514106182E+03-0.1579502165798E+03-0.1579489931872E+03-0.1579477397183E+03 +-0.1579464554328E+03-0.1579451395728E+03-0.1579437913610E+03-0.1579424100015E+03 +-0.1579409946785E+03-0.1579395445562E+03-0.1579380587781E+03-0.1579365364667E+03 +-0.1579349767226E+03-0.1579333786246E+03-0.1579317412285E+03-0.1579300635669E+03 +-0.1579283446486E+03-0.1579265834577E+03-0.1579247789535E+03-0.1579229300695E+03 +-0.1579210357128E+03-0.1579190947638E+03-0.1579171060750E+03-0.1579150684705E+03 +-0.1579129807456E+03-0.1579108416658E+03-0.1579086499660E+03-0.1579064043500E+03 +-0.1579041034893E+03-0.1579017460228E+03-0.1578993305558E+03-0.1578968556590E+03 +-0.1578943198678E+03-0.1578917216813E+03-0.1578890595617E+03-0.1578863319330E+03 +-0.1578835371802E+03-0.1578806736485E+03-0.1578777396420E+03-0.1578747334228E+03 +-0.1578716532102E+03-0.1578684971791E+03-0.1578652634594E+03-0.1578619501347E+03 +-0.1578585552411E+03-0.1578550768222E+03-0.1578515128205E+03-0.1578478610750E+03 +-0.1578441194207E+03-0.1578402856386E+03-0.1578363574551E+03-0.1578323325399E+03 +-0.1578282085050E+03-0.1578239829031E+03-0.1578196532264E+03-0.1578152169044E+03 +-0.1578106713033E+03-0.1578060137234E+03-0.1578012413983E+03-0.1577963514927E+03 +-0.1577913411007E+03-0.1577862072445E+03-0.1577809470055E+03-0.1577755571770E+03 +-0.1577700345276E+03-0.1577643757734E+03-0.1577585775491E+03-0.1577526364055E+03 +-0.1577465488076E+03-0.1577403111321E+03-0.1577339196658E+03-0.1577273706027E+03 +-0.1577206600420E+03-0.1577137839855E+03-0.1577067385477E+03-0.1576995193700E+03 +-0.1576921221432E+03-0.1576845424581E+03-0.1576767757952E+03-0.1576688175226E+03 +-0.1576606628923E+03-0.1576523070378E+03-0.1576437449708E+03-0.1576349717847E+03 +-0.1576259821557E+03-0.1576167706613E+03-0.1576073317998E+03-0.1575976599318E+03 +-0.1575877492766E+03-0.1575775939086E+03-0.1575671878293E+03-0.1575565250172E+03 +-0.1575455989077E+03-0.1575344029647E+03-0.1575229304883E+03-0.1575111746105E+03 +-0.1574991282908E+03-0.1574867846551E+03-0.1574741361192E+03-0.1574611750757E+03 +-0.1574478937583E+03-0.1574342842056E+03-0.1574203384132E+03-0.1574060482034E+03 +-0.1573914048171E+03-0.1573763994720E+03-0.1573610231648E+03-0.1573452669121E+03 +-0.1573291213600E+03-0.1573125766577E+03-0.1572956228738E+03-0.1572782499202E+03 +-0.1572604477677E+03-0.1572422054534E+03-0.1572235120224E+03-0.1572043563837E+03 +-0.1571847274214E+03-0.1571646130786E+03-0.1571440012727E+03-0.1571228800350E+03 +-0.1571012367049E+03-0.1570790581451E+03-0.1570563312646E+03-0.1570330427816E+03 +-0.1570091784576E+03-0.1569847241709E+03-0.1569596657006E+03-0.1569339877274E+03 +-0.1569076751802E+03-0.1568807126050E+03-0.1568530835951E+03-0.1568247722562E+03 +-0.1567957615233E+03-0.1567660341642E+03-0.1567355729643E+03-0.1567043593938E+03 +-0.1566723756960E+03-0.1566396024249E+03-0.1566060206698E+03-0.1565716104739E+03 +-0.1565363517412E+03-0.1565002237584E+03-0.1564632053931E+03-0.1564252749724E+03 +-0.1563864104510E+03-0.1563465889974E+03-0.1563057877912E+03-0.1562639827107E+03 +-0.1562211500424E+03-0.1561772645999E+03-0.1561323012804E+03-0.1560862343066E+03 +-0.1560390369737E+03-0.1559906825450E+03-0.1559411433434E+03-0.1558903910321E+03 +-0.1558383969474E+03-0.1557851318150E+03-0.1557305653557E+03-0.1556746669684E+03 +-0.1556174053980E+03-0.1555587488825E+03-0.1554986646174E+03-0.1554371193821E+03 +-0.1553740792868E+03-0.1553095097635E+03-0.1552433755617E+03-0.1551756407445E+03 +-0.1551062687074E+03-0.1550352226138E+03-0.1549624629060E+03-0.1548879531185E+03 +-0.1548116511978E+03-0.1547335196465E+03-0.1546535153287E+03-0.1545715984376E+03 +-0.1544877255826E+03-0.1544018539360E+03-0.1543139407156E+03-0.1542239399562E+03 +-0.1541318078613E+03-0.1540374979333E+03-0.1539409631856E+03-0.1538421572625E+03 +-0.1537410316491E+03-0.1536375372568E+03-0.1535316248875E+03-0.1534232451236E+03 +-0.1533123461999E+03-0.1531988766810E+03-0.1530827844241E+03-0.1529640165470E+03 +-0.1528425197242E+03-0.1527182392516E+03-0.1525911201486E+03-0.1524611067784E+03 +-0.1523281428151E+03-0.1521921712513E+03-0.1520531344044E+03-0.1519109739201E+03 +-0.1517656307746E+03-0.1516170452755E+03-0.1514651570607E+03-0.1513099050966E+03 +-0.1511512276786E+03-0.1509890625573E+03-0.1508233465506E+03-0.1506540158952E+03 +-0.1504810061311E+03-0.1503042521451E+03-0.1501236880283E+03-0.1499392471303E+03 +-0.1497508620765E+03-0.1495584647089E+03-0.1493619859398E+03-0.1491613562144E+03 +-0.1489565051058E+03-0.1487473607335E+03-0.1485338507639E+03-0.1483159023582E+03 +-0.1480934413371E+03-0.1478663925641E+03-0.1476346799442E+03-0.1473982264162E+03 +-0.1471569539212E+03-0.1469107832291E+03-0.1466596346668E+03-0.1464034272057E+03 +-0.1461420788412E+03-0.1458755064403E+03-0.1456036259039E+03-0.1453263525748E+03 +-0.1450436007918E+03-0.1447552841345E+03-0.1444613155159E+03-0.1441616072891E+03 +-0.1438560702304E+03-0.1435446158502E+03-0.1432271558036E+03-0.1429036007486E+03 +-0.1425738606480E+03-0.1422378480034E+03-0.1418954729924E+03-0.1415466483866E+03 +-0.1411912863201E+03-0.1408293006143E+03-0.1404606048369E+03-0.1400851152400E+03 +-0.1397027472136E+03-0.1393134170642E+03-0.1389170449157E+03-0.1385135506571E+03 +-0.1381028531396E+03-0.1376848724928E+03-0.1372595290695E+03-0.1368267430897E+03 +-0.1363864342697E+03-0.1359385214514E+03-0.1354829222552E+03-0.1350195527809E+03 +-0.1345483273833E+03-0.1340691585486E+03-0.1335819560731E+03-0.1330866243468E+03 +-0.1325830730302E+03-0.1320712077676E+03-0.1315509329988E+03-0.1310221460300E+03 +-0.1304847566818E+03-0.1299386700632E+03-0.1293837886208E+03-0.1288200293679E+03 +-0.1282472988620E+03-0.1276655234350E+03-0.1270746208949E+03-0.1264745270942E+03 +-0.1258651762268E+03-0.1252465127315E+03-0.1246184896503E+03-0.1239810651086E+03 +-0.1233341995641E+03-0.1226778750660E+03-0.1220120709154E+03-0.1213367704633E+03 +-0.1206519655684E+03-0.1199576531505E+03-0.1192538345002E+03-0.1185405146067E+03 +-0.1178177015393E+03-0.1170854059196E+03-0.1163436405212E+03-0.1155924168024E+03 +-0.1148317422781E+03-0.1140616400328E+03-0.1132821291668E+03-0.1124932257562E+03 +-0.1116949418408E+03-0.1108873143605E+03-0.1100703654863E+03-0.1092441273127E+03 +-0.1084086492487E+03-0.1075639698257E+03-0.1067101593744E+03-0.1058472800054E+03 +-0.1049754161934E+03-0.1040946640841E+03-0.1032051201533E+03-0.1023069039252E+03 +-0.1014001398286E+03-0.1004849569847E+03-0.9956147274760E+02-0.9862985136389E+02 +-0.9769023185633E+02-0.9674274429459E+02-0.9578752697343E+02-0.9482471834943E+02 +-0.9385445865955E+02-0.9287689251598E+02-0.9189217222742E+02-0.9090043828694E+02 +-0.8990187458654E+02-0.8889668181383E+02-0.8788507587937E+02-0.8686728013487E+02 +-0.8584353171632E+02-0.8481412886471E+02-0.8377935054793E+02-0.8273948591542E+02 +-0.8169487125050E+02-0.8064578099587E+02-0.7959256074316E+02-0.7853547084528E+02 +-0.7747482170476E+02-0.7641085968945E+02-0.7534381140764E+02-0.7427390370042E+02 +-0.7320131057375E+02-0.7212617738473E+02-0.7104864230807E+02-0.6996880223150E+02 +-0.6888680438126E+02-0.6780276765594E+02-0.6671679728140E+02-0.6562904521125E+02 +-0.6453969820228E+02-0.6344898570382E+02-0.6235718437426E+02-0.6126461824676E+02 +-0.6017164571874E+02-0.5907861559994E+02-0.5798597704325E+02-0.5689414737164E+02 +-0.5580352000055E+02-0.5471441077687E+02-0.5362722779155E+02-0.5254228429562E+02 +-0.5145987443935E+02-0.5038039618124E+02-0.4930417017926E+02-0.4823170791496E+02 +-0.4716349456039E+02-0.4610018889047E+02-0.4504251516734E+02-0.4399126998622E+02 +-0.4294738759992E+02-0.4191183518635E+02-0.4088562773889E+02-0.3986980034802E+02 +-0.3886534176324E+02-0.3787331763323E+02-0.3689467202425E+02-0.3593024273905E+02 +-0.3498079353082E+02-0.3404697458995E+02-0.3312931222339E+02-0.3222820339281E+02 +-0.3134388127849E+02-0.3047638743661E+02-0.2962580308793E+02-0.2879203362064E+02 +-0.2797488789835E+02-0.2717390063641E+02-0.2638887909651E+02-0.2561945278770E+02 +-0.2486507555264E+02-0.2412551875864E+02-0.2340026800231E+02-0.2268910289660E+02 +-0.2199161178392E+02-0.2130764338279E+02-0.2063688998638E+02-0.1997920613254E+02 +-0.1933440806782E+02-0.1870228216228E+02-0.1808261625023E+02-0.1747519683093E+02 +-0.1687974950997E+02-0.1629597686140E+02-0.1572356807054E+02-0.1516221408266E+02 +-0.1461162507623E+02-0.1407154743908E+02-0.1354177765616E+02-0.1302216063535E+02 +-0.1251261234356E+02-0.1201311914445E+02-0.1152370675115E+02-0.1104445389318E+02 +-0.1057548197636E+02-0.1011692845315E+02-0.9668947072183E+01-0.9231699274782E+01 +-0.8805330967911E+01-0.8389971981578E+01-0.7985733297765E+01-0.7592696505290E+01 +-0.7210910498387E+01-0.6840388634392E+01-0.6481113256360E+01-0.6133035303348E+01 +-0.5796074311278E+01-0.5470116194374E+01-0.5155034086191E+01-0.4850670194820E+01 +-0.4556858646660E+01-0.4273419900468E+01-0.4000166911439E+01-0.3736910617322E+01 +-0.3483463548708E+01-0.3239643420548E+01-0.3005276411389E+01-0.2780196774145E+01 +-0.2564257374078E+01-0.2357320810989E+01-0.2159271307744E+01-0.1970008717618E+01 +-0.1789452665670E+01-0.1617543058277E+01-0.1454240734886E+01-0.1299527880203E+01 +-0.1153408231930E+01-0.1015909586914E+01-0.8870825887660E+00-0.7670027940329E+00 +-0.6557718699380E+00-0.5535188900939E+00-0.4604019386691E+00-0.3766095458764E+00 +-0.3023626993334E+00 + 0.0000000000000E+00 muffin-tin zero + 505 0.3300000000000E+02 0.7819433444120E+02 + 0.1036149577238E-02 0.1084733986847E-02 0.1135662250028E-02 0.1189049154424E-02 + 0.1245015152664E-02 0.1303686643020E-02 0.1365196263997E-02 0.1429683203552E-02 + 0.1497293523657E-02 0.1568180500978E-02 0.1642504984464E-02 0.1720435770700E-02 + 0.1802149997879E-02 0.1887833559356E-02 0.1977681537724E-02 0.2071898660446E-02 + 0.2170699778118E-02 0.2274310366479E-02 0.2382967053352E-02 0.2496918171758E-02 + 0.2616424340494E-02 0.2741759073554E-02 0.2873209419815E-02 0.3011076634497E-02 + 0.3155676883974E-02 0.3307341985598E-02 0.3466420184268E-02 0.3633276967587E-02 + 0.3808295921496E-02 0.3991879628434E-02 0.4184450610101E-02 0.4386452317071E-02 + 0.4598350167564E-02 0.4820632637832E-02 0.5053812406717E-02 0.5298427557078E-02 + 0.5555042836914E-02 0.5824250983151E-02 0.6106674111203E-02 0.6402965173582E-02 + 0.6713809490989E-02 0.7039926359485E-02 0.7382070737532E-02 0.7741035016862E-02 + 0.8117650881359E-02 0.8512791258310E-02 0.8927372366630E-02 0.9362355866883E-02 + 0.9818751118150E-02 0.1029761754706E-01 0.1080006713457E-01 0.1132726702632E-01 + 0.1188044227276E-01 0.1246087870543E-01 0.1306992595624E-01 0.1370900062675E-01 + 0.1437958961509E-01 0.1508325360812E-01 0.1582163074722E-01 0.1659644047638E-01 + 0.1740948758137E-01 0.1826266642992E-01 0.1915796542244E-01 0.2009747166420E-01 + 0.2108337586957E-01 0.2211797751017E-01 0.2320369021889E-01 0.2434304746255E-01 + 0.2553870849657E-01 0.2679346461566E-01 0.2811024571519E-01 0.2949212717884E-01 + 0.3094233710858E-01 0.3246426391410E-01 0.3406146427957E-01 0.3573767152640E-01 + 0.3749680439186E-01 0.3934297624401E-01 0.4128050475483E-01 0.4331392205428E-01 + 0.4544798538912E-01 0.4768768831168E-01 0.5003827242496E-01 0.5250523971157E-01 + 0.5509436547573E-01 0.5781171192866E-01 0.6066364244948E-01 0.6365683655508E-01 + 0.6679830561434E-01 0.7009540934363E-01 0.7355587312246E-01 0.7718780617009E-01 + 0.8099972062586E-01 0.8500055157825E-01 0.8919967808976E-01 0.9360694526722E-01 + 0.9823268742946E-01 0.1030642042826E+00 0.1081108934555E+00 0.1134044378876E+00 + 0.1189570128429E+00 0.1247813965978E+00 0.1308910003687E+00 0.1372998997272E+00 + 0.1440228675753E+00 0.1510754087598E+00 0.1584737964060E+00 0.1662351100568E+00 + 0.1743772757067E+00 0.1829191078254E+00 0.1918803534696E+00 0.2012817385858E+00 + 0.2111450166154E+00 0.2214930195140E+00 0.2322940591880E+00 0.2436060266078E+00 + 0.2554660559464E+00 0.2679010471776E+00 0.2809392293170E+00 0.2946102263030E+00 + 0.3089451261492E+00 0.3239765535303E+00 0.3397387459724E+00 0.3562676338270E+00 + 0.3736009242159E+00 0.3917781891452E+00 0.4107527210438E+00 0.4306338951290E+00 + 0.4514692712051E+00 0.4733054555777E+00 0.4961913515504E+00 0.5201782731474E+00 + 0.5453200644797E+00 0.5716732250347E+00 0.5992970411834E+00 0.6281681819675E+00 + 0.6583862100209E+00 0.6900404841869E+00 0.7232007970851E+00 0.7579403721125E+00 + 0.7943360330764E+00 0.8324683822408E+00 0.8723907551372E+00 0.9141070622479E+00 + 0.9577857706361E+00 0.1003521717691E+01 0.1051414389547E+01 0.1101568150613E+01 + 0.1154092484476E+01 0.1208960503346E+01 0.1266369260559E+01 0.1326454480594E+01 + 0.1389344882644E+01 0.1455175479622E+01 0.1524022991455E+01 0.1595958016308E+01 + 0.1671210972548E+01 0.1749940492414E+01 0.1832312940659E+01 0.1918401477569E+01 + 0.2008344428558E+01 0.2102395125922E+01 0.2200748883245E+01 0.2303572190505E+01 + 0.2410912083456E+01 0.2523100867387E+01 0.2640367667583E+01 0.2762895657100E+01 + 0.2890773804945E+01 0.3024369285471E+01 0.3163950655851E+01 0.3309628778024E+01 + 0.3461682694641E+01 0.3620469482801E+01 0.3786208218379E+01 0.3959074082248E+01 + 0.4139506996571E+01 0.4327772301807E+01 0.4524045800464E+01 0.4728816180388E+01 + 0.4942336036642E+01 0.5164874313663E+01 0.5396944013842E+01 0.5638706561922E+01 + 0.5890648265940E+01 0.6153172122277E+01 0.6426520039870E+01 0.6711306016083E+01 + 0.7007672746601E+01 0.7316277370753E+01 0.7637394991928E+01 0.7971548814809E+01 + 0.8319143743835E+01 0.8680653105383E+01 0.9056536411395E+01 0.9447282900944E+01 + 0.9853341676336E+01 0.1027528101435E+02 0.1071346890160E+02 0.1116858811117E+02 + 0.1164095592282E+02 0.1213126314964E+02 0.1263995500823E+02 0.1316752660990E+02 + 0.1371461656983E+02 0.1428165179210E+02 0.1486919899941E+02 0.1547784933333E+02 + 0.1610808449771E+02 0.1676038301904E+02 0.1743536779996E+02 0.1813353630783E+02 + 0.1885538488292E+02 0.1960134734169E+02 0.2037199574270E+02 0.2116779419661E+02 + 0.2198919170724E+02 0.2283662050015E+02 0.2371049251742E+02 0.2461119591072E+02 + 0.2553908297840E+02 0.2649430222631E+02 0.2747772452266E+02 0.2848876442585E+02 + 0.2952861505156E+02 0.3059646651539E+02 0.3169328597257E+02 0.3281849239928E+02 + 0.3397250538575E+02 0.3515520755264E+02 0.3636614540069E+02 0.3760568124807E+02 + 0.3887306286625E+02 0.4016817282127E+02 0.4149072059245E+02 0.4283985434987E+02 + 0.4421513946269E+02 0.4561598816516E+02 0.4704150219302E+02 0.4849051574706E+02 + 0.4996230915635E+02 0.5145569414815E+02 0.5296936836356E+02 0.5450192288002E+02 + 0.5605174599074E+02 0.5761731841929E+02 0.5919685689315E+02 0.6078843787111E+02 + 0.6239000739708E+02 0.6399937837513E+02 0.6561422915089E+02 0.6723210332454E+02 + 0.6885041073113E+02 0.7046642952265E+02 0.7207730929119E+02 0.7368007517785E+02 + 0.7527163246988E+02 0.7684875526122E+02 0.7840815689338E+02 0.7994643341344E+02 + 0.8146009362588E+02 0.8294558513725E+02 0.8439929972147E+02 0.8581756397791E+02 + 0.8719672985500E+02 0.8853310370530E+02 0.8982314870031E+02 0.9106302218907E+02 + 0.9224908764252E+02 0.9337819216742E+02 0.9444693706410E+02 0.9545187330191E+02 + 0.9638999921370E+02 0.9725851424436E+02 0.9805484714747E+02 0.9877668840437E+02 + 0.9942203488754E+02 0.9998930441821E+02 0.1004770135143E+03 0.1008842647378E+03 + 0.1012105938521E+03 0.1014560818292E+03 0.1016213215629E+03 0.1017072822360E+03 + 0.1017155657779E+03 0.1016483481478E+03 0.1015083950918E+03 0.1012990738028E+03 + 0.1010248468897E+03 0.1006903341656E+03 0.1003006402777E+03 0.9986189086190E+02 + 0.9938114141016E+02 0.9886511870799E+02 0.9832221692380E+02 0.9776058357070E+02 + 0.9718938237791E+02 0.9661794046823E+02 0.9605638485301E+02 0.9551465817740E+02 + 0.9500355619708E+02 0.9453398968439E+02 0.9411627420718E+02 0.9376112229936E+02 + 0.9347948366224E+02 0.9328185106972E+02 0.9317833637367E+02 0.9317855529226E+02 + 0.9329150279453E+02 0.9352542282273E+02 0.9388767608684E+02 0.9438460869951E+02 + 0.9502142551768E+02 0.9580207042925E+02 0.9672902481518E+02 0.9780282308117E+02 + 0.9902316904855E+02 0.1003878274506E+03 0.1018927635831E+03 0.1035310345379E+03 + 0.1052956761938E+03 0.1071769637451E+03 0.1091624435156E+03 0.1112399005842E+03 + 0.1133927369796E+03 0.1156055493676E+03 0.1178587316007E+03 0.1201336493325E+03 + 0.1224088879167E+03 0.1246627521035E+03 0.1268730924368E+03 0.1290169488632E+03 + 0.1310707304984E+03 0.1330130650404E+03 0.1348216939949E+03 0.1364752084738E+03 + 0.1379539150622E+03 0.1392398290156E+03 0.1403170308449E+03 0.1411720116590E+03 + 0.1417939965461E+03 0.1421752361326E+03 0.1423112576723E+03 0.1422014352293E+03 + 0.1418495942452E+03 0.1412617561162E+03 0.1404483572730E+03 0.1394244056054E+03 + 0.1382097083501E+03 0.1368246363471E+03 0.1352957078005E+03 0.1336518857860E+03 + 0.1319232922682E+03 0.1301448901872E+03 0.1283501310782E+03 0.1265758882627E+03 + 0.1248577365968E+03 0.1232308818783E+03 0.1217300582720E+03 0.1203867085450E+03 + 0.1192298200892E+03 0.1182847260640E+03 0.1175719903776E+03 0.1171071250027E+03 + 0.1169003603011E+03 0.1169553030669E+03 0.1172690372577E+03 0.1178318989839E+03 + 0.1186274359809E+03 0.1196325338568E+03 0.1208177149482E+03 0.1221453112789E+03 + 0.1235750776548E+03 0.1250623266897E+03 0.1265583637581E+03 0.1280107021435E+03 + 0.1293651704920E+03 0.1305713628989E+03 0.1315772327389E+03 0.1323338239390E+03 + 0.1327983787743E+03 0.1329314493174E+03 0.1327021818590E+03 0.1320869923830E+03 + 0.1310709144613E+03 0.1296489267487E+03 0.1278264470052E+03 0.1256176173910E+03 + 0.1230478340939E+03 0.1201525725245E+03 0.1169754064617E+03 0.1135696490069E+03 + 0.1099913584369E+03 0.1063040070020E+03 0.1025765856025E+03 0.9887833494122E+02 + 0.9527834060698E+02 0.9184359863267E+02 0.8863711495400E+02 0.8571610880731E+02 + 0.8313048863014E+02 0.8092097399134E+02 0.7911687440875E+02 0.7773745064903E+02 + 0.7678995598186E+02 0.7626697365177E+02 0.7615361145140E+02 0.7642049747075E+02 + 0.7702636062651E+02 0.7792561582767E+02 0.7905712381707E+02 0.8036205130125E+02 + 0.8176953989907E+02 0.8321222583692E+02 0.8461999554738E+02 0.8592308940495E+02 + 0.8705931786863E+02 0.8796954945106E+02 0.8860215777508E+02 0.8891418487381E+02 + 0.8887254486716E+02 0.8845405728214E+02 0.8764526731350E+02 0.8644386935655E+02 + 0.8485760254076E+02 0.8290366316305E+02 0.8060781933371E+02 0.7800331520490E+02 + 0.7513055187839E+02 0.7203637747868E+02 0.6876657674764E+02 0.6537089256716E+02 + 0.6189979117377E+02 0.5840666168594E+02 0.5493617131793E+02 0.5153295107528E+02 + 0.4823923049145E+02 0.4508862914886E+02 0.4211200780332E+02 0.3933279424064E+02 + 0.3676874043077E+02 0.3443171735321E+02 0.3232689971045E+02 0.3045459721291E+02 + 0.2880978422961E+02 0.2738235609099E+02 0.2615850796757E+02 0.2512175799834E+02 + 0.2425265906389E+02 0.2353010528251E+02 0.2293205038355E+02 0.2243618623090E+02 + 0.2202056564528E+02 0.2166415968656E+02 0.2134733101867E+02 0.2105208922808E+02 + 0.2076265056525E+02 0.2046557951226E+02 0.2014966019622E+02 0.1980615922135E+02 + 0.1942873602420E+02 0.1901326555044E+02 0.1855774050141E+02 0.1806200934002E+02 + 0.1752757493345E+02 0.1695733267220E+02 0.1635526997424E+02 0.1572622517641E+02 + 0.1507563593868E+02 0.1440930652584E+02 0.1373316665566E+02 0.1305309253181E+02 + 0.1237475633369E+02 0.1170350564714E+02 0.1104421159184E+02 0.1040125476869E+02 + 0.9778433576256E+01 0.9178962344824E+01 0.8605460591691E+01 0.8059966450144E+01 + 0.7543965583779E+01 0.7058430737576E+01 0.6603869270690E+01 0.6180366547480E+01 + 0.5787662820968E+01 0.5425179137063E+01 0.5092103064489E+01 0.4787418239720E+01 + 0.4509963351201E+01 0.4258476131838E+01 0.4031633263013E+01 0.3828084406693E+01 + 0.3646480291152E+01 0.3485504764848E+01 0.3343886652725E+01 0.3220419744793E+01 + 0.3113974408761E+01 0.3023505173786E+01 0.2948055492879E+01 0.2886758454556E+01 + 0.2838836398812E+01 0.2804881644036E+01 0.2785230294779E+01 0.2778903258937E+01 + 0.2785069393544E+01 + 19 505 + 1 0 -1 -5366.16027 1 s core state + 2 0 -1 -894.94060 2 s core state + 2 1 1 -860.61248 2 p core state + 2 1 -2 -860.61248 2 p core state + 3 0 -1 -208.81808 3 s core state + 3 1 1 -192.58274 3 p core state + 3 1 -2 -192.58274 3 p core state + 3 2 2 -162.19225 3 d core state + 3 2 -3 -162.19225 3 d core state + 4 0 -1 -43.32948 4 s core state + 4 1 1 -36.33021 4 p core state + 4 1 -2 -36.33021 4 p core state + 4 2 2 -23.43700 4 d core state + 4 2 -3 -23.43700 4 d core state + 4 3 3 -6.14774 4 f semi-core state + 4 3 -4 -6.14774 4 f semi-core state + 5 0 -1 -5.40446 5 s semi-core state + 5 1 1 -3.18058 5 p semi-core state + 5 1 -2 -3.18058 5 p semi-core state + 0.1036091445995E-02 0.1084673135534E-02 0.1135598547569E-02 0.1188982463327E-02 + 0.1244945328705E-02 0.1303613534909E-02 0.1365119713034E-02 0.1429603043256E-02 + 0.1497209579386E-02 0.1568092589521E-02 0.1642412913624E-02 0.1720339338843E-02 + 0.1802048993475E-02 0.1887727760486E-02 0.1977570711566E-02 0.2071782562737E-02 + 0.2170578152591E-02 0.2274182944264E-02 0.2382833552359E-02 0.2496778296018E-02 + 0.2616277779475E-02 0.2741605501441E-02 0.2873048494753E-02 0.3010907997796E-02 + 0.3155500159278E-02 0.3307156778010E-02 0.3466226079433E-02 0.3633073530726E-02 + 0.3808082696398E-02 0.3991656136393E-02 0.4184216348803E-02 0.4386206759423E-02 + 0.4598092760470E-02 0.4820362800901E-02 0.5053529530913E-02 0.5298131003301E-02 + 0.5554731934510E-02 0.5823925028349E-02 0.6106332365471E-02 0.6402606861908E-02 + 0.6713433800071E-02 0.7039532435832E-02 0.7381657685472E-02 0.7740601896454E-02 + 0.8117196706194E-02 0.8512314993208E-02 0.8926872925224E-02 0.9361832109077E-02 + 0.9818201847456E-02 0.1029704150780E-01 0.1079946300894E-01 0.1132663343132E-01 + 0.1187977775693E-01 0.1246018174549E-01 0.1306919495350E-01 0.1370823390339E-01 + 0.1437878541021E-01 0.1508241007361E-01 0.1582074594350E-01 0.1659551236776E-01 + 0.1740851403138E-01 0.1826164519621E-01 0.1915689415163E-01 0.2009634788626E-01 + 0.2108219699212E-01 0.2211674081235E-01 0.2320239284500E-01 0.2434168641538E-01 + 0.2553728063034E-01 0.2679196662864E-01 0.2810867414199E-01 0.2949047838226E-01 + 0.3094060727108E-01 0.3246244902885E-01 0.3405956014105E-01 0.3573567372055E-01 + 0.3749470828567E-01 0.3934077697471E-01 0.4127819721844E-01 0.4331150089364E-01 + 0.4544544498131E-01 0.4768502275484E-01 0.5003547552441E-01 0.5250230496527E-01 + 0.5509128605901E-01 0.5780848067821E-01 0.6066025184651E-01 0.6365327870767E-01 + 0.6679457223890E-01 0.7009149174544E-01 0.7355176217523E-01 0.7718349229455E-01 + 0.8099519376722E-01 0.8499580118260E-01 0.8919469307926E-01 0.9360171401407E-01 + 0.9822719772854E-01 0.1030584446289E+00 0.1081048518184E+00 0.1133981004743E+00 + 0.1189503651804E+00 0.1247744234990E+00 0.1308836858973E+00 0.1372922271604E+00 + 0.1440148193649E+00 0.1510669664914E+00 0.1584649407559E+00 0.1662258207468E+00 + 0.1743675314573E+00 0.1829088863059E+00 0.1918696312456E+00 0.2012704910653E+00 + 0.2111332179907E+00 0.2214806427018E+00 0.2322810788367E+00 0.2435924141539E+00 + 0.2554517807592E+00 0.2678860771230E+00 0.2809235306829E+00 0.2945937637215E+00 + 0.3089278625148E+00 0.3239584499139E+00 0.3397197615313E+00 0.3562477257100E+00 + 0.3735800474643E+00 0.3917562965884E+00 0.4107297688896E+00 0.4306098328990E+00 + 0.4514440457470E+00 0.4732790111462E+00 0.4961636296795E+00 0.5201492125165E+00 + 0.5452896007728E+00 0.5716412907928E+00 0.5992635656497E+00 0.6281330951782E+00 + 0.6583494366762E+00 0.6900019441399E+00 0.7231604062945E+00 0.7578980424511E+00 + 0.7942916721308E+00 0.8324218930999E+00 0.8723420379212E+00 0.9140560168786E+00 + 0.9577322876177E+00 0.1003465682237E+01 0.1051355681320E+01 0.1101506643453E+01 + 0.1154028046117E+01 0.1208893003204E+01 0.1266298556958E+01 0.1326380424276E+01 + 0.1389267317171E+01 0.1455094241031E+01 0.1523937911590E+01 0.1595868923175E+01 + 0.1671117681213E+01 0.1749842809105E+01 0.1832210662338E+01 0.1918294397258E+01 + 0.2008232331635E+01 0.2102277783569E+01 0.2200626055773E+01 0.2303443629014E+01 + 0.2410777536675E+01 0.2522960065391E+01 0.2640220327706E+01 0.2762741486581E+01 + 0.2890612506339E+01 0.3024200540749E+01 0.3163774132062E+01 0.3309444136365E+01 + 0.3461489580835E+01 0.3620267522648E+01 0.3785997025710E+01 0.3958853261476E+01 + 0.4139276127511E+01 0.4327530949752E+01 0.4523793521376E+01 0.4728552502922E+01 + 0.4942060475799E+01 0.5164586369886E+01 0.5396643158989E+01 0.5638392259626E+01 + 0.5890319952885E+01 0.6152829213084E+01 0.6426161936478E+01 0.6710932086274E+01 + 0.7007282351393E+01 0.7315869834636E+01 0.7636969625132E+01 0.7971104898969E+01 + 0.8318680538980E+01 0.8680169846043E+01 0.9056032307527E+01 0.9446757136304E+01 + 0.9852793410959E+01 0.1027470937770E+02 0.1071287300445E+02 0.1116796702729E+02 + 0.1164030871059E+02 0.1213058883027E+02 0.1263925258037E+02 0.1316679504693E+02 + 0.1371385481156E+02 0.1428085875749E+02 0.1486837357854E+02 0.1547699038619E+02 + 0.1610719086090E+02 0.1675945350626E+02 0.1743440119373E+02 0.1813253136718E+02 + 0.1885434034381E+02 0.1960026192103E+02 0.2037086813047E+02 0.2116662306265E+02 + 0.2198797570251E+02 0.2283535825816E+02 0.2370918265580E+02 0.2460983703309E+02 + 0.2553767367705E+02 0.2649284109463E+02 0.2747621011270E+02 0.2848719533567E+02 + 0.2952698981969E+02 0.3059478374077E+02 0.3169154420927E+02 0.3281669024929E+02 + 0.3397064143999E+02 0.3515328042362E+02 0.3636415374300E+02 0.3760362371074E+02 + 0.3887093815944E+02 0.4016597967928E+02 0.4148845778485E+02 0.4283752071504E+02 + 0.4421273388371E+02 0.4561350958003E+02 0.4703894961410E+02 0.4848788827534E+02 + 0.4995960595730E+02 0.5145291447903E+02 0.5296651158049E+02 0.5449898844621E+02 + 0.5604873348485E+02 0.5761422753228E+02 0.5919368744277E+02 0.6078518980943E+02 + 0.6238668081747E+02 0.6399597351871E+02 0.6561074641191E+02 0.6722854325659E+02 + 0.6884677405240E+02 0.7046271712004E+02 0.7207352222382E+02 0.7367621467944E+02 + 0.7526769995067E+02 0.7684475230410E+02 0.7840408525596E+02 0.7994229502347E+02 + 0.8145589057570E+02 0.8294131967183E+02 0.8439497422853E+02 0.8581318098105E+02 + 0.8719229198470E+02 0.8852861368726E+02 0.8981860930719E+02 0.9105843626601E+02 + 0.9224445804903E+02 0.9337352167895E+02 0.9444222838383E+02 0.9544712903575E+02 + 0.9638522177286E+02 0.9725370578181E+02 0.9805000948685E+02 0.9877182296058E+02 + 0.9941714257695E+02 0.9998438554774E+02 0.1004720677117E+03 0.1008792908067E+03 + 0.1012055896424E+03 0.1014510440830E+03 0.1016162457575E+03 0.1017021624390E+03 + 0.1017103944534E+03 0.1016431159612E+03 0.1015030907003E+03 0.1012936836308E+03 + 0.1010193548027E+03 0.1006847212892E+03 0.1002948847490E+03 0.9985596746059E+02 + 0.9937502117565E+02 0.9885876879082E+02 0.9831559999803E+02 0.9775365762686E+02 + 0.9718210020078E+02 0.9661024928621E+02 0.9604822577441E+02 0.9550596586797E+02 + 0.9499425819754E+02 0.9452400586742E+02 0.9410551648543E+02 0.9374949393319E+02 + 0.9346687848176E+02 0.9326815290971E+02 0.9316341840649E+02 0.9316227932825E+02 + 0.9327371855043E+02 0.9350596715881E+02 0.9386637221324E+02 0.9436126535249E+02 + 0.9499583610646E+02 0.9577401215427E+02 0.9669825772579E+02 0.9776908901129E+02 + 0.9898619093847E+02 0.1003473082417E+03 0.1018483852021E+03 0.1034824567059E+03 + 0.1052425357480E+03 0.1071188734518E+03 0.1090989910006E+03 0.1111706475534E+03 + 0.1133172181206E+03 0.1155232715982E+03 0.1177691732584E+03 0.1200362594476E+03 + 0.1223030857352E+03 0.1245479267040E+03 0.1267486024477E+03 0.1288821226225E+03 + 0.1309248669208E+03 0.1328554330914E+03 0.1346515343386E+03 0.1362917354203E+03 + 0.1377563183418E+03 0.1390272764054E+03 0.1400886714465E+03 0.1409269798914E+03 + 0.1415314169203E+03 0.1418942288679E+03 0.1420109452088E+03 0.1418809507069E+03 + 0.1415080920783E+03 0.1408984172942E+03 0.1400624018884E+03 0.1390151068887E+03 + 0.1377764089760E+03 0.1363667544540E+03 0.1348127606975E+03 0.1331435052272E+03 + 0.1313892374741E+03 0.1295850758187E+03 0.1277646353936E+03 0.1259649823078E+03 + 0.1242218979853E+03 0.1225708137778E+03 0.1210467137096E+03 0.1196813004892E+03 + 0.1185038390914E+03 0.1175399543298E+03 0.1168105170862E+03 0.1163313357209E+03 + 0.1161129504871E+03 0.1161592767958E+03 0.1164676948596E+03 0.1170288216215E+03 + 0.1178264633586E+03 0.1188377345923E+03 0.1200333492742E+03 0.1213757757079E+03 + 0.1228248577103E+03 0.1243359368424E+03 0.1258602762822E+03 0.1273452565898E+03 + 0.1287364761148E+03 0.1299832482582E+03 0.1310331118351E+03 0.1318365896987E+03 + 0.1323503524919E+03 0.1325341948887E+03 0.1323565100334E+03 0.1317927692427E+03 + 0.1308270665976E+03 0.1294533149204E+03 0.1276758225072E+03 0.1255076118418E+03 + 0.1229729174483E+03 0.1201060641052E+03 0.1169495219547E+03 0.1135555562204E+03 + 0.1099792849535E+03 0.1062833468604E+03 0.1025360533657E+03 0.9880615345759E+02 + 0.9516244953943E+02 0.9167188652690E+02 0.8839767072261E+02 0.8539748745222E+02 + 0.8272200135400E+02 0.8041305849744E+02 0.7850121189453E+02 0.7700728710647E+02 + 0.7594037675719E+02 0.7529524307077E+02 0.7505910626795E+02 0.7520501511195E+02 + 0.7569423682430E+02 0.7648360425006E+02 0.7751457462797E+02 0.7873062235074E+02 + 0.8006319733670E+02 0.8144692237202E+02 0.8281343828438E+02 0.8409437423111E+02 + 0.8522854140375E+02 0.8615737163000E+02 0.8682931552416E+02 0.8720098369636E+02 + 0.8723823597667E+02 0.8691660972946E+02 0.8622071120042E+02 0.8514574227897E+02 + 0.8369657858485E+02 0.8188716928807E+02 0.7973971249770E+02 0.7728362903218E+02 + 0.7455526455603E+02 0.7159724324227E+02 0.6845138246224E+02 0.6516341040369E+02 + 0.6177993594066E+02 0.5835062141363E+02 0.5491700328472E+02 0.5152086107061E+02 + 0.4820207507821E+02 0.4499245032799E+02 0.4192171019143E+02 0.3901261478593E+02 + 0.3628311527177E+02 0.3374565621620E+02 0.3140691776349E+02 0.2926910131416E+02 + 0.2732972598556E+02 0.2558183615499E+02 0.2401519884669E+02 0.2261722411102E+02 + 0.2137275650990E+02 0.2026520783577E+02 0.1927719083326E+02 0.1839110865918E+02 + 0.1758969047340E+02 0.1685646087305E+02 0.1617613225919E+02 0.1553479814313E+02 + 0.1492039134364E+02 0.1432278086514E+02 0.1373361577992E+02 0.1314653008144E+02 + 0.1255704050714E+02 0.1196235214704E+02 0.1136126243718E+02 0.1075392984383E+02 + 0.1014166558229E+02 0.9526714234397E+01 0.8912009640542E+01 0.8300963732617E+01 + 0.7697263384503E+01 0.7104685454692E+01 0.6526924572959E+01 0.5967462120510E+01 + 0.5429460175041E+01 0.4915676037373E+01 0.4428403641855E+01 0.3969439968829E+01 + 0.3540078894877E+01 0.3141114634575E+01 0.2772863117433E+01 0.2435199476377E+01 + 0.2127602513822E+01 0.1849205679049E+01 0.1598851789426E+01 0.1375137670440E+01 + 0.1176496529627E+01 0.1001213347056E+01 0.8475093718041E+00 0.7135602873433E+00 + 0.5975459398289E+00 0.4976821419082E+00 0.4122475629664E+00 0.3396043781990E+00 + 0.2782131303566E+00 0.2266491872781E+00 0.1836031900326E+00 0.1478871833508E+00 + 0.1184340827987E+00 0.9429433792906E-01 0.7463086188421E-01 0.5871189801862E-01 + 0.4590372023947E-01 0.3567951860383E-01 0.2757960926186E-01 0.2119786226693E-01 + 0.1619816828197E-01 diff --git a/src/SingleSite/SingleSiteScattering.cpp b/src/SingleSite/SingleSiteScattering.cpp index b5cfa8b7a..e55b31e73 100644 --- a/src/SingleSite/SingleSiteScattering.cpp +++ b/src/SingleSite/SingleSiteScattering.cpp @@ -106,7 +106,7 @@ void calculateSingleScattererSolution(LSMSSystemParameters &lsms, AtomData &atom solution.fz.resize(iprpts, 2, kmymax); solution.fj.resize(iprpts, 2, kmymax); - printf("iprpts = %d ir = %d\n", iprpts, ir); + // printf("iprpts = %d ir = %d\n", iprpts, ir); // exit(1); single_scatterer_rel_(&energy, &psq, &atom.lmax, &kmymax, @@ -119,6 +119,7 @@ void calculateSingleScattererSolution(LSMSSystemParameters &lsms, AtomData &atom &solution.nuz[0], &solution.indz(0,0), &iflag, &soscal, &iprpts, &lsms.global.iprint,lsms.global.istop,32); +/* for(int j1=0; j1