钻井水眼面积计算_百度文库
两大类热门资源免费畅读
续费一年阅读会员，立省24元！
评价文档：
&&¥80.00
喜欢此文档的还喜欢
钻井水眼面积计算
阅读已结束，如果下载本文需要使用
想免费下载本文？
把文档贴到Blog、BBS或个人站等：
普通尺寸(450*500pix)
较大尺寸(630*500pix)
你可能喜欢上传时间：
56官方微信
扫一扫发现精彩北航数值分析计算实习大作业2_百度文库
两大类热门资源免费畅读
续费一年阅读会员，立省24元！
评价文档：
喜欢此文档的还喜欢
北航数值分析计算实习大作业2
北​航​数​值​分​析​计​算​实​习​大​作​业​2
阅读已结束，如果下载本文需要使用
想免费下载本文？
把文档贴到Blog、BBS或个人站等：
普通尺寸(450*500pix)
较大尺寸(630*500pix)
你可能喜欢&& 查看话题
Openmx中Graphene的优化问题！
小虾最近在学习Openmx，想用它来优化Graphene，元胞中使用了两个C原子，但是一直得不到收敛的结果，程序一直运行到最大循环数然后被强制终止，请大侠帮忙看看问题出在哪了？谢谢！下面是我的输入文件：
# Graphene written by XXX
System.CurrrentDirectory /home/wangdb/Workspaces/Openmx/work-wdb/ # default=./
System.Name graphene
DATA.PATH /home/wangdb/Workspaces/Openmx/DFT_DATA06 # default=../DFT_DATA11
level.of.stdout 1 # default=1 (1-3)
level.of.fileout 1 # default=1 (1-3)
# Definition of Atomic Species
Species.Number 1
&Definition.of.Atomic.Species
C C6.0-s2p2d1 C_PBE
Definition.of.Atomic.Species&
Atoms.Number 2
Atoms.SpeciesAndCoordinates.Unit FRAC # Ang|AU
&Atoms.SpeciesAndCoordinates # Unit=Ang.
1 C 0. 0. 0. 2.0 2.0
2 C 0. -0. 0. 2.0 2.0
Atoms.SpeciesAndCoordinates&
Atoms.UnitVectors.Unit AU # Ang|AU
&Atoms.UnitVectors # unit=Ang.
2. 3.38 0.
2. -3.38 0.
0. 0. 20.6
Atoms.UnitVectors&
# SCF or Electronic System
scf.XcType LDA # LDA|LSDA-CA|LSDA-PW
scf.SpinPolarization off # On|Off
scf.ElectronicTemperature 300.0 # default=300 (K)
scf.energycutoff 700.0 # default=150 (Ry)
scf.maxIter 100 # default=40
scf.EigenvalueSolver band # Recursion|Cluster|Band
scf.lapack.dste dstevx # dstegr|dstedc|dstevx, default=dstegr
scf.Kgrid 24 24 1 # means 4x4x4
scf.Mixing.Type rmm-diisk # Simple|Rmm-Diis|Gr-Pulay
scf.Init.Mixing.Weight 0.010 # default=0.30
scf.Min.Mixing.Weight 0.001 # default=0.001
scf.Max.Mixing.Weight 0.200 # default=0.40
scf.Mixing.History 15 # default=5
scf.Mixing.StartPulay 5 # default=6
scf.criterion 1.0e-6 # default=1.0e-6 (Hartree)
1DFFT.NumGridK 900 # default=900
1DFFT.NumGridR 900 # default=900
1DFFT.EnergyCutoff 3600.0 # default=3600 (Ry)
# SCF Order-N
orderN.HoppingRanges 8.389 # default=5.0 (Ang)
orderN.NumHoppings 2 # default=2
orderN.KrylovH.order 800 # default=400
# MD or Geometry Optimization
MD.Type opt # Nomd|Opt|DIIS|NVE|NVT_VS|NVT_NH
MD.Opt.DIIS.History 7 # default=7
MD.Opt.StartDIIS 5 # default=5
MD.maxIter 100 # default=1
MD.TimeStep 1.0 # default=0.5 (fs)
MD.Opt.criterion 1.0e-7 # default=1.0e-4 (Hartree/bohr)
我覺得你直接去 OPENMX 官方發問會比較快
http://www.openmx-square.org/forum/patio.cgi
小木蟲上似乎很很很少人用 OPENMX
OPENMX.... 不好上手 scf.energycutoff 200.0 # default=150 (Ry)
&Definition.of.Atomic.Species
C C6.0-s2p2d1 C_CA13
Definition.of.Atomic.Species&
改动2处，可以收敛
*******************************************************
& && && && && & Total Energy (Hartree) at MD = 4& && &&&
*******************************************************
&&Uele&&=& && &-3.
&&Ukin&&=& && & 8.
&&UH0& &=& & -115.
&&UH1& &=& && & 0.
&&Una& &=& && &-5.
&&Unl& &=& && &-1.
&&Uxc0&&=& && &-1.
&&Uxc1&&=& && &-1.
&&Ucore =& &&&106.
&&Uhub&&=& && & 0.
&&Ucs& &=& && & 0.
&&Uzs& &=& && & 0.
&&Uzo& &=& && & 0.
&&Uef& &=& && & 0.
&&UvdW&&=& && & 0.
&&Utot&&=& &&&-11.
&&Utot = Ukin+UH0+UH1+Una+Unl+Uxc0+Uxc1+Ucore+Uhub+Ucs+Uzs+Uzo+Uef+UvdW
&&Uele:& &band energy
&&Ukin:& &kinetic energy
&&UH0:& & electric part of screened Coulomb energy
&&UH1:& & difference electron-electron Coulomb energy
&&Una:& & neutral atom potential energy
&&Unl:& & non-local potential energy
&&Uxc0:& &exchange-correlation energy for alpha spin
&&Uxc1:& &exchange-correlation energy for beta spin
&&Ucore:&&core-core Coulomb energy
&&Uhub:& &LDA+U energy
&&Ucs:& & constraint energy for spin orientation
&&Uzs:& & Zeeman term for spin magnetic moment
&&Uzo:& & Zeeman term for orbital magnetic moment
&&Uef:& & electric energy by electric field
&&UvdW:& &semi-empirical vdW energy
&&(see also PRB 72, 5) for the energy contributions)
*******************************************************
& && && &&&Computational times (s) at MD = 4& && && && &
*******************************************************
&&DFT in total& && &=& &12.12758
&&Set_OLP_Kin& && & =& & 1.44986
&&Set_Nonlocal& && &=& & 2.16248
&&Set_ProExpn_VNA& &=& & 2.26113
&&Set_Hamiltonian& &=& & 1.70408
&&Poisson& && && &&&=& & 0.00246
&&diagonalization& &=& & 0.50167
&&Mixing_DM& && && &=& & 0.01076
&&Force& && && && & =& & 1.85790
&&Total_Energy& && &=& & 0.92433
&&Set_Aden_Grid& &&&=& & 0.02858
&&Set_Orbitals_Grid =& & 0.27365
&&Set_Density_Grid&&=& & 0.89989
&&RestartFileDFT& & =& & 0.02197
&&Mulliken_Charge& &=& & 0.00213
&&FFT(2D)_Density& &=& & 0.00316
*******************************************************
& && && && & MD or geometry opt. at MD = 4& && && && &&&
*******************************************************
&Steepest_Descent&&&SD_scaling= 0.
&Steepest_Descent&&&|Maximum force| (Hartree/Bohr) = 0.
&Steepest_Descent&&&Criterion& && & (Hartree/Bohr) = 0.
&&atom=&&1, XYZ(ang) Fxyz(a.u.)=&&-0.0000& &-0.0000& & 0.0000& & -0.0000& & 0.0000& & 0.0000
&&atom=&&2, XYZ(ang) Fxyz(a.u.)=& &1.2195& &-0.7041& &-0.0000& &&&0.0000& &-0.0000& &-0.0000
outputting data on grids to files...
***********************************************************
***********************************************************
& && && && && &Computational Time (second)& && && && && &&&
***********************************************************
***********************************************************
& && && && && && && && && && & Min_ID& &Min_Time& && & Max_ID& &Max_Time
& &Total Computational Time =& &&&4& && & 67.461& && && & 0& && & 68.100
& &readfile& && && && && &&&=& &&&7& && &&&6.883& && && & 0& && &&&6.886
& &truncation& && && && && &=& &&&5& && &&&1.936& && && & 1& && &&&2.428
& &MD_pac& && && && && && & =& &&&0& && &&&0.028& && && & 5& && &&&0.036
& &OutData& && && && && && &=& &&&1& && &&&1.956& && && & 0& && &&&2.593
& &DFT& && && && && && && & =& &&&1& && & 56.048& && && & 5& && & 56.540
*** In DFT ***
& &Set_OLP_Kin& && && && &&&=& &&&5& && &&&1.045& && && & 0& && &&&5.790
& &Set_Nonlocal& && && && & =& &&&0& && &&&8.654& && && & 5& && & 13.400
& &Set_ProExpn_VNA& && && & =& &&&5& && &&&6.023& && && & 0& && &&&8.983
& &Set_Hamiltonian& && && & =& &&&1& && & 11.370& && && & 6& && & 11.370
& &Poisson& && && && && && &=& &&&7& && &&&0.014& && && & 1& && &&&0.015
& &Diagonalization& && && & =& &&&0& && &&&3.395& && && & 1& && &&&3.395
& &Mixing_DM& && && && && & =& &&&7& && &&&0.070& && && & 0& && &&&0.070
& &Force& && && && && && &&&=& &&&1& && &&&7.363& && && & 6& && &&&7.363
& &Total_Energy& && && && & =& &&&6& && &&&3.596& && && & 0& && &&&3.693
& &Set_Aden_Grid& && && && &=& &&&0& && &&&0.116& && && & 5& && &&&3.067
& &Set_Orbitals_Grid& && &&&=& &&&3& && &&&0.000& && && & 0& && &&&1.102
& &Set_Density_Grid& && && &=& &&&4& && &&&5.358& && && & 0& && &&&5.364
& &RestartFileDFT& && && &&&=& &&&5& && &&&0.066& && && & 7& && &&&0.895
& &Mulliken_Charge& && && & =& &&&1& && &&&0.009& && && & 7& && &&&0.010
& &FFT(2D)_Density& && && & =& &&&7& && &&&0.068& && && & 3& && &&&0.068
& &Others& && && && && && & =& &&&1& && &&&0.091& && && & 5& && &&&1.694
The calculation was normally finished. 谢谢y1ding，问题已解决，不过我发现在3.5上仍然不收敛，在3.7上运行到MD=13时计算正常结束了。以后还请多多指教！
var cpro_id = 'u1216994';
欢迎监督和反馈：本帖内容由
提供，小木虫仅提供交流平台，不对该内容负责。欢迎协助我们监督管理，共同维护互联网健康，如果您对该内容有异议，请立即发邮件到
联系通知管理员，也可以通过QQ周知，我们的QQ号为：8835100
我们保证在1个工作日内给予处理和答复，谢谢您的监督。
小木虫，学术科研第一站，为中国学术科研研究提供免费动力
广告投放请联系QQ： &
违规贴举报删除请联系邮箱： 或者 QQ:8835100
Copyright &
eMuch.net, All Rights Reserved. 小木虫 版权所有.0.0.0.-保定找情人-帅哥找情人-情人网&&&&下一个我感兴趣的情人：你现在的位置：&&>>&&&&>>&& 昵称：诚信值:0 普通会员性别：男身高：179年龄：19出生年月：199508婚姻状况：未婚最高学历：本科从事职业：个体老板年收入：_50_100万元现居住地：河北&&&&蠡县寻找：异性情人个性标签：他的联系方式您需要 / 会员，才能查看！他的交友密语我要寻找：寂寞 我能提供：不提供帮助 我希望得到：不提供帮助 我的爱情观念：您需要 / 会员，才能查看！ 我的外貌：英俊 我的个性：成熟 兴趣爱好：拥抱 自我介绍：我是一位来自河北保定的帅哥，与大家一样的目的相约玫瑰情人交友网 期望的情人：在玫瑰情人网寻觅一个漂亮温柔的女情人 他收到的礼物我和他之间的对话写信回信&&&&&& 请选择一句话推荐实力会员john2012leisanshui成熟男士爱情的美凯迪拉无敌孤独人情人网提示举报投诉投诉对象：.0.0.0.&&&&&ID:2565409投诉原因：酒托代孕钱财欺诈昵称资料不雅联系方式虚假其它补充说明：　|　　|　　|　　|　　|　　|　　|　　|　客服电话：400-000-7969（周一至周六 9:00-18:00) &&&&&&官方客服QQ：&&&&&&&&&&&&网站备案号：增值电信业务经营许可证：| 版权所有：& 湖北爱爱网络有限公司&& 法律声明：禁止未成年人、色情服务者、婚姻幸福者注册，禁止在本站发布任何色情信息。本站所表述的“情人”是指“情侣”。