&& 查看话题
加了一个电子算垂直亲和势，总是收敛失败，已经试过各种方法，求大神帮忙看看！
因为工作里要计算甲基修饰的GC碱基对的垂直亲和势，在已经优化好的结构，我加上了一个电子算单点能，基组稍微大了一些，算了好多遍，都告诉我不收敛TAT，但这个碱基对感觉在生化领域十分基本，应该不太可能算不出来，难道说是这样计算垂直亲和势的方法本身就有错？本人年龄小，见识少，能力有限，金币不够实在是不好意思TAT，请教高手帮忙看看怎么能算下来
这是log文件
------------------------------------------------------------------
#p ub3lyp/6-311++g(d,p) pop=reg scf=(maxcycle=600,tight) guess=mix
------------------------------------------------------------------
2/17=6,18=5,40=1/2;
3/5=4,6=6,7=,16=1,25=1,30=1,74=-5/1,2,3;
4/7=2,9=2,13=1/1;
5/5=2,7=600,32=2,38=5/2;
99/5=1,9=1/99;
Leave Link& & 1 at Sat Feb&&7 13:33:06 2015, MaxMem=
cpu:& &&&172.3
(Enter /opt/APP//g03/l101.exe)
-------------------
Title Card Required
-------------------
Symbolic Z-matrix:
Charge = -1 Multiplicity = 2
C& && && && && && &&&-4.95996& &0.56748& &0.00364
H& && && && && && &&&-6.04276& &0.59368& &0.00588
C& && && && && && &&&-4.21082& &1.69729&&-0.00032
H& && && && && && &&&-4.67491& &2.67323&&-0.00142
C& && && && && && &&&-2.78312& &1.53566&&-0.00295
N& && && && && && &&&-1.97118& &2.59902&&-0.00677
H& && && && && && &&&-0.94447& &2.47785&&-0.00556
H& && && && && && &&&-2.35612& &3.5283& &-0.00593
N& && && && && && &&&-2.22645& &0.32201&&-0.0019
C& && && && && && &&&-2.98083&&-0.80415& &0.00232
O& && && && && && &&&-2.50721&&-1.94317& &0.00383
C& && && && && && &&&-5.19402&&-1.89151& &0.0094
H& && && && && && &&&-6.249& & -1.61876& &0.01151
H& && && && && && &&&-4.96461&&-2.48761& &0.89316
H& && && && && && &&&-4.96895&&-2.49088&&-0.87326
N& && && && && && &&&-4.39119&&-0.6678& & 0.00517
C& && && && && && && &4.956& &&&1.02489& &0.0076
H& && && && && && && &5.98275& &1.36194& &0.01099
N& && && && && && && &3.90339& &1.79511& &0.00616
C& && && && && && && &2.83466& &0.91972& &0.00123
C& && && && && && && &1.42549& &1.15149&&-0.00212
O& && && && && && && &0.81175& &2.22371&&-0.00113
N& && && && && && && &0.70705&&-0.06226&&-0.00686
H& && && && && && &&&-0.3204& & 0.03598&&-0.00732
C& && && && && && && &1.25481&&-1.32168&&-0.00857
N& && && && && && && &0.39068&&-2.3598& &-0.01649
H& && && && && && &&&-0.62272&&-2.23661&&-0.00774
H& && && && && && && &0.78727&&-3.28353&&-0.00748
N& && && && && && && &2.55822&&-1.54352&&-0.00484
C& && && && && && && &3.28063&&-0.4051& &-0.00032
C& && && && && && && &5.57116&&-1.45208& &0.00665
H& && && && && && && &5.40583&&-2.08274&&-0.86878
H& && && && && && && &5.43402&&-2.05767& &0.90467
H& && && && && && && &6.59106&&-1.06805&&-0.01497
N& && && && && && && &4.64986&&-0.32869& &0.00343
结果如下：
&&&&&&&&&& Convergence criterion not met.
SCF Done:&&E(UB+HF-LYP) =&&-2& &&&A.U. after&&601 cycles
& && && && & Convg&&=& & 0.4271D-03& && && && & -V/T =&&2.0045
& && && && & S**2& &=& &0.7578
KE= 1.D+03 PE=-5.D+03 EE= 1.D+03
Annihilation of the first spin contaminant:
S**2 before annihilation& &&&0.7578,& &after& &&&0.7500
Convergence failure -- run terminated.
Error termination via Lnk1e in /opt/APP//g03/l502.exe at Sun Feb&&8 09:45:48 2015.
Job cpu time:&&6 days&&9 hours 34 minutes 57.8 seconds.
File lengths (MBytes):&&RWF=& & 512 Int=& && &0 D2E=& && &0 Chk=& && &8 Scr=& && &8
请各位帮忙！！！:cry:
1. 已经尝试过改变初始猜测为core，或直接去掉guess指令，都仍然收敛失败；
2. 试过单纯加到maxcycle=800，收敛失败；
3. 试过加上novarint指令，收敛失败；
4. 论坛中有帖子提到降低收敛的条件，不要用tight那么苛刻的也行，试了一下收敛到5确实可以，但问过老师说在这个工作里不能精确度这么低，还是需要tight所以这条路是行不通的……
5. 有试过先跑小基组然后guess=read，不知道是不是我写的有问题，还是算不过去TAT，输入文件和结果如下：
输入文件：
%chk=MGMC3.chk
#p ub3lyp/6-311+g(d,p) pop=reg scf=(maxcycle=600,tight) guess=mix
Title Card Required
&&C& & & && && &&&-4.959961& & 0.567481& & 0.003638
&&H& & & && && &&&-6.042758& & 0.593683& & 0.005877
&&C& & & && && &&&-4.210823& & 1.697286& &-0.000320
&&H& & & && && &&&-4.674910& & 2.673230& &-0.001416
&&C& & & && && &&&-2.783118& & 1.535657& &-0.002947
&&N& & & && && &&&-1.971178& & 2.599023& &-0.006769
&&H& & & && && &&&-0.944468& & 2.477845& &-0.005562
&&H& & & && && &&&-2.356124& & 3.528297& &-0.005932
&&N& & & && && &&&-2.226450& & 0.322013& &-0.001902
&&C& & & && && &&&-2.980831& &-0.804154& & 0.002323
&&O& & & && && &&&-2.507206& &-1.943169& & 0.003829
&&C& & & && && &&&-5.194020& &-1.891506& & 0.009403
&&H& & & && && &&&-6.248998& &-1.618761& & 0.011510
&&H& & & && && &&&-4.964606& &-2.487608& & 0.893161
&&H& & & && && &&&-4.968948& &-2.490879& &-0.873261
&&N& & & && && &&&-4.391194& &-0.667797& & 0.005174
&&C& & & && && && &4.956001& & 1.024888& & 0.007602
&&H& & & && && && &5.982750& & 1.361941& & 0.010986
&&N& & & && && && &3.903393& & 1.795105& & 0.006164
&&C& & & && && && &2.834662& & 0.919721& & 0.001225
&&C& & & && && && &1.425485& & 1.151489& &-0.002118
&&O& & & && && && &0.811752& & 2.223710& &-0.001133
&&N& & & && && && &0.707050& &-0.062262& &-0.006862
&&H& & & && && &&&-0.320396& & 0.035975& &-0.007320
&&C& & & && && && &1.254809& &-1.321684& &-0.008569
&&N& & & && && && &0.390676& &-2.359796& &-0.016494
&&H& & & && && &&&-0.622717& &-2.236612& &-0.007742
&&H& & & && && && &0.787273& &-3.283532& &-0.007483
&&N& & & && && && &2.558222& &-1.543516& &-0.004840
&&C& & & && && && &3.280632& &-0.405102& &-0.000315
&&C& & & && && && &5.571163& &-1.452079& & 0.006646
&&H& & & && && && &5.405829& &-2.082737& &-0.868784
&&H& & & && && && &5.434015& &-2.057667& & 0.904665
&&H& & & && && && &6.591056& &-1.068051& &-0.014971
&&N& & & && && && &4.649857& &-0.328691& & 0.003429
--Link1--& && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && &
%nproc=8& && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && &
%mem=14GB& && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && &&&
%chk=MGMC3.chk& && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && &
#p ub3lyp/6-311++g(d,p) pop=reg scf=(maxcycle=500,tight) geom=AllCheck guess=Read& && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && &
&&&&&&&&&& Convergence criterion not met.
SCF Done:&&E(UB+HF-LYP) =&&-8& &&&A.U. after&&501 cycles
& && && && & Convg&&=& & 0.4118D-03& && && && & -V/T =&&2.0037
& && && && & S**2& &=& &0.7539
KE= 1.D+03 PE=-5.D+03 EE= 2.D+03
Annihilation of the first spin contaminant:
S**2 before annihilation& &&&0.7539,& &after& &&&0.7500
Convergence failure -- run terminated.
Error termination via Lnk1e in /opt/APP//g03/l502.exe at Mon Feb&&9 19:21:10 2015.
Job cpu time:&&4 days 22 hours 19 minutes 22.7 seconds.
File lengths (MBytes):&&RWF=& & 512 Int=& && &0 D2E=& && &0 Chk=& &&&28 Scr=& && &8
综上，我实在是已经脑力枯竭，无力回天，再拜！万请各位不吝赐教！:hand:
楼主用的方法是计算电子亲和势的一种。
不妨试试OVGF方法在优化结构基础上直接计算电子亲和势。此时电荷和自旋多重度仍然为0和1，应该比较容易收敛。
OVGF结果中最低空轨道对应的数值就是电子亲和势。
请注意，该数值不是自洽场迭代完之后给出的轨道能级，而是要找到如下的信息（以我之前的一个计算为例）：
& & Alpha HOMO is orbital&&12
这个12意思是第12个轨道，在你的分子中应该是一个挺大的数值。
那个我的这个计算电子亲和能就要找第13个轨道的信息
Summary of results for alpha spin-orbital&&13&&OVGF:
Koopmans theorem:& && && && &-0.18900D+00 au& &-5.143 eV
Converged second order pole: -0.23296D+00 au& &-6.339 eV&&0.951 (PS)
Converged third&&order pole: -0.22027D+00 au& &-5.994 eV&&0.959 (PS)
Outer Valence Approximation: -0.22284D+00 au& &-6.064 eV&&0.957 (PS)
其中最下面一行的 Outer Valence Approximation对应的-6.064 eV就是电子亲和能
计算方法很简单
# OVGF/你的基组 : Originally posted by hairan at
楼主用的方法是计算电子亲和势的一种。
不妨试试OVGF方法在优化结构基础上直接计算电子亲和势。此时电荷和自旋多重度仍然为0和1，应该比较容易收敛。
OVGF结果中最低空轨道对应的数值就是电子亲和势。
请注意，该 ... 呃，首先万分感谢您赐教！我写了一下试了试，所以花了比较久，不好意思。但不知道是不是因为我有哪里写出了差错，计算还是有错误……如果您愿意的话请再帮忙看看TAT多谢
输出文件：
SCF Done:&&E(UHF) =&&-6& &&&A.U. after& &60 cycles
& && && && & Convg&&=& & 0.1287D-08& && && && & -V/T =&&2.0020
& && && && & S**2& &=& &0.7502
Annihilation of the first spin contaminant:
S**2 before annihilation& &&&0.7502,& &after& &&&0.7500
ExpMin= 3.60D-02 ExpMax= 8.59D+03 ExpMxC= 1.30D+03 IAcc=3 IRadAn=& && && &5 AccDes= 0.00D+00
HarFok:&&IExCor= 205 AccDes= 0.00D+00 IRadAn=& && && &5 IDoV=1
ScaDFX=&&1....000000
Range of M.O.s used for correlation:& & 22& &560
NBasis=& &560 NAE=& & 77 NBE=& & 76 NFC=& & 21 NFV=& &&&0
NROrb=& & 539 NOA=& & 56 NOB=& & 55 NVA=& &483 NVB=& &484
**** Warning!!: The largest alpha MO coefficient is&&0.
**** Warning!!: The smallest alpha delta epsilon is&&0.
**** Warning!!: The largest beta MO coefficient is&&0.
Semi-Direct transformation.
ModeAB=& && && &&&4 MOrb=& && && && &56 LenV=& &
LASXX=& & LTotXX=&& LenRXX=&&
LTotAB=&& MaxLAS=&& LenRXY=& && && &&&0
LnScr1=&& LExtra=& && && &&&0 Total=&&
MaxDsk=& && && & -1 SrtSym=& && && &&&T ITran=& && && && &4
File extend in NtrExt1 failed.
You are probably out of disk space.
File extend in NtrExt1 probably out of disk space.
输入文件：
%chk=MGMC5.chk
# OVGF/6-311++g(d,p) pop=reg scf=(maxcycle=500,tight) guess=mix
Title Card Required
C& & & && && &&&-4.959961& & 0.567481& & 0.003638&&
H& & & && && &&&-6.042758& & 0.593683& & 0.005877&&
C& & & && && &&&-4.210823& & 1.697286& &-0.000320&&
H& & & && && &&&-4.674910& & 2.673230& &-0.001416&&
C& & & && && &&&-2.783118& & 1.535657& &-0.002947&&
N& & & && && &&&-1.971178& & 2.599023& &-0.006769&&
H& & & && && &&&-0.944468& & 2.477845& &-0.005562&&
H& & & && && &&&-2.356124& & 3.528297& &-0.005932&&
N& & & && && &&&-2.226450& & 0.322013& &-0.001902&&
C& & & && && &&&-2.980831& &-0.804154& & 0.002323&&
O& & & && && &&&-2.507206& &-1.943169& & 0.003829&&
C& & & && && &&&-5.194020& &-1.891506& & 0.009403&&
H& & & && && &&&-6.248998& &-1.618761& & 0.011510&&
H& & & && && &&&-4.964606& &-2.487608& & 0.893161&&
H& & & && && &&&-4.968948& &-2.490879& &-0.873261&&
N& & & && && &&&-4.391194& &-0.667797& & 0.005174&&
C& & & && && && &4.956001& & 1.024888& & 0.007602&&
H& & & && && && &5.982750& & 1.361941& & 0.010986&&
N& & & && && && &3.903393& & 1.795105& & 0.006164&&
C& & & && && && &2.834662& & 0.919721& & 0.001225&&
C& & & && && && &1.425485& & 1.151489& &-0.002118&&
O& & & && && && &0.811752& & 2.223710& &-0.001133&&
N& & & && && && &0.707050& &-0.062262& &-0.006862&&
H& & & && && &&&-0.320396& & 0.035975& &-0.007320&&
C& & & && && && &1.254809& &-1.321684& &-0.008569&&
N& & & && && && &0.390676& &-2.359796& &-0.016494&&
H& & & && && &&&-0.622717& &-2.236612& &-0.007742&&
H& & & && && && &0.787273& &-3.283532& &-0.007483&&
N& & & && && && &2.558222& &-1.543516& &-0.004840&&
C& & & && && && &3.280632& &-0.405102& &-0.000315&&
C& & & && && && &5.571163& &-1.452079& & 0.006646&&
H& & & && && && &5.405829& &-2.082737& &-0.868784&&
H& & & && && && &5.434015& &-2.057667& & 0.904665&&
H& & & && && && &6.591056& &-1.068051& &-0.014971&&
N& & & && && && &4.649857& &-0.328691& & 0.003429&&
再次感谢啊！ You are probably out of disk space.
这句话是主要的出错信息，意思是磁盘需求超过你的磁盘容量。
如果是在Windows下计算的，请用64位Linux来算，因为windows大多数是32位的系统，有单个文件2G，总读写文件16G的限制。
如果已经在64位linux下面做的计算，请换成一台硬盘大的机器来算。 : Originally posted by hairan at
You are probably out of disk space.
这句话是主要的出错信息，意思是磁盘需求超过你的磁盘容量。
如果是在Windows下计算的，请用64位Linux来算，因为windows大多数是32位的系统，有单个文件2G，总读写文件16G的 ... 这位不好意思哈，一直在外头这么久才爬回来回复。我看到out of disk space也觉得是计算受限制了，但本人只是个在实验室帮师兄师姐做给定课题的小小砖工，平时基本就是长个经验那个地步，所以没啥能力和资源换计算平台和硬盘什么的TAT就是老师给了啥就是啥了……这种情况下请问还能有别的可能的方法吗？我查了相关文献发现目前算过同样的东西的文章里基组要不就是6-31++G**，要不就是6-311+G**……可能我用的这个基组太大了？但因为其他数据都是这个基组，一旦换掉其他所有的都得从头再算，未免太烧了…………所以如果能有别的方法，真是跪求啊TAT再次感谢之前提供的建议！ 这么久了发现还是没人理我TAT只好厚着脸皮再踢一脚，希望有大神能帮帮忙:cry:
以及重新看了一下楼上那位高人的指点，发现有一个地方似乎和我的需要有点出入：我是在算单点的时候无法收敛，所以这个构型没优化过也不打算优化，因为想算一个瞬时的结构的性质，但好像提供的OVGF方法是优化好的构型？如果理解错了请指正，再次感谢，虽然不知道您能不能看到……牡丹江买的二手房，68平16万，满五年为首房，求大神帮忙算一下详细的过户费！_百度知道
牡丹江买的二手房，68平16万，满五年为首房，求大神帮忙算一下详细的过户费！
按房产证已满五年算，过户税费约为3000元。过户税费具体如下：二手房房产证过五年不用交个人所得税及营业税，只需交契税1.5%，其它过户税费约数百。即总共1.8%左右。报税价是你们自己定的。房管局的电脑系统对当地每个区域有一个最低的评估价。你们报的报税价高于那评估价就按你们的价计税。低于的话就按系统的评估价计算。对这评估价最了解的是当地同区域的房产中介，因为他们经常去过户，知道能最低报多少可以通过，所以建议你去中介那了解下。注：以上只为房产过户的房管局及税局收取的费用，未含房产中介费、银行贷款费及小区房的住房维修基金等在内。
来自团队：
其他类似问题
为您推荐：
过户费的相关知识
等待您来回答
下载知道APP
随时随地咨询
出门在外也不愁&& 查看话题
gibbs程序算热力学性质输出数据不合理是怎么回事，求大神帮我看看
在用gibbs的时候，输入文件如下，输出文件中没有体弹模量，Gmin(&&P=0)& & =& &&&498.60 kJ/mol&&为正值，而且德拜温度不合理，是什么原因呢，请高手帮我找找原因。，金币不多，只有这些了
输入文件如下：
sr2tio4 rutile phase. aiPI-uCHF calculation.
sample.out
11 0.0 2.0 4.0 6.0 8.0 10.0 12.0 14.0 16.0 18.0 20.0
6 0.0 100.0 200.0 300.0 400.0 500.0
& && &712.2588795& & & & -
& && &705.5076105& & & & -
& && &698.7563415& & & & -
& && &692.0050725& & & & -
& && &685.2538035& & & & -
& && &678.5025345& & & & -
& && &671.7512655& & & & -
& && &664.9999965& & & & -
& && &658.2487275& & & & -
& && &644.7461895& & & & -
& && &637.9949205& & & & -
& && &631.2436515& & & & -
& && &624.4923825& & & & -
& && &617.7411135& & & & -
& && &610.9898445& & & & -
部分输出文献如下：
sr2tio4 rutile phase. aiPI-uCHF calculation.
Number of data points:& && && & 15
Static EOS calculation - Numerical results
Vmin(&&P=0)& & =& & 657.48 bohr^3
Gmin(&&P=0)& & =& &&&498.60 kJ/mol
NUMERICAL EQUILIBRIUM PROPERTIES
================================
& &P(GPa) G(kJ/mol) V(bohr3)& &V/V0& &B(GPa) rel.err.
-----------------------------------------------------
& &&&0.00& & 498.60& &657.48 1.00000 ******* 0.237630
& &&&2.00& & 615.94& &657.46 0.99998 ******* 0.237199
& &&&4.00& & 733.28& &657.45 0.99996 ******* 0.236675
& &&&6.00& & 850.62& &657.43 0.99993 ******* 0.236171
& &&&8.00& & 967.95& &657.42 0.99991 ******* 0.235842
& & 10.00& &1085.29& &657.41 0.99989 ******* 0.235497
& & 12.00& &1202.62& &657.39 0.99987 ******* 0.234708
& & 14.00& &1319.94& &657.38 0.99985 ******* 0.234250
& & 16.00& &1437.27& &657.36 0.99982 ******* 0.234085
& & 18.00& &1554.59& &657.35 0.99980 ******* 0.233554
& & 20.00& &1671.91& &657.33 0.99978 ******* 0.232847
& & Debye temperature - numerical derivatives
& & Poisson coefficient: 0.250000&&Poisson ratio function: 0.859949
& && && && && & 610.99& & 14533.69& &&&17966.53
& && && && && & 617.74& & 14533.69& &&&17500.91
& && && && && & 624.49& & 14533.69& &&&17028.44
& && && && && & 631.24& & 14533.69& &&&16548.37
& && && && && & 637.99& & 14533.69& &&&16059.88
& && && && && & 644.75& & 14533.69& &&&15562.02
& && && && && & 658.25& & 14533.69& &&&14533.69
& && && && && & 665.00& & 14533.69& &&&14000.51
& && && && && & 698.76& & 14533.69& &&&11063.04
& && && && && & 705.51& & 14533.69& &&&10398.51
不好意思 借帖子请问gibbs程序在哪个网页可以下载到？ : Originally posted by gump_813 at
不好意思 借帖子请问gibbs程序在哪个网页可以下载到？ http://cpc.cs.qub.ac.uk/summaries/ADSY
CPC程序库，开源的 : Originally posted by happy小海 at
http://cpc.cs.qub.ac.uk/summaries/ADSY
CPC程序库，开源的... 谢了 你们是用这个做相变？
我画了一下你的V-E曲线 我觉得出现的问题可能是因为V-E曲线有个跳跃 导致优化出了问题 : Originally posted by gump_813 at
谢了 你们是用这个做相变？
我画了一下你的V-E曲线 我觉得出现的问题可能是因为V-E曲线有个跳跃 导致优化出了问题... v变化范围大的曲线就判断了一下一下机械稳定性，也算是相变吧，
这个v变化范围小的E-V打算用来算热学性质，主要是想算热膨胀
曲线我按你指出的问题又算了一下，确实是跳跃的地方有错误，数据多弄乱了，但是改正后还是这样的。。对输出结果影响不大，我再看看 用最新版本试试。这个软件出了gibbs2. : Originally posted by lz830617 at
用最新版本试试。这个软件出了gibbs2. 谢谢帮助，我去看了一下，gibbs2要用Octave，安上以后看着很晕呢，没怎么学过编程软件，matlab都不熟，貌似以我的水平很难搞定。。:cry: : Originally posted by happy小海 at
谢谢帮助，我去看了一下，gibbs2要用Octave，安上以后看着很晕呢，没怎么学过编程软件，matlab都不熟，貌似以我的水平很难搞定。。:cry:... gibbs2不需要matlab，fortran语言编写的 我还是不想放弃gibbs1，这个包输入输出命令是师兄帮我改好的，自己不会，我还想问一下B算不出来可能是由于E-V值范围或者密度取的不恰当吗，或者是我的体系比较大，单包7个原子，对数据选取更敏感？看了一下介绍gibbs的那篇文章有些这些内容。还有我有通过第一性原理弹性常数得来的最稳结构的B，gibbs可以直接用我自己的体弹模量值去算后面的热学性质吗，这样不就省掉了对E-V曲线敏感的近似求值过程，而且结果更准确？ : Originally posted by lz830617 at
gibbs2不需要matlab，fortran语言编写的
... 数据怎么修改都不出结果，表示要放弃gibbs1，请问你是用fortran软件运行的gibbs2是吗？因为gibbs2库里写的是用Octave，在fortran上运行也可以是吗 虽然本人很不推荐用Gibbs程序做热力学计算，事实上做CALPHAD的人没有一个推荐用Debye法的～不过抛开这个不谈，LZ问问题的方式本身就值得商榷。
我们知道，软件是最忠实的，你给它什么命令，它就干什么活。换句话说，LZ给我们一套输入，然后告诉我们输出很烂，让我们找原因，这个是很难的。因为LZ给的输入已经不是个“第一性”的东西了，我们没有人有能力去判断LZ的输入文件的数据是合理数据，还是garbage。今天老板还在说，如果input是garbage，那么output得到的不是garbage的几率几乎就没有。
所以，LZ如果真想figure out为什么Gibbs得到的结果有问题，就该先show出E-V曲线来。另外，对于LZ计算的东西，结构是rutile的，也就是tetragonal结构的。为什么连个poisson ratio都不给出呢？
顺便，LZ既然有Gibbs1的下载地址，就应该知道那个软件包里还有个gibbs.doc，不知道LZ有没有仔细看过，如果看过的话，不至于这个输入文件都要找个师兄给调试。人家又不用负责你的input是否合理，只要能运转，就可以了啊，人家也没必要为你的结果准确性负责不是么？ 可以的。linux下编译，算起来超快。 : Originally posted by souledge at
虽然本人很不推荐用Gibbs程序做热力学计算，事实上做CALPHAD的人没有一个推荐用Debye法的～不过抛开这个不谈，LZ问问题的方式本身就值得商榷。
我们知道，软件是最忠实的，你给它什么命令，它就干什么活。换句话说 ... 师兄给我调试的是fortran主程序中的文件读入和文件输出，在里面加了两个语句，如果不改就只有一个小黑框啊，我这编程能力让高手见笑了。输入文件有固定格式，很简单这个没问题，我算的是四方结构，由弹性常数算得的泊松比就是0.251，且gibbbs文章中说过，默认值0.25可以得到合理值。师兄调试的当然只管运转。。。我输出不合理当然是因为输入不合理。。。经验丰富的人知道从哪里找问题，是不是我的数据取值范围、密度不合理或者其他我不知道的原因，这是我的意图，亲。 请问楼主，你输入文件的第五行-1433.8是什么意思呢？具体如何获得呢？可能这里有问题！ : Originally posted by Brilliantly at
请问楼主，你输入文件的第五行-1433.8是什么意思呢？具体如何获得呢？可能这里有问题！ 问题已经解决，是我输入数据的单位跟gibbs要求的不一样，转换一下就可以了，可能是你说的那个数据，过去太久了记不清了。谢谢关注！ 楼主不容易，值得赞赏不是，顺便问问，Gibbs2知道如何使用了吗？如果知道，请分享分享。。。急！！这个考勤表要怎么做出一个月的迟到次数统计，总共几次。求高手解答！！9：30以后算迟到，计算的时间要10：00以内9：30以后的时间总次数，大神求帮忙！~_百度作业帮
急！！这个考勤表要怎么做出一个月的迟到次数统计，总共几次。求高手解答！！9：30以后算迟到，计算的时间要10：00以内9：30以后的时间总次数，大神求帮忙！~
急！！这个考勤表要怎么做出一个月的迟到次数统计，总共几次。求高手解答！！9：30以后算迟到，计算的时间要10：00以内9：30以后的时间总次数，大神求帮忙！~
问题：1，如何知道已经晚了？早退？ 2，什么是天数休息了一个月的基础上计算出来的？
和我们的比较类似，使用则科考勤系统，建议还是上个专业的收费系统，我说的也很值得考察