24小时热门版块排行榜

返回列表

本帖产生 1 个模拟EPI ，点击这里进行查看

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

zyj8119

木虫 (著名写手)

模拟EPI: 10
应助: 65 (初中生)
贵宾: 0.003
金币: 915.1
散金: 1440
红花: 35
帖子: 2936
在线: 1329.4小时
虫号: 664177
注册: 2008-11-29
性别: GG
专业: 理论和计算化学

[交流] 【原创】MCM-41的孔径分布的计算已有4人参与

接上一个帖子，MCM-41建好的模型会有孔径分布，这个在定义中是很复杂的，应该怎么计算呢，可以导出MCM-41中所有的H原子，然后自己写程序去计算孔径分布，我使用的是FORTRAN语言，导出的car结构中，可以使用EXCEL找到所有的H原子，这些H原子的 car文件可以保存为一个名叫H.car的文件，这个文件就可以进行孔径分布的计算，car文件的优点是带有电荷。读入这个car文件，可以输出一个名叫output.car的文件。

CODE:

        integer hno

        parameter(hno=2356)

        character a(hno)*5,fx(hno)*4,fft(hno)*5,natom(hno)*1

        double precision xc(hno),yc(hno),zc(hno)

        integer occupation(hno)

        real charge(hno)

      open(10,file='H.car',status='old')

        do 20 i=1,hno

        read(10,*)a(i),xc(i),yc(i),zc(i),fx(i),occupation(i),fft(i),

     &          natom(i),charge(i)

20    continue

      close(10)

        jmm=0

      open(30,file='output.car',access='append')

*40        amm=RAN2(IDUM)

40        imm=int(RAN2(IDUM)*hno)

        if(charge(imm).lt.1.0)then

        jmm=jmm+1

        write(30,888)a(imm),xc(imm),yc(imm),zc(imm),fx(imm),

     &             occupation(imm),fft(imm),natom(imm),charge(imm)

        charge(imm)=2.0

        endif

        if(jmm.eq.50.or.jmm.eq.100.or.jmm.eq.150.or.jmm.eq.200.or.

     &          jmm.eq.250.or.jmm.eq.300.or.jmm.eq.350.or.jmm.eq.400)then

        write(30,*)'********************',jmm,'**************************'

        write(30,*)'********************',jmm,'**************************'

        endif

        if(jmm.lt.415)then

      goto 40

        else

        goto 50

        endif

50    close(30)

888   format(A5,3X,F12.9,2X,F13.9,3X,F12.9,1X,A4,1X,I1,6X,A5,3X,A1,F8.3)   

        end 

*&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&

*     Random Number Generator                

*&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&

                          FUNCTION RAN2(IDUM)

      INTEGER idum,IM1,IM2,IMM1,IA1,IA2,IQ1,IQ2,IR1,IR2,NTAB,NDIV

      REAL ran2,AM,EPS,RNMX

      PARAMETER (IM1=2147483563,IM2=2147483399,AM=1./IM1,IMM1=IM1-1,

     *IA1=40014,IA2=40692,IQ1=53668,IQ2=52774,IR1=12211,IR2=3791,

     *NTAB=32,NDIV=1+IMM1/NTAB,EPS=1.2e-7,RNMX=1.-EPS)

        INTEGER idum2,j,k,iv(NTAB),iy

      SAVE iv,iy,idum2

      DATA idum2/123456789/, iv/NTAB*0/, iy/0/

      if (idum.le.0) then

        idum=max(-idum,1)

        idum2=idum

        do 11 j=NTAB+8,1,-1

          k=idum/IQ1

          idum=IA1*(idum-k*IQ1)-k*IR1

          if (idum.lt.0) idum=idum+IM1

          if (j.le.NTAB) iv(j)=idum

11      continue

        iy=iv(1)

      endif

      k=idum/IQ1

      idum=IA1*(idum-k*IQ1)-k*IR1

      if (idum.lt.0) idum=idum+IM1

      k=idum2/IQ2

      idum2=IA2*(idum2-k*IQ2)-k*IR2

      if (idum2.lt.0) idum2=idum2+IM2

      j=1+iy/NDIV

      iy=iv(j)-idum2

      iv(j)=idum

      if(iy.lt.1)iy=iy+IMM1

      ran2=min(AM*iy,RNMX)

      return

      END

其中的nho是你的H.car中H原子的个数，可以按照你的文件的不同随意更改成你所需要的值，得到的output.car可以做下一步的计算。

[ Last edited by zyj8119 on 2010-12-2 at 14:14 ]

回复此楼

» 猜你喜欢

博士读完未来一定会好吗已经有6人回复
小论文投稿已经有3人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有9人回复
心脉受损已经有3人回复
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有8人回复
申请2026年博士已经有6人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有5人回复

» 本主题相关价值贴推荐，对您同样有帮助:

【求助】氮吸附曲线已经有17人回复
【求助】为什么我在不同温度下制备的MCM-41的小角XRD和孔径都一样，但BET和孔容差别很已经有5人回复
【讨论】关于N2吸附脱附的问题问题. 已经有9人回复
【求助】BET 孔径计算已经有10人回复
【求助】如何计算Cu-MCM-41的晶胞常数已经有3人回复
【求助】mcm-41的xrd，求高手指点已经有32人回复
【求助成功】MS模拟MCM-41过程中forcite filed计算的问题已经有12人回复

好好学习，天天向上。

1楼2010-12-02 00:24:27

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

zyj8119

木虫 (著名写手)

模拟EPI: 10
应助: 65 (初中生)
贵宾: 0.003
金币: 915.1
散金: 1440
红花: 35
帖子: 2936
在线: 1329.4小时
虫号: 664177
注册: 2008-11-29
性别: GG
专业: 理论和计算化学

引用回帖:

Originally posted by zyj8119 at 2010-12-02 10:01:19:
统计完成以后，接下来就是计算孔径分布了，源程序是这样的：

[code]
*             main program
            real cx,cy
            integer HNO,lattice
            open(1000,file='para ...

版主请评阅下。。。

回复此楼