各种计算中INCAR设置汇总
--------------------------------------------
Band
--------------------------------------------
ICHARG = 11
ENCUT = 600
PREC = Accuracy precision:Medium/High/Low default medium . the same time, ROPT accuracy of approximately 1meV/atom (10e-3). recommend to use PREC Normal for everyday calculations in VASP.4.5
LWAVE = .FALSE.
LCHARG= .FALSE.
Band计算要读取静态（精度提高后比如提高PREC,ENCUT,）的CHGCAR，除了ICHARG=11以及KPOINTS改为线性模式以外，其他的设置与静态计算的设置一致。
--------------------------------------------
bader
--------------------------------------------
ISTART=1
ICHARG=0
LCHARG= .TURE.
LAECHG=.TURE.
NGX=132
NGY=188
NGZ=268
~
Bader计算读取优化后或静态（与优化一致）的WAVECAR（ISTART=1）,要保存Bader计算后的CHGCAR(LCHARG=.TRUE.),要设置LAECHG=.TRUE.，NGX/Y/Z的设置是优化计算中OUTCAR里面推荐的NGX/Y/Z值的2倍（比如，优化时OUTCAR中WARNING: aliasing errors must be expected set NGX to 54 to avoid them，此时bader计算中应设置NGX=108）.其他设置与优化时一致。
-------------------------------------------
GGA+U
-------------------------------------------
LMAXMIX=4! for d elements increase LMAXMIX to 4, f: LMAXMIX = 6
LDAU=.TURE.
LDAUTYPE=2
LDAUU=3 0 0
LDAUL=2 -1 -1
GGA+U设置，可从OUTCAR中看出，
LDA+U is selected, type is set to LDAUTYPE = 2
angular momentum for each species LDAUL = 2 -1 -1
U (eV) for each species LDAUU = 3.0 0.0 0.0
J (eV) for each species LDAUJ = 0.0 0.0 0.0
---------------------------------------------
HSE HF
----------------------------------------------
ISTART=1
ICHARG=1
GGA=PE RP stands for RPBE, PE for PBE, 91 for PW91
LHFCALC = .TRUE. Hartree Fock is set to；The flag specifies, whether Hartree-Fock type calculations are performed
HFSCREEN = 0.207screening length (either q_TF or 0.3 A-1)
ALGO = Damped ;
TIME = 0.3timestep for ELM
PRECFOCK= FastNormal, Fast or Accurate (Low or Medium for compatibility)
NKREDZ=3reduce k-point grid by
--------------------------------------------
Output Control:
--------------------------------------------
LWAVE=.FALSE.
LCHARG=.FALSE.
LORBIT=11
EMIN=-10.0
EMAX=10.0
NEDOS=2001
HSE计算时由于计算量较大，一般只做单点校正，顺便会一起做DOS计算。做HSE计算时由于只做单点优化，因此最好读取PBE或RPBE优化后的WAVECAR/CHGCAR.具体设置只需在PBE 优化中加入HF部分即可，HF的计算一般设置LHFCALC = .TRUE.,ALGO = Damped, TIME = 0.3(一般推荐设为0.4，但不易收敛时可将其调小),具体设置详见VASP5.3.5手册119-121，123页，也可参见网页https://www.wendangku.net/doc/c015420786.html,/blog-567091-732988.html
--------------------------------------------
DOS Related:
--------------------------------------------
ISTART =1
ICHARG =0
PREC = Accurate precision:Medium/High/Low default medium . the same time, ROPT accuracy of approximately
ISMEAR = 0
SIGMA = 0.02
LORBIT =11
NEDOS=2101
EMAX=6
EMIN=-15
DOS计算要读取静态或优化后的WAVECAR,然后提高精度ENCUT,PREC=Accurate,KPOINTS （abc）设置比优化的要高即a+2,b+2,c+2,加2是为了使取点为基数，也可以+4/6/8等，但计算量会很大，有时候服务器承受不了，计算不动，因此不宜取太大；或者读取提高精度后(ENCUT,PREC,KPOINTS)静态的CHARG,然后将SIGMA设小一点，之后就是能窗的设置，一般是优化后的E-Fermi值加减10就可以，比如上面所设的值，优化后的fermi能级应该在4-5之间，然后EMAX设为6，EMIN设为-15，此时NEDOS=(EMAX-EMIN)*100+1，加1是由于加上本身，比如0-100之间总共有100+1个个数。其他的设置与一般优化设置相同，最后不必保留计算后的WAVECAR,CHGCAR,我们需要的是DOSCAR
-------------------------------------------
PARCHG
-------------------------------------------
LPARD=.TURE.Evaluate partial (band and/or k-point) decomposed charge density.
NBMOD=-3 ! > 0 Number of values in the array IBAND . If IBAND is specified, NBMOD is set automatically to the correct value (in thatcase NBMOD should not be set manually in the INCAR file) 0 Take all bands to calculate the charge density, even unoccupied bands are taken into account.
-1 Calculate the total charge density as usual. This is the default value if nothing else is given.
-2 Calculate the partial charge density for electrons with there eigenvalues in the range specified by EINT .
-3 The same as before, but the energy range is given vs. the Fermi energy
EINT= -1.2
LSEPB=.FALSE.Specifies whether the charge density is calculated for every band separately and written to a file PARCHG.nb.?(TRUE) or whether charge density is merged for all selected bands and write to the file PARCHG.ALLB.? or PARCHG.Default is FALSE.
LSEPK=.FALSE.Specifies whether the charge density of every k-point is write to the files PARCHG.?.nk (TRUE) or whether itis merged (FALSE) to a single file. If the merged file is written, then the weight of each k-point is determined from theKPOINTS file, otherwise the kpoints weights of one are chosen.
PARCHG分析主要是可以分析DOS图中具体的轨道，从而看其组成。PARCHG计算主要是基于DOS计算之后选取某一能区（能量区间）范围内的轨道进行分析，相当于是从DOSCAR中取出单独的信息，屏蔽不需要的轨道信息，因此不需要另外单独的计算，只是分解DOSCAR而已，孤儿其他的设置应该与DOS计算的设置一致，有时候KPOINTS的设置也可以小点，主要是精度会上不去，但都是想要单独分析某一能区的信息。ENIT设置区间。详见VASP5.3.5手册105页
--------------------------------------------
TS opt
--------------------------------------------
POTIM =0.05 default is 0.5 for frequency calculations a smaller values should be set ,0.02 example.time-step for ionic-motion
EDIFFG =-0.1 Converge criterion(positive for energy,negetive for force)
NSW = 300 number of steps for IOM
IBRION = 2 -1:no update 0:MD 1:quasi-New 2:CG
IMAGES=4time-step for ionic-motion
SPRING=-5recommended setting; The nudged elastic band method is applied when SPRING is set to a negative value
LCLIMB=.FALSE.
NEB方法计算过渡态，在00,05（XX）中间插入4个点（IMAGES=4，XX=number of images+1）,通过SPRING设置力常数( In addition to the IMAGES tag, a spring constant can be supplied in the SPRING tag)，SPRING为负数时表示采用NEB方法计算力常数，一般默认为-5，LCLIMB=.FALSE. 表示采用NEB方法计算，如果LCLIMB=.TRUE.则表示采用CINEB方法计算过渡态，一般最好不要用CINEB，有时候结果会比较诡异，一般难以收敛，除非过渡态很简单；过渡态优化时只能用POTIM=0.05，无论哪种计算方法，QN或CG;刚开始计算TS时可取精度低一点比如力收敛标准为0.1，此时可以做单点（能量最高的点）优化，或者继续NEB优化，刚开始切不可取过大的NSW，一般300步以内力不到0.1左右，则需改方法或改变初始猜测。也可参见网页https://www.wendangku.net/doc/c015420786.html,/vtsttools/
--------------------------------------------
optmization for freq
--------------------------------------------
LWAVE = .TURE.
LCHARG= .TURE.
POTIM =0.05 default is 0.5 for frequency calculations a smaller values should be set ,0.02 example.
IBRION = 1 -1:no update 0:MD 1:quasi-New 2:CG
EDIFFG =-0.03 Converge criterion(positive for energy,negetive for force)
NSW = 1000 number of steps for IOM
NEB插入点做过渡态后取能量最高的点然后只能用QN的方法优化得到过渡态，此时POTIM可取范围0.05~0.2，最好刚开始用0.05，当力接近0.1但一直在0.1附近徘徊时，可以重新建POTIM 取大一点比如0.1或0.2，之后再继续优化，但当力趋近0.05又一直不能更小时，这时可以将能量的收敛精度调小比如EDIFF=1e-5，然后继续优化，一般过渡态收敛到力≤0.05，能量标准1e-5就可以啦（有时候过渡态比较难收敛时做到0.05就可以啦）。最好保留WAVECAR或CHARG，下一步频率计算需要用到。
--------------------------------------------
frequency
--------------------------------------------
ISTART =1
ICHARG =1
POTIM =0.02 default is 0.5 for frequency calculations a smaller values should be
set ,0.02 example.
NSW = 1 number of steps for IOM
IBRION = 5 -1:no update 0:MD 1:quasi-New 2:CG
# NPAR = 4
NFREE=2
频率计算时要读上一步优化的WAVECAR或CHGCAR,或2者同时读取也可以，有时候也可以不用读，但一般是磁性体系或是自选变化比较太的体系其态不易收敛的必须要读，主要是设置POTIM=0.02或0.01，NFREE=2或4，NSW=1，IBRION=5，这4个参数是计算频率的必设，其中NFREE=2表示是在我们所取点的左右2边各取一个点来拟合曲线做其二阶导求力常数k，继而求得频率，详见VASP5.3.5手册60-61页（NFREE =2 uses central difference, i.e. each ion is displaced in each direction by a small positive and negative displacement，±POTIM × ? x , ± POTIM × ? y, ± POTIM × ? z ；For NFREE =4, four displacement are used，±POTIM × ? x and ± 2 POTIM × ? x，±POTIM × ? y and ± 2 POTIM × ? y…，For NFREE =1, only a single displacement is applied (it is strongly recommend to avoid NFREE =1).）；For IBRION =5, all atoms are displaced in all three Cartesian directions, resulting in a significant computational effort even formoderately sized high symmetry systems. For IBRION =6, however only symmetry inequivalent displacements are considered,and the remainder of the Hessian matrix is filled using symmetry considerations.Selective dynamics are presently only supported for IBRION =5; The parameter NFREE determines how many displacements are used for each direction and ion, and POTIM determines the step size. The step size is defaulted to 0.015 ? A，我们推荐的使用0.01或0.02
其中对于我们的服务器cluster，不知道出于什么原因要注释掉NPAR,不然设为1或4都会报错，然后死掉。
-------------------------------------------
MD Ionic Relaxation:
-------------------------------------------
EDIFFG = -0.03 Converge criterion(positive for energy,negetive for force) POTIM = 2 default is 0.5 for frequency calculations a smaller values should be set ,0.02 example. time step in fs
NSW = 2000
IBRION = 0 -1:no update 0:MD 1:quasi-New 2:CG
SMASS = 3controls the velocities during an ab-initio molecular dynamics.
TEBEG = 300start temperature,temperature in K
TEEND = 300final temperature
NBLOCK = 1 !After NBLOCK ionic steps the pair correlation function and the DOS are calculated and the ionic configuration will be writtento the XDATCAR-file,The CPU costs for these tasks are quite small so use NBLOCK =1.After KBLOCK*NBLOCK main loops the averaged pair correlation function and DOS are written to the files PCDAT and
DOSCAR .
NWRITE = 0write down eigenvectors of dynamical matrix after division by SQRT(mass)
--------------------------------------------
MD计算我做的比较少，不是很熟悉,详见VASP5.3.5手册52，65-66页。TEBEG and TEEND
control the temperature during an ab-initio molecular dynamics run; This means that the real simulation temperature isT = TEBEG×N ions /(N ions ?1).TEBEG = T requested ×(N ions ?1)/N ions ,
Vasp中一些网站，能量含义的https://www.wendangku.net/doc/c015420786.html,/blog-779486-671812.html
上述只是我个人（pwang）的理解和一些经验，不足之处欢迎大家来更正，补充啊
锅炉耗水量计算
§2 锅炉基本特性的表示 为了区别各类锅炉构造、燃用燃料、燃烧方式、容量大小、参数高低以及运行经济性等特点，经常用到如下参数： 一、锅炉额定出力 锅炉额定出力是指锅炉在额定参数（压力、温度）和保证一定效率下的最大连续出力。对于蒸汽锅炉，叫额定蒸发量，单位为吨/小时；对于热水锅炉，叫额定产热量。单位为MW（老单位为万大卡/小时）。 产热量与蒸发量之间的关系： Q=D（iq-igs）×1000 千焦/小时 式中：D----锅炉蒸发量，吨/小时 iq----蒸汽焓，千焦/公斤 igs----锅炉给水焓，千焦/公斤 对于热水锅炉： Q=G（irs “-irs…）×1000 千焦/小时 式中：G----热水锅炉循环水量，吨/小时 irs “---锅炉出水焓，千焦/公斤 irs …---锅炉进水焓，千焦/公斤 注：1千卡（kcal）=4.1868千焦(KJ) 二、蒸汽（或热水）参数 锅炉产生蒸汽的参数，是指锅炉出口处蒸汽的额定压力（表压）和温度。对生产饱和蒸汽的锅炉来说，一般只标明蒸汽压力；对生产过热蒸汽的锅炉，则需标明压力和过热蒸汽温度；对热水锅炉来说，则需标明出水压力和温度。 工业锅炉的容量、参数，既要满足生产工艺上对蒸汽的要求，又要便于锅炉房的设计，
锅炉配套设备的供应以及锅炉本身的标准化，因而要求有一定的锅炉参数系列。见 GB1921-88《工业蒸汽锅炉参数系列》及GB3166-88《热水锅炉参数系列》GB1921-88《工业蒸汽锅炉参数系列》 额定蒸发量 t/h 额定出口蒸汽压力MPa （表压） 0.4 0.7 1.0 1.25 1.6 2.5 额定出口蒸汽温度℃ 饱和饱和饱和饱和250 350 饱和350 饱和350 400 0.1 ★ 0.2 ★ 0.5 ★★ 1 ★★★ 2 ★★★★ 4 ★★★★★ 6 ★★★★★★★ 8 ★★★★★★★ 10 ★★★★★★★★★ 15 ★★★★★★★★ 20 ★★★★★★★ 35 ★★★★★★ 65 ★★ 本表中的额定蒸发量,对于<6t/h的饱和蒸汽锅炉是20℃给水温度下锅炉额定蒸发量，对
VASP参数设置详解
VASP参数设置详解 计算材料2010-11-30 20:11:32 阅读197 评论0 字号：大中小订阅 转自小木虫，略有增减 软件主要功能： 采用周期性边界条件（或超原胞模型）处理原子、分子、团簇、纳米线（或管）、薄膜、晶体、准晶和无定性材料，以及表面体系和固体 l 计算材料的结构参数（键长、键角、晶格常数、原子位置等）和构型 l 计算材料的状态方程和力学性质（体弹性模量和弹性常数） l 计算材料的电子结构（能级、电荷密度分布、能带、电子态密度和ELF） l 计算材料的光学性质 l 计算材料的磁学性质 l 计算材料的晶格动力学性质（声子谱等） l 表面体系的模拟（重构、表面态和STM模拟） l 从头分子动力学模拟 l 计算材料的激发态（GW准粒子修正） 计算主要的四个参数文件：INCAR ,POSCAR,POTCAR ,KPOINTS,下面简要介绍，详细权威的请参照手册 INCAR文件： 该文件控制VASP进行何种性质的计算，并设置了计算方法中一些重要的参数，这些参数主要包括以下几类： 对所计算的体系进行注释：SYSTEM
●定义如何输入或构造初始的电荷密度和波函数：ISTART，ICHARG，INIWAV ●定义电子的优化 –平面波切断动能和缀加电荷时的切断值：ENCUT，ENAUG –电子部分优化的方法：ALGO，IALGO，LDIAG –电荷密度混合的方法：IMIX，AMIX，AMIN，BMIX，AMIX_MAG，BMIX_MAG，WC，INIMIX，MIXPRE，MAXMIX –自洽迭代步数和收敛标准：NELM，NELMIN，NELMDL，EDIFF ●定义离子或原子的优化 –原子位置优化的方法、移动的步长和步数：IBRION，NFREE，POTIM，NSW –分子动力学相关参数：SMASS，TEBEG，TEEND，POMASS，NBLOCK，KBLOCK，PSTRESS –离子弛豫收敛标准：EDIFFG ●定义态密度积分的方法和参数 –smearing方法和参数：ISMEAR，SIGMA –计算态密度时能量范围和点数：EMIN，EMAX，NEDOS –计算分波态密度的参数：RWIGS，LORBIT ●其它 –计算精度控制：PREC –磁性计算：ISPIN，MAGMOM，NUPDOWN –交换关联函数：GGA，VOSKOWN –计算ELF和总的局域势：LELF，LVTOT –结构优化参数：ISIF –等等。 主要参数说明如下： ?SYSTEM：该输入文件所要执行的任务的名字。取值：字符串，缺省值：SYSTEM ?NWRITE：输出内容详细程度。取值：0~4，缺省值：2
c15111期权中希腊字母100分
一、单项选择题 1. （）度量了期权价格对期权存续期的敏感性。 A. Delta B. Theta C. Vega D. Rho 您的答案：B 题目分数：10 此题得分：10.0 2. 其他条件不变的情况下，存续期越长，利率的影响越明显，期 权的Rho的绝对值（）。 A. 越大 B. 越小 C. 不变 D. 以上都不正确 您的答案：A 题目分数：10 此题得分：10.0 3. 当期权处于（）状态时，时间价值最大。 A. 实值 B. 虚值 C. 平值 D. 极端实值 您的答案：C 题目分数：10 此题得分：10.0 4. 当认沽期权处于极端实值状态，股票价格变动1元，内在价值 （）变动趋向1元。 A. 横向 B. 正向 C. 纵向 D. 反向 您的答案：D
题目分数：10 此题得分：10.0 二、多项选择题 5. 下列关于期权的Vega值的说法正确的是（）。 A. 随着到期日的临近，Vega值逐渐减小 B. 波动率越大，实值和虚值期权的Vega值越大，平值期 权的Vega值相对较平稳 C. 当期权处于平值时，Vega值最大 D. 当期权处于平值时，Vega值最小 您的答案：B,A,C 题目分数：10 此题得分：10.0 6. 下列关于期权的Theta值说法正确的是（）。 A. 当期权处于平值时，Theta值最小 B. 平值期权Theta逐渐减小，时间价值加速流失 C. 波动率越高，Theta值越小 D. 一般情况下期权的Theta值均为负 您的答案：A,C,D,B 题目分数：10 此题得分：10.0 三、判断题 7. 期权为平值时Gamma最大，Delta变化最快。（） 您的答案：正确 题目分数：10 此题得分：10.0 8. 当认购期权处于极端实值时，时间价值将变小，Delta趋向0。 （） 您的答案：错误 题目分数：10 此题得分：10.0 9. 波动率越高，时间价值越大，期权的Theta值越小，相同存续 期内时间价值流失越快。（）
趣谈期权有关的希腊字母
趣谈期权有关的希腊字母！Delta, Gamma, Vega和Theta 当我们理解期权价值与其影响因素的敏感性时，可以作这样比喻。股票期权作为股票的“孩子”，其脾气秉性自然受三方面的影响：一是自身“基因”的制约，比如：权利属性（认购还是认沽）、行权价（K）、到期时间（T）；二是“父母亲”的言传身教：股价（S）、股价的波动率（Sigma）；三是社会大环境的熏陶：无风险收益率（r）。 那么一份股票期权的价格（V）究竟是如何被这些因素所影响的呢？换而言之，股票价格上涨1％，或者股价波动率上升1%，作为孩子的期权的“脾气”变化多少呢？为了回答这个问题，我们就必须认识五个“希腊字母”了。毫不夸张地说，这五个希腊字母就是期权价格变化的生命源泉，也是“孩子”与“父母”的纽带。这五个希腊字母就叫做Delta，Gamma，Vega，Theta和Rho。 先让我们来认识第一个希腊字母—— Delta。 1． Delta是什么？ 期权是标的资产的衍生产品。两者之间就像是“父子”一样，父亲的一举一动无时无刻不在影响着孩子的行为。父亲的这种影响力就是Delta。 以50ETF为例，当ETF价格发生变化时，期权价格也会随之改变。ETF与期权之间的价格关系可以用Delta来刻画：当ETF价格变化0.001元时，对期权价格的影响就是0.001*Delta元。 认购期权是“乖孩子”，当“父亲”ETF价格上涨的时候，认购期权价格也会上涨，认购期权的Delta大于零；而“坏孩子”认沽期权则恰恰相反，当ETF 价格上涨时，认沽期权的价格反而是下跌的，它的Delta小于零。 2． Delta在投资中的两个简单应用 一个是对冲作用。如果我们有着如下对冲组合：由Delta份ETF空头和1份期权多头组成。当ETF价格变化0.001元时，Delta份ETF空头价格会变化-0.001*Delta元，1份期权合约价格会变化0.001*Delta元。两者相互抵消，对冲组合的整体价格几乎不变。因此，我们可以用Delta份ETF空头去对冲1份期权。 另一个是计算杠杆。我们知道期权具有一定的杠杆性。比如ETF上涨1%，期权上涨10%，那么期权的杠杆就是10倍。那么通过Delta，我们可以计算期权的杠杆倍数。假设目前50ETF的价格是3.00元，有一份1个月后到期行权价为3.20的认购期权，现在的价格是0.10元，Delta为0.33。如果ETF上涨1%，也就是0.030元，期权价格就会上涨0.030*Delta，等于0.1元。从涨幅来看，期权合约上涨了10%。因此，这个期权合约的杠杆大概是10倍。 （1）Delta与标的价格的变动关系 无巧不巧，不论是认购还是认沽期权，Delta的绝对值都介于0与1之间，而且越实值的期权Delta越接近于1，越虚值的期权Delta越趋近于0，平值期权的Delta恰好是0.5。因此我们也可以把Delta想象成期权到期实值的概率。
VASP控制参数文件INCAR的简单介绍
限于能力，只对部分最基本的一些参数(>,没有这个标志的参数都是可以不出现的) 详细说明，在这里只是简单介绍这些参数的设置，详细的问题在后文具体示例中展开。 部分可能会干扰VASP运行的参数在这里被刻意隐去了，需要的同学还是请查看VASP自带的帮助文档原文。 参数列表如下： >SYSTEM name of System 任务的名字*** >NWRITE verbosity write-flag (how much is written) 输出内容详细程度0-3 缺省2 如果是做长时间动力学计算的话最好选0或1(首末步/每步核运动输出) 据说也可以结合shell的tail或grep命令手动输出 >ISTART startjob: restart选项0-3 缺省0/1 for 无/有前次计算的WAVECAR(波函数) 1 'restart with constant energy cut-off' 2 'restart with constant basis set' 3 'full restart including wave function and charge prediction' ICHARG charge: 1-file 2-atom 10-const Default:if ISTART=0 2 else 0 ISPIN spin polarized calculation (2-yes 1-no) default 2 MAGMOM initial mag moment / atom Default NIONS*1 INIWAV initial electr wf. : 0-lowe 1-rand Default 1 only used for start jobs (ISTART=0) IDIPOL calculate monopole/dipole and quadrupole corrections 1-3 只计算第一/二/三晶矢方向适于slab的计算 4 全部计算尤其适于就算孤立分子 >PREC precession: medium, high or low(VASP.4.5+ also: normal, accurate) Default: Medium VASP4.5+采用了优化的accurate来替代high,所以一般不推荐使用 high。不过high可以确保'绝对收敛'，作为参考值有时也是必要的。 同样受推荐的是normal，作为日常计算选项，可惜的是说明文档提供的信息不足。 受PREC影响的参数有四类：ENCUT; NGX,NGY,NGZ; NGXF, NGYF, NGZF; ROPT 如果设置了PREC，这些参数就都不需要出现了 当然直接设置相应的参数也是同样效果的，这里不展开了，随后详释
锅炉计算简介
锅炉计算简介 锅炉热力计算结果是由锅炉制造厂家设计计算完成的。厂家为用户提供热力计算说明书或热力计算汇总表，其结果往往和实际运行的参数有差别。为便于电厂锅炉技术人员对锅炉技术问题进行分析，需要了解锅炉热力计算的基本思路。对受热面改造要依据《锅炉热力计算标准》进行设计。 第一节 燃烧计算 一、化学燃烧方程式 单位数量的燃料[固体及液体用千克（kg ）计，气体燃料用标准立方米（Nm 3）计]完全燃烧时所需要的空气量称为理论空气量。单位为标准立方米每千克（Nm 3/kg ），在此情况下空气中的氧全部与燃料中的可燃元素化合，烟气中没有自由氧存在，即1kg 燃料中的可燃成分100ar C kg ，100 ar H kg ，100ar S kg ，完全燃烧时所需空气量之和。 1.碳（C ） 碳完全燃烧时，化学反应式为 C +O 2 → CO 2 12kgC + 22.4 Nm 3 O 2 → 22.4 Nm 3CO 2 或 1kgC + 1.866 Nm 3 O 2 → 1.866Nm 3CO 2 1kg 收到基燃料中含有100 ar C kg 碳，因而1kg 燃料中的碳完全燃烧时所需要的氧量为1.866100 ar C Nm 3。 2.氢（H ） 氢完全燃烧时，化学反应式为 2H 2 + O 2 → 2H 2O 4.032kgH 2 + 22.4Nm 3O 2 → 44．8 Nm 3H 2O 或 1kgH 2 +5.56Nm 3O 2 →11.1Nm 3H 2O 1kg 收到基燃料中含有100 ar H kg 氢，因而1kg 燃料中的氢完全燃烧时所需要的氧量为
5.56100 ar H Nm 3。 3.硫（S ） 硫完全燃烧时，化学反应式为 S + O 2 → SO 2 32kgS + 22.4Nm 3O 2 → 22．4Nm 3SO 2 或 1kgS + 0.7Nm 3O 2 →0．7Nm 3SO 2 1kg 收到基燃料中含有100 ar S kg 硫，因而1kg 燃料中的硫完全燃烧时所需要的氧量为0．7100 ar S Nm 3。 二、空气量的计算 (一)理论空气需要量用0 V 表示―单位燃料中的可燃质完全燃烧，而且空气中没有剩余氧时所需要空气的体积。 由于1kg 燃料本身含有的氧量为100ar O kg ，相当于100 7.0100324.22ar ar O O =?Nm 3。 所以，1kg 收到基燃料燃烧所需的理论氧量02O V 为式（2-1） 20ar ar ar ar O C S H O V 1.8660.7 5.560.7100100100100=++- Nm 3/kg （2—1） 在干空气中氧的容积含量为21%，所以，1kg 收到基固体或液体燃料燃烧时所需的理论空气量V 0 为 0ar ar ar ar ar ar ar ar C S H O 1V (1.8660.7 5.560.7)0.21100100100100 0.0889C 0.0333S 0.265H 0.0333O =++-=++- Nm 3/kg （2—2） 或写成式（2-3）的形式 0ar ar ar ar V 0.0889C 0.265H 0.0333(S O )=++- Nm 3/kg （2—3） 式中 C AR ——碳的收到基百分含量，%； ar S ——硫的收到基百分含量，%； ar H ——氢的收到基百分含量，%； ar 0——氧的收到基百分含量，%。
期权价值敏感性——希腊字母汇总
第三章 期权敏感性（希腊字母） 顾名思义，期权敏感性是指期权价格受某些定价参数的变动而变动的敏感 程度，本章主要介绍期权价格对其四个参数（标的资产市场价格、到期时间、 波动率和无风险利率）的敏感性指标，这些敏感性指标也称作希腊值（Greeks ）。 每一个希腊值刻画了某个特定风险，如果期权价格对某一参数的敏感性为 零，可以想见，该参数变化时给期权带来的价格风险就为零。实际上，当我们 运用期权给其标的资产或其它期权进行套期保值时，一种较常用的方法就是分 别算出保值工具与保值对象两者的价值对一些共同的变量（如标的资产价格、 时间、标的资产价格的波动率、无风险利率等）的敏感性，然后建立适当数量 的证券头寸，组成套期保值组合，使组合中的保值工具与保值对象的价格变动 能相互抵消，也就是说让套期保值组合对该参数变化的敏感性变为零，这样就 能起到消除相应风险的套期保值的目的。 本章将主要介绍 Delta 、Gamma 、Vega 、Theta 、Rho 五个常用希腊字母。 符号 风险因素 量化公式 Gamma Γ 标的证券价格变化 Delta 变化/标的证券价格变化 Vega ν 波动率变化 权利金变化/波动率变化 Theta Θ 到期时间变化 权利金变化/到期时间变化 本章符号释义： T 为期权到期时间 S 为标的证券价格， S 0 为标的证券现价， S T 为标的证券行权时价格 K 为期权行权价格 σ 为标的证券波动率 r 为无风险利率 π t 为资产组合在 t 时刻的价值 N () 为标准正态分布的累积密度函数，可以查表或用计算机(如 Excel)求得 1
VASP-INCAR参数设置
V A S P-I N C A R参数设置-CAL-FENGHAI-(2020YEAR-YICAI)_JINGBIAN
1. 结构优化 (Opt) SYSTEM = opt ISTART = 0 INIWAV = 1 ICHARG = 2 ISPIN = 2 LREAL = Auto ENCUT = 400 PREC = high NSW= 600 NELM = 60 IBRION = 2 ISIF = 2 POTIM = 0.1 ALGO= Fast LVDW = .TRUE. EDIFF = 1E-5 EDIFFG = 1E-4 or -0.05 # 体系需计算TS时，全部结构优化EDIFFG均设置为-0.05 ISMEAR = 0 SIGMA = 0.2 LCHARG = .FALSE. LWAVE = .FALSE.
2. 过渡态搜索 (TS): 计算时先进行低精度计算，再进行高精度计算 SYSTEM= TS ISTART = 0 INIWAV = 1 ICHARG = 2 ISPIN = 2 LREAL = Auto ENCUT = 400 PREC = high NSW = 600 NELMIN = 6 IBRION = 3 or 1 # 过渡态计算低精度为3，高精度为1 ISIF = 2 POTIM = 0.01 ALGO = Fast LVDW = .TRUE. EDIFF = 1E-5 EDIFFG = -1 or -0.05 # 过渡态计算低精度为-1，高精度为-0.05 ISMEAR = 0 SIGMA = 0.05 LCHARG= .FALSE. LWAVE= .FALSE. IMAGES=8 # TS专属设置 SPRING=-5 # TS专属设置 LCLIMB=.TRUE. # TS专属设置
锅炉废气,计算方式总结
1、燃煤锅炉废气污染源强 ①过渡期燃煤锅炉废气污染源强 燃煤锅炉额定煤用量可根据下式计算： B＝(D×L)/(Q dw×η1) 式中：B——锅炉额定煤用量，t/h； D——锅炉每小时产汽量(根据型号为4t/h)； L——锅炉锅炉工作压力下饱和蒸气焓(查有关锅炉手册为659.9kcal/kg)； Q dw——燃煤的低位发热值，(取5800kcal/kg)； η 1——锅炉热效率(取80％)。 根据项目拟采用的龙岩无烟煤的煤质分析报告相关参数，含硫量0.72%，灰份21.8%，挥发份4.12%，低位发热值5487～6007kcal/kg。再由上式计算出一台4t/h蒸汽锅炉额定耗煤量为569kg/h。 a. 锅炉烟气排放量计算 燃煤锅炉的烟气量与锅炉型号、燃料的热值、燃烧方式以及配置的引风机型号均有密切关系。根据国家环境保护局科技标准司编写的《工业污染物产生和排 的产生量进行估算。放系数手册》计算公式对该项目燃煤烟气量、烟尘和SO 2 锅炉烟气量计算公式： Vo=1.01（QyL/1000）+0.5 （Nm3/ kg） Vy=0.89（QyL/1000）+1.65+（α-1）Vo （Nm3/ kg）式中：Vo——燃料燃烧所需理论空气量，Nm3/kg； Vy——实际烟气量，Nm3/kg； QyL——燃煤的低位发热值，(取5800kcal/kg)； α——炉膛过剩空气系数，α=αo+Δα，αo取1.3，Δα取0.5。
由上式可计算出Vo 为6.358Nm 3/kg ，Vy 为11.898Nm 3/kg ，4t/h 蒸汽锅炉额定耗煤量为569kg/h ，烟气排放量理论值为6770Nm 3/h ，即2.11×107 Nm 3/a 。 b. 烟尘产生量计算 G d ＝1000×A y ×a fh ×K C 11fh ?-）（ 式中：G d ——烟尘产污系数，kg/t-煤； A y ——煤中含灰量，%； a fh ——烟尘中飞灰占灰分总量的百分比； C fh ——烟尘中含碳量； K ——锅炉出力影响系数。 公式中各项参数，按调查资料和实测数据确定如下：A y =21.8%，a fh =10%， C fh =30%，K=1，计算出烟尘产污系数G d 为31.1 kg/t 。项目燃煤蒸汽锅炉每天运行12小时，每年运行260日，预计煤用量为1775t/a （569kg/h ），则烟尘产生量为55.2 t/a ，产生浓度为2616mg/Nm 3。 c. SO 2产生量计算 G SO 2 ＝1.6×B×S 式中：G SO 2 ——SO 2产生量，kg/h ； B ——耗煤量，kg/h ； S ——煤中全硫份。 经计算项目燃煤的SO 2产生量为20.45t/a ，产生浓度为969mg/Nm 3。 对照《锅炉大气污染物排放标准》（GB13271-2001）可知，燃煤烟气中SO 2和烟尘排放浓度均超过相应标准（900mg/Nm 3、200mg/Nm 3），项目拟采用文丘里—水膜（添加碱液）除尘设施对燃煤烟气进行净化，脱硫率按60%、除尘率按94%计，则SO 2和烟尘实际排放量分别为8.18t/a 和3.32t/a ，实际排放浓度分别为388mg/Nm 3和157mg/Nm 3。 ② 过渡期之后的运营期燃气锅炉废气污染源强
vasp计算参数设置
软件主要功能： 采用周期性边界条件（或超原胞模型）处理原子、分子、团簇、纳米线（或管）、薄膜、晶体、准晶和无定性材料，以及表面体系和固体 l 计算材料的结构参数（键长、键角、晶格常数、原子位置等）和构型 l 计算材料的状态方程和力学性质（体弹性模量和弹性常数） l 计算材料的电子结构（能级、电荷密度分布、能带、电子态密度和ELF） l 计算材料的光学性质 l 计算材料的磁学性质 l 计算材料的晶格动力学性质（声子谱等） l 表面体系的模拟（重构、表面态和STM模拟） l 从头分子动力学模拟 l 计算材料的激发态（GW准粒子修正） 计算主要的四个参数文件：INCAR ,POSCAR,POTCAR ,KPOINTS,下面简要介绍，详细权威的请参照手册 INCAR文件： 该文件控制VASP进行何种性质的计算，并设置了计算方法中一些重要的参数，这些参数主要包括以下几类： l 对所计算的体系进行注释：SYSTEM l 定义如何输入或构造初始的电荷密度和波函数：ISTART，ICHARG，INIWA V l 定义电子的优化 –平面波切断动能和缀加电荷时的切断值：ENCUT，ENAUG –电子部分优化的方法：ALGO，IALGO，LDIAG –电荷密度混合的方法：IMIX，AMIX，AMIN，BMIX，AMIX_MAG，BMIX_MAG，WC，INIMIX，MIXPRE，MAXMIX –自洽迭代步数和收敛标准：NELM，NELMIN，NELMDL，EDIFF l 定义离子或原子的优化 –原子位置优化的方法、移动的步长和步数：IBRION，NFREE，POTIM，NSW –分子动力学相关参数：SMASS，TEBEG，TEEND，POMASS，NBLOCK，KBLOCK，PSTRESS –离子弛豫收敛标准：EDIFFG l 定义态密度积分的方法和参数 –smearing方法和参数：ISMEAR，SIGMA –计算态密度时能量范围和点数：EMIN，EMAX，NEDOS –计算分波态密度的参数：RWIGS，LORBIT l 其它 –计算精度控制：PREC –磁性计算：ISPIN，MAGMOM，NUPDOWN –交换关联函数：GGA，VOSKOWN –计算ELF和总的局域势：LELF，LVTOT –结构优化参数：ISIF –等等。 主要参数说明如下： ? SYSTEM：该输入文件所要执行的任务的名字。取值：字符串，缺省值：SYSTEM
初学VASP中电子态密度计算设置参考
初学VASP中电子态密度计算基本设置参考主要分成三步：一、结构优化；二、静态自洽计算；三、非自洽计算以Al－FCC为例子 第一步结构优化 输入文件（INCAR, POTCAR, POSCAR, KPOINT） INCAR文件 System=Al ISTART=0 ISMEAR=1 SIGMA=0.2 ISPIN=2 GGA=91; VOSKOWN=1; EDIFF=0.1E-05; EDIFFG=-0.01 IBRION=2 NSW=50 ISIF=2 (OR 3) NPAR=10 POTCAR 文件直接在势库中拷贝 POSCAR文件 Al 4.05 1.0 0.0 0.0 0.0 1.0 0.0
0.0 0.0 1.0 4 Direct 0.0 0.0 0.0 0.5 0.5 0.0 0.5 0.0 0.5 0.0 0.5 0.5 KPOINT 文件 Automatic generation Mohkorst Pack 15 15 15 0.0 0.0 0.0 第二步静态自洽计算 INCAR：PREC = Medium，ISTART = 0，ICHARG = 2，ISMEAR = -5输入文件（INCAR, POTCAR, POSCAR, KPOINT） INCAR文件 System=Al ISTART=0 ISMEAR=1 SIGMA=0.2 ISPIN=2
GGA=91; VOSKOWN=1; EDIFF=0.1E-05; EDIFFG=-0.01 #IBRION=2 #NSW=50 #ISIF=2 (OR 3) NPAR=10 POTCAR 文件直接在势库中拷贝 POSCAR文件 Al 4.05 1.0 0.0 0.0 0.0 1.0 0.0 0.0 0.0 1.0 4 Selective Dynamic Direct 0.0 0.0 0.0 T T T 0.5 0.5 0.0 T T T 0.5 0.0 0.5 T T T 0.0 0.5 0.5 T T T KPOINT 文件 Automatic generation
VASP遇到小总结问题
VASP 计算的过程遇到的问题 01、第一原理计算的一些心得 （1）第一性原理其实是包括基于密度泛函的从头算和基于Hartree-Fock自洽计算的从头算，前者以电子密度作为基本变量（霍亨伯格-科洪定理），通过求解Kohn-Sham方程，迭代自洽得到体系的基态电子密度，然后求体系的基态性质；后者则通过自洽求解Hartree-Fock方程，获得体系的波函数，求基态性质； 评述：K-S方程的计算水平达到了H-F水平，同时还考虑了电子间的交换关联作用。 （2）关于DFT中密度泛函的Functional，其实是交换关联泛函 包括LDA，GGA，杂化泛函等等 一般LDA为局域密度近似，在空间某点用均匀电子气密度作为交换关联泛函的唯一变量，多数为参数化的CA-PZ方案； GGA为广义梯度近似，不仅将电子密度作为交换关联泛函的变量，也考虑了密度的梯度为变量，包括PBE,PW,RPBE等方案，BL YP泛函也属于GGA； 此外还有一些杂化泛函，B3L YP等。 （3）关于赝势 在处理计算体系中原子的电子态时，有两种方法，一种是考虑所有电子，叫做全电子法，比如WIEN2K中的FLAPW方法(线性缀加平面波)；此外还有一种方法是只考虑价电子，而把芯电子和原子核构成离子实放在一起考虑，即赝势法，一般赝势法是选取一个截断半径，截断半径以内，波函数变化较平滑，和真实的不同，截断半径以外则和真实情况相同，而且赝势法得到的能量本征值和全电子法应该相同。 赝势包括模守恒和超软，模守恒较硬，一般需要较大的截断能，超软势则可以用较小的截断能即可。另外，模守恒势的散射特性和全电子相同，因此一般红外，拉曼等光谱的计算需要用模守恒势。 赝势的测试标准应是赝势与全电子法计算结果的匹配度，而不是赝势与实验结果的匹配度，因为和实验结果的匹配可能是偶然的。 （4）关于收敛测试 （a）Ecut，也就是截断能，一般情况下，总能相对于不同Ecut做计算，当Ecut增大时总能变化不明显了即可；然而，在需要考虑体系应力时，还需对应力进行收敛测试，而且应力相对于Ecut的收敛要比总能更为苛刻，也就是某个截断能下总能已经收敛了，但应力未必收敛。 （b）K-point，即K网格，一般金属需要较大的K网格，采用超晶胞时可以选用相对较小的K网格，但实际上还是要经过测试。 （5）关于磁性 一般何时考虑自旋呢？举例子，例如BaTiO3中，Ba、Ti和O分别为+2，+4和-2价，离子全部为各个轨道满壳层的结构，就不必考虑自旋了；对于BaMnO3中，由于Mn+3价时d 轨道还有电子，但未满，因此需考虑Mn的自旋，至于Ba和O则不必考虑。其实设定自旋就是给定一个原子磁矩的初始值，只在刚开始计算时作为初始值使用，具体的可参照磁性物理。 （6）关于几何优化 包括很多种了，比如晶格常数和原子位置同时优化，只优化原子位置，只优化晶格常数，还有晶格常数和原子位置分开优化等等。
希腊字母在期权中的应用
希腊字母在期权中的应用 在衡量期权组合风险的时候，若用希腊字母来表示期权的风险指标，原本繁多复杂的 期权交易和持仓就会显得简洁明了。在交易中，投资者不仅要关注做多做空多少手不同的 期权合约，而且还要注意所有持仓的Delta、Gamma等参数。 选择策略 以最简单的买入标的和单腿策略为例，预计标的价格上涨，想要做多Delta，有买入 期货、买入看涨期权和卖出看跌期权三种方法，但预计标的价格上涨的同时波动率下跌， 即需要做多Delta、做空Vega，那么卖出看跌期权则是相对有利的策略。 对冲期权 对于同一个品种的期货和期权，希腊字母都可以直接相加减。当投资者利用跨式策略、价差策略、蝶式策略等多腿策略来交易期权时，有时候固定的策略并不能完美贴合投资者 的交易需求，此时就可以根据叠加后的希腊字母总和去对冲存在风险的部分。 例如，当预计市场有重大消息披露、标的价格可能有大幅变化、波动率将会变大时， 通常可以利用买入平值跨式期权策略来做多波动率。比如说，当豆粕期货1901合约价格 为3111元/吨时，同时买入行权价为3100元/吨的看涨期权和看跌期权构建买入跨式期 权策略。 可以看到这个策略中，两个期权的Delta并没有完全对冲掉，还存在一小部分方向上 的风险，当标的价格下跌时，会对这个跨式组合策略造成不利影响。此时可以做空0.073 倍的期货，得到-0.073个Delta，使得期权部位的总Delta为零。 管理持仓 由于希腊字母可以直接相加减，当持有的期权合约类型、行权价、数量等各不相同时，可以通过计算持仓部位的希腊字母来管理持仓风险。因此，即使持仓的头寸繁多复杂，利 用希腊字母的叠加，持仓的风险状况就会变得更直观明了，分析起来也更方便。 下面以铜期权2018年9月21日收盘时的风险参数为例，假设同时持有数量不一、 行权价不同的若干期权，结果如下表所示：
初学VASP最重要的INCAR参数
初学VASP(六) 最重要的INCAR参数 初学VASP(六) 最重要的INCAR参数 INCAR是决定how to do 的文件 限于能力，只对部分最基本的一些参数(>,没有这个标志的参数都是可以不出现的) 详细说明，在这里只是简单介绍这些参数的设置，详细的问题在后文具体示例中展开。 部分可能会干扰VASP运行的参数在这里被刻意隐去了，需要的同学还是请查看VASP自带 的帮助文档原文。 参数列表如下： >SYSTEM name of System 任务的名字 *** >NWRITE verbosity write-flag (how much is written) 输出内容详细程度 0-3 缺省2 如果是做长时间动力学计算的话最好选0或1(首末步/每步核运动输出) 据说也可以结合shell的tail或grep命令手动输出 >ISTART startjob: restart选项 0-3 缺省0/1 for 无/有前次计算的WAVECAR(波函数) 1 'restart with constant energy cut-off' 2 'restart with constant basis set' 3 'full restart including wave function and charge prediction' ICHARG charge: 1-file 2-atom 10-const Default:if ISTART=0 2 else 0 ISPIN spin polarized calculation (2-yes 1-no) default 2 MAGMOM initial mag moment / atom Default NIONS*1 INIWAV initial electr wf. : 0-lowe 1-rand Default 1 only used for start jobs (ISTART=0) IDIPOL calculate monopole/dipole and quadrupole corrections 1-3 只计算第一/二/三晶矢方向适于slab的计算 4 全部计算尤其适于就算孤立分子 >PREC precession: medium, high or low(VASP.4.5+ also: normal, accurate)
VASP-INCAR参数设置
1. 结构优化(Opt) SYSTEM = opt ISTART = 0 INIWAV = 1 ICHARG = 2 ISPIN = 2 LREAL = Auto ENCUT = 400 PREC = high NSW= 600 NELM = 60 IBRION = 2 ISIF = 2 POTIM = 0.1 ALGO= Fast LVDW = .TRUE. EDIFF = 1E-5 EDIFFG = 1E-4 or -0.05 # 体系需计算TS时，全部结构优化EDIFFG均设置为-0.05 ISMEAR = 0 SIGMA = 0.2 LCHARG = .FALSE. LWAVE = .FALSE.
2. 过渡态搜索(TS): 计算时先进行低精度计算，再进行高精度计算 SYSTEM= TS ISTART = 0 INIWAV = 1 ICHARG = 2 ISPIN = 2 LREAL = Auto ENCUT = 400 PREC = high NSW = 600 NELMIN = 6 IBRION = 3 or 1 # 过渡态计算低精度为3，高精度为1 ISIF = 2 POTIM = 0.01 ALGO = Fast LVDW = .TRUE. EDIFF = 1E-5 EDIFFG = -1 or -0.05 # 过渡态计算低精度为-1，高精度为-0.05 ISMEAR = 0 SIGMA = 0.05 LCHARG= .FALSE. LWAVE= .FALSE. IMAGES=8 # TS专属设置 SPRING=-5 # TS专属设置 LCLIMB=.TRUE. # TS专属设置
锅炉热力计算
锅炉热力计算 ●计算依据 燃煤热值按4500千卡/公斤、醇基燃料热值按6500千卡/公斤、柴油热值按10200千卡/公斤，燃煤价格按750元/吨、醇基燃料按3500元/吨、柴油价格按7500元/吨，煤锅炉的效率按45%、油气锅炉的效率按95%计算： ●4吨燃油蒸汽锅炉 4吨燃油蒸汽锅炉的热功率为248万大卡/小时， * 使用燃煤蒸汽锅炉，使用成本为： 248×104÷4500÷45%=1225公斤/小时×0.75=919元/小时*换装燃醇蒸汽锅炉使用醇基燃料使用成本为： 248×104÷6500÷95%=401公斤/小时×3.5=1404元/小时*换装油气蒸汽锅炉使用柴油作为燃料的使用成本为： 248×104÷10200÷95%=256公斤/小时×7.5=1920元/小时 ●300万大卡导热油锅炉 *使用燃煤导热油锅炉，使用成本为：
300×104÷4500÷45%=1482公斤/小时×0.75=1112元/小时*换装燃醇导热油锅炉使用醇基燃料使用成本为： 300×104÷6500÷95%=486公斤/小时×3.5=1700元/小时*换装油气导热油锅炉使用柴油作为燃料的使用成本为： 300×104÷10200÷95%=310公斤/小时×7.5=2325元/小时 三、综合效益计算 1、设备成本 ●4吨蒸汽锅炉 沿用现有的燃煤锅炉使用醇基燃料，每小时使用成本为： 248×104÷6500÷95%×3.5=1404元/小时 每天按8小时计算，则每天为11232元。 若更换同等功率的燃油燃气蒸汽锅炉约需55万元，每小时使用成本为1920元，每天按8小时计算，则每天为15360元，每天节省燃料费3984元，约130天即可收回设备投入。 ●300万大卡导热油锅炉
vasp在计算磁性的实例和讨论
兄弟，问3个问题 1，vasp在计算磁性的时候，oszicar中得到的磁矩和outcar中得到各原子磁矩之和不一致，在投稿的是否曾碰到有审稿人质疑，对于这个不一致你们一般是怎么解释的了？ 2，另外，磁性计算应该比较负责。你应该还使用别的程序计算过磁性，与vasp结果比较是否一致，对磁性计算采用的程序有什么推荐。 ps：由于曾使用vasp和dmol算过非周期体系磁性，结构对磁性影响非常大，因此使用这两个程序计算的磁性要一致很麻烦。还不敢确定到底是哪个程序可能不可靠。 3，如果采用vasp计算磁性，对采用的方法和设置有什么推荐。 1，OSZICAR中得到的磁矩是OUTCAR中最后一步得到的总磁矩是相等的。总磁矩和各原子的磁矩(RMT球内的磁矩)之和之差就是间隙区的磁矩。因为有间隙区存在,不一致是正常的。 2，如果算磁性，全电子的结果更精确，我的一些计算结果显示磁性原子对在最近邻的位置时，PAW与FPLAW给出的能量差不一致，在长程时符合的很好。虽然并没有改变定性结论。感觉PAW似乎不能很好地描述较强耦合。我试图在找出原因，主要使用exciting和vasp做比较。计算磁性推荐使用FP-LAPW, FP-LMTO, FPLO很吸引人(不过是商业的)，后者是O(N)算法。 3，使用vasp计算磁性，注意不同的初始磁矩是否收敛为同一个磁矩。倒没有特别要注意的地方，个人认为。 归根结底，需要一个优秀的交换关联形式出现 VASP计算是否也是像计算DOS和能带一样要进行三步（结构优化，静态自洽计算，非自洽计算），然后看最后一步的出的磁矩呢？ 一直想计算固体中某个原子的磁矩，根据OUTCAR的结果似乎不能分析，因为它里面总磁矩跟OSZICAR的值有一定的差别，据说是OUTCAR中只考虑WS半径内磁矩造成的。最近看到一个帖子说是可以用bader电荷分析方法分析原子磁矩。如法炮制之后发现给出的总磁矩与OSZICAR的结果符合的甚好，可是觉得没有根据，有谁知道这样做的依据吗，欢迎讨论！ 设置ISPIN=2计算得到的态密度成为自旋态密度。 设置ISPIN=2就可以计算磁性，铁磁和反铁磁在MAG里设置。最后得到的DOS是分up和down的。 磁性计算 (2006-12-03 21:02) 标签： - 分类：Vasp ·磁性计算
锅炉计算公式汇总
锅炉计算公式： 1、 x k 1= x 为干度 2、a a o ) 1(212-= 3、a o o ?=-- -' 2 ' '2 21212121 4、 b B Q B =2 .29271. 5、y y y y O H S C V 0333.0265.0)375.0(0889.00-++= 6、风机每小时耗电量KWh/锅炉实际负荷t/h 7、)/(08 .01kwh g b b g -=发 8、%.) ()(1 11 net r gs o Q B h h D h h D -+-=η 9、28.024 3600.?= net Q B P ≈200MW 10、% 3076.0538.12.05 .7211110 %2.05=?=?= q 11、)%(10065432q q q q q ++++-=η 12、%100.100.32866.4 ??? ? ????-+-=LZ LZ LZ fh fh fh net ar c c a c c a Q A q 44.q B B = （t/h ） 13、%79.89.)(=-= net gs o gL Q B t h D η 14、kwh t P w 4410144007201020.?=??== 15、 zf g f x B D p ηηηη+++=∑)( 16、 % ..3600net gL jx Q B b p ∑-= ηη
%100271.293600 ??= b η )1(n b b ?-=发 % 100271.29??= 发 b q gL η 17、)/() (h kg Q D B net η 给水焓过热蒸汽焓 -= 1 2 11 2 2121212121o o a a a -- -= -=? BV V a ?=?（M 3/h ） 18、D PW =PD 19、1mm Hg=133.3P a P g =P-P 0 20、H 2=150+60=210mm ∵22 11υυh h = ∴mm h h 2800012 .00016.02102 1 21=?= ?= υυ h=280-（150+50）=80mm 21、2 1 2111/) (m w t t q δλ-= 2 2 21 1212/) (m w t t q λδλδ+ -= 22、求热伸长值？ )(mm t La L L ?=? 23、 24、 [ ] kg kj s o H c Q y y y y gr /186.4)(2630081?--+= [ ]186.4) 654(?+-=y y gr net M H Q Q ) /(.3 h m V Q B q net v =