[Smeagol-discuss] About MD calculation
Jian ZHOU
zhoujian at gmail.com
Tue Apr 7 13:34:08 IST 2009
hi,
First of all, I think this is not a siesta maillist.
Then, it seems that your energy does not converge at all. And your force and
stress are very large. Your unit cell is also a little strange. So, I think
you should check your structure first.
btw, are you sure the pure B-nanotube stable? not BN nanotube?
Best wishes,
jian
On Tue, Apr 7, 2009 at 8:02 PM, plijun2 <plijun2 at student.cityu.edu.hk>wrote:
> Dear All,
> when i do some MD calculations about boron nanotubes, my problem is the
> calculation stopped after two MD runs though task being in process, If
> you know where is the problems, how can I avoid them? the parameters is
> below:
> SystemLabel B-tube
> NumberOfAtoms 80
> NumberOfSpecies 1
> %block ChemicalSpeciesLabel
> 1 5 B
> %endblock ChemicalSpeciesLabel
> %block kgrid_Monkhorst_Pack
> 1 0 0 0.0
> 0 1 0 0.0
> 0 0 4 0.5
> %endblock kgrid_Monkhorst_Pack
> MeshCutoff 210. Ry
> SolutionMethod diagon
> %block PAO.BasisSizes
> B DZP
> %endblock PAO.BasisSizes
> MaxSCFIterations 200
> DM.MixingWeight 0.08
> DM.Tolerance 1.d-4
> DM.NumberPulay 1
> ElectronicTemperature 300.0 K
> DM.UseSaveDM true
> UseSaveData true
> xc.functional LDA
> xc.authors CA
> MD.TypeOfRun Nose
> MD.InitialTimeStep 1
> MD.FinalTimeStep 200
> MD.LengthTimeStep 3.0 fs
> MD.InitialTemperature 300.0 K
> MD.TargetTemperature 900.0 K
> MD.NoseMass 10.0 Ry*fs**2
> WriteMullikenPop 1
> SaveDeltaRho .true.
> SaveRho .true.
> SaveElectrostaticPotential .true.
> SaveTotalPotential .true.
>
> the tail of the *.out file is below:
> Begin MD step = 2
> ==============================
>
> superc: Internal auxiliary supercell: 1 x 1 x 4 = 4
> superc: Number of atoms, orbitals, and projectors: 320 4160 5120
>
> outcell: Unit cell vectors (Ang):
> 23.380649 0.228425 0.000028
> 0.228425 23.015078 -0.000118
> 0.000006 -0.000026 5.079132
>
> outcell: Cell vector modules (Ang) : 23.381765 23.016212 5.079132
> outcell: Cell angles (23,13,12) (deg): 90.0006 89.9999 88.8716
> outcell: Cell volume (Ang**3) : 2732.8538
>
> siesta: iscf Eharris(eV) E_KS(eV) FreeEng(eV) dDmax Ef(eV)
> siesta: 1 -7051.9838 -7050.6406 -7050.6408 0.0368 -4.2754
> siesta: 2 -7211.8922 -7043.2105 -7043.2539 0.4719 -4.7185
> siesta: 3 -7517.7718 -7017.4702 -7017.5268 0.9271 -5.1648
> siesta: 4 -8487.9698 -6948.9505 -6948.9798 1.8118 -6.8120
> siesta: 5 -9970.2379 -6787.3928 -6787.4027 1.9907-12.3450
> siesta: 6 -9238.3203 -6702.1485 -6702.1926 1.1823-10.2885
> siesta: 7 -8487.7483 -6665.2915 -6665.3395 1.0848 -7.1378
> siesta: 8 -10053.7822 -6528.6350 -6528.6907 2.6354-12.4869
> siesta: 9 -9522.1997 -6430.2782 -6430.3076 0.8784-11.9571
> siesta: 10 -8142.5899 -6445.8784 -6445.9326 0.7887 -6.4984
> siesta: 11 -9609.2443 -6360.4433 -6360.5007 2.7786-11.0113
> siesta: 12 -9905.1171 -6235.2577 -6235.2667 0.7563-13.0335
> siesta: 13 -8406.6404 -6263.0145 -6263.0605 1.2971 -7.7136
> siesta: 14 -9402.4740 -6203.3521 -6203.4003 3.9235-10.2040
> siesta: 15 -10286.2003 -6054.9228 -6054.9603 1.2908-14.8732
> siesta: 16 -8664.2461 -6095.9931 -6096.0067 1.5461 -8.7113
> siesta: 17 -8752.0660 -6078.6655 -6078.6988 2.3421 -8.2695
> siesta: 18 -10382.3825 -5938.9577 -5938.9996 1.3597-15.0789
> siesta: 19 -8839.0386 -5974.2922 -5974.3135 1.5707 -9.5967
> siesta: 20 -8492.4826 -5982.0717 -5982.1096 1.3146 -7.5641
> siesta: 21 -10040.9585 -5881.6632 -5881.7100 1.2289-13.9986
> siesta: 22 -8918.5609 -5902.2913 -5902.3343 1.2351 -9.9892
> siesta: 23 -8570.6836 -5917.3222 -5917.3667 1.2807 -8.0055
> siesta: 24 -10239.2188 -5814.4202 -5814.4479 1.7368-14.3132
> siesta: 25 -9069.4759 -5826.5595 -5826.6002 1.3656-10.7125
> siesta: 26 -8317.3863 -5863.6814 -5863.7054 0.7638 -7.4809
> siesta: 27 -9910.8094 -5795.2108 -5795.2854 1.9631-12.9554
> siesta: 28 -9315.0563 -5778.5191 -5778.5462 1.1572-11.5062
> siesta: 29 -8437.6514 -5820.4458 -5820.4890 1.0494 -7.7454
> siesta: 30 -9869.7187 -5765.1493 -5765.1930 2.2396-12.6481
> siesta: 31 -9509.0257 -5726.9072 -5726.9652 0.7567-12.3101
> siesta: 32 -8443.3986 -5783.1722 -5783.2106 1.0727 -8.0071
> siesta: 33 -9579.2344 -5751.8445 -5751.8919 2.4912-11.3123
> siesta: 34 -9927.1264 -5673.2859 -5673.3125 0.7882-13.3716
> siesta: 35 -8512.0068 -5740.6950 -5740.7455 1.0929 -8.4899
> siesta: 36 -9154.6754 -5735.6999 -5735.7420 2.5554 -9.4506
> siesta: 37 -10458.4781 -5613.6417 -5613.6722 1.5036-15.6415
> siesta: 38 -8742.4789 -5686.0833 -5686.1391 1.6101 -9.2708
> siesta: 39 -8569.4632 -5711.1614 -5711.1774 1.6164 -7.9883
> siesta: 40 -10088.9488 -5627.1143 -5627.1633 1.0975-14.2494
> siesta: 41 -8891.9769 -5672.6126 -5672.6331 1.2706-10.0050
> siesta: 42 -8573.9096 -5704.5143 -5704.5323 1.2055 -7.9981
> siesta: 43 -10148.6987 -5623.4262 -5623.4638 1.6959-14.2753
> siesta: 44 -9030.5692 -5654.0008 -5654.0411 1.2873-10.5518
> siesta: 45 -8397.0494 -5700.3125 -5700.3501 0.8491 -7.7440
> siesta: 46 -10020.0781 -5635.6614 -5635.7207 1.6199-13.3542
> siesta: 47 -9231.9917 -5641.0768 -5641.1264 1.0854-11.2867
> siesta: 48 -8437.3208 -5691.9775 -5691.9852 0.7581 -7.6764
> siesta: 49 -9876.3726 -5645.0079 -5645.0733 2.2266-12.8718
> siesta: 50 -9459.7964 -5622.4120 -5622.4566 0.8199-12.0511
> siesta: 51 -8444.9542 -5684.3599 -5684.3859 1.0647 -8.0058
> siesta: 52 -9678.9889 -5654.6371 -5654.6810 2.3202-11.6211
> siesta: 53 -9804.6276 -5594.9503 -5595.0130 0.7755-13.1750
> siesta: 54 -8438.6655 -5669.9969 -5670.0403 1.0148 -8.1112
> siesta: 55 -9315.5151 -5663.2986 -5663.3324 2.5546-10.1888
> siesta: 56 -10304.6819 -5558.4286 -5558.4344 1.1718-14.9908
> siesta: 57 -8661.4703 -5637.6553 -5637.6639 1.3469 -9.0713
> siesta: 58 -8693.3974 -5660.9386 -5660.9734 1.8276 -8.0799
> siesta: 59 -10273.8334 -5565.5763 -5565.6223 1.2809-14.8710
> siesta: 60 -8875.2551 -5622.6722 -5622.7029 1.4816 -9.7118
> siesta: 61 -8485.4394 -5660.3015 -5660.3456 1.2566 -7.8253
> siesta: 62 -10100.6419 -5583.1351 -5583.1817 1.2248-14.0941
> siesta: 63 -8963.3176 -5624.1649 -5624.2079 1.2065-10.3598
> siesta: 64 -8496.1470 -5665.3158 -5665.3445 1.0190 -7.8931
> siesta: 65 -10121.5453 -5592.2553 -5592.2829 1.7926-13.9691
> siesta: 66 -9139.8110 -5613.7295 -5613.7435 1.3358-10.9342
> siesta: 67 -8340.5437 -5668.4117 -5668.4501 0.6632 -7.5965
> siesta: 68 -9820.0768 -5624.0804 -5624.1286 1.7509-12.5305
> siesta: 69 -9374.0870 -5613.5330 -5613.5575 0.9935-11.7692
> siesta: 70 -8500.5623 -5668.6437 -5668.6668 0.9460 -7.9523
> siesta: 71 -9769.5023 -5633.6301 -5633.6570 2.4725-12.2741
> siesta: 72 -9664.0188 -5589.5647 -5589.5980 0.8247-12.6557
> siesta: 73 -8446.1700 -5660.8125 -5660.8473 1.0900 -8.1173
> siesta: 74 -9462.8175 -5646.0491 -5646.1041 2.3941-10.7612
> siesta: 75 -10084.0359 -5561.8936 -5561.9381 0.8592-14.0664
> siesta: 76 -8632.6975 -5635.9174 -5635.9312 1.2849 -8.9571
> siesta: 77 -8870.5318 -5652.3468 -5652.3818 2.2110 -8.6815
> siesta: 78 -10441.0242 -5542.7073 -5542.7399 1.3971-15.4385
> siesta: 79 -8818.0001 -5611.8445 -5611.9165 1.7031 -9.5951
> siesta: 80 -8490.2921 -5648.5548 -5648.5805 1.1770 -7.8317
> siesta: 81 -9994.4247 -5580.6352 -5580.6845 1.1858-13.7552
> siesta: 82 -9012.7419 -5615.9920 -5616.0282 1.1936-10.4131
> siesta: 83 -8501.0126 -5659.0213 -5659.0415 1.0583 -8.0669
> siesta: 84 -10175.1166 -5585.6331 -5585.6771 1.7894-14.0007
> siesta: 85 -9163.0081 -5603.7544 -5603.7762 1.2157-11.0853
> siesta: 86 -8333.7909 -5661.4548 -5661.4921 0.7171 -7.4583
> siesta: 87 -9812.5028 -5619.7493 -5619.7946 2.0637-12.5758
> siesta: 88 -9460.6295 -5600.3247 -5600.3547 0.9102-11.9832
> siesta: 89 -8442.0079 -5662.7298 -5662.7704 1.1241 -7.8775
> siesta: 90 -9770.0714 -5628.0612 -5628.0880 2.3146-12.1339
> siesta: 91 -9689.3057 -5582.1804 -5582.2216 0.8544-12.9087
> siesta: 92 -8448.8938 -5655.0787 -5655.1090 0.9965 -8.1075
> siesta: 93 -9352.9123 -5647.7952 -5647.8545 2.5212-10.3005
> siesta: 94 -10216.8661 -5551.4901 -5551.5222 0.8879-14.6251
> siesta: 95 -8598.8635 -5632.4368 -5632.4679 1.2999 -8.8197
> siesta: 96 -8800.8676 -5651.4554 -5651.4929 2.0217 -8.5035
> siesta: 97 -10360.8276 -5547.9517 -5548.0024 1.3200-15.2337
> siesta: 98 -8829.5452 -5613.8166 -5613.8520 1.5192 -9.5220
> siesta: 99 -8513.6081 -5649.5709 -5649.6270 1.3969 -7.8735
> siesta: 100 -10094.1715 -5572.0353 -5572.0654 1.1595-14.0921
> siesta: 101 -8950.5556 -5615.8910 -5615.9298 1.1935-10.3083
> siesta: 102 -8519.4557 -5655.7401 -5655.7739 1.0787 -7.8919
> siesta: 103 -10138.3168 -5581.2062 -5581.2631 1.7966-14.1184
> siesta: 104 -9111.8582 -5607.1235 -5607.1669 1.3913-10.8401
> siesta: 105 -8344.6234 -5661.2576 -5661.2731 0.7307 -7.6292
> siesta: 106 -9854.8474 -5614.2426 -5614.2796 1.7136-12.6909
> siesta: 107 -9338.3898 -5608.4909 -5608.5098 1.0269-11.6893
> siesta: 108 -8496.4376 -5662.8954 -5662.9439 0.9528 -7.9510
> siesta: 109 -9839.6555 -5623.3793 -5623.4281 2.4983-12.5959
> siesta: 110 -9618.0945 -5584.5849 -5584.6215 0.8310-12.5922
> siesta: 111 -8404.5112 -5658.2189 -5658.2645 1.0881 -7.9402
> siesta: 112 -9531.0784 -5639.7592 -5639.8007 2.3884-11.0256
> siesta: 113 -9993.6055 -5564.0348 -5564.0843 0.7799-13.7309
> siesta: 114 -8618.1172 -5636.4979 -5636.5445 1.3880 -8.9027
> siesta: 115 -8944.3640 -5649.4603 -5649.4826 2.3275 -8.9307
> siesta: 116 -10427.0081 -5541.1351 -5541.1681 1.3490-15.4498
> siesta: 117 -8811.6590 -5610.8592 -5610.8936 1.7509 -9.5572
> siesta: 118 -8512.2832 -5646.0211 -5646.0650 1.2241 -7.9360
> siesta: 119 -9990.4845 -5578.9927 -5579.0448 1.2197-13.7620
> siesta: 120 -9003.0016 -5615.2493 -5615.2994 1.1505-10.3402
> siesta: 121 -8496.6365 -5658.4066 -5658.4348 1.1237 -8.0059
> siesta: 122 -10194.9398 -5582.9178 -5582.9591 1.8324-13.9858
> siesta: 123 -9164.9489 -5601.7604 -5601.8090 1.1487-11.2049
> siesta: 124 -8250.2455 -5664.0925 -5664.1309 0.6696 -7.2722
> siesta: 125 -9770.0737 -5624.9340 -5624.9884 2.0483-12.4082
> siesta: 126 -9467.3059 -5604.4251 -5604.4628 0.9205-11.9826
> siesta: 127 -8499.1795 -5663.3066 -5663.3632 1.1328 -8.0966
> siesta: 128 -9694.2260 -5633.0959 -5633.1111 2.1704-11.9076
> siesta: 129 -9682.3252 -5588.8938 -5588.9232 0.8870-12.7905
> siesta: 130 -8511.7061 -5655.6994 -5655.7154 0.9326 -8.3077
> siesta: 131 -9339.1331 -5648.7871 -5648.8437 2.3368-10.3841
> siesta: 132 -10109.7778 -5563.9281 -5563.9696 0.9060-14.1235
> siesta: 133 -8665.5106 -5634.2819 -5634.3224 1.1900 -9.1929
> siesta: 134 -8800.8342 -5655.2940 -5655.3198 1.9376 -8.4173
> siesta: 135 -10364.2135 -5553.5667 -5553.6034 1.6020-15.1670
> siesta: 136 -8877.6050 -5613.0720 -5613.1265 1.6799 -9.7454
> siesta: 137 -8447.8961 -5653.3424 -5653.3783 1.1751 -7.8689
> siesta: 138 -9967.3923 -5588.2606 -5588.3004 1.2282-13.5290
> siesta: 139 -9031.8385 -5620.2823 -5620.3149 1.0557-10.5385
> siesta: 140 -8534.8153 -5661.0192 -5661.0765 0.8680 -8.1744
> siesta: 141 -10110.3659 -5593.0037 -5593.0482 2.0396-13.8243
> siesta: 142 -9206.3829 -5605.4008 -5605.4137 1.2884-11.1055
> siesta: 143 -8363.2891 -5662.4557 -5662.4815 0.7471 -7.4611
> siesta: 144 -9918.0387 -5615.0808 -5615.1223 2.0175-12.8427
> siesta: 145 -9435.2649 -5597.3685 -5597.3943 0.7634-11.9433
> siesta: 146 -8420.7043 -5662.0572 -5662.0812 0.9966 -7.7666
> siesta: 147 -9747.3248 -5629.9638 -5629.9994 2.5051-11.9821
> siesta: 148 -9791.9127 -5572.8811 -5572.9067 0.8278-13.0885
> siesta: 149 -8411.5010 -5651.4348 -5651.4675 1.0861 -8.0797
> siesta: 150 -9396.2395 -5641.2163 -5641.2756 2.5596-10.4623
> siesta: 151 -10256.7472 -5540.7172 -5540.7655 1.1142-14.7729
> siesta: 152 -8644.6528 -5622.3526 -5622.4095 1.3989 -8.9542
> siesta: 153 -8745.9825 -5643.7612 -5643.7776 2.0422 -8.4550
> siesta: 154 -10345.7717 -5542.8067 -5542.8596 1.2601-15.1544
> siesta: 155 -8838.0808 -5607.7562 -5607.7976 1.4295 -9.7080
> siesta: 156 -8512.9481 -5644.1875 -5644.2096 1.2699 -7.7640
> siesta: 157 -10087.6217 -5568.4544 -5568.5240 1.2199-14.1067
> siesta: 158 -8968.1590 -5610.5704 -5610.6093 1.3025-10.2910
> siesta: 159 -8495.3279 -5653.0157 -5653.0594 1.0972 -7.9836
> siesta: 160 -10123.5319 -5580.6851 -5580.7265 1.5457-13.9417
> siesta: 161 -9116.2330 -5605.0277 -5605.0641 1.2344-10.9117
> siesta: 162 -8399.8987 -5656.6344 -5656.6694 0.7486 -7.6923
> siesta: 163 -9997.8441 -5600.5165 -5600.5653 2.0473-13.2600
> siesta: 164 -9350.7867 -5594.7649 -5594.8189 1.0836-11.7375
> siesta: 165 -8352.2045 -5658.3977 -5658.4387 1.0885 -7.5467
> siesta: 166 -9816.8144 -5619.8906 -5619.9466 2.1152-12.4275
> siesta: 167 -9558.6754 -5588.5679 -5588.6056 0.7678-12.3928
> siesta: 168 -8466.4716 -5655.5422 -5655.5794 1.1128 -8.1787
> siesta: 169 -9492.6878 -5639.6990 -5639.7414 2.3984-11.0129
> siesta: 170 -9989.1579 -5565.0925 -5565.1284 0.8116-13.6840
> siesta: 171 -8542.8955 -5640.3986 -5640.4395 1.0920 -8.6481
> siesta: 172 -9073.3683 -5647.1447 -5647.1746 2.3477 -9.1690
> siesta: 173 -10515.9076 -5529.6179 -5529.6622 1.7811-15.8275
> siesta: 174 -8792.6658 -5605.2299 -5605.2462 1.6987 -9.5061
> siesta: 175 -8482.9981 -5641.3583 -5641.3926 1.3605 -7.8837
> siesta: 176 -10029.8888 -5568.7356 -5568.7717 1.1072-14.0000
> siesta: 177 -8922.8528 -5614.4376 -5614.4794 1.2029-10.1590
> siesta: 178 -8567.0041 -5651.9414 -5651.9595 1.1367 -7.9747
> siesta: 179 -10173.6715 -5574.6622 -5574.7248 1.8152-14.2713
> siesta: 180 -9079.0915 -5604.4193 -5604.4605 1.2566-10.7918
> siesta: 181 -8324.4587 -5659.6636 -5659.6970 0.7113 -7.5972
> siesta: 182 -9886.9403 -5610.1979 -5610.2540 1.6876-12.7899
> siesta: 183 -9330.4231 -5606.0511 -5606.0953 1.0345-11.6460
> siesta: 184 -8468.4507 -5661.3789 -5661.4291 0.9341 -7.8577
> siesta: 185 -9886.0365 -5619.1852 -5619.2231 2.5593-12.7148
> siesta: 186 -9612.5213 -5580.8300 -5580.8651 0.8532-12.6494
> siesta: 187 -8370.4775 -5656.9422 -5656.9846 1.0919 -7.8323
> siesta: 188 -9575.1917 -5636.3021 -5636.3506 2.4405-11.2189
> siesta: 189 -10020.8409 -5556.9178 -5556.9509 0.7811-13.8235
> siesta: 190 -8568.9803 -5634.8735 -5634.9185 1.3820 -8.7310
> siesta: 191 -8992.3038 -5645.1658 -5645.2020 2.3696 -9.0535
> siesta: 192 -10472.4081 -5530.8721 -5530.8935 1.3187-15.6447
> siesta: 193 -8778.0448 -5607.1254 -5607.1482 1.8138 -9.4213
> siesta: 194 -8529.1101 -5640.4489 -5640.4816 1.2937 -7.9593
> siesta: 195 -9980.9570 -5572.6494 -5572.6976 1.1523-13.8323
> siesta: 196 -8971.8495 -5613.1909 -5613.2205 1.1310-10.2730
> siesta: 197 -8541.3931 -5653.0380 -5653.0662 1.1923 -7.9950
> siesta: 198 -10218.2106 -5573.5170 -5573.5667 1.7006-14.2400
> siesta: 199 -9078.7862 -5602.9174 -5602.9796 1.0555-10.8879
>
> siesta: E_KS(eV) = -5659.8460
>
> siesta: Atomic forces (eV/Ang):
> ----------------------------------------
> Tot -0.169126 0.024447 0.006474
> ----------------------------------------
> Max 17.095007
> Res 5.357858 sqrt( Sum f_i^2 / 3N )
> ----------------------------------------
> Max 17.095007 constrained
>
> Stress-tensor-Voigt (kbar): -194.33 -401.10 -622.89 -6.85
> 0.26 0.42
> Target enthalpy (eV/cell) -5659.8880
>
> siesta: Stress tensor (static) (eV/Ang**3):
> -0.121169 -0.004594 0.000215
> -0.004299 -0.250340 0.000060
> 0.000224 0.000169 -0.388761
>
> siesta: Pressure (static): 406.03298931 kBar
>
> siesta: Stress tensor (total) (eV/Ang**3):
> -0.121288 -0.004573 0.000252
> -0.004278 -0.250343 0.000053
> 0.000261 0.000162 -0.388772
>
> siesta: Pressure (total): 406.10448906 kBar
>
> mulliken: Atomic and Orbital Populations:
>
> Species: B
> Atom Qatom Qorb
> 2s 2s 2py 2pz 2px 2py 2pz 2px
> 2Pdxy 2Pdyz 2Pdz2 2Pdxz 2Pdx2-y2
> 1 3.028 0.151 0.509 0.103 0.077 0.299 0.494 0.521 0.442
> 0.068 0.094 0.094 0.073 0.103
> 2 3.068 0.115 0.517 0.188 0.073 0.284 0.487 0.494 0.458
> 0.096 0.091 0.107 0.080 0.078
> 3 3.175 0.133 0.514 0.227 0.137 0.218 0.476 0.492 0.513
> 0.080 0.104 0.098 0.102 0.083
> 4 3.121 0.185 0.562 0.324 0.079 0.145 0.408 0.542 0.498
> 0.064 0.068 0.086 0.086 0.073
> 5 2.952 0.215 0.594 0.341 0.124 0.147 0.295 0.536 0.498
> 0.032 0.048 0.050 0.054 0.018
> 6 2.921 0.097 0.524 0.284 0.164 0.147 0.421 0.411 0.462
> 0.068 0.080 0.090 0.107 0.066
> 7 2.916 0.112 0.480 0.225 0.071 0.273 0.442 0.464 0.426
> 0.083 0.088 0.104 0.084 0.065
> 8 2.927 0.112 0.510 0.153 0.071 0.268 0.461 0.480 0.452
> 0.072 0.086 0.099 0.068 0.094
> 9 2.939 0.127 0.507 0.054 0.112 0.273 0.492 0.489 0.455
> 0.067 0.103 0.092 0.079 0.089
> 10 2.932 0.140 0.458 0.185 0.064 0.236 0.463 0.482 0.472
> 0.091 0.086 0.104 0.083 0.068
> 11 2.901 0.125 0.492 0.299 0.046 0.147 0.450 0.473 0.453
> 0.071 0.073 0.102 0.094 0.076
> 12 2.984 0.155 0.640 0.339 0.238 0.179 0.281 0.446 0.490
> 0.043 0.039 0.044 0.068 0.023
> 13 3.133 0.183 0.560 0.337 0.094 0.115 0.402 0.537 0.507
> 0.069 0.072 0.094 0.089 0.075
> 14 3.147 0.129 0.498 0.241 0.082 0.219 0.517 0.492 0.514
> 0.091 0.086 0.105 0.095 0.076
> 15 3.123 0.163 0.497 0.067 0.131 0.286 0.531 0.495 0.484
> 0.085 0.113 0.101 0.086 0.086
> 16 3.019 0.161 0.504 0.086 0.072 0.308 0.514 0.511 0.435
> 0.069 0.087 0.099 0.073 0.099
> 17 3.060 0.122 0.512 0.184 0.053 0.275 0.484 0.519 0.464
> 0.098 0.092 0.103 0.079 0.077
> 18 3.173 0.137 0.512 0.225 0.147 0.216 0.476 0.485 0.514
> 0.079 0.103 0.097 0.100 0.081
> 19 3.122 0.174 0.569 0.314 0.090 0.135 0.417 0.537 0.509
> 0.063 0.069 0.087 0.085 0.072
> 20 2.948 0.220 0.591 0.341 0.114 0.155 0.293 0.541 0.488
> 0.032 0.048 0.050 0.054 0.019
> 21 2.915 0.101 0.524 0.274 0.147 0.146 0.432 0.423 0.465
> 0.069 0.078 0.087 0.105 0.065
> 22 2.913 0.099 0.480 0.203 0.076 0.248 0.462 0.455 0.457
> 0.088 0.087 0.111 0.082 0.067
> 23 2.920 0.117 0.507 0.147 0.058 0.268 0.462 0.491 0.453
> 0.072 0.085 0.098 0.069 0.094
> 24 2.945 0.133 0.507 0.086 0.132 0.282 0.478 0.472 0.440
> 0.064 0.099 0.090 0.075 0.088
> 25 2.935 0.140 0.459 0.200 0.065 0.238 0.458 0.476 0.469
> 0.093 0.085 0.103 0.081 0.069
> 26 2.921 0.116 0.502 0.292 0.063 0.178 0.457 0.458 0.435
> 0.073 0.074 0.106 0.095 0.074
> 27 2.966 0.146 0.642 0.334 0.233 0.156 0.285 0.448 0.503
> 0.043 0.039 0.045 0.068 0.022
> 28 3.129 0.195 0.551 0.353 0.080 0.106 0.392 0.549 0.509
> 0.069 0.070 0.090 0.090 0.074
> 29 3.159 0.146 0.499 0.273 0.061 0.218 0.483 0.515 0.502
> 0.091 0.093 0.101 0.104 0.074
> 30 3.127 0.173 0.496 0.083 0.135 0.290 0.520 0.486 0.474
> 0.086 0.111 0.102 0.086 0.085
> 31 3.031 0.131 0.524 0.070 0.151 0.300 0.511 0.469 0.444
> 0.069 0.103 0.097 0.076 0.086
> 32 3.127 0.148 0.488 0.217 0.072 0.264 0.485 0.508 0.492
> 0.093 0.094 0.104 0.090 0.072
> 33 3.196 0.122 0.532 0.319 0.071 0.174 0.497 0.506 0.496
> 0.085 0.084 0.122 0.098 0.090
> 34 3.041 0.144 0.643 0.351 0.221 0.163 0.310 0.457 0.504
> 0.050 0.045 0.052 0.074 0.029
> 35 2.900 0.178 0.540 0.312 0.077 0.112 0.350 0.516 0.471
> 0.063 0.062 0.077 0.078 0.065
> 36 2.911 0.134 0.461 0.248 0.074 0.189 0.449 0.456 0.475
> 0.082 0.084 0.100 0.095 0.065
> 37 2.970 0.162 0.467 0.064 0.124 0.271 0.498 0.475 0.461
> 0.081 0.108 0.096 0.083 0.081
> 38 2.941 0.147 0.497 0.088 0.057 0.274 0.489 0.507 0.448
> 0.069 0.090 0.095 0.072 0.108
> 39 2.932 0.123 0.483 0.168 0.065 0.261 0.467 0.484 0.440
> 0.101 0.091 0.101 0.076 0.074
> 40 2.906 0.125 0.468 0.194 0.144 0.195 0.455 0.431 0.472
> 0.070 0.094 0.090 0.094 0.073
> 41 2.913 0.165 0.564 0.292 0.085 0.132 0.383 0.506 0.462
> 0.052 0.060 0.076 0.080 0.057
> 42 3.011 0.226 0.609 0.350 0.124 0.165 0.297 0.546 0.486
> 0.033 0.047 0.051 0.053 0.024
> 43 3.193 0.127 0.551 0.341 0.155 0.178 0.442 0.469 0.488
> 0.073 0.088 0.092 0.118 0.072
> 44 3.098 0.104 0.524 0.242 0.063 0.281 0.471 0.507 0.466
> 0.084 0.091 0.108 0.086 0.071
> 45 3.015 0.120 0.529 0.175 0.076 0.297 0.465 0.501 0.439
> 0.073 0.086 0.098 0.069 0.086
> 46 3.040 0.143 0.515 0.049 0.137 0.301 0.528 0.482 0.449
> 0.070 0.102 0.097 0.077 0.088
> 47 3.124 0.140 0.488 0.201 0.069 0.249 0.499 0.512 0.503
> 0.096 0.095 0.108 0.090 0.073
> 48 3.192 0.133 0.529 0.323 0.054 0.172 0.490 0.528 0.492
> 0.083 0.084 0.114 0.101 0.089
> 49 3.040 0.148 0.641 0.350 0.206 0.152 0.313 0.466 0.511
> 0.051 0.045 0.052 0.075 0.029
> 50 2.891 0.182 0.537 0.318 0.077 0.109 0.346 0.512 0.473
> 0.061 0.061 0.077 0.075 0.065
> 51 2.903 0.133 0.466 0.255 0.058 0.199 0.445 0.465 0.471
> 0.080 0.079 0.095 0.091 0.065
> 52 2.975 0.172 0.464 0.071 0.122 0.283 0.500 0.473 0.450
> 0.080 0.105 0.093 0.080 0.082
> 53 2.944 0.150 0.496 0.085 0.068 0.276 0.492 0.499 0.445
> 0.068 0.090 0.096 0.072 0.107
> 54 2.924 0.102 0.493 0.168 0.067 0.248 0.470 0.474 0.458
> 0.100 0.089 0.105 0.075 0.077
> 55 2.910 0.115 0.469 0.208 0.145 0.201 0.452 0.428 0.470
> 0.069 0.093 0.092 0.094 0.074
> 56 2.897 0.186 0.550 0.297 0.079 0.134 0.374 0.503 0.459
> 0.051 0.059 0.076 0.074 0.056
> 57 3.018 0.218 0.616 0.343 0.117 0.152 0.307 0.552 0.502
> 0.032 0.048 0.051 0.057 0.023
> 58 3.171 0.101 0.554 0.334 0.149 0.139 0.448 0.481 0.521
> 0.073 0.091 0.091 0.117 0.074
> 59 3.096 0.107 0.519 0.244 0.091 0.280 0.465 0.483 0.470
> 0.081 0.092 0.113 0.084 0.069
> 60 3.024 0.115 0.527 0.150 0.066 0.279 0.482 0.508 0.464
> 0.076 0.090 0.101 0.072 0.094
> 61 3.131 0.082 0.554 0.187 0.092 0.272 0.493 0.526 0.476
> 0.097 0.094 0.097 0.086 0.076
> 62 3.127 0.118 0.567 0.280 0.014 0.079 0.497 0.562 0.551
> 0.089 0.082 0.103 0.105 0.080
> 63 2.826 0.108 0.621 0.241 0.004 -0.035 0.401 0.579 0.584
> 0.056 0.061 0.072 0.074 0.059
> 64 2.914 0.072 0.528 0.211 0.072 0.166 0.479 0.481 0.485
> 0.088 0.081 0.098 0.082 0.070
> 65 2.963 0.087 0.549 0.057 0.062 0.253 0.526 0.518 0.466
> 0.078 0.092 0.094 0.074 0.106
> 66 2.861 0.104 0.506 0.052 -0.004 0.190 0.521 0.532 0.506
> 0.093 0.093 0.103 0.084 0.080
> 67 2.843 0.062 0.640 0.251 0.079 0.082 0.413 0.499 0.511
> 0.049 0.057 0.067 0.077 0.057
> 68 2.916 0.104 0.667 0.275 0.014 -0.040 0.382 0.623 0.669
> 0.043 0.051 0.049 0.064 0.016
> 69 3.050 0.098 0.548 0.103 -0.013 0.192 0.548 0.558 0.539
> 0.100 0.093 0.110 0.100 0.074
> 70 2.931 0.079 0.556 0.020 0.002 0.204 0.551 0.558 0.511
> 0.078 0.095 0.103 0.077 0.096
> 71 3.045 0.094 0.547 0.063 0.095 0.305 0.548 0.519 0.441
> 0.076 0.093 0.093 0.075 0.095
> 72 2.981 0.096 0.548 0.053 -0.002 0.207 0.544 0.551 0.526
> 0.090 0.096 0.104 0.086 0.081
> 73 3.049 0.055 0.659 0.289 0.096 0.104 0.443 0.524 0.531
> 0.055 0.065 0.075 0.088 0.065
> 74 2.851 0.101 0.654 0.279 0.002 -0.037 0.373 0.629 0.648
> 0.042 0.049 0.044 0.056 0.012
> 75 2.852 0.088 0.503 0.092 -0.010 0.168 0.523 0.515 0.515
> 0.102 0.086 0.107 0.090 0.072
> 76 2.838 0.088 0.542 0.027 0.002 0.200 0.517 0.527 0.496
> 0.074 0.093 0.102 0.076 0.095
> 77 2.941 0.090 0.510 0.145 0.069 0.222 0.489 0.507 0.468
> 0.101 0.089 0.099 0.080 0.073
> 78 2.837 0.085 0.541 0.219 0.001 0.047 0.493 0.512 0.519
> 0.080 0.075 0.092 0.098 0.074
> 79 3.044 0.149 0.622 0.299 0.014 -0.029 0.432 0.583 0.602
> 0.064 0.066 0.084 0.086 0.072
> 80 3.147 0.078 0.578 0.262 0.101 0.210 0.480 0.499 0.506
> 0.092 0.083 0.097 0.089 0.072
>
> mulliken: Qtot = 240.000
>
>
> Thank you all in advance !
> PAN Lijun
>
> _______________________________________________
> Smeagol-discuss mailing list
> Smeagol-discuss at lists.tchpc.tcd.ie
> http://lists.tchpc.tcd.ie/listinfo/smeagol-discuss
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.tchpc.tcd.ie/pipermail/smeagol-discuss/attachments/20090407/ab3829a2/attachment-0001.html
More information about the Smeagol-discuss
mailing list